34 #include "SteppingAction.hh"
35 #include "RunAction.hh"
36 #include "HistoManager.hh"
74 analysisManager->FillH1(
id,energy);
79 analysisManager->FillH1(
id,costeta);
84 const std::vector<const G4Track*>* secondary
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); }
G4double GetStepLength() const
void SumeTransf(G4double energy)
G4StepStatus GetStepStatus() const
G4TrackStatus GetTrackStatus() const
void CountProcesses(G4String procName)
void UserSteppingAction(const G4Step *)
const G4ThreeVector & GetMomentumDirection() const
double precision function energy(A, Z)
ExG4HbookAnalysisManager G4AnalysisManager
const std::vector< const G4Track * > * GetSecondaryInCurrentStep() const
const G4String & GetProcessName() const
void SumTrack(G4double track)
static G4RunManager * GetRunManager()
const G4VProcess * GetProcessDefinedStep() const
G4StepPoint * GetPostStepPoint() const
virtual void AbortEvent()
G4double GetKineticEnergy() const
G4Track * GetTrack() const