From ecec9e94eaa9b73dd7a572b41a72a7df76d8ea51 Mon Sep 17 00:00:00 2001
From: Richard Jones <rjones30@gmail.com>
Date: Thu, 29 Jul 2021 18:07:32 -0400
Subject: [PATCH] * extend GlueXBernardConversionProcess.cc to support nuclear
 targets,   like the 209Pb target to be used in CPP. [rtj]

---
 src/GlueXBernardConversionProcess.cc | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/src/GlueXBernardConversionProcess.cc b/src/GlueXBernardConversionProcess.cc
index e93e521..cbd2632 100644
--- a/src/GlueXBernardConversionProcess.cc
+++ b/src/GlueXBernardConversionProcess.cc
@@ -28,6 +28,7 @@
 // pairs inside, otherwise the standard pair conversion
 // probabilities apply.
 #define FORCED_LIH2_PAIR_CONVERSION 0
+#define FORCED_TGT0_PAIR_CONVERSION 0
 
 #include <G4SystemOfUnits.hh>
 #include "G4PhysicalConstants.hh"
@@ -204,6 +205,15 @@ G4double GlueXBernardConversionProcess::PostStepGetPhysicalInteractionLength(
       *condition = Forced;
       return 100*cm;
    }
+   else if (track.GetTrackID() == 1 && pvol && pvol->GetName() == "TGT0" &&
+       (FORCED_TGT0_PAIR_CONVERSION || 
+        fStopBeamAfterTarget ))
+   {
+      fPIL = G4VEmProcess::PostStepGetPhysicalInteractionLength(
+                             track, previousStepSize, condition);
+      *condition = Forced;
+      return 100*cm;
+   }
    *condition = NotForced;
    return 1e99;
 }