91 anEvent->GetTrajectoryContainer();
93 G4int n_trajectories = 0;
94 if(trajectoryContainer)
95 n_trajectories = trajectoryContainer->entries();
98 if(G4VVisManager::GetConcreteInstance())
100 for(G4int i = 0; i < n_trajectories; ++i)
102 auto trj = (
LXeTrajectory*) ((*(anEvent->GetTrajectoryContainer()))[i]);
103 if(trj->GetParticleName() ==
"opticalphoton")
108 trj->DrawTrajectory();
132 size_t n_hit = scintHC->entries();
133 G4ThreeVector eWeightPos(0.);
135 G4double edepMax = 0;
137 for(
size_t i = 0; i < n_hit; ++i)
139 edep = (*scintHC)[i]->GetEdep();
142 (*scintHC)[i]->GetPos() * edep;
146 G4ThreeVector posMax = (*scintHC)[i]->GetPos();
152 G4AnalysisManager::Instance()->FillH1(7,
fTotE);
157 G4cout <<
"No hits in the scintillator this event." << G4endl;
166 G4cout <<
"\tEnergy weighted position of hits in LXe : "
167 << eWeightPos / mm << G4endl;
172 G4cout <<
"\tTotal energy deposition in scintillator : " <<
fTotE / keV
173 <<
" (keV)" << G4endl;
179 G4ThreeVector reconPos(0., 0., 0.);
180 size_t pmts = pmtHC->entries();
182 for(
size_t i = 0; i < pmts; ++i)
184 fHitCount += (*pmtHC)[i]->GetPhotonCount();
185 reconPos += (*pmtHC)[i]->GetPMTPos() * (*pmtHC)[i]->GetPhotonCount();
192 (*pmtHC)[i]->SetDrawit(
false);
196 G4AnalysisManager::Instance()->FillH1(1,
fHitCount);
204 G4cout <<
"\tReconstructed position of hits in LXe : " << reconPos / mm
209 pmtHC->DrawAllHits();
220 G4cout <<
"\tNumber of photons that hit PMTs in this event : " <<
fHitCount
222 G4cout <<
"\tNumber of PMTs above threshold(" <<
fPMTThreshold
224 G4cout <<
"\tNumber of photons produced by scintillation in this event : "
226 G4cout <<
"\tNumber of photons produced by cerenkov in this event : "
228 G4cout <<
"\tNumber of photons absorbed (OpAbsorption) in this event : "
230 G4cout <<
"\tNumber of photons absorbed at boundaries (OpBoundary) in "
232 G4cout <<
"Unaccounted for photons in this event : "
239 auto run =
static_cast<LXeRun*
>(
240 G4RunManager::GetRunManager()->GetNonConstCurrentRun());
253 G4RunManager::GetRunManager()->rndmSaveThisEvent();