35 #include "DetectorConstruction.hh"
36 #include "EventAction.hh"
47 fEdeposit(0.), fEdeposit2(0.),
48 fTrackLen(0.), fTrackLen2(0.),
49 fProjRange(0.), fProjRange2(0.),
50 fNbOfSteps(0), fNbOfSteps2(0),
51 fStepSize(0.), fStepSize2(0.)
69 std::ios::fmtflags mode =
G4cout.flags();
70 G4cout.setf(std::ios::fixed,std::ios::floatfield);
75 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
79 <<
"\n Total Energy deposited = " <<
G4BestUnit(fEdeposit,
"Energy")
87 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
91 <<
"\n Track length of primary track = " <<
G4BestUnit(fTrackLen,
"Length")
100 <<
" (from full dE/dx)" <<
G4endl;
107 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
110 <<
"\n Projected range = " <<
G4BestUnit(fProjRange,
"Length")
119 rms = fNbSteps2 - fNbSteps*fNbSteps;
120 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
123 G4cout <<
"\n Nb of steps of primary track = " << fNbSteps <<
" +- " << rms;
127 if (rms>0.) rms = std::sqrt(rms);
else rms = 0.;
131 <<
"\t Step size= " <<
G4BestUnit(fStepSize,
"Length")
143 <<
"\n absorbed = " << absorbed <<
" %"
144 <<
" transmit = " << transmit <<
" %"
145 <<
" reflected = " << reflected <<
" %" <<
G4endl;
148 G4cout.setf(mode,std::ios::floatfield);
156 const Run* localRun =
static_cast<const Run*
>(run);
171 for (
G4int i=1; i<=nbOfAbsor; ++i) {
virtual void Merge(const G4Run *)
#define G4BestUnit(a, b)
#define G4_USE_G4BESTUNIT_FOR_VERBOSE 1
G4GLOB_DLL std::ostream G4cout
DetectorConstruction * fDetector
Detector construction class to demonstrate various ways of placement.
G4double fCsdaRange[MaxAbsor]
G4double fXfrontNorm[MaxAbsor]
virtual void Merge(const G4Run *)
Run(DetectorConstruction *)