33 #include "MyMollerBhabhaModel.hh"
68 if (kineticEnergy < th) {
88 dedx = log(2.0*(tau + 2.0)/eexc2) - 1.0 - beta2
89 + log((tau-d)*d) + tau/(tau-
d)
90 + (0.5*d*d + (2.0*tau + 1.)*log(1. - d/tau))/gamma2;
99 dedx = log(2.0*(tau + 2.0)/eexc2) + log(tau*d)
100 - beta2*(tau + 2.0*d - y*(3.0*d2
101 + y*(d - d3 + y*(d2 - tau*d3 + d4))))/tau;
119 if (dedx < 0.0) dedx = 0.0;
125 if (kineticEnergy >=
lowLimit) dedx *= sqrt(tkin/kineticEnergy);
126 else dedx *= sqrt(tkin*kineticEnergy)/
lowLimit;