33 #include "SteppingAction.hh"
35 #include "HistoManager.hh"
91 if (target) targetName = target->
GetName();
92 nuclearChannel +=
" + " + targetName +
" --> ";
101 analysis->
FillH1(ih,energy);
107 nuclearChannel += partName +
" + ";
112 const std::vector<const G4Track*>* secondary
114 for (
size_t lp=0; lp<(*secondary).size(); lp++) {
115 particle = (*secondary)[lp]->GetDefinition();
118 G4double energy = (*secondary)[lp]->GetKineticEnergy();
127 else if (type ==
"nucleus") ih = 7;
128 else if (type ==
"meson") ih = 8;
129 else if (type ==
"baryon") ih = 9;
130 if (ih > 0) analysis->
FillH1(ih,energy);
132 if (type ==
"nucleus") {
152 analysis->
FillH1(ih,Pbal);
156 const G4String conver[] = {
"0",
"",
"2 ",
"3 ",
"4 ",
"5 ",
"6 ",
"7 ",
"8 ",
"9 ",
157 "10 ",
"11 ",
"12 ",
"13 ",
"14 ",
"15 ",
"16 "};
158 std::map<G4ParticleDefinition*,G4int>::iterator ip;
160 particle = ip->first;
162 G4int nb = ip->second;
163 if (nb > kMax) nb =
kMax;
168 name =
"gamma or e-";
171 nuclearChannel += Nb +
name;