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
LET =3D deposited_energy / particle_travel_length / material_density
is not the way FLUKA calculates, if i dont misunderstand, isnt it?
On Mon, 2008-11-03 at 16:37 +0100, Francesco Cerutti wrote:
> Dear Ercan,
> 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 =
> LET spectrum given in counts/[GeV keV/(micrometer g/cm3)] per primary in=
> the first case and in counts/[GeV/n keV/(micrometer g/cm3)] per primary i=
> the second one. Since your kinetic energy intervals are 0.8 GeV and 0.02=
> 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=
> and 0.02 GeV/n, respectively, and give the two (identical) LET spectra in=
> counts/[keV/(micrometer g/cm3)] per primary.
> With respect to your usreou.f implementation, note that LET values are=20
> indeed gotten, for particles crossing a given boundary, as unrestricted=
> stopping power in any material you wish (specified again in WHAT(6) of th=
> USRYIELD continuation card - as you did - provided that it is assigned to=
> some region). In principle this is not the energy lost up to that boundar=
> divided by the travelled path (the stopping power depends on energy).
> Francesco Cerutti
> CH-1211 Geneva 23
> tel. ++41 22 7678962
> fax ++41 22 7668854
()".'."() Ercan Pilicer
( (T) ) Uludag University, Turkey
(=C3=B6) (=C3=B6) High Energy Physics Department
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