80 lowEnergyLimit(1.022000*
MeV),
81 highEnergyLimit(100*
GeV),
82 intrinsicLowEnergyLimit(1.022000*
MeV),
83 intrinsicHighEnergyLimit(100*
GeV),
90 G4Exception(
"G4LowEnergyGammaConversion::G4LowEnergyGammaConversion()",
92 "Energy limit outside intrinsic process validity range!");
123 G4String crossSectionFile =
"pair/pp-cs-";
166 G4cout <<
"G4LowEnergyGammaConversion::PostStepDoIt - element = 0" <<
G4endl;
171 G4cout <<
"G4LowEnergyGammaConversion::PostStepDoIt - ionisation = 0" <<
G4endl;
176 if (photonEnergy > 50. *
MeV) fZ += 8. * (element->
GetfCoulomb());
180 G4double screenMax = std::exp ((42.24 - fZ)/8.368) - 0.952 ;
184 G4double epsilon1 = 0.5 - 0.5 * std::sqrt(1. - screenMin / screenMax) ;
186 G4double epsilonRange = 0.5 - epsilonMin ;
200 epsilon = 0.5 - epsilonRange * std::pow(
G4UniformRand(), 0.3333) ;
201 screen = screenFactor / (epsilon * (1. -
epsilon));
207 screen = screenFactor / (epsilon * (1 -
epsilon));
221 electronTotEnergy = (1. -
epsilon) * photonEnergy;
222 positronTotEnergy = epsilon * photonEnergy;
226 positronTotEnergy = (1. -
epsilon) * photonEnergy;
227 electronTotEnergy = epsilon * photonEnergy;
253 G4double dxEle= std::sin(thetaEle)*std::cos(phi),dyEle= std::sin(thetaEle)*std::sin(phi),dzEle=std::cos(thetaEle);
254 G4double dxPos=-std::sin(thetaPos)*std::cos(phi),dyPos=-std::sin(thetaPos)*std::sin(phi),dzPos=std::cos(thetaPos);
273 electronDirection.
rotateUz(photonDirection);
282 localEnergyDeposit += electronKineEnergy ;
291 localEnergyDeposit += positronKineEnergy ;
292 positronKineEnergy = 0. ;
296 positronDirection.
rotateUz(photonDirection);
300 positronDirection, positronKineEnergy);
326 size_t materialIndex = couple->
GetIndex();
341 if (screenVariable > 1.)
342 value = 42.24 - 8.368 * std::log(screenVariable + 0.952);
344 value = 42.392 - screenVariable * (7.796 - 1.961 * screenVariable);
355 if (screenVariable > 1.)
356 value = 42.24 - 8.368 * std::log(screenVariable + 0.952);
358 value = 41.405 - screenVariable * (5.828 - 0.8945 * screenVariable);
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4RDVEMDataSet * BuildMeanFreePathForMaterials(const G4DataVector *energyCuts=0)
const G4double intrinsicHighEnergyLimit
static constexpr double MeV
G4LowEnergyGammaConversion(const G4String &processName="LowEnConversion")
G4RDVRangeTest * rangeTest
const G4ThreeVector & GetMomentumDirection() const
virtual G4bool Escape(const G4ParticleDefinition *particle, const G4MaterialCutsCouple *couple, G4double energy, G4double safety) const =0
G4double GetMeanFreePath(const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
Hep3Vector & rotateUz(const Hep3Vector &)
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
static constexpr double epsilon0
const G4double smallEnergy
void AddSecondary(G4Track *aSecondary)
G4ParticleChange aParticleChange
const G4double intrinsicLowEnergyLimit
void LoadData(const G4String &dataFile)
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
const XML_Char int const XML_Char * value
virtual void Initialize(const G4Track &)
virtual G4double FindValue(G4double x, G4int componentId=0) const =0
static constexpr double electron_mass_c2
G4RDVCrossSectionHandler * crossSectionHandler
static constexpr double twopi
const G4String & GetProcessName() const
const G4Element * SelectRandomElement(const G4MaterialCutsCouple *material, G4double e) const
static G4Positron * Positron()
G4StepPoint * GetPostStepPoint() const
void ProposeEnergy(G4double finalEnergy)
G4VParticleChange * PostStepDoIt(const G4Track &aTrack, const G4Step &aStep)
static G4Electron * Electron()
G4bool IsApplicable(const G4ParticleDefinition &photon)
double epsilon(double density, double temperature)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
void BuildPhysicsTable(const G4ParticleDefinition &photon)
~G4LowEnergyGammaConversion()
G4double ScreenFunction2(G4double screenVariable)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4double GetSafety() const
G4IonisParamElm * GetIonisation() const
void Initialise(G4RDVDataSetAlgorithm *interpolation=0, G4double minE=250 *CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int numberOfBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
G4RDVEMDataSet * meanFreePathTable
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4double GetfCoulomb() const
G4double GetlogZ3() const
G4double ScreenFunction1(G4double screenVariable)
void ProposeTrackStatus(G4TrackStatus status)
static constexpr double GeV
const G4DynamicParticle * GetDynamicParticle() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
void SetNumberOfSecondaries(G4int totSecondaries)