Merge branch 'develop' into feature/munjung-calohit_tmatch

Author: bear-is-asleep

CMakeLists.txt

@@ -15,12 +15,12 @@
 
 cmake_minimum_required(VERSION 3.20 FATAL_ERROR)
 
-set(${PROJECT_NAME}_CMAKE_PROJECT_VERSION_STRING 09.93.01.02.01)
+set(${PROJECT_NAME}_CMAKE_PROJECT_VERSION_STRING 10.03.01)
 find_package(cetmodules REQUIRED)
 project(sbndcode LANGUAGES CXX)
 
 # for CI:
-# project(sbndcode VERSION 09.93.01.02.01)
+# project(sbndcode VERSION 10.03.01)
 
 message(STATUS "\n")
 message(STATUS "=================================   ${PROJECT_NAME}   =================================")
@@ -80,6 +80,7 @@ find_package( larcore REQUIRED )
 find_package( larcorealg REQUIRED )
 find_package( larcoreobj REQUIRED )
 find_package( larreco REQUIRED )
+find_package( larrecodnn REQUIRED )
 find_package( larana REQUIRED )
 find_package( PostgreSQL REQUIRED )
 find_package( hep_concurrency REQUIRED )

sbndcode/AnalysisTree/AnalysisTree_module.cc

@@ -3677,17 +3677,15 @@ bool sbnd::AnalysisTree::TrackIdToMCTruth( Int_t const trkID, art::Ptr<simb::MCT
 double sbnd::AnalysisTree::bdist(const recob::Track::Point_t& pos)
 {
   // Get geometry.
-  art::ServiceHandle<geo::Geometry> geom;
-
+  geo::TPCGeo const& tpc = art::ServiceHandle<geo::Geometry>()->TPC({0, 0});
   double d1 = pos.X();                             // Distance to right side (wires).
-  double d2 = 2.*geom->DetHalfWidth() - pos.X();   // Distance to left side (cathode).
-  double d3 = pos.Y() + geom->DetHalfHeight();     // Distance to bottom.
-  double d4 = geom->DetHalfHeight() - pos.Y();     // Distance to top.
+  double d2 = 2.*tpc.HalfWidth() - pos.X();   // Distance to left side (cathode).
+  double d3 = pos.Y() + tpc.HalfHeight();     // Distance to bottom.
+  double d4 = tpc.HalfHeight() - pos.Y();     // Distance to top.
   double d5 = pos.Z();                             // Distance to front.
-  double d6 = geom->DetLength() - pos.Z();           // Distance to back.
+  double d6 = tpc.Length() - pos.Z();           // Distance to back.
 
-  double result = std::min(std::min(std::min(std::min(std::min(d1, d2), d3), d4), d5), d6);
-  return result;
+  return std::min({d1, d2, d3, d4, d5, d6});
 }
 
 
@@ -3701,16 +3699,14 @@ double sbnd::AnalysisTree::length(const recob::Track& track)
 // Length of MC particle, trajectory by trajectory.
 double sbnd::AnalysisTree::length(const simb::MCParticle& part, TVector3& start, TVector3& end)
 {
-  // Get geometry.
-  art::ServiceHandle<geo::Geometry> geom;
-
   // Get active volume boundary.
-  double xmin = -2.0 * geom->DetHalfWidth() - 1e-8;
-  double xmax = 2.0 * geom->DetHalfWidth() + 1e-8;
-  double ymin = -geom->DetHalfHeight() -1e-8;
-  double ymax = geom->DetHalfHeight() + 1e-8;
+  geo::TPCGeo const& tpc = art::ServiceHandle<geo::Geometry>()->TPC({0, 0});
+  double xmin = -2.0 * tpc.HalfWidth() - 1e-8;
+  double xmax = 2.0 * tpc.HalfWidth() + 1e-8;
+  double ymin = -tpc.HalfHeight() -1e-8;
+  double ymax = tpc.HalfHeight() + 1e-8;
   double zmin = 0. -1e-8;
-  double zmax = geom->DetLength() + 1e-8;
+  double zmax = tpc.Length() + 1e-8;
 
   // Get number traj points
   int n = part.NumberTrajectoryPoints();

sbndcode/AnalysisTree/CMakeLists.txt

@@ -13,7 +13,6 @@ set(      MODULE_LIBRARIES larcorealg::Geometry
                            lardata::RecoObjects
                            larpandora::LArPandoraInterface
                            sbndcode_CRTUtils
-                           sbndcode_CRT
                            sbnobj::Common_CRT
                            nusimdata::SimulationBase
                            art::Framework_Core

sbndcode/CMakeLists.txt

@@ -42,3 +42,7 @@ add_subdirectory(WireCell)
 
 add_subdirectory(Decoders)
 add_subdirectory(ChannelMaps)
+
+# for CVN
+# Disable for now
+#add_subdirectory(SBNDCVN)

sbndcode/CRT/CMakeLists.txt

@@ -6,48 +6,37 @@ add_subdirectory(CRTSimulation)
 add_subdirectory(CRTTPCMatching)
 add_subdirectory(CRTUtils)
 
-art_make(
-         EXCLUDE
-                CRTChannelMapAlg.cxx
-                CRTGeometryHelper_service.cc
-         LIBRARY_NAME
-                sbndcode_CRTData
-         LIB_LIBRARIES
-                larcorealg::Geometry
-                larcore::Geometry_Geometry_service
-                lardataobj::RecoBase
-                art::Framework_Core
-                art::Framework_IO_Sources
-                art::Framework_Principal
-                canvas::canvas
-                art::Framework_Services_Registry
-                messagefacility::MF_MessageLogger
-
-                fhiclcpp::fhiclcpp
-                ROOT::Core
-                ROOT::Geom
-                cetlib::cetlib
-)
-
 art_make_library(
-                 SOURCE       CRTChannelMapAlg.cxx
-                 LIBRARIES
-                        larcorealg::Geometry
-                        sbndcode_CRTData
-                        cetlib_except::cetlib_except
-                        messagefacility::MF_MessageLogger
-                        fhiclcpp::fhiclcpp
-                        ROOT::Core
-                        ROOT::Geom
-                        cetlib::cetlib
+    SOURCE
+        CRTGeoObjectSorter.cxx
+    LIBRARY_NAME
+        sbndcode_CRTData
+    LIBRARIES
+        larcorealg::Geometry
+        larcore::Geometry_Geometry_service
+        lardataobj::RecoBase
+        art::Framework_Core
+        art::Framework_IO_Sources
+        art::Framework_Principal
+        canvas::canvas
+        art::Framework_Services_Registry
+        messagefacility::MF_MessageLogger
+        fhiclcpp::fhiclcpp
+        ROOT::Core
+        ROOT::Geom
+        cetlib::cetlib
 )
-simple_plugin(
-              CRTGeometryHelper service
-              sbndcode_CRT
-              art::Framework_Core
-              art::Framework_IO_Sources
-              art::Framework_Principal
-              canvas::canvas
-              art::Framework_Services_Registry
-              messagefacility::MF_MessageLogger
+
+cet_build_plugin( CRTAuxDetInitializer art::tool
+    LIBRARIES REG
+        larcorealg::Geometry
+        sbndcode_CRTData
+        cetlib_except::cetlib_except
+        messagefacility::MF_MessageLogger
+        fhiclcpp::fhiclcpp
+        cetlib::cetlib
 )
+
+install_headers()
+install_fhicl()
+install_source()

sbndcode/CRT/CRTAna/CRTAnalysis_module.cc

@@ -39,6 +39,7 @@
 #include "sbndcode/CRT/CRTBackTracker/CRTBackTrackerAlg.h"
 #include "sbndcode/CRT/CRTUtils/CRTCommonUtils.h"
 #include "sbndcode/Decoders/PTB/sbndptb.h"
+#include "sbndcode/Timing/SBNDRawTimingObj.h"
 
 namespace sbnd::crt {
   class CRTAnalysis;
@@ -94,7 +95,7 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
   std::string fMCParticleModuleLabel, fSimDepositModuleLabel, fFEBDataModuleLabel, fCRTStripHitModuleLabel,
     fCRTClusterModuleLabel, fCRTSpacePointModuleLabel, fCRTTrackModuleLabel, fTPCTrackModuleLabel,
     fCRTSpacePointMatchingModuleLabel, fCRTTrackMatchingModuleLabel, fPFPModuleLabel, fPTBModuleLabel,
-    fTDCModuleLabel;
+    fTDCModuleLabel, fTimingReferenceModuleLabel;
   bool fDebug, fDataMode, fNoTPC, fHasPTB, fHasTDC;
 
   TTree* fTree;
@@ -104,6 +105,8 @@ class sbnd::crt::CRTAnalysis : public art::EDAnalyzer {
   int _run;
   int _subrun;
   int _event;
+  int _crt_timing_reference_type;
+  int _crt_timing_reference_channel;
 
   //mc truth
   std::vector<int16_t>              _mc_trackid;
@@ -331,6 +334,7 @@ sbnd::crt::CRTAnalysis::CRTAnalysis(fhicl::ParameterSet const& p)
   fPFPModuleLabel                   = p.get<std::string>("PFPModuleLabel", "pandora");
   fPTBModuleLabel                   = p.get<std::string>("PTBModuleLabel", "ptbdecoder");
   fTDCModuleLabel                   = p.get<std::string>("TDCModuleLabel", "tdcdecoder");
+  fTimingReferenceModuleLabel       = p.get<std::string>("TimingReferenceModuleLabel", "crtstrips");
   fDebug                            = p.get<bool>("Debug", false);
   fDataMode                         = p.get<bool>("DataMode", false);
   fNoTPC                            = p.get<bool>("NoTPC", false);
@@ -346,6 +350,8 @@ sbnd::crt::CRTAnalysis::CRTAnalysis(fhicl::ParameterSet const& p)
   fTree->Branch("run", &_run);
   fTree->Branch("subrun", &_subrun);
   fTree->Branch("event", &_event);
+  fTree->Branch("crt_timing_reference_type", &_crt_timing_reference_type);
+  fTree->Branch("crt_timing_reference_channel", &_crt_timing_reference_channel);
 
   if(!fDataMode)
     {
@@ -621,6 +627,17 @@ void sbnd::crt::CRTAnalysis::analyze(art::Event const& e)
 
   if(fDebug) std::cout << "This is event " << _run << "-" << _subrun << "-" << _event << std::endl;
 
+  _crt_timing_reference_type    = -1;
+  _crt_timing_reference_channel = -1;
+
+  art::Handle<raw::TimingReferenceInfo> TimingReferenceHandle;
+  e.getByLabel(fTimingReferenceModuleLabel, TimingReferenceHandle);
+  if(TimingReferenceHandle.isValid())
+    {
+      _crt_timing_reference_type    = TimingReferenceHandle->timingType;
+      _crt_timing_reference_channel = TimingReferenceHandle->timingChannel;
+    }
+
   if(fHasPTB)
     {
       // Get PTBs

sbndcode/CRT/CRTAuxDetInitializer_tool.cc

@@ -0,0 +1,59 @@
+#include "larcorealg/Geometry/AuxDetGeo.h"
+#include "larcorealg/Geometry/AuxDetSensitiveGeo.h"
+#include "larcorealg/Geometry/AuxDetGeometryCore.h"
+
+#include "art/Utilities/ToolMacros.h"
+
+namespace {
+  class CRTAuxDetInitializer : public geo::AuxDetInitializer {
+  public:
+    explicit CRTAuxDetInitializer(fhicl::ParameterSet const&) {}
+
+  private:
+    geo::AuxDetReadoutInitializers
+    initialize(std::vector<geo::AuxDetGeo> const& adgeo) const override
+    {
+      geo::AuxDetReadoutInitializers result;
+      // Map the AuxDetGeo names to their position in the sorted vector
+      //
+      // Each tagger is composed of scintillator modules, composed of 16 strips.
+      // In the geometry, CRTStripArrays are AuxDets and CRTStrips are the
+      // AuxDetSensitives. Each strip has two SiPM channels, one per optical
+      // fiber (not in the geometry).
+      //
+
+      for (size_t a=0; a<adgeo.size(); a++){
+        std::string volName(adgeo[a].TotalVolume()->GetName());
+
+        long unsigned int number_scintillating_strips = 0;
+
+        if (strncmp(((adgeo[a].TotalVolume())->GetShape())->GetName(), "CRTstripMINOSArray", 18) == 0) {
+          number_scintillating_strips = 20;    //To account for the MINOS modules.
+        }
+        else {number_scintillating_strips = 16;}
+
+        size_t nsv = adgeo[a].NSensitiveVolume();
+        if (nsv != number_scintillating_strips) {
+          throw cet::exception("CRTChannelMap")
+            << "Wrong number of sensitive volumes for CRT volume "
+            << volName << " (got " << nsv << ", expected 16)" << std::endl;
+        }
+
+        result.ADGeoToName[a] = volName;
+        result.NameToADGeo[volName] = a;
+
+        if (volName.find("CRTStripArray") != std::string::npos) {
+          for (size_t svID=0; svID<number_scintillating_strips; svID++) {
+            for (size_t ich=0; ich<2; ich++) {
+              size_t chID = 2 * svID + ich;
+              result.ADGeoToChannelAndSV[a].push_back(std::make_pair(chID, svID));
+            }
+          }
+        }
+      }
+      return result;
+    }
+  };
+}
+
+DEFINE_ART_CLASS_TOOL(CRTAuxDetInitializer)

sbndcode/CRT/CRTBackTracker/CMakeLists.txt

@@ -2,7 +2,7 @@ art_make_library(
                  SOURCE       CRTBackTrackerAlg.cc
                  LIBRARIES
                         sbnobj::SBND_CRT
-                        sbndcode_GeoWrappers
+                        sbndcode::GeoWrappers
                         larsim::MCCheater_ParticleInventoryService_service
 )