34 #include "SteppingAction.hh"
35 #include "PrimaryGeneratorAction.hh"
36 #include "RunAction.hh"
37 #include "HistoManager.hh"
43 :fPrimary(prim),fRunAction(RuAct)
75 analysisManager->FillH1(
id,energy);
80 analysisManager->FillH1(
id,costeta);
86 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
87 G4double charge = (*secondary)[lp]->GetDefinition()->GetPDGCharge();
88 if (charge != 0.) {
id = 3; }
else {
id = 5; }
89 G4double energy = (*secondary)[lp]->GetKineticEnergy();
90 analysisManager->FillH1(
id,energy);
93 G4ThreeVector direction = (*secondary)[lp]->GetMomentumDirection();
95 analysisManager->FillH1(
id,costeta);
98 if (charge != 0.) { fRunAction->
SumeTransf(energy); }