#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.
#define E_PURITY 2*CLHEP::m |
#define GRID_DENSITY 8.03*(CLHEP::g/CLHEP::cm3) |
G4double GetGasElectronDriftSpeed |
( |
G4double |
efieldinput, |
|
|
G4double |
density |
|
) |
| |
Definition at line 315 of file G4S2Light.cc.
318 double gas1a=395.50266631436,gas1b=-357384143.004642,gas1c=0.518110447340587;
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;
324 G4double edrift=0, gasdep=efieldinput/density, gas1fix=0, gas2fix=0;
326 if ( gasdep < 1.2e-19 && gasdep >= 0 ) edrift = 4e22*gasdep;
327 if ( gasdep < 3.5e-19 && gasdep >= 1.2
e-19 ) {
328 gas1fix = gas1b*
pow(gasdep,gas1c); edrift = gas1a*
pow(gasdep,gas1fix);
330 if ( gasdep < 3.8e-17 && gasdep >= 3.5
e-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));
335 if ( gasdep >= 3.8
e-17 ) edrift = 6e21*gasdep-32279;
static constexpr double cm
G4double GetLiquidElectronDriftSpeed |
( |
double |
T, |
|
|
double |
F, |
|
|
G4bool |
M, |
|
|
G4int |
Z |
|
) |
| |