Re: [fluka-discuss]: How to distinguish between inelastically scattered ejectile protons: primary or target proton

From: George Kharashvili <georgek_at_jlab.org>
Date: Thu, 3 Aug 2017 14:10:29 -0400 (EDT)

Dear Tabea,

The AUXSCORE cards have A=0 (except of the one applied to oxygen).
Also, COALESCEnce is not turned on - WHAT(1) in the PHYSICS card is left blank.

Regards,
George



----- Original Message -----
From: "Tabea Pfuhl" <t.pfuhl_at_gsi.de>
To: "FLUKA Discussion List" <fluka-discuss_at_fluka.org>
Sent: Thursday, August 3, 2017 5:32:39 AM
Subject: Re: [fluka-discuss]: How to distinguish between inelastically scattered ejectile protons: primary or target proton

Dear Anton and other FLUKA experts,


thank you for the answer to my questions.


Concerning the second question I added your advised options to my input
file. In my experiment I have a beam of 200 MeV protons passing through
a water target. My goal is to score energy spectra of heavy recoils but
until now I have still experienced some difficulties while trying to do so.

The input file (which I attached to this mail) includes the following cards:

-Physics: Evaporat; New Evap with heavy fragments

-Physics: Coalesce

-Transport: Iontrans; Heavyion (whenever I set Full Transport in flair,
save the file and reopen it with flair the setting mysteriously goes
back to Heavyion??)

-Part-Thr: Energy; 4-He; 1e-5

If I use USRTRACK to score neutrons, it works, so there must be
reactions that create heavy recoils, e.g. O-15. Also if I score
RESNUCLEi I can see that heavy recoils are created with all kinds of
combinations of Z and A.

However, if I try to score the energy spectra of these heavy fragments
(USRTRACK card with Heavyion combined with AUXSCORE card with definition
of Z and A) there is nothing scored. Only for Z=8 and A=16 I can see a
spectrum which I think originates from elastic scattering processes with
the target nuclei. I also included the executable flukadpm3, but this
did not change any of the results.

So my overall question is: What is the key to score energy spectra of
heavy recoils?

There was a similar question asked a while ago but, there has not been a
final answer to that
(https://urldefense.proofpoint.com/v2/url?u=http-3A__www.fluka.org_FLUKA_web-5Farchive_earchive_new-2Dfluka-2Ddiscuss_10544.html&d=DwID-g&c=lz9TcOasaINaaC3U7FbMev2lsutwpI4--09aP8Lu18s&r=VEsjWcjHR_RwJLerQJyK4g&m=disWsN_PhhujzVHTP4a9IAP7qC85bQNmoHXf1w2TQhU&s=-eapZie-0ogBWIIhoG6I5E3lhRn8Vw8f6s68Z43_pms&e= )
so I decided to ask you for help concerning this question.


I would be very thankful for any kind of hint or comment.

Best regards,

Tabea


> On 07/14/2017 02:46 PM, Anton Lechner wrote:
>> 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 athttps://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 Thu Aug 03 2017 - 21:29:24 CEST

This archive was generated by hypermail 2.3.0 : Thu Aug 03 2017 - 21:29:28 CEST