Mikhail,
Thank you so much for your answer, it was very helpful. I just have a question about the example you directed me to:
https://indico.cern.ch/event/334606/contribution/35/attachments/653364/898409/AdvancedSources2014.pdf (pg. 38)
I don't quite understand the plot in this example. Ok, they got t, (energy in my case) and they sample it N times but what do they do with that to get the plot (MC for 3 different N-s)?
Based on how you wrote your source routine:
* Kinetic energy of the particle (GeV)
c Kinetic energy in MeV
EKIN = (EMIN1 + (EMAX1 - EMIN1) * FLRNDM(XDUMMY))
C WRITE( LUNOUT, *) 'DEBUG EKIN: ', EKIN
C Final energy in GeV
EKIN = (C10 * DEXP(-(EKIN - C20) / C30)) * EMVGEV
C WRITE( LUNOUT, *) 'DEBUG EKIN1: ', EKIN
TKEFLK (NPFLKA) = EKIN
where, i guess, your EKIN would be analogous to the t in the example. So, based on the way you wrote your code, I would say they plugged t back into f(x) function with t in place of x but in that case I don't see why MC should deviate from the theoretical y=-exp(-x) curve.
Again, my question is how they obtained MC plot.
My other, perhaps more important question is whether this is the only way to model energy distribution in fluka? For example, in MCNPX one would enter numerical values of energies and corresponding values of probabilities in source card. I saw some people did it in fluka too. Like this:
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/att-0768/source.f
from this thread:
http://www.fluka.org/web_archive/earchive/new-fluka-discuss/0768.html
Hope my questions make sense.
thanks,
Roman
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Aug 04 2015 - 03:50:10 CEST