Macros | Functions | Variables
NestAlg.cxx File Reference
#include "Geant4/G4Ions.hh"
#include "Geant4/G4OpticalPhoton.hh"
#include "Geant4/G4VProcess.hh"
#include "larsim/LegacyLArG4/NestAlg.h"
#include "larsim/LegacyLArG4/G4ThermalElectron.hh"
#include "CLHEP/Random/RandGauss.h"
#include "CLHEP/Random/RandFlat.h"

Go to the source code of this file.

Macros

#define AVO   6.022e23
 
#define EMASS   9.109e-31*CLHEP::kg
 
#define MillerDriftSpeed   true
 
#define GASGAP   0.25*CLHEP::cm
 
#define BORDER   0*CLHEP::cm
 
#define QE_EFF   1
 
#define phe_per_e   1
 
#define WIN   0*CLHEP::mm
 
#define TOP   0
 
#define ANE   0
 
#define SRF   0
 
#define GAT   0
 
#define CTH   0
 
#define BOT   0
 
#define PMT   0
 
#define MIN_ENE   -1*CLHEP::eV
 
#define MAX_ENE   1.*CLHEP::TeV
 
#define HIENLIM   5*CLHEP::MeV
 
#define R_TOL   0.2*CLHEP::mm
 
#define R_MAX   1*CLHEP::km
 
#define Density_LXe   2.9
 
#define Density_LAr   1.393
 
#define Density_LNe   1.207
 
#define Density_LKr   2.413
 

Functions

G4double UnivScreenFunc (G4double E, G4double Z, G4double A)
 

Variables

G4bool diffusion = true
 
G4bool SinglePhase =false
 
G4bool ThomasImelTail =true
 
G4bool OutElectrons =true
 
G4double biExc = 0.77
 

Macro Definition Documentation

#define ANE   0

Definition at line 14 of file NestAlg.cxx.

#define AVO   6.022e23

Definition at line 1 of file NestAlg.cxx.

#define BORDER   0*CLHEP::cm

Definition at line 6 of file NestAlg.cxx.

#define BOT   0

Definition at line 18 of file NestAlg.cxx.

#define CTH   0

Definition at line 17 of file NestAlg.cxx.

#define Density_LAr   1.393

Definition at line 26 of file NestAlg.cxx.

#define Density_LKr   2.413

Definition at line 28 of file NestAlg.cxx.

#define Density_LNe   1.207

Definition at line 27 of file NestAlg.cxx.

#define Density_LXe   2.9

Definition at line 25 of file NestAlg.cxx.

#define EMASS   9.109e-31*CLHEP::kg

Definition at line 2 of file NestAlg.cxx.

#define GASGAP   0.25*CLHEP::cm

Definition at line 5 of file NestAlg.cxx.

#define GAT   0

Definition at line 16 of file NestAlg.cxx.

#define HIENLIM   5*CLHEP::MeV

Definition at line 22 of file NestAlg.cxx.

#define MAX_ENE   1.*CLHEP::TeV

Definition at line 21 of file NestAlg.cxx.

#define MillerDriftSpeed   true

Definition at line 3 of file NestAlg.cxx.

#define MIN_ENE   -1*CLHEP::eV

Definition at line 20 of file NestAlg.cxx.

#define phe_per_e   1

Definition at line 9 of file NestAlg.cxx.

#define PMT   0

Definition at line 19 of file NestAlg.cxx.

#define QE_EFF   1

Definition at line 8 of file NestAlg.cxx.

#define R_MAX   1*CLHEP::km

Definition at line 24 of file NestAlg.cxx.

#define R_TOL   0.2*CLHEP::mm

Definition at line 23 of file NestAlg.cxx.

#define SRF   0

Definition at line 15 of file NestAlg.cxx.

#define TOP   0

Definition at line 13 of file NestAlg.cxx.

#define WIN   0*CLHEP::mm

Definition at line 12 of file NestAlg.cxx.

Function Documentation

G4double UnivScreenFunc ( G4double  E,
G4double  Z,
G4double  A 
)

Definition at line 1253 of file NestAlg.cxx.

1253  {
1254  G4double a_0 = 5.29e-11*CLHEP::m; G4double a = 0.626*a_0*pow(Z,(-1./3.));
1255  G4double epsilon_0 = 8.854e-12*(CLHEP::farad/CLHEP::m);
1256  G4double epsilon = (a*E*2.*CLHEP::twopi*epsilon_0)/(2*pow(CLHEP::eplus,2.)*pow(Z,2.));
1257  G4double zeta_0 = pow(Z,(1./6.)); G4double m_N = A*1.66e-27*CLHEP::kg;
1258  G4double hbar = 6.582e-16*CLHEP::eV*CLHEP::s;
1259  if ( Z == 54 ) {
1260  epsilon *= 1.068; //zeta_0 = 1.63;
1261  } //special case for LXe from Bezrukov et al. 2011
1262  G4double s_n = log(1+1.1383*epsilon)/(2.*(epsilon +
1263  0.01321*pow(epsilon,0.21226) +
1264  0.19593*sqrt(epsilon)));
1265  G4double s_e = (a_0*zeta_0/a)*hbar*sqrt(8*epsilon*2.*CLHEP::twopi*epsilon_0/
1266  (a*m_N*pow(CLHEP::eplus,2.)));
1267  return 1.38e5*0.5*(1+tanh(50*epsilon-0.25))*epsilon*(s_e/s_n);
1268 }
hbar
const double hbar
Definition: CommonIMapIncludes.h:25
cet::pow
constexpr T pow(T x)
Definition: pow.h:72
genie::units::kg
static constexpr double kg
Definition: Units.h:143
genie::units::eV
static constexpr double eV
Definition: Units.h:127
a
const double a
Definition: gUpMuFluxGen.cxx:164
wirecell.units.twopi
int twopi
Definition: units.py:12
cache_state.m
m
Definition: cache_state.py:415
s
static QCString * s
Definition: config.cpp:1042

Variable Documentation

G4double biExc = 0.77

Definition at line 44 of file NestAlg.cxx.

G4bool diffusion = true

Definition at line 40 of file NestAlg.cxx.

G4bool OutElectrons =true

Definition at line 42 of file NestAlg.cxx.

G4bool SinglePhase =false

Definition at line 42 of file NestAlg.cxx.

G4bool ThomasImelTail =true

Definition at line 42 of file NestAlg.cxx.

Generated by   doxygen 1.8.11