Re: How to combine data of different FLUKA runs

From: <rhaelg_at_phys.ethz.ch>
Date: Fri, 24 Apr 2009 13:57:29 +0200

--=-OBEALYLedQlC8yLJmXI8
Content-Type: text/plain
Content-Transfer-Encoding: 7bit

Hi Vasilis

Attached you find all the files I use for this simulation.

Thanks for your help!
Greetz

Roger

On Fri, 2009-04-24 at 11:46 +0200, Vasilis Vlachoudis wrote:
> Hi Roger,
>
> could you please send the .inp and .flair file (and maybe what else is
> needed) just to check.
> This feature of flair of running in subdirectories is something new and
> maybe there is a bug somewhere.
>
> Cheers
> Vasilis
>
> Roger =3D?ISO-8859-1?Q?H=3DE4lg?=3D wrote:
> > =09 <200904071923.n37JNg5J024581_at_smtp1.mi.infn.it>
> > =09 <1239975668.17275.71.camel_at_rvext.aria>
> > =09 <Pine.LNX.4.58.0904191342160.17652_at_lxplus219.cern.ch>
> > =09 <1240502424.25525.32.camel_at_rvext.aria>
> > =09 <Pine.LNX.4.58.0904231814170.4927_at_lxplus244.cern.ch>
> > Content-Type: text/plain; charset=3D"UTF-8"
> > Message-Id: <1240504220.25525.46.camel_at_rvext.aria>
> > Mime-Version: 1.0
> > X-Mailer: Evolution 2.24.5=20
> > Content-Transfer-Encoding: 8bit
> > X-Greylist: Sender IP whitelisted, not delayed by milter-greylist-4.0 (smtp=
> > 1.mi.infn.it [192.84.138.69]); Thu, 23 Apr 2009 18:30:28 +0200 (CEST)
> > Sender: owner-fluka-discuss_at_mi.infn.it
> >
> > Hi Stefan
> >
> > The goal is to run several fluka runs in parallel, in order to use a
> > multiprocessor machine in an optimal way. The only difference of the
> > runs will be the initial random seed.=20
> >
> > To use the new possibility of recursively combine the data of different
> > runs in flair (0.7.6, "data" tab, "type" field), Vasilis recommended to
> > start the parallel runs in subfolders of the main directory. This means
> > that they can have the same name and with activated option "recursive"
> > the data from all runs would be combined in the summary files. This
> > would be an easy and clearly arranged way of handling parallel runs.=20
> >
> > So I wonder, how can I start a run in a subfolder?
> > Thanks for your help!
> >
> > Greetz
> >
> > Roger
> >
> > On Thu, 2009-04-23 at 18:15 +0200, Stefan Roesler wrote:
> >
> >> Hi Roger
> >> =20
> >> If I understand your file structure correctly then the input should be in=
> >>
> > =20
> >
> >> =20
> >> ~/flukawork/test
> >> =20
> >> and not in the subdirectory run1.
> >> =20
> >> Cheers
> >> Stefan
> >> =20
> >> =20
> >> =20
> >> On Thu, 23 Apr 2009, Roger H=C3=A4lg wrote:
> >> =20
> >>
> >>> Hi Stefan
> >>> =20
> >>> Thanks for your hints.
> >>> Yes, I did check these files. The problem is, that fort.11, fort.12,
> >>> fort.15 and fort.2 are dead links:
> >>> =20
> >>> ~/flukawork/test/fluka_10140 $ ls -lh
> >>> total 0
> >>> lrwxrwxrwx 1 rhaelg rhaelg 27 Apr 22 18:08 elasct.bin
> >>> -> /usr/local/fluka/elasct.bin
> >>> lrwxrwxrwx 1 rhaelg rhaelg 27 Apr 22 18:08 fluodt.dat
> >>> -> /usr/local/fluka/fluodt.dat
> >>> lrwxrwxrwx 1 rhaelg rhaelg 45 Apr 22 18:08 fort.1
> >>> -> /home/rhaelg/flukawork/test/ranneutron_rw3001
> >>> lrwxrwxrwx 1 rhaelg rhaelg 18 Apr 22 18:08 fort.11 -> neutron_rw3001.ou=
> >>>
> > t
> >
> >>> lrwxrwxrwx 1 rhaelg rhaelg 44 Apr 22 18:08 fort.12
> >>> -> /usr/local/fluka/libec_thihecufealw_10t.pemf
> >>> lrwxrwxrwx 1 rhaelg rhaelg 18 Apr 22 18:08 fort.15 -> neutron_rw3001.er=
> >>>
> > r
> >
> >>> lrwxrwxrwx 1 rhaelg rhaelg 17 Apr 22 18:08 fort.2 -> ranneutron_rw3002
> >>> lrwxrwxrwx 1 rhaelg rhaelg 35 Apr 22 18:08 neuxsc.bin
> >>> -> /usr/local/fluka/neuxsc-ind_260.bin
> >>> lrwxrwxrwx 1 rhaelg rhaelg 28 Apr 22 18:08 nuclear.bin
> >>> -> /usr/local/fluka/nuclear.bin
> >>> lrwxrwxrwx 1 rhaelg rhaelg 28 Apr 22 18:08 sigmapi.bin
> >>> -> /usr/local/fluka/sigmapi.bin
> >>> lrwxrwxrwx 1 rhaelg rhaelg 27 Apr 22 18:08 xnloan.dat
> >>> -> /usr/local/fluka/xnloan.dat
> >>> =20
> >>> nohup.out has the following content when just running the run in the
> >>> subfolder:
> >>> =20
> >>> $TARGET_MACHINE =3D Linux
> >>> $FLUPRO =3D /usr/local/fluka
> >>> $PEMF =3D /usr/local/fluka/libec_thihecufealw_10t.pemf
> >>> =20
> >>> Initial seed copied from /usr/local/fluka
> >>> Running fluka in /home/rhaelg/flukawork/test/fluka_10140cp: cannot stat
> >>> `run1/neutron_rw3.inp': No such file or directory
> >>> =20
> >>> =20
> >>> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D R=
> >>>
> > unning FLUKA for cycle # 1
> >
> >>> =3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
> >>> /usr/local/fluka/flutil/rfluka: line 280: neutron_rw3.inp: No such file
> >>> or directory
> >>> =20
> >>> But the stated file exists:
> >>> =20
> >>> ~/flukawork/test $ ls -lh run1/
> >>> total 8.0K
> >>> -rw-r--r-- 1 rhaelg rhaelg 7.4K Apr 22 18:08 neutron_rw3.inp
> >>> =20
> >>> I guess it is a problem with the paths. Any further suggestions?
> >>> =20
> >>> Greetz
> >>> =20
> >>> Roger
> >>> =20
> >>> On Sun, 2009-04-19 at 13:45 +0200, Stefan Roesler wrote:
> >>>
> >>>> Hi Roger,
> >>>> =20
> >>>> Did you check the content of fort.11 (the standard FLUKA output file)=
> >>>>
> > and=20
> >
> >>>> fort.15 (the error file)? Also see nohup.out which should have been=
> >>>>
> > =20
> >
> >>>> created in ~/flukawork/test and which contains the printouts from rfl=
> >>>>
> > uka.=20
> >
> >>>> In case of doubt please send me those files.
> >>>> =20
> >>>> Cheers
> >>>> Stefan
> >>>> =20
> >>>> =20
> >>>> On Fri, 17 Apr 2009, Roger H=C3=A4lg wrote:
> >>>> =20
> >>>>
> >>>>> Hi Vasilis
> >>>>> =20
> >>>>> Thank you for your answer. Your second way works fine. But I was no=
> >>>>>
> > t
> >
> >>>>> able to get your first suggestion to work.=20
> >>>>> =20
> >>>>> I can duplicate my standard run in the Run Fluka tab. I save the in=
> >>>>>
> > put
> >
> >>>>> file in an already created subfolder of the standard run. I enter a
> >>>>> Number for the random seed and select the executable.=20
> >>>>> When I start the run by clicking Run, flair tries to attach to the =
> >>>>>
> > run.
> >
> >>>>> After a few minutes the Status changes to Finished OK, without star=
> >>>>>
> > ting
> >
> >>>>> to calculate anything. The standard run in the main directory finis=
> >>>>>
> > hes
> >
> >>>>> successfully. So I guess it is not because of the input file.=20
> >>>>> =20
> >>>>> I tried to copy some eventually needed files to the subdirectory,
> >>>>> without any success.
> >>>>> =20
> >>>>> I see the following commands in the flair output:
> >>>>> =20
> >>>>> Cmd: /usr/bin/nohup /usr/local/fluka/flutil/rfluka -e myfluka -M 3 =
> >>>>>
> > neutron_rw3=20
> >
> >>>>> /usr/bin/nohup: ignoring input and appending output to `nohup.out'
> >>>>> =20
> >>>>> Cmd: /usr/bin/nohup /usr/local/fluka/flutil/rfluka -e myfluka -M 3 =
> >>>>>
> > run1/neutron_rw3=20
> >
> >>>>> /usr/bin/nohup: ignoring input and appending output to `nohup.out'
> >>>>> =20
> >>>>> A temporary directory is created for the run in the subfolder and i=
> >>>>>
> > t isn't deleted after finishing the run. It contains the following files:
> >
> >>>>> =20
> >>>>> ~/flukawork/test $ ls fluka_17313
> >>>>> elasct.bin fort.1 fort.12 fort.2 nuclear.bin xnloan.dat
> >>>>> fluodt.dat fort.11 fort.15 neuxsc.bin sigmapi.bin
> >>>>> =20
> >>>>> The manually created subdirectory only contains the input file.=20
> >>>>> =20
> >>>>> Any suggestions what I am doing wrong?
> >>>>> =20
> >>>>> Best Regards
> >>>>> =20
> >>>>> Roger
> >>>>> =20
> >>>>> On Tue, 2009-04-07 at 18:10 +0200, Vasilis Vlachoudis wrote:
> >>>>>
> >>>>>> Hi Roger,
> >>>>>> =20
> >>>>>> there are two ways in flair.
> >>>>>> =20
> >>>>>> 1. Create the runs in sub-folders. From version 0.7.6 you have th=
> >>>>>>
> > is=20
> >
> >>>>>> possibility
> >>>>>> and then in the Data Process frame you mark the "Recursive"=
> >>>>>>
> > =20
> >
> >>>>>> checkbox. Then it will
> >>>>>> scan all sub-folders for the same pattern.
> >>>>>> =20
> >>>>>> 2. In the Data frame select the master run and select also the de=
> >>>>>>
> > tector(s).
> >
> >>>>>> You click the "fitler" button right to the file display.
> >>>>>> Flair will popup a dialog to enter the pattern to search. It m=
> >>>>>>
> > ay look=20
> >
> >>>>>> a bit confusing in the begging but the most important rules are w=
> >>>>>>
> > ritten=20
> >
> >>>>>> in the syntax text box.
> >>>>>> Lets say that the auxiliary runs you created them by adding a=
> >>>>>>
> > =20
> >
> >>>>>> numerical suffix 2,3,... then modify the rule to
> >>>>>> +^\I\d\d\d\d_fort.\U$
> >>>>>> The first +/- is to include or remove a file,
> >>>>>> ^ means that the pattern should match the beginning of the string
> >>>>>> \I is substituted by flair with the filename of the input file (w=
> >>>>>>
> > ithout=20
> >
> >>>>>> the .inp)
> >>>>>> \d is any digit from 0-9 (I've added an extra one from the defaul=
> >>>>>>
> > t three=20
> >
> >>>>>> digits corresponding to the cycle added by fluka)
> >>>>>> \U is the unit number
> >>>>>> $ means the pattern should match also the end of the string
> >>>>>> =20
> >>>>>> click ok and process
> >>>>>> =20
> >>>>>> WARNING you can have different patterns for each Usrxx detector. =
> >>>>>>
> > But if=20
> >
> >>>>>> you select several then only from the pattern of the first select=
> >>>>>>
> > ed=20
> >
> >>>>>> detector will appear.
> >>>>>> =20
> >>>>>> Best Regards
> >>>>>> Vasilis
> >>>>>>
> >>>>> =20
> >>>>> =20
> >>>>> =20
> >>>>> =20
> >>>>>
> >>>> =20
> >>>>
> >>> =20
> >>> =20
> >>> =20
> >>>
> >> =20
> >>
> >
> >

--=-OBEALYLedQlC8yLJmXI8
Content-Disposition: attachment; filename="neutron_rw3.flair"
Content-Type: text/plain; name="neutron_rw3.flair"; charset="UTF-8"
Content-Transfer-Encoding: 7bit

# FLUKA Project file
Title: Photon beam distributed by 15MV spectrum impinging a RW3 phantom
Input: neutron_rw3.inp
Exec: myfluka
Debug: "Region #1" 25.0 -1.0 -25.0 30.0 1.0 30.0 100 100 1000
Source: source.f
LinkPrg: lfluka
F77bound: True
F77dline: False
FlukaVar: FLUPRO

# Run information
Run: <default>
Prev: 0
Last: 3
Status: 3
Pid: 0
StartRun: 1239974696

# USRxxx data file
Data: neutron_rw3_usrbin_neutron_eq-dose
Unit: 30
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_total_eq-dose
Unit: 31
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_energy_phantom
Unit: 32
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_energy_total
Unit: 33
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_photon_fluence_phantom
Unit: 40
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_photon_fluence_total
Unit: 41
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_electron_fluence_phantom
Unit: 42
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_electron_fluence_total
Unit: 43
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_neutron_fluence_phantom
Unit: 44
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_neutron_fluence_total
Unit: 45
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_helium4_fluence_phantom
Unit: 46
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_helium4_fluence_total
Unit: 47
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_proton_fluence_phantom
Unit: 48
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbin_proton_fluence_total
Unit: 49
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_photon_sprectrum_generated
Unit: 50
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_top_phantom_enter
Unit: 51
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_top_phantom_exit
Unit: 52
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-5_down
Unit: 53
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-5_up
Unit: 54
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-10_down
Unit: 55
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-10_up
Unit: 56
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-15_down
Unit: 57
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-15_up
Unit: 58
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-20_down
Unit: 59
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-20_up
Unit: 60
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-25_down
Unit: 61
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_z-25_up
Unit: 62
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_bottom_phantom_exit
Unit: 63
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_neutron_spectrum_bottom_phantom_enter
Unit: 64
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_photon_spectrum_top_phantom_enter
Unit: 65
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_photon_spectrum_phantom_z-15
Unit: 66
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: neutron_rw3_usrbdx_photon_spectrum_bottom_phantom_exit
Unit: 67
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData
EndRun

# Run information
Run: run1/neutron_rw3
Rnd: 82365
Exe: myfluka
Prev: 0
Last: 3
Status: 3
Pid: 0
StartRun: 1240416494

# USRxxx data file
Data: run1/neutron_rw3_usrbin_30
Unit: 30
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_31
Unit: 31
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_32
Unit: 32
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_33
Unit: 33
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_40
Unit: 40
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_41
Unit: 41
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_42
Unit: 42
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_43
Unit: 43
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_44
Unit: 44
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_45
Unit: 45
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_46
Unit: 46
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_47
Unit: 47
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_48
Unit: 48
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbin_49
Unit: 49
Type: b
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_50
Unit: 50
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_51
Unit: 51
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_52
Unit: 52
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_53
Unit: 53
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_54
Unit: 54
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_55
Unit: 55
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_56
Unit: 56
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_57
Unit: 57
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_58
Unit: 58
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_59
Unit: 59
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_60
Unit: 60
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_61
Unit: 61
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_62
Unit: 62
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_63
Unit: 63
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_64
Unit: 64
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_65
Unit: 65
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_66
Unit: 66
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData

# USRxxx data file
Data: run1/neutron_rw3_usrbdx_67
Unit: 67
Type: x
Recursive: 0
Rule: +,^\I\d\d\d_fort.\U$
Rule: +,^\I\d\d\d_ftn.\U$
EndData
EndRun

# Geometry plot "Geometry Photon beam distributed by 15MV spectrum impinging a RW3 phantom"
Plot: neutron_rw3_geometry
Title: Geometry Photon beam distributed by 15MV spectrum impinging a RW3 phantom
Format:
Type: Geometry
Keys: 1
AutoAspect: 1
origin: 0 0 0
showcenter: 0
subtype: Geometry
labels: 1
coord: X-Y
boundaries: 0
vectscale: 0.1
extends: 100 100
basisu: 1 0 0
basisv: 0 1 0
nu: 200
nv: 200
EndPlot

# USRBIN plot "Neutron equivalent dose in the phantom"
Plot: neutron_eq-dose
Title: Neutron equivalent dose in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Neutron dose equivalent / pSv per primary
Commands:
set size square
zmax: -0.2
int: 0.0120769586
zmin: -30
lt.0: 1
lw.0: 1
ls.0: 0
xmin: -15
cbcolors: 30
errors: 0
log: 0
min: 3.01340536E-10
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 1.70019598E-06
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_neutron_eq-dose
det: 1
cbpalette: FLUKA
pt.0: 1
xmax: 15
EndPlot

# USRBIN plot "Total equivalent dose in the phantom"
Plot: total_eq-dose
Title: Total equivalent dose in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Dose equivalent / pSv per primary
Commands:
set size square
zmax: -0.2
int: 237.483278
zmin: -30
lt.0: 1
lw.0: 1
ls.0: 0
xmin: -15
cbcolors: 30
errors: 0
log: 0
min: 0.00012058665
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 0.0337481759
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_total_eq-dose
det: 1
cbpalette: FLUKA
pt.0: 1
xmax: 15
EndPlot

# USRBIN plot "Energy deposition in the phantom"
Plot: energy_phantom
Title: Energy deposition in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Energy / GeV/cm^3 per primary
int: 0.00138390476
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 1
min: 2.52319304E-10
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 2.13201602E-07
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_energy_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Energy deposition"
Plot: energy_total
Title: Energy deposition
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Energy / GeV/cm^3 per primary
int: 0.00138333999
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 1
min: 9.57776322E-15
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 1.47663556E-07
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_energy_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Photon fluence in the phantom"
Plot: photon_fluence_phantom
Title: Photon fluence in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Photon fluence / particles/cm^2 per primary
Commands:
set size square
int: 35.234015
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 7.94249572E-05
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 0.00427424256
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_photon_fluence_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Photon fluence"
Plot: photon_fluence_total
Title: Photon fluence
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Photon fluence / particles/cm^2 per primary
int: 145.719607
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 4.50745802E-06
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 0.00323674991
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_photon_fluence_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Electron fluence in the phantom"
Plot: electron_fluence_phantom
Title: Electron fluence in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Electron fluence / particles/cm^2 per primary
int: 0.623472785
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 8.82563143E-08
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 9.09695591E-05
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_electron_fluence_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Electron fluence"
Plot: electron_fluence_total
Title: Electron fluence
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Electron fluence / particles/cm^2 per primary
int: 0.899659276
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 2.45344389E-09
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 6.76936543E-05
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_electron_fluence_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Neutron fluence in the phantom"
Plot: neutron_fluence_phantom
Title: Neutron fluence in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Neutron fluence / particles/cm^2 per primary
Commands:
set size square
int: 0.000105111312
lt.0: 1
lw.0: 1
ls.0: 0
ymin: -15
cbcolors: 30
errors: 0
log: 0
min: 4.83709557E-12
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 1.17880612E-08
ymax: 14.75
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_neutron_fluence_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Neutron fluence"
Plot: neutron_fluence_total
Title: Neutron fluence
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Neutron fluence / particles/cm^2 per primary
int: 0.00013460031
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 1
min: 2.6458758E-14
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 8.33944025E-09
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_neutron_fluence_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Helium-4 fluence in the phantom"
Plot: helium4_fluence_phantom
Title: Helium-4 fluence in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Helium4 fluence / particles/cm^2 per primary
int: 3.66632824E-09
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 1.86110452E-14
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 2.97974275E-12
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_helium4_fluence_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Helium-4 fluence"
Plot: helium4_fluence_total
Title: Helium-4 fluence
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Helium4 fluence / particles/cm^2 per primary
int: 8.1993217E-08
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 1
min: 1.13822958E-14
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 1.02742835E-10
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_helium4_fluence_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Proton fluence in the phantom"
Plot: proton_fluence_phantom
Title: Proton fluence in the phantom
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
XLabel: z / cm
YLabel: x / cm
CBLabel: Proton fluence / particles/cm^2 per primary
int: 1.33342436E-08
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 0
min: 2.80435774E-15
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 1.97092134E-11
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_proton_fluence_phantom
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USRBIN plot "Proton fluence"
Plot: proton_fluence_total
Title: Proton fluence
Format: .eps
Type: USRBIN
Keys: 1
AutoAspect: 1
Ratio: 0.296
XLabel: z / cm
YLabel: x / cm
CBLabel: Proton fluence / particles/cm^2 per primary
int: 9.13955468E-07
lt.0: 1
lw.0: 1
ls.0: 0
cbcolors: 30
errors: 0
log: 1
min: 3.27645282E-15
xrebin: 1
axes: Auto
hist: 2D Histogram
proj: Y
swap: 0
with.0: steps
zrebin: 1
max: 4.13004936E-10
yrebin: 1
geo: -Auto-
axes.0: x1y1
ps.0: 1
datafile: neutron_rw3_usrbin_proton_fluence_total
det: 1
cbpalette: FLUKA
pt.0: 1
EndPlot

# USR-1D plot "Photon spectrum entering treatment room"
Plot: photon_spectrum_generated
Title: Photon spectrum entering treatment room
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_photon_sprectrum_generated_tab.lis
file.1: neutron_rw3_usrbdx_50_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum entering the phantom from top"
Plot: neutron_spectrum_top_phantom_enter
Title: Neutron spectrum entering the phantom from top
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_top_phantom_enter_tab.lis
file.1: neutron_rw3_usrbdx_51_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum leaving the phantom at top"
Plot: neutron_spectrum_top_phantom_exit
Title: Neutron spectrum leaving the phantom at top
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_top_phantom_exit_tab.lis
file.1: neutron_rw3_usrbdx_52_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -5 cm beam direction"
Plot: neutron_spectrum_z-5_down
Title: Neutron spectrum in the phantom at z = -5 cm beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-5_down_tab.lis
file.1: neutron_rw3_usrbdx_52_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -5 cm negative beam direction"
Plot: neutron_spectrum_z-5_up
Title: Neutron spectrum in the phantom at z = -5 cm negative beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-5_up_tab.lis
file.1: neutron_rw3_usrbdx_53_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -10 cm beam direction"
Plot: neutron_spectrum_z-10_down
Title: Neutron spectrum in the phantom at z = -10 cm beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-10_down_tab.lis
file.1: neutron_rw3_usrbdx_54_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -10 cm negative beam direction"
Plot: neutron_spectrum_z-10_up
Title: Neutron spectrum in the phantom at z = -10 cm negative beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-10_up_tab.lis
file.1: neutron_rw3_usrbdx_55_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -15 cm beam direction"
Plot: neutron_spectrum_z-15_down
Title: Neutron spectrum in the phantom at z = -15 cm beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-15_down_tab.lis
file.1: neutron_rw3_usrbdx_56_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -15 cm negative beam direction"
Plot: neutron_spectrum_z-15_up
Title: Neutron spectrum in the phantom at z = -15 cm negative beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-15_up_tab.lis
file.1: neutron_rw3_usrbdx_57_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -20 cm beam direction"
Plot: neutron_spectrum_z-20_down
Title: Neutron spectrum in the phantom at z = -20 cm beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-20_down_tab.lis
file.1: neutron_rw3_usrbdx_58_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -20 cm negative beam direction"
Plot: neutron_spectrum_z-20_up
Title: Neutron spectrum in the phantom at z = -20 cm negative beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-20_up_tab.lis
file.1: neutron_rw3_usrbdx_59_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -25 cm beam direction"
Plot: neutron_spectrum_z-25_down
Title: Neutron spectrum in the phantom at z = -25 cm beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-25_down_tab.lis
file.1: neutron_rw3_usrbdx_61_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum in the phantom at z = -25 cm negative beam direction"
Plot: neutron_spectrum_z-25_up
Title: Neutron spectrum in the phantom at z = -25 cm negative beam direction
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_z-25_up_tab.lis
file.1: neutron_rw3_usrbdx_62_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum leaving the phantom at bottom"
Plot: neutron_spectrum_bottom_phantom_exit
Title: Neutron spectrum leaving the phantom at bottom
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_bottom_phantom_exit_tab.lis
file.1: neutron_rw3_usrbdx_63_tab.lis
EndPlot

# USR-1D plot "Neutron spectrum entering the phantom from bottom"
Plot: neutron_spectrum_bottom_phantom_enter
Title: Neutron spectrum entering the phantom from bottom
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_neutron_spectrum_bottom_phantom_enter_tab.lis
file.1: neutron_rw3_usrbdx_64_tab.lis
EndPlot

# USR-1D plot "Photon spectrum entering the phantom"
Plot: photon_spectrum_top_phantom_enter
Title: Photon spectrum entering the phantom
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
axes.0: x1y1
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x.0: 2
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_photon_spectrum_top_phantom_enter_tab.lis
file.1: neutron_rw3_usrbdx_60_tab.lis
EndPlot

# USR-1D plot "Photon spectrum inside the photon at z = -15 cm"
Plot: photon_spectrum_phantom_z-15
Title: Photon spectrum inside the photon at z = -15 cm
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
Commands:
set grid
set xtics add ("0.000511" 0.000511)
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_photon_spectrum_phantom_z-15_tab.lis
file.1: neutron_rw3_usrbdx_61_tab.lis
EndPlot

# USR-1D plot "Photon spectrum leaving the phantom"
Plot: photon_spectrum_bottom_phantom_exit
Title: Photon spectrum leaving the phantom
Format: .eps
Type: USR-1D
Keys: 1
AutoAspect: 1
XLabel: Energy / GeV
YLabel: dN/dE/dO / 1/GeV/cm2/sr per primary
Commands:
set grid
set xtics add ("0.000511" 0.000511)
xlog: 0
ylog: 0
lt.0: 1
name.1: err.
lw.0: 2
ls.0: 0
det.0: 0
det.1: 0
ndetectors: 2
block.0: 0
y2log: 0
axes.0: x1y1
show.0: 1
with.0: steps
with.1: errorbars
key.0: 1
key.1: 0
x2tics: 0
x.0: 2
x.1: 0
ps.0: 1
name.0: Detector 1
y2tics: 0
pt.0: 1
x2log: 0
y.0: 2
y.1: 2
file.0: neutron_rw3_usrbdx_photon_spectrum_bottom_phantom_exit_tab.lis
file.1: neutron_rw3_usrbdx_62_tab.lis
EndPlot

--=-OBEALYLedQlC8yLJmXI8
Content-Disposition: attachment; filename="neutron_rw3.geo"
Content-Type: text/plain; name="neutron_rw3.geo"; charset="UTF-8"
Content-Transfer-Encoding: 7bit

    0 0 Rectangular phantom (30x30x30 cm^3) at SSD100 from top z-end of the treatment room
*
* Bodies
* ------
*
* General definitions:
* blackhole to include geometry
SPH BLK 0.0 0.0 0.0 10000.0
* void
RPP VAC -1000.0 1000.0 -1000.0 1000.0 -1000.0 1000.0
* Treatment room air
RPP ROOM -200.0 200.0 -200.0 200.0 -200.0 100.0
* Body for phantom
RPP PHANTOM -15.0 15.0 -15.0 15.0 -30.0 0.0
XYP WATERP01 -1.0
XYP WATERP02 -2.0
XYP WATERP03 -3.0
XYP WATERP04 -4.0
XYP WATERP05 -5.0
XYP WATERP06 -6.0
XYP WATERP07 -7.0
XYP WATERP08 -8.0
XYP WATERP09 -9.0
XYP WATERP10 -10.0
XYP WATERP11 -11.0
XYP WATERP12 -12.0
XYP WATERP13 -13.0
XYP WATERP14 -14.0
XYP WATERP15 -15.0
XYP WATERP16 -16.0
XYP WATERP17 -17.0
XYP WATERP18 -18.0
XYP WATERP19 -19.0
XYP WATERP20 -20.0
XYP WATERP21 -21.0
XYP WATERP22 -22.0
XYP WATERP23 -23.0
XYP WATERP24 -24.0
XYP WATERP25 -25.0
XYP WATERP26 -26.0
XYP WATERP27 -27.0
XYP WATERP28 -28.0
XYP WATERP29 -29.0
END
*
* Regions
* -------
*
* Blackhole
BLKHOLE 5 +BLK -VAC
* Vacuum outside treatment room
VACUUM 5 +VAC -ROOM
* Air around target
AIR 5 +ROOM -PHANTOM
* Water phantom
WPHAN01 5 +PHANTOM -WATERP01
WPHAN02 5 +PHANTOM +WATERP01 -WATERP02
WPHAN03 5 +PHANTOM +WATERP02 -WATERP03
WPHAN04 5 +PHANTOM +WATERP03 -WATERP04
WPHAN05 5 +PHANTOM +WATERP04 -WATERP05
WPHAN06 5 +PHANTOM +WATERP05 -WATERP06
WPHAN07 5 +PHANTOM +WATERP06 -WATERP07
WPHAN08 5 +PHANTOM +WATERP07 -WATERP08
WPHAN09 5 +PHANTOM +WATERP08 -WATERP09
WPHAN10 5 +PHANTOM +WATERP09 -WATERP10
WPHAN11 5 +PHANTOM +WATERP10 -WATERP11
WPHAN12 5 +PHANTOM +WATERP11 -WATERP12
WPHAN13 5 +PHANTOM +WATERP12 -WATERP13
WPHAN14 5 +PHANTOM +WATERP13 -WATERP14
WPHAN15 5 +PHANTOM +WATERP14 -WATERP15
WPHAN16 5 +PHANTOM +WATERP15 -WATERP16
WPHAN17 5 +PHANTOM +WATERP16 -WATERP17
WPHAN18 5 +PHANTOM +WATERP17 -WATERP18
WPHAN19 5 +PHANTOM +WATERP18 -WATERP19
WPHAN20 5 +PHANTOM +WATERP19 -WATERP20
WPHAN21 5 +PHANTOM +WATERP20 -WATERP21
WPHAN22 5 +PHANTOM +WATERP21 -WATERP22
WPHAN23 5 +PHANTOM +WATERP22 -WATERP23
WPHAN24 5 +PHANTOM +WATERP23 -WATERP24
WPHAN25 5 +PHANTOM +WATERP24 -WATERP25
WPHAN26 5 +PHANTOM +WATERP25 -WATERP26
WPHAN27 5 +PHANTOM +WATERP26 -WATERP27
WPHAN28 5 +PHANTOM +WATERP27 -WATERP28
WPHAN29 5 +PHANTOM +WATERP28 -WATERP29
WPHAN30 5 +PHANTOM +WATERP29
END

--=-OBEALYLedQlC8yLJmXI8
Content-Disposition: attachment; filename="neutron_rw3.inp"
Content-Type: text/plain; name="neutron_rw3.inp"; charset="UTF-8"
Content-Transfer-Encoding: 7bit

* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7...
TITLE
Photon beam distributed by 15MV spectrum impinging a RW3 phantom
*
* use names everywhere and free format for geometry
DEFAULTS PRECISIO
RANDOMIZ 1. 0.0
START 500.0 0.0
*
* beam definitions
BEAM -.015 10000.0 1.0PHOTON
BEAMPOS 0.0 0.0 100.1 0.0 0.0 NEGATIVE
SOURCE
EMFRAY 4. BLKHOLE @LASTREG 1.0
PHOTONUC 1. AIR
PHOTONUC 1. RW3
DELTARAY 0.00001 100.0 1.04 AIR NOPRINT
DELTARAY 0.00001 100.0 1.04 RW3 NOPRINT
GEOBEGIN 21. COMBNAME
neutron_rw3.geo
GEOEND
MATERIAL 0.001225 AIR
COMPOUND -.9256 NITROGEN -.2837 OXYGEN -.01572 ARGONAIR
* 226 Polystyrene
* Chemical H-C = C-H |
* Formula / \ |
* H - C C -- C - H
* |C H | \\ // |
* | 8 8|n C - C H-C-H
* | | |
* H H |
* The variation in density ranges from 0.98 to 1.075 g/cm3.
MATERIAL 1.06 Polystyr
COMPOUND 8.0 HYDROGEN 8.0 CARBON Polystyr
* 652 Titanium oxide
MATERIAL 4.23 Titanium
COMPOUND 2.0 OXYGEN 1.0 TITANIUM Titanium
MATERIAL 1.045 RW3
COMPOUND -.98 Polystyr -0.02 Titanium RW3
* Material assignments
ASSIGNMA BLCKHOLE BLKHOLE
ASSIGNMA VACUUM VACUUM
ASSIGNMA RW3 WPHAN01 WPHAN30 1.0
ASSIGNMA AIR AIR
EMFFIX AIR 0.1 RW3 0.05
LAM-BIAS 0.0 0.01 PHOTON
EMFCUT -1.0D-06 1.0D-06 0.0 AIR PROD-CUT
EMFCUT -1.0D-06 1.0D-06 0.0 RW3 PROD-CUT
USRBIN 10. DOSE-EQ -30. 15.0 15.0 0.0NDOSPHAN
USRBIN -15.0 -15.0 -30.0 150.0 150.0 150.0&
AUXSCORE USRBIN NEUTRON NDOSPHAN NDOSPHAN AMB74
USRBIN 10. DOSE-EQ -31. 15.0 15.0 0.0DOSPHAN
USRBIN -15.0 -15.0 -30.0 150.0 150.0 150.0&
AUXSCORE USRBIN ALL-PART DOSPHAN DOSPHAN AMB74
USRBIN 10. ENERGY -32. 15.0 15.0 0.0ENPHAN
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. ENERGY -33. 20.0 20.0 100.0ENTOT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBIN 10. PHOTON -40. 15.0 15.0 0.0PHOFLU
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. PHOTON -41. 20.0 20.0 100.0PHOFLUT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBIN 10. ELECTRON -42. 15.0 15.0 0.0ELEFLU
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. ELECTRON -43. 20.0 20.0 100.0ELEFLUT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBIN 10. NEUTRON -44. 15.0 15.0 0.0NEUFLU
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. NEUTRON -45. 20.0 20.0 100.0NEUFLUT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBIN 10. 4-HELIUM -46. 15.0 15.0 0.0HE4FLU
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. 4-HELIUM -47. 20.0 20.0 100.0HE4FLUT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBIN 10. PROTON -48. 15.0 15.0 0.0PROFLU
USRBIN -15.0 -15.0 -30.0 120.0 120.0 120.0&
USRBIN 10. PROTON -49. 20.0 20.0 100.0PROFLUT
USRBIN -20.0 -20.0 -35.0 80.0 80.0 270.0&
USRBDX 99. PHOTON -50. VACUUM AIR AIRSPEC
USRBDX 0.015 0.000001 600.0 &
USRBDX 99. NEUTRON -51. AIR WPHAN01 NSPEC0D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -52. WPHAN01 AIR NSPEC0U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -53. WPHAN05 WPHAN06 NSPEC5D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -54. WPHAN06 WPHAN05 NSPEC5U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -55. WPHAN10 WPHAN11 NSPEC10D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -56. WPHAN11 WPHAN10 NSPEC10U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -57. WPHAN15 WPHAN16 NSPEC15D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -58. WPHAN16 WPHAN15 NSPEC15U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -59. WPHAN20 WPHAN21 NSPEC20D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -60. WPHAN21 WPHAN20 NSPEC20U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -61. WPHAN25 WPHAN26 NSPEC25D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -62. WPHAN26 WPHAN25 NSPEC25U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -63. WPHAN30 AIR NSPEC30D
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. NEUTRON -64. AIR WPHAN30 NSPEC30U
USRBDX 0.005 0.000001 600.0 &
USRBDX 99. PHOTON -65. AIR WPHAN01 PHSPC0D
USRBDX 0.015 0.000001 600.0 &
USRBDX 99. PHOTON -66. WPHAN15 WPHAN16 PHSPC15D
USRBDX 0.015 0.000001 600.0 &
USRBDX 99. PHOTON -67. WPHAN30 AIR PHSPC30D
USRBDX 0.015 0.000001 600.0 &
STOP

--=-OBEALYLedQlC8yLJmXI8
Content-Disposition: attachment; filename="source.f"
Content-Type: text/x-fortran; name="source.f"; charset="UTF-8"
Content-Transfer-Encoding: 7bit

*$ CREATE SOURCE.FOR
*COPY SOURCE
*
*=== source ===========================================================*
*
      SUBROUTINE SOURCE ( NOMORE )

      INCLUDE '(DBLPRC)'
      INCLUDE '(DIMPAR)'
      INCLUDE '(IOUNIT)'
*
*----------------------------------------------------------------------*
* *
* Copyright (C) 1990-2006 by Alfredo Ferrari & Paola Sala *
* All Rights Reserved. *
* *
* *
* New source for FLUKA9x-FLUKA200x: *
* *
* Created on 07 january 1990 by Alfredo Ferrari & Paola Sala *
* Infn - Milan *
* *
* Last change on 03-mar-06 by Alfredo Ferrari *
* *
* This is just an example of a possible user written source routine. *
* note that the beam card still has some meaning - in the scoring the *
* maximum momentum used in deciding the binning is taken from the *
* beam momentum. Other beam card parameters are obsolete. *
* *
*----------------------------------------------------------------------*
*
      INCLUDE '(BEAMCM)'
      INCLUDE '(FHEAVY)'
      INCLUDE '(FLKSTK)'
      INCLUDE '(IOIOCM)'
      INCLUDE '(LTCLCM)'
      INCLUDE '(PAPROP)'
      INCLUDE '(SOURCM)'
      INCLUDE '(SUMCOU)'
*
      LOGICAL LFIRST
*
c defining and saving spectrum arrays
      DIMENSION ENEPOI(0:1000),ENEPRO(0:1000),ENECUM(0:1000)
      SAVE ENEPOI, ENEPRO, ENECUM
c saving spectrum dimension
      SAVE IMAX
*
      SAVE LFIRST
      DATA LFIRST / .TRUE. /
*======================================================================*
* *
* BASIC VERSION *
* *
*======================================================================*
      NOMORE = 0
* +-------------------------------------------------------------------*
* | First call initializations:
      IF ( LFIRST ) THEN
* | *** The following 3 cards are mandatory ***
         TKESUM = ZERZER
         LFIRST = .FALSE.
         LUSSRC = .TRUE.
* | *** User initialization ***
         CALL OAUXFI('spectrum.dat',LUNRDB,'OLD',IERR)
c reading spectrum
         DO I=0,1000
           READ(LUNRDB,*,END=1972) ENEPOI(I),ENEPRO(I)
           IMAX=I
         ENDDO
         STOP ' spectrum reading uncomplete!'
 1972 CONTINUE
         ENECUM(0)=ZERZER
c building cumulative spectrum
         DO I=1,IMAX
           ENECUM(I)=(ENEPRO(I-1)+ENEPRO(I))*(ENEPOI(I)-ENEPOI(I-1))
     & +ENECUM(I-1)
         ENDDO
c normalizing cumulative spectrum
         DO I=1,IMAX
           ENECUM(I)=ENECUM(I)/ENECUM(IMAX)
         ENDDO
      END IF
* |
* +-------------------------------------------------------------------*
* Push one source particle to the stack. Note that you could as well
* push many but this way we reserve a maximum amount of space in the
* stack for the secondaries to be generated
* Npflka is the stack counter: of course any time source is called it
* must be =0
      NPFLKA = NPFLKA + 1
* Wt is the weight of the particle
      WTFLK (NPFLKA) = ONEONE
      WEIPRI = WEIPRI + WTFLK (NPFLKA)
* Particle type (1=proton.....). Ijbeam is the type set by the BEAM
* card
* +-------------------------------------------------------------------*
* | (Radioactive) isotope:
      IF ( IJBEAM .EQ. -2 .AND. LRDBEA ) THEN
         IARES = IPROA
         IZRES = IPROZ
         IISRES = IPROM
         CALL STISBM ( IARES, IZRES, IISRES )
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
* |
* +-------------------------------------------------------------------*
* | Heavy ion:
      ELSE IF ( IJBEAM .EQ. -2 ) THEN
         IJHION = IPROZ * 1000 + IPROA
         IJHION = IJHION * 100 + KXHEAV
         IONID = IJHION
         CALL DCDION ( IONID )
         CALL SETION ( IONID )
         ILOFLK (NPFLKA) = IJHION
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
* |
* +-------------------------------------------------------------------*
* | Normal hadron:
      ELSE
         IONID = IJBEAM
         ILOFLK (NPFLKA) = IJBEAM
* | Flag this is prompt radiation
         LRADDC (NPFLKA) = .FALSE.
      END IF
* |
* +-------------------------------------------------------------------*
* From this point .....
* Particle generation (1 for primaries)
      LOFLK (NPFLKA) = 1
* User dependent flag:
      LOUSE (NPFLKA) = 0
* User dependent spare variables:
      DO 100 ISPR = 1, MKBMX1
         SPAREK (ISPR,NPFLKA) = ZERZER
 100 CONTINUE
* User dependent spare flags:
      DO 200 ISPR = 1, MKBMX2
         ISPARK (ISPR,NPFLKA) = 0
 200 CONTINUE
* Save the track number of the stack particle:
      ISPARK (MKBMX2,NPFLKA) = NPFLKA
      NPARMA = NPARMA + 1
      NUMPAR (NPFLKA) = NPARMA
      NEVENT (NPFLKA) = 0
      DFNEAR (NPFLKA) = +ZERZER
* ... to this point: don't change anything
* Particle age (s)
      AGESTK (NPFLKA) = +ZERZER
      AKNSHR (NPFLKA) = -TWOTWO
* Group number for "low" energy neutrons, set to 0 anyway
      IGROUP (NPFLKA) = 0
c
c sampling from the normalized cumulative spectrum
      XYZ=FLRNDM(XYZ)
      DO I=1,IMAX
        IF(XYZ.LT.ENECUM(I)) GOTO 1973
      ENDDO
      STOP ' I did a big mistake'
 1973 CONTINUE
c the sampled energy lies in the bin I (between ENEPOI(I-1) and ENEPOI(I))
c now determining the energy inside the bin I according to a linear spectrum
      XYZ=(XYZ-ENECUM(I-1))/(ENECUM(I)-ENECUM(I-1))
      IF(ENEPRO(I).EQ.ENEPRO(I-1)) THEN
c indeed, inside the bin I the spectrum is flat!
        ESAMPLE=XYZ
      ELSE
        ESAMPLE=XYZ*(ENEPRO(I)-ENEPRO(I-1))*(ENEPRO(I)+ENEPRO(I-1))
        ESAMPLE=SQRT(ESAMPLE+ENEPRO(I-1)**2)-ENEPRO(I-1)
        ESAMPLE=ESAMPLE/(ENEPRO(I)-ENEPRO(I-1))
      ENDIF
      ESAMPLE=ESAMPLE*(ENEPOI(I)-ENEPOI(I-1))+ENEPOI(I-1)
* Kinetic energy of the particle (GeV)
      TKEFLK (NPFLKA) = ESAMPLE
* Particle momentum
      PMOFLK (NPFLKA) = SQRT ( TKEFLK (NPFLKA) * ( TKEFLK (NPFLKA)
     & + TWOTWO * AM (IONID) ) )
* Cosines (tx,ty,tz)
      TXFLK (NPFLKA) = UBEAM
      TYFLK (NPFLKA) = VBEAM
      TZFLK (NPFLKA) = WBEAM
* TZFLK (NPFLKA) = SQRT ( ONEONE - TXFLK (NPFLKA)**2
* & - TYFLK (NPFLKA)**2 )
* Polarization cosines:
      TXPOL (NPFLKA) = -TWOTWO
      TYPOL (NPFLKA) = +ZERZER
      TZPOL (NPFLKA) = +ZERZER
* Particle coordinates
      XSIZE = 10.D0
      YSIZE = 10.D0
      XFLK (NPFLKA) = -XSIZE / TWOTWO + XSIZE * FLRNDM(XXX)
      YFLK (NPFLKA) = -YSIZE / TWOTWO + YSIZE * FLRNDM(XXX)
      ZFLK (NPFLKA) = ZBEAM
* Calculate the total kinetic energy of the primaries: don't change
      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
      RADDLY (NPFLKA) = ZERZER
* Here we ask for the region number of the hitting point.
* NREG (NPFLKA) = ...
* The following line makes the starting region search much more
* robust if particles are starting very close to a boundary:
      CALL GEOCRS ( TXFLK (NPFLKA), TYFLK (NPFLKA), TZFLK (NPFLKA) )
      CALL GEOREG ( XFLK (NPFLKA), YFLK (NPFLKA), ZFLK (NPFLKA),
     & NRGFLK(NPFLKA), IDISC )
* Do not change these cards:
      CALL GEOHSM ( NHSPNT (NPFLKA), 1, -11, MLATTC )
      NLATTC (NPFLKA) = MLATTC
      CMPATH (NPFLKA) = ZERZER
      CALL SOEVSV
      RETURN
*=== End of subroutine Source =========================================*
      END

--=-OBEALYLedQlC8yLJmXI8
Content-Disposition: attachment; filename="spectrum.dat"
Content-Type: text/plain; name="spectrum.dat"; charset="UTF-8"
Content-Transfer-Encoding: 7bit

0.00012550 2.29425000
0.00015530 364.14740000
0.00018510 575.07386000
0.00021490 701.02856000
0.00024470 775.61164000
0.00027450 817.64147000
0.00030430 838.40187000
0.00033410 845.01066000
0.00036390 842.13286000
0.00039370 832.91564000
0.00042350 819.52837000
0.00045330 803.48951000
0.00048310 785.87232000
0.00051290 767.43874000
0.00054270 748.72899000
0.00057250 730.12311000
0.00060230 711.88407000
0.00063210 694.18856000
0.00066190 677.14936000
0.00069170 660.83184000
0.00072150 645.26625000
0.00075130 630.45694000
0.00078110 616.38950000
0.00081090 603.03606000
0.00084070 590.35952000
0.00087050 578.31677000
0.00090030 566.86117000
0.00093010 555.94447000
0.00095990 545.51837000
0.00098970 535.53562000
0.00101950 525.95091000
0.00104930 516.72155000
0.00107910 507.80789000
0.00110890 499.17374000
0.00113870 490.78650000
0.00116850 482.61733000
0.00119830 474.64118000
0.00122810 466.83674000
0.00125790 459.18639000
0.00128770 451.67603000
0.00131750 444.29494000
0.00134730 437.03557000
0.00137710 429.89331000
0.00140690 422.86621000
0.00143670 415.95476000
0.00146650 409.16155000
0.00149630 402.49096000
0.00152610 395.94887000
0.00155590 389.54232000
0.00158570 383.27916000
0.00161550 377.16767000
0.00164530 371.21628000
0.00167510 365.43313000
0.00170490 359.82573000
0.00173470 354.40057000
0.00176450 349.16279000
0.00179430 344.11573000
0.00182410 339.26058000
0.00185390 334.59597000
0.00188370 330.11762000
0.00191350 325.81787000
0.00194330 321.68537000
0.00197310 317.70461000
0.00200290 313.85556000
0.00203270 310.11322000
0.00206250 306.44730000
0.00209230 302.82172000
0.00212210 299.19426000
0.00215190 295.51615000
0.00218170 291.73163000
0.00221150 287.77757000
0.00224130 283.58306000
0.00227110 281.05120000
0.00230090 277.06673000
0.00233070 273.17922000
0.00236050 269.38506000
0.00239030 265.68081000
0.00242010 262.06318000
0.00244990 258.52908000
0.00247970 255.07555000
0.00250950 251.69976000
0.00253930 248.39902000
0.00256910 245.17077000
0.00259890 242.01256000
0.00262870 238.92207000
0.00265850 235.89706000
0.00268830 232.93539000
0.00271810 230.03503000
0.00274790 227.19404000
0.00277770 224.41054000
0.00280750 221.68274000
0.00283730 219.00894000
0.00286710 216.38750000
0.00289690 213.81684000
0.00292670 211.29545000
0.00295650 208.82188000
0.00298630 206.39475000
0.00301610 204.01272000
0.00304590 201.67450000
0.00307570 199.37886000
0.00310550 197.12461000
0.00313530 194.91062000
0.00316510 192.73577000
0.00319490 190.59903000
0.00322470 188.49936000
0.00325450 186.43579000
0.00328430 184.40737000
0.00331410 182.41318000
0.00334390 180.45236000
0.00337370 178.52404000
0.00340350 176.62742000
0.00343330 174.76170000
0.00346310 172.92611000
0.00349290 171.11993000
0.00352270 169.34243000
0.00355250 167.59293000
0.00358230 165.87077000
0.00361210 164.17529000
0.00364190 162.50588000
0.00367170 160.86194000
0.00370150 159.24287000
0.00373130 157.64812000
0.00376110 156.07713000
0.00379090 154.52939000
0.00382070 153.00437000
0.00385050 151.50157000
0.00388030 150.02052000
0.00391010 148.56074000
0.00393990 147.12179000
0.00396970 145.70322000
0.00399950 144.30460000
0.00402930 142.92552000
0.00405910 141.56558000
0.00408890 140.22438000
0.00411870 138.90155000
0.00414850 137.59671000
0.00417830 136.30952000
0.00420810 135.03961000
0.00423790 133.78665000
0.00426770 132.55032000
0.00429750 131.33028000
0.00432730 130.12623000
0.00435710 128.93787000
0.00438690 127.76490000
0.00441670 126.60703000
0.00444650 125.46398000
0.00447630 124.33548000
0.00450610 123.22127000
0.00453590 122.12108000
0.00456570 121.03467000
0.00459550 119.96179000
0.00462530 118.90219000
0.00465510 117.85566000
0.00468490 116.82195000
0.00471470 115.80085000
0.00474450 114.79214000
0.00477430 113.79561000
0.00480410 112.81105000
0.00483390 111.83827000
0.00486370 110.87706000
0.00489350 109.92724000
0.00492330 108.98862000
0.00495310 108.06101000
0.00498290 107.14424000
0.00501270 106.23813000
0.00504250 105.34252000
0.00507230 104.45723000
0.00510210 103.58210000
0.00513190 102.71698000
0.00516170 101.86170000
0.00519150 101.01612000
0.00522130 100.18009000
0.00525110 99.35347000
0.00528090 98.53610000
0.00531070 97.72785000
0.00534050 96.92859000
0.00537030 96.13818000
0.00540010 95.35649000
0.00542990 94.58339000
0.00545970 93.81876000
0.00548950 93.06248000
0.00551930 92.31442000
0.00554910 91.57446000
0.00557890 90.84250000
0.00560870 90.11842000
0.00563850 89.40211000
0.00566830 88.69345000
0.00569810 87.99235000
0.00572790 87.29870000
0.00575770 86.61239000
0.00578750 85.93333000
0.00581730 85.26142000
0.00584710 84.59656000
0.00587690 83.93865000
0.00590670 83.28761000
0.00593650 82.64334000
0.00596630 82.00576000
0.00599610 81.37476000
0.00602590 80.75028000
0.00605570 80.13222000
0.00608550 79.52050000
0.00611530 78.91503000
0.00614510 78.31575000
0.00617490 77.72256000
0.00620470 77.13540000
0.00623450 76.55418000
0.00626430 75.97883000
0.00629410 75.40928000
0.00632390 74.84545000
0.00635370 74.28727000
0.00638350 73.73467000
0.00641330 73.18759000
0.00644310 72.64596000
0.00647290 72.10970000
0.00650270 71.57876000
0.00653250 71.05307000
0.00656230 70.53256000
0.00659210 70.01717000
0.00662190 69.50685000
0.00665170 69.00152000
0.00668150 68.50114000
0.00671130 68.00564000
0.00674110 67.51496000
0.00677090 67.02905000
0.00680070 66.54785000
0.00683050 66.07131000
0.00686030 65.59936000
0.00689010 65.13197000
0.00691990 64.66906000
0.00694970 64.21060000
0.00697950 63.75653000
0.00700930 63.30680000
0.00703910 62.86135000
0.00706890 62.42015000
0.00709870 61.98313000
0.00712850 61.55026000
0.00715830 61.12149000
0.00718810 60.69676000
0.00721790 60.27604000
0.00724770 59.85927000
0.00727750 59.44641000
0.00730730 59.03742000
0.00733710 58.63225000
0.00736690 58.23087000
0.00739670 57.83321000
0.00742650 57.43926000
0.00745630 57.04895000
0.00748610 56.66226000
0.00751590 56.27914000
0.00754570 55.89955000
0.00757550 55.52345000
0.00760530 55.15080000
0.00763510 54.78156000
0.00766490 54.41569000
0.00769470 54.05316000
0.00772450 53.69393000
0.00775430 53.33796000
0.00778410 52.98521000
0.00781390 52.63565000
0.00784370 52.28924000
0.00787350 51.94595000
0.00790330 51.60573000
0.00793310 51.26857000
0.00796290 50.93441000
0.00799270 50.60323000
0.00802250 50.27499000
0.00805230 49.94966000
0.00808210 49.62721000
0.00811190 49.30760000
0.00814170 48.99079000
0.00817150 48.67677000
0.00820130 48.36549000
0.00823110 48.05693000
0.00826090 47.75105000
0.00829070 47.44782000
0.00832050 47.14721000
0.00835030 46.84920000
0.00838010 46.55374000
0.00840990 46.26081000
0.00843970 45.97038000
0.00846950 45.68242000
0.00849930 45.39691000
0.00852910 45.11380000
0.00855890 44.83308000
0.00858870 44.55471000
0.00861850 44.27867000
0.00864830 44.00493000
0.00867810 43.73345000
0.00870790 43.46422000
0.00873770 43.19720000
0.00876750 42.93237000
0.00879730 42.66970000
0.00882710 42.40916000
0.00885690 42.15073000
0.00888670 41.89438000
0.00891650 41.64008000
0.00894630 41.38781000
0.00897610 41.13754000
0.00900590 40.88924000
0.00903570 40.64289000
0.00906550 40.39847000
0.00909530 40.15594000
0.00912510 39.91529000
0.00915490 39.67649000
0.00918470 39.43950000
0.00921450 39.20432000
0.00924430 38.97091000
0.00927410 38.73925000
0.00930390 38.50931000
0.00933370 38.28107000
0.00936350 38.05451000
0.00939330 37.82960000
0.00942310 37.60632000
0.00945290 37.38464000
0.00948270 37.16454000
0.00951250 36.94600000
0.00954230 36.72900000
0.00957210 36.51350000
0.00960190 36.29949000
0.00963170 36.08695000
0.00966150 35.87584000
0.00969130 35.66616000
0.00972110 35.45787000
0.00975090 35.25095000
0.00978070 35.04539000
0.00981050 34.84115000
0.00984030 34.63822000
0.00987010 34.43657000
0.00989990 34.23618000
0.00992970 34.03703000
0.00995950 33.83909000
0.00998930 33.64235000
0.01001910 33.44679000
0.01004890 33.25237000
0.01007870 33.05909000
0.01010850 32.86691000
0.01013830 32.67582000
0.01016810 32.48579000
0.01019790 32.29680000
0.01022770 32.10884000
0.01025750 31.92188000
0.01028730 31.73589000
0.01031710 31.55087000
0.01034690 31.36678000
0.01037670 31.18360000
0.01040650 31.00132000
0.01043630 30.81992000
0.01046610 30.63936000
0.01049590 30.45964000
0.01052570 30.28073000
0.01055550 30.10261000
0.01058530 29.92526000
0.01061510 29.74866000
0.01064490 29.57279000
0.01067470 29.39762000
0.01070450 29.22315000
0.01073430 29.04934000
0.01076410 28.87617000
0.01079390 28.70363000
0.01082370 28.53170000
0.01085350 28.36035000
0.01088330 28.18957000
0.01091310 28.01933000
0.01094290 27.84962000
0.01097270 27.68041000
0.01100250 27.51168000
0.01103230 27.34342000
0.01106210 27.17560000
0.01109190 27.00821000
0.01112170 26.84122000
0.01115150 26.67461000
0.01118130 26.50837000
0.01121110 26.34247000
0.01124090 26.17690000
0.01127070 26.01163000
0.01130050 25.84664000
0.01133030 25.68192000
0.01136010 25.51744000
0.01138990 25.35319000
0.01141970 25.18914000
0.01144950 25.02527000
0.01147930 24.86157000
0.01150910 24.69802000
0.01153890 24.53459000
0.01156870 24.37127000
0.01159850 24.20804000
0.01162830 24.04487000
0.01165810 23.88174000
0.01168790 23.71865000
0.01171770 23.55556000
0.01174750 23.39246000
0.01177730 23.22933000
0.01180710 23.06614000
0.01183690 22.90289000
0.01186670 22.73954000
0.01189650 22.57608000
0.01192630 22.41249000
0.01195610 22.24876000
0.01198590 22.08485000
0.01201570 21.92076000
0.01204550 21.75645000
0.01207530 21.59192000
0.01210510 21.42714000
0.01213490 21.26209000
0.01216470 21.09676000
0.01219450 20.93112000
0.01222430 20.76515000
0.01225410 20.59883000
0.01228390 20.43215000
0.01231370 20.26509000
0.01234350 20.09762000
0.01237330 19.92973000
0.01240310 19.76139000
0.01243290 19.59259000
0.01246270 19.42331000
0.01249250 19.25352000
0.01252230 19.08321000
0.01255210 18.91236000
0.01258190 18.74094000
0.01261170 18.56895000
0.01264150 18.39635000
0.01267130 18.22314000
0.01270110 18.04928000
0.01273090 17.87476000
0.01276070 17.69956000
0.01279050 17.52367000
0.01282030 17.34705000
0.01285010 17.16969000
0.01287990 16.99158000
0.01290970 16.81269000
0.01293950 16.63299000
0.01296930 16.45248000
0.01299910 16.27114000
0.01302890 16.08893000
0.01305870 15.90585000
0.01308850 15.72187000
0.01311830 15.53697000
0.01314810 15.35113000
0.01317790 15.16434000
0.01320770 14.97657000
0.01323750 14.78780000
0.01326730 14.59802000
0.01329710 14.40720000
0.01332690 14.21532000
0.01335670 14.02237000
0.01338650 13.82831000
0.01341630 13.63315000
0.01344610 13.43684000
0.01347590 13.23938000
0.01350570 13.04074000
0.01353550 12.84090000
0.01356530 12.63985000
0.01359510 12.43756000
0.01362490 12.23401000
0.01365470 12.02919000
0.01368450 11.82306000
0.01371430 11.61562000
0.01374410 11.40684000
0.01377390 11.19670000
0.01380370 10.98518000
0.01383350 10.77226000
0.01386330 10.55793000
0.01389310 10.34215000
0.01392290 10.12491000
0.01395270 9.90619000
0.01398250 9.68597000
0.01401230 9.46423000
0.01404210 9.24095000
0.01407190 9.01610000
0.01410170 8.78967000
0.01413150 8.56164000
0.01416130 8.33198000
0.01419110 8.10068000
0.01422090 7.86772000
0.01425070 7.63306000
0.01428050 7.39670000
0.01431030 7.15861000
0.01434010 6.91878000
0.01436990 6.67717000
0.01439970 6.43378000
0.01442950 6.18857000
0.01445930 5.94154000
0.01448910 5.69265000
0.01451890 5.44188000
0.01454870 5.18923000
0.01457850 4.93466000
0.01460830 4.67815000
0.01463810 4.41968000
0.01466790 4.15924000
0.01469770 3.89680000
0.01472750 3.63234000
0.01475730 3.36583000
0.01478710 3.09727000
0.01481690 2.82662000
0.01484670 2.55386000
0.01487650 2.27898000
0.01490630 2.00195000
0.01493610 1.72276000
0.01496590 1.44137000
0.01499570 1.15777000

--=-OBEALYLedQlC8yLJmXI8--
Received on Fri Apr 24 2009 - 18:44:47 CEST

This archive was generated by hypermail 2.2.0 : Fri Apr 24 2009 - 18:44:47 CEST