51 const auto localRun =
static_cast<const OpNoviceRun*
>(run);
66 fMie2 += localRun->fMie2;
76 if(numberOfEvent == 0)
78 auto TotNbofEvents = G4double(numberOfEvent);
84 rmsCerenkov = std::sqrt(rmsCerenkov);
93 rmsScint = std::sqrt(rmsScint);
101 rmsRayleigh = std::sqrt(rmsRayleigh);
107 G4double rmsAbsorption =
109 if(rmsAbsorption > 0.)
110 rmsAbsorption = std::sqrt(rmsAbsorption);
115 fMie2 /= TotNbofEvents;
118 rmsMie = std::sqrt(rmsMie);
126 rmsBoundary = std::sqrt(rmsBoundary);
130 G4int prec = G4cout.precision(3);
131 G4cout <<
"\n ======================== run summary ======================\n";
133 G4cout <<
"Primary particle was: " <<
fParticle->GetParticleName()
135 G4cout <<
"Number of events: " << numberOfEvent << G4endl;
137 G4cout <<
"Average number of Cerenkov photons created per event: "
139 G4cout <<
"Average number of scintillation photons created per event: "
142 G4cout <<
"Average number of optical Rayleigh interactions per event: "
144 G4cout <<
"Average number of optical absorption interactions per event: "
146 G4cout <<
"Average number of optical Mie interactions per event: "
148 G4cout <<
"Average number of optical boundary interactions per event: "
152 G4cout.precision(prec);