# Re: USRYIELD card normalization

From: ercan pilicer <ercanpilicer_at_gmail.com>
Date: Wed, 05 Nov 2008 09:12:48 +0200

hi all,

as Francesco pointed that the multiplication factor overlapped two plots
very well. thank you.

as a second point about LET calculation with usreou.f routine by using

[MeV/cm/(g/cm3)]
LET =3D deposited_energy / particle_travel_length / material_density

is not the way FLUKA calculates, if i dont misunderstand, isnt it?

ciao
e.

On Mon, 2008-11-03 at 16:37 +0100, Francesco Cerutti wrote:
> Dear Ercan,
>=20
> according to your (meaningful) choice for WHAT(6) of the continuation=20
> card, USRYIELD is giving d2N/dx1dx2, as yourself wrote in the card=20
> comment. x1 is LET, expressed in keV/(micrometer g/cm3), as in your=20
> plot label. x2 is ion energy (GeV in one case and GeV/n in the other).
> So your y label is wrong: we are not dealing with fluence, but just with =
a=20
> LET spectrum given in counts/[GeV keV/(micrometer g/cm3)] per primary in=
=20
> the first case and in counts/[GeV/n keV/(micrometer g/cm3)] per primary i=
n=20
> the second one. Since your kinetic energy intervals are 0.8 GeV and 0.02=
=20
> GeV/n wide, respectively (both choices are correct), your two curves=20
> are supposed to differ by a factor of 40. Let us multiply them by 0.8 GeV=
=20
> and 0.02 GeV/n, respectively, and give the two (identical) LET spectra in=
=20
> counts/[keV/(micrometer g/cm3)] per primary.
>=20
> With respect to your usreou.f implementation, note that LET values are=20
> indeed gotten, for particles crossing a given boundary, as unrestricted=
=20
> stopping power in any material you wish (specified again in WHAT(6) of th=
e=20
> USRYIELD continuation card - as you did - provided that it is assigned to=
=20
> some region). In principle this is not the energy lost up to that boundar=
y=20
> divided by the travelled path (the stopping power depends on energy).
>=20
> Cheers
>=20
> Francesco
>=20
> **************************************************
> Francesco Cerutti
> CERN-AB
> CH-1211 Geneva 23
> Switzerland
> tel. ++41 22 7678962
> fax ++41 22 7668854
--=20
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
()".'."() Ercan Pilicer
( (T) ) Uludag University, Turkey
(=C3=B6) (=C3=B6) High Energy Physics Department
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Received on Wed Nov 05 2008 - 08:54:39 CET

This archive was generated by hypermail 2.2.0 : Wed Nov 05 2008 - 08:54:41 CET