Re: Use of Deq

Hi Sunil,

Any entry into USRTRACK bins is multplied with a conversion factor in
units of pSv * cm^2. Thus, fluence (spectral fluence and integrals) are
in units of pSv/primary (=p/cmq/pr * pSv *cmq) and tracklength is in units
of pSv * cm^3/primary (=cm/pr * pSv *cmq). The former is your required
result.

Best regards
Stefan

On Wed, 2 Apr 2008, Sunil C wrote:

> Hi Stefan
>
> On processing the USRTRACK outputs using ustsuw.f I get two results in the
> *_sum.lis.
> one has unit p/cmq/pr and the other cm/pr. While computing dose equivalent
> using the deq99.f the correct one will the integral one as I suppose the deq.f
> will be folding the fluence with the conversion coefficients bin by bin?
> So the cm/pr should be the one giving the pSv/primary?
> Am I right?
> What does the p/cmq/pr give in terms of pSv?
>
> Best Regards
> Sunil
> Stefan Roesler wrote:
> > Hi Sunil,
> >
> > I don't have NCRP 144 and IAEA 188; thus can only comment on your input:
> >
> > - Photonuclear interaction are off by default. Therefore, you'll have to
> > uncomment PHOTONUC. - deq99c.f does not apply to USRYIELD.
> >
> > - Thus, the only card which folds with dose-equivalent factors in your input
> > is USRTRACK. Is this the dose output which you are refering to?
> >
> > - You are correct, results are in pSv/primary.
> >
> > - I'd recommend using EMFCUT (both for production and transport thresholds)
> > so that you are sure of your cutoffs.
> >
> > - I guess that you have verified in the output file that deq99c folding was
> > indeed applied?
> >
> >
> > I take the liberty to cc fluka-discuss as the discussiuon might also be of
> > interest to others (or know a better answer).
> >
> > Regards
> > Stefan
> >
> >
> > On Mon, 31 Mar 2008, Sunil C wrote:
> >
> >
> > > Hi Stefan
> > > I am having a normalization problem in a work of mine.
> > > Attached is a typical input where I use a cylindrical detector after a
> > > certain
> > > shield thickness to score photon dose using deq99.f routine.The primary
> > > beam
> > > is electron incident on 10 radiation lengths thickness of several targets.
> > > I
> > > vary the shield thickness and score the dose. On plotting the dose rate
> > > after
> > > shield against the shield thickness, the well known transmission curve
> > > can be
> > > obtained. I get a neat graph except of course at the higher thickness.
> > > Extrapolating the graph to zero thickness should give the unshielded dose
> > > rate. I am comparing this with the NCRP 144 (figure 3.5 page 45) and
> > > Swanson's
> > > IAEA tecdoc 188. I am getting results 100 times lesser unshielded dose
> > > rate
> > > with FLUKA calculations compared to the Swanson's rule of thumb given in
> > > IAEA
> > > 188 (page 53).
> > >
> > > Am I making a mistake while calculating the dose equivalent?
> > > I understand the units are in pSv/primary.
> > >
> > > Your comments will help me debug the problem.
> > >
> > > Thanks a lot
> > > Best Regards
> > > Sunil
> > >
> > >
> >
> >
>
>

-- 
___________________________________
Stefan Roesler
CERN, SC/RP
CH-1211 Geneva 23
Switzerland
Phone:  +41-22-7679891
Fax:    +41-22-7669639
E-mail: Stefan.Roesler@cern.ch