Detailed Description

A stack of materials in which successive energy losses of charged particles can be calculated ,.

A KVMaterialStack is basically a copy of an existing detector which can be used to explore the effect of different incident angles of particles on the effective thickness of the different layers and therefore the maximum energy loss that may be observed for a given particle in the active layer of the detector.

Author: John Frankland

Date: Wed Nov 30 16:04:14 2022

Definition at line 21 of file KVMaterialStack.h.

#include <KVMaterialStack.h>

Inheritance diagram for KVMaterialStack:

[legend]

Public Member Functions
	KVMaterialStack ()=default

	KVMaterialStack (const KVDetector *, double=0.0)
	Copy all materials to stack, adjust thicknesses according to incidence angle. More...

Double_t	GetDeltaE (Int_t Z, Int_t A, Double_t Einc)
	Returns energy loss of given nucleus in the active layer of the detector. More...

Double_t	GetEIncOfMaxDeltaE (Int_t Z, Int_t A)

TF1 &	GetELossFunction (Int_t Z, Int_t A)

Double_t	GetERes (Int_t Z, Int_t A, Double_t Einc)

TF1 &	GetEResFunction (Int_t Z, Int_t A)

Double_t	GetIncidentEnergy (Int_t Z, Int_t A, Double_t delta_e, enum KVMaterial::SolType type=KVMaterial::SolType::kEmax)

Double_t	GetIncidentEnergyFromERes (Int_t Z, Int_t A, Double_t Eres)

Double_t	GetMaxDeltaE (Int_t Z, Int_t A)

TF1 &	GetMaxELossFunction (Int_t Z, Int_t A)

Double_t	GetMinimumIncidentAngleForDEMax (Int_t Z, Int_t A, Double_t Emax)

Double_t	GetTotalDeltaE (Int_t Z, Int_t A, Double_t Einc)

void	SetIncidenceAngle (double psi)

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

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

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

Private Member Functions
Double_t	ELossActive (Double_t x, Double_t par)

Double_t	EResDet (Double_t x, Double_t par)

Double_t	MaxELossActive (Double_t x, Double_t par)

Private Attributes
std::vector< KVMaterial >	fAbsorbers

Int_t	fActiveLayer = 0

const KVDetector *	fDetector = nullptr

TF1	fELossF

TF1	fEResF

Double_t	fIncAngle = 0.

Double_t	fIncAngleCosine = 1.

std::vector< double >	fLayThick
	layer thicknesses in cm More...

TF1	fmaxELossF

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 }

Public Types inherited from TObject
enum	EDeprecatedStatusBits

enum	EStatusBits

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

Constructor & Destructor Documentation

◆ KVMaterialStack() [1/2]

KVMaterialStack::KVMaterialStack ( )

default

◆ KVMaterialStack() [2/2]

KVMaterialStack::KVMaterialStack	(	const KVDetector *	D,
		double	incidence = `0.0`
	)

Copy all materials to stack, adjust thicknesses according to incidence angle.

Definition at line 9 of file KVMaterialStack.cpp.

Member Function Documentation

◆ ELossActive()

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

private

Calculates energy loss (in MeV) in active layer of stack, taking into account preceding layers

Parameters

[in]	x[0]	is incident energy in MeV
[in]	par[0]	Z of ion
[in]	par[1]	A of ion

Definition at line 78 of file KVMaterialStack.cpp.

◆ EResDet()

Double_t KVMaterialStack::EResDet	(	Double_t *	x,
		Double_t *	par
	)

private

Calculates residual energy (in MeV) of particle after traversing all layers of detector.

Warning: Returned value is -1000 if particle stops in one of the layers of the detector.

Parameters

[in]	x[0]	is incident energy in MeV
[in]	par[0]	Z of ion
[in]	par[1]	A of ion

Definition at line 140 of file KVMaterialStack.cpp.

◆ GetDeltaE()

Double_t KVMaterialStack::GetDeltaE	(	Int_t	Z,
		Int_t	A,
		Double_t	Einc
	)

Returns energy loss of given nucleus in the active layer of the detector.

Definition at line 298 of file KVMaterialStack.cpp.

◆ GetEIncOfMaxDeltaE()

Double_t KVMaterialStack::GetEIncOfMaxDeltaE	(	Int_t	Z,
		Int_t	A
	)

inline

Returns: incident energy corresponding to maximum energy loss GetMaxDeltaE() in active layer for given $(Z,A)$ and incident angle (see SetIncidenceAngle()).

Definition at line 60 of file KVMaterialStack.h.

◆ GetELossFunction()

TF1 & KVMaterialStack::GetELossFunction	(	Int_t	Z,
		Int_t	A
	)

Definition at line 38 of file KVMaterialStack.cpp.

◆ GetERes()

Double_t KVMaterialStack::GetERes	(	Int_t	Z,
		Int_t	A,
		Double_t	Einc
	)

Returns residual energy of given nucleus after the detector.

Returns 0 if Einc<=0

Definition at line 265 of file KVMaterialStack.cpp.

◆ GetEResFunction()

TF1 & KVMaterialStack::GetEResFunction	(	Int_t	Z,
		Int_t	A
	)

Definition at line 50 of file KVMaterialStack.cpp.

◆ GetIncidentEnergy()

Double_t KVMaterialStack::GetIncidentEnergy	(	Int_t	Z,
		Int_t	A,
		Double_t	delta_e,
		enum KVMaterial::SolType	type = `KVMaterial::SolType::kEmax`
	)

Returns incident energy corresponding to energy loss delta_e in active layer for a given nucleus.

By default the solution corresponding to the highest incident energy is returned This is the solution found for Einc greater than the maximum of the dE(Einc) curve. If you want the low energy solution set SolType = KVIonRangeTable::kEmin.

Warning: calculating the incident energy of a particle using only the dE in a detector is ambiguous, as in general (and especially for very heavy ions) the maximum of the dE curve occurs for Einc greater than the punch-through energy, therefore it is not always true to assume that if the particle does not stop in the detector the required solution is that for type=KVIonRangeTable::kEmax. For a range of energies between punch-through and dE_max, the required solution is still that for type=KVIonRangeTable::kEmin. If the residual energy of the particle is unknown, there is no way to know which is the correct solution.; If the given energy loss in the active layer is greater than the maximum theoretical dE for given Z & A, (dE > GetMaxDeltaE(Z,A)) then we return a NEGATIVE incident energy corresponding to the maximum, GetEIncOfMaxDeltaE(Z,A)

Definition at line 207 of file KVMaterialStack.cpp.

◆ GetIncidentEnergyFromERes()

Double_t KVMaterialStack::GetIncidentEnergyFromERes	(	Int_t	Z,
		Int_t	A,
		Double_t	Eres
	)

Calculate incident energy of nucleus from residual energy.

Warning: Returns -1 if Eres is out of defined range of values

Definition at line 169 of file KVMaterialStack.cpp.

◆ GetMaxDeltaE()

Double_t KVMaterialStack::GetMaxDeltaE	(	Int_t	Z,
		Int_t	A
	)

inline

Returns: maximum energy loss in active layer for given $(Z,A)$ and incident angle (see SetIncidenceAngle()).

See also: GetEIncOfMaxDeltaE()

Definition at line 52 of file KVMaterialStack.h.

◆ GetMaxELossFunction()

TF1 & KVMaterialStack::GetMaxELossFunction	(	Int_t	Z,
		Int_t	A
	)

Definition at line 62 of file KVMaterialStack.cpp.

◆ GetMinimumIncidentAngleForDEMax()

Double_t KVMaterialStack::GetMinimumIncidentAngleForDEMax	(	Int_t	Z,
		Int_t	A,
		Double_t	Emax
	)

inline

Returns: minimum value of incident angle $\psi$ (i.e. deviation from normal to detector) that would be required for $E_{max}$ to be the maximum energy loss possible for given $(Z,A)$

Definition at line 67 of file KVMaterialStack.h.

◆ GetTotalDeltaE()

Double_t KVMaterialStack::GetTotalDeltaE	(	Int_t	Z,
		Int_t	A,
		Double_t	Einc
	)

Returns calculated total energy loss of ion in ALL layers of the detector. This is just (Einc - GetERes(Z,A,Einc))

Definition at line 285 of file KVMaterialStack.cpp.

◆ MaxELossActive()

Double_t KVMaterialStack::MaxELossActive	(	Double_t *	x,
		Double_t *	par
	)

private

Calculates maximum energy loss (in MeV) in active layer of stack, taking into account preceding layers

Parameters

[in]	x[0]	is incident angle [deg] of particle measured wrt normal to entrance window
[in]	par[0]	Z of ion
[in]	par[1]	A of ion

Definition at line 112 of file KVMaterialStack.cpp.

◆ SetIncidenceAngle()

void KVMaterialStack::SetIncidenceAngle ( double psi )

inline

Parameters

[in] psi angle of incidence $\psi$ [deg.] of particle w.r.t. normal to detector i.e. axis along which detector has nominal thickness

Definition at line 41 of file KVMaterialStack.h.

Member Data Documentation

◆ fAbsorbers

std::vector<KVMaterial> KVMaterialStack::fAbsorbers

private

Definition at line 27 of file KVMaterialStack.h.

◆ fActiveLayer

Int_t KVMaterialStack::fActiveLayer = 0

private

Definition at line 29 of file KVMaterialStack.h.

◆ fDetector

const KVDetector* KVMaterialStack::fDetector = nullptr

private

Definition at line 26 of file KVMaterialStack.h.

◆ fELossF

TF1 KVMaterialStack::fELossF

private

Definition at line 23 of file KVMaterialStack.h.

◆ fEResF

TF1 KVMaterialStack::fEResF

private

Definition at line 24 of file KVMaterialStack.h.

◆ fIncAngle

Double_t KVMaterialStack::fIncAngle = 0.

private

Definition at line 30 of file KVMaterialStack.h.

◆ fIncAngleCosine

Double_t KVMaterialStack::fIncAngleCosine = 1.

private

Definition at line 31 of file KVMaterialStack.h.

◆ fLayThick

std::vector<double> KVMaterialStack::fLayThick

private

layer thicknesses in cm

Definition at line 28 of file KVMaterialStack.h.

◆ fmaxELossF

TF1 KVMaterialStack::fmaxELossF

private

Definition at line 25 of file KVMaterialStack.h.

Detailed Description

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ KVMaterialStack() [1/2]

◆ KVMaterialStack() [2/2]

Member Function Documentation

◆ ELossActive()

◆ EResDet()

◆ GetDeltaE()

◆ GetEIncOfMaxDeltaE()

◆ GetELossFunction()

◆ GetERes()

◆ GetEResFunction()

◆ GetIncidentEnergy()

◆ GetIncidentEnergyFromERes()

◆ GetMaxDeltaE()

◆ GetMaxELossFunction()

◆ GetMinimumIncidentAngleForDEMax()

◆ GetTotalDeltaE()

◆ MaxELossActive()

◆ SetIncidenceAngle()

Member Data Documentation

◆ fAbsorbers

◆ fActiveLayer

◆ fDetector

◆ fELossF

◆ fEResF

◆ fIncAngle

◆ fIncAngleCosine

◆ fLayThick

◆ fmaxELossF