41 #ifndef G4MuonDecayChannelWithSpin_hh
42 #define G4MuonDecayChannelWithSpin_hh 1
44 #include <CLHEP/Units/PhysicalConstants.h>
92 f_c = (5.+17.*x-34.*x*
x)*(omega+std::log(x))-22.*x+34.*x*x;
93 f_c = (1.-
x)/(3.*x*x)*f_c;
94 f_c = (6.-4.*
x)*
R_c(x,omega)+(6.-6.*
x)*std::log(x) + f_c;
95 f_c = (CLHEP::fine_structure_const/
CLHEP::twopi) * (x*x-x0*x0) * f_c;
104 f_theta = (1.+x+34*x*
x)*(omega+std::log(x))+3.-7.*x-32.*x*x;
105 f_theta = f_theta + ((4.*(1.-
x)*(1.-x))/x)*std::log(1.-x);
106 f_theta = (1.-
x)/(3.*x*x) * f_theta;
107 f_theta = (2.-4.*
x)*
R_c(x,omega)+(2.-6.*
x)*std::log(x)-f_theta;
108 f_theta = (CLHEP::fine_structure_const/
CLHEP::twopi) * (x*x-x0*x0) * f_theta;
virtual ~G4MuonDecayChannelWithSpin()
G4double R_c(G4double x, G4double omega)
G4double F_c(G4double x, G4double x0, G4double omega)
virtual G4DecayProducts * DecayIt(G4double)
G4MuonDecayChannelWithSpin & operator=(const G4MuonDecayChannelWithSpin &)
static const double twopi
G4MuonDecayChannelWithSpin()
G4double F_theta(G4double x, G4double x0, G4double omega)
const G4double x[NPOINTSGL]