Detailed Description

Modify masses used in multi-gauss fits of an id-grid to be conform to those of a 'reference' grid.

This class will take one KVIDZAFromZGrid as "reference" and modify the masses of the all the KVMultiGaussIsotopeFit fits in a second KVIDZAFromZGrid so that they are consistent with the "reference" grid.

Example

KVIDZAFromZGrid* ref_gr; // reference grid pointer
KVIDZAFromZGrid* mod_gr; // pointer to grid to check/modify
 
KVIDZAFromZGridMassCorrector corr(ref_gr,mod_gr);
 
auto corrections = corr.DeterminMassCorrections(0.8); // use all isotopes with >80% of max yield
 
TH2* hdata; // 2d identification map used with grid to check/modify
auto pid_dist = mod_gr->LinearizeHistoToPID(hdata); // generate 1d PID distribution for check/modify grid
 
corr.ApplyCorrections(corrections, pid_dist, KVItvFinderDialog::mass_fit_parameters); // change masses of all isotopes according to determined corrections

Note: For reasons beyond our ken, the default parameters for fitting the gaussians to PID spectra are part of the GUI class KVItvFinderDialog. They should obviously be part of KVMultiGaussIsotopeFit but this causes problems with ROOT/Cling which are for the moment incomprehensible.

Definition at line 35 of file KVIDZAFromZGridMassCorrector.h.

#include <KVIDZAFromZGridMassCorrector.h>

Inheritance diagram for KVIDZAFromZGridMassCorrector:

[legend]

Public Types
using	mass_correction_t = std::pair< int, int >

Public Types inherited from KVBase
enum	EKaliVedaBits { kIsKaliVedaObject = BIT(23) }

enum class	GetX_status { ok , above_maximum , below_minimum , at_maximum , at_minimum }

Public Types inherited from TObject
enum	EDeprecatedStatusBits

enum	EStatusBits

Public Member Functions
	KVIDZAFromZGridMassCorrector (KVIDZAFromZGrid _reference, KVIDZAFromZGrid _to_correct)

void	ApplyCorrections (const std::vector< KVIDZAFromZGridMassCorrector::mass_correction_t > &, TH1 *pid_dist, const KVNameValueList &mass_fit_parameters)

std::vector< mass_correction_t >	DetermineMassCorrections (double min_yield)

Public Member Functions inherited from KVBase
	KVBase ()
	Default constructor. More...

	KVBase (const Char_t name, const Char_t title="")
	Ctor for object with given name and type. More...

	KVBase (const KVBase &)
	copy ctor More...

void	Clear (Option_t *opt="") override
	Clear object properties : name, type/title, number, label. More...

void	Copy (TObject &) const override
	Make a copy of this object. More...

const Char_t *	GetLabel () const

UInt_t	GetNumber () const

virtual TObject *	GetObject () const

virtual const Char_t *	GetType () const

Bool_t	HasLabel () const

Bool_t	IsCalled (const Char_t *name, bool case_insensitive=false) const

Bool_t	IsLabelled (const Char_t *l, bool case_insensitive=false) const

Bool_t	IsType (const Char_t *typ, bool case_insensitive=false) const

virtual void	List ()

KVBase &	operator= (const KVBase &)
	copy assignment operator More...

void	Print (Option_t *option="") const override

void	SetLabel (const Char_t *lab)

virtual void	SetNumber (UInt_t num)

virtual void	SetType (const Char_t *str)

Public Member Functions inherited from TNamed
	TNamed ()

	TNamed (const char name, const char title)

	TNamed (const TNamed &named)

	TNamed (const TString &name, const TString &title)

virtual	~TNamed ()

TObject *	Clone (const char *newname="") const override

Int_t	Compare (const TObject *obj) const override

virtual void	FillBuffer (char *&buffer)

const char *	GetName () const override

const char *	GetTitle () const override

ULong_t	Hash () const override

TClass *	IsA () const override

Bool_t	IsSortable () const override

void	ls (Option_t *option="") const override

TNamed &	operator= (const TNamed &rhs)

virtual void	SetName (const char *name)

virtual void	SetNameTitle (const char name, const char title)

virtual void	SetTitle (const char *title="")

virtual Int_t	Sizeof () const

void	Streamer (TBuffer &) override

void	StreamerNVirtual (TBuffer &ClassDef_StreamerNVirtual_b)

Public Member Functions inherited from TObject
	TObject ()

	TObject (const TObject &object)

virtual	~TObject ()

void	AbstractMethod (const char *method) const

virtual void	AppendPad (Option_t *option="")

virtual void	Browse (TBrowser *b)

ULong_t	CheckedHash ()

virtual const char *	ClassName () const

virtual void	Delete (Option_t *option="")

virtual Int_t	DistancetoPrimitive (Int_t px, Int_t py)

virtual void	Draw (Option_t *option="")

virtual void	DrawClass () const

virtual TObject *	DrawClone (Option_t *option="") const

virtual void	Dump () const

virtual void	Error (const char method, const char msgfmt,...) const

virtual void	Execute (const char method, const char params, Int_t *error=nullptr)

virtual void	Execute (TMethod method, TObjArray params, Int_t *error=nullptr)

virtual void	ExecuteEvent (Int_t event, Int_t px, Int_t py)

virtual void	Fatal (const char method, const char msgfmt,...) const

virtual TObject *	FindObject (const char *name) const

virtual TObject *	FindObject (const TObject *obj) const

virtual Option_t *	GetDrawOption () const

virtual const char *	GetIconName () const

virtual char *	GetObjectInfo (Int_t px, Int_t py) const

virtual Option_t *	GetOption () const

virtual UInt_t	GetUniqueID () const

virtual Bool_t	HandleTimer (TTimer *timer)

Bool_t	HasInconsistentHash () const

virtual void	Info (const char method, const char msgfmt,...) const

virtual Bool_t	InheritsFrom (const char *classname) const

virtual Bool_t	InheritsFrom (const TClass *cl) const

virtual void	Inspect () const

void	InvertBit (UInt_t f)

Bool_t	IsDestructed () const

virtual Bool_t	IsEqual (const TObject *obj) const

virtual Bool_t	IsFolder () const

R__ALWAYS_INLINE Bool_t	IsOnHeap () const

R__ALWAYS_INLINE Bool_t	IsZombie () const

void	MayNotUse (const char *method) const

virtual Bool_t	Notify ()

void	Obsolete (const char method, const char asOfVers, const char *removedFromVers) const

void	operator delete (void *ptr)

void	operator delete (void ptr, void vp)

void	operator delete[] (void *ptr)

void	operator delete[] (void ptr, void vp)

void *	operator new (size_t sz)

void *	operator new (size_t sz, void *vp)

void *	operator new[] (size_t sz)

void *	operator new[] (size_t sz, void *vp)

TObject &	operator= (const TObject &rhs)

virtual void	Paint (Option_t *option="")

virtual void	Pop ()

virtual Int_t	Read (const char *name)

virtual void	RecursiveRemove (TObject *obj)

void	ResetBit (UInt_t f)

virtual void	SaveAs (const char filename="", Option_t option="") const

virtual void	SavePrimitive (std::ostream &out, Option_t *option="")

void	SetBit (UInt_t f)

void	SetBit (UInt_t f, Bool_t set)

virtual void	SetDrawOption (Option_t *option="")

virtual void	SetUniqueID (UInt_t uid)

void	StreamerNVirtual (TBuffer &ClassDef_StreamerNVirtual_b)

virtual void	SysError (const char method, const char msgfmt,...) const

R__ALWAYS_INLINE Bool_t	TestBit (UInt_t f) const

Int_t	TestBits (UInt_t f) const

virtual void	UseCurrentStyle ()

virtual void	Warning (const char method, const char msgfmt,...) const

virtual Int_t	Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0)

virtual Int_t	Write (const char *name=nullptr, Int_t option=0, Int_t bufsize=0) const

Static Public Attributes
static int	fOnlyShowZ = 0

Additional Inherited Members
Static Public Member Functions inherited from KVBase
static std::optional< TString >	AbsoluteUnixPath (const TString &)

static Bool_t	AreEqual (Double_t x, Double_t y, Long64_t maxdif=1)
	Comparison between two 64-bit floating-point values. More...

static void	BackupFileWithDate (const Char_t *path)

static void	CombineFiles (const Char_t file1, const Char_t file2, const Char_t *newfilename, Bool_t keep=kTRUE)

static void	Deprecated (const char method, const char advice)

static Bool_t	FindClassSourceFiles (const KVString &class_name, KVString &imp_file, KVString &dec_file, const KVString &dir_name=".")

static Bool_t	FindExecutable (TString &exec, const Char_t *path="$(PATH)")

static const Char_t *	FindFile (const Char_t *search, TString &wfil)

static const Char_t *	GetBINDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetDATABASEFilePath ()

static const Char_t *	GetDATADIRFilePath (const Char_t *namefile="")

template<typename ValType >
static ValType	GetDataSetEnv (const KVString &dataset, const KVString &type, const ValType &defval={})

static const Char_t *	GetETCDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetExampleFilePath (const Char_t library, const Char_t namefile)
	Return full path to example file for given library (="KVMultiDet", "BackTrack", etc.) More...

static const Char_t *	GetINCDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetKVBuildDate ()
	Returns KaliVeda build date. More...

static const Char_t *	GetKVBuildDir ()
	Returns top-level directory used for build. More...

static const Char_t *	GetKVBuildTime ()
	Returns KaliVeda build time. More...

static const Char_t *	GetKVBuildType ()
	Returns KaliVeda build type (cmake build: Release, Debug, RelWithDebInfo, ...) More...

static const Char_t *	GetKVBuildUser ()
	Returns username of person who performed build. More...

static int	GetKVMajorVersion ()

static int	GetKVMinorVersion ()

static int	GetKVPatchVersion ()

static const Char_t *	GetKVSourceDir ()
	Returns top-level directory of source tree used for build. More...

static const Char_t *	GetKVVersion ()
	Returns KaliVeda version string. More...

static const Char_t *	GetLIBDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetListOfPlugins (const Char_t *base)

static const Char_t *	GetListOfPluginURIs (const Char_t *base)

static const Char_t *	GetPluginURI (const Char_t base, const Char_t plugin)

static void	GetTempFileName (TString &base)

static const Char_t *	GetTEMPLATEDIRFilePath (const Char_t *namefile="")

static const Char_t *	GetWORKDIRFilePath (const Char_t *namefile="")

static const Char_t *	gitBranch ()
	Returns git branch of sources. More...

static const Char_t *	gitCommit ()
	Returns last git commit of sources. More...

static void	InitEnvironment ()

static bool	is_gnuinstall ()

static Bool_t	IsEnvInit ()

static Bool_t	IsThisAPlugin (const TString &uri, TString &base)

static TPluginHandler *	LoadPlugin (const Char_t base, const Char_t uri="0")

static Bool_t	OpenContextMenu (const char method, TObject obj, const char *alt_method_name="")

static void	OpenTempFile (TString &base, std::ofstream &fp)

static void	PrintSplashScreen ()
	Prints welcome message and infos on version etc. More...

static Double_t	ProtectedGetX (const TF1 &func, Double_t val, GetX_status &status, std::optional< Double_t > xmin={}, std::optional< Double_t > xmax={})

static Double_t	ProtectedGetX (const TF1 *func, Double_t val, GetX_status &status, std::optional< Double_t > xmin={}, std::optional< Double_t > xmax={})

static Bool_t	SearchAndOpenKVFile (const Char_t name, KVSQLite::database &dbfile, const Char_t kvsubdir="")

static Bool_t	SearchAndOpenKVFile (const Char_t name, std::ifstream &file, const Char_t kvsubdir="", KVLockfile *locks=0)

static Bool_t	SearchAndOpenKVFile (const Char_t name, std::ofstream &file, const Char_t kvsubdir="", KVLockfile *locks=0)

static Bool_t	SearchKVFile (const Char_t name, TString &fullpath, const Char_t kvsubdir="")

static const Char_t *	WorkingDirectory ()

Static Public Member Functions inherited from TNamed
static TClass *	Class ()

static const char *	Class_Name ()

static constexpr Version_t	Class_Version ()

static const char *	DeclFileName ()

Static Public Member Functions inherited from TObject
static TClass *	Class ()

static const char *	Class_Name ()

static constexpr Version_t	Class_Version ()

static const char *	DeclFileName ()

static Longptr_t	GetDtorOnly ()

static Bool_t	GetObjectStat ()

static void	SetDtorOnly (void *obj)

static void	SetObjectStat (Bool_t stat)

Public Attributes inherited from TObject
	kBitMask

	kCanDelete

	kCannotPick

	kHasUUID

	kInconsistent

	kInvalidObject

	kIsOnHeap

	kIsReferenced

	kMustCleanup

	kNoContextMenu

	kNotDeleted

	kObjInCanvas

	kOverwrite

	kSingleKey

	kWriteDelete

	kZombie

Protected Member Functions inherited from TObject
virtual void	DoError (int level, const char location, const char fmt, va_list va) const

void	MakeZombie ()

Protected Attributes inherited from TNamed
TString	fName

TString	fTitle

Protected Attributes inherited from TObject
	kOnlyPrepStep

Member Typedef Documentation

◆ mass_correction_t

using KVIDZAFromZGridMassCorrector::mass_correction_t = std::pair<int,int>

Definition at line 51 of file KVIDZAFromZGridMassCorrector.h.

Constructor & Destructor Documentation

◆ KVIDZAFromZGridMassCorrector()

KVIDZAFromZGridMassCorrector::KVIDZAFromZGridMassCorrector	(	KVIDZAFromZGrid *	_reference,
		KVIDZAFromZGrid *	_to_correct
	)

inline

Definition at line 46 of file KVIDZAFromZGridMassCorrector.h.

Member Function Documentation

◆ ApplyCorrections()

void KVIDZAFromZGridMassCorrector::ApplyCorrections	(	const std::vector< KVIDZAFromZGridMassCorrector::mass_correction_t > &	corrections,
		TH1 *	pid_dist,
		const KVNameValueList &	mass_fit_parameters
	)

Takes the mass corrections previously determined by DetermineMassCorrections() and applies them to each fit using KVGaussFitMassModifier

Parameters

[in]	corrections	see DetermineMassCorrections()
[in]	pid_dist	a PID distribution corresponding to the grid being used
[in]	mass_fit_parameters	parameters to control fit; by default KVMultiGaussIsotopeFit::mass_fit_parameters

Note: For any KVIDZAGrid, a 1-dimensional PID distribution to be used with multi-gaussian fits can be generated from a 2D identification map (e.g. $\Delta E-E$ matrix) using KVIDZAFromZGrid::LinearizeHistoToPID()

Definition at line 148 of file KVIDZAFromZGridMassCorrector.cpp.

◆ DetermineMassCorrections()

std::vector< KVIDZAFromZGridMassCorrector::mass_correction_t > KVIDZAFromZGridMassCorrector::DetermineMassCorrections ( double _min_yield )

Determine the mass-shifts to apply to each fitted element $Z$ by comparing the mass distributions with those of the reference fits.

Parameters

[in] _min_yield minimum yield (expressed as a fraction of the largest yield in the reference fit) above which to consider isotopes in mass distributions

Note: For $Z=6$ _min_yield is fixed at 0.5

Up to 3 isotope masses will be compared for the 2 fits of each element, depending on the threshold minimum yield parameter.

If the sequence of masses for the reference & other fit can be related with a single offset (i.e. masses "24,25,26" for one, masses "23,24,25" for the other implies an offset of $\pm 1$), then this is taken to be the required mass correction, unless such a correction would lead to assign a spurious value corresponding to an undetectable resonance etc.

Note: set static member variable KVIDZAFromZGridMassCorrector::fOnlyShowZ to a non-zero value in order to only consider the given $Z$. Reset KVIDZAFromZGridMassCorrector::fOnlyShowZ=0 in order to treat all elements.

Returns: a vector containing the mass correction required for each $Z$

Definition at line 27 of file KVIDZAFromZGridMassCorrector.cpp.

Member Data Documentation

◆ fOnlyShowZ

int KVIDZAFromZGridMassCorrector::fOnlyShowZ = 0

static

Definition at line 44 of file KVIDZAFromZGridMassCorrector.h.