PEN: add resolutionscale

Without this Geant4 will abort when trying to scintillate in a PEN volume.
legend-exp · May 7, 2024 · 78b6cc4 · 78b6cc4
1 parent eede31f
commit 78b6cc4
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/legendoptics/pen.py b/src/legendoptics/pen.py
@@ -189,6 +189,10 @@ def pyg4_pen_attach_scintillation(mat, reg) -> None:
 
     mat.addConstPropertyPint("SCINTILLATIONTIMECONSTANT1", pen_scint_timeconstant())
 
+    # We do not know the PEN fano factor. Geant4 calculates σ = RESOLUTIONSCALE × √mean,
+    # so we set it to 1; otherwise Geant4 will crash.
+    mat.addConstPropertyPint("RESOLUTIONSCALE", 1)
+
     # setup scintillation response just for electrons.
     scint_config = ScintConfig(
         flat_top=pen_scint_light_yield(),