Dear FLUKA users,
I'm trying to simulate muons interacting with a block of water. I sample
the muon kinetic energy spectrum, using the SOURCE card and source.f.
After the run, I check the spectrum from Entry SODRAW in mgdraw.f. I find
that the TKESUM seems to be the sum of sampled TOTAL energy of primaries,
instead of kinetic.
I checked this block in source.f:
IF ( ILOFLK (NPFLKA) .EQ. -2 .OR. ILOFLK (NPFLKA) .GT. 100000 ) THEN
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
ELSE IF ( ILOFLK (NPFLKA) .NE. 0 ) THEN
TKESUM = TKESUM + ( TKEFLK (NPFLKA) + AMDISC (ILOFLK(NPFLKA)) ) *
WTFLK (NPFLKA)
ELSE
TKESUM = TKESUM + TKEFLK (NPFLKA) * WTFLK (NPFLKA)
END IF
I found that the middle formula is used, where TKESUM has an mass term on
top of the sampled kinetic energy. Did I make some mistakes in the
sampling, or I misunderstood the variable TKESUM?
Thanks a lot!
Shirley
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
- application/octet-stream attachment: data
Received on Mon Mar 28 2016 - 21:57:02 CEST