75{
76 if(numberOfEvent == 0)
77 return;
78 auto TotNbofEvents = G4double(numberOfEvent);
79
83 if(rmsCerenkov > 0.)
84 rmsCerenkov = std::sqrt(rmsCerenkov);
85 else
86 rmsCerenkov = 0.;
87
90 G4double rmsScint =
92 if(rmsScint > 0.)
93 rmsScint = std::sqrt(rmsScint);
94 else
95 rmsScint = 0.;
96
100 if(rmsRayleigh > 0.)
101 rmsRayleigh = std::sqrt(rmsRayleigh);
102 else
103 rmsRayleigh = 0.;
104
107 G4double rmsAbsorption =
109 if(rmsAbsorption > 0.)
110 rmsAbsorption = std::sqrt(rmsAbsorption);
111 else
112 rmsAbsorption = 0.;
113
115 fMie2 /= TotNbofEvents;
117 if(rmsMie > 0.)
118 rmsMie = std::sqrt(rmsMie);
119 else
120 rmsMie = 0.;
121
125 if(rmsBoundary > 0.)
126 rmsBoundary = std::sqrt(rmsBoundary);
127 else
128 rmsBoundary = 0.;
129
130 G4int prec = G4cout.precision(3);
131 G4cout << "\n ======================== run summary ======================\n";
132
133 G4cout <<
"Primary particle was: " <<
fParticle->GetParticleName()
135 G4cout << "Number of events: " << numberOfEvent << G4endl;
136
137 G4cout << "Average number of Cerenkov photons created per event: "
139 G4cout << "Average number of scintillation photons created per event: "
141
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: "
150
151 G4cout << G4endl;
152 G4cout.precision(prec);
153}