Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
G4NeutronDecay Class Reference

#include <G4NeutronDecay.hh>

Inheritance diagram for G4NeutronDecay:
Collaboration diagram for G4NeutronDecay:

Public Member Functions

 G4NeutronDecay (const G4ParticleDefinition *theParentNucleus, const G4double &theBR, const G4double &Qvalue, const G4double &excitation, const G4Ions::G4FloatLevelBase &flb)
 
virtual ~G4NeutronDecay ()
 
virtual G4DecayProductsDecayIt (G4double)
 
virtual void DumpNuclearInfo ()
 
- Public Member Functions inherited from G4NuclearDecay
 G4NuclearDecay (const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
 
virtual ~G4NuclearDecay ()
 
G4RadioactiveDecayMode GetDecayMode ()
 
G4double GetDaughterExcitation ()
 
G4Ions::G4FloatLevelBase GetFloatingLevel ()
 
G4ParticleDefinitionGetDaughterNucleus ()
 
void SetHLThreshold (G4double HLT)
 
G4double GetHLThreshold ()
 
- Public Member Functions inherited from G4VDecayChannel
 G4VDecayChannel (const G4String &aName, G4int Verbose=1)
 
 G4VDecayChannel (const G4String &aName, const G4String &theParentName, G4double theBR, G4int theNumberOfDaughters, const G4String &theDaughterName1, const G4String &theDaughterName2="", const G4String &theDaughterName3="", const G4String &theDaughterName4="")
 
virtual ~G4VDecayChannel ()
 
G4int operator== (const G4VDecayChannel &right) const
 
G4int operator!= (const G4VDecayChannel &right) const
 
G4int operator< (const G4VDecayChannel &right) const
 
const G4StringGetKinematicsName () const
 
G4double GetBR () const
 
G4int GetNumberOfDaughters () const
 
G4ParticleDefinitionGetParent ()
 
G4ParticleDefinitionGetDaughter (G4int anIndex)
 
G4int GetAngularMomentum ()
 
const G4StringGetParentName () const
 
const G4StringGetDaughterName (G4int anIndex) const
 
G4double GetParentMass () const
 
G4double GetDaughterMass (G4int anIndex) const
 
void SetParent (const G4ParticleDefinition *particle_type)
 
void SetParent (const G4String &particle_name)
 
void SetBR (G4double value)
 
void SetNumberOfDaughters (G4int value)
 
void SetDaughter (G4int anIndex, const G4ParticleDefinition *particle_type)
 
void SetDaughter (G4int anIndex, const G4String &particle_name)
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
void DumpInfo ()
 
G4double GetRangeMass () const
 
void SetRangeMass (G4double val)
 
virtual G4bool IsOKWithParentMass (G4double parentMass)
 
void SetPolarization (const G4ThreeVector &)
 
const G4ThreeVectorGetPolarization () const
 

Additional Inherited Members

- Protected Member Functions inherited from G4VDecayChannel
void ClearDaughtersName ()
 
void CheckAndFillDaughters ()
 
void CheckAndFillParent ()
 
 G4VDecayChannel ()
 
 G4VDecayChannel (const G4VDecayChannel &)
 
G4VDecayChanneloperator= (const G4VDecayChannel &)
 
G4double DynamicalMass (G4double massPDG, G4double width, G4double maxDev=+1.) const
 
- Protected Attributes inherited from G4NuclearDecay
const G4RadioactiveDecayMode theMode
 
- Protected Attributes inherited from G4VDecayChannel
G4String kinematics_name
 
G4double rbranch
 
G4int numberOfDaughters
 
G4Stringparent_name
 
G4String ** daughters_name
 
G4double rangeMass
 
G4ThreeVector parent_polarization
 
G4ParticleTableparticletable
 
G4int verboseLevel
 
G4ParticleDefinitionG4MT_parent
 
G4ParticleDefinition ** G4MT_daughters
 
G4double G4MT_parent_mass
 
G4doubleG4MT_daughters_mass
 
G4doubleG4MT_daughters_width
 
G4Mutex daughtersMutex
 
G4Mutex parentMutex
 
- Static Protected Attributes inherited from G4VDecayChannel
static const G4String noName = " "
 

Detailed Description

Definition at line 44 of file G4NeutronDecay.hh.

Constructor & Destructor Documentation

G4NeutronDecay::G4NeutronDecay ( const G4ParticleDefinition theParentNucleus,
const G4double theBR,
const G4double Qvalue,
const G4double excitation,
const G4Ions::G4FloatLevelBase flb 
)

Definition at line 45 of file G4NeutronDecay.cc.

49  : G4NuclearDecay("neutron decay", Neutron, excitationE, flb), transitionQ(Qvalue)
50 {
51  SetParent(theParentNucleus); // Store name of parent nucleus, delete G4MT_parent
52  SetBR(branch);
53 
55  G4IonTable* theIonTable =
57  G4int daughterZ = theParentNucleus->GetAtomicNumber();
58  G4int daughterA = theParentNucleus->GetAtomicMass() - 1;
59  SetDaughter(0, theIonTable->GetIon(daughterZ, daughterA, excitationE, flb) );
60  SetDaughter(1, "neutron");
61 }
void SetBR(G4double value)
G4ParticleDefinition * GetIon(G4int Z, G4int A, G4int lvl=0)
Definition: G4IonTable.cc:503
int G4int
Definition: G4Types.hh:78
G4int GetAtomicNumber() const
G4IonTable * GetIonTable() const
void SetNumberOfDaughters(G4int value)
G4int GetAtomicMass() const
G4NuclearDecay(const G4String &channelName, const G4RadioactiveDecayMode &mode, const G4double &excitation, const G4Ions::G4FloatLevelBase &floatingLevel)
static G4ParticleTable * GetParticleTable()
void SetParent(const G4ParticleDefinition *particle_type)
void SetDaughter(G4int anIndex, const G4ParticleDefinition *particle_type)

Here is the call graph for this function:

G4NeutronDecay::~G4NeutronDecay ( )
virtual

Definition at line 64 of file G4NeutronDecay.cc.

65 {}

Member Function Documentation

G4DecayProducts * G4NeutronDecay::DecayIt ( G4double  )
virtual

Implements G4VDecayChannel.

Definition at line 68 of file G4NeutronDecay.cc.

69 {
70  // Fill G4MT_parent with theParentNucleus (stored by SetParent in ctor)
72 
73  // Fill G4MT_daughters with neutron and residual nucleus (stored by SetDaughter)
75 
76  G4double neutronMass = G4MT_daughters[1]->GetPDGMass();
77  // Excitation energy included in PDG mass
78  G4double nucleusMass = G4MT_daughters[0]->GetPDGMass();
79 
80  // Q value was calculated from atomic masses.
81  // Use it to get correct neutron energy.
82  G4double cmMomentum = std::sqrt(transitionQ*(transitionQ + 2.*neutronMass)*
83  (transitionQ + 2.*nucleusMass)*
84  (transitionQ + 2.*neutronMass + 2.*nucleusMass) )/
85  (transitionQ + neutronMass + nucleusMass)/2.;
86 
87  // Set up final state
88  // parentParticle is set at rest here because boost with correct momentum
89  // is done later
90  G4DynamicParticle parentParticle(G4MT_parent, G4ThreeVector(0,0,0), 0.0);
91  G4DecayProducts* products = new G4DecayProducts(parentParticle);
92 
93  G4double costheta = 2.*G4UniformRand()-1.0;
94  G4double sintheta = std::sqrt(1.0 - costheta*costheta);
96  G4ThreeVector direction(sintheta*std::cos(phi),sintheta*std::sin(phi),
97  costheta);
98 
99  G4double KE = std::sqrt(cmMomentum*cmMomentum + neutronMass*neutronMass)
100  - neutronMass;
101  G4DynamicParticle* daughterparticle =
102  new G4DynamicParticle(G4MT_daughters[1], direction, KE, neutronMass);
103  products->PushProducts(daughterparticle);
104 
105  KE = std::sqrt(cmMomentum*cmMomentum + nucleusMass*nucleusMass) - nucleusMass;
106  daughterparticle =
107  new G4DynamicParticle(G4MT_daughters[0], -1.0*direction, KE, nucleusMass);
108  products->PushProducts(daughterparticle);
109 
110  // Energy conservation check
111  // For neutron decays, do final energy check against reaction Q value
112  // which is well-measured using atomic mass differences. Nuclear masses
113  // should not be used since they are not usually directly measured and we
114  // always decay atoms and not fully stripped nuclei.
115  /*
116  G4int nProd = products->entries();
117  G4DynamicParticle* temp = 0;
118  G4double Esum = 0.0;
119  for (G4int i = 0; i < nProd; i++) {
120  temp = products->operator[](i);
121  Esum += temp->GetKineticEnergy();
122  }
123  G4double eCons = (transitionQ - Esum)/keV;
124  if (eCons > 1.e-07) G4cout << " Neutron decay check: Ediff (keV) = " << eCons << G4endl;
125  */
126  return products;
127 }
void CheckAndFillDaughters()
CLHEP::Hep3Vector G4ThreeVector
G4int PushProducts(G4DynamicParticle *aParticle)
G4ParticleDefinition * G4MT_parent
G4ParticleDefinition ** G4MT_daughters
static constexpr double rad
Definition: G4SIunits.hh:149
static constexpr double twopi
Definition: G4SIunits.hh:76
#define G4UniformRand()
Definition: Randomize.hh:97
G4double GetPDGMass() const
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

void G4NeutronDecay::DumpNuclearInfo ( )
virtual

Implements G4NuclearDecay.

Definition at line 130 of file G4NeutronDecay.cc.

131 {
132  G4cout << " G4NeutronDecay for parent nucleus " << GetParentName() << G4endl;
133  G4cout << " decays to " << GetDaughterName(0) << " + " << GetDaughterName(1)
134  << " with branching ratio " << GetBR() << "% and Q value "
135  << transitionQ << G4endl;
136 }
G4double GetBR() const
const G4String & GetParentName() const
G4GLOB_DLL std::ostream G4cout
const G4String & GetDaughterName(G4int anIndex) const
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:


The documentation for this class was generated from the following files: