Dear Tabea,
It is indeed true that the particle generation number of a proton remains unchanged after an elastic nuclear interaction, and that all products of inelastic interactions (protons, alphas etc.) have a generation number which is one higher than the generation number of the incoming proton.
Concerning your questions:
1) no, there is a priori no such distinction made (which can become arbitrary in particular at high energies). If you are interested in more detail about the inelastic interaction products (particle type, energy, direction etc.) you can access this information in the USDRAW entry of the MGDRAW routine (to be activated with USERDUMP card). I include a brief code snippet below. I know this was not your primary question, but it might be of help anyway.
2) yes they are explicitly transported if you enable ion transport via the IONTRANS card, e.g. with WHAT(2)=-2.0. However, the transport only takes place if their energy is above the transport thresholds which can be set via the PART-THR card. Note that the transport thresholds of ions heavier than alphas are derived from alpha cuts (see Note 2 under the PART-THR card in the manual). Hence if you set for example the following then the transport cut for ions is 0.25 keV/u (which is the minimum recommended value):
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
PART-THR -1.0E-6 4-HELIUM
Cheers,
Anton
* Code snippet for printing information about inel. collision products
* to a file on unit 79 (note: evaporation products are not included in the print-out
* but can be accesses via the FHEAVY common)
* ICODE for inelastic interactions is 101:
IF (ICODE .EQ. 101) THEN
* print particle ID, energy and generation number of projectile:
WRITE( 79, * ) "ParentID=", JTRACK, ", E(GeV)=",ETRACK,
& ", Generation=", LTRACK
* print mass number/atomic number of residual nuclei:
WRITE( 79, * ) "Res. Nuclei=", ICRES,"/", IBRES
* print particle ID and kinetic energy of secondaries:
WRITE( 79, * ) " Secondaries (",NP," ):"
DO IP = 1, NP
WRITE( 79, * ) KPART(IP),TKI(IP)
END DO
END IF
________________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Tabea Pfuhl [t.pfuhl_at_gsi.de]
Sent: 14 July 2017 10:05
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: How to distinguish between inelastically scattered ejectile protons: primary or target proton
Dear FLUKA experts,
I study the depth dose profile of 200 MeV protons in a water block target.
I would like to distinguish between dose which originates from primary
and secondary protons. Therefore I use USRBIN as a scorer with an
additional AUXSCORE card (with BEAMPART or PROTON). As I learned from
older discussions and the manual, if I score the dose of BEAMPART I only
get the dose of incident protons and the dose of protons which have
undergone elastic scattering processes. As soon as the proton
experiences inelastic nuclear processes it is no longer counted as a
BEAMPART. However, in inelastic nuclear reactions it is also possible
that the primary proton is one of the ejectiles. But in that case the
proton's generation number would increase so that it is not counted as
BEAMPART anymore.
So here are my questions:
1) Is there a possibility to distinguish between protons produced in
inelastic nuclear reactions which originate from the beam and actual
target protons that originate from the target nucleus?
2) Also I wondered if the heavy recoils (e.g. O-15, N-15) which are
produced in these collisions (protons on water) are explicitly transported?
Dear FLUKA experts, I would really appreciate some help concerning my
questions.
Best regards,
Tabea
--
Tabea Pfuhl
Biophysics - BIO
Phone / Telefon: +49 6159 71 1893
t.pfuhl_at_gsi.de
GSI Helmholtzzentrum für Schwerionenforschung GmbH
Planckstraße 1, 64291 Darmstadt, Germany, www.gsi.de
Commercial Register / Handelsregister: Amtsgericht Darmstadt, HRB 1528
Managing Directors / Geschäftsführung:
Professor Dr. Paolo Giubellino, Ursula Weyrich, Jörg Blaurock
Chairman of the Supervisory Board / Vorsitzender des GSI-Aufsichtsrats:
State Secretary / Staatssekretär Dr. Georg Schütte
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id¬c_info
Received on Fri Jul 14 2017 - 15:49:43 CEST