Monochromatic photon beam source routine

Advanced Features and User Routines
1

Dear FLUKA and Flair experts,

Following the discussion on the Marinelli Beaker source, I have reconstructed a Marinelli beaker model to simulate a uniformly distributed 40K source. However, the results I obtained are not consistent with the experimental data. To address this, I would like to replace the 40K source with a photon source emitting at 1460 keV.

The issue I am encountering is that, after modifying the beam using the BEAM card, the energy deposition scored with USRBIN changes significantly, as shown in the comparison between Figure 1 and Figure 2. Should I modify any settings in the source routine or adjust the input parameters in Flair to resolve this discrepancy?

Thank you in advance,
Cheers

Matteo

Fig. 1: 40K source
Fig. 1: 40K source1049×849 52 KB

Fig. 2: 1460keV photons source
Fig. 2: 1460keV photons source994×837 91.2 KB

marinelli.f (9.8 KB)
USER.flair (3.8 KB)

3

Dear Matteo,

if you are using a source routine for a non-isotope source, the beam divergence specified on the BEAM card is notused. You need to sample isotropic beam direction manually in the routine. You can use the RACO() subroutine to do so.
You can find more information on the available sampling routines here.

Cheers,
David