GLG4Scint Class Reference

GLG4Scint is an extremely modified version of the G4Scintillation process, so much so that it's not even a G4Process anymore! Features include arbitrary scintillation light time profile and spectra, Birks' law, particle-dependent specification of all parameters, and reemission of optical photons killed by other processes. More...

#include <GLG4Scint.hh>

List of all members.

Public Member Functions
	GLG4Scint (const G4String &tableName="", G4double lowerMassLimit=0.0)
G4VParticleChange *	PostPostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
G4double	GetLowerMassLimit (void) const
void	DumpInfo () const
MyPhysicsTable *	GetMyPhysicsTable (void) const
G4int	GetVerboseLevel (void) const
void	SetVerboseLevel (G4int level)
void	SetNewValue (G4UIcommand *command, G4String newValues)
G4String	GetCurrentValue (G4UIcommand *command)
Static Public Member Functions
G4VParticleChange *	GenericPostPostStepDoIt (const G4Step *pStep)
void	ResetTotEdep ()
G4double	GetTotEdep ()
G4double	GetTotEdepQuenched ()
G4bool	GetDoScintillation ()
G4ThreeVector	GetScintCentroid ()
Protected Attributes
int	verboseLevel
MyPhysicsTable *	myPhysicsTable
G4double	myLowerMassLimit
G4ParticleChange	aParticleChange
Static Protected Attributes
int	quenchtype = 0
G4std::vector< GLG4Scint * >	masterVectorOfGLG4Scint
G4UIdirectory *	GLG4ScintDir = 0
G4int	maxTracksPerStep = 180000
G4double	meanPhotonsPerSecondary = 3.0
G4bool	doScintillation = true
G4bool	doReemission = true
G4double	totEdep = 0.0
G4double	totEdep_quenched = 0.0
G4ThreeVector	scintCentroidSum = 0.0

---

Detailed Description

GLG4Scint is an extremely modified version of the G4Scintillation process, so much so that it's not even a G4Process anymore! Features include arbitrary scintillation light time profile and spectra, Birks' law, particle-dependent specification of all parameters, and reemission of optical photons killed by other processes.

• Has a GenericPostPostStepDoIt() function (note two "Post"s) instead of a PostStepDoIt() function. GenericPostPostStepDoIt() should be called by user in UserSteppingAction. This guarantees that GLG4Scint will absolutely be the last process considered, and will definitely see the energy loss by charged particles accurately.

• Modified to allow specification of absolute yield spectra, resolution scale, Birk's-law coefficient, and digitized waveform, customized for medium and (optionally) particle type.

• No longer calls G4MaterialPropertiesTable::GetProperty() in [Post]PostStepDoit() -- all needed data can be found quickly in the internal physics table.

• Uses poisson random distribution for number of photons if mean number of photons <= 12.

• The total scintillation yield is now found implicitly from the integral of the scintillation spectrum, which must now be in units of photons per photon energy.

• The above feature has been modified by Dario Motta: a scintillation yield CAN be defined and -if found- used instead of the implicit integral of the scintillation spectrum. This allows having scintillators with the same spectrum, but different light yields.

• The materials property tables used are SCINTILLATION == scintillation spectrum SCINTWAVEFORM == scintillation waveform or time constant SCINTMOD == resolution scale, Birk's constant, reference dE/dx

• SCINTILLATION is required in each scintillating medium. (Okay to omit if you don't want the medium to scintillate.)

• If SCINTWAVEFORM is missing, uses exponential waveform with default ScintillationTime. If SCINTWAVEFORM contains negative "Momentum"'s then each "Momentum" is the decay time and its corresponding value is the relative strength of that exponential decay. Otherwise, the "PhotonEnergy" of each element is a time, and the Value of each element is the relative strength.

• Default values of resolution scale (=1.0), Birk's constant (=0.0) and reference dE/dx (=0.0) are used if all or part of SCINTMOD is is missing. SCINTMOD "PhotonEnergy" values should be set to the index number (0.0, 1.0, 2.0, with no units).

• Birk's law (see 1998 Particle Data Booklet eq. 25.1) is implemented as yield(dE/dx) = yield_ref * dE/dx * (1 + kb*(dE/dx)_ref) / (1 + kb*(dE/dx)). I.e., the scintillation spectrum given in SCINTILLATION is measured for particles with dE/dx = (dE/dx)_ref. The usual formula is recovered if (dE/dx)_ref = 0.0 (the default). This is useful if you have an empirically-measured spectrum for some densely-ionizing particle (like an alpha).

• The constructor now accepts an additional string argument, tablename, which allows selection of alternate property tables. E.g, tablename = "neutron" might be used to allow specification of a different waveform for scintillation due to neutron energy deposition. The code then searches for tables with names of the form "SCINTILLATIONneutron" If it finds such a table, that table is used in preference to the default (un-suffixed) table when stepping particles of that type.

• The process generates at most maxTracksPerStep secondaries per step. If more "real" photons are needed, it increases the weight of the tracked opticalphotons. Opticalphotons are thus macro-particles in the high-scintillation case. The code preserves an integer number of real photons per macro-particle.

---

The documentation for this class was generated from the following files:

• GLG4Scint.hh
• GLG4Scint.cc

---