54 : fDetector(det),fScintCollID(-1),fPMTCollID(-1),fVerbose(0),
55 fPMTThreshold(1),fForcedrawphotons(
false),fForcenophotons(
false)
105 G4int n_trajectories = 0;
106 if (trajectoryContainer) n_trajectories = trajectoryContainer->
entries();
110 for (
G4int i=0; i<n_trajectories; i++){
137 int n_hit = scintHC->
entries();
142 for(
int i=0;i<n_hit;i++){
143 edep=(*scintHC)[i]->GetEdep();
145 eWeightPos += (*scintHC)[i]->GetPos()*
edep;
164 G4cout <<
"\tEnergy weighted position of hits in LXe : "
169 G4cout <<
"\tTotal energy deposition in scintillator : "
178 for(
G4int i=0;i<pmts;i++){
179 fHitCount += (*pmtHC)[i]->GetPhotonCount();
180 reconPos+=(*pmtHC)[i]->GetPMTPos()*(*pmtHC)[i]->GetPhotonCount();
185 (*pmtHC)[i]->SetDrawit(
false);
195 G4cout <<
"\tReconstructed position of hits in LXe : "
210 G4cout <<
"\tNumber of photons that hit PMTs in this event : "
214 G4cout <<
"\tNumber of photons produced by scintillation in this event : "
216 G4cout <<
"\tNumber of photons produced by cerenkov in this event : "
218 G4cout <<
"\tNumber of photons absorbed (OpAbsorption) in this event : "
220 G4cout <<
"\tNumber of photons absorbed at boundaries (OpBoundary) in "
222 G4cout <<
"Unaccounted for photons in this event : "