60 #include "G4HCofThisEvent.hh"
61 #include "G4VHitsCollection.hh"
77 : runAct(0),genAction(0),hitsfile(0),pmtfile(0)
219 G4cout <<
" Number of hits in LXe: " << S_hits <<
G4endl;
254 G4cout <<
" Total energy in LXe: "
275 G4cout <<
" Average light collection time: "
277 G4cout <<
" Number of PMT hits (photoelectron equivalent): "
323 std::stringstream
sss;
325 filename = sss.str();
331 (*hitsfile) <<
"Evt Eprim Etot LXe LXeTime PMT PMTTime Seed1 Seed2 First Flags"
333 (*hitsfile) <<
"# MeV MeV hits ns hits ns hit"
343 (*hitsfile) << std::setiosflags(std::ios::fixed)
344 << std::setprecision(4)
351 << std::setiosflags(std::ios::scientific)
352 << std::setprecision(2)
355 << std::setiosflags(std::ios::fixed)
356 << std::setprecision(4)
359 << *(
seeds+1) <<
"\t"
369 G4cout <<
" Event summary in file " << filename <<
G4endl;
373 G4int firstparticleIndex = 0;
379 firstparticleIndex = 5;
391 long seed2 = *(
seeds+1);
432 std::stringstream
sss;
434 filename = sss.str();
446 (*pmtfile) <<
"Hit# X, mm Y, mm Z, mm" <<
G4endl;
447 (*pmtfile) << std::setiosflags(std::ios::fixed)
448 << std::setprecision(3)
453 x = ((*hits)[i]->GetPos()).
x()/
mm;
454 y = ((*hits)[i]->GetPos()).
y()/
mm;
455 z = ((*hits)[i]->GetPos()).
z()/
mm;
456 (*pmtfile) << i <<
"\t"
476 G4cout <<
" " << P_hits <<
" PMT hits in " << filename <<
G4endl;
489 G4int n_trajectories = 0;
491 if(trajContainer) n_trajectories = trajContainer->
entries();
492 for (
G4int i=0; i<n_trajectories; i++) {