Loading...
Searching...
No Matches
G4XrayRayleighModel Class Reference
#include <G4XrayRayleighModel.hh>
Inheritance diagram for G4XrayRayleighModel:
Additional Inherited Members | |
| Protected Member Functions inherited from G4VEmModel | |
| G4ParticleChangeForLoss * | GetParticleChangeForLoss () |
| G4ParticleChangeForGamma * | GetParticleChangeForGamma () |
| virtual G4double | MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kineticEnergy) |
| const G4MaterialCutsCouple * | CurrentCouple () const |
| void | SetCurrentElement (const G4Element *) |
| Protected Attributes inherited from G4VEmModel | |
| G4ElementData * | fElementData = nullptr |
| G4VParticleChange * | pParticleChange = nullptr |
| G4PhysicsTable * | xSectionTable = nullptr |
| const G4Material * | pBaseMaterial = nullptr |
| const std::vector< G4double > * | theDensityFactor = nullptr |
| const std::vector< G4int > * | theDensityIdx = nullptr |
| G4double | inveplus |
| G4double | pFactor = 1.0 |
| std::size_t | currentCoupleIndex = 0 |
| std::size_t | basedCoupleIndex = 0 |
| G4bool | lossFlucFlag = true |
Detailed Description
Definition at line 42 of file G4XrayRayleighModel.hh.
Constructor & Destructor Documentation
◆ G4XrayRayleighModel() [1/2]
|
explicit |
Definition at line 47 of file G4XrayRayleighModel.cc.
50{
51 fParticleChange = nullptr;
52 lowEnergyLimit = 250*eV;
53 highEnergyLimit = 10.*MeV;
54 fFormFactor = 0.0;
55
56 // SetLowEnergyLimit(lowEnergyLimit);
57 SetHighEnergyLimit(highEnergyLimit);
58 //
59 verboseLevel= 0;
60 // Verbosity scale:
61 // 0 = nothing
62 // 1 = warning for energy non-conservation
63 // 2 = details of energy budget
64 // 3 = calculation of cross sections, file openings, sampling of atoms
65 // 4 = entering in methods
66
67 if(verboseLevel > 0)
68 {
70 << "Energy range: "
71 << lowEnergyLimit / eV << " eV - "
72 << highEnergyLimit / MeV << " MeV"
73 << G4endl;
74 }
75}
G4GLOB_DLL std::ostream G4cout
◆ ~G4XrayRayleighModel()
|
overridedefault |
◆ G4XrayRayleighModel() [2/2]
|
delete |
Member Function Documentation
◆ ComputeCrossSectionPerAtom()
|
overridevirtual |
Reimplemented from G4VEmModel.
Definition at line 102 of file G4XrayRayleighModel.cc.
107{
108 if (verboseLevel > 3)
109 {
111 }
112 if (gammaEnergy < lowEnergyLimit || gammaEnergy > highEnergyLimit)
113 {
114 return 0.0;
115 }
116 G4double k = gammaEnergy/hbarc;
117 k *= Bohr_radius;
118 G4double p0 = 0.680654;
119 G4double p1 = -0.0224188;
120 G4double lnZ = std::log(Z);
121
122 G4double lna = p0 + p1*lnZ;
123
125
126 G4double fo = std::pow(k, alpha);
127
128 p0 = 3.68455;
129 p1 = -0.464806;
130 lna = p0 + p1*lnZ;
131
132 fo *= 0.01*std::exp(lna);
133
134 fFormFactor = fo;
135
136 G4double b = 1. + 2.*fo;
137 G4double b2 = b*b;
138 G4double b3 = b*b2;
139
140 G4double xsc = fCofR*Z*Z/b3;
141 xsc *= fo*fo + (1. + fo)*(1. + fo);
142
143
144 return xsc;
145
146}
◆ Initialise()
|
overridevirtual |
Implements G4VEmModel.
Definition at line 83 of file G4XrayRayleighModel.cc.
85{
86 if (verboseLevel > 3)
87 {
89 }
90
91 InitialiseElementSelectors(particle,cuts);
92
93
94 if(isInitialised) return;
95 fParticleChange = GetParticleChangeForGamma();
96 isInitialised = true;
97
98}
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition G4VEmModel.cc:120
void InitialiseElementSelectors(const G4ParticleDefinition *, const G4DataVector &)
Definition G4VEmModel.cc:135
◆ operator=()
|
delete |
◆ SampleSecondaries()
|
overridevirtual |
Implements G4VEmModel.
Definition at line 150 of file G4XrayRayleighModel.cc.
155{
156 if ( verboseLevel > 3)
157 {
159 }
161
163
164
165 // Sample the angle of the scattered photon
166 // according to 1 + cosTheta*cosTheta distribution
167
168 G4double cosDipole, cosTheta, sinTheta;
169 G4double c, delta, cofA, signc = 1., a, power = 1./3.;
170
171 c = 4. - 8.*G4UniformRand();
172 a = c;
173
174 if( c < 0. )
175 {
176 signc = -1.;
177 a = -c;
178 }
179 delta = std::sqrt(a*a+4.);
180 delta += a;
181 delta *= 0.5;
182 cofA = -signc*std::pow(delta, power);
183 cosDipole = cofA - 1./cofA;
184
185 // select atom
187 photonEnergy0,aDPGamma->GetLogKineticEnergy());
189
190 G4double k = photonEnergy0/hbarc;
191 k *= Bohr_radius;
192 G4double p0 = 0.680654;
193 G4double p1 = -0.0224188;
194 G4double lnZ = std::log(Z);
195
196 G4double lna = p0 + p1*lnZ;
197
199
200 G4double fo = std::pow(k, alpha);
201
202 p0 = 3.68455;
203 p1 = -0.464806;
204 lna = p0 + p1*lnZ;
205
206 fo *= 0.01*pi*std::exp(lna);
207
208
209 G4double beta = fo/(1 + fo);
210
211 cosTheta = (cosDipole + beta)/(1. + cosDipole*beta);
212
213
214 if( cosTheta > 1.) cosTheta = 1.;
215 if( cosTheta < -1.) cosTheta = -1.;
216
217 sinTheta = std::sqrt( (1. - cosTheta)*(1. + cosTheta) );
218
219 // Scattered photon angles. ( Z - axis along the parent photon)
220
222 G4double dirX = sinTheta*std::cos(phi);
223 G4double dirY = sinTheta*std::sin(phi);
224 G4double dirZ = cosTheta;
225
226 // Update G4VParticleChange for the scattered photon
227
228 G4ThreeVector photonDirection1(dirX, dirY, dirZ);
229 photonDirection1.rotateUz(photonDirection0);
230 fParticleChange->ProposeMomentumDirection(photonDirection1);
231
232 fParticleChange->SetProposedKineticEnergy(photonEnergy0);
233}
Definition ThreeVector.h:36
const G4ThreeVector & GetMomentumDirection() const
G4double GetLogKineticEnergy() const
const G4ParticleDefinition * GetParticleDefinition() const
G4double GetKineticEnergy() const
Definition G4Element.hh:91
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition G4ParticleChangeForGamma.hh:106
void ProposeMomentumDirection(const G4ThreeVector &Pfinal)
Definition G4ParticleChangeForGamma.hh:120
const G4Element * SelectTargetAtom(const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition G4VEmModel.hh:576
The documentation for this class was generated from the following files:
