From 24caa3967b7b34f329972992cc3bee5a2e3cf009 Mon Sep 17 00:00:00 2001
From: michaelmackenzie <michaelmackenzie@u.northwestern.edu>
Date: Tue, 3 Feb 2026 14:15:01 -0600
Subject: [PATCH 1/2] Add fcl to simulate the 1.8 MeV muon capture photon

---
 JobConfig/primary/MuCap1809keV.fcl     | 24 ++++++++++++++++++++++++
 JobConfig/primary/MuCap1809keVCalo.fcl | 10 ++++++++++
 Tests/MuCap1809keVCaloSteps.fcl        |  4 ++++
 Tests/MuCap1809keVSteps.fcl            |  4 ++++
 4 files changed, 42 insertions(+)
 create mode 100644 JobConfig/primary/MuCap1809keV.fcl
 create mode 100644 JobConfig/primary/MuCap1809keVCalo.fcl
 create mode 100644 Tests/MuCap1809keVCaloSteps.fcl
 create mode 100644 Tests/MuCap1809keVSteps.fcl

diff --git a/JobConfig/primary/MuCap1809keV.fcl b/JobConfig/primary/MuCap1809keV.fcl
new file mode 100644
index 00000000..6530a88e
--- /dev/null
+++ b/JobConfig/primary/MuCap1809keV.fcl
@@ -0,0 +1,24 @@
+#
+# Generate 1809 keV muon capture photons
+#
+# Original author: Michael MacKenzie, 2026
+#
+#include "Production/JobConfig/primary/TargetStopParticle.fcl"
+
+physics.producers.generate : {
+  module_type : SingleProcessGenerator
+  inputSimParticles: TargetStopResampler
+  stoppingTargetMaterial : "Al"
+  decayProducts :
+  {
+    tool_type : MuCap1809keVGammaGenerator
+    fireAll   : true # Insert a photon into every event
+  }
+  verbosity : 0
+}
+
+physics.filters.PrimaryFilter.MinimumPartMom : 0.1
+physics.filters.PrimaryFilter.MinimumSumCaloStepE : 0.1
+
+physics.producers.FindMCPrimary.PrimaryProcess : "mu2eMuonCaptureAtRest"
+outputs.PrimaryOutput.fileName: "dts.owner.MuCap1809keV.version.sequencer.art"
diff --git a/JobConfig/primary/MuCap1809keVCalo.fcl b/JobConfig/primary/MuCap1809keVCalo.fcl
new file mode 100644
index 00000000..f6feb0c0
--- /dev/null
+++ b/JobConfig/primary/MuCap1809keVCalo.fcl
@@ -0,0 +1,10 @@
+#
+# Generate 1809 keV muon capture photons, aimed at the calorimeter
+#
+# Original author: Michael MacKenzie, 2026
+#
+#include "Production/JobConfig/primary/MuCap1809keV.fcl"
+
+physics.producers.generate.decayProducts.czmin : 0.992
+physics.producers.generate.decayProducts.czmax : 1.
+outputs.PrimaryOutput.fileName: "dts.owner.MuCap1809keVCalo.version.sequencer.art"
diff --git a/Tests/MuCap1809keVCaloSteps.fcl b/Tests/MuCap1809keVCaloSteps.fcl
new file mode 100644
index 00000000..9552afc9
--- /dev/null
+++ b/Tests/MuCap1809keVCaloSteps.fcl
@@ -0,0 +1,4 @@
+#include "Production/JobConfig/primary/MuCap1809keVCalo.fcl"
+#include "Production/Tests/MuonStopConfig.fcl"
+
+services.TFileService.fileName: "nts.owner.test.version.sequencer.root"
diff --git a/Tests/MuCap1809keVSteps.fcl b/Tests/MuCap1809keVSteps.fcl
new file mode 100644
index 00000000..2e9e83ce
--- /dev/null
+++ b/Tests/MuCap1809keVSteps.fcl
@@ -0,0 +1,4 @@
+#include "Production/JobConfig/primary/MuCap1809keV.fcl"
+#include "Production/Tests/MuonStopConfig.fcl"
+
+services.TFileService.fileName: "nts.owner.test.version.sequencer.root"

From 7bcd0c64b6d64be6554f763811570de906dc260d Mon Sep 17 00:00:00 2001
From: michaelmackenzie <michaelmackenzie@u.northwestern.edu>
Date: Wed, 4 Feb 2026 08:53:33 -0600
Subject: [PATCH 2/2] Make czmin a bit more conservative

---
 JobConfig/primary/MuCap1809keVCalo.fcl | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/JobConfig/primary/MuCap1809keVCalo.fcl b/JobConfig/primary/MuCap1809keVCalo.fcl
index f6feb0c0..f3ac094e 100644
--- a/JobConfig/primary/MuCap1809keVCalo.fcl
+++ b/JobConfig/primary/MuCap1809keVCalo.fcl
@@ -5,6 +5,6 @@
 #
 #include "Production/JobConfig/primary/MuCap1809keV.fcl"
 
-physics.producers.generate.decayProducts.czmin : 0.992
+physics.producers.generate.decayProducts.czmin : 0.99
 physics.producers.generate.decayProducts.czmax : 1.
 outputs.PrimaryOutput.fileName: "dts.owner.MuCap1809keVCalo.version.sequencer.art"