fcl to filter or save choppy events

Author: Mun Jung Jung

sbndcode/Calibration/CMakeLists.txt

@@ -1,2 +1,3 @@
 add_subdirectory(OpReco)
 add_subdirectory(CRT)
+add_subdirectory(DQM)

sbndcode/Calibration/DQM/CMakeLists.txt

@@ -1,28 +1,3 @@
-# cet_build_plugin(FlagChoppy art::module 
-#         lardataobj::RecoBase
-#         ROOT::Core
-#         ROOT::Tree
-#         art::Framework_Core
-#         art::Framework_Principal
-#         art::Utilities
-#         art_root_io::tfile_support
-#         art_root_io::TFileService_service
-#         canvas::canvas
-#         fhiclcpp::fhiclcpp
-#         cetlib::cetlib
-#         cetlib_except::cetlib_except
-#         messagefacility::MF_MessageLogger
-#         lardataobj::RawData
-#         larcore::Geometry_Geometry_service
-#         fhiclcpp::fhiclcpp
-#         cetlib::cetlib
-#         CLHEP::Random
-#         ROOT::Geom
-#         ROOT::XMLIO
-#         ROOT::Gdml
-#         ROOT::FFTW
-# )
-
 cet_build_plugin(FilterChoppy art::module 
         lardataobj::RecoBase
         ROOT::Core

sbndcode/Calibration/DQM/FilterChoppy_module.cc

@@ -98,7 +98,6 @@ namespace filt{
     npeak = 0;
     delta_tmaxmin = 0;
 
-
     // get the timestamps
     art::Handle<std::vector<raw::RDTimeStamp>> timestamp_handle;
     evt.getByLabel("daq", timestamp_handle);
@@ -110,16 +109,14 @@ namespace filt{
     }
     art::fill_ptr_vector(raw_timestamps_handle, timestamp_handle);
 
-
-    // analyze timestamp distributions
-    // collect timestamps
+    // collect TS
     unsigned timeindex = 0;
     for (auto const& timestamps: raw_timestamps_handle) {
       long this_ts = timestamps->GetTimeStamp();
       ts_vec.push_back(this_ts);
       timeindex++;
     }
-    // make binned hist
+    // find min, max TS values
     long ts_min = ts_vec[0];
     long ts_max = ts_vec[0];
     for (long ts : ts_vec) {
@@ -132,39 +129,42 @@ namespace filt{
     }
     delta_tmaxmin = std::abs(ts_max - ts_min);
 
-    float nbin = 100.;
-    float bin_width = (ts_max - ts_min)/nbin;
-    for (int i = 0; i < 100; ++i) {
-      float bin_low = ts_min + i*bin_width;
-      float bin_high = ts_min + (i+1)*bin_width;
-      int bin_count = 0;
-      for (long ts : ts_vec) {
-        if ((bin_low < ts) & (ts <= bin_high)) {
-          bin_count += 1;
-        }
-      }
-      ts_hist_vec.push_back(bin_count);
-    }
-    // find peak bins
-    for (int bc : ts_hist_vec) {
-      //TODO: make threshold fcl parameter
-      if (bc > 1000) {
-        npeak += 1;
-      }
-    }
-
-    // If the energy deposit within the beam time is greater than some limit then trigger the event
-    if (npeak > 1){
-      std::cout << "Event " << event << " has " << npeak << " peaks" << std::endl;
-    }
-    if (delta_tmaxmin > 100000){
-      std::cout << "Event " << event << " has a delta_tmaxmin of " << delta_tmaxmin << std::endl;
-    }
-    //return npeak < 2;
+    // if (delta_tmaxmin > 100000){
+    //   std::cout << "Event " << event << " has a delta_tmaxmin of " << delta_tmaxmin << std::endl;
+    // }
     return delta_tmaxmin < 100000;
 
-  }
 
+    // // Identify choppy events by looking at the # of peaks in the TS distribution
+    // make binned hist
+    // float nbin = 100.;
+    // float bin_width = (ts_max - ts_min)/nbin;
+    // for (int i = 0; i < 100; ++i) {
+    //   float bin_low = ts_min + i*bin_width;
+    //   float bin_high = ts_min + (i+1)*bin_width;
+    //   int bin_count = 0;
+    //   for (long ts : ts_vec) {
+    //     if ((bin_low < ts) & (ts <= bin_high)) {
+    //       bin_count += 1;
+    //     }
+    //   }
+    //   ts_hist_vec.push_back(bin_count);
+    // }
+    // // find peak bins
+    // for (int bc : ts_hist_vec) {
+    //   //TODO: make threshold fcl parameter
+    //   if (bc > 1000) {
+    //     npeak += 1;
+    //   }
+    // }
+
+    // // If the energy deposit within the beam time is greater than some limit then trigger the event
+    // if (npeak > 1){
+    //   std::cout << "Event " << event << " has " << npeak << " peaks" << std::endl;
+    // }
+    // return npeak < 2;
+    
+  }
 
   void FilterChoppy::beginJob() {
   }

sbndcode/Calibration/DQM/filterchoppy.fcl

@@ -31,7 +31,7 @@ outputs:
  out1:
  {
    module_type: RootOutput
-   fileName:    "%ifb_filtered.root"
+   fileName:    "%ifb_choppyfiltered.root"
    dataTier:    "filtered"
    fastCloning: true
    SelectEvents: [ filter ]

sbndcode/Calibration/DQM/savechoppy.fcl

@@ -0,0 +1,55 @@
+#include "services_sbnd.fcl"
+#include "TPCChannelMapService.fcl"
+
+process_name: FilterChoppy
+
+services:
+{
+  TFileService: { fileName: "reco_hist.root"}
+  MemoryTracker:     {}
+  TimeTracker:       {}
+  RandomNumberGenerator: {} 
+  message:              @local::sbnd_message_services_prod_debug
+  FileCatalogMetadata:  @local::art_file_catalog_mc
+  TPCChannelMapService: @local::SBNDTPCChannelMapServiceDefaults
+  AuxDetExptGeoHelperInterface: { service_provider: "sbndcode/CRT/CRTGeometryHelper" }
+  ExptGeoHelperInterface:    @local::sbnd_geometry_helper
+  Geometry:                  @local::sbnd_geo
+  GeometryConfigurationWriter:    {}
+  @table::sbnd_basic_services
+  @table::sbnd_random_services
+}
+
+source:
+{
+  module_type: RootInput
+  maxEvents:  -1
+}
+
+outputs:
+{
+ out1:
+ {
+   module_type: RootOutput
+   fileName:    "%ifb_choppy.root"
+   dataTier:    "filtered"
+   fastCloning: true
+   SelectEvents: [ savechoppy ]
+ }
+}
+
+physics:
+{
+  filters:
+  {
+    filterchoppy: {
+      module_type: FilterChoppy
+    }
+  }
+
+  savechoppy: [ "!filterchoppy"]
+  stream1:  [ out1 ]
+  trigger_paths: [ savechoppy ]
+  end_paths:     [ stream1 ]
+}
+