**From:** Mr. Bhushankumar Jagnnath Patil. (*bjp@physics.unipune.ernet.in*)

**Date:** Sat Jan 12 2008 - 12:57:44 CET

**Previous message:**Alberto Fasso': "Re: usrbdx + flair"**In reply to:**Alberto Fasso': "Re: usrbdx + flair"**Next in thread:**Alberto Fasso': "Re: usrbdx + flair"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

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.

**Previous message:**Alberto Fasso': "Re: usrbdx + flair"**In reply to:**Alberto Fasso': "Re: usrbdx + flair"**Next in thread:**Alberto Fasso': "Re: usrbdx + flair"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Sat Jan 12 2008 - 23:29:07 CET
*