How to insert zenith angle info on a source.f from txt file

[horvathd (Dávid Horváth)]

Dear Georgios,

please take a look the source routine posted here:

Source routine - Reading external spectrum

Dear Hassan, you can use the following source routine to read your spectrum file: source_spectrum.f (10.5 KB) A couple things to note: You need to change the file name on line 68. If the energy unit in the spectrum file is not GeV, then enable and modify line 74 to convert it to GeV The comments (first line) in your spectrum file has to be removed The routine uses a linear interpolation between the specified energy points Cheers, David

1. You can replicate the same reading and sampling parts for the zenith angle, then modify the direction cosine variables as needed.
   The SFECFE subroutine will return a random direction in 2D.

   CALL SFECFE(SINT, COST)
   

2. Unfortunately the volumetric beam parameters are not inherited, only the point defined on the BEAMPOS card via the XBEAM, YBEAM and ZBEAM variables. If you want a rectangle as a source plane, you will have to sample it by yourself.

   You can use the FLRNMD(XDUMMY) function to get a random number between [0,1) sampled uniformly.

Cheers,
David