From 981b9cb744e4739fd24af6c0ff4fc37a876ed9f9 Mon Sep 17 00:00:00 2001
From: michaelmackenzie <michaelmackenzie@u.northwestern.edu>
Date: Fri, 10 Jul 2026 13:55:03 -0500
Subject: [PATCH 1/2] Add neutron tail primaries for calo cluster studies, add
 missing product line

---
 JobConfig/primary/MuCapNeutronTail.fcl     | 26 ++++++++++++++++++++++
 JobConfig/primary/MuCapNeutronTailCalo.fcl | 10 +++++++++
 JobConfig/primary/prolog.fcl               |  3 ++-
 Tests/MuCapNeutronTailSteps.fcl            |  4 ++++
 4 files changed, 42 insertions(+), 1 deletion(-)
 create mode 100644 JobConfig/primary/MuCapNeutronTail.fcl
 create mode 100644 JobConfig/primary/MuCapNeutronTailCalo.fcl
 create mode 100644 Tests/MuCapNeutronTailSteps.fcl

diff --git a/JobConfig/primary/MuCapNeutronTail.fcl b/JobConfig/primary/MuCapNeutronTail.fcl
new file mode 100644
index 00000000..8c65b4fe
--- /dev/null
+++ b/JobConfig/primary/MuCapNeutronTail.fcl
@@ -0,0 +1,26 @@
+#
+# Generate high KE muon capture neutrons
+#
+# Original author: Michael MacKenzie, 2026
+#
+#include "Production/JobConfig/primary/TargetStopParticle.fcl"
+
+physics.producers.generate : {
+  module_type : SingleProcessGenerator
+  inputSimParticles: TargetStopResampler
+  stoppingTargetMaterial : "Al"
+  decayProducts :
+  {
+    tool_type : MuCapNeutronGenerator
+    spectrumVariable : kineticEnergy
+    spectrum : {
+      spectrumShape: tabulated
+      spectrumFileName: "Offline/EventGenerator/data/SchroederNeutronSpectrum.txt"
+      elow : 60.
+    }
+  }
+  verbosity : 0
+}
+
+physics.producers.FindMCPrimary.PrimaryProcess : "mu2eMuonCaptureAtRest"
+outputs.PrimaryOutput.fileName: "dts.owner.MuCapNeutronTail.version.sequencer.art"
diff --git a/JobConfig/primary/MuCapNeutronTailCalo.fcl b/JobConfig/primary/MuCapNeutronTailCalo.fcl
new file mode 100644
index 00000000..2aeb1843
--- /dev/null
+++ b/JobConfig/primary/MuCapNeutronTailCalo.fcl
@@ -0,0 +1,10 @@
+#
+# Generate high energy muon capture neutrons, aimed at the calorimeter
+#
+# Original author: Michael MacKenzie, 2026
+#
+#include "Production/JobConfig/primary/MuCapNeutronTail.fcl"
+
+physics.producers.generate.decayProducts.czmin : 0.99
+physics.producers.generate.decayProducts.czmax : 1.
+outputs.PrimaryOutput.fileName: "dts.owner.MuCapNeutronCalo.version.sequencer.art"
diff --git a/JobConfig/primary/prolog.fcl b/JobConfig/primary/prolog.fcl
index eeedf0c1..9700a68e 100644
--- a/JobConfig/primary/prolog.fcl
+++ b/JobConfig/primary/prolog.fcl
@@ -188,7 +188,8 @@ Primary: {
     "keep mu2e::EventWeight_*_*_*",
     "keep art::TriggerResults_*_*_*",
     "keep mu2e::StepPointMCs_CosmicResampler_crvStage1_*", # this is needed to keep track of low/all CRV split
-    "keep mu2e::PrimaryParticle_*_*_*"
+    "keep mu2e::PrimaryParticle_*_*_*",
+    "keep mu2e::SpectrumConfig_*_*_*"
   ]
 
   StepProducts : [
diff --git a/Tests/MuCapNeutronTailSteps.fcl b/Tests/MuCapNeutronTailSteps.fcl
new file mode 100644
index 00000000..1d5c6a16
--- /dev/null
+++ b/Tests/MuCapNeutronTailSteps.fcl
@@ -0,0 +1,4 @@
+#include "Production/JobConfig/primary/MuCapNeutronTailCalo.fcl"
+#include "Production/Tests/MuonStopConfig.fcl"
+
+services.TFileService.fileName: "nts.owner.test.version.sequencer.root"

From 33b832117ac281139473ad70a0913e67fe731fa7 Mon Sep 17 00:00:00 2001
From: michaelmackenzie <michaelmackenzie@u.northwestern.edu>
Date: Fri, 10 Jul 2026 13:58:23 -0500
Subject: [PATCH 2/2] Fix test fcl

---
 Tests/MuCapNeutronTailSteps.fcl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/Tests/MuCapNeutronTailSteps.fcl b/Tests/MuCapNeutronTailSteps.fcl
index 1d5c6a16..0c5aa181 100644
--- a/Tests/MuCapNeutronTailSteps.fcl
+++ b/Tests/MuCapNeutronTailSteps.fcl
@@ -1,4 +1,4 @@
-#include "Production/JobConfig/primary/MuCapNeutronTailCalo.fcl"
+#include "Production/JobConfig/primary/MuCapNeutronTail.fcl"
 #include "Production/Tests/MuonStopConfig.fcl"
 
 services.TFileService.fileName: "nts.owner.test.version.sequencer.root"