Changeset 1315 for trunk/source/processes/electromagnetic/xrays/include/G4Scintillation.hh

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/processes/electromagnetic/xrays/include/G4Scintillation.hh (modified) (6 diffs)

trunk/source/processes/electromagnetic/xrays/include/G4Scintillation.hh

@@ -25 +25 @@
 //
 //
-// $Id: G4Scintillation.hh,v 1.16 2009/07/29 23:45:20 gum Exp $
-// GEANT4 tag $Name: geant4-09-03 $
+// $Id: G4Scintillation.hh,v 1.17 2010/05/27 20:48:35 gum Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 // 
@@ -146 +146 @@
         // have been tracked, the tracking of the primary resumes.
 
+        void SetFiniteRiseTime(const G4bool state);
+        // If set, the G4Scintillation process expects the user to have
+        // set the constant material property FAST/SLOWSCINTILLATIONRISETIME.
+
         G4bool GetTrackSecondariesFirst() const;
         // Returns the boolean flag for tracking secondaries first.
+
+        G4bool GetFiniteRiseTime() const;
+        // Returns the boolean flag for a finite scintillation rise time.
         
         void SetScintillationYieldFactor(const G4double yieldfactor);
@@ -198 +205 @@
 
         G4bool fTrackSecondariesFirst;
+        G4bool fFiniteRiseTime;
 
         G4double YieldFactor;
@@ -204 +212 @@
 
 private:
+
+        G4double single_exp(G4double t, G4double tau2);
+        G4double bi_exp(G4double t, G4double tau1, G4double tau2);
+
+        // emission time distribution when there is a finite rise time
+        G4double sample_time(G4double tau1, G4double tau2);
 
         G4EmSaturation* emSaturation;
@@ -229 +243 @@
 
 inline
+void G4Scintillation::SetFiniteRiseTime(const G4bool state)
+{
+        fFiniteRiseTime = state;
+}
+
+inline
 G4bool G4Scintillation::GetTrackSecondariesFirst() const
 {
         return fTrackSecondariesFirst;
+}
+
+inline 
+G4bool G4Scintillation::GetFiniteRiseTime() const
+{
+        return fFiniteRiseTime;
 }
 
@@ -296 +322 @@
 }
 
+inline
+G4double G4Scintillation::single_exp(G4double t, G4double tau2)
+{
+         return exp(-1.0*t/tau2)/tau2;
+}
+
+inline
+G4double G4Scintillation::bi_exp(G4double t, G4double tau1, G4double tau2)
+{
+         return exp(-1.0*t/tau2)*(1-exp(-1.0*t/tau1))/tau2/tau2*(tau1+tau2);
+}
+
 #endif /* G4Scintillation_h */