diff --git a/fcl/CMakeLists.txt b/fcl/CMakeLists.txt
index 0d1ef85ed..f48ba48e4 100644
--- a/fcl/CMakeLists.txt
+++ b/fcl/CMakeLists.txt
@@ -9,6 +9,7 @@ add_subdirectory(caf)
 add_subdirectory(services)
 add_subdirectory(utilities)
 add_subdirectory(compatibility)
+add_subdirectory(numi-anaA)
 add_subdirectory(syst_variations)
 
 install_fhicl()
diff --git a/fcl/caf/SystTools/configs/config_test_Multisigma.fcl b/fcl/caf/SystTools/configs/config_test_Multisigma.fcl
new file mode 100644
index 000000000..b952eaa40
--- /dev/null
+++ b/fcl/caf/SystTools/configs/config_test_Multisigma.fcl
@@ -0,0 +1,44 @@
+GENIEReWeight_Multisigma_Tool_Config: {
+
+    tool_type: "GENIEReWeight"
+
+    # Updated 25 July 2023 by S. Gardiner
+    instance_name: "ICARUS_v2_multisigma"
+
+    # Keep one GENIEReWeight instance per discrete dial tweak in memory
+    # Requires no reconfigures within the event loop
+    UseFullHERG: true
+
+    ignore_parameter_dependence: true
+
+    ################################## CCQE Parameters
+
+    # This dial varies the strength of the RPA corrections that are applied in
+    # the Valencia CCQE model. The implementation was contributed to GENIE by
+    # MicroBooNE and is documented in Phys. Rev. D 105, 072001 (2022).
+    RPA_CCQE_central_value: 0
+    RPA_CCQE_variation_descriptor: "[-1,+1,-2,+2,-3,+3]"
+
+    # This dial varies the strength of the electromagnetic potential used by
+    # the Valencia CCQE model to apply Coulomb corrections to the cross
+    # section. The implementation was contributed to GENIE by MicroBooNE and is
+    # documented in Phys. Rev. D 105, 072001 (2022).
+    CoulombCCQE_central_value: 0
+    CoulombCCQE_variation_descriptor: "[-1,+1,-2,+2,-3,+3]"
+
+    # Photon angular distribution in Delta -> photon N decays
+    # ThetaDelta2NRad = 0 gives the default isotropic angular distribution
+    # ThetaDelta2NRad = 1 gives an alternative distribution proportional to
+    # cos^2(theta). See also Phys. Rev. D 105, 072001 (2022).
+    ThetaDelta2NRad_central_value: 0
+    ThetaDelta2NRad_variation_descriptor: "[1.0]"
+
+    NormCCCOH_central_value: 0
+    NormCCCOH_variation_descriptor: "[-1,+1,-2,+2,-3,+3]"
+
+    NormNCCOH_central_value: 0
+    NormNCCOH_variation_descriptor: "[-1,+1,-2,+2,-3,+3]"
+
+}
+
+syst_providers: [GENIEReWeight_Multisigma_Tool_Config]
diff --git a/fcl/caf/cafmaker_defs.fcl b/fcl/caf/cafmaker_defs.fcl
index 5c8e0e554..e5817d351 100644
--- a/fcl/caf/cafmaker_defs.fcl
+++ b/fcl/caf/cafmaker_defs.fcl
@@ -14,6 +14,7 @@
 #include "eventweight_genie_sbn.fcl"
 #include "eventweight_genie_systtools.fcl"
 #include "eventweight_flux_sbn.fcl"
+#include "eventweight_geant4_sbn.fcl"
 #include "mcreco.fcl"
 #include "mcsproducer.fcl"
 #include "rangeproducer.fcl"
@@ -119,6 +120,7 @@ recoana_caf_preprocess_producers: {
   genieweight: @local::sbn_eventweight_genie
   systtools: @local::sbn_systtools
   fluxweight: @local::sbn_eventweight_flux
+  geant4weight: @local::sbn_eventweight_geant4
 }
 
 # Overwrite labels
@@ -138,6 +140,7 @@ recoana_caf_preprocess_producers.pandoraPidGausCryoW.CalorimetryModuleLabel: "pa
 
 recoana_caf_preprocess_producers.genieweight.weight_functions: @local::recoana_caf_preprocess_producers.genieweight.weight_functions_genie
 recoana_caf_preprocess_producers.fluxweight.weight_functions: @local::recoana_caf_preprocess_producers.fluxweight.weight_functions_flux
+recoana_caf_preprocess_producers.geant4weight.weight_functions: @local::recoana_caf_preprocess_producers.geant4weight.weight_functions_reint
 
 # Producers for making SCE products
 recoana_sce_producers: {
@@ -272,7 +275,7 @@ caf_preprocess_data_sequence: [
     ShowerCosmicDistCryoE, ShowerCosmicDistCryoW
     ]
 
-caf_preprocess_evtw_sequence: [@sequence::caf_preprocess_sequence, rns, genieweight, fluxweight]
+caf_preprocess_evtw_sequence: [@sequence::caf_preprocess_sequence, rns, genieweight, fluxweight, geant4weight]
 
 caf_preprocess_sce_sequence: [ mcreco,
   # Run the SCE correction
@@ -294,7 +297,7 @@ caf_preprocess_sce_sequence: [ mcreco,
   # TODO: rns??
 ]
 
-caf_preprocess_sce_evtw_sequence: [@sequence::caf_preprocess_sce_sequence, rns, genieweight, fluxweight]
+caf_preprocess_sce_evtw_sequence: [@sequence::caf_preprocess_sce_sequence, rns, genieweight, fluxweight, geant4weight]
 
 # CAFMaker config
 cafmaker: @local::standard_cafmaker
@@ -329,7 +332,7 @@ cafmaker.OpFlashLabel: "opflash"
 cafmaker.TriggerLabel: "emuTrigger"
 cafmaker.FlashTrigLabel: "" # unavailable
 cafmaker.SimChannelLabel: "largeant"
-cafmaker.SystWeightLabels: ["genieweight", "fluxweight"]
+cafmaker.SystWeightLabels: ["genieweight", "fluxweight", "geant4weight"]
 cafmaker.TPCPMTBarycenterMatchLabel: "tpcpmtbarycentermatch"
 cafmaker.TrackHitFillRREndCut: 30 # include entire PID region
 
diff --git a/fcl/caf/cafmakerjob_icarus_data_rereco_numi.fcl b/fcl/caf/cafmakerjob_icarus_data_rereco_numi.fcl
new file mode 100644
index 000000000..bf467fd7d
--- /dev/null
+++ b/fcl/caf/cafmakerjob_icarus_data_rereco_numi.fcl
@@ -0,0 +1,68 @@
+#include "stage1_icarus_defs.fcl"
+#include "cafmakerjob_icarus_data.fcl"
+
+process_name: stage1REDO
+
+services:
+{
+  TFileService:           { }
+  IICARUSChannelMap:      @local::icarus_channelmappinggservice
+                          @table::icarus_wirecalibration_services
+  SpaceChargeService: @local::icarus_spacecharge
+  IFBeam: {}
+}
+
+physics.producers: {
+  @table::physics.producers
+  @table::icarus_stage1_producers
+}
+
+physics.filters: {
+  @table::icarus_stage1_filters
+}
+
+physics.analyzers: {
+  @table::icarus_stage1_analyzers
+}
+
+physics.runprod: [
+  @sequence::icarus_filter_cluster3D,
+  @sequence::icarus_pandora_Gauss,
+  # @sequence::icarus_reco_fm,
+  # @sequence::icarus_crttrack,
+  # @sequence::icarus_crtt0match,
+  caloskimCalorimetryCryoE, caloskimCalorimetryCryoW,
+  
+  @sequence::physics.runprod
+]
+physics.outana:            [ caloskimE, caloskimW]
+physics.stream1: []
+
+physics.trigger_paths:     [ runprod ]
+physics.end_paths:         [ outana, stream1 ]
+
+source.inputCommands: [
+  "keep *_*_*_*",
+  "drop *_*_*_stage1"
+]
+
+physics.analyzers.caloskimE.SelectEvents: [ runprod ]
+physics.analyzers.caloskimW.SelectEvents: [ runprod ]
+
+services.SpaceChargeService: {
+    EnableCalEfieldSCE: false
+    EnableCalSpatialSCE: false
+    EnableCorrSCE: false
+    EnableSimEfieldSCE: false
+    EnableSimSpatialSCE: false
+    InputFilename: "SCEoffsets/SCEoffsets_ICARUS_E500_voxelTH3.root"
+    RepresentationType: "Voxelized_TH3"
+    service_provider: "SpaceChargeServiceICARUS"
+}
+
+### adding t0 info in the calorimetry variables
+physics.producers.pandoraCaloGausCryoE.T0ModuleLabel: "pandoraGausCryoE"
+physics.producers.pandoraCaloGausCryoW.T0ModuleLabel: "pandoraGausCryoW" 
+
+# physics.producers.pandoraCaloGausCryoE.ChargeMethod: 2 # SumADC
+# physics.producers.pandoraCaloGausCryoW.ChargeMethod: 2 # SumADC
diff --git a/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run1.fcl b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run1.fcl
new file mode 100644
index 000000000..cfcc0095b
--- /dev/null
+++ b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run1.fcl
@@ -0,0 +1,3 @@
+#include "cafmakerjob_icarus_data_rereco_numi.fcl"
+
+services.Geometry: @local::icarus_geometry_services_no_overburden_legacy_icarus_v3.Geometry
diff --git a/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OffBeam.fcl b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OffBeam.fcl
new file mode 100644
index 000000000..540458d17
--- /dev/null
+++ b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OffBeam.fcl
@@ -0,0 +1,4 @@
+#include "cafmakerjob_icarus_data_rereco_numi.fcl"
+
+# Get EXT
+physics.runprod: [numiextinfo, @sequence::physics.runprod]
diff --git a/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OnBeam.fcl b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OnBeam.fcl
new file mode 100644
index 000000000..ba8162933
--- /dev/null
+++ b/fcl/caf/cafmakerjob_icarus_data_rereco_numi_Run2OnBeam.fcl
@@ -0,0 +1,5 @@
+#include "cafmakerjob_icarus_data_rereco_numi.fcl"
+
+# Get POT
+# physics.runprod: [numiinfo, @sequence::physics.runprod]
+
diff --git a/fcl/caf/cafmakerjob_icarus_systtools_and_fluxwgt.fcl b/fcl/caf/cafmakerjob_icarus_systtools_and_fluxwgt.fcl
index 7ba1d48a3..75f7d6250 100644
--- a/fcl/caf/cafmakerjob_icarus_systtools_and_fluxwgt.fcl
+++ b/fcl/caf/cafmakerjob_icarus_systtools_and_fluxwgt.fcl
@@ -1,6 +1,6 @@
 #include "cafmakerjob_icarus.fcl"
 
-physics.runprod: [ @sequence::caf_preprocess_sequence, rns, systtools, fluxweight, cafmaker ]
+physics.runprod: [ @sequence::caf_preprocess_sequence, rns, systtools, fluxweight, geant4weight, cafmaker ]
 
-physics.producers.cafmaker.SystWeightLabels: ["systtools", "fluxweight"]
+physics.producers.cafmaker.SystWeightLabels: ["systtools", "fluxweight", "geant4weight"]
 
diff --git a/fcl/caf/cafmakerjob_rereco_numi_icarus_systtools_and_fluxwgt.fcl b/fcl/caf/cafmakerjob_rereco_numi_icarus_systtools_and_fluxwgt.fcl
new file mode 100644
index 000000000..e16f69c6b
--- /dev/null
+++ b/fcl/caf/cafmakerjob_rereco_numi_icarus_systtools_and_fluxwgt.fcl
@@ -0,0 +1,60 @@
+#include "stage1_icarus_defs.fcl"
+#include "services_icarus_simulation.fcl"
+#include "cafmakerjob_icarus.fcl"
+
+process_name: stage1REDO
+
+services:
+{
+  TFileService:           { }
+  IICARUSChannelMap:      @local::icarus_channelmappinggservice
+                          @table::icarus_wirecalibration_services
+  SpaceChargeService: @local::icarus_spacecharge
+  IFBeam: {}
+}
+services.ParticleInventoryService:  @local::standard_particleinventoryservice
+
+
+physics.producers: {
+  @table::physics.producers
+  @table::icarus_stage1_producers
+}
+
+physics.filters: {
+  @table::icarus_stage1_filters
+}
+
+physics.analyzers: {
+  @table::icarus_stage1_analyzers
+}
+
+physics.runprod: [
+  rns, 
+  @sequence::icarus_filter_cluster3D,
+  @sequence::icarus_pandora_Gauss,
+  # @sequence::icarus_reco_fm,
+  # @sequence::icarus_crttrack,
+  # @sequence::icarus_crtt0match,
+  caloskimCalorimetryCryoE, caloskimCalorimetryCryoW,
+  systtools, fluxweight, geant4weight,
+  
+  @sequence::physics.runprod
+]
+
+physics.outana:            [ caloskimE, caloskimW]
+physics.stream1: []
+
+physics.trigger_paths:     [ runprod ]
+physics.end_paths:         [ outana, stream1 ]
+
+source.inputCommands: [
+  "keep *_*_*_*",
+  "drop *_*_*_stage1"
+]
+
+physics.analyzers.caloskimE.SelectEvents: [ runprod ]
+physics.analyzers.caloskimW.SelectEvents: [ runprod ]
+
+physics.producers.cafmaker.SystWeightLabels: ["systtools", "fluxweight", "geant4weight"]
+#include "cafmaker_add_detsim2d_icarus.fcl"
+
diff --git a/fcl/gen/MultiVertex/cohlikevertex_icarus.fcl b/fcl/gen/MultiVertex/cohlikevertex_icarus.fcl
new file mode 100644
index 000000000..0e6501bc9
--- /dev/null
+++ b/fcl/gen/MultiVertex/cohlikevertex_icarus.fcl
@@ -0,0 +1,91 @@
+#include "services_icarus_simulation.fcl"
+#include "beamgates_icarus.fcl"
+
+process_name: CohLikeVertexGen
+
+CohLikeVertex: {
+    module_type  : "MultiPartVertex"
+    DebugMode    : 1
+    G4Time       : 1500
+    G4TimeJitter : 0
+    TPCRange     : [[0,0], [0,1], [0,2], [0,3], [1,0], [1,1], [1,2], [1,3]]
+    XRange       : [10]
+    YRange       : [10]
+    ZRange       : [10, 100]
+    MultiMax     : 2
+    MultiMin     : 2
+    ParticleParameter: {
+        PDGCode      : [   [13, -13], [211,-211]]
+        MinMulti     : [          1,           1]
+        MaxMulti     : [          1,           1]
+        ProbWeight   : [          1,           1]
+        KERange      : [ [0.2,4],   [0.2,4]]
+        MomRange     : []
+    }
+}
+
+services:
+{
+  scheduler:    { defaultExceptions: false }    # Make all uncaught exceptions fatal.
+  # Load the service that manages root files for histograms.
+  TFileService: { fileName: "multivertex_hist_uboone.root" }
+  TimeTracker:  {}
+  MemoryTracker:     { ignoreTotal: 1 } # default is one
+  RandomNumberGenerator: {} #ART native random number generator
+  #FileCatalogMetadata:  @local::art_file_catalog_mc
+  @table::icarus_gen_services
+}
+
+#Start each new event with an empty event.
+source:
+{
+  module_type: EmptyEvent
+  timestampPlugin: { plugin_type: "GeneratedEventTimestamp" }
+  maxEvents:   10          # Number of events to create
+  firstRun:    1           # Run number to use for this file
+  firstEvent:  1           # number of first event in the file
+}
+
+# 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:
+ {
+   generator: @local::CohLikeVertex
+   beamgate:  @local::icarus_standardbeamgate_NuMI  # from beamgates_icarus.fcl
+   rns:       { module_type: "RandomNumberSaver" }
+ }
+
+ #define the producer and filter modules for this path, order matters, 
+ #filters reject all following items.  see lines starting physics.producers below
+ simulate: [ rns, generator, beamgate ]
+
+ #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: [simulate] 
+
+ #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:     [stream1]  
+}
+
+#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:    "prod_multivertex_icarus_%tc_gen.root"
+   dataTier:    "generated"
+   compressionLevel: 1
+   saveMemoryObjectThreshold: 0
+ }
+}
diff --git a/fcl/gen/mevprtl/alp/alp_gen_icarus_M350.fcl b/fcl/gen/mevprtl/alp/alp_gen_icarus_M350.fcl
new file mode 100644
index 000000000..12a8568ee
--- /dev/null
+++ b/fcl/gen/mevprtl/alp/alp_gen_icarus_M350.fcl
@@ -0,0 +1,4 @@
+#include "alp_corsika_p_gen_icarus.fcl"
+
+alpM: 0.35
+#include "set_alpM.fcl"
diff --git a/fcl/gen/mevprtl/alp/alp_gen_icarus_M400.fcl b/fcl/gen/mevprtl/alp/alp_gen_icarus_M400.fcl
new file mode 100644
index 000000000..a3bf280d0
--- /dev/null
+++ b/fcl/gen/mevprtl/alp/alp_gen_icarus_M400.fcl
@@ -0,0 +1,4 @@
+#include "alp_corsika_p_gen_icarus.fcl"
+
+alpM: 0.4
+#include "set_alpM.fcl"
diff --git a/fcl/gen/mevprtl/alp/alp_gen_icarus_M450.fcl b/fcl/gen/mevprtl/alp/alp_gen_icarus_M450.fcl
new file mode 100644
index 000000000..c142a4e62
--- /dev/null
+++ b/fcl/gen/mevprtl/alp/alp_gen_icarus_M450.fcl
@@ -0,0 +1,4 @@
+#include "alp_corsika_p_gen_icarus.fcl"
+
+alpM: 0.45
+#include "set_alpM.fcl"
diff --git a/fcl/gen/mevprtl/alp/alp_gen_icarus_M650.fcl b/fcl/gen/mevprtl/alp/alp_gen_icarus_M650.fcl
new file mode 100644
index 000000000..8497209ff
--- /dev/null
+++ b/fcl/gen/mevprtl/alp/alp_gen_icarus_M650.fcl
@@ -0,0 +1,6 @@
+#include "alp_corsika_p_gen_icarus.fcl"
+
+alpM: 0.65
+#include "set_alpM.fcl"
+alpfa: 1e6
+#include "set_alpfa.fcl"
diff --git a/fcl/gen/mevprtl/alp/alp_gen_icarus_M750.fcl b/fcl/gen/mevprtl/alp/alp_gen_icarus_M750.fcl
new file mode 100644
index 000000000..d4cc78e35
--- /dev/null
+++ b/fcl/gen/mevprtl/alp/alp_gen_icarus_M750.fcl
@@ -0,0 +1,6 @@
+#include "alp_corsika_p_gen_icarus.fcl"
+
+alpM: 0.75
+#include "set_alpM.fcl"
+alpfa: 1e6
+#include "set_alpfa.fcl"
diff --git a/fcl/gen/mevprtl/higgs/dissH_M280.fcl b/fcl/gen/mevprtl/higgs/dissH_M280.fcl
new file mode 100644
index 000000000..e6fbc4b64
--- /dev/null
+++ b/fcl/gen/mevprtl/higgs/dissH_M280.fcl
@@ -0,0 +1,13 @@
+#include "dissonant_higgs_gen.fcl"
+
+physics.producers.generator.Decay.AllowPionDecay: false
+physics.producers.generator.Decay.AllowPi0Decay: false
+
+physics.producers.generator.Decay.ReferenceHiggsMixing: 3e-5
+physics.producers.generator.Flux.MixingAngle: 3e-5
+
+physics.producers.generator.RayTrace.MaxWeightFudge: 1.
+physics.producers.generator.RayTrace.ReferenceLabSolidAngle: 1e-4
+
+higgsM: 0.28
+#include "set_higgsM.fcl"
diff --git a/fcl/gen/mevprtl/higgs/dissH_M300.fcl b/fcl/gen/mevprtl/higgs/dissH_M300.fcl
index 46854f064..dc75aff32 100644
--- a/fcl/gen/mevprtl/higgs/dissH_M300.fcl
+++ b/fcl/gen/mevprtl/higgs/dissH_M300.fcl
@@ -1,4 +1,13 @@
 #include "dissonant_higgs_gen.fcl"
 
+physics.producers.generator.Decay.AllowPionDecay: false
+physics.producers.generator.Decay.AllowPi0Decay: false
+
+physics.producers.generator.Decay.ReferenceHiggsMixing: 3e-5
+physics.producers.generator.Flux.MixingAngle: 3e-5
+
+physics.producers.generator.RayTrace.MaxWeightFudge: 1.
+physics.producers.generator.RayTrace.ReferenceLabSolidAngle: 1e-4
+
 higgsM: 0.3
 #include "set_higgsM.fcl"
diff --git a/fcl/gen/mevprtl/higgs/dissH_M330.fcl b/fcl/gen/mevprtl/higgs/dissH_M330.fcl
new file mode 100644
index 000000000..328b34f7a
--- /dev/null
+++ b/fcl/gen/mevprtl/higgs/dissH_M330.fcl
@@ -0,0 +1,13 @@
+#include "dissonant_higgs_gen.fcl"
+
+physics.producers.generator.Decay.AllowPionDecay: false
+physics.producers.generator.Decay.AllowPi0Decay: false
+
+physics.producers.generator.Decay.ReferenceHiggsMixing: 2e-5
+physics.producers.generator.Flux.MixingAngle: 2e-5
+
+physics.producers.generator.RayTrace.MaxWeightFudge: 1.
+physics.producers.generator.RayTrace.ReferenceLabSolidAngle: 1e-4
+
+higgsM: 0.33
+#include "set_higgsM.fcl"
diff --git a/fcl/gen/mevprtl/higgs/dissonant_higgs_gen.fcl b/fcl/gen/mevprtl/higgs/dissonant_higgs_gen.fcl
index 424c56c06..018783075 100644
--- a/fcl/gen/mevprtl/higgs/dissonant_higgs_gen.fcl
+++ b/fcl/gen/mevprtl/higgs/dissonant_higgs_gen.fcl
@@ -53,6 +53,12 @@ outputs: {
 
 }
 
+# NuMI-like-timing
+physics.producers.generator.Flux.SpillTimeConfig: @local::FNAL_NuMI_default # from beamspilltimings.fcl
+
+physics.producers.generator.Verbose: false
+phyiscs.producers.generator.RayTrace.Verbose: false
+
 outputs.out.fileName: "simulation_higgs_icarus_numi_%tc-%p.root"
 
 services.NuRandomService.policy: "random"
diff --git a/fcl/numi-anaA/CMakeLists.txt b/fcl/numi-anaA/CMakeLists.txt
new file mode 100644
index 000000000..13355789a
--- /dev/null
+++ b/fcl/numi-anaA/CMakeLists.txt
@@ -0,0 +1 @@
+install_fhicl()
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom.fcl
new file mode 100644
index 000000000..23bdd23e0
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom.fcl
@@ -0,0 +1,3 @@
+#include "detsim_2d_icarus.fcl"
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR_P0nom
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom_gain0-8x.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom_gain0-8x.fcl
new file mode 100644
index 000000000..7faaf7b48
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_P0nom_gain0-8x.fcl
@@ -0,0 +1,5 @@
+#include "detsim_2d_icarus.fcl"
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR_P0nom
+
+physics.producers.daq.wcls_main.structs.gain0: 13.641696
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainhi.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainhi.fcl
new file mode 100644
index 000000000..34302c68e
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainhi.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 11.8031771
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainlo.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainlo.fcl
new file mode 100644
index 000000000..ff199a1e0
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainlo.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 8.724087455
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvhi.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvhi.fcl
new file mode 100644
index 000000000..279636478
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvhi.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 12.316358759999998
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvlo.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvlo.fcl
new file mode 100644
index 000000000..52d60efa2
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0gainvlo.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 8.21090584
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0var.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0var.fcl
new file mode 100644
index 000000000..0ed634359
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind0var.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, turned up noise on each plane
+#
+#include "detsim_2d_icarus.fcl"
+
+# Add in the tuned field responses
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+# variation on ind 0
+physics.producers.daq.wcls_main.params.files_fields: "icarus_fnal_fit_ks_ind0var.json.bz2"
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin0.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin0.fcl
new file mode 100644
index 000000000..7bfa4bd01
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin0.fcl
@@ -0,0 +1,13 @@
+# Fit field responses, turned up noise on each plane
+#
+#include "detsim_2d_icarus.fcl"
+
+# Add in the tuned field responses
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+# variation on ind 0
+physics.producers.daq.wcls_main.params.files_fields: "icarus_fnal_fit_ks_ind1bin0.json.bz2"
+
+# Set gains to MCNuPhase2
+physics.producers.daq.wcls_main.structs.gain0: 10.2636323
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin14.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin14.fcl
new file mode 100644
index 000000000..192ba6ced
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_ind1bin14.fcl
@@ -0,0 +1,13 @@
+# Fit field responses, turned up noise on each plane
+#
+#include "detsim_2d_icarus.fcl"
+
+# Add in the tuned field responses
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+# variation on ind 0
+physics.producers.daq.wcls_main.params.files_fields: "icarus_fnal_fit_ks_ind1bin14.json.bz2"
+
+# Set gains to MCNuPhase2
+physics.producers.daq.wcls_main.structs.gain0: 10.2636323
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-1x.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-1x.fcl
new file mode 100644
index 000000000..266577e83
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-1x.fcl
@@ -0,0 +1,11 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 10.2636323
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
+
+physics.producers.daq.wcls_main.structs.int_noise_scale: 1.10
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-2x.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-2x.fcl
new file mode 100644
index 000000000..32e048ea7
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_noise1-2x.fcl
@@ -0,0 +1,11 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 10.2636323
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
+
+physics.producers.daq.wcls_main.structs.int_noise_scale: 1.20
diff --git a/fcl/numi-anaA/detsim_2d_icarus_fitFR_oldgain.fcl b/fcl/numi-anaA/detsim_2d_icarus_fitFR_oldgain.fcl
new file mode 100644
index 000000000..fa40ba0db
--- /dev/null
+++ b/fcl/numi-anaA/detsim_2d_icarus_fitFR_oldgain.fcl
@@ -0,0 +1,9 @@
+# Fit field responses, same gains as used in the MCNuPhase2 2023A monte carlo
+#
+#include "detsim_2d_icarus.fcl"
+
+
+physics.producers.daq: @local::icarus_simwire_wirecell_fitSR
+
+physics.producers.daq.wcls_main.structs.gain0: 10.2636323
+physics.producers.daq.wcls_main.structs.gain1: 12.1420344
diff --git a/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl b/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl
index 64734aa13..ebff398e2 100644
--- a/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl
+++ b/icaruscode/TPC/ICARUSWireCell/detsimmodules_wirecell_ICARUS.fcl
@@ -70,6 +70,7 @@ icarus_simwire_wirecell:
     }
 }
 
+# TUNED FR CONFIG
 icarus_simwire_wirecell_fitSR: @local::icarus_simwire_wirecell
 # Add in the ER tail
 icarus_simwire_wirecell_fitSR.wcls_main.params.file_rcresp: "icarus_fnal_rc_tail.json"
@@ -80,4 +81,16 @@ icarus_simwire_wirecell_fitSR.wcls_main.structs.gain0: 11.9918701 # mV/fC
 icarus_simwire_wirecell_fitSR.wcls_main.structs.gain1: 12.6181926 # mV/fC
 icarus_simwire_wirecell_fitSR.wcls_main.structs.gain2: 13.0261362 # mV/fC
 icarus_simwire_wirecell_fitSR.wcls_main.structs.shaping1: 1.45 # us
+
+# TUNED FR CONFIG w/ Front Induction Nominal
+icarus_simwire_wirecell_fitSR_P0nom: @local::icarus_simwire_wirecell
+# Add in the ER tail
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.params.file_rcresp: "icarus_fnal_rc_tail.json"
+# Add in the tuned field responses (except for P0)
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.params.files_fields: "icarus_fnal_fit_ks_P0nom.json.bz2"
+# futz with shaping+gain values (note these are really just scale factors and should not be taken literally)
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.structs.gain0: 17.05212 # mV/fC 
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.structs.gain1: 12.6181926 # mV/fC
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.structs.gain2: 13.0261362 # mV/fC
+icarus_simwire_wirecell_fitSR_P0nom.wcls_main.structs.shaping1: 1.45 # us
 END_PROLOG