Add event displays for questionable cases

Author: henrylay97

sbncode/Calibration/CMakeLists.txt

@@ -25,6 +25,7 @@ cet_build_plugin( TrackCaloSkimmer art::module
                             sbncode_CAFMaker
                             sbnobj::Common_Calibration_dict
                             larevt::SpaceCharge
+			    sbndcode::CRT_CRTEventDisplay
 )
 
 cet_build_plugin(TrackCaloSkimmerSelectStoppingTrack art::tool

sbncode/Calibration/TrackCaloSkimmer.h

@@ -66,6 +66,8 @@
 #include "sbnobj/SBND/CRT/CRTTrack.hh"
 #include "sbnobj/SBND/CRT/CRTSpacePoint.hh"
 
+#include "sbndcode/CRT/CRTEventDisplay/CRTEventDisplayAlg.h"
+
 #include "ITCSSelectionTool.h"
 
 namespace sbn {
@@ -240,5 +242,7 @@ class sbn::TrackCaloSkimmer : public art::EDAnalyzer {
   // Fitting info
   TFitter fFitExp;
   TFitter fFitConst;
+
+  sbnd::crt::CRTEventDisplayAlg fCRTEventDisplayAlg;
 };
 #endif

sbncode/Calibration/TrackCaloSkimmer_module.cc

@@ -64,6 +64,7 @@ sbn::TrackCaloSkimmer::TrackCaloSkimmer(fhicl::ParameterSet const& p)
   : EDAnalyzer{p},
     fFitExp(2),
     fFitConst(1)
+  , fCRTEventDisplayAlg(p.get<fhicl::ParameterSet>("CRTEventDisplayAlg"))
 {
   // Grab config
   fPFPproducer  = p.get< art::InputTag > ("PFPproducer","pandoraGausCryo0");
@@ -260,6 +261,7 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
 
   // The raw digit list is not sorted, so make it into a map on the WireID
   std::map<geo::WireID, art::Ptr<raw::RawDigit>> rawdigits;
+  /*
   for (const art::Ptr<raw::RawDigit> &d: rawdigitlist) {
     
     std::vector<geo::WireID> wids;
@@ -279,7 +281,7 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
 
     rawdigits[wids[0]] = d;
   } 
-
+  */
   // Collect all hits
   art::ValidHandle<std::vector<recob::Hit>> allhit_handle = e.getValidHandle<std::vector<recob::Hit>>(fHITproducer);
   std::vector<art::Ptr<recob::Hit>> allHits;
@@ -296,10 +298,12 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
   const cheat::BackTrackerService *bt = NULL;
 
   if (simchannels.size()) {
+    /*
     art::ServiceHandle<cheat::BackTrackerService> bt_serv;
     id_to_ide_map = caf::PrepSimChannels(simchannels, *wireReadout);
     id_to_truehit_map = caf::PrepTrueHits(allHits, clock_data, *bt_serv.get());
     bt = bt_serv.get();
+    */
   }
 
   // service data
@@ -408,6 +412,23 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
 
     hasT0 = hasPFPT0 || hasCRTTrackT0 || hasCRTHitT0 || hasCRTSpacePointT0;
 
+    if(hasPFPT0)
+      {
+	fCRTEventDisplayAlg.SetPrint(true);
+	fCRTEventDisplayAlg.SetMinTime(t0PFP - 3e3);
+	fCRTEventDisplayAlg.SetMaxTime(t0PFP + 3e3);
+
+	auto const clockData = art::ServiceHandle<detinfo::DetectorClocksService const>()->DataFor(e);
+
+	if(hasCRTTrackT0)
+	  {
+	    std::cout << t0CRTTrack << std::endl;
+	    fCRTEventDisplayAlg.Draw(clockData, e, Form("/exp/sbnd/data/users/hlay/crt_fall_production_2025/t0_event_displays/good/crtEventDisplayRun%iSubrun%iEvent%iPFP%li", e.run(), e.subRun(), e.event(), p_pfp.key()));
+	  }
+	else
+	  fCRTEventDisplayAlg.Draw(clockData, e, Form("/exp/sbnd/data/users/hlay/crt_fall_production_2025/t0_event_displays/bad/crtEventDisplayRun%iSubrun%iEvent%iPFP%li", e.run(), e.subRun(), e.event(), p_pfp.key()));
+      }
+    continue;
     // "whicht0" should reflect the T0 used for the reconstruction of the drift coordinate.
     if(!hasT0) whicht0 = -1;
     // In this way, if a track is T0 tagged from PFP and CRT tagged, which T0 reflects the PFP Tag.

sbncode/Calibration/fcl/sbnd/sbnd_trackcalo_skimmer.fcl

@@ -1,3 +1,5 @@
+#include "crteventdisplayalg_sbnd.fcl"
+
 BEGIN_PROLOG
 
 caloskim_calorimetry: {
@@ -56,7 +58,7 @@ stopping_selection: {
   EndMediandQdxCut: 1600
   NumberTimeSamples: 3400
   MediandQdxRRMax: 5.
-  CheckFiducialX: true
+  CheckFiducialX: false
 }
 
 exiting_selection: {
@@ -74,7 +76,7 @@ exiting_selection: {
   EndMediandQdxCut: -1
   RequireDownwards: false
   NumberTimeSamples: 3400
-  CheckFiducialX: true
+  CheckFiducialX: false
 }
 
 a2c_selection: {
@@ -90,10 +92,11 @@ throughgoing_selection: {
   FVInsetMaxY: 5
   FVInsetMinZ: 5
   FVInsetMaxZ: 5
-  CheckFiducialX: true
+  CheckFiducialX: false
 }
 
 caloskim_nodigits_goldentracks: @local::caloskim_nodigits
 caloskim_nodigits_goldentracks.SelectionTools: [@local::stopping_selection, @local::a2c_selection, @local::throughgoing_selection, @local::exiting_selection]
+caloskim_nodigits_goldentracks.CRTEventDisplayAlg: @local::crteventdisplayalg_sbnd_temp
 
 END_PROLOG