Hi Shirley,
You are right, TKESUM seems to be integrating the total energy for the muons- (including the 105.1474 MeV rest mass energy). By comparison with the other cases and also from its name, I'd concur with you in that the intention probably was that it would report just the sum of kinetic energies. Anyway, I guess it is not a big deal as long as you know what the variable means.
-Mario
From: <owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it>> on behalf of Shirley Li <shirleyli.phys_at_gmail.com<mailto:shirleyli.phys_at_gmail.com>>
Date: Monday, March 28, 2016 at 11:12 AM
To: fluka-discuss <fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>>
Subject: [fluka-discuss]: question about sampling source kinetic energy
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
Received on Thu Mar 31 2016 - 21:11:26 CEST