Message: Force a nuclear decay Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Force a nuclear decay 

Forum: Hadronic Processes
Date: 23 Jul, 2012
From: Kareem Kazkaz <Kareem Kazkaz>

Hello everyone.

I am using the Radioactive Decay Manager in a simulation to create the particles I need. I'm manually setting the global time of the event when said event is generated. Unfortunately, that's only when the radioisotope is itself created--it then decays only after a random time has passed (based on that radioisotope's half-life, of course).

I thought I might be able to handle this by subtracting the time of the first step of whatever track has a parent ID of 1 (which would be the actual decay time of the primary radioisotope), but then I run into the fact that the particle not always a radioisotope (e.g., it might be a neutron), and sometimes I might have stable primary particles. In cases where the primary particle isn't a radioisotope, I don't want to subtract the time of the secondary track / step 1.

A second problem I'm running into is the limitations of the double data type. There are only about 15 decimal places of accuracy in the double, which means if the primary radioisotope has a half-life of a year or more, the accuracy is exceeded (because the time is stored internally in units of nanoseconds, and a year is roughly 1e16 nanoseconds).

One solution to my conundrum is to force the primary nucleus to decay instantly. But I can't figure out how to force a nuclear decay. The class G4DynamicParticle has a SetPreAssignedDecayProperTime method, but it doesn't seem to be applicable to radioisotopes.

Does anybody have any suggestions on how to force a nuclear decay at a specific time? Or any other ideas on how to handle the timing of a delayed decay?

Thanks for any and all help.

          Kareem

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Force a nuclear decay   (Susana Cebrian - 28 Aug, 2012)
 Add Message Add Message
to: "Force a nuclear decay"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]