RE: [fluka-discuss]: FW: FLUKA

From: Joachim Vollaire <>
Date: Fri, 12 Aug 2016 07:19:28 +0000

You are missing the following cards for the activation of the decay of the “source” isotope (semi-analogue mode)


And to associate the detectors (two USRBIN) to score from decay particles

DCYSCORE -1.0 0. 0. doses doses USRBIN

DCYSCORE -1.0 0. 0. photon photon USRBIN

It should work adding those cards and you should be able to check that it is consistent with the “mono-energetic” photon source

From: James Wetzel []
Sent: 11 August 2016 21:33
To: Joachim Vollaire <>
Subject: Re: [fluka-discuss]: FW: FLUKA

Alright many thanks for your help this far -

The geometry makes sense and looks good in the viewer without any errors.

I was able to get dose values on a map with the extended source of mono energetic photons, however, when trying to set up the isotope, my plots are empty! And I think there are errors, but I’m not seeing any on the output.

I’ve attached my updated inp file - I assume I am probably initializing the isotope incorrectly.

On Aug 11, 2016, at 1:10 PM, James Wetzel <<>> wrote:

Thank you very much!

In this case, the entire radiator is exposed to the scintillator, there is no collimation, which should create a 1/r effect, compared to a 1/r^2 effect, naively, so I think I need to simulate the entire cylinder as being an activated source.

Will this still accomplish that?

On Aug 11, 2016, at 11:56 AM, Joachim Vollaire <<>> wrote:

Hi James

Looking at your input a few things.

The geometry errors:

You define your blakchole external and internal surface with the same sphere radius (use a big number for the blakchole)

SPH blkbody 0.0 0.0 0.0 50.0
* Air sphere
SPH air 0.0 0.0 0.0 50.0

To define the blackbody replace :

BLKBODY 5 +blkbody -(air+Plastic)
By (“Plastic” is already contained inside air…)
BLKBODY 5 +blkbody -air

The error problem should be solved with the above.

For scoring, you can use a region based USRBIN Scoring for DOSE (GeV/g) in the region describing your scintillator.

For accuracy of the results, did you consider simulating the collimation system of the irradiator and a point source and isotropic source inside ? I have no experience with irradiator but I am not sure how you can derive from the source activity the photon fluence out of the collimator. I guess that using only solid angle does not work due to scattering effect ?

For the source description, note that instead of using mono-energetic 667 keV photon you can use an isotope as the source and benefit from the decay data in FLUKA (of interest for more complicated cases and to get accurate branching ratio which for example is not 100 % for decay of Cs-137 in Ba-137m emitting the 661.6 keV gamma…).
To do so look at the course material related to definition of sources slide 21 in the following presentation (don’t forget to activate the decay of the particle by calling the RADDECAY card:

For normalization, the results (deposited energy) will be by isotope decay so you just need to multiply the results by the source activity ( in Bq)


From: James Wetzel []
Sent: 11 August 2016 17:26
To: Joachim Vollaire <<>>
Cc: Giuseppe Battistoni <<>>; Paola Sala <<>>
Subject: Re: FLUKA

Thank you very much!

Attached is my .inp file I am working with.

I have a cylindrical gamma radiator of cesium 137, with dimensions Diameter: 2.54 cm and Length: 30.48 cm, and 6 Ci activity.

I have a sample of plastic scintillator with dimensions 10cm x 10cm x 0.1cm, which is placed a distance 5 cm from the radiator as in the picture below, where the 10x10cm scintillator plane is parallel with the cylinder, and the center of the plastic is aligned with the center of the cylinder, as below:

Radiator: ===============


Plastic Scint: ———

What I need at the end of the day is the total dose absorbed by the plastic scintillator, in Gy or Mrad.

I am getting Geometry errors (I think), and I’m not sure how to get the dose out properly with scoring.

Thank you very much for your help!


You can manage unsubscription from this mailing list at
Received on Fri Aug 12 2016 - 11:08:40 CEST

This archive was generated by hypermail 2.3.0 : Fri Aug 12 2016 - 11:08:44 CEST