update slice selection (ignore clear cosmics), and add scaling factor

Author: lynnt20

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_data.fcl

@@ -5,14 +5,18 @@ BEGIN_PROLOG
 sbnd_opt0_finder_data: @local::sbnd_opt0_finder_one_to_many
 
 sbnd_opt0_finder_data.CaloProducer:     "pandoraCaloData"
-sbnd_opt0_finder_data.SelectNeutrino:   false # for data, unsure about trusting pandora to select neutrinos
+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, # OFF/uncalibrated PMTs
-                                        134,135,150,151,152,153,154,155,156,157,158,159,160,161,176,177] # APSIA channels
+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
 
-# set to arbitrary high value to ignore saturation correction until we know it better
+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