Add Fano factor and ionization energy to GXe properties

next-exp · Sep 13, 2024 · e23d806 · e23d806
1 parent e4a157e
commit e23d806
Show file tree

Hide file tree

Showing 2 changed files with 12 additions and 12 deletions.
diff --git a/source/materials/MaterialProperties.cc b/source/materials/MaterialProperties.cc
@@ -647,12 +647,14 @@ namespace materialprops {
 
     // CONST PROPERTIES
     mpt->AddConstProperty("SCINTILLATIONYIELD", sc_yield);
-    mpt->AddConstProperty("RESOLUTIONSCALE",    1.0);
-    mpt->AddConstProperty("SCINTILLATIONTIMECONSTANT1",   4.5  * ns);
-    mpt->AddConstProperty("SCINTILLATIONTIMECONSTANT2",   100. * ns);
+    mpt->AddConstProperty("RESOLUTIONSCALE",         1.0);
+    mpt->AddConstProperty("SCINTILLATIONTIMECONSTANT1",   4.5 * ns);
+    mpt->AddConstProperty("SCINTILLATIONTIMECONSTANT2", 100.  * ns);
     mpt->AddConstProperty("SCINTILLATIONYIELD1", .1);
     mpt->AddConstProperty("SCINTILLATIONYIELD2", .9);
     mpt->AddConstProperty("ATTACHMENT",         e_lifetime, 1);
+    mpt->AddConstProperty("IONIZATIONENERGY",    22.4 * eV, 1);
+    mpt->AddConstProperty("FANOFACTOR",                .15, 1);
 
     return mpt;
   }

diff --git a/source/physics/IonizationClustering.cc b/source/physics/IonizationClustering.cc
@@ -109,15 +109,13 @@ namespace nexus {
     // and N is the average number of charges.
 
     // Fetch the W_i and F from the material properties table
-    //G4MaterialPropertiesTable* mpt =
-    //  track.GetMaterial()->GetMaterialPropertiesTable();
-    //if (!mpt)
-    //  return G4VRestDiscreteProcess::PostStepDoIt(track, step);
-
-    //G4double ioni_energy = mpt->GetConstProperty("IONIZATIONENERGY");
-    //G4double fano_factor = mpt->GetConstProperty("FANOFACTOR");
-    G4double ioni_energy = 22.4 * eV;
-    G4double fano_factor = .15;
+    G4MaterialPropertiesTable* mpt =
+      track.GetMaterial()->GetMaterialPropertiesTable();
+    if (!mpt)
+      return G4VRestDiscreteProcess::PostStepDoIt(track, step);
+
+    G4double ioni_energy = mpt->GetConstProperty("IONIZATIONENERGY");
+    G4double fano_factor = mpt->GetConstProperty("FANOFACTOR");
 
     G4double mean = energy_dep / ioni_energy;