# Re: usrbdx + flair

From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Fri Jan 11 2008 - 19:23:50 CET

• Next message: Alberto Fasso': "Re: usrbdx + flair"

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

• Next message: Alberto Fasso': "Re: usrbdx + flair"

This archive was generated by hypermail 2.1.6 : Fri Jan 11 2008 - 21:18:57 CET