Fix bug for opdettype in light propagation time

Author: asanchezcastillo

sbndcode/LightPropagationCorrection/LightPropagationCorrection_module.cc

@@ -36,6 +36,11 @@ sbnd::LightPropagationCorrection::LightPropagationCorrection(fhicl::ParameterSet
         fOpDetX.push_back(pdCenter.X());
         fOpDetY.push_back(pdCenter.Y());
         fOpDetZ.push_back(pdCenter.Z());
+        if(fPDSMap.pdType(opch)=="pmt_coated") fOpDetType.push_back(0);
+        else if(fPDSMap.pdType(opch)=="pmt_uncoated") fOpDetType.push_back(1);
+        else if(fPDSMap.pdType(opch)=="xarapuca_vuv") fOpDetType.push_back(2);
+        else if(fPDSMap.pdType(opch)=="xarapuca_vis") fOpDetType.push_back(3);
+        else fOpDetType.push_back(-1);
     }
 
     auto const& tpc = art::ServiceHandle<geo::Geometry>()->TPC();
@@ -442,16 +447,18 @@ 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(fOpDetType[opdet]==0)
+                lightPropTime = std::min(lightPropTimeVIS, lightPropTimeVUV);
+            else if(fOpDetType[opdet]==1)
+                lightPropTime = lightPropTimeVIS;
+            else
+                throw art::Exception(art::errors::LogicError) << " OpDet Type " <<  fOpDetType[opdet] << " not supported ." << std::endl;
             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;

sbndcode/LightPropagationCorrection/LightPropagationCorrection_module.hh

@@ -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;
@@ -162,6 +162,7 @@ private:
     std::vector<double> fOpDetX;
     std::vector<double> fOpDetY;
     std::vector<double> fOpDetZ;
+    std::vector<int> fOpDetType;
 
     std::string fReco2Label;
     std::string fOpT0FinderModuleLabel;