Add a model that prevents PN reactions

Author: EinarElen

include/SimCore/PhotoNuclearModels/NoPhotoNuclearModel.h

@@ -0,0 +1,22 @@
+#ifndef NOPHOTONUCLEARMODEL_H
+#define NOPHOTONUCLEARMODEL_H
+
+#include <G4ProcessManager.hh>
+
+#include "Framework/Configure/Parameters.h"
+#include "SimCore/PhotoNuclearModel.h"
+
+namespace simcore {
+
+class NoPhotoNuclearModel : public PhotoNuclearModel {
+ public:
+  NoPhotoNuclearModel(const std::string& name,
+                      const framework::config::Parameters& parameters)
+      : PhotoNuclearModel{name, parameters} {}
+  virtual ~NoPhotoNuclearModel() = default;
+  void ConstructGammaProcess(G4ProcessManager* processManager) override;
+};
+
+}  // namespace simcore
+
+#endif /* NOPHOTONUCLEARMODEL_H */

python/photonuclear_models.py

@@ -125,3 +125,16 @@ def kaon(min_products = 1, hard_particle_threshold=200.):
         ]
         model.min_products = min_products
         return model
+
+
+class NoPhotoNuclearModel(simcfg.PhotoNuclearModel):
+    """A PhotoNuclear model that disables the photonuclear process entirely.
+
+    Make sure that no biasing operators for photonuclear reactions are enabled
+    when using this model.
+
+    """
+    def __init__(self):
+        super().__init__('NoPhotoNuclearModel',
+                         'simcore::NoPhotoNuclearModel',
+                         'SimCore_PhotoNuclearModels')

src/SimCore/PhotoNuclearModels/NoPhotoNuclearModel.cxx

@@ -0,0 +1,11 @@
+#include "SimCore/PhotoNuclearModels/NoPhotoNuclearModel.h"
+
+namespace simcore {
+
+void NoPhotoNuclearModel::ConstructGammaProcess(
+    G4ProcessManager* processManager) {
+  // Do nothing
+}
+}  // namespace simcore
+
+DECLARE_PHOTONUCLEAR_MODEL(simcore::NoPhotoNuclearModel);