33 #include "RunAction.hh"
34 #include "DetectorConstruction.hh"
35 #include "PhysicsList.hh"
37 #include "PrimaryGeneratorAction.hh"
52 fAnalysisManager(0), fDetector(det), fPhysics(phys), fKinematic(kin),
54 fEdeptot(0.), fEniel(0.), fNbPrimarySteps(0), fRange(0)
95 if (nbBins < nbmin) nbBins = nbmin;
105 if (nbofEvents == 0)
return;
115 G4cout <<
"\n The run consists of " << nbofEvents <<
" "<< particle <<
" of "
118 << material->
GetName() <<
" (density: "
128 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
136 G4cout <<
" Mean number of primary steps = "<< nstep <<
G4endl;
150 if (tallyNumber > 0) {
152 G4cout <<
"\n---------------------------------------------------------\n";
153 G4cout <<
" Cumulated Doses : \tEdep \tEdep/Ebeam \tDose" <<
G4endl;
154 for (
G4int j=0; j < tallyNumber; ++j) {
158 G4cout <<
" tally " << j <<
": \t \t"
163 G4cout <<
"\n---------------------------------------------------------\n";
170 for (
G4int j=1; j<3; ++j) {
203 const G4String id[] = {
"h0",
"h1",
"h2",
"h3" };
206 "Edep (MeV/mm) along absorber ",
207 "Edep (MeV/mm) along absorber zoomed",
221 if (k == 1) activ =
true;