clean up configurations to only keep cases that have implementations in geant4reweight

Author: pgreen135

sbncode/SBNEventWeight/jobs/geant4/eventweight_geant4_sbn.fcl

@@ -61,11 +61,11 @@ sbn_eventweight_geant4: {
     material: @local::LAr
     random_seed: 60
     parameters: @local::ProtonParameters
-    mode: multisigma
-    number_of_sigmas: 3
-    #mode: multisim
-    #number_of_multisims: 1000
-    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root"
+    #mode: multisigma
+    #number_of_sigmas: 3
+    mode: multisim
+    number_of_multisims: 1000
+    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root" # dummy file, cross section pulled directly from Geant4
     makeoutputtree: true
     pdg_to_reweight: 2212
     debug: true
@@ -74,13 +74,13 @@ sbn_eventweight_geant4: {
   reinteractions_neutron: {
     type: Geant4
     material: @local::LAr
-    random_seed: 60
+    random_seed: 61
     parameters: @local::NeutronParameters
-    mode: multisigma
-    number_of_sigmas: 3
-    #mode: multisim
-    #number_of_multisims: 1000
-    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root"
+    #mode: multisigma
+    #number_of_sigmas: 3
+    mode: multisim
+    number_of_multisims: 1000
+    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root" # dummy file, cross section pulled directly from Geant4
     makeoutputtree: true
     pdg_to_reweight: 2112
     debug: true
@@ -89,12 +89,12 @@ sbn_eventweight_geant4: {
   reinteractions_kplus: {
     type: Geant4
     material: @local::LAr
-    random_seed: 60
+    random_seed: 62
     parameters: @local::KPlusParameters
     #mode: pm1sigma
     mode: multisim
     number_of_multisims: 1000
-    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root"
+    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root" # dummy file, cross section pulled directly from Geant4
     makeoutputtree: true
     pdg_to_reweight: 321
     debug: true
@@ -103,21 +103,19 @@ sbn_eventweight_geant4: {
   reinteractions_kminus: {
     type: Geant4
     material: @local::LAr
-    random_seed: 60
+    random_seed: 63
     parameters: @local::KMinusParameters
     #mode: pm1sigma
     mode: multisim
     number_of_multisims: 1000
-    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root"
+    fracsfile: "$SBNDATA_DIR/systematics/reint/g4_fracs_proton.root" # dummy file, cross section pulled directly from Geant4
     makeoutputtree: true
     pdg_to_reweight: -321
     debug: true
   }
 
 }
 
-# reinteractions_neutron
-
 ###
 # Reweighting parameters should be defined as
 #
@@ -139,15 +137,15 @@ sbn_eventweight_geant4: {
 #   "nominal" value in that case, just the set value)
 # - Sigma defines the 1-sigma range for multisims and unisims. Currently the
 #   code can only accept a single value for sigma, giving symmetric
-#   uncertainty bounds: nominal-sigma and nominal+sigma
+#   uncertainty bounds: nominal-sigma and nominal+sigma 
+#   [Only in use for fitting purposes, instead range defined through input files - Jul 2025]
 # - Cut must be one of the following:
-#      "reac" <- total inelastic scattering cross section
+#      "total" <- total scattering cross section, pulled from Geant4
 #      "abs" <- absorption cross section (pi+ and pi- only at the moment -- Oct 2019)
 #      "cex" <- charge exchange cross section (pi+ and pi- only at the moment -- Oct 2019)
 #      "dcex" <- double charge exchange cross section (pi+ and pi- only at the moment -- Oct 2019)
 #      "prod" <- pion production cross section (pi+ and pi- only at the moment -- Oct 2019)
 #      "inel" <- quasi-elastic inelastic scattering cross section (pi+ and pi- only at the moment -- Oct 2019)
-#      "elast" <- total elastic scattering cross section
 # - Name can be anything you like (but should uniquely identify that parameter/
 #   variation)
 ###

sbncode/SBNEventWeight/jobs/geant4/kminus_reweight_parameters.fcl

@@ -1,18 +1,12 @@
 BEGIN_PROLOG
 
-KMinusParameters: [
-  {
-    Name:       "fKMinusReac"
-    Cut:        "reac"
-    Range:      [10., 2005.]
-    Nominal:    1
-    Sigma:      0.2
-  },
+## KMinus available process keys: "total"
 
+KMinusParameters: [
   {
-    Name:       "fKMinusElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
+    Name:       "fKMinusTotal"
+    Cut:        "total"
+    Range:      [10.0, 2005.0]
     Nominal:    1.0
     Sigma:      0.2
   }

sbncode/SBNEventWeight/jobs/geant4/kplus_reweight_parameters.fcl

@@ -1,18 +1,12 @@
 BEGIN_PROLOG
 
-KPlusParameters: [
-  {
-    Name:       "fKPlusReac"
-    Cut:        "reac"
-    Range:      [10., 2005.]
-    Nominal:    1
-    Sigma:      0.2
-  },
+## KPlus available process keys: "total"
 
+KPlusParameters: [
   {
-    Name:       "fKPlusElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
+    Name:       "fKPlusTotal"
+    Cut:        "total"
+    Range:      [10.0, 2005.0]
     Nominal:    1.0
     Sigma:      0.2
   }

sbncode/SBNEventWeight/jobs/geant4/neutron_reweight_parameters.fcl

@@ -1,30 +1,15 @@
 BEGIN_PROLOG
 
-NeutronParameters: [
-  {
-    Name:       "fNeutronReac"
-    Cut:        "reac"
-    Range:      [10., 2005.]
-    Nominal:    1
-    Sigma:      0.2
-  },
-
-  {
-    Name:       "fNeutronElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
-    Nominal:    1.0
-    Sigma:      0.2
-  },
+## Neutron available process keys: "total"
 
+NeutronParameters: [ 
   {
     Name:       "fNeutronTotal"
     Cut:        "total"
-    Range:      [10.0, 2005.00]
+    Range:      [10.0, 2005.0]
     Nominal:    1.0
     Sigma:      0.2
   }
-
 ]
 
 END_PROLOG 

sbncode/SBNEventWeight/jobs/geant4/piminus_reweight_parameters.fcl

@@ -1,6 +1,6 @@
 BEGIN_PROLOG
 
-# From MicroBooNE, ubsim UBOONE_SUITE_v08_00_00_69
+## PiMinus available process keys: "abs", "cex", "prod", "dcex", "inel"
 
 PiMinusParameters: [ 
   {
@@ -52,15 +52,8 @@ PiMinusParameters: [
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
-  },
-
-  {
-    Name:       "fPiMinusElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
-    Nominal:    1.0
-    Sigma:      0.2
   }
+    
 ]
 
 END_PROLOG

sbncode/SBNEventWeight/jobs/geant4/piplus_reweight_parameters.fcl

@@ -1,13 +1,13 @@
 BEGIN_PROLOG
 
-# From MicroBooNE, ubsim UBOONE_SUITE_v08_00_00_69
+## PiPlus available process keys: "abs", "cex", "prod", "dcex", "inel"
 
 PiPlusParameters: [
   {
     Name:       "fPiPlusReacLow"
     Cut:        "reac"
     Range:      [10., 200.]
-    Nominal:    1
+    Nominal:    1.0
     Sigma:      0.2
   } ,
   {
@@ -52,15 +52,8 @@ PiPlusParameters: [
     Range:      [10.0, 2005.00]
     Nominal:    1.0
     Sigma:      0.2
-  },
-
-  {
-    Name:       "fPiPlusElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
-    Nominal:    1.0
-    Sigma:      0.2
   }
+
 ]
 
 END_PROLOG

sbncode/SBNEventWeight/jobs/geant4/proton_reweight_parameters.fcl

@@ -1,29 +1,13 @@
 BEGIN_PROLOG
 
-# From MicroBooNE, ubsim UBOONE_SUITE_v08_00_00_69
+## Proton available process keys: "total"
 
 ProtonParameters: [
-  {
-    Name:       "fProtonReac"
-    Cut:        "reac"
-    Range:      [10., 2005.]
-    Nominal:    1
-    Sigma:      0.2
-  },
-
-  {
-    Name:       "fProtonElast"
-    Cut:        "elast"
-    Range:      [10.0, 2005.00]
-    Nominal:    1.0
-    Sigma:      0.2
-  },
-
   {
     Name:       "fProtonTotal"
     Cut:        "total"
-    Range:      [10., 2005.]
-    Nominal:    1
+    Range:      [10.0, 2005.0]
+    Nominal:    1.0
     Sigma:      0.2
   }
 ]