Message: Problems with Radioactive Decay: no decay lines Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Problems with Radioactive Decay: no decay lines 

Forum: Hadronic Processes
Date: 05 Jun, 2014
From: Kseniia Rezynkina <Kseniia Rezynkina>

Hello,

I'm trying to simulate the decay of an isomer in 210Ra implanted into a silicon strip detector with use of Geant4.10.00.p01. For this I define the ions through particle gun like this:

    fParticleGun  = new G4ParticleGun(nofParticles);

    G4int Z = 88, A = 210;                  //210Ra
    G4double ionCharge   = 0.*eplus;
    G4double excitEnergy = 2100.*keV;

    G4ParticleDefinition* ion = G4IonTable::GetIonTable()->GetIon(Z,A,excitEnergy);
    ion->SetPDGStable(false);
    fParticleGun->SetParticleDefinition(ion);
    fParticleGun->SetParticleCharge(ionCharge);
 
    fParticleGun->GeneratePrimaryVertex(event);

For the physics list I use, for a start, the one from extended/rdecay02 example.

The code compiles without any problems, but the spectrum in the output is not what I expect.

I would have thought to see the transitions to the ground state and then the alpha-decay from the g.s. => I want to see the alpha-lines in the spectrum. But I don't. Same happens (see attachment) when I set excitEnergy to zero (directly produce the ground state).

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2014/06/05/02.02-8651--06-05_at_11.00.10_AM.png

In the output I get a lot of this warning:

G4WT1 > 
-------- WWWW ------- G4Exception-START -------- WWWW -------
*** G4Exception : HAD_RDM_011
      issued by : G4RadioactiveDecay::GetMeanFreePath()
Ion has negative lifetime -1001 but is not stable.  Setting mean free path to DBL_MAX

*** This is just a warning message. ***

Could you please tell me what may be wrong with my code?

Thanks in advance, Kseniia

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

1 None: Re: Problems with Radioactive Decay: no decay lines   (michel maire - 05 Jun, 2014)
(_ Question: Re: Problems with Radioactive Decay: no decay lines   (Kseniia Rezynkina - 16 Jun, 2014)
(_ None: Re: Problems with Radioactive Decay: no decay lines   (michel maire - 17 Jun, 2014)
 Add Message Add Message
to: "Problems with Radioactive Decay: no decay lines"

 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 ]