60 const LXeRun* localRun =
static_cast<const LXeRun*
>(run);
84 G4cout <<
"\n ======================== run summary ======================\n";
89 G4cout <<
"The run was " << n_evt <<
" events." <<
G4endl;
95 if (rms_hits > 0.) rms_hits = std::sqrt(rms_hits/n_evt);
97 G4cout <<
"Number of hits per event:\t " << hits <<
" +- " << rms_hits
102 G4double rms_hitsAbove = hitsAbove2 - hitsAbove*hitsAbove;
103 if (rms_hitsAbove > 0.) rms_hitsAbove = std::sqrt(rms_hitsAbove/n_evt);
104 else rms_hitsAbove = 0.;
106 G4cout <<
"Number of hits per event above threshold:\t " << hitsAbove
107 <<
" +- " << rms_hitsAbove <<
G4endl;
111 G4double rms_scint = scint2 - scint*scint;
112 if (rms_scint > 0.) rms_scint = std::sqrt(rms_scint/n_evt);
115 G4cout <<
"Number of scintillation photons per event :\t " << scint <<
" +- "
120 G4double rms_ceren = ceren2 - ceren*ceren;
121 if (rms_ceren > 0.) rms_ceren = std::sqrt(rms_ceren/n_evt);
124 G4cout <<
"Number of Cerenkov photons per event:\t " << ceren <<
" +- "
129 G4double rms_absorb = absorb2 - absorb*absorb;
130 if (rms_absorb > 0.) rms_absorb = std::sqrt(rms_absorb/n_evt);
131 else rms_absorb = 0.;
133 G4cout <<
"Number of absorbed photons per event :\t " << absorb <<
" +- "
138 G4double rms_bdry = bdry2 - bdry*bdry;
139 if (rms_bdry > 0.) rms_bdry = std::sqrt(rms_bdry/n_evt);
142 G4cout <<
"Number of photons absorbed at boundary per event:\t " << bdry
143 <<
" +- " << rms_bdry <<
G4endl;
149 if (rms_en > 0.) rms_en = std::sqrt(rms_en/n_evt);
152 G4cout <<
"Total energy deposition in scintillator per event:\t " << en/
keV
153 <<
" +- " << rms_en/
keV <<
" keV." <<
G4endl;
G4int fPhotonCount_Scint2
G4int fBoundaryAbsorptionCount2
virtual void Merge(const G4Run *run)
static constexpr double keV
G4int fPMTsAboveThreshold2
G4int fBoundaryAbsorptionCount
G4int fPhotonCount_Ceren2
G4GLOB_DLL std::ostream G4cout
virtual void Merge(const G4Run *)
G4int fPMTsAboveThreshold