75 lowEnergyLimit(250*
eV),
76 highEnergyLimit(100*
GeV),
77 intrinsicLowEnergyLimit(10*
eV),
78 intrinsicHighEnergyLimit(100*
GeV)
83 G4Exception(
"G4LowEnergyCompton::G4LowEnergyCompton()",
85 "Energy outside intrinsic process validity range!");
91 G4String scatterFile =
"comp/ce-sf-";
124 G4String crossSectionFile =
"comp/ce-cs-";
170 G4double alpha1 = -std::log(epsilon0);
171 G4double alpha2 = 0.5 * (1. - epsilon0Sq);
193 epsilonSq = epsilon0Sq + (1. - epsilon0Sq) *
G4UniformRand();
194 epsilon = std::sqrt(epsilonSq);
197 oneCosT = (1. -
epsilon) / ( epsilon * e0m);
198 sinT2 = oneCosT * (2. - oneCosT);
199 G4double x = std::sqrt(oneCosT/2.) / (wlPhoton/
cm);
201 gReject = (1. - epsilon * sinT2 / (1. + epsilonSq)) * scatteringFunction;
206 G4double sinTheta = std::sqrt (sinT2);
208 G4double dirX = sinTheta * std::cos(phi);
209 G4double dirY = sinTheta * std::sin(phi);
218 G4int maxDopplerIterations = 1000;
220 G4double photonEoriginal = epsilon * photonEnergy0;
231 eMax = photonEnergy0 - bindingE;
237 G4double pDoppler2 = pDoppler * pDoppler;
239 G4double var3 = var2*var2 - pDoppler2;
240 G4double var4 = var2 - pDoppler2 * cosTheta;
241 G4double var = var4*var4 - var3 + pDoppler2 * var3;
247 if (
G4UniformRand() < 0.5) photonE = (var4 - varSqrt) * scale;
248 else photonE = (var4 + varSqrt) * scale;
254 }
while ( iteration <= maxDopplerIterations &&
255 (photonE < 0. || photonE > eMax || photonE < eMax*
G4UniformRand()) );
259 if (iteration >= maxDopplerIterations)
261 photonE = photonEoriginal;
268 photonDirection1.
rotateUz(photonDirection0);
274 if (photonEnergy1 > 0.)
285 G4double eKineticEnergy = photonEnergy0 - photonEnergy1 - bindingE;
288 G4double electronE = photonEnergy0 * (1. -
epsilon) + electron_mass_c2;
294 cosThetaE = (eTotalEnergy + photonEnergy1 )* (1. - epsilon) / std::sqrt(electronP2);
295 sinThetaE = -1. * std::sqrt(1. - cosThetaE * cosThetaE);
298 G4double eDirX = sinThetaE * std::cos(phi);
299 G4double eDirY = sinThetaE * std::sin(phi);
309 eDirection.
rotateUz(photonDirection0);
338 size_t materialIndex = couple->
GetIndex();
void BuildPhysicsTable(const G4ParticleDefinition &definition)
G4RDVEMDataSet * BuildMeanFreePathForMaterials(const G4DataVector *energyCuts=0)
G4int SelectRandomAtom(const G4MaterialCutsCouple *couple, G4double e) const
static constexpr double keV
const G4ThreeVector & GetMomentumDirection() const
virtual G4bool Escape(const G4ParticleDefinition *particle, const G4MaterialCutsCouple *couple, G4double energy, G4double safety) const =0
Hep3Vector & rotateUz(const Hep3Vector &)
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
static constexpr double epsilon0
static constexpr double h_Planck
void AddSecondary(G4Track *aSecondary)
G4ParticleChange aParticleChange
void LoadData(const G4String &dataFile)
G4VParticleChange * PostStepDoIt(const G4Track &track, const G4Step &step)
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
virtual void Initialize(const G4Track &)
G4RDDopplerProfile profileData
virtual G4double FindValue(G4double x, G4int componentId=0) const =0
const G4double intrinsicLowEnergyLimit
static constexpr double electron_mass_c2
const G4double intrinsicHighEnergyLimit
static constexpr double twopi
virtual G4bool LoadData(const G4String &fileName)=0
G4double RandomSelectMomentum(G4int Z, G4int shellIndex) const
const G4String & GetProcessName() const
static constexpr double eV
G4StepPoint * GetPostStepPoint() const
void ProposeEnergy(G4double finalEnergy)
static G4Electron * Electron()
double epsilon(double density, double temperature)
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
static constexpr double c_light
G4bool IsApplicable(const G4ParticleDefinition &definition)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4RDVEMDataSet * meanFreePathTable
G4double GetSafety() const
G4double GetKineticEnergy() const
static constexpr double cm
G4GLOB_DLL std::ostream G4cout
G4RDVEMDataSet * scatterFunctionData
G4RDVCrossSectionHandler * crossSectionHandler
G4RDVRangeTest * rangeTest
G4double BindingEnergy(G4int Z, G4int shellIndex) const
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4double GetMeanFreePath(const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
static constexpr double fine_structure_const
void LoadData(const G4String &fileName)
G4int SelectRandomShell(G4int Z) const
void ProposeTrackStatus(G4TrackStatus status)
static constexpr double GeV
const G4DynamicParticle * GetDynamicParticle() const
G4LowEnergyCompton(const G4String &processName="LowEnCompton")
void SetNumberOfSecondaries(G4int totSecondaries)