53 G4Track& photon = *aStep->GetTrack();
56 if (photon.GetDefinition() != G4OpticalPhoton::OpticalPhotonDefinition())
return false;
66 G4ThreeVector localPosition = photon.GetTouchableHandle()->GetHistory()->GetTopTransform().TransformPoint(photon.GetPosition());
67 double xLocal = localPosition.x() *
Unit::mm;
68 double yLocal = localPosition.y() *
Unit::mm;
73 auto* info =
dynamic_cast<Simulation::TrackInfo*
>(photon.GetUserInformation());
75 applyQE = info->getStatus() < 3;
76 fraction = info->getFraction();
79 double qeffi =
m_topgp->getPMTEfficiency(energy, moduleID, pmtID, xLocal, yLocal);
80 if (qeffi == 0 or gRandom->Uniform() * fraction > qeffi) {
81 photon.SetTrackStatus(fStopAndKill);
86 info->setFraction(qeffi);
91 const G4ThreeVector& g_detPoint = photon.GetPosition();
92 const G4ThreeVector& g_detMomDir = photon.GetMomentumDirection();
93 ROOT::Math::XYZPoint detPoint(g_detPoint.x(), g_detPoint.y(), g_detPoint.z());
94 ROOT::Math::XYZVector detMomDir(g_detMomDir.x(), g_detMomDir.y(), g_detMomDir.z());
95 double detTime = photon.GetGlobalTime();
96 double length = photon.GetTrackLength();
99 const G4ThreeVector& g_emiPoint = photon.GetVertexPosition();
100 const G4ThreeVector& g_emiMomDir = photon.GetVertexMomentumDirection();
101 ROOT::Math::XYZPoint emiPoint(g_emiPoint.x(), g_emiPoint.y(), g_emiPoint.z());
102 ROOT::Math::XYZVector emiMomDir(g_emiMomDir.x(), g_emiMomDir.y(), g_emiMomDir.z());
103 double emiTime = photon.GetGlobalTime() - photon.GetLocalTime();
112 auto* simHit =
m_simHits.appendNew(moduleID, pmtID, xLocal, yLocal,
115 int parentID = photon.GetParentID();
116 if (parentID == 0) parentID = photon.GetTrackID();
119 const auto* geo =
m_topgp->getGeometry();
120 if (geo->isModuleIDValid(moduleID)) {
122 const auto&
module = geo->getModule(moduleID);
123 emiPoint =
module.pointToLocal(emiPoint);
124 detPoint =
module.pointToLocal(detPoint);
125 emiMomDir =
module.momentumToLocal(emiMomDir);
126 detMomDir =
module.momentumToLocal(detMomDir);
128 B2ERROR(
"TOP::SensitivePMT: undefined module ID."
129 <<
LogVar(
"moduleID", moduleID));
132 emiPoint, emiMomDir, emiTime,
133 detPoint, detMomDir, detTime,
136 simHit->addRelationTo(simPhoton);
139 photon.SetTrackStatus(fStopAndKill);