Merge pull request #647 from SBNSoftware/feature/lynnt_opt0finder_data

Add fcls to run OpT0Finder on data

Author: bear-is-asleep

sbndcode/JobConfigurations/standard/reco/reco2_data.fcl

@@ -1,5 +1,6 @@
 #include "crt_channel_map_service.fcl"
 #include "crt_calib_service.fcl"
+#include "opt0finder_sbnd_data.fcl"
 #include "standard_reco2_sbnd.fcl"
 
 services:
@@ -16,10 +17,11 @@ physics.producers:
     crtclustering:  @local::crtclusterproducer_data_sbnd
     crtspacepoints: @local::crtspacepointproducer_data_sbnd
     crttracks:      @local::crttrackproducer_data_sbnd
+    opt0finder:     @local::sbnd_opt0_finder_data
 }
 
 physics.reco2: [ pandora, pandoraTrack, pandoraShower, pandoraShowerSBN, pandoraCaloData, pandoraPidData, caloskimCalorimetry,
-  crtclustering, crtspacepoints, crttracks, cvn ]
+  crtclustering, crtspacepoints, crttracks, cvn, opt0finder]
 
 physics.analyzers.caloskim.G4producer: ""
 physics.analyzers.caloskim.SimChannelproducer: ""

sbndcode/OpT0Finder/SBNDOpT0Finder_module.cc

@@ -631,6 +631,8 @@ bool SBNDOpT0Finder::ConstructLightClusters(art::Event& e, unsigned int tpc) {
   art::FindManyP<recob::SpacePoint> shw_to_spacepoints(shw_h, e, _shw_producer);
   art::FindManyP<recob::Hit> spacepoint_to_hits (spacepoint_h, e, _slice_producer);
 
+  art::FindOneP<larpandoraobj::PFParticleMetadata> pfp_to_metadata(pfp_h, e, _slice_producer);
+
   // Loop over the Slices
   for (size_t n_slice = 0; n_slice < slice_h->size(); n_slice++) {
     flashmatch::QCluster_t light_cluster;
@@ -653,15 +655,17 @@ bool SBNDOpT0Finder::ConstructLightClusters(art::Event& e, unsigned int tpc) {
     std::vector<art::Ptr<recob::PFParticle>> pfp_v = slice_to_pfps.at(n_slice);
 
     if (_select_nus){
-      bool nu_pfp = false;
+      bool found_clear_cosmic = false;
       for (size_t n_pfp = 0; n_pfp < pfp_v.size(); n_pfp++) {
         auto pfp = pfp_v[n_pfp];
-        unsigned pfpPDGC = std::abs(pfp->PdgCode());
-        if ((pfpPDGC == 12) || (pfpPDGC == 14) || (pfpPDGC == 16))
-          nu_pfp = true;
+        const auto meta  = pfp_to_metadata.at(pfp.key());
+        const auto props = meta->GetPropertiesMap();
+        if (props.count("IsClearCosmic")){
+          found_clear_cosmic = true;
+          break;
+        }
       }
-      if (nu_pfp == false)
-        break;
+      if (found_clear_cosmic) continue;
     }
 
     for (size_t n_pfp = 0; n_pfp < pfp_v.size(); n_pfp++) {

sbndcode/OpT0Finder/job/opt0finder_sbnd.fcl

@@ -70,17 +70,6 @@ sbnd_opt0_finder_one_to_many.FlashMatchConfig.QLLMatch.ChiErrorWidth:      0.1
 sbnd_opt0_finder_one_to_many.FlashMatchConfig.QLLMatch.UseMinuit:          false   # do not perform minimiziation!!; Minuit not needed for one-to-many matching
 sbnd_opt0_finder_one_to_many.FlashMatchConfig.QLLMatch.SaturatedThreshold: 2.5e3   # fix for saturated PMTs, determines when to evaluate non-linearities 
 
-# ... for TPC 0
-sbnd_opt0_finder_one_to_many_tpc0: @local::sbnd_opt0_finder_one_to_many
-sbnd_opt0_finder_one_to_many_tpc0.TPC: 0
-sbnd_opt0_finder_one_to_many_tpc0.OpFlashProducer: ["opflashtpc0","opflashtpc0xarapuca"]
-
-# ... for TPC 1
-sbnd_opt0_finder_one_to_many_tpc1: @local::sbnd_opt0_finder_one_to_many
-sbnd_opt0_finder_one_to_many_tpc1.TPC: 1
-sbnd_opt0_finder_one_to_many_tpc1.OpFlashProducer: ["opflashtpc1","opflashtpc1xarapuca"]
-
-
 # WARNING: many-to-many matching is untested, performance unknown! 
 #
 # Configuration to run the flash matching in many flashes to many slices configuration

sbndcode/OpT0Finder/job/opt0finder_sbnd_data.fcl

@@ -0,0 +1,22 @@
+#include "opt0finder_sbnd.fcl"
+
+BEGIN_PROLOG
+
+sbnd_opt0_finder_data: @local::sbnd_opt0_finder_one_to_many
+
+sbnd_opt0_finder_data.CaloProducer:     "pandoraCaloData"
+sbnd_opt0_finder_data.SelectNeutrino:   true 
+sbnd_opt0_finder_data.CalAreaConstants: [ 0.0204 , 0.0210, 0.0193 ]
+sbnd_opt0_finder_data.OpChannelsToSkip: [39, 66, 67, 69, 71, 85, 86, 87, 92, 115, 138, 141, 170, 197, 217, 218, 221, 222, 223, 226, 245, 248, 249, 302]
+# since the flash time is w.r.t the TDC ETRIG, this may move with light arrival; need to extend the flash window
+sbnd_opt0_finder_data.FlashVetoTimeStart: -1
+sbnd_opt0_finder_data.FlashVetoTimeEnd:    5
+
+sbnd_opt0_finder_data.FlashMatchConfig.QLLMatch.ChiErrorWidth: 0.25 # increase the error width from 0.1
+sbnd_opt0_finder_data.FlashMatchConfig.PhotonLibHypothesis.GlobalQE:     0.25 # effectively a global scaling factor
+sbnd_opt0_finder_data.FlashMatchConfig.PhotonLibHypothesis.GlobalQERefl: 0.25 # effectively a global scaling factor
+
+# set to arbitrary high value to ignore non-linearity correction until we know it better
+sbnd_opt0_finder_data.FlashMatchConfig.QLLMatch.SaturatedThreshold: 1e9  
+
+END_PROLOG

sbndcode/OpT0Finder/job/run_opt0finder_sbnd_data.fcl

@@ -0,0 +1,83 @@
+#include "services_sbnd.fcl"
+#include "simulationservices_sbnd.fcl"
+#include "messages_sbnd.fcl"
+#include "sam_sbnd.fcl"
+
+#include "opt0finder_sbnd_data.fcl"
+#include "opt0finderana_sbnd.fcl"
+
+process_name: OpT0Finder
+
+services:
+{
+  # Load the service that manages root files for histograms.
+  TFileService:              { fileName: "opt0finder_tree.root" }
+  RandomNumberGenerator:     {}
+  message:                   @local::sbnd_message_services_prod # from messages_sbnd.fcl
+  FileCatalogMetadata:       @local::sbnd_file_catalog_mc       # from sam_sbnd.fcl
+                             @table::sbnd_services              # from services_sbnd.fcl
+                             @table::sbnd_g4_services
+}
+
+
+#source is now a root file
+source:
+{
+  module_type: RootInput
+  maxEvents:  -1  # Number of events to create
+}
+
+
+#block to define where the output goes.  if you defined a filter in the physics
+#block and put it in the trigger_paths then you need to put a SelectEvents: {SelectEvents: [XXX]}
+#entry in the output stream you want those to go to, where XXX is the label of the filter module(s)
+outputs:
+{
+  out1:
+  {
+    module_type: RootOutput
+    fileName:    "opt0finder-art.root"
+    dataTier:    "reconstructed"
+    compressionLevel: 1
+    SelectEvents: ["reco"]
+ }
+}
+
+# Define and configure some modules to do work on each event.
+# First modules are defined; they are scheduled later.
+# Modules are grouped by type.
+physics:
+{
+  producers:
+  {
+    opt0finder:    @local::sbnd_opt0_finder_data
+  }
+
+  analyzers:
+  {
+    opt0finderana:     @local::sbnd_opt0_finder_ana
+  }
+
+
+  #define the producer and filter modules for this path, order matters,
+  #filters reject all following items.  see lines starting physics.producers below
+  reco: [
+    opt0finder
+  ]
+
+  ana: [
+    opt0finderana
+  ]
+
+  # define the output stream, there could be more than one if using filters
+  stream1: [out1]
+
+  # trigger_paths is a keyword and contains the paths that modify the art::event,
+  # ie filters and producers
+  trigger_paths: [reco]
+
+  # end_paths is a keyword and contains the paths that do not modify the art::Event,
+  # ie analyzers and output streams.  these all run simultaneously
+  end_paths: [ana, stream1]
+}
+