65 if (theParentName ==
"neutron") {
72 }
else if (theParentName ==
"anti_neutron") {
82 G4cout <<
"G4NeutronBetaDecayChannel:: constructor :";
83 G4cout <<
" parent particle is not neutron but ";
103 if (
this != &right) {
147 for (
G4int index=0; index<3; index++){
149 sumofdaughtermass += daughtermass[index];
151 G4double xmax = parentmass-sumofdaughtermass;
159 delete parentparticle;
171 const size_t MAX_LOOP=10000;
172 for (
size_t loop_counter=0; loop_counter <MAX_LOOP; ++loop_counter){
174 p = std::sqrt(x*(x+2.0*dm));
176 r = p*(x+dm)*(xmax-x)*(xmax-
x)*(1.0+
aENuCorr*p/(x+dm)*
w);
191 daughtermomentum[0] = p;
193 direction0 = rm * direction0;
200 eNu = (parentmass-daughtermass[2])*(parentmass+daughtermass[2])+(dm*dm)-2.*parentmass*(x+dm);
201 eNu /= 2.*(parentmass+p*w-(x+dm));
204 G4double sinn = std::sqrt((1.0-cosn)*(1.0+cosn));
206 G4ThreeVector direction1(sinn*std::cos(phin), sinn*std::sin(phin), cosn);
207 direction1 = rm * direction1;
214 eP = parentmass-eNu-(x+dm)-daughtermass[2];
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);
227 G4cout <<
"G4NeutronBetaDecayChannel::DecayIt ";
228 G4cout <<
" create decay products in rest frame " <<
G4endl;
void ClearDaughtersName()
void SetBR(G4double value)
G4NeutronBetaDecayChannel & operator=(const G4NeutronBetaDecayChannel &)
CLHEP::Hep3Vector G4ThreeVector
CLHEP::HepRotation G4RotationMatrix
G4int PushProducts(G4DynamicParticle *aParticle)
G4ParticleDefinition * G4MT_parent
G4ParticleDefinition ** G4MT_daughters
const G4double w[NPOINTSGL]
G4GLOB_DLL std::ostream G4cout
void SetNumberOfDaughters(G4int value)
static const double twopi
G4NeutronBetaDecayChannel()
G4int GetVerboseLevel() const
G4double GetPDGMass() const
const G4double x[NPOINTSGL]
void SetParent(const G4ParticleDefinition *particle_type)
void SetDaughter(G4int anIndex, const G4ParticleDefinition *particle_type)
virtual ~G4NeutronBetaDecayChannel()
G4String ** daughters_name
virtual G4DecayProducts * DecayIt(G4double)