RE: LET from USRYIELD: units

From: Reddell, Brandon D (brandon.d.reddell@boeing.com)
Date: Tue May 23 2006 - 17:09:33 CEST

Next message: Dipanwita Dutta: "RE: LET from USRYIELD: units"

Previous message: Chancellor, Jeffery C. \(JSC-SF\)[UH]: "RE: LET from USRYIELD: units"
Maybe in reply to: Reddell, Brandon D: "LET from USRYIELD: units"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This is a multi-part message in MIME format.

------_=_NextPart_001_01C67E7A.E6002514
Content-Type: text/plain;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

Jeff,
=20
Thanks for the lengthy description. Sounds a bit more involved then just
using the rdyld program.
=20
Yes... I was looking at a GCR iron nuclei on a silicon target behind
aluminum shielding (the reason for going up to 50 GeV). Ultimately, I
want to look at the LET spectra from this run -- flux of particles vs.
LET . Based on the current state of FLUKA (with the low energy neutron
kerma issues) it doesn't sound like this can be done then (if the LET
curve is cutoff)???
=20
If this task can be done, does usysuw.f produce this type of output?
=20
Brandon

________________________________

From: Chancellor, Jeffery C. (JSC-SF)[UH]
[mailto:jeffery.c.chancellor@nasa.gov]=20
Sent: Monday, May 22, 2006 8:08 PM
To: Reddell, Brandon D; fluka-discuss@fluka.org
Subject: RE: LET from USRYIELD: units

Hi Brandon,

Looks like you are putting my input card to use. Let me clarify a few
things for you.

USRYIELD -123.0 201.0 -51.0 3.0 4.0 1.0LET

The first what(1) indicates which two quantities you wish FLUKA to
score. I assume here that the first quantity you want to score is LET
and you want it binned logarithmically. The second quantity is the
particles kinetic energy, which only has one interval. This input does
not request FLUKA to score a double differential yield (d2N/dLETdE).
This is done on the next line;

USRYIELD 1000.0 0.0 100.0 50.0 0.0 1403.0 &

Here you have asked FLUKA to bin the LET in 100 logarithmic bins from 0
to 1000 GeV/cm, and score the particles kinetic energy from 0 to 50 GeV
(which seems awfully high). The last what(6), 1403.0, is where you
indicate what type of output you are requesting. If you are requesting
d2N/dLETdE, then I assume here that the FLUKA material id number is 14.
Which would give for what(6) 14*100+3=3D1403. Right?

However, looking at your output, it does not look like you have LET and
energy. The first two columns look like the first 4 of the 100 LET bins,
I verified this. The next output looks like it is binning (one bin) some
type of angular distribution. Which is not what you indicated was the
desired output. Double check your calculation of what(6) and try using
the program usysuw.f located in the flutil directory. This will take
your binary file and give you two outputs. One is a summary of your run
for each detector. Part of it would look like this;

Detector n: 1 LET+0.1 =20
     (User norm: 1.
      sigma: 1. mb
      distr. scored: 201 ,
      from reg. 5 to reg. 6 )
      logar. 1st variable binning from 1.0000E-03 to 1.0000E+03 100
bins (ratio : 1.1482E+00)
      2nd variable ranges from 0.0000E+00 to 1.5000E+00
      1st variable is: Particle LET =20
      2nd variable is: Laboratory Kinetic Energy=20

The second output will be in a list format, which is much easier to
manipulate.

Finally, you can easily check your output by running a simulation
without a source card. Use protons as your primaries and use water as
your target material. Pick a more reasonable energy, say 100 MeV. This
should give you an output with the majority of the scoring (if not all)
centered around the 7-10 MeV/cm bin. Remember the output in normalized
per event.

One final note. There is a known problem with the scoring of the low
energy neutron kerma, especially at depth where most of the protons have
ranged out. From what I understand nuclear interaction with A>1 are not
being taking into account when scoring LET. You will see this if you
plot your distribution, there will be a noticeable cutoff of the higher
values of LET. Alfredo is aware of the problem. I do not know the status
of the solution.

Regards,

Jeff Chancellor

Space Radiation Analysis Group, Johnson Space Center
University of Houston, Department of Physics
jeffery.c.chancellor@nasa.gov
jchancellor@uh.edu

-----Original Message-----
From: owner-fluka-discuss@fisica.unimi.it on behalf of Reddell, Brandon
D
Sent: Mon 5/22/2006 8:20 AM
To: fluka-discuss@fluka.org
Subject: LET from USRYIELD: units

Dear FLUKA Users,

Below, I have the following syntax in my .inp file to get LET using
USRYIELD.

USRYIELD -123.0 201.0 -51.0 3.0 4.0 1.0
LET
USRYIELD 1000.0 0.0 100.0 50.0 0.0 1403.0 &

If one uses the rdyld program from the manual to read the binary output
file, the data is presented as follows:

   Yield n. 1 "LET " , generalised particle n. 201, from region n.
3 to region n. 4
      user normalisation: 1.00000000E+00, adopted cross section (if
any): 1.00000000E+00 mb
      logar. energy binning from 1.0000E-03 to 1.0000E+03 GeV, 100
bins (ratio : 1.1482E+00)
      One angular bin from 0.0000E+00 to 5.0000E+01 rad , ( 5.0000E+01
sr wide )

Differential Integral
Yield Yield
Lower energy Upper energy sr**-1 GeV**-1 number of particles

1.0000E-03 1.1482E-03 0.0000000E+00 0.0000000E+00
1.1482E-03 1.3183E-03 0.0000000E+00 0.0000000E+00
1.3183E-03 1.5136E-03 0.0000000E+00 0.0000000E+00
1.5136E-03 1.7378E-03 0.0000000E+00 0.0000000E+00
.
.
.

I thought by using code -123 I am getting d2N/d(LET)/dE. Are the first
two columns energy or LET, is the 3rd column the flux of particles in
that energy range or LET range? If energy range where is the LET
information contained?

To sum it all up: can someone tell me exactly how to read LET from this
output file (and clarify the units)?

Thanks,

Brandon

------_=_NextPart_001_01C67E7A.E6002514
Content-Type: text/html;
charset="us-ascii"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD><TITLE>RE: LET from USRYIELD: units</TITLE>
<META http-equiv=3DContent-Type content=3D"text/html; =
charset=3Dus-ascii">
<META content=3D"MSHTML 6.00.2900.2873" name=3DGENERATOR></HEAD>
<BODY>
<DIV dir=3Dltr align=3Dleft>Jeff,</DIV>
<DIV dir=3Dltr align=3Dleft> </DIV>
<DIV dir=3Dltr align=3Dleft>Thanks for the lengthy description. Sounds a =
bit more=20
involved then just using the rdyld program.</DIV>
<DIV dir=3Dltr align=3Dleft> </DIV>
<DIV dir=3Dltr align=3Dleft>Yes... I was looking at a GCR iron nuclei on a =
silicon=20
target behind aluminum shielding (the reason for going up to 50 =
GeV).=20
Ultimately, I want to look at the LET spectra from this run -- flux of =
particles=20
vs. LET . Based on the current state of FLUKA (with the low energy =
neutron kerma=20
issues) it doesn't sound like this can be done then (if the LET curve is =

cutoff)???</DIV>
<DIV dir=3Dltr align=3Dleft> </DIV>
<DIV dir=3Dltr align=3Dleft>If this task can be done, does usysuw.f =
produce this=20
type of output?</DIV>
<DIV dir=3Dltr align=3Dleft> </DIV>
<DIV dir=3Dltr align=3Dleft>Brandon</DIV> 
<DIV class=3DOutlookMessageHeader lang=3Den-us dir=3Dltr align=3Dleft>
<HR tabIndex=3D-1>
From: Chancellor, Jeffery C. =
(JSC-SF)[UH]=20
[mailto:jeffery.c.chancellor@nasa.gov] Sent: Monday, May 22, =
2006=20
8:08 PM To: Reddell, Brandon D;=20
fluka-discuss@fluka.org Subject: RE: LET from USRYIELD:=20
units </DIV>
<DIV></DIV>
Hi Brandon, Looks like you are putting my input =
card to=20
use. Let me clarify a few things for=20
you. USRYIELD     =
-123.0    =20
201.0     -51.0      =20
3.0       =
4.0      =20
1.0LET The first what(1) indicates which two quantities you wish =
FLUKA to=20
score. I assume here that the first quantity you want to score is LET =
and you=20
want it binned logarithmically. The second quantity is the particles =
kinetic=20
energy, which only has one interval. This input does not request FLUKA =
to score=20
a double differential yield (d2N/dLETdE). This is done on the next=20
line; USRYIELD    =20
1000.0       0.0    =20
100.0      =
50.0      =20
0.0    1403.0 & Here you have asked FLUKA to =
bin the=20
LET in 100 logarithmic bins from 0 to 1000 GeV/cm, and score the =
particles=20
kinetic energy from 0 to 50 GeV (which seems awfully high). The last =
what(6),=20
1403.0, is where you indicate what type of output you are requesting. If =
you are=20
requesting d2N/dLETdE, then I assume here that the FLUKA material id =
number is=20
14. Which would give for what(6) 14*100+3=3D1403. Right? However, =
looking=20
at your output, it does not look like you have LET and energy. The first =
two=20
columns look like the first 4 of the 100 LET bins, I verified this. The =
next=20
output looks like it is binning (one bin) some type of angular =
distribution.=20
Which is not what you indicated was the desired output. Double check =
your=20
calculation of what(6) and try using the program usysuw.f located in the =
flutil=20
directory. This will take your binary file and give you two outputs. One =
is a=20
summary of your run for each detector. Part of it would look like=20
this; Detector n:  1 =20
LET+0.1        (User=20
norm:        =20
1.      =20
sigma:           &=
nbsp;=20
1. mb       distr. scored: 201  =20
,       from reg. 5 to reg.  6=20
)       logar. 1st variable binning =
from =20
1.0000E-03 to  1.0000E+03   100 bins (ratio :=20
1.1482E+00)       2nd variable ranges=20
from   0.0000E+00 to =20
1.5000E+00       1st variable is: Particle=20
LET           &nbs=
p;          =20
2nd variable is: Laboratory Kinetic Energy  The second =
output will=20
be in a list format, which is much easier to manipulate. Finally, =
you can=20
easily check your output by running a simulation without a source card. =
Use=20
protons as your primaries and use water as your target material. Pick a =
more=20
reasonable energy, say 100 MeV. This should give you an output with the =
majority=20
of the scoring (if not all) centered around the 7-10 MeV/cm bin. =
Remember the=20
output in normalized per event. One final note. There is a known =
problem=20
with the scoring of the low energy neutron kerma, especially at depth =
where most=20
of the protons have ranged out. From what I understand nuclear =
interaction with=20
A>1 are not being taking into account when scoring LET. You will see =
this if=20
you plot your distribution, there will be a noticeable cutoff of the =
higher=20
values of LET. Alfredo is aware of the problem. I do not know the status =
of the=20
solution. Regards, Jeff Chancellor Space Radiation =
Analysis=20
Group, Johnson Space Center University of Houston, Department of=20
Physics jeffery.c.chancellor@nasa.gov jchancellor@uh.edu <BR=
>-----Original=20
Message----- From: owner-fluka-discuss@fisica.unimi.it on behalf of =
Reddell,=20
Brandon D Sent: Mon 5/22/2006 8:20 AM To:=20
fluka-discuss@fluka.org Subject: LET from USRYIELD: units Dear =
FLUKA=20
Users, Below, I have the following syntax in my .inp file to get =
LET=20
using USRYIELD. USRYIELD    =20
-123.0     201.0    =20
-51.0      =20
3.0       =
4.0      =20
1.0 LET USRYIELD    =20
1000.0       0.0    =20
100.0      =
50.0      =20
0.0    1403.0 & If one uses the rdyld program =
from the=20
manual to read the binary output file, the data is presented as=20
follows:    Yield n.   1  "LET  " , =

generalised particle n.  201, from region n. 3 to region=20
n.      4       =
user=20
normalisation:  1.00000000E+00, adopted cross section =
(if any): =20
1.00000000E+00 mb       logar. energy =
binning=20
from  1.0000E-03 to  1.0000E+03 GeV,   100 bins =
(ratio :=20
1.1482E+00)       One angular bin from  =

0.0000E+00 to  5.0000E+01 rad , ( 5.0000E+01 sr wide=20
)           &nbs=
p;           &nbsp=
;         =20
Differential        =20
Integral           &=
nbsp;           &n=
bsp;           &nb=
sp;=20
Yield           &n=
bsp; =20
Yield Lower energy    Upper energy    =
sr**-1=20
GeV**-1   number of=20
particles  1.0000E-03     =20
1.1482E-03     =20
0.0000000E+00     =20
0.0000000E+00  1.1482E-03     =20
1.3183E-03     =20
0.0000000E+00     =20
0.0000000E+00  1.3183E-03     =20
1.5136E-03     =20
0.0000000E+00     =20
0.0000000E+00  1.5136E-03     =20
1.7378E-03     =20
0.0000000E+00     =20
0.0000000E+00 . . . I thought by using code -123 I am =
getting=20
d2N/d(LET)/dE. Are the first two columns energy or LET, is the 3rd =
column the=20
flux of particles in that energy range or LET range? If energy range =
where is=20
the LET information contained? To sum it all up: can someone =
tell me=20
exactly how to read LET from this output file (and clarify the=20
units)? Thanks, Brandon </BODY></HTML=
>

------_=_NextPart_001_01C67E7A.E6002514--

Next message: Dipanwita Dutta: "RE: LET from USRYIELD: units"

Previous message: Chancellor, Jeffery C. \(JSC-SF\)[UH]: "RE: LET from USRYIELD: units"
Maybe in reply to: Reddell, Brandon D: "LET from USRYIELD: units"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This archive was generated by hypermail 2.1.6 : Tue May 23 2006 - 17:47:36 CEST