#include "Geant4/Randomize.hh"
#include "Geant4/G4Ions.hh"
#include "Geant4/G4OpticalPhoton.hh"
#include "Geant4/G4VProcess.hh"
#include "larsim/LArG4/NestAlg.h"
#include "larsim/LArG4/G4ThermalElectron.hh"
#include "CLHEP/Random/RandGauss.h"
#include "CLHEP/Random/RandFlat.h"

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.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define AVO 6.022e23

Definition at line 1 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define BORDER 0*CLHEP::cm

Definition at line 6 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define BOT 0

Definition at line 18 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define CTH 0

Definition at line 17 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define Density_LAr 1.393

Definition at line 26 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#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.

Referenced by NestAlg::CalculateIonizationAndScintillation().

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

Definition at line 2 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation(), and NestAlg::GetLiquidElectronDriftSpeed().

#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.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define HIENLIM 5*CLHEP::MeV

Definition at line 22 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#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.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define MIN_ENE -1*CLHEP::eV

Definition at line 20 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define phe_per_e 1

Definition at line 9 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define PMT 0

Definition at line 19 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define QE_EFF 1

Definition at line 8 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#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.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define SRF 0

Definition at line 15 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define TOP 0

Definition at line 13 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

#define WIN 0*CLHEP::mm

Definition at line 12 of file NestAlg.cxx.

Referenced by NestAlg::CalculateIonizationAndScintillation().

Function Documentation

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

Definition at line 1254 of file NestAlg.cxx.

References s.

Referenced by NestAlg::EnergyDeposition().

                                                                {
   G4double a_0 = 5.29e-11*CLHEP::m; G4double a = 0.626*a_0*pow(Z,(-1./3.));
   G4double epsilon_0 = 8.854e-12*(CLHEP::farad/CLHEP::m);
   G4double epsilon = (a*E*2.*CLHEP::twopi*epsilon_0)/(2*pow(CLHEP::eplus,2.)*pow(Z,2.));
   G4double zeta_0 = pow(Z,(1./6.)); G4double m_N = A*1.66e-27*CLHEP::kg;
   G4double hbar = 6.582e-16*CLHEP::eV*CLHEP::s;
   if ( Z == 54 ) {
     epsilon *= 1.068; //zeta_0 = 1.63;
   } //special case for LXe from Bezrukov et al. 2011
   G4double s_n = log(1+1.1383*epsilon)/(2.*(epsilon +
                                             0.01321*pow(epsilon,0.21226) +
                                             0.19593*sqrt(epsilon)));
   G4double s_e = (a_0*zeta_0/a)*hbar*sqrt(8*epsilon*2.*CLHEP::twopi*epsilon_0/
                                           (a*m_N*pow(CLHEP::eplus,2.)));
   return 1.38e5*0.5*(1+tanh(50*epsilon-0.25))*epsilon*(s_e/s_n);
 }