#include "Geant4/G4Ions.hh"
#include "Geant4/G4OpticalPhoton.hh"
#include "Geant4/G4VProcess.hh"
#include "larsim/LegacyLArG4/G4ThermalElectron.hh"
#include "larsim/LegacyLArG4/NestAlg.h"
#include "CLHEP/Random/RandFlat.h"
#include "CLHEP/Random/RandGauss.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 larg4::ISCalcNESTLAr::CalcIonAndScint(), and 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 1284 of file NestAlg.cxx.

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);
 }