41 #ifndef G4MuonDecayChannelWithSpin_hh
42 #define G4MuonDecayChannelWithSpin_hh 1
98 parent_polarization = polar;
103 return parent_polarization;
108 G4double omega = std::log(EMMU/EMASS);
112 f_c = (5.+17.*x-34.*x*
x)*(omega+std::log(x))-22.*x+34.*x*x;
113 f_c = (1.-
x)/(3.*x*x)*f_c;
114 f_c = (6.-4.*
x)*R_c(x)+(6.-6.*
x)*std::log(x) + f_c;
115 f_c = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_c;
122 G4double omega = std::log(EMMU/EMASS);
126 f_theta = (1.+x+34*x*
x)*(omega+std::log(x))+3.-7.*x-32.*x*x;
127 f_theta = f_theta + ((4.*(1.-
x)*(1.-x))/x)*std::log(1.-x);
128 f_theta = (1.-
x)/(3.*x*x) * f_theta;
129 f_theta = (2.-4.*
x)*R_c(x)+(2.-6.*
x)*std::log(x)-f_theta;
130 f_theta = (CLHEP::fine_structure_const/CLHEP::twopi) * (x*x-x0*x0) * f_theta;
virtual ~G4MuonDecayChannelWithSpin()
void SetPolarization(G4ThreeVector)
virtual G4DecayProducts * DecayIt(G4double)
G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &)
const G4ThreeVector & GetPolarization() const