Skip to content

Fix bug for opdettype in light propagation time #898

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
asanchezcastillo wants to merge 3 commits into
base: develop
Choose a base branch
from
+6 −6
Open

Fix bug for opdettype in light propagation time #898

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
6686745
Fix bug for opdettype in light propagation time
asanchezcastillo Jan 20, 2026
306b771
Add continue for non-PMT channels
asanchezcastillo Jan 20, 2026
c461ba5
Remove vector and use XAs
asanchezcastillo Jan 26, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
10 changes: 5 additions & 5 deletions sbndcode/LightPropagationCorrection/LightPropagationCorrection_module.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -442,16 +442,16 @@ void sbnd::LightPropagationCorrection::GetPropagationTimeCorrectionPerChannel()
double dy = fSpacePointY[sp] - _opDetY;
double dz = fSpacePointZ[sp] - _opDetZ;
double distanceToOpDet = std::sqrt(dx*dx + dy*dy + dz*dz);
//double spToCathode = abs(fSpacePointX[sp]); // Distance from space point to cathode in mm
//double cathodeToOpDet = std::sqrt(_opDetX*_opDetX + dy*dy + dz*dz); // Distance from cathode to OpDet in mm
//float lightPropTimeVIS = spToCathode/fVGroupVUV + cathodeToOpDet/fVGroupVIS; // Speed

double cathodeToOpDet = std::sqrt(_opDetX*_opDetX + (dy/2)*(dy/2) + (dz/2)*(dz/2)); // Distance from cathode to OpDet in mm
double spToCathode = std::sqrt( fSpacePointX[sp]*fSpacePointX[sp] + (dy/2)*(dy/2) + (dz/2)*(dz/2)); // Distance from space point to cathode in mm

float lightPropTimeVIS = spToCathode/fVGroupVUV + cathodeToOpDet/fVGroupVIS; // Speed
float lightPropTimeVUV = distanceToOpDet / fVGroupVUV; // Speed of light in mm/ns for VUV
float lightPropTime = std::min(lightPropTimeVIS, lightPropTimeVUV);
float lightPropTime = 0;
if(fPDSMap.pdType(opch)=="pmt_coated" || fPDSMap.pdType(opch)=="xarapuca_vuv")
lightPropTime = std::min(lightPropTimeVIS, lightPropTimeVUV);
else if(fPDSMap.pdType(opch)=="pmt_uncoated" || fPDSMap.pdType(opch)=="xarapuca_vis")
lightPropTime = lightPropTimeVIS;
float partPropTime = std::sqrt((fSpacePointX[sp]-fRecoVx)*(fSpacePointX[sp]-fRecoVx) + (fSpacePointY[sp]-fRecoVy)*(fSpacePointY[sp]-fRecoVy) + (fSpacePointZ[sp]-fRecoVz)*(fSpacePointZ[sp]-fRecoVz))/fSpeedOfLight;
float PropTime = lightPropTime + partPropTime;
if(PropTime < minPropTime) minPropTime = PropTime;
Expand Down
2 changes: 1 addition & 1 deletion sbndcode/LightPropagationCorrection/LightPropagationCorrection_module.hh
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -145,7 +145,7 @@ private:


geo::WireReadoutGeom const& fWireReadout = art::ServiceHandle<geo::WireReadout>()->Get();

opdet::sbndPDMapAlg fPDSMap;
//Flash finder manager
::lightana::FlashFinderManager _mgr;
::lightana::FlashFinderManager _mgr_tpc0;
Expand Down