Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
G4NeutronBetaDecayChannel Class Reference

#include <G4NeutronBetaDecayChannel.hh>

Inheritance diagram for G4NeutronBetaDecayChannel:
Inheritance graph
[legend]
Collaboration diagram for G4NeutronBetaDecayChannel:
Collaboration graph
[legend]

Public Member Functions

 G4NeutronBetaDecayChannel (const G4String &theParentName, G4double theBR)
 
virtual ~G4NeutronBetaDecayChannel ()
 
virtual G4DecayProducts * DecayIt (G4double)
 
- 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 G4String & GetKinematicsName () const
 
G4double GetBR () const
 
G4int GetNumberOfDaughters () const
 
G4ParticleDefinition * GetParent ()
 
G4ParticleDefinition * GetDaughter (G4int anIndex)
 
G4int GetAngularMomentum ()
 
const G4String & GetParentName () const
 
const G4String & GetDaughterName (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 G4ThreeVector & GetPolarization () const
 

Protected Member Functions

 G4NeutronBetaDecayChannel (const G4NeutronBetaDecayChannel &)
 
G4NeutronBetaDecayChannel & operator= (const G4NeutronBetaDecayChannel &)
 
 G4NeutronBetaDecayChannel ()
 
- Protected Member Functions inherited from G4VDecayChannel
void ClearDaughtersName ()
 
void CheckAndFillDaughters ()
 
void CheckAndFillParent ()
 
 G4VDecayChannel ()
 
 G4VDecayChannel (const G4VDecayChannel &)
 
G4VDecayChannel & operator= (const G4VDecayChannel &)
 
G4double DynamicalMass (G4double massPDG, G4double width, G4double maxDev=+1.) const
 

Protected Attributes

const G4double aENuCorr
 
- Protected Attributes inherited from G4VDecayChannel
G4String kinematics_name
 
G4double rbranch
 
G4int numberOfDaughters
 
G4String * parent_name
 
G4String ** daughters_name
 
G4double rangeMass
 
G4ThreeVector parent_polarization
 
G4ParticleTable * particletable
 
G4int verboseLevel
 
G4ParticleDefinition * G4MT_parent
 
G4ParticleDefinition ** G4MT_daughters
 
G4double G4MT_parent_mass
 
G4double * G4MT_daughters_mass
 
G4double * G4MT_daughters_width
 
G4Mutex daughtersMutex
 
G4Mutex parentMutex
 

Additional Inherited Members

- Static Protected Attributes inherited from G4VDecayChannel
static const G4String noName = " "
 

Detailed Description

Definition at line 43 of file G4NeutronBetaDecayChannel.hh.

Constructor & Destructor Documentation

G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel ( const G4String &  theParentName,
G4double  theBR 
)

Definition at line 58 of file G4NeutronBetaDecayChannel.cc.

61  :G4VDecayChannel("Neutron Decay"),
62  aENuCorr(-0.102)
63 {
64  // set names for daughter particles
65  if (theParentName == "neutron") {
66  SetBR(theBR);
67  SetParent("neutron");
68  SetNumberOfDaughters(3);
69  SetDaughter(0, "e-");
70  SetDaughter(1, "anti_nu_e");
71  SetDaughter(2, "proton");
72  } else if (theParentName == "anti_neutron") {
73  SetBR(theBR);
74  SetParent("anti_neutron");
75  SetNumberOfDaughters(3);
76  SetDaughter(0, "e+");
77  SetDaughter(1, "nu_e");
78  SetDaughter(2, "anti_proton");
79  } else {
80 #ifdef G4VERBOSE
81  if (GetVerboseLevel()>0) {
82  G4cout << "G4NeutronBetaDecayChannel:: constructor :";
83  G4cout << " parent particle is not neutron but ";
84  G4cout << theParentName << G4endl;
85  }
86 #endif
87  }
88 }
G4VDecayChannel::SetBR
void SetBR(G4double value)
Definition: G4VDecayChannel.hh:294
G4cout
G4GLOB_DLL std::ostream G4cout
G4VDecayChannel::SetNumberOfDaughters
void SetNumberOfDaughters(G4int value)
Definition: G4VDecayChannel.cc:257
G4VDecayChannel::GetVerboseLevel
G4int GetVerboseLevel() const
Definition: G4VDecayChannel.hh:303
G4VDecayChannel::SetParent
void SetParent(const G4ParticleDefinition *particle_type)
Definition: G4VDecayChannel.cc:485
G4VDecayChannel::SetDaughter
void SetDaughter(G4int anIndex, const G4ParticleDefinition *particle_type)
Definition: G4VDecayChannel.cc:325
G4endl
#define G4endl
Definition: G4ios.hh:61

Here is the call graph for this function:

G4NeutronBetaDecayChannel::~G4NeutronBetaDecayChannel ( )
virtual

Definition at line 90 of file G4NeutronBetaDecayChannel.cc.

91 {
92 }
G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel ( const G4NeutronBetaDecayChannel &  right)
protected

Definition at line 94 of file G4NeutronBetaDecayChannel.cc.

95  : G4VDecayChannel(right),
96  aENuCorr(-0.102)
97 {
98 }
G4NeutronBetaDecayChannel::G4NeutronBetaDecayChannel ( )
protected

Definition at line 52 of file G4NeutronBetaDecayChannel.cc.

53  :G4VDecayChannel(),
54  aENuCorr(-0.102)
55 {
56 }

Member Function Documentation

G4DecayProducts * G4NeutronBetaDecayChannel::DecayIt ( G4double  )
virtual

Implements G4VDecayChannel.

Definition at line 128 of file G4NeutronBetaDecayChannel.cc.

129 {
130  // This class describes free neutron beta decay kinemtics.
131  // This version neglects neutron/electron polarization
132  // without Coulomb effect
133 
134 #ifdef G4VERBOSE
135  if (GetVerboseLevel()>1) G4cout << "G4NeutronBetaDecayChannel::DecayIt ";
136 #endif
137 
138  CheckAndFillParent();
139  CheckAndFillDaughters();
140 
141  // parent mass
142  G4double parentmass = G4MT_parent->GetPDGMass();
143 
144  //daughters'mass
145  G4double daughtermass[3];
146  G4double sumofdaughtermass = 0.0;
147  for (G4int index=0; index<3; index++){
148  daughtermass[index] = G4MT_daughters[index]->GetPDGMass();
149  sumofdaughtermass += daughtermass[index];
150  }
151  G4double xmax = parentmass-sumofdaughtermass;
152 
153  //create parent G4DynamicParticle at rest
154  G4ThreeVector dummy;
155  G4DynamicParticle * parentparticle = new G4DynamicParticle( G4MT_parent, dummy, 0.0);
156 
157  //create G4Decayproducts
158  G4DecayProducts *products = new G4DecayProducts(*parentparticle);
159  delete parentparticle;
160 
161  // calculate daughter momentum
162  G4double daughtermomentum[3];
163 
164  // calcurate electron energy
165  G4double x; // Ee
166  G4double p; // Pe
167  G4double dm = daughtermass[0]; //Me
168  G4double w; // cosine of e-nu angle
169  G4double r;
170  G4double r0;
171  const size_t MAX_LOOP=10000;
172  for (size_t loop_counter=0; loop_counter <MAX_LOOP; ++loop_counter){
173  x = xmax*G4UniformRand();
174  p = std::sqrt(x*(x+2.0*dm));
175  w = 1.0-2.0*G4UniformRand();
176  r = p*(x+dm)*(xmax-x)*(xmax-x)*(1.0+aENuCorr*p/(x+dm)*w);
177  r0 = G4UniformRand()*(xmax+dm)*(xmax+dm)*xmax*xmax*(1.0+aENuCorr);
178  if (r > r0) break;
179  }
180 
181  //create daughter G4DynamicParticle
182  // rotation materix to lab frame
183  G4double costheta = 2.*G4UniformRand()-1.0;
184  G4double theta = std::acos(costheta)*rad;
185  G4double phi = twopi*G4UniformRand()*rad;
186  G4RotationMatrix rm;
187  rm.rotateY(theta);
188  rm.rotateZ(phi);
189 
190  // daughter 0 (electron) in Z direction
191  daughtermomentum[0] = p;
192  G4ThreeVector direction0(0.0, 0.0, 1.0);
193  direction0 = rm * direction0;
194  G4DynamicParticle * daughterparticle0
195  = new G4DynamicParticle( G4MT_daughters[0], direction0*daughtermomentum[0]);
196  products->PushProducts(daughterparticle0);
197 
198  // daughter 1 (nutrino) in XZ plane
199  G4double eNu; // Enu
200  eNu = (parentmass-daughtermass[2])*(parentmass+daughtermass[2])+(dm*dm)-2.*parentmass*(x+dm);
201  eNu /= 2.*(parentmass+p*w-(x+dm));
202  G4double cosn = w;
203  G4double phin = twopi*G4UniformRand()*rad;
204  G4double sinn = std::sqrt((1.0-cosn)*(1.0+cosn));
205 
206  G4ThreeVector direction1(sinn*std::cos(phin), sinn*std::sin(phin), cosn);
207  direction1 = rm * direction1;
208  G4DynamicParticle * daughterparticle1
209  = new G4DynamicParticle( G4MT_daughters[1], direction1*eNu);
210  products->PushProducts(daughterparticle1);
211 
212  // daughter 2 (proton) at REST
213  G4double eP; // Eproton
214  eP = parentmass-eNu-(x+dm)-daughtermass[2];
215  G4double pPx = -eNu*sinn;
216  G4double pPz = -p-eNu*cosn;
217  G4double pP = std::sqrt(eP*(eP+2.*daughtermass[2]));
218  G4ThreeVector direction2(pPx/pP*std::cos(phin), pPx/pP*std::sin(phin), pPz/pP);
219  G4DynamicParticle * daughterparticle2
220  = new G4DynamicParticle( G4MT_daughters[2], direction2);
221  products->PushProducts(daughterparticle2);
222 
223 
224  // output message
225 #ifdef G4VERBOSE
226  if (GetVerboseLevel()>1) {
227  G4cout << "G4NeutronBetaDecayChannel::DecayIt ";
228  G4cout << " create decay products in rest frame " <<G4endl;
229  products->DumpInfo();
230  }
231 #endif
232  return products;
233 }
G4VDecayChannel::CheckAndFillDaughters
void CheckAndFillDaughters()
Definition: G4VDecayChannel.hh:324
G4DecayProducts::PushProducts
G4int PushProducts(G4DynamicParticle *aParticle)
Definition: G4DecayProducts.cc:153
G4VDecayChannel::G4MT_parent
G4ParticleDefinition * G4MT_parent
Definition: G4VDecayChannel.hh:192
p
const char * p
Definition: xmltok.h:285
G4VDecayChannel::G4MT_daughters
G4ParticleDefinition ** G4MT_daughters
Definition: G4VDecayChannel.hh:193
test.x
tuple x
Definition: test.py:50
CLHEP::HepRotation::rotateY
HepRotation & rotateY(double delta)
Definition: Rotation.cc:79
rad
static constexpr double rad
Definition: G4SIunits.hh:149
G4int
int G4int
Definition: G4Types.hh:78
twopi
static constexpr double twopi
Definition: G4SIunits.hh:76
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:97
G4cout
G4GLOB_DLL std::ostream G4cout
G4DecayProducts::DumpInfo
void DumpInfo() const
Definition: G4DecayProducts.cc:311
G4VDecayChannel::GetVerboseLevel
G4int GetVerboseLevel() const
Definition: G4VDecayChannel.hh:303
G4ParticleDefinition::GetPDGMass
G4double GetPDGMass() const
Definition: G4ParticleDefinition.hh:122
CLHEP::HepRotation::rotateZ
HepRotation & rotateZ(double delta)
Definition: Rotation.cc:92
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76

Here is the call graph for this function:

G4NeutronBetaDecayChannel & G4NeutronBetaDecayChannel::operator= ( const G4NeutronBetaDecayChannel &  right)
protected

Definition at line 101 of file G4NeutronBetaDecayChannel.cc.

102 {
103  if (this != &right) {
104  kinematics_name = right.kinematics_name;
105  verboseLevel = right.verboseLevel;
106  rbranch = right.rbranch;
107 
108  // copy parent name
109  parent_name = new G4String(*right.parent_name);
110 
111  // clear daughters_name array
112  ClearDaughtersName();
113 
114  // recreate array
115  numberOfDaughters = right.numberOfDaughters;
116  if ( numberOfDaughters >0 ) {
117  if (daughters_name !=0) ClearDaughtersName();
118  daughters_name = new G4String*[numberOfDaughters];
119  //copy daughters name
120  for (G4int index=0; index < numberOfDaughters; index++) {
121  daughters_name[index] = new G4String(*right.daughters_name[index]);
122  }
123  }
124  }
125  return *this;
126 }
G4int
int G4int
Definition: G4Types.hh:78
G4VDecayChannel::kinematics_name
G4String kinematics_name
Definition: G4VDecayChannel.hh:165
G4VDecayChannel::parent_name
G4String * parent_name
Definition: G4VDecayChannel.hh:171
G4VDecayChannel::daughters_name
G4String ** daughters_name
Definition: G4VDecayChannel.hh:173

Here is the call graph for this function:

Member Data Documentation

const G4double G4NeutronBetaDecayChannel::aENuCorr
protected

Definition at line 70 of file G4NeutronBetaDecayChannel.hh.

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