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

**Date:** Sat Jan 12 2008 - 18:51:32 CET

**Previous message:**Alberto Fasso': "Re: usrbdx + flair"**Maybe in reply to:**B H U S H A N: "usrbdx + flair"**Next in thread:**Alberto Fasso': "Re: usrbdx + flair"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

Dear Bhushan,

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

Not completely. The X 2pi part is already done by the usxsuw program

(if you look at the output you will see a line such as:

# Detector n: 1 piFluenUD (integrated over solid angle)

-----------------------------

# N. of energy intervals 50

=2E..........................................................

so, you need only to do

IIIrd column X (IInd column - Ist column)

The lines above are reproduced form the Fluka manual, Chap. 2

(A FLUKA beginner's guide).

I suggest that you study carefully that Chapter: there are clear

answers to most of your questions.

Alberto

*>
*

*> 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=3D
*

*>> 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=3D
*

*>> hose
*

*>> two energies. If the two energies are the the limits of a single
*

*> bin,
*

*>> you get the fluence in that bin.=3D20
*

*>> 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=3D
*

*>> ,
*

*>> 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 =3D3D Ist column, "High" xh =3D3D IInd column, Y =
*

=3D3D II=3D

*>> 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 ma=
*

ke

*>>>> 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 d=
*

at=3D

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

*>>>>> =3D09*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*=
*

!*!

*> =3D09!* =3D09Mr. PATIL BHUSHANKUMAR JAGNNATH. =3D09=3D09!*
*

*> =3D09*!=3D09 Junior Reasearch Fellow *!
*

*>>>>> !*=3D09=3D09DEPARTMENT OF PHYSICS,=3D09 =
*

!*

*>>>>> =3D09!*=3D09=3D09UNIVERSITY OF PUNE,=3D09=3D09=3D09=3D09!*
*

*>>>>> =3D09*!=3D09=3D09GANESHKHIND, PUNE=3D09=3D09=3D09=3D09*!
*

*>>>>> =3D09!*=3D09=3D09MAHARASHTRA=3D09=3D09=3D09=3D09=3D09!*
*

*>>>>> =3D09*!=3D09=3D09INDIA.=3D09=3D09=3D09=3D09=3D09=3D09*!
*

*>>>>> =3D09!*=3D09=3D09411007=3D09=3D09=3D09=3D09=3D09=3D09!*
*

*>>>>> =3D09*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*!*=
*

!*!

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

*>>
*

*>> --=3D20
*

*>> Alberto Fass=3DF2
*

*>> 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=3D3D: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
*

*>
*

*>
*

*>
*

*>
*

*>
*

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

--1334196310-1161952472-1200160292=3D:30799--

**Previous message:**Alberto Fasso': "Re: usrbdx + flair"**Maybe in reply to:**B H U S H A N: "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:04 CET
*