From: Alberto Fasso' (fasso@SLAC.Stanford.EDU)
Date: Fri Nov 16 2007 - 18:29:01 CET
V?
Message-ID: <Pine.LNX.4.64.0711160928100.580@noric08.slac.stanford.edu>
MIME-Version: 1.0
Content-Type: MULTIPART/Mixed; boundary=3D"1334195286-823803711-1195080046=
=3D:23609"
Content-ID: <Pine.LNX.4.64.0711160928101.580@noric08.slac.stanford.edu>
Sender: owner-fluka-discuss@smtp2.mi.infn.it
This message is in MIME format. The first part should be readable text,
while the remaining parts are likely unreadable without MIME-aware tools.
--1334195286-823803711-1195080046=3D:23609
Content-Type: TEXT/PLAIN; format=3Dflowed; charset=3DISO-8859-1
Content-Transfer-Encoding: 8BIT
Content-ID: <Pine.LNX.4.64.0711141442341.23609@noric02.slac.stanford.edu>
Dear Niels,
indeed you can hack fluscw.f, but you must use it in a non-standard way.
That means that you don't assign any value to FLUSCW (leave it at the
default value =3D ONEONE), but exploit the fact that fluscw.f is called at
tracklength scoring to manipulate the energy of the ion.
Put the following lines in fluscw.f:
=2E..................................
INCLUDE '(FHEAVY)'
INCLUDE '(PAPROP)'
=2E..................................
FLUSCW =3D ONEONE
IF (-6 .LE. IJ ) THEN
IA =3D IBARCH(IJ)
ELSE IF (IJ .LT. -6) THEN
IA =3D IBHEAV(-IJ)
END IF
PLA =3D PLA/DBLE(IA)
Of course, in your USRTRACK commands must set the maximum and minimum
energy of the spectrum in a way consistent with the fact that it will
be a spectrum of E/n and not of E.
Alberto
On Tue, 13 Nov 2007, Niels Bassler wrote:
> Hi Alberto,
>=20
> thanks for your reply!
> I don't mind hacking fluscw.f if it is possible. But even that I dont
> know how to do:
> Sure I can modify what is being scored (by changing FLUSCW =3D
> WHATEVER_I_WANT ), but I cannot see any possibility to change _where_
> (meaning: in which energy bin) it will be scored.
>=20
> Example:
>=20
> I want to use the USRTRACK & fluscw.f results with my own stoppingpower
> calculations.
>=20
> If I run fluka for some hypothetical situation and get
>=20
> 1 x 100 MeV proton
> 1 x 200 MeV deuteron
> 1 x 300 MeV triton
>=20
> I know how to get three output files with the reults.
> And of course I could bin all in one single output file, but then the
> important information is lost of what particle had one energy:
>=20
> Energy bin:=09=09Fluence:
> 0-100=09=09=091
> 101-200=09=09=091
> 201-300=09=09=091
> 301-400=09=09=090
>=20
> this I cannot use for stopping power calculations, because I have lost
> the information of the isotope composition for Z =3D 1.
>=20
> But if I somehow could bin it in energy/nucleon instead It would be
> cool:
>=20
> Energy/amu bin:=09=09Fluence:
> 0-100=09=09=093
> 101-200=09=09=090
> 201-300=09=09=090
> 301-400=09=09=090
>=20
> because, the stopping power per nucleon for equal Z particles is the
> same.
>=20
> It would then be sufficent for me, to just generate one output table for
> every ion-charge Z.
>=20
> I think it is possible to do it with HBOOK, but there must be an easier
> way??
>=20
>=20
> I hope what I am writing can be understood ... thanks a lot!
>=20
> Cheers
>=20
> Niels
>=20
>=20
>=20
> On Mon, 2007-11-12 at 08:43 -0800, Alberto Fasso' wrote:
>> Hi Niels,
>>=20
>> as far as I know, the only way to change the USRTRACK energy
>> spectrum scoring to GeV/nucleon is to write a fluscw.f user routine.
>> But that is what you are trying to avoid, if I have well understood
>> what you said.
>>=20
>> Regards,
>>=20
>> Alberto
>>=20
>> On Mon, 12 Nov 2007, Niels Bassler wrote:
>>=20
>>> Hi,
>>>=20
>>> I want to score the fluence particle-energy spectrum of heavy ions
>>> (Z>2), and currently I use the USRTRACK card (scoring particle -2)
>>> together with a customized fluscw.f routine (which identifies the ions
>>> via USRDCI). This works fine enough.
>>> However, it is a bit unelegant to have one USRTRACK card for each
>>> possible isotope of each ion species, as both the .f and .inp files get
>>> rather bloated.
>>> Instead, I could directly score all isotopes with one Z into one logica=
l
>>> unit (defined by a single USRTRACK card), but this makes only sense for
>>> my work, if I can score it in terms of GeV/nucleon, and not GeV.
>>>=20
>>> So the question is: is it somehow possible, to change the USRTRACK
>>> energy spectrum scoring to GeV/nucleon instead of GeV? (and if yes, how
>>> do I do it?)
>>>=20
>>> Cheers
>>>=20
>>> Niels
>>>=20
>>>=20
>>>=20
>>>=20
>>=20
>=20
>
--=20
Alberto Fass=F2
SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025
Phone: (1 650) 926 4762 Fax: (1 650) 926 3569
fasso@slac.stanford.edu
--1334195286-823803711-1195080046=3D:23609--
This archive was generated by hypermail 2.1.6 : Sun Nov 18 2007 - 22:47:31 CET