RE: LET from USRYIELD: units

From: Chancellor, Jeffery C. \(JSC-SF\)[UH] (jeffery.c.chancellor@nasa.gov)
Date: Tue May 23 2006 - 03:07:33 CEST

Next message: Reddell, Brandon D: "RE: LET from USRYIELD: units"

Previous message: Katherine Harine: "Re: LET from USRYIELD: units"
Maybe in reply to: Reddell, Brandon D: "LET from USRYIELD: units"
Next in thread: Dipanwita Dutta: "RE: LET from USRYIELD: units"
Reply: Dipanwita Dutta: "RE: LET from USRYIELD: units"
Messages sorted by: [ date ] [ thread ] [ subject ] [ author ] [ attachment ]

This is a multi-part message in MIME format.

------_=_NextPart_001_01C67E05.45ADB613
Content-Type: text/plain;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

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?=20

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.=20
      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.=20

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.=20

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
=20
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=20
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?=20

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_01C67E05.45ADB613
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable

<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV=3D"Content-Type" CONTENT=3D"text/html; =
charset=3Diso-8859-1">
<META NAME=3D"Generator" CONTENT=3D"MS Exchange Server version =
6.5.7638.1">
<TITLE>RE: LET from USRYIELD: units</TITLE>
</HEAD>
<BODY>

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   
     (User =
norm:         1. 
      =
sigma:           &=
nbsp; 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           &nbs=
p;     
      2nd variable is: Laboratory Kinetic =
Energy  
 
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 ) 
 
            &=
nbsp;           &n=
bsp;        =
Differential         =
Integral 
            &=
nbsp;           &n=
bsp;           =
Yield           &n=
bsp;  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

</BODY>
</HTML>
------_=_NextPart_001_01C67E05.45ADB613--

Next message: Reddell, Brandon D: "RE: LET from USRYIELD: units"

Previous message: Katherine Harine: "Re: LET from USRYIELD: units"
Maybe in reply to: Reddell, Brandon D: "LET from USRYIELD: units"
Next in thread: Dipanwita Dutta: "RE: LET from USRYIELD: units"
Reply: Dipanwita Dutta: "RE: 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 - 08:08:00 CEST