35 #include "TrackingAction.hh"
36 #include "PrimaryGeneratorAction.hh"
40 #include "RunAction.hh"
41 #include "HistoManager.hh"
49 :fPrimary(prim),fHistoManager(histo)
68 if(selectedGeneratorAction==2)
75 else if(selectedGeneratorAction==3)
84 if (dOmega > 0.) fHistoManager->
FillHisto(
id, alpha, 1./dOmega);
93 G4double psi = std::atan2(u2*um, u1*um);
94 if (psi < 0.) psi +=
twopi;
97 dOmega = (1. - std::cos(alphaMax))*dpsi;
98 if (dOmega > 0.) fHistoManager->
FillHisto(
id, psi, 1./dOmega);
100 else if(selectedGeneratorAction==4)
104 if (r <= 0.0)
return;
113 if (dv > 0.) fHistoManager->
FillHisto(
id, r, 1./dv);
120 dv = fDeltaR3*
twopi*std::sin(theta)*dteta;
121 if (dv > 0.) fHistoManager->
FillHisto(
id, theta, 1./dv);
127 if (phi < 0.) phi +=
twopi;
129 dv = 2*fDeltaR3*dphi;
130 if (dv > 0.) fHistoManager->
FillHisto(
id, phi, 1./dv);
139 if (dOmega > 0.) fHistoManager->
FillHisto(
id, alpha, 1./dOmega);
148 G4double psi = std::atan2(u2*um, u1*um);
149 if (psi < 0.) psi +=
twopi;
151 dOmega = (fCosAlphaMin - fCosAlphaMax)*dpsi;
152 if (dOmega > 0.) fHistoManager->
FillHisto(
id, psi, 1./dOmega);