66 {fEnergyDeposit +=
edep; fEnergyDeposit2 += edep*
edep;};
69 {fTrakLenCharged += length; fTrakLenCharged2 += length*length;};
72 {fTrakLenNeutral += length; fTrakLenNeutral2 += length*length;};
75 {
if (std::abs(theta) <= fMscThetaCentral) { fMscEntryCentral++;
76 fMscProjecTheta += theta; fMscProjecTheta2 += theta*theta;}
80 {fNbStepsCharged += nSteps; fNbStepsCharged2 += nSteps*nSteps;};
83 {fNbStepsNeutral += nSteps; fNbStepsNeutral2 += nSteps*nSteps;};
91 {
if (flag == 1) fTransmit[0]++;
92 else if (flag == 2) {fTransmit[0]++; fTransmit[1]++; }};
95 {
if (flag == 1) fReflect[0]++;
96 else if (flag == 2) {fReflect[0]++; fReflect[1]++; }};
101 { fEnergyLeak[
index] += eleak; fEnergyLeak2[
index] += eleak*eleak;};
104 G4double fEnergyDeposit, fEnergyDeposit2;
105 G4double fTrakLenCharged, fTrakLenCharged2;
106 G4double fTrakLenNeutral, fTrakLenNeutral2;
107 G4double fNbStepsCharged, fNbStepsCharged2;
108 G4double fNbStepsNeutral, fNbStepsNeutral2;
109 G4double fMscProjecTheta, fMscProjecTheta2;
112 G4int fNbGamma, fNbElect, fNbPosit;
113 G4int fTransmit[2], fReflect[2];
114 G4int fMscEntryCentral;
116 G4double fEnergyLeak[2], fEnergyLeak2[2];