41 #ifndef G4MuonDecayChannelWithSpin_hh
42 #define G4MuonDecayChannelWithSpin_hh 1
44 #include <CLHEP/Units/PhysicalConstants.h>
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;
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()
CLHEP::Hep3Vector G4ThreeVector
void SetPolarization(G4ThreeVector)
G4double F_theta(G4double x, G4double x0)
virtual G4DecayProducts * DecayIt(G4double)
G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &)
G4ThreeVector parent_polarization
G4MuonDecayChannelWithSpin()
G4double F_c(G4double x, G4double x0)
const G4ThreeVector & GetPolarization() const