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


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

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