33 #include "SteppingAction.hh"
35 #include "HistoManager.hh"
56 Run* run =
static_cast<Run*
>(
77 analysisManager->
FillH1(
id,energy);
82 analysisManager->
FillH1(
id,costeta);
87 const std::vector<const G4Track*>* secondary
89 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
90 G4double charge = (*secondary)[lp]->GetDefinition()->GetPDGCharge();
91 if (charge != 0.) {
id = 3; }
else {
id = 5; }
92 G4double energy = (*secondary)[lp]->GetKineticEnergy();
93 analysisManager->
FillH1(
id,energy);
96 G4ThreeVector direction = (*secondary)[lp]->GetMomentumDirection();
98 analysisManager->
FillH1(
id,costeta);
101 if (charge != 0.) { run->
SumeTransf(energy); }