Message: Re: Alpha energies in Radioactive Decay are wrong! Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

None Re: Alpha energies in Radioactive Decay are wrong! 

Forum: Hadronic Processes
Re: Warning Alpha energies in Radioactive Decay are wrong! (Kseniia Rezynkina)
Date: 18 Jul, 2014
From: Laurent Desorgher <Laurent Desorgher>

Dear Ksenilia,

The following code to add directly after lines 248 of G4NuclearDecayChannel::DecayIt(G4double)

 else if (decayMode == Alpha){
    //Compute Q from nuclear mass, correct for daughterExcitation as it is not
    // in daughter mass at this level and as Qtransition does contain it
    //daughterExcitation is corrected again afterwards when calling  SetParentMass()
    G4double Q_nuc = G4MT_parent->GetPDGMass()-daughterNucleus->GetPDGMass()
                                -G4Alpha::Alpha()->GetPDGMass()-daughterExcitation;
    deltaM=Q_nuc-Qtransition;
  }

will produce the correct alpha energy.

As such it does use the Q value given in the radioactive decay file. So you can also modify this value if needed in you own database. I first try to work with a solution where Q is directly computed from atomic mass but it is at the moment, as far as I know, impossible to get this atomic mass as public method in Geant4. That this why I rely on the Q of the database for the mass correction deltaM. Note that the Q values in the database are taken from ENSDF or from Qcalc program and these values should be correctly computed from atomic masses.

Best regards

Laurent

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

1 None: Re: Alpha energies in Radioactive Decay are wrong!   (Laurent Desorgher - 18 Jul, 2014)
2 None: Re: Alpha energies in Radioactive Decay are wrong!   (Karl Hauschild - 09 Oct, 2014)
 Add Message Add Message
to: "Re: Alpha energies in Radioactive Decay are wrong!"

 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 ]