Macros | Functions | Variables
G4S2Light.cc File Reference
#include "G4ParticleTypes.hh"
#include "G4EmProcessSubType.hh"
#include "G4S2Light.hh"
#include "G4S1Light.hh"
#include "G4SystemOfUnits.hh"
#include "G4PhysicalConstants.hh"

Go to the source code of this file.

Macros

#define E_PURITY   2*CLHEP::m
 
#define GRID_DENSITY   8.03*(CLHEP::g/CLHEP::cm3)
 

Functions

G4double GetGasElectronDriftSpeed (G4double efieldinput, G4double density)
 
G4double GetLiquidElectronDriftSpeed (double T, double F, G4bool M, G4int Z)
 

Variables

G4double thr [100]
 
G4double E_eV [100]
 
G4double tau1 [100]
 
G4double tau3 [100]
 
G4double MolarMass [100]
 
G4double ConvertEff [100]
 

Macro Definition Documentation

#define E_PURITY   2*CLHEP::m

Definition at line 56 of file G4S2Light.cc.

#define GRID_DENSITY   8.03*(CLHEP::g/CLHEP::cm3)

Definition at line 57 of file G4S2Light.cc.

Function Documentation

G4double GetGasElectronDriftSpeed ( G4double  efieldinput,
G4double  density 
)

Definition at line 315 of file G4S2Light.cc.

316 {
317  //Gas equation one coefficients(E/N of 1.2E-19 to 3.5E-19)
318  double gas1a=395.50266631436,gas1b=-357384143.004642,gas1c=0.518110447340587;
319  //Gas equation two coefficients(E/N of 3.5E-19 to 3.8E-17)
320  double gas2a=-592981.611357632,gas2b=-90261.9643716643,
321  gas2c=-4911.83213989609,gas2d=-115.157545835228, gas2f=-0.990440443390298,
322  gas2g=1008.30998933704,gas2h=223.711221224885;
323 
324  G4double edrift=0, gasdep=efieldinput/density, gas1fix=0, gas2fix=0;
325 
326  if ( gasdep < 1.2e-19 && gasdep >= 0 ) edrift = 4e22*gasdep;
327  if ( gasdep < 3.5e-19 && gasdep >= 1.2e-19 ) {
328  gas1fix = gas1b*pow(gasdep,gas1c); edrift = gas1a*pow(gasdep,gas1fix);
329  }
330  if ( gasdep < 3.8e-17 && gasdep >= 3.5e-19 ) {
331  gas2fix = log(gas2g*gasdep);
332  edrift = (gas2a+gas2b*gas2fix+gas2c*pow(gas2fix,2)+gas2d*pow(gas2fix,3)+
333  gas2f*pow(gas2fix,4))*(gas2h*exp(gasdep));
334  }
335  if ( gasdep >= 3.8e-17 ) edrift = 6e21*gasdep-32279;
336 
337  return 0.5*EMASS*pow(edrift*(CLHEP::cm/CLHEP::s),2.);
338 }
static constexpr double cm
Definition: Units.h:68
#define EMASS
Definition: G4S1Light.hh:25
constexpr T pow(T x)
Definition: pow.h:72
const double e
static QCString * s
Definition: config.cpp:1042
G4double GetLiquidElectronDriftSpeed ( double  T,
double  F,
G4bool  M,
G4int  Z 
)

Variable Documentation

G4double ConvertEff[100]

Definition at line 61 of file G4S2Light.cc.

G4double E_eV[100]

Definition at line 60 of file G4S2Light.cc.

G4double MolarMass[100]

Definition at line 61 of file G4S2Light.cc.

G4double tau1[100]

Definition at line 61 of file G4S2Light.cc.

G4double tau3[100]

Definition at line 61 of file G4S2Light.cc.

G4double thr[100]

Definition at line 59 of file G4S2Light.cc.