Added a CRT Veto Producer module. Stores boolean values based on CRT hit information

Author: aantonakis

sbndcode/CRT/CMakeLists.txt

@@ -2,6 +2,7 @@ add_subdirectory(CRTAna)
 add_subdirectory(CRTBackTracker)
 add_subdirectory(CRTEventDisplay)
 add_subdirectory(CRTReco)
+add_subdirectory(CRTVeto)
 add_subdirectory(CRTSimulation)
 add_subdirectory(CRTTPCMatching)
 add_subdirectory(CRTUtils)

sbndcode/CRT/CRTReco/CMakeLists.txt

@@ -34,11 +34,5 @@ simple_plugin(
               Eigen3::Eigen
 )
 
-simple_plugin(
-              CRTVetoProducer module
-              sbnobj::SBND_CRT
-              sbndcode_CRTUtils
-              sbndcode_GeoWrappers
-)
 
 install_fhicl()

sbndcode/CRT/CRTReco/crtrecoproducers_sbnd.fcl

@@ -83,17 +83,4 @@ crttrackproducer_data_sbnd.CRTGeoAlg:     @local::crtgeoalg_data_sbnd
 crttrackproducer_data_sbnd.UseTs0:        true
 crttrackproducer_data_sbnd.MaskedTaggers: [ ]
 
-crtvetoproducer_sbnd:
-{
-   CRTGeoAlg:                  @local::crtgeoalg_sbnd
-   CRTSpacePointModuleLabel:   "crtspacepoints"
-   WindowStart:                0.4
-   WindowEnd:                  2.0
-   module_type:                "CRTVetoProducer"
-}
-
-crtvetoproducer_data_sbnd: @local::crtvetoproducer_sbnd
-crtvetoproducer_data_sbnd.CRTGeoAlg:     @local::crtgeoalg_data_sbnd
-
-
 END_PROLOG

sbndcode/CRT/CRTVeto/CMakeLists.txt

@@ -0,0 +1,14 @@
+art_make_library(
+                 LIBRARIES
+                        sbnobj::SBND_CRT
+                        sbndcode_GeoWrappers
+)
+
+simple_plugin(
+              CRTVetoProducer module
+              sbnobj::SBND_CRT
+              sbndcode_CRTUtils
+              sbndcode_GeoWrappers
+)
+
+install_fhicl()

…de/CRT/CRTReco/CRTVetoProducer_module.cc …de/CRT/CRTVeto/CRTVetoProducer_module.ccsbndcode/CRT/CRTReco/CRTVetoProducer_module.cc renamed to sbndcode/CRT/CRTVeto/CRTVetoProducer_module.cc sbndcode/CRT/CRTReco/CRTVetoProducer_module.cc renamed to sbndcode/CRT/CRTVeto/CRTVetoProducer_module.cc

@@ -108,7 +108,21 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
     double x = CRTSpacePointVec[i]->X();
     double y = CRTSpacePointVec[i]->Y();
     double z = CRTSpacePointVec[i]->Z();
-    double t = CRTSpacePointVec[i]->Ts0();
+    double t = CRTSpacePointVec[i]->Ts0()/1000; // Convert to us
+    
+    // Debugging statement 
+    
+    /*
+    if ( (t <= 2) && (t >= 0.4) ) {
+      std::cout << std::endl;
+      std::cout << std::endl;
+      std::cout << "Found a CRT Hit within the South Wall Top Hat Window!" << std::endl;
+      std::cout << std::endl;
+      std::cout << std::endl;
+
+    } 
+    */
+
     // check if in the beam window --> Make configureable
     if ((t > fWindowEnd) || (t < fWindowStart)) {
       continue;
@@ -148,6 +162,7 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
   bool v1 = false;
   bool v2 = false;
   bool v3 = false;
+  bool v4 = false;
 
   // V0
   if ((Nwalls_all + Ntop) > 0) {
@@ -165,9 +180,12 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
   if ( (Nwalls_noNorth > 0) || ((TL > 0) && (TH > 0)) ) {
     v3 = true;
   }
+  if (S > 0)  {
+    v4 = true;
+  }
 
   //crtveto =  CRTVeto(v0, v1, v2, v3); 
-  auto crtveto            = std::make_unique<CRTVeto>(v0, v1, v2, v3);
+  auto crtveto            = std::make_unique<CRTVeto>(v0, v1, v2, v3, v4);
 
   e.put(std::move(crtveto));
   //e.put(std::move(vetoSpacePointAssn));

sbndcode/CRT/CRTVeto/crtvetoproducer_sbnd.fcl

@@ -0,0 +1,20 @@
+#include "crtgeoalg_sbnd.fcl"
+#include "crtsimmodules_sbnd.fcl"
+
+BEGIN_PROLOG
+
+
+crtvetoproducer_sbnd:
+{
+   CRTGeoAlg:                  @local::crtgeoalg_sbnd
+   CRTSpacePointModuleLabel:   "crtspacepoints"
+   WindowStart:                0.4
+   WindowEnd:                  2.0
+   module_type:                "CRTVetoProducer"
+}
+
+crtvetoproducer_data_sbnd: @local::crtvetoproducer_sbnd
+crtvetoproducer_data_sbnd.CRTGeoAlg:     @local::crtgeoalg_data_sbnd
+
+
+END_PROLOG

sbndcode/CRT/CRTReco/run_crtveto.fcl sbndcode/CRT/CRTVeto/run_crtveto.fclsbndcode/CRT/CRTReco/run_crtveto.fcl renamed to sbndcode/CRT/CRTVeto/run_crtveto.fcl

@@ -1,6 +1,6 @@
 #include "services_sbnd.fcl"
 #include "rootoutput_sbnd.fcl"
-#include "crtrecoproducers_sbnd.fcl"
+#include "crtvetoproducer_sbnd.fcl"
 
 process_name: CRTVeto
 
@@ -9,6 +9,8 @@ services:
   @table::sbnd_services
 }
 
+services.DetectorClocksService.InheritClockConfig: false
+
 source:
 {
   module_type: RootInput