Re: [fluka-discuss]: Not all isotopes decays

From: Paola Sala <paola.sala_at_mi.infn.it>
Date: Sun, 23 Feb 2020 18:34:47 +0100

Dear Alice
sorry for the late answer.

When you record products from radioactive decays in mgdraw (usdraw with
icode=110), you will get in a single call all the products from the full
decay chain, not from a single isotope. The decay products carry a flag
that identifies the specific isotope: infext(kp). Did you check with this
one?
Be careful also with isomeric states (flagged in infext(kp) )

 However * some of the decays do not emit particles, even though the decay
chain is foolwed **. this happens when the information on the particle
energy is not known in the decay database. Thus, for some of them, you'll
get a residual without having any product. One example is the isomeric
state of 16N.

If you wish to intercept rellay each residual nucleus, you can use the
usrrnc.f routine, that is called at each resnucle filling (as the comscw.f
for energy deposition) provided USERWEIG has been requested with WHAT(5)
> 0.

As for thresholds: decay products are generated without checking on
production/transport thresholds. Tranpost thresholds are applied as oon as
their transport starts. So yes, you'll have electrons generated below
threshold

Regards
Paola


> Dear Paola,
>
> thanks for your help. I am using semi-analogue REDDECAY (what(1)=2) and
> no DCYSCORE, in my mgdraw.f (USRDRAW part) I record all the particles
> doing an interaction with iCode=110, if I remember correctly this is the
> code for isotope decays. I also have some EMF cuts (2 Mev) but in the
> spectra of the isotope decay products it doesn't seem to impact, meaning
> that I have beta decays with products with energy lower than the
> threshold I put in EMFCUT. Is this expected?
>
> I attach the input file.
>
> Thnaks again.
>
> Best,
>
> Alice
>
> Il 14-Feb-20 10:16 AM, Paola Sala ha scritto:
>> Dear Alice,
>> it would be useful if you could send me your input file.
>> In any case, may I ask you which type of RADDECAY do you set?
>> (semi-analogue or activation study, value of what(1)).
>> And: do you associate a DCYSCORE to your scorings?
>> Regards
>> Paola
>>> Dear all,
>>>
>>> I am running a simulation of muons doing spallation in water, I am
>>> using
>>> RESNUCLEI card to check isotopes produced and REDDECAY to check their
>>> decay products. I found that the isotope yield is pretty reasonable and
>>> the energy spectra of the decay products too, however, not all my
>>> isotopes are decaying, meaning that the number of isotopes is always
>>> greater than the number of decays. I notice also that this is not
>>> uniform for different type of isotopes (for example only 22% of 16N do
>>> decay,  44% for 18N, 75 %17N, 81% for 15C, while for other like 8He or
>>> 8Li, 12B, 13O, 12C ... is more than 95%). Moreover i checked that is
>>> not
>>> related to the geometry of my simulation nor the thresholds I put. Does
>>> anyone know if I am missing something in the physics or in the
>>> simulation? Any comment would be appreciated.
>>>
>>> Thanks a lot.
>>>
>>> Best,
>>>
>>> Alice
>>>
>>> __________________________________________________________________________
>>> You can manage unsubscription from this mailing list at
>>> https://www.fluka.org/fluka.php?id=acc_info
>>>
>>
>> Paola Sala
>> INFN Milano
>> tel. Milano +39-0250317374
>> tel. CERN +41-227679148
>>
>


Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148

__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Sun Feb 23 2020 - 20:18:30 CET

This archive was generated by hypermail 2.3.0 : Sun Feb 23 2020 - 20:18:38 CET