Merge pull request #205 from paolafer/g4-11-1

Use total drift length of EL gap to generate EL photons
next-exp · Apr 26, 2023 · 4960b3b · 4960b3b
2 parents 8ae6374 + 5f504f6
commit 4960b3b
Show file tree

Hide file tree

Showing 4 changed files with 13 additions and 1 deletion.
diff --git a/source/physics/BaseDriftField.h b/source/physics/BaseDriftField.h
@@ -40,6 +40,7 @@ namespace nexus {
       GeneratePointAlongDriftLine(const G4LorentzVector&, const G4LorentzVector&) = 0;
 
     virtual G4double LightYield() const;
+    virtual G4double GetTotalDriftLength() const;
 
   private:
     void Print() const;
@@ -53,6 +54,8 @@ namespace nexus {
 
   inline G4double BaseDriftField::LightYield() const {return 0.;}
 
+  inline G4double BaseDriftField::GetTotalDriftLength() const {return 0.;}
+
   inline void BaseDriftField::Print() const {}
 
 } // end namespace nexus

diff --git a/source/physics/Electroluminescence.cc b/source/physics/Electroluminescence.cc
@@ -79,7 +79,7 @@ Electroluminescence::PostStepDoIt(const G4Track& track, const G4Step& step)
 
   // Get the light yield from the field
   const G4double yield = field->LightYield();
-  G4double step_length = step.GetStepLength();
+  G4double step_length = field->GetTotalDriftLength();
 
   if (yield <= 0.)
     return G4VDiscreteProcess::PostStepDoIt(track, step);

diff --git a/source/physics/UniformElectricDriftField.cc b/source/physics/UniformElectricDriftField.cc
@@ -115,4 +115,11 @@ namespace nexus {
   }
 
 
+
+  G4double UniformElectricDriftField::GetTotalDriftLength() const
+  {
+    return std::abs(anode_pos_ - cathode_pos_);
+  }
+
+
 } // end namespace nexus
diff --git a/source/physics/UniformElectricDriftField.h b/source/physics/UniformElectricDriftField.h
@@ -64,6 +64,8 @@ namespace nexus {
     void SetNumberOfPhotons(G4double);
     G4double GetNumberOfPhotons() const;
 
+    virtual G4double GetTotalDriftLength() const;
+
   private:
     /// Returns true if coordinate is between anode and cathode
     G4bool CheckCoordinate(G4double);