Re: problem with USRYIELD

From: Sunil C <csunil11_at_gmail.com>
Date: Tue, 11 Jan 2011 14:56:21 +0530

Hi
It is a bit sketchy how you got the surface area from sr-1.

The units of usryield is given as GeV-1sr-1primary-1. (Don't confuse this
with the boundary crossing USRBDX estimator, which has a cm-2 also in its
unit.) Thus to get a cm-2 into your unit, you have to use the relation
mentioned before. Once you do this, you have area in your units and no solid
angle (cm-2 GeV-1 primary-1) but now at a distance (which is 'r'). Multiply
the results by the bin width (GeV) and you have cm-2 primary-1. The USRYIELD
results are now ready to be multiplied with the conversion coefficients
(pSvcm2).

The final result will be in pSv primary-1 at distance 'r'.

Hope it helps.

Sunil

On Thu, Jan 6, 2011 at 9:52 PM, Sunil C <csunil11_at_gmail.com> wrote:

> Hi
> Perhaps you should try the LAM-BIAS card also along with photonuc card.
> There is an example in the manual.
>
> The USRYIELD card scores the yield in unit solid angle(sr-1) at the given
> direction. The half solid angle is defined by the upper and lower angle that
> are input as what(4) and what(5) in the second card. See the answer
> by Giuseppe.
> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/0439.html
>
> The steradian can be converted to unit area at a distance (solid angle =
> dA*cos(theta)/r2), theta=0 if normal to surface. ONce you convert the sr to
> cm-2, you can use the ICRP or Maurizio Pelliccioni's fluence to dose
> conversion coefficients given in units of pSvcm2 and obtain dose equvalent.
>
> Good luck!
>
> Sunil
>
>
> On Mon, Jan 3, 2011 at 1:16 PM, Hamideh Jalali <jalali.hb_at_gmail.com>wrote:
>
>> Dear FLUKA users
>>
>> In my problem a 3GeV electron beam strike a target of iron and I'm
>> going to score gamma and neutron dose equivalent in angles 0 and 90
>> degree with respect to beam direction around the wall of concrete. So
>> I used USRYIELD for this aim but in 0 degree, the results are one or
>> two order less than what I expected and in 90 degree, the results are
>> zero. Please guide me where I am wrong?
>> The other question, as I've got I can set angles in USRYIELD in radian
>> or degree but the results are in sr and if I normalize results to
>> surface area I will obtain dose equivalent in pSv/GeV.sr/primary. Now
>> I want to find dose equivalent in pSv/primary, is it correct if I
>> multiply cumulative spectrum by angular interval width in radian?
>> the input file is attached.
>>
>> Best Regards
>> Hamideh
>>
>
>

--0016e65b52e471eeb804998ead6e
Content-Type: text/html; charset=ISO-8859-1
Content-Transfer-Encoding: quoted-printable

Hi<br>It is a bit sketchy how you got the surface area from sr-1. <br><br>T=
he units of usryield is given as GeV-1sr-1primary-1. (Don&#39;t confuse thi=
s with the boundary crossing USRBDX estimator, which has a cm-2 also in its=
  unit.) Thus to get a cm-2 into your unit, you have to use the relation
men=
tioned before. Once you do this, you have area in your units and no solid a=
ngle (cm-2 GeV-1 primary-1) but now at a distance (which is &#39;r&#39;). M=
ultiply the results by the bin width (GeV) and you have cm-2 primary-1. The=
  USRYIELD results are now ready to be multiplied with the conversion
coeffi=
cients (pSvcm2).<br>
<br>The final result will be in pSv primary-1 at distance &#39;r&#39;.<br><=
br>Hope it helps.<br><br>Sunil<br><br><br><div class=3D"gmail_quote">On Thu=
, Jan 6, 2011 at 9:52 PM, Sunil C <span dir=3D"ltr">&lt;<a href=3D"mailto:c=
sunil11_at_gmail.com">csunil11_at_gmail.com</a>&gt;</span> wrote:<br>
<blockquote class=3D"gmail_quote" style=3D"margin: 0pt 0pt 0pt 0.8ex; borde=
r-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;"><div>Hi</div>Perh=
aps you should try the LAM-BIAS card also along with photonuc card. There i=
s an example in the manual.<div>
<br>The USRYIELD card scores the yield in unit solid angle(sr-1) at the giv=
en direction. The half solid angle is defined by the upper and lower angle =
that are input as what(4) and what(5) in the second card. =A0See the answer=
  by=A0Giuseppe.=A0</div>

<div><a href=3D"http://www.fluka.org/web_archive/earchive/new-fluka-discuss=
/0439.html" target=3D"_blank">http://www.fluka.org/web_archive/earchive/new=
-fluka-discuss/0439.html</a></div><div><br></div><div>The steradian can be =
converted to unit area at a distance (solid angle =3D dA*cos(theta)/r2), th=
eta=3D0 if normal to surface. ONce you convert the sr to cm-2, you can use =
the ICRP or Maurizio Pelliccioni&#39;s fluence to dose conversion=A0coeffic=
ients=A0given in units of pSvcm2 and obtain dose equvalent.</div>

<div><br></div><div>Good luck!</div><div><br></div><font color=3D"#888888">=
<div>Sunil</div></font><div><div></div><div class=3D"h5"><div><br></div><di=
v><br></div><div><div class=3D"gmail_quote">On Mon, Jan 3, 2011 at 1:16 PM,=
  Hamideh Jalali <span dir=3D"ltr">&lt;<a
href=3D"mailto:jalali.hb_at_gmail.com=
" target=3D"_blank">jalali.hb_at_gmail.com</a>&gt;</span> wrote:<br>

<blockquote class=3D"gmail_quote" style=3D"margin: 0pt 0pt 0pt 0.8ex; borde=
r-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">Dear FLUKA users<=
br>
<br>
In my problem a 3GeV electron beam strike a target of iron and I&#39;m<br>
going to score gamma and neutron dose equivalent in angles 0 and 90<br>
degree with respect to beam direction around the wall of concrete. So<br>
I used USRYIELD for this aim but in 0 degree, the results are one or<br>
two order less than what I expected and =A0in 90 degree, the results are<br=
>
zero. Please guide me where I am wrong?<br>
The other question, as I&#39;ve got I can set angles in USRYIELD in radian<=
br>
or degree but the results are in sr and if I normalize results to<br>
surface area I will obtain dose equivalent in pSv/GeV.sr/primary. Now<br>
I want to find dose equivalent in pSv/primary, is it correct if I<br>
multiply cumulative spectrum by angular interval width in radian?<br>
the input file is attached.<br>
<br>
Best Regards<br>
<font color=3D"#888888">Hamideh<br>
</font></blockquote></div><br></div>
</div></div></blockquote></div><br>

--0016e65b52e471eeb804998ead6e--
Received on Wed Jan 12 2011 - 00:38:17 CET

This archive was generated by hypermail 2.2.0 : Wed Jan 12 2011 - 00:38:18 CET