Introduce Barycenter Flash Match

fcl/caf/cafmaker_defs.fcl

@@ -326,6 +326,7 @@ cafmaker.TriggerLabel: "emuTrigger"
 cafmaker.FlashTrigLabel: "" # unavailable
 cafmaker.SimChannelLabel: "largeant"
 cafmaker.SystWeightLabels: ["genieweight", "fluxweight"]
+cafmaker.TPCPMTBarycenterMatchLabel: "tpcpmtbarycentermatch"
 cafmaker.TrackHitFillRREndCut: 30 # include entire PID region
 
 # Add CAFMaker to the list of producers

fcl/reco/Definitions/stage1_icarus_defs.fcl

@@ -21,6 +21,7 @@
 ## The below can be found from the softlink to Supera in sbncode
 #include "supera_modules.fcl"
 #include "crtpmtmatching_parameters.fcl"
+#include "tpcpmtbarycentermatch_config.fcl"
 
 BEGIN_PROLOG
 
@@ -65,6 +66,9 @@ icarus_stage1_producers:
   CRTT0Matching:     @local::standard_crtt0producer
   CRTT0MatchingW:    @local::standard_crtt0producerW
   CRTT0MatchingE:    @local::standard_crtt0producerE
+
+  tpcpmtbarycentermatchCryoE: @local::data_tpcpmtbarycentermatchproducer_east
+  tpcpmtbarycentermatchCryoW: @local::data_tpcpmtbarycentermatchproducer_west
 }
 
 icarus_stage1_filters:
@@ -191,6 +195,11 @@ icarus_pandora_Gauss:              [
 icarus_reco_fm:                    [ fmatchCryoE,
                                      fmatchCryoW ]
 
+icarus_tpcpmtbarycentermatch:      [
+                                     tpcpmtbarycentermatchCryoE, 
+                                     tpcpmtbarycentermatchCryoW
+                                   ]
+
 icarus_crttrack:                   [crttrack]
 
 icarus_crtt0match:                 [CRTT0Matching]

fcl/reco/Stage1/Run2/stage1_run2_icarus.fcl

@@ -5,6 +5,7 @@ process_name: stage1
 physics.reco: [ @sequence::icarus_filter_cluster3D,
                 @sequence::icarus_pandora_Gauss,
                 @sequence::icarus_reco_fm,
+                @sequence::icarus_tpcpmtbarycentermatch,
                 @sequence::icarus_crttrack,
                 @sequence::icarus_crtt0match,
                 caloskimCalorimetryCryoE, caloskimCalorimetryCryoW]

fcl/reco/Stage1/Run2/stage1_run2_icarus_MC.fcl

@@ -28,10 +28,15 @@ physics.producers: {
             mcassociationsGausCryoW:        @local::standard_mcparticlehitmatching
 }
 
+#Redefine TPCPMTBarycenterMatch producers with MC parameters
+physics.producers.tpcpmtbarycentermatchCryoE: @local::mc_tpcpmtbarycentermatchproducer_east
+physics.producers.tpcpmtbarycentermatchCryoW: @local::mc_tpcpmtbarycentermatchproducer_west
+
 physics.reco: [ 
                 @sequence::icarus_reco_Gauss_CryoE ,
                 @sequence::icarus_reco_Gauss_CryoW ,
 		            @sequence::icarus_reco_fm,
+                @sequence::icarus_tpcpmtbarycentermatch,
                 caloskimCalorimetryCryoE, caloskimCalorimetryCryoW,
                 mcassociationsGausCryoE,  mcassociationsGausCryoW
               ]

icaruscode/PMT/OpReco/CMakeLists.txt

@@ -16,6 +16,7 @@ cet_build_plugin(ICARUSOpHitFinder art::module
 
 
 set(	MODULE_LIBRARIES
+                        icarusalg::Utilities
 			sbnobj::Common_Trigger
 			larcorealg::Geometry
 			larcore::Geometry_Geometry_service
@@ -59,6 +60,7 @@ cet_build_plugin(ICARUSOpFlashAna art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(ICARUSOpHitAna art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(ICARUSOpHitTuple art::module LIBRARIES ${MODULE_LIBRARIES})
 cet_build_plugin(ICARUSParticleAna art::module LIBRARIES ${MODULE_LIBRARIES})
+cet_build_plugin(TPCPMTBarycenterMatchProducer art::module LIBRARIES ${MODULE_LIBRARIES})
 
 
 install_headers()

icaruscode/PMT/OpReco/fcl/tpcpmtbarycentermatch_config.fcl

@@ -0,0 +1,98 @@
+BEGIN_PROLOG
+
+#Common paramters for the barycenter matching
+tpcpmtbarycentermatch_common_params:
+{
+    OpFlashLabel:     "opflash"
+    PandoraLabel:     "pandoraGaus"
+    CollectionOnly:   true
+    UseTimeRange:     true
+    Verbose:          false
+    FillMatchTree:    false
+    TriggerTolerance: 0.15
+    TimeRangeMargin:  35.
+}
+
+
+#Parameters specific to data and Monte Carlo separately
+######
+
+tpcpmtbarycentermatch_data_params:
+{   
+    TriggerLabel:  "daqTrigger"
+    TriggerDelay:  0.6
+}
+
+tpcpmtbarycentermatch_mc_params:
+{
+    TriggerLabel:  "emuTrigger"
+    TriggerDelay:  0.02
+}
+
+######
+
+
+#Separate configurations for data or MC for running the matching over BOTH cryostats
+######
+
+data_tpcpmtbarycentermatchproducer:
+{ 
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_data_params
+    InputTags:         ["CryoE","CryoW"]
+}
+
+mc_tpcpmtbarycentermatchproducer:
+{
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_mc_params
+    InputTags:         ["CryoE","CryoW"]
+}
+
+######
+
+#Separate configurations for data or MC for running the matching over the east cryostat
+######
+
+data_tpcpmtbarycentermatchproducer_east:
+{
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_data_params
+    InputTags:         ["CryoE"]
+}
+
+mc_tpcpmtbarycentermatchproducer_east:
+{ 
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_mc_params
+    InputTags:         ["CryoE"]
+}
+
+######
+
+#Separate configurations for data or MC for running the matching over the west cryostat
+######
+
+data_tpcpmtbarycentermatchproducer_west:
+{
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_data_params
+    InputTags:         ["CryoW"]
+}
+
+mc_tpcpmtbarycentermatchproducer_west:
+{
+  module_type: "icaruscode/PMT/OpReco/TPCPMTBarycenterMatchProducer"
+    @table::tpcpmtbarycentermatch_common_params
+    @table::tpcpmtbarycentermatch_mc_params
+    InputTags:         ["CryoW"]
+}
+
+######
+
+END_PROLOG

icaruscode/PMT/OpReco/fcl/tpcpmtbarycentermatch_data.fcl

@@ -0,0 +1,59 @@
+# A script to run the producer module: TPCPMTBarycenterMatchProducer
+
+#include "tpcpmtbarycentermatch_config.fcl"
+
+#include "services_icarus_simulation.fcl"
+#include "services_common_icarus.fcl"
+
+process_name: TPCPMTBarycenterMatchProducer
+
+services:
+{
+  # Load the service that manages root files for histograms.
+  TFileService:           { fileName: "tpcpmtmatchingana.root" }
+ # This constrols the display in the output of how long each job step takes for each event. 
+  TimeTracker:            {}
+
+
+
+  # This parameter controls the level of descriptive output from various LArSoft modules.
+  message:                @local::icarus_message_services_prod_debug
+                          @table::icarus_common_services
+
+} # services
+
+
+outputs:
+{
+  out1:
+  {
+    module_type: RootOutput
+    fileName:    "%ifb_%tc_TPCPMTBarycenterMatch.root"
+    dataTier: "reconstructed"
+    saveMemoryObjectThreshold: 0
+    compressionLevel: 1
+    fastCloning: true
+  }
+}
+
+# The 'physics' section defines and configures some modules to do work on each event.
+physics:
+{
+   producers:
+  {
+    tpcpmtbarycentermatchCryoE: @local::data_tpcpmtbarycentermatchproducer_east
+    tpcpmtbarycentermatchCryoW: @local::data_tpcpmtbarycentermatchproducer_west
+  }
+
+  # Schedule job step(s) for execution by defining the analysis module for this job.
+  reco: [ tpcpmtbarycentermatchCryoE, tpcpmtbarycentermatchCryoW ]
+
+  stream1: [ out1 ]
+
+  trigger_paths: [reco]
+
+  # "end_paths" is a keyword and contains the modules that do not modify the art::Event;
+  # i.e., analyzers and output streams.
+  end_paths: [ stream1 ]
+
+}

icaruscode/PMT/OpReco/fcl/tpcpmtbarycentermatch_mc.fcl

@@ -0,0 +1,58 @@
+# A script to run the producer module: TPCPMTBarycenterMatchProducer
+
+#include "tpcpmtbarycentermatch_config.fcl"
+
+#include "services_icarus_simulation.fcl"
+#include "services_common_icarus.fcl"
+
+process_name: TPCPMTBarycenterMatchProducer
+
+services:
+{
+  # Load the service that manages root files for histograms.
+  TFileService:           { fileName: "tpcpmtmatchingana.root" }
+ # This constrols the display in the output of how long each job step takes for each event. 
+  TimeTracker:            {}
+
+
+
+  # This parameter controls the level of descriptive output from various LArSoft modules.
+  message:                @local::icarus_message_services_prod_debug
+                          @table::icarus_common_services
+
+} # services
+
+outputs:
+{
+  out1:
+  {
+    module_type: RootOutput
+    fileName:    "%ifb_%tc_TPCPMTBarycenterMatch.root"
+    dataTier: "reconstructed"
+    saveMemoryObjectThreshold: 0
+    compressionLevel: 1
+    fastCloning: true
+  }
+}
+
+# The 'physics' section defines and configures some modules to do work on each event.
+physics:
+{
+   producers:
+  {
+    tpcpmtbarycentermatchCryoE: @local::mc_tpcpmtbarycentermatchproducer_east
+    tpcpmtbarycentermatchCryoW: @local::mc_tpcpmtbarycentermatchproducer_west
+  }
+
+  # Schedule job step(s) for execution by defining the analysis module for this job.
+  reco: [ tpcpmtbarycentermatchCryoE, tpcpmtbarycentermatchCryoW ]
+
+  stream1: [ out1 ]
+
+  trigger_paths: [reco]
+
+  # "end_paths" is a keyword and contains the modules that do not modify the art::Event;
+  # i.e., analyzers and output streams.
+  end_paths: [ stream1 ]
+
+}