31 #ifndef G4LivermorePolarizedPhotoElectricGDModel_h
32 #define G4LivermorePolarizedPhotoElectricGDModel_h 1
49 const G4String& nam =
"LivermorePolarizedPhotoElectric");
96 static std::vector<G4double>*
fParam[99];
std::vector< G4double > fSandiaCof
static G4int fNShells[99]
std::vector< ExP01TrackerHit * > a
G4double SetPhi(G4double, G4double, G4double)
void ReadData(G4int Z, const char *path=0)
G4ThreeVector GetRandomPolarization(G4ThreeVector &direction0)
static const G4double emax
G4LivermorePolarizedPhotoElectricGDModel & operator=(const G4LivermorePolarizedPhotoElectricGDModel &right)
void SetLimitNumberOfShells(G4int)
static G4LPhysicsFreeVector * fCrossSectionLE[99]
void SystemOfRefChange(G4ThreeVector &direction0, G4ThreeVector &direction1, G4ThreeVector &polarization0)
static G4Material * fWater
G4ParticleDefinition * theGamma
G4ThreeVector SetPerpendicularVector(G4ThreeVector &a)
static std::vector< G4double > * fParam[99]
double A(double temperature)
G4ParticleChangeForGamma * fParticleChange
G4double SetCosTheta(G4double)
G4LivermorePolarizedPhotoElectricGDModel(const G4String &nam="LivermorePolarizedPhotoElectric")
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0, G4double cut=0, G4double emax=DBL_MAX)
static G4ElementData * fShellCrossSection
double epsilon(double density, double temperature)
static G4int fNShellsUsed[99]
G4ParticleDefinition * theElectron
G4ThreeVector GetPerpendicularPolarization(const G4ThreeVector &direction0, const G4ThreeVector &polarization0) const
G4VAtomDeexcitation * fAtomDeexcitation
G4ThreeVector SetNewPolarization(G4double epsilon, G4double sinSqrTheta, G4double phi, G4double cosTheta)
G4bool fDeexcitationActive
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
static G4double fWaterEnergyLimit
virtual ~G4LivermorePolarizedPhotoElectricGDModel()
virtual void InitialiseForElement(const G4ParticleDefinition *, G4int Z)
static G4LPhysicsFreeVector * fCrossSection[99]
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)