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 Fri Feb 12 2016 - 11:59:47 CET