34 #include "SteppingAction.hh"
36 #include "HistoManager.hh"
90 if (target) targetName = target->
GetName();
91 nuclearChannel +=
" + " + targetName +
" --> ";
100 analysis->FillH1(ih,energy);
104 Pbalance += momentum;
106 nuclearChannel += partName +
" + ";
111 const std::vector<const G4Track*>* secondary
113 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
114 particle = (*secondary)[lp]->GetDefinition();
117 G4double energy = (*secondary)[lp]->GetKineticEnergy();
126 else if (type ==
"nucleus") ih = 7;
127 else if (type ==
"meson") ih = 8;
128 else if (type ==
"baryon") ih = 9;
129 if (ih > 0) analysis->FillH1(ih,energy);
131 if (type ==
"nucleus") {
133 analysis->FillH1(12, A);
138 Pbalance += momentum;
147 analysis->FillH1(ih,Q);
149 analysis->FillH1(ih,Pbal);
152 const G4int kMax = 16;
153 const G4String conver[] = {
"0",
"",
"2 ",
"3 ",
"4 ",
"5 ",
"6 ",
"7 ",
"8 ",
"9 ",
154 "10 ",
"11 ",
"12 ",
"13 ",
"14 ",
"15 ",
"16 "};
155 std::map<G4ParticleDefinition*,G4int>::iterator ip;
157 particle = ip->first;
159 G4int nb = ip->second;
160 if (nb > kMax) nb = kMax;
167 nuclearChannel += Nb +
name;
G4ParticleDefinition * GetDefinition() const
void CountProcesses(G4String procName)
const G4String & GetName() const
void SetTargetXXX(G4bool)
CLHEP::Hep3Vector G4ThreeVector
G4double GetStepLength() const
std::map< G4ParticleDefinition *, G4int > fParticleFlag
G4StepStatus GetStepStatus() const
AnalysisManager * analysis
G4TrackStatus GetTrackStatus() const
G4ThreeVector GetMomentum() const
void UserSteppingAction(const G4Step *)
const G4Isotope * GetTargetIsotope()
const G4String & GetParticleName() const
G4StepPoint * GetPreStepPoint() const
void CountNuclearChannel(G4String, G4double)
void SumTrack(G4double track)
static G4Proton * Proton()
const G4String & GetParticleType() const
ExG4HbookAnalysisManager G4AnalysisManager
static G4Neutron * Neutron()
const std::vector< const G4Track * > * GetSecondaryInCurrentStep() const
static const G4double A[nN]
static G4Deuteron * Deuteron()
G4int GetAtomicMass() const
static G4RunManager * GetRunManager()
const G4VProcess * GetProcessDefinedStep() const
G4double energy(const ThreeVector &p, const G4double m)
G4StepPoint * GetPostStepPoint() const
virtual void AbortEvent()
G4double GetKineticEnergy() const
G4Track * GetTrack() const
G4Run * GetNonConstCurrentRun() const
void ParticleCount(G4String, G4double)