Fixed Timing bug for RWM referencing in BNBLight

Author: aantonakis

sbndcode/CRT/CRTVeto/CRTVetoProducer_module.cc

@@ -109,7 +109,7 @@ sbnd::crt::CRTVetoProducer::CRTVetoProducer(fhicl::ParameterSet const& p)
   , fCRTClusterModuleLabel(p.get<std::string>("CRTClusterModuleLabel"))
   , fCRTSpacePointModuleLabel(p.get<std::string>("CRTSpacePointModuleLabel"))
   , fCRTTimingReferenceInfoLabel(p.get<std::string>("CRTTimingReferenceInfoLabel"))
-  , fSPECTDCModuleLabel(p.get<std::string>("SPECTDCModuleLabel"))
+  , fSPECTDCModuleLabel(p.get<std::string>("SPECTDCModuleLabel", ""))
   , fIsData(p.get<bool>("IsData"))
   , fDAQHeaderModuleLabel(p.get<std::string>("DAQHeaderModuleLabel", ""))
   , fDAQHeaderInstanceLabel(p.get<std::string>("DAQHeaderInstanceLabel", ""))
@@ -137,8 +137,9 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
   std::cout << std::endl;
 
   // May only be needed for BNB+Light Data --> TODO
-  uint64_t raw_ts = 0, ref_time = 0;
-  uint32_t ref_time_ns = 0;
+  uint64_t raw_ts = 0;
+  int64_t ref_time = 0; 
+  int64_t ref_time_ns = 0;
 
   // Data needs to be handled differently than MC
 
@@ -197,7 +198,7 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
       std::cout << std::endl;
       std::cout << std::endl;
     }
-    ref_time_ns = ref_time % static_cast<uint64_t>(1e9);
+    ref_time_ns = ref_time % static_cast<int64_t>(1e9);
 
   } // end of check for Data
 
@@ -278,7 +279,10 @@ void sbnd::crt::CRTVetoProducer::produce(art::Event& e)
       }
     } // end loop over space points
   } // end of loop over clusters
- 
+  
+  std::cout << "Finished Loop over Clusters --> Produce CRTVeto products" << std::endl;  
+  std::cout << std::endl; 
+
   int Nwalls_all = S + N + E + W;
   int Nwalls_noNorth = S + E + W;
   int Ntop = TL + TH;

sbndcode/CRT/CRTVeto/crtvetoproducer_sbnd.fcl

@@ -23,13 +23,27 @@ crtvetoproducer_data_sbnd.IsData:                    true
 #crtvetoproducer_data_sbnd.WindowEnd:                 2.065
 
 # Values for testing
-crtvetoproducer_data_sbnd.WindowStart:               -2 
-crtvetoproducer_data_sbnd.WindowEnd:                 2
+crtvetoproducer_data_sbnd.WindowStart:               -1.135 
+crtvetoproducer_data_sbnd.WindowEnd:                 0.565
 crtvetoproducer_data_sbnd.SPECTDCModuleLabel:        "tdcdecoder"
 
 crtvetoproducer_data_sbnd.DAQHeaderModuleLabel:      "daq"
 crtvetoproducer_data_sbnd.DAQHeaderInstanceLabel:    "RawEventHeader"
 crtvetoproducer_data_sbnd.RawTSCorrection:           367000
 crtvetoproducer_data_sbnd.MaxAllowedRefTimeDiff:     3000000
 
+# Top Hat Configuration for Debugging
+crtvetoproducer_data_sbnd_tophat:                    @local::crtvetoproducer_sbnd
+crtvetoproducer_data_sbnd_tophat.CRTGeoAlg:          @local::crtgeoalg_data_sbnd
+crtvetoproducer_data_sbnd_tophat.IsData:             true
+
+# Values for testing
+crtvetoproducer_data_sbnd_tophat.SPECTDCModuleLabel:        "tdcdecoder"
+crtvetoproducer_data_sbnd_tophat.DAQHeaderModuleLabel:      "daq"
+crtvetoproducer_data_sbnd_tophat.DAQHeaderInstanceLabel:    "RawEventHeader"
+crtvetoproducer_data_sbnd_tophat.RawTSCorrection:           367000
+crtvetoproducer_data_sbnd_tophat.MaxAllowedRefTimeDiff:     3000000
+crtvetoproducer_data_sbnd_tophat.WindowStart:               -2 
+crtvetoproducer_data_sbnd_tophat.WindowEnd:                 2
+
 END_PROLOG

sbndcode/CRT/CRTVeto/run_crtveto_data_tophat.fcl

@@ -0,0 +1,10 @@
+#include "crt_channel_map_service.fcl"
+#include "crt_calib_service.fcl"
+#include "run_crtveto.fcl"
+
+services.message:              @local::sbnd_message_services
+services.CRTChannelMapService: @local::crt_channel_map_standard
+services.CRTCalibService:      @local::crt_calib_service
+
+physics.producers.crtveto:     @local::crtvetoproducer_data_sbnd_tophat
+