91 G4cout <<
"Calling G4XrayRayleighModel::Initialise()" <<
G4endl;
113 G4cout <<
"Calling CrossSectionPerAtom() of G4XrayRayleighModel" <<
G4endl;
135 fo *= 0.01*std::exp(lna);
144 xsc *= fo*fo + (1. + fo)*(1. + fo);
161 G4cout <<
"Calling SampleSecondaries() of G4XrayRayleighModel" <<
G4endl;
171 G4double cosDipole, cosTheta, sinTheta;
172 G4double c, delta, cofA, signc = 1., a, power = 1./3.;
182 delta = std::sqrt(a*a+4.);
185 cofA = -signc*std::pow(delta, power);
186 cosDipole = cofA - 1./cofA;
209 fo *= 0.01*
pi*std::exp(lna);
214 cosTheta = (cosDipole +
beta)/(1. + cosDipole*
beta);
217 if( cosTheta > 1.) cosTheta = 1.;
218 if( cosTheta < -1.) cosTheta = -1.;
220 sinTheta = std::sqrt( (1. - cosTheta)*(1. + cosTheta) );
225 G4double dirX = sinTheta*std::cos(phi);
226 G4double dirY = sinTheta*std::sin(phi);
232 photonDirection1.
rotateUz(photonDirection0);
static const G4double alpha
static constexpr double pi2
static constexpr double twopi
static constexpr double eV
static constexpr double MeV
static constexpr double pi
G4GLOB_DLL std::ostream G4cout
Hep3Vector & rotateUz(const Hep3Vector &)
const G4ThreeVector & GetMomentumDirection() const
G4double GetLogKineticEnergy() const
const G4ParticleDefinition * GetParticleDefinition() const
G4double GetKineticEnergy() const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
void SetHighEnergyLimit(G4double)
G4ParticleChangeForGamma * GetParticleChangeForGamma()
const G4Element * SelectTargetAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
static const G4double fCofR
~G4XrayRayleighModel() override
G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0, G4double cut=0, G4double emax=DBL_MAX) override
G4XrayRayleighModel(const G4ParticleDefinition *p=0, const G4String &nam="XrayRayleigh")
void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static const G4double fCofA
G4ParticleChangeForGamma * fParticleChange
int classic_electr_radius