|Message: Re: Energy deposited is larger than the energy of the primary gamma beam||Not Logged In (login)|
Click on the Forum title, e.g. on the "Forums by Category" page, to read a sequence of postings to the Forum and its threads all in one page. If you are only interested in one thread or the thread following a specific posting, click the thread or the posting, which takes you to a smaller page, which contains only the part you are interested in and may be easier to navigate.
Messages are "chained" if there are only replies at the first level, i.e. 1/1.html, 1/1/1.html etc. In case of "chained" messages the message number is replaced by the icon and there is no indentation.
Inline: Display the subject line only or also the text of the posting(s); for the choice "All" the "Outline" choices are switched off.
|1||0||1||no text / full text of posting|
|2||1||All||text for level 1 only / text for All postings|
Outline: Choose the depth of the posting thread, successive toggle controls provide increasing detail.
|1||2||1||2 levels / 1 level (original posting)|
|2||3||2||3 levels / 2 levels|
|3||3||All||3 levels / all levels (all postings)|
> -I am using unix version of geant4.9.6.
Can you please use geant4.9.6.p2 if you don't already use the 2nd correction patch.
> (i) The energy deposited in a NaI scintillator is larger than the energy > of the primary gamma beam. As an example the energy of the gamma beam > was set to 1MeV shooting straight into the NaI scintillator (defined in > PrimaryGenerationAction.cc) however the energy deposited occured most > often at about 1.01965MeV. I would expect the maximum energy deposited > at 1MeV. Not sure where the stray energy comes from...
I have no explanation for this. I suggest that you run one of the /examples/extended/electromagnetic/TestEmX examples that deals with tallying the total energy deposit per event.
Your code is sufficiently complex that I can't spend the time to carefully verify how you are collecting energy deposit etc. is correct.
> (ii) The number of optical photons detected at the photon detector was > ~6800 which is less than expected ~8000. The scintillation yield was set > to 40000 photons per 1MeV. The efficiency of the PMT was 20% so the > expected photon hits should be around 8000 per 1MeV.
You should verify that the number of scintillation photons created corresponds to your yield and the energy deposition. The number of detected photons can be less because you have absorption in the medium or at the walls or you have some leak.
If all else fails, you need to simplify your code to the bare minimum - just the NaI, the Al housing and the detector volume.
> -The surfaces of NaI cylinder and photon detector (metallic disc) are > defined via G4LogicalBorderSurface. Refractive index is defined for > AlHousing,
Why do you define a refractive index for the AlHousing? Define this surface as dielectric_metal and give 100% reflectivity.
> Reflectivity of the NaI > surface is set to 94%.
Bingo, this explains why you see fewer than expected photons. And although you know that each time a photon hits the wall 6% will be absorbed. What you don't know is how often, on average, your photons hit a wall before arriving at your detector.
OK, you did define the G4BorderSurface to be dielectric_metal.
I agree with everything else you explained in your posting.
|Inline Depth:||Outline Depth:||Add message:|