|Message: Re: How to explain energy spectra of emmited photons from a scintillator, in my code?||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)|
> 1. How translate Geant4 these values in such array? Only one by one, or > not? I mean, for example, if energy of a photon be 2.035*eV (very close > to 2.034*eV), do Geant4 accept it as a value of the array? From C++, I > know the answer in NO.
The user provides the (relataive) strength of the emmission spectrum at a finite number of photon energies. You can pick those out from your continuous spectrum as best and accurate as you can. The program will then numerically integrate this array, normalize it to one, and store it in a new array. When a scintillation photon energy is to be sampled (throwing a number between 0->1) the corresponding photon energy is found by interpolating the integral array. The interpolation used to be linear but is now quadratic. This will result in a photon energy not identical to any of your array input values.
The key word here is 'Inverse Transform Sampling' of an empirical density function at discrete points.
> 2. Perhaps because matter of origin of these photons is a scintillator, > so only "some energies" will release; Is that correct?
No, the program will produce a continuous spectrum that smoothly approaches your input the more points you provide.
> 3. Please observe attached pic, about NE102A scintillator The emission > spectra is "continuous". How to explain wavelengths of emitted photons, > in my code as an array?
As I said, just pick as many points from the spectrum as you feel necessary to describe it in sufficient detail, put the energy and emission strength in an array and feed it to the program. To test how well the program reproduces your spectrum you can histogram the scintillation photon energy the program has sampled. If you are not happy with the representation you will have to increase your number of points or where you choose to pick out the points. Obviously, you'll want points at peaks and valleys of your spectrum
|Inline Depth:||Outline Depth:||Add message:|