35 #include "DetectorConstruction.hh"
79 std::vector<G4DynamicParticle*>* fvect,
113 epsilon = sqrt(epsilonsq);
116 onecost = (1.-
epsilon)/(epsilon*E0_m);
117 sint2 = onecost*(2.-onecost);
118 greject = 1. - epsilon*sint2/(1.+ epsilonsq);
129 G4double dirx = sinTeta*cos(Phi), diry = sinTeta*sin(Phi), dirz = cosTeta;
137 gamDirection1.
rotateUz(gamDirection0);
138 G4double gamEnergy1 = epsilon*gamEnergy0;
146 G4double eKinEnergy = gamEnergy0 - gamEnergy1;
150 = gamEnergy0*gamDirection0 - gamEnergy1*gamDirection1;
151 eDirection = eDirection.
unit();
156 fvect->push_back(dp);