33 #include "SteppingAction.hh"
34 #include "RunAction.hh"
57 if (process == 0)
return;
61 if (processName !=
"GammaToMuPair")
return;
70 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
72 Eplus = (*secondary)[lp]->GetTotalEnergy();
73 Pplus = (*secondary)[lp]->GetMomentum();
75 Eminus = (*secondary)[lp]->GetTotalEnergy();
76 Pminus = (*secondary)[lp]->GetMomentum();
80 G4double xPlus = Eplus/EGamma, xMinus = Eminus/EGamma;
87 if(0.0 == thetaPlus || 0.0 == thetaMinus) {
88 G4cout <<
"SteppingAction: "
89 <<
"thetaPlus= " << thetaPlus <<
" thetaMinus= " << thetaMinus
90 <<
" gamPlus= " << GammaPlus <<
" gamMinus= " << GammaMinus
91 <<
" " << thetaPlus *GammaPlus - thetaMinus*GammaMinus <<
G4endl;
94 analysisManager->
FillH1(1,1./(1.+std::pow(thetaPlus*GammaPlus,2)));
95 analysisManager->
FillH1(2,std::log10(thetaPlus*GammaPlus));
97 analysisManager->
FillH1(3,std::log10(thetaMinus*GammaMinus));
98 analysisManager->
FillH1(4,std::log10(std::fabs(thetaPlus *GammaPlus
99 -thetaMinus*GammaMinus)));
101 analysisManager->
FillH1(5,xPlus);
102 analysisManager->
FillH1(6,xMinus);