update reco1 stage

Author: Brinden Carlson

sbndcode/JobConfigurations/README.md

@@ -6,14 +6,22 @@ At the time of writing, the core workflow (for BNB + Dirt + Cosmics) is the foll
 
 - `prodoverlay_corsika_cosmics_proton_genie_rockbox_sce.fcl`
 - `g4_sce_dirt_filter_lite.fcl` (does not include TPC electron drift simulation)
-- `detsim_sce_lite.fcl` (includes TPC drift simulation, TPC electronics simulation, and signal processing)
+- `detsim_sce_overlay_lite.fcl` (includes TPC drift simulation, TPC electronics simulation, and signal processing)
 - `reco1_sce_lite.fcl`
 - `reco2_sce.fcl`
 
+For single generator workflows (like intrinsic neutrino samples) the workflow is the follow:
+
+- `<your-gen>.fcl`
+- `g4_sce_lite.fcl`
+- `detsim_sce_lite.fcl`
+- `reco1_sce_no_overlay_lite.fcl`
+- `reco2_sce.fcl`
+
 This may well change over the coming months, and this README should be updated to reflect this.
 
 ** ALERT **
 
 Due to changes implemented in sbndcode PRs #408 and #409 the 1D simulation fcls will not work out of the box, they will need editing by experts!
 
-Henry Lay & Lynn Tung - Feb 2024
+Henry Lay, Lynn Tung, Bear Carlson - Feb 2024

sbndcode/JobConfigurations/base/reco_drops.fcl

@@ -19,7 +19,10 @@ BEGIN_PROLOG
 
 reco_drops: [ @sequence::detsim_drops,
               "drop raw::RawDigits_*_*_*",
-              "drop raw::OpDetWaveforms_*_*_*" ]
+              "drop raw::OpDetWaveforms_*_*_*",
+              "drop *_cluster3d_*_*", #drop all mlreco output
+              "drop *_simplemerge_*_*" #drop all mlreco output
+               ]
 
 END_PROLOG
 

sbndcode/JobConfigurations/base/reco_sbnd.fcl

@@ -48,6 +48,10 @@
 #include "scecorrections.fcl"
 #include "flashmatch_sbnd.fcl"
 
+#include "simplemerge.fcl" 
+#include "mlreco_sbnd.fcl"
+
+
 ###include "hitparticleassociations.fcl"
 
 #include "sbnd_trackcalo_skimmer.fcl"
@@ -189,6 +193,10 @@ physics:
   ### Uncalibrated calorimetry producer for calibration caloskimmer
   caloskimCalorimetry: @local::caloskim_calorimetry
 
+  #ML Reco
+  cluster3d: @local::standard_cluster3d
+  simplemerge: @local::simplemerge
+
  }
 
  #define the producer and filter modules for this path, order matters,
@@ -211,6 +219,8 @@ physics:
    #      , corner
    #      , fuzzycluster
          , crtstrips
+         , simplemerge #merges dropped particle collections into nominal
+         , cluster3d  
  ]
 
  reco2_no_opt0finder: [
@@ -245,9 +255,11 @@ physics:
 
  analyzers:{
   caloskim: @local::caloskim_nodigits_goldentracks
+  supera: @local::sbnd_supera
  }
 
  caloskimana: [ caloskim ]
+ superaana: [ supera ] #outputs larcv.root for use in mlreco
 
  #define the output stream, there could be more than one if using filters
  stream1:  [ out1 ]
@@ -276,7 +288,9 @@ outputs:
          #"drop sim::SimChannels_*_*_*",
          "drop sim::OpDetBacktrackerRecords_*_*_*",
          "drop *_linecluster_*_*",
-         "drop *_fasthit_*_*"
+         "drop *_fasthit_*_*",
+         "drop *_cluster3d_*_*", #drop all mlreco output
+         "drop *_simplemerge_*_*" #drop all mlreco output
      ]
 
 

sbndcode/JobConfigurations/standard/reco/reco1_sce_no_overlay_lite.fcl

@@ -0,0 +1,13 @@
+#-------------------------------------------------------------------
+#
+# Name: reco1_sce_no_overlay_lite.fcl
+#
+# Purpose: Handles simplemerge labels from single source generators
+#
+# Created: 28-Feb-2024 B. Carlson (bcarlson1@ufl.edu)
+#
+#-------------------------------------------------------------------
+
+#include "reco1_sce_lite.fcl"
+
+physics.producers.simplemerge.InputSourcesLabels: ["largeant","largeant:droppedMCParticles"]

sbndcode/JobConfigurations/standard/standard_reco1_sbnd.fcl

@@ -6,11 +6,16 @@
 # 20200923 (d.brailsford@lancster.ac.uk) [v1.0]
 #   Initial definition
 #
+# Changes:
+# 20240228 (bcarlson1@ufl.edu) [v1.1]
+#    Add supera output using an analyzer end path
 # job configuration
 
 #include "reco_sbnd.fcl"
 
 process_name: Reco1
 
 physics.trigger_paths: [ reco1 ]
+physics.end_paths: [stream1, superaana]
+
 

sbndcode/JobConfigurations/standard/standard_reco1reco2_sbnd.fcl

@@ -13,3 +13,4 @@
 process_name: Reco1Reco2
 
 physics.trigger_paths: [ fullreco ]
+physics.end_paths: [stream1, superaana]

sbndcode/LArSoftConfigurations/mlreco_sbnd.fcl

@@ -7,4 +7,24 @@ sbnd_sedlite: @local::sbn_largeant_info_reducer
 sbnd_sedlite.SimEnergyDepositLabel: "largeant:LArG4DetectorServicevolTPCActive"
 sbnd_sedlite.useOrigTrackID: true #needed since TrackID overwritten by SBND modules
 
+#supera for mpvmpr
+sbnd_supera_mpvmpr: {
+    module_type:     "LArSoftSuperaDriver"
+    supera_params:   "supera_sbnd_mpvmpr.fcl"
+    out_filename:    "larcv.root"
+    unique_filename: false
+    stream:          "mc"
+    StrictDataLoading: false
+}
+
+#supera for everything else
+sbnd_supera: {
+    module_type:     "LArSoftSuperaDriver"
+    supera_params:   "supera_sbnd_corsika.fcl"
+    out_filename:    "larcv.root"
+    unique_filename: false
+    stream:          "mc"
+    StrictDataLoading: false
+}
+
 END_PROLOG