Detailed Description

Silicon PHD described by Moulton formula.

Created by KVClassFactory on Mon Feb 19 17:32:55 2007 Author: franklan

The static TF1 "fMoulton" calculates the PHD of a particle of incident energy E and charge Z calculated according to Moulton's formula:

log_10(PHD) = b(Z) + a(Z)*log_10(Ed)  for Z > 2,<br>
PHD = 0  for Z <= Zmin<br>
 with  a(Z) = a_1*(Z**2/1000) + a_2<br>
         b(Z) = b_1*(100/Z) + b_2<br>
           Ed = energy lost by particle in detector (=E if particle stops)

The five parameters (a_1, a_2, b_1, b_2) can be set using method
SetParameters(a_1,a_2,b_1,b_2). It should be noted that all the PHD of all INDRA silicon detectors are calibrated with the fixed values of a_1=0.0223 and a_2=0.5682.

Definition at line 39 of file KVPulseHeightDefect.h.

#include <KVPulseHeightDefect.h>

Inheritance diagram for KVPulseHeightDefect:

[legend]

Public Member Functions
	KVPulseHeightDefect ()
	Default constructor. More...

	KVPulseHeightDefect (KVDetector *)
	Associate PHD calculation to detector. More...

virtual	~KVPulseHeightDefect ()
	Destructor. More...

virtual Double_t	Compute (Double_t, const KVNameValueList &="") const

Double_t	ELossActive (Double_t x, Double_t par)

TF1 *	GetELossFunction (Int_t Z, Int_t A, Bool_t Wrong=kFALSE)

TF1 *	GetMoultonPHDFunction (Int_t Z)
	Create TF1* fMoulton if not already done. More...

Int_t	GetZ () const

void	init ()
	default initialisations More...

virtual Double_t	Invert (Double_t, const KVNameValueList &="") const

Double_t	PHDMoulton (Double_t x, Double_t par)

void	SetZ (Int_t z)

Public Member Functions inherited from KVCalibrator
	KVCalibrator ()

	KVCalibrator (const TString &formula, const TString &type)

virtual	~KVCalibrator ()

KVDetector *	GetDetector () const

TGraph *	GetGraphOfCalibration (int npts, KVValueRange< double > input_range, const KVNameValueList &par="")

TString	GetInputSignalType () const

Int_t	GetNumberParams () const

TString	GetOutputSignalType () const

Double_t	GetParameter (int i) const

Bool_t	GetStatus () const

Bool_t	InversionFailure () const

virtual Bool_t	IsAvailableFor (const KVNameValueList &) const

Bool_t	IsUseInverseFunction () const

Double_t	operator() (Double_t x, const KVNameValueList &par="")

virtual void	Print (Option_t *opt="") const
	Print a description of the calibration object, including a list of its parameters. More...

virtual void	Reset ()

void	SetDetector (KVDetector *d)

void	SetInputSignalType (const TString &type)

virtual void	SetOptions (const KVNameValueList &)

void	SetOutputSignalType (const TString &type)

void	SetParameter (int i, Double_t par_val) const

void	SetStatus (Bool_t ready)

void	SetUseInverseFunction (Bool_t yes=kTRUE)

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...

virtual	~ KVBase ()

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

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

const Char_t *	GetLabel () const

UInt_t	GetNumber () const

UInt_t	GetNumberOfObjects () const

virtual TObject *	GetObject () const

virtual const Char_t *	GetType () const

Bool_t	HasLabel () const

virtual Bool_t	IsCalled (const Char_t *name) const

Bool_t	IsLabelled (const Char_t *l) const

virtual Bool_t	IsType (const Char_t *typ) const

virtual void	List ()

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

Double_t	ProtectedGetX (const TF1 *func, Double_t val, int &status, Double_t xmin=0.0, Double_t xmax=0.0) const

void	SetLabel (const Char_t *lab)

virtual void	SetNumber (UInt_t num)

virtual void	SetType (const Char_t *str)

Private Attributes
TF1 *	fDeltaEphd
	deltaE calculated including PHD More...

TF1 *	fMoulton
	Moulton formula for PHD = f(E,Z) More...

Int_t	fZ
	Z of nucleus to be calibrated. More...

Bool_t	wrong

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

Static Public Member Functions inherited from KVCalibrator
static KVCalibrator *	MakeCalibrator (const Char_t *type)

Static Public Member Functions inherited from KVBase
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 Char_t class_name, KVString &imp_file, KVString &dec_file, const Char_t 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="")

static Bool_t	GetDataSetEnv (const Char_t dataset, const Char_t type, Bool_t defval)

static const Char_t *	GetDataSetEnv (const Char_t dataset, const Char_t type, const Char_t *defval)

static Double_t	GetDataSetEnv (const Char_t dataset, const Char_t type, Double_t 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 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	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 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

◆ KVPulseHeightDefect() [1/2]

KVPulseHeightDefect::KVPulseHeightDefect ( )

Default constructor.

Definition at line 79 of file KVPulseHeightDefect.cpp.

◆ KVPulseHeightDefect() [2/2]

KVPulseHeightDefect::KVPulseHeightDefect ( KVDetector * d )

Associate PHD calculation to detector.

Definition at line 90 of file KVPulseHeightDefect.cpp.

◆ ~KVPulseHeightDefect()

KVPulseHeightDefect::~KVPulseHeightDefect ( )

virtual

Destructor.

Definition at line 102 of file KVPulseHeightDefect.cpp.

Member Function Documentation

◆ Compute()

Double_t KVPulseHeightDefect::Compute	(	Double_t	energy,
		const KVNameValueList &	= `""`
	)		const

virtual

Calculate the pulse height defect (in MeV) for a given energy loss in the detector and a given Z (the Z of the particle should be set first using method SetZ) The Moulton formula is used:

log_10(PHD) = b(Z) + a(Z)*log_10(E)

PHD = 0, for Z <= Zmin

with a(Z) = a_1*(Z**2/1000) + a_2 b(Z) = b_1*(100/Z) + b_2 E = energy loss of particle

Reimplemented from KVCalibrator.

Definition at line 125 of file KVPulseHeightDefect.cpp.

◆ ELossActive()

Double_t KVPulseHeightDefect::ELossActive	(	Double_t *	x,
		Double_t *	par
	)

Calculate energy lost in active layer of detector minus the Moulton PHD x[0] = incident energy par[0] = Z par[1] = A

Definition at line 185 of file KVPulseHeightDefect.cpp.

◆ GetELossFunction()

TF1 * KVPulseHeightDefect::GetELossFunction	(	Int_t	Z,
		Int_t	A,
		Bool_t	Wrong = `kFALSE`
	)

Return pointer to TF1 giving energy loss in active layer of detector minus the pulse height defect for a given nucleus (Z,A).

If Wrong=kTRUE (default:kFALSE) this will be calculated incorrectly (if the particle does not stop in the detector) by using the Moulton formula with the incident energy of the particle instead of the calculated energy loss of the particle.

Definition at line 153 of file KVPulseHeightDefect.cpp.

◆ GetMoultonPHDFunction()

TF1 * KVPulseHeightDefect::GetMoultonPHDFunction ( Int_t Z )

Create TF1* fMoulton if not already done.

Definition at line 216 of file KVPulseHeightDefect.cpp.

◆ GetZ()

Int_t KVPulseHeightDefect::GetZ ( ) const

inline

Definition at line 59 of file KVPulseHeightDefect.h.

◆ init()

void KVPulseHeightDefect::init ( void )

default initialisations

Definition at line 66 of file KVPulseHeightDefect.cpp.

◆ Invert()

Double_t KVPulseHeightDefect::Invert	(	Double_t	energy,
		const KVNameValueList &	= `""`
	)		const

virtual

Given the PHD (in MeV) of a particle of charge Z (set using SetZ() method), this method inverts the Moulton formula in order to find the energy loss of the particle in the detector.

Reimplemented from KVCalibrator.

Definition at line 237 of file KVPulseHeightDefect.cpp.

◆ PHDMoulton()

Double_t KVPulseHeightDefect::PHDMoulton	(	Double_t *	x,
		Double_t *	par
	)

Returns Moulton PHD for given E and Z:

log_10(PHD) = b(Z) + a(Z)*log_10(E), for Z > Zmin

PHD = 0, for Z <= Zmin

with a(Z) = a_1*(Z**2/1000) + a_2 b(Z) = b_1*(100/Z) + b_2 E = energy lost by particle in detector

x[0] = E (MeV) par[0] = Z

Definition at line 32 of file KVPulseHeightDefect.cpp.

◆ SetZ()

void KVPulseHeightDefect::SetZ ( Int_t z )

inline

Definition at line 55 of file KVPulseHeightDefect.h.

Member Data Documentation

◆ fDeltaEphd

TF1* KVPulseHeightDefect::fDeltaEphd

private

deltaE calculated including PHD

Definition at line 41 of file KVPulseHeightDefect.h.

◆ fMoulton

TF1* KVPulseHeightDefect::fMoulton

private

Moulton formula for PHD = f(E,Z)

Definition at line 40 of file KVPulseHeightDefect.h.

◆ fZ

Int_t KVPulseHeightDefect::fZ

private

Z of nucleus to be calibrated.

Definition at line 42 of file KVPulseHeightDefect.h.

◆ wrong

Bool_t KVPulseHeightDefect::wrong

private

Definition at line 43 of file KVPulseHeightDefect.h.

Detailed Description

Public Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ KVPulseHeightDefect() [1/2]

◆ KVPulseHeightDefect() [2/2]

◆ ~KVPulseHeightDefect()

Member Function Documentation

◆ Compute()

◆ ELossActive()

◆ GetELossFunction()

◆ GetMoultonPHDFunction()

◆ GetZ()

◆ init()

◆ Invert()

◆ PHDMoulton()

◆ SetZ()

Member Data Documentation

◆ fDeltaEphd

◆ fMoulton

◆ fZ

◆ wrong