Re: [fluka-discuss]: activate radioactive decays

From: Shirley Li <shirleyli.phys_at_gmail.com>
Date: Wed, 22 Jan 2014 11:05:40 -0500

Dear Francesco,

Thank you for your reply! It worked~
Is there a time variable to show when the interaction happens? Since ATRACK
is associated with JTRACK, I assume it doesn't describe the decay time.

Thanks again!


On Mon, Jan 20, 2014 at 10:50 AM, Francesco Cerutti <
Francesco.Cerutti_at_cern.ch> wrote:

>
> Dear Shirley,
>
> yes, RADDECAY asking for semi-analogue decay of nuclear residues, together
> with the implementation of a customized usdraw scoring, look like a viable
> approach for your event-by-event analysis.
>
> If so, what's the ICODE for radioactive decay and does the decaying
>> isotope
>> have a corresponding JTRACK value? If not, what's the right way to do
>> this?
>>
>
> As written in the USDRAW header, the radioactive decay code is 110.
> Despite the fact that you asked for ion transport, the decaying isotope
> will hardly be above the transport threshold (you are taking the default
> value of 25keV/n, which in principle can be lowered by a factor 100 with a
> PART-THR card setting the 4-HELIUM threshold to the minimum value of 1keV)
> and anyway their decay will be performed only once they are at rest (no
> in-flight radioactive decay is yet available).
> In USDRAW, under the proper ICODE condition, you can identify the isotope
> at the origin of the decay chain (Z->ICHTAR,A->IBTAR, after inclusion of
> the (RESNUC) common) and the emitted radiation (looping over the NP
> secondaries stored in (GENSTK)). In case of subsequent decays, the
> INFEXT(IP=1,NP) variable identifies the respective parent nucleus, coded as
> (A + Z * 1000) (+ 1000000 in case of isomers).
>
> Kind regards
>
> Francesco
>
> **************************************************
> Francesco Cerutti
> CERN-EN/STI
> CH-1211 Geneva 23
> Switzerland
> tel. ++41 22 7678962
> fax ++41 22 7668854
>
> On Sun, 19 Jan 2014, Shirley Li wrote:
>
> Dear FLUKA users,
>> I'm simulating cosmic-ray muons interact in detectors, producing
>> backgrounds. I'm not sure what card(s) to use to activate radioactive
>> decays. So far I have RADDECAY card. I want to analyse the decay in an
>> event-by-event basis in mgdraw.f (USDRAW). Is this the right way to do it?
>> If so, what's the ICODE for radioactive decay and does the decaying
>> isotope
>> have a corresponding JTRACK value? If not, what's the right way to do
>>
> this?
>
>>
>> Thank you very much!
>> Shirley
>>
>>
Received on Wed Jan 22 2014 - 18:13:14 CET