To verify whether it has read the file properly or not, I added a write statement [write (,)] inside it so that energy values and position coordinates can be printed in the log file (see the attached tet001.txt, the first column represents energy value).
In my understanding, it should print exactly same energy values from spectrum.txt. But I found out that there is a small difference. Is it expected ?
Hi @riya,
you didn’t attach your source routine, therefore it is not evident how to compare the two files you posted.
But please consider that inside each energy bin the spectrum is assumed flat.
Cheers, luigi
Thank you for your response. Here I am attaching a sample input which I prepared as a verification. My source in this case is a spherical source whose centre and radius are written in another text file. Also I am attaching the spectrum.txt file. My questions are as follows:
From the output of log file, it seems that the energy is sampled in the form of histogram, not as a discrete spectra. Is it so?
The energy spectra I have added is for beta spectra from K - 40. Since beta is a continuous spectra, then it might be helpful if we go for histogram sampling, instead of discrete sampling. But for gamma having discrete spectra (for example Co - 60), can we use this same source.f file ?
Hi @riya,
as I anticipated in the previous email, inside each energy bin the spectrum is assumed flat.
This is correct in case of continuous distributions.
In the case of discrete spectra, you need to modify the logic in two points:
building the cumulative spectrum by adding the weight of each gamma line ENECUM(I)=ENEPRO(I)+ENECUM(I-1)
no need to sample within the bin, just use ESAMPLE = ENEPOI(I)
