**From:** Alberto Fasso' (*fasso@SLAC.Stanford.EDU*)

**Date:** Fri Jan 11 2008 - 19:23:50 CET

**Previous message:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Maybe in reply to:**B H U S H A N: "usrbdx + flair"**Next in thread:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Reply:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

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

**Previous message:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Maybe in reply to:**B H U S H A N: "usrbdx + flair"**Next in thread:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Reply:**Mr. Bhushankumar Jagnnath Patil.: "Re: usrbdx + flair"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

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