From: Mr. Bhushankumar Jagnnath Patil. (bjp@physics.unipune.ernet.in)
Date: Sat Jan 12 2008 - 12:57:44 CET
Respected sir,
I am very thankful to you.
Just check this what i understand from your help.
When i am scoring in linear binning in energy and for only one bin solid
angle (one way), to get fluence (part/cm-2/pri) i will do
flux (pat/Gev/cm-2/pr) X (Max. energy of the bin - Min. energy of the bin)
X 2pi.
i.e form *_tab.lis file
IIIrd column X (IInd column - Ist column) X 2pi
for total fluence addition of the above quantity for all the bins.
is it correct?
and what if i am scoring energy in log binning, in case of neutron i am
using log binning.
the same eq i have to use ?
Also I want to know why to multiply by 2 pi, as i am using only single
binning of angular distribution.
I will be very thankful to you.
Thanking you.
Patil Bhushan J.
> Dear Bhushan
> in Flair, the option you want is Y x DX. But be careful: this will give
you the INTEGRATED fluence in each interval: in general, this
> representatio=
> n
> is not good for plotting, unless each interval has the same width (i.e.,
can be used for linear binning but is meaningless for logarithmic
binning).
>
> I would suggest that you try the utility program usxsuw (you find it in
$FLUPRO/flutil). It gives you:
> 1) the total fluence, integrated over energy and solid angle
> 2) the fluence differential in energy, but integrated over solid angle
> (from what I understand you are not interested in this)
> 3) the fluence CUMULATIVE in energy and integrated over solid angle.
> This is exactly what you are looking for: if you take the difference
between any two energies you get the total integrated fluence
between
> t=
> hose
> two energies. If the two energies are the the limits of a single
bin,
> you get the fluence in that bin.=20
> Only, remember that usxsuw works only with binary output files.
>
> Finally, you can obtain what you want by a simple multiplication:
multiply
> each bin by 2 pi (or 4 pi, depending if you have a 1-way or 2-way
scoring)
> and by the difference of energy limits of the bin. The latter is a
constant=
> ,
> since you are doing a linear scoring. This constant is equal to
> (Max. energy requested - Min. energy)/number of bins.
>
> Alberto
>
>
> On Fri, 11 Jan 2008, Mr. Bhushankumar Jagnnath Patil. wrote:
>
>> Respected sir
>> Thank you for your kind help. The information is very useful for
>> me.
>> I want to ask one more question that which option of flair i will use
to calculate the y axis result in (cm-2 per incident particle) unit.
Just i want to remove the per GeV and per Sr terms or how i will do it
manually.
>> waiting for your reply
>> Thanking you in advance.
>>> Dear Bhushan,
>>> From the manual:
>>> "Note: The results of a USRBDX boundary crossing estimator are
>>> always given as DOUBLE DIFFERENTIAL distributions of
fluence (or current) in energy and solid angle, in units of
cm-2 GeV-1 sr-1 per incident primary, EVEN WHEN ONLY 1
INTERVAL (BIN) HAS BEEN REQUESTED, which is often the case
for angular distributions."
>>> The output of the USRBDX card is fluence. Now if you do what you are
asking "multiply by the difference of the Ist and IInd column" what
you
>>> will get is a histogram of weighted counts for every energy bin. For
flair please look in the last fluka courses the flair-2.pdf file. You
have several options depending on what you will like to plot, both on
the X-axis and the Y-axis.
>>> Variables: "Low" xl =3D Ist column, "High" xh =3D IInd column, Y =3D II=
> Ird
>>> column
>>> X-axis
>>> 1. GeoMean: sqrt(xl*xh) <Xgeo>
>>> The geometrical mean of the two values. Useful if you want to plot
>>> point in the geometrical center of each bin. To be used only if you make
>>> a logX plot. Be careful this is only for plotting it doesn't care what
you have scored in fluka linear or log in energy.
>>> 2. Mean (xl+xh)/2 <X>
>>> Normal mean. To be used for linear plots
>>> 3. Low (xl)
>>> Low value of the bin. To be used when plotting like a histogram
the
>>> data, together with the "steps" style
>>> 4. High (xh)
>>> High value of the bin. Same as 3 but for "histeps"
>>> Y-axis
>>> There you have several choices on the way you want to represent your dat=
> a
>>> i) Y
>>> The value of the IIIrd column as given by fluka
>>> ii) Y x <X>
>>> The value multiplied by the mean of the first two columns. To plot in
isolethargic way with linear x-axis
>>> iii) Y x <Xgeo>
>>> Multiplied by the geometric mean. Isolethargic way with log x-axis iv)
Y x Xl
>>> ...
>>> v) Y x Xh
>>> ...
>>> vi) Y x DX
>>> Number of "counts" per bin.
>>> Examples:
>>> * To plot normal fluence like a histogram (log or linear) you have to
choose
>>> X: Low [xl]
>>> Y: Y
>>> With: steps
>>> * To plot normal fluence with error bars (linear)
>>> X: Mean (xl+xh)/2
>>> Y: Y
>>> With: errorbars
>>> * To plot normal fluence with error bars (log)
>>> X: GeoMean sqrt(xl*xh)
>>> Y: Y
>>> With: errorbars
>>> * To plot isolethargic fluence as histogram (log)
>>> X: Low [xl]
>>> Y: X*<Ygeo>
>>> With: steps
>>> ...
>>> Regards
>>> Vasilis
>>> B H U S H A N wrote:
>>>> Dear fluka users and experts,
>>>> USRBDX 101.0 7.0 -31.0 20.0 3.0
>>>> 10.5768PhFluenUD
>>>> USRBDX 0.006 0.00 100.0 &
>>>> This is my usrbdx card for the detection of photon fluence. I want to
calculate fluence of photon passing through the region 20 to 3. how
can
>>>> i
>>>> get it?
>>>> I was doing fluence calulation by this way
>>>> IIIrd column in *_tab.lis (i.e photon/cm^2/GeV/primary
particle)
>>>> multiply this by the differece of Ist and IInd column of the same
row.
>>>> Is this right? if not please reply me the correct way. does
>>>> the
>>>> same way is useful for neutron fluence calculation?
>>>> Also I have installed flair, in flair for plotting USR-1D and in
>>>> detector information of this there are 6 different option for Y axix
selection. like
>>>> i) Y
>>>> ii) Y x <X>
>>>> iii) Y x <Xgeo>
>>>> iv) Y x Xl
>>>> v) Y x Xh
>>>> vi) Y x DX
>>>> I want to know that which one is correct for fluence calculation(i
was
>>>> using last one) and where other option will be useful. what are its
meaning.
>>>> THank you very much.
>>>> Bhushan
>>>> =09*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!
=09!* =09Mr. PATIL BHUSHANKUMAR JAGNNATH. =09=09!*
=09*!=09 Junior Reasearch Fellow *!
>>>> !*=09=09DEPARTMENT OF PHYSICS,=09 !*
>>>> =09!*=09=09UNIVERSITY OF PUNE,=09=09=09=09!*
>>>> =09*!=09=09GANESHKHIND, PUNE=09=09=09=09*!
>>>> =09!*=09=09MAHARASHTRA=09=09=09=09=09!*
>>>> =09*!=09=09INDIA.=09=09=09=09=09=09*!
>>>> =09!*=09=09411007=09=09=09=09=09=09!*
>>>> =09*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!
>>>> Alternate mail-id amolbhushan2003@yahoo.com
>>>> Web-page http://physics.unipune.ernet.in/~bjp/
>>> --
>>> This message has been scanned for viruses and
>>> dangerous content by MailScanner, and is
>>> believed to be clean.
>> Mr. Bhushankumar Jagnnath Patil.
>> Microtron Accelerator Lab.
>> Department of Physics,
>> University of Pune,
>> Ganeshkhind, Pune 411007
>> Ph. No. Off. 020-25692678 Ext. 421
>> Mob. 9823968377
>> Webpage http://physics.unipune.ernet.in/~bjp
>> Alternate e-mail amolbhushan@gmail.com
>
> --=20
> Alberto Fass=F2
> SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025
> Phone: (1 650) 926 4762 Fax: (1 650) 926 3569
> fasso@slac.stanford.edu
> --1334196822-2031075482-1200075830=3D:32695--
>
>
> --
> This message has been scanned for viruses and
> dangerous content by MailScanner, and is
> believed to be clean.
>
Mr. Bhushankumar Jagnnath Patil.
Microtron Accelerator Lab.
Department of Physics,
University of Pune,
Ganeshkhind, Pune 411007
Ph. No. Off. 020-25692678 Ext. 421
Mob. 9823968377
Webpage http://physics.unipune.ernet.in/~bjp
Alternate e-mail amolbhushan@gmail.com
-- This message has been scanned for viruses and dangerous content by MailScanner, and is believed to be clean.
This archive was generated by hypermail 2.1.6 : Sat Jan 12 2008 - 23:29:07 CET