34 #include "TrackingAction.hh"
35 #include "PrimaryGeneratorAction.hh"
39 #include "RunAction.hh"
40 #include "HistoManager.hh"
62 if(selectedGeneratorAction==2)
70 else if(selectedGeneratorAction==0)
77 analysis->
FillH1(
id, cosalpha);
83 if (psi < 0.) psi +=
twopi;
86 else if(selectedGeneratorAction==3)
93 analysis->
FillH1(
id, cosalpha);
102 G4double psi = std::atan2(u2*um, u1*um);
103 if (psi < 0.) psi +=
twopi;
104 analysis->
FillH1(
id, psi);
107 else if(selectedGeneratorAction==4)
111 if (r <= 0.0)
return;
120 if (dv > 0.) analysis->
FillH1(
id, r, 1./dv);
126 analysis->
FillH1(
id, costheta);
132 if (phi < 0.) phi +=
twopi;
133 analysis->
FillH1(
id, phi);
140 analysis->
FillH1(
id, cosalpha);
149 G4double psi = std::atan2(u2*um, u1*um);
150 if (psi < 0.) psi +=
twopi;
151 analysis->
FillH1(
id, psi);