**From:** Giuseppe Battistoni (*Giuseppe.Battistoni@mi.infn.it*)

**Date:** Tue Sep 19 2006 - 21:21:37 CEST

**Previous message:**Alberto Fasso': "Re: Plotting USRBDX with one angular bin"**In reply to:**Alberto Fasso': "Re: Plotting USRBDX with one angular bin"**Next in thread:**Nicole Patricia Lee Pratt-Boyden: "Re: RE: Plotting USRBDX with one angular bin"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

I can answer:

the *_tab file was created in order to make life easier to users.

It replicates the same numbers available in

the *_sum.lis file, but written in a way that is very very easy to be

readout from any analysis/graphic program.

For instance, some of us developed a few root-based tools to readout the

_tab files in order to produce immediately graphics.

Giuseppe Battistoni

On Tue, 19 Sep 2006, Alberto Fasso' wrote:

*> Date: Tue, 19 Sep 2006 09:21:58 -0700 (PDT)
*

*> From: Alberto Fasso' <fasso@SLAC.Stanford.EDU>
*

*> To: Sebastien WURTH <wurth@ipno.in2p3.fr>
*

*> Cc: fluka-discuss@fluka.org, Paola R. Sala <Paola.Sala@cern.ch>
*

*> Subject: Re: Plotting USRBDX with one angular bin
*

*>
*

*> I know for sure for what concerns the *_sum.lis file, but I don't know about
*

*> the *_tab.lis file, which is a recent addition (by Paola, I think).
*

*> Paola, can you answer this question?
*

*>
*

*> Alberto
*

*>
*

*> On Tue, 19 Sep 2006, Sebastien WURTH wrote:
*

*>
*

*> > All right, apparently, it was not very clear for me, although I saw it
*

*> > (the fact that you must not multiply again by 2 pi or 4 pi) by
*

*> > experience on comparing some USRTRACK an "similar" USRBDX results...
*

*> > But then with the same example (one angular bin), in the *_tab.lis file
*

*> > when it says "integrated over solid angle" at the first lign, it refers
*

*> > to the differential Flux as a function of energy in (Part/GeV/cmq/pr) in
*

*> > *_sum.lis file, I do have to multiply by each energy bin value to obtain
*

*> > integrated flux in part/cmq/pr. Am I right this time ?
*

*> >
*

*> > Thank you.
*

*> > Regards.
*

*> > Sebastien.
*

*> >
*

*> > Alberto Fasso' a écrit :
*

*> >
*

*> > >Be careful: the manual says that you must multiply by 2 pi or 4 pi,
*

*> > >but assuming that you do the integration yourself. The usxsuw program
*

*> > >takes already care of that, so you must not do it again
*

*> > >
*

*> > >Alberto
*

*> > >
*

*> > >On Tue, 19 Sep 2006, Sebastien WURTH wrote:
*

*> > >
*

*> > >
*

*> > >
*

*> > >>Hello,
*

*> > >>
*

*> > >>In your example, you cannot have a neutron fluence, you put WHAT(2) =
*

*> > >>7.0 in your first USRBDX card, this would give you photon fluence (maybe
*

*> > >>you choose only the wrong example to illustrate).
*

*> > >>To answer your question (I was asking me the same one once), I quote the
*

*> > >>manual from USRBDX part :
*

*> > >>*********
*

*> > >>Notes
*

*> > >>1. IMPORTANT! 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.
*

*> > >>Thus, for example, when requesting a fluence or current energy spectrum,
*

*> > >>with no angular distribution, to
*

*> > >>obtain integral binned results (fluence or current in cm−2 per energy
*

*> > >>bin per primary) one must multiply the
*

*> > >>value of each energy bin by the width of the bin (even for logarithmic
*

*> > >>binning), and by 2 pi or 4 pi (depending
*

*> > >>on whether one-way or two-way scoring has been requested).[...]
*

*> > >>*********
*

*> > >>
*

*> > >>Best Regards.
*

*> > >>Sebastien.
*

*> > >>
*

*> > >>
*

*> > >>
*

*> > >>
*

*> > >>Nicole Patricia Lee Pratt-Boyden a écrit :
*

*> > >>
*

*> > >>
*

*> > >>
*

*> > >>>Hello,
*

*> > >>>
*

*> > >>>I'm using a USRBDX card to score neutron fluence (linear angular and
*

*> > >>>energy bins). The programme seems to be running fine, and I'm compiling
*

*> > >>>the data with the usxsuw.f programme.
*

*> > >>>
*

*> > >>>I asked for eight energy bins, and one angular bin.
*

*> > >>>
*

*> > >>>USRBDX 101.0 7.0 -48.0 3.0 4.0 400.0
*

*> > >>>PH-199m
*

*> > >>>USRBDX 6.5E-4 2.5E-4 8.0 0.0 1. &
*

*> > >>>
*

*> > >>>The *_tab file only returns the differential fluence integrated over
*

*> > >>>solid angle; there are no double differential values returned at all,
*

*> > >>>no mention of solid angle boundaries.
*

*> > >>>
*

*> > >>>Assuming this is fine, do I still need to multiply the differential
*

*> > >>>values by 2PI (one-way scoring) in order to get diff. fluence in
*

*> > >>>particles/GeV/cm^2/primary, or do I only do this if I want fluence
*

*> > >>>(/cm^2/pr)?
*

*> > >>>
*

*> > >>>Apologies if this is extremely simple, I was just thrown by the lack of
*

*> > >>>any information for the double differential! Best regards,
*

*> > >>>
*

*> > >>>Nicole Pratt-Boyden.
*

*> > >>>
*

*> > >>>
*

*> > >>>
*

*> > >>>
*

*> > >>
*

*> > >>
*

*> > >>
*

*> > >
*

*> > >
*

*> > >
*

*> >
*

*> >
*

*>
*

*>
*

**Previous message:**Alberto Fasso': "Re: Plotting USRBDX with one angular bin"**In reply to:**Alberto Fasso': "Re: Plotting USRBDX with one angular bin"**Next in thread:**Nicole Patricia Lee Pratt-Boyden: "Re: RE: Plotting USRBDX with one angular bin"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ] [ attachment ]

*
This archive was generated by hypermail 2.1.6
: Tue Sep 19 2006 - 23:14:14 CEST
*