Message: G4Scintillation - optical photon generation Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question G4Scintillation - optical photon generation 

Forum: Particles
Date: 07 Jun, 2012
From: Bruno Arsenali <Bruno Arsenali>

Hi,

I have problem with number of generated optical photons from scintillation process (G4Scintillation). 

To demonstrate whats my problem lets start wit definition of scintillation crystal properties. I will make it as simple as I can, so that other scintillation parameters such as RESOLUTIONSCALE dont affect number of generated photons.

Here is definition of scintillation crystal properties:

// *** ScintillationCrystal properties
G4MaterialPropertiesTable* ScintillationCrystalProperties = new G4MaterialPropertiesTable();
const G4int N = 1;
G4double Energies[N] = {1.0*eV}; 
G4double Fastcomp[N] = {1.0}; 
G4double Rindices[N] = {1.0};
ScintillationCrystalProperties-&gt;AddProperty("FASTCOMPONENT", Energies, Fastcomp, N);
ScintillationCrystalProperties-&gt;AddProperty("RINDEX",        Energies, Rindices, N); 
ScintillationCrystalProperties-&gt;AddConstProperty("SCINTILLATIONYIELD",10./keV); 
ScintillationCrystalProperties-&gt;AddConstProperty("RESOLUTIONSCALE", 0.0);
ScintillationCrystalProperties-&gt;AddConstProperty("FASTTIMECONSTANT", 100.0*ns);
ScintillationCrystalProperties-&gt;AddConstProperty("YIELDRATIO", 1.0);


Number of generated photons follow normal distribution. Standard deviation and expected photon number of normal distribution can be calculated with:

ExpectedPhotonNumber=SCINTILLATIONYIELD*EnergyDeposit
StandardDeviation=RESOLUTIONSCALE*sqrt(ExpectedPhotonNumber)

Base on given ScintillationCrystal properties we can calculate that for energy deposit of 122keV (Gamma ray) we should get 1220 optical photons (standard deviation should be 0).
But instead of this in my simulation for energy deposit of 122keV I'm getting 332 photons (standard deviation 0).

EnergyDeposit is calculated with:
step-&gt;GetPreStepPoint()-&gt;GetTotalEnergy()-step-&gt;GetPostStepPoint()-&gt;GetTotalEnergy()

Number of generated photons is calculated correct.

Please help.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 Feedback: Re: G4Scintillation - optical photon generation   (Gumplinger Peter - 07 Jun, 2012)
2 None: Re: G4Scintillation - optical photon generation   (Bruno Arsenali - 08 Jun, 2012)
1 Feedback: Re: G4Scintillation - optical photon generation   (Gumplinger Peter - 08 Jun, 2012)
1 None: Re: G4Scintillation - optical photon generation   (Bruno Arsenali - 09 Jun, 2012)
2 None: Re: G4Scintillation - optical photon generation   (Cejo George - 15 Jun, 2012)
... 2 Message(s)
 Add Message Add Message
to: "G4Scintillation - optical photon generation"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]