Detailed Description

Check & correct if necessary that isotopic masses used in different multi-gauss fits are consistent.

Use to determine whether or not attributed masses are consistent for a set of KVIDZAFromZGrid using multi-gauss isotope fits, and correct where necessary.

Example

KVIDZAFromZGridMassHomogenizer gmh("my_grids.dat", "my_data.root");
 
gmh.Analyse(0.8);
gmh.Print();
 
...[example output]...
Z = 12 :
      25 : 1    {ID_SI1_SI2_443}
        26 : 2  {ID_SI1_SI2_112, ID_SI1_SI2_143}
        25-26 : 161 {ID_SI1_SI2_1011, ID_SI1_SI2_1012, ID_SI1_SI2_1013, ID_SI1_SI2_1014, ID_SI1_SI2_1021, ... }
Z = 13 :
    27-28 : 71  {ID_SI1_SI2_1012, ID_SI1_SI2_1021, ID_SI1_SI2_1032, ID_SI1_SI2_1034, ID_SI1_SI2_111, ... }
    26-27 : 92  {ID_SI1_SI2_1011, ID_SI1_SI2_1013, ID_SI1_SI2_1014, ID_SI1_SI2_1022, ID_SI1_SI2_1023, ... }
Z = 14 :
    29-30 : 20  {ID_SI1_SI2_1114, ID_SI1_SI2_411, ID_SI1_SI2_412, ID_SI1_SI2_442, ID_SI1_SI2_511, ... }
    30 : 142    {ID_SI1_SI2_1011, ID_SI1_SI2_1012, ID_SI1_SI2_1013, ID_SI1_SI2_1014, ID_SI1_SI2_1021, ... }
...
 
gmh.CorrectMasses(13, "26-27", 1, KVItvFinderDialog::mass_fit_parameters); // increase A by 1 for Z=13 isotopes in all grids in list {ID_SI1_SI2_1011, ID_SI1_SI2_1013, ... }

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 44 of file KVIDZAFromZGridMassHomogenizer.h.

#include <KVIDZAFromZGridMassHomogenizer.h>

Inheritance diagram for KVIDZAFromZGridMassHomogenizer:

[legend]

Public Member Functions
	KVIDZAFromZGridMassHomogenizer (const TString &gridfile, const TString &datafile)

void	Analyse (double min_yield)

void	Clear ()

void	CorrectMasses (int Z, const std::string &alist, int deltaA, const KVNameValueList &mass_fit_parameters)

void	CorrectMassesForTelescope (int Z, const TString &idtel, int deltaA, const KVNameValueList &mass_fit_parameters)

KVNumberList	GetCurrentRunList () const

KVIDZAFromZGrid *	GetGrid (const TString &gridname)

void	Print () const

void	SaveGrids (const TString &filename)
	Save all grids: overwrites original file if no filename given (filename=""), otherwise writes a new file. More...

void	SetCurrentRunList (const TString &runs)

void	ShowFits (int Z, const std::string &alist)

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)

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

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

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 ()

Constructor & Destructor Documentation

◆ KVIDZAFromZGridMassHomogenizer()

KVIDZAFromZGridMassHomogenizer::KVIDZAFromZGridMassHomogenizer	(	const TString &	gridfile,
		const TString &	datafile
	)

Prepare mass homogeneity analysis

Parameters

[in]	gridfile	name/path to file containing grids to analyse
[in]	datafile	name/path to ROOT file containing 2D identification maps corresponding to the grids

Note: We assume that the histogram and grid for each telescope have the same name, which allows to associate them easily!

If more than 1 'runlist' parameter is used for the grids in the file, we print the list of runlists and inform that 1 runlist must be chosen before proceeding.

Definition at line 53 of file KVIDZAFromZGridMassHomogenizer.cpp.

Member Function Documentation

◆ Analyse()

void KVIDZAFromZGridMassHomogenizer::Analyse ( double min_yield )

inline

Parameters

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

Analyse fits for all grids according to the provide minimum (relative) yield cut-off

Definition at line 96 of file KVIDZAFromZGridMassHomogenizer.h.

◆ Clear()

void KVIDZAFromZGridMassHomogenizer::Clear ( )

inline

Reset statistics on isotope masses. After calling Clear(), Analyse() can be called with a different minimum yield threshold.

Definition at line 131 of file KVIDZAFromZGridMassHomogenizer.h.

◆ CorrectMasses()

void KVIDZAFromZGridMassHomogenizer::CorrectMasses	(	int	Z,
		const std::string &	alist,
		int	deltaA,
		const KVNameValueList &	mass_fit_parameters
	)

Apply the given mass-shift to the isotopes for element $Z$ in all grids corresponding to the given 'mass list' (see Print()).

Parameters

[in]	Z	element atomic number
[in]	alist	'mass list' attributed to grids to modify (see Print())
[in]	deltaA	mass modification $\Delta A$ to apply to all isotopes of element, can be positive or negative
[in]	mass_fit_parameters	parameters to control fit; by default KVMultiGaussIsotopeFit::mass_fit_parameters

Definition at line 268 of file KVIDZAFromZGridMassHomogenizer.cpp.

◆ CorrectMassesForTelescope()

void KVIDZAFromZGridMassHomogenizer::CorrectMassesForTelescope	(	int	Z,
		const TString &	idtel,
		int	deltaA,
		const KVNameValueList &	mass_fit_parameters
	)

Apply the given mass-shift to the isotopes for element $Z$ of the given telescope/grid

Parameters

[in]	Z	element atomic number
[in]	idtel	name of grid/telescope to modify
[in]	deltaA	mass modification $\Delta A$ to apply to all isotopes of element, can be positive or negative
[in]	mass_fit_parameters	parameters to control fit; by default KVMultiGaussIsotopeFit::mass_fit_parameters

Definition at line 297 of file KVIDZAFromZGridMassHomogenizer.cpp.

◆ GetCurrentRunList()

KVNumberList KVIDZAFromZGridMassHomogenizer::GetCurrentRunList ( ) const

inline

Definition at line 92 of file KVIDZAFromZGridMassHomogenizer.h.

◆ GetGrid()

KVIDZAFromZGrid* KVIDZAFromZGridMassHomogenizer::GetGrid ( const TString & gridname )

inline

If multiple runlists are used in the grid file, we only return a valid pointer if:

the currentRunList has been set
the grid is valid for all runs of the current runlist

Definition at line 118 of file KVIDZAFromZGridMassHomogenizer.h.

◆ Print()

void KVIDZAFromZGridMassHomogenizer::Print ( ) const

Print in the terminal the mass distribution statistics obtained from the analysis

For each element $Z$ for which a fit exists in at least 1 grid, this prints out each associated 'mass list' (list of $A$-values, KVNumberList format) and a (partial) list of the telescopes/grids which are concerned

Definition at line 150 of file KVIDZAFromZGridMassHomogenizer.cpp.

◆ SaveGrids()

void KVIDZAFromZGridMassHomogenizer::SaveGrids ( const TString & filename )

Save all grids: overwrites original file if no filename given (filename=""), otherwise writes a new file.

Definition at line 314 of file KVIDZAFromZGridMassHomogenizer.cpp.

◆ SetCurrentRunList()

void KVIDZAFromZGridMassHomogenizer::SetCurrentRunList ( const TString & runs )

inline

Definition at line 88 of file KVIDZAFromZGridMassHomogenizer.h.

◆ ShowFits()

void KVIDZAFromZGridMassHomogenizer::ShowFits	(	int	Z,
		const std::string &	alist
	)

For a given element and 'mass list' (see Print()), display the associated multi-gaussian fits for different grids/telescopes (up to 20 can be displayed).

Definition at line 187 of file KVIDZAFromZGridMassHomogenizer.cpp.

Detailed Description

Public Member Functions

Additional Inherited Members

Constructor & Destructor Documentation

◆ KVIDZAFromZGridMassHomogenizer()

Member Function Documentation

◆ Analyse()

◆ Clear()

◆ CorrectMasses()

◆ CorrectMassesForTelescope()

◆ GetCurrentRunList()

◆ GetGrid()

◆ Print()

◆ SaveGrids()

◆ SetCurrentRunList()

◆ ShowFits()