64 if (theParentName ==
"mu+") {
71 }
else if (theParentName ==
"mu-") {
81 G4cout <<
"G4MuonDecayChannel:: constructor :";
82 G4cout <<
" parent particle is not muon but ";
100 if (
this != &right) {
139 const int N_DAUGHTER=3;
144 for (
G4int index=0; index<N_DAUGHTER; index++){
146 sumofdaughtermass += daughtermass[index];
154 delete parentparticle;
157 G4double daughtermomentum[N_DAUGHTER];
159 G4double xmax = (1.0+daughtermass[0]*daughtermass[0]/parentmass/parentmass);
165 G4double EMax=parentmass/2-daughtermass[0];
167 const size_t MAX_LOOP=1000;
169 for (
size_t loop1=0; loop1 <MAX_LOOP; ++loop1){
171 for (
size_t loop2 =0; loop2<MAX_LOOP; ++loop2){
174 if (gam <= x*(1.-x))
break;
178 if ( Ene >= (1.-Ee))
break;
186 G4double costheta,sintheta,rphi,rtheta,rpsi;
187 costheta= 1.-2./Ee-2./Ene+2./Ene/Ee;
188 sintheta=std::sqrt(1.-costheta*costheta);
196 rot.
set(rphi,rtheta,rpsi);
199 daughtermomentum[0]=std::sqrt(Ee*Ee*EMax*EMax+2.0*Ee*EMax * daughtermass[0]);
210 daughtermomentum[1]=std::sqrt(Ene*Ene*EMax*EMax+2.0*Ene*EMax * daughtermass[1]);
220 daughtermomentum[2]=std::sqrt(Enm*Enm*EMax*EMax +2.0*Enm*EMax*daughtermass[2]);
221 G4ThreeVector direction2(-Ene/Enm*sintheta,0,-Ee/Enm-Ene/Enm*costheta);
226 direction2 * daughtermomentum[2]);
232 G4cout <<
"G4MuonDecayChannel::DecayIt ";
233 G4cout <<
" create decay products in rest frame " <<
G4endl;
G4int PushProducts(G4DynamicParticle *aParticle)
G4String ** daughters_name
HepRotation & set(const Hep3Vector &axis, double delta)
G4double GetPDGMass() const
G4MuonDecayChannel & operator=(const G4MuonDecayChannel &)
void SetBR(G4double value)
virtual ~G4MuonDecayChannel()
void ClearDaughtersName()
G4ParticleDefinition * G4MT_parent
static constexpr double twopi
G4ParticleDefinition ** G4MT_daughters
static constexpr double rad
void SetNumberOfDaughters(G4int value)
void CheckAndFillParent()
void SetParent(const G4ParticleDefinition *particle_type)
G4int GetVerboseLevel() const
G4GLOB_DLL std::ostream G4cout
virtual G4DecayProducts * DecayIt(G4double)
void SetDaughter(G4int anIndex, const G4ParticleDefinition *particle_type)
void CheckAndFillDaughters()