Re: [fluka-discuss]: activate radioactive decays

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Tue, 28 Jan 2014 15:06:41 +0100

---------- Forwarded message ----------
Date: Sun, 26 Jan 2014 19:55:33 +0100 (CET)
From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
To: Shirley Li <shirleyli.phys_at_gmail.com>
Cc: "fluka-discuss_at_fluka.org" <fluka-discuss_at_fluka.org>
Subject: Re: [fluka-discuss]: activate radioactive decays


Dear Shirley,

when an unstable nucleus is created by a reaction in the prompt phase
(ICODE=101 in USDRAW), ATRACK - age of the projectile from the history
beginning - gives you the interaction time. Then, when its decay takes place
(ICODE=110) in semi-analogue mode, in (GENSTK) you have also AGESEC(IP=1,NP)
giving the (cumulative) decay time associated with the decay radiation. Anyway,
when the latter (electrons, positrons, photons, and - upon request - neutrinos)
is transported, ATRACK still carries its absolute time.
ATRACK and AGESEC are in seconds.

Cheers

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854

On Wed, 22 Jan 2014, Shirley Li wrote:

> 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 Tue Jan 28 2014 - 17:20:26 CET

This archive was generated by hypermail 2.3.0 : Tue Jan 28 2014 - 17:20:29 CET