34 #include "SteppingAction.hh"
35 #include "RunAction.hh"
58 if (process == 0)
return;
62 if (processName !=
"GammaToMuPair")
return;
71 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
73 Eplus = (*secondary)[lp]->GetTotalEnergy();
74 Pplus = (*secondary)[lp]->GetMomentum();
76 Eminus = (*secondary)[lp]->GetTotalEnergy();
77 Pminus = (*secondary)[lp]->GetMomentum();
81 G4double xPlus = Eplus/EGamma, xMinus = Eminus/EGamma;
83 G4double GammaPlus=EGamma*xPlus/fMuonMass;
84 G4double GammaMinus=EGamma*xMinus/fMuonMass;
88 analysisManager->FillH1(1,1./(1.+std::pow(thetaPlus*GammaPlus,2)));
89 analysisManager->FillH1(2,std::log10(thetaPlus*GammaPlus));
91 analysisManager->FillH1(3,std::log10(thetaMinus*GammaMinus));
92 analysisManager->FillH1(4,std::log10(std::fabs(thetaPlus *GammaPlus
93 -thetaMinus*GammaMinus)));
95 analysisManager->FillH1(5,xPlus);
96 analysisManager->FillH1(6,xMinus);