All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
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.

print PMT<< EOF< solids >< tubename="PMTVolume"rmax="(6*2.54)"z="(11.0*2.54)"deltaphi="2*(3.1415926535897)"aunit="rad"lunit="cm"/>< tubename="PMT_TPBCoating"rmax="(6.0*2.54)"z="0.01"deltaphi="2*(3.1415926535897)"aunit="rad"lunit="cm"/>< tubename="PMT_AcrylicPlate"rmax="(6.0*2.54)"z="(0.2)"deltaphi="2*(3.1415926535897)"aunit="rad"lunit="cm"/>< tubename="PMT_Stalk"rmax="(1.25*2.54)"z="(3.0*2.54)"deltaphi="2*(3.1415926535897)"aunit="rad"lunit="cm"/>< tubename="PMT_SteelBase"rmax="(6.0*2.54)"z="(1.5*2.54)"deltaphi="2*(3.1415926535897)"aunit="rad"lunit="cm"/>< tubename="PMT_Underside"rmax="2.54*4.0"z="2.54*2.5"deltaphi="2*3.1415926535897"aunit="rad"lunit="cm"/>< tubename="PMT_Lens"rmax="2.54*4.0"z="2.54*2.5"deltaphi="2*3.1415926535897"aunit="rad"lunit="cm"/></solids >< structure >< volumename="vol_PMT_TPBCoating">< materialrefref="TPB"/>< solidrefref="PMT_TPBCoating"/></volume >< volumename="vol_PMT_AcrylicPlate">< materialrefref="Acrylic"/>< solidrefref="PMT_AcrylicPlate"/></volume >< volumename="vol_PMT_Stalk">< materialrefref="Glass"/>< solidrefref="PMT_Stalk"/></volume >< volumename="vol_PMT_SteelBase">< materialrefref="STEEL_STAINLESS_Fe7Cr2Ni"/>< solidrefref="PMT_SteelBase"/></volume >< volumename="vol_PMT_Underside">< materialrefref="Glass"/>< solidrefref="PMT_Underside"/></volume >< volumename="volOpDetSensitive">< materialrefref="LAr"/>< solidrefref="PMT_Lens"/></volume >< volumename="volPMT">< materialrefref="LAr"/>< solidrefref="PMTVolume"/>< physvol >< volumerefref="vol_PMT_TPBCoating"/>< positionname="pos_PMT_TPBCoating"unit="cm"x="0"y="0"z="(5.5*2.54)-(0.5*0.005)"/></physvol >< physvol >< volumerefref="vol_PMT_AcrylicPlate"/>< positionname="pos_PMT_AcrylicPlate"unit="cm"x="0"y="0"z="(5.5*2.54)-0.01-(0.5*0.2)"/></physvol >< physvol >< volumerefref="vol_PMT_Stalk"/>< positionname="pos_PMT_Stalk"unit="cm"x="0"y="0"z="(3.0*2.54)-(5.5*2.54)"/></physvol >< physvol >< volumerefref="vol_PMT_SteelBase"/>< positionname="pos_PMT_SteelBase"unit="cm"x="0"y="0"z="(0.75*2.54)-(5.5*2.54)"/></physvol >< physvol >< volumerefref="volOpDetSensitive"/>< positionname="pos_PMT_Lens"unit="cm"x="0"y="0"z="(7.0*2.54)-(5.5*2.54)"/></physvol >< physvol >< volumerefref="vol_PMT_Underside"/>< positionname="pos_PMT_Underside"unit="cm"x="0"y="0"z="(7.0*2.54)-(5.5*2.54)"/></physvol ></volume ></structure > EOF sub gen_cryostat()   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 }
E
process_name E
Definition: rawtorootconvert.fcl:26
gen_crt_frags.m
tuple m
now if test mode generate materials, CRT shell, world, gdml header else just generate CRT shell for u...
Definition: gen_crt_frags.py:1347
a
process_name gaushit a
Definition: decode_signalprocess_icarus.fcl:37
s
then echo File list $list not found else cat $list while read file do echo $file sed s
Definition: file_to_url.sh:60
lib.dedx.A
float A
Definition: dedx.py:137

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