diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..3465aed69
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.125
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..19d090358
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.140
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..8eda885b0
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.035
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..b5c5f4877
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.050
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..d4718fa5d
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.075
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..45d116d78
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.140
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..39ae4e5fb
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.165
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..bce3da7c8
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.190
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..3536a27c1
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.215
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..2e88040c2
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.240
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..0deec68eb
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodMeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodMeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.260
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..b86268374
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_125_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.125
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..2c0c188ea
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_140_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.140
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..63aeced23
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_35_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.035
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..9d0169841
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_50_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.050
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..c76c05dc1
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nuee_m_75_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.075
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_e_e"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..e7aa5ccfb
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_140_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.140
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..9a3ded599
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_165_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.165
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..cb8188429
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_190_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.190
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..166886b9d
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_215_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.215
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..e924fe2b6
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_240_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.240
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"
diff --git a/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl
new file mode 100644
index 000000000..0c823b8da
--- /dev/null
+++ b/sbndcode/JobConfigurations/standard/gen/MeVPrtl/prodoverlay_corsika_MeVPrtl_hnl_nupi0_m_260_Um4_1p00e7_TPC_sbnd.fcl
@@ -0,0 +1,24 @@
+#include "prodoverlay_corsika_MeVPrtl_hnl_sbnd.fcl"
+
+###----------Set mass of HNL [GeV]----------###
+hnlM: 0.260
+physics.producers.generator.Flux.M: @local::hnlM
+physics.producers.generator.RayTrace.ReferencePrtlMass: @local::hnlM
+physics.producers.generator.Decay.ReferenceHNLMass: @local::hnlM
+
+###----------Set mixing angle of HNL ----------###
+Ume4: 0
+Umu4: 1.00e-7
+Umt4: 0
+physics.producers.generator.Flux.MagUe4: @local::Ume4
+physics.producers.generator.Flux.MagUm4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUE4: @local::Ume4
+physics.producers.generator.Decay.ReferenceUM4: @local::Umu4
+physics.producers.generator.Decay.ReferenceUT4: @local::Umt4
+physics.producers.generator.Decay.Decays: ["nu_pi0"]
+
+###-----------Verbose--------------###
+#include "set_MeVPrtl_Verbose.fcl" 
+                                                                      
+###----------Set geometry ----------###
+#include "set_MeVPrtl_GeoTPC.fcl"