[fluka-discuss]: question about sampling source kinetic energy

From: Shirley Li <shirleyli.phys_at_gmail.com>
Date: Mon, 28 Mar 2016 14:12:07 -0400

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 Mon Mar 28 2016 - 21:57:02 CEST