LDMX-Software · tomeichlersmith · Aug 18, 2023 · Mar 27, 2023 · May 18, 2023 · May 22, 2023
diff --git a/G4DarkBreM b/G4DarkBreM
diff --git a/include/SimCore/UserEventInformation.h b/include/SimCore/UserEventInformation.h
@@ -24,6 +24,23 @@ class UserEventInformation : public G4VUserEventInformation {
   /// Decrease the number of brem candidates in an event.
   void decBremCandidateCount() { bremCandidateCount_ -= 1; }
 
+  /**
+   * Set the Z of the element in which the dark brem ocurred
+   *
+   * @param[in] z atomic Z of element in which the dark brem ocurred
+   */
+  void setDarkBremMaterialZ(double z) { db_material_z_ = z; }
+
+  /**
+   * Get the Z of the element in which the dark brem ocurred
+   *
+   * @note This will return -1 if no dark brem ocurred within
+   * this event.
+   *
+   * @param[in] z atomic Z of element in which the dark brem ocurred
+   */
+  double getDarkBremMaterialZ() const { return db_material_z_; }
+
   /**
    * Set the event weight.
    *
@@ -140,6 +157,14 @@ class UserEventInformation : public G4VUserEventInformation {
    * Was the most recent step a electron-nuclear interaction?
    */
   bool last_step_en_{false};
+
+  /**
+   * atomic Z of the element in which dark brem occurred (-1 if didn't happen)
+   *
+   * The default is -1. and so will provide unphysical results if the
+   * dark brem did not occur within the event in question.
+   */
+  double db_material_z_{-1.};
 };
 }  // namespace simcore
 

diff --git a/src/SimCore/APrimePhysics.cxx b/src/SimCore/APrimePhysics.cxx
@@ -9,16 +9,34 @@
 
 #include "G4DarkBreM/G4APrime.h"
 #include "G4DarkBreM/G4DarkBreMModel.h"
+#include "SimCore/UserEventInformation.h"
 
 // Geant4
 #include "G4Electron.hh"
+#include "G4EventManager.hh"
 #include "G4ProcessManager.hh"
 
 namespace simcore {
 
 const std::string APrimePhysics::NAME = "APrime";
 
-APrimePhysics::APrimePhysics(const framework::config::Parameters &params)
+/**
+ * Store the atomic Z for the element in which the dark brem occurred.
+ *
+ * This function is registered with G4DarkBremsstrahlung and will be
+ * called with the element that the dark brem will occurr off of. We
+ * store the element Z in UserEventInformation for later serialization
+ * into the EventHeader
+ *
+ * @param[in] element G4Element off-which the dark brem occurred
+ */
+static void store_element_z(const G4Element& element) {
+  static_cast<UserEventInformation*>(
+      G4EventManager::GetEventManager()->GetUserInformation())
+      ->setDarkBremMaterialZ(element.GetZ());
+}
+
+APrimePhysics::APrimePhysics(const framework::config::Parameters& params)
     : G4VPhysicsConstructor(APrimePhysics::NAME),
       parameters_{params},
       process_{nullptr} {
@@ -74,6 +92,7 @@ void APrimePhysics::ConstructProcess() {
           parameters_.getParameter<bool>("only_one_per_event"),
           1., /* global bias - should use bias operator instead */
           parameters_.getParameter<bool>("cache_xsec"));
+      process_->RegisterStorageMechanism(store_element_z);
     } else {
       EXCEPTION_RAISE("BadConf",
                       "Unrecognized model name '" + model_name + "'.");

diff --git a/src/SimCore/DetectorConstruction.cxx b/src/SimCore/DetectorConstruction.cxx
@@ -26,7 +26,8 @@ static bool isInEcal(G4LogicalVolume* vol, const std::string& vol_to_bias) {
            volumeName.contains("PCB") || volumeName.contains("strongback") ||
            volumeName.contains("Glue") || volumeName.contains("CFMix") ||
            volumeName.contains("Al") || volumeName.contains("C")) &&
-          volumeName.contains("volume"));
+          volumeName.contains("volume")) ||
+         (volumeName.contains("nohole_motherboard"));
 }
 
 /**

diff --git a/src/SimCore/SimulatorBase.cxx b/src/SimCore/SimulatorBase.cxx
@@ -24,6 +24,8 @@ void SimulatorBase::updateEventHeader(ldmx::EventHeader& eventHeader) const {
                                 event_info->getPNEnergy());
   eventHeader.setFloatParameter("total_electronuclear_energy",
                                 event_info->getENEnergy());
+  eventHeader.setFloatParameter("db_material_z",
+                                event_info->getDarkBremMaterialZ());
 }
 void SimulatorBase::onProcessEnd() {
   runManager_->TerminateEventLoop();
+2 −1		CMakeLists.txt
+1 −1		app/extract_library.cxx
+10 −2		app/scale.cxx
+2 −2		app/simulate.cxx
+53 −19		app/xsec_calc.cxx
+177 −0		include/G4DarkBreM/ElementXsecInterpolation.h
+11 −7		include/G4DarkBreM/G4DarkBreMModel.h
+38 −1		include/G4DarkBreM/G4DarkBremsstrahlung.h
+25 −14		include/G4DarkBreM/ParseLibrary.h
+3 −1		include/G4DarkBreM/PrototypeModel.h
+154 −0		src/G4DarkBreM/ElementXsecInterpolation.cxx
+45 −23		src/G4DarkBreM/G4DarkBreMModel.cxx
+53 −6		src/G4DarkBreM/G4DarkBremsstrahlung.cxx
+61 −33		src/G4DarkBreM/ParseLibrary.cxx