41 #ifndef G4MuonDecayChannelWithSpin_hh
42 #define G4MuonDecayChannelWithSpin_hh 1
44 #include <CLHEP/Units/PhysicalConstants.h>
96 f_c = (5.+17.*x-34.*x*
x)*(omega+std::log(x))-22.*x+34.*x*x;
97 f_c = (1.-
x)/(3.*x*x)*f_c;
98 f_c = (6.-4.*
x)*
R_c(x)+(6.-6.*
x)*std::log(x) + f_c;
99 f_c = (CLHEP::fine_structure_const/
CLHEP::twopi) * (x*x-x0*x0) * f_c;
110 f_theta = (1.+x+34*x*
x)*(omega+std::log(x))+3.-7.*x-32.*x*x;
111 f_theta = f_theta + ((4.*(1.-
x)*(1.-x))/x)*std::log(1.-x);
112 f_theta = (1.-
x)/(3.*x*x) * f_theta;
113 f_theta = (2.-4.*
x)*
R_c(x)+(2.-6.*
x)*std::log(x)-f_theta;
114 f_theta = (CLHEP::fine_structure_const/
CLHEP::twopi) * (x*x-x0*x0) * f_theta;
virtual ~G4MuonDecayChannelWithSpin()
G4double F_theta(G4double x, G4double x0)
virtual G4DecayProducts * DecayIt(G4double)
G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &)
static const double twopi
G4MuonDecayChannelWithSpin()
G4double F_c(G4double x, G4double x0)
const G4double x[NPOINTSGL]