Hi Vasilis,
and how do I get the number of incoming particles flux? Do I need a separate scoring card or this value can be found in the output file?
thanks,
Roman
________________________________
From: Vasilis Vlachoudis <Vasilis.Vlachoudis_at_cern.ch>
Sent: Friday, February 12, 2016 1:17 AM
To: Roman Savinov; fluka-discuss_at_fluka.org
Subject: RE: Conversion factor to get dose from USRBIN output
Hi Roman,
In this case you have to multiply (not divide) by the number of incoming particles flux times the surface
of your emission sphere
Cheers
Vasilis
________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Roman Savinov [rsavinov_at_calpoly.edu]
Sent: Friday, February 12, 2016 02:05
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: Re: Conversion factor to get dose from USRBIN output
dear Vasilis,
I am a bit confused about USRBDX. I am not using this card for dose scoring.
The value that I use with units of particles/cm2/pr is "Flux" found under "Output During Transport" in the output file:
*** NASA - O'Neill (Pinsky) flux ***
*** Samp. transition energy: 0.00100000005
*** Samp. spectral index: 0.
*** Real spectral index: 2.75
*** Minimum energy: 0.00100000005 GeV ***
*** Maximum energy: 10000. GeV ***
*** Spnorm: 0.522643089 nuc/sr/cm^2/s
*** Fluxst: 4.67721367 ions/cm^2/s
*** for no shadowing ***
*** <Z_ion>: 1.15238881 ***
*** Emission radius: 50. cm ***
*** Flux: 0.000127323961 Part/pr/cmq ***
*** Geomagnetic latitude : 0. (deg) ***
*** Scoring range +/-: 90. (deg) ***
*** Geomagnetic longitude: 0. (deg) ***
*** Scoring range +/-: 180. (deg) ***
*** Sampling latitude +/-: 90. -90. (deg) ***
*** Sampling longit. +/-: 180. -180. (deg) ***
*** P: R_Earth=R_Stormer : 59.6580887 (GV) ***
*** Vertical rigidity : 0. (GV) ***
so when you say forget about USRBDX information and multiply Gy/g by flux I am confused ...
thank you,
Roman
________________________________
From: Vasilis Vlachoudis <Vasilis.Vlachoudis_at_cern.ch>
Sent: Thursday, February 11, 2016 7:29 AM
To: Roman Savinov; fluka-discuss_at_fluka.org
Subject: RE: Conversion factor to get dose from USRBIN output
Hi Roman,
the result is not correct because the units are correct :)
Please give a look on the scoring lecture from our last course:
https://indico.cern.ch/event/442634/contribution/12/attachments/1184625/1716758/07_Scoring_2015.pdf
The output of FLUKA is normalized per primary particle (since you are also asking NO-NORM) in your GCR-SPE card.
So you have a USRBIN giving Gy/pr and a USRBDX giving the particles/cm2/pr
Dividing these two quantities just cancels the primary beam information as you correctly say you get rid of it.
However then it gives a quantity of Dose vs fluence, which tells nothing about your incoming spectrum!
Forget the information from the USRBDX and simply take the dose Gy/pr and multiply with with the
incoming flux of particles then you will get Gy per the total number of particles of your event.
Cheers
Vasilis
________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Roman Savinov [rsavinov_at_calpoly.edu]
Sent: Thursday, February 11, 2016 07:00
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: Conversion factor to get dose from USRBIN output
dear Fluka experts,
I am scoring dose with USRBIN. It has output of GeV/g/pr.
Multiplying this by 1.062e-7 gives Gy/pr (manual-USRBIN-note 5).
Now, to get rid of the primary, I divide by the flux (found in the output file) to get (Gy)*(cm^2)
Dividing by the area of the detector should give me Gy but the numbers Im getting is in the order of 10^-13 which is obviously too small given the GCR source and the target (in space) is a cylinder of water 30 cm long shielded by 3.7 cm of Aluminum.
What do I do wrong?
thank you,
Roman
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Sat Feb 13 2016 - 00:35:14 CET