59 PBound(0), Parameterisation(0), HMaker(0)
69 PBound(0), Parameterisation(0), HMaker(0)
148 CrossShower = DirectionShower.
cross(OrthoShower);
152 G4int CosPhi[4] = {1,0,-1,0};
153 G4int SinPhi[4] = {0,1,0,-1};
156 G4int NlateralInside=0;
159 for(
int i=0; i<4 ;i++)
162 Position = InitialPositionShower +
164 R*CosPhi[i]*OrthoShower +
165 R*SinPhi[i]*CrossShower ;
172 if(NlateralInside==4) filter=
true;
211 CrossShower = DirectionShower.
cross(OrthoShower);
242 G4int NSpotDeposited =0;
252 if(Bound < stepLength)
262 ZEndStep=ZEndStep+
Dz;
267 LastEneIntegral = EneIntegral;
270 (EneIntegral-LastEneIntegral)*Energy);
271 LastNspIntegral = NspIntegral;
273 DNsp =
std::max(1., std::floor( (NspIntegral-LastNspIntegral)
280 DNsp =
std::max(1., std::floor( (1.- NspIntegral)
283 EnergyNow = EnergyNow - DEne;
296 StepLenght = StepLenght + Dz/2.00;
297 NewPositionShower = NewPositionShower +
298 StepLenght*DirectionShower;
299 StepLenght = Dz/2.00;
302 for (
int i = 0; i < DNsp; i++)
317 Dz/DNsp*DirectionShower*(i+1/2.-DNsp/2.) +
318 RSpot*std::cos(PhiSpot)*OrthoShower +
319 RSpot*std::sin(PhiSpot)*CrossShower;
326 while(EnergyNow > 0.0 && Bound> 0.0);