120 for (
G4int index=0; index<3; index++){
122 sumofdaughtermass += daughtermass[index];
132 delete parentparticle;
146 G4double W_mue = (EMMU*EMMU+EMASS*EMASS)/(2.*EMMU);
159 const size_t MAX_LOOP=10000;
160 for (
size_t loop_count =0; loop_count <MAX_LOOP; ++loop_count){
166 x = x0 + rndm*(1.-x0);
172 F_IS = 1./6.*(-2.*x_squared+3.*x-x0_squared);
173 F_AS = 1./6.*std::sqrt(x_squared-x0_squared)*(2.*x-2.+std::sqrt(1.-x0_squared));
175 G_IS = 2./9.*(michel_rho-0.75)*(4.*x_squared-3.*x-x0_squared);
176 G_IS = G_IS + michel_eta*(1.-
x)*x0;
178 G_AS = 3.*(michel_xsi-1.)*(1.-
x);
179 G_AS = G_AS+2.*(michel_xsi*michel_delta-0.75)*(4.*x-4.+std::sqrt(1.-x0_squared));
180 G_AS = 1./9.*std::sqrt(x_squared-x0_squared)*G_AS;
186 const G4double omega = std::log(EMMU/EMASS);
189 G4double F = 6.*F_IS + R_IS/std::sqrt(x_squared-x0_squared);
199 G4double G = 6.*F_AS - R_AS/std::sqrt(x_squared-x0_squared);
201 FG = std::sqrt(x_squared-x0_squared)*F*(1.+(G/F)*ctheta);
204 G4cout<<
"***Problem in Muon Decay *** : FG > FG_max"<<
G4endl;
210 if (FG >= rndm*FG_max)
break;
219 if(energy < EMASS) energy = EMASS;
224 daughtermomentum[0] = std::sqrt(energy*energy - EMASS*EMASS);
226 G4double stheta = std::sqrt(1.-ctheta*ctheta);
248 G4double energy2 = parentmass*(1.0 - x/2.0);
249 G4double vmass = std::sqrt((energy2-daughtermomentum[0])*(energy2+daughtermomentum[0]));
252 G4double sinthetan = std::sqrt((1.0-costhetan)*(1.0+costhetan));
257 G4ThreeVector direction1(sinthetan*cosphin,sinthetan*sinphin,costhetan);
279 G4cout <<
"G4MuonDecayChannelWithSpin::DecayIt ";
280 G4cout <<
" create decay products in rest frame " <<
G4endl;
291 if(n_max<10)n_max=10;
296 L2 += std::pow(x,
n)/(
n*
n);
301 r_c = 2.*L2-(
pi*
pi/3.)-2.;
302 r_c = r_c + omega * (1.5+2.*std::log((1.-x)/x));
303 r_c = r_c - std::log(x)*(2.*std::log(x)-1.);
304 r_c = r_c + (3.*std::log(x)-1.-1./
x)*std::log(1.-x);
G4int PushProducts(G4DynamicParticle *aParticle)
G4String ** daughters_name
G4double F_theta(G4double x, G4double x0, G4double omega)
G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &)
virtual G4DecayProducts * DecayIt(G4double)
Hep3Vector & rotateUz(const Hep3Vector &)
G4MuonDecayChannelWithSpin()
G4double GetPDGMass() const
void ClearDaughtersName()
G4ParticleDefinition * G4MT_parent
static constexpr double twopi
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
void Set4Momentum(const G4LorentzVector &momentum)
G4LorentzVector Get4Momentum() const
G4ParticleDefinition ** G4MT_daughters
static constexpr double rad
G4ThreeVector parent_polarization
void CheckAndFillParent()
G4int GetVerboseLevel() const
G4double R_c(G4double x, G4double omega)
G4GLOB_DLL std::ostream G4cout
virtual ~G4MuonDecayChannelWithSpin()
G4double GetTotalMomentum() const
static constexpr double pi
G4double F_c(G4double x, G4double x0, G4double omega)
HepLorentzVector & boost(double, double, double)
void CheckAndFillDaughters()