Loading...
Searching...
No Matches
G4BraggIonModel Class Reference
#include <G4BraggIonModel.hh>
Inheritance diagram for G4BraggIonModel:
Detailed Description
Definition at line 63 of file G4BraggIonModel.hh.
Constructor & Destructor Documentation
◆ G4BraggIonModel() [1/2]
|
explicit |
Definition at line 83 of file G4BraggIonModel.cc.
85 : G4BraggModel(p, nam)
86{
87 HeMass = 3.727417*CLHEP::GeV;
88 massFactor = 1000.*CLHEP::amu_c2/HeMass;
89}
G4BraggModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="Bragg")
Definition G4BraggModel.cc:91
Referenced by G4BraggIonModel(), G4BraggNoDeltaModel::G4BraggNoDeltaModel(), and operator=().
◆ ~G4BraggIonModel()
|
override |
Definition at line 93 of file G4BraggIonModel.cc.
94{
95 if(isFirstAlpha) {
96 delete fASTAR;
97 fASTAR = nullptr;
98 }
99}
◆ G4BraggIonModel() [2/2]
|
delete |
Member Function Documentation
◆ ComputeCrossSectionPerAtom()
|
overridevirtual |
Reimplemented from G4VEmModel.
Definition at line 136 of file G4BraggIonModel.cc.
142{
143 G4double sigma =
144 Z*ComputeCrossSectionPerElectron(p,kinEnergy,cutEnergy,maxEnergy);
145 if(isAlpha) {
146 sigma *= (HeEffChargeSquare(Z, kinEnergy/CLHEP::MeV)/chargeSquareRatio);
147 }
148 return sigma;
149}
G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
Definition G4BraggModel.cc:200
◆ ComputeDEDXPerVolume()
|
overridevirtual |
Reimplemented from G4VEmModel.
Reimplemented in G4BraggNoDeltaModel.
Definition at line 172 of file G4BraggIonModel.cc.
176{
180 G4double dedx = 0.0;
181
182 if(kineticEnergy < tlim) {
183 dedx = HeDEDX(material, tlim)*std::sqrt(kineticEnergy/tlim);
184 } else {
185 dedx = HeDEDX(material, kineticEnergy);
186
187 if (tmin < tmax) {
190
191 G4double del =
192 (G4Log(x)*(tau + 1.)*(tau + 1.)/(tau * (tau + 2.0)) + 1.0 - x) *
193 CLHEP::twopi_mc2_rcl2*material->GetElectronDensity();
194 if(isAlpha) {
196 material->GetTotNbOfAtomsPerVolume();
197 heChargeSquare = HeEffChargeSquare(zeff, kineticEnergy/CLHEP::MeV);
198 del *= heChargeSquare;
199 }
200 dedx += del;
201 }
202 }
203 dedx = std::max(dedx, 0.0);
204 /*
205 G4cout << "BraggIon: " << material->GetName()
206 << " E(MeV)=" << kineticEnergy/MeV
207 << " Tmin(MeV)=" << tmin << " dedx(MeV*cm^2/g)="
208 << dedx*gram/(MeV*cm2*material->GetDensity())
209 << " q2=" << chargeSquare << G4endl;
210 */
211 return dedx;
212}
G4double MaxSecondaryEnergy(const G4ParticleDefinition *, G4double kinEnergy) final
Definition G4BraggModel.cc:368
G4double GetTotNbOfAtomsPerVolume() const
Definition G4Material.hh:197
G4double GetTotNbOfElectPerVolume() const
Definition G4Material.hh:199
Referenced by G4BraggNoDeltaModel::ComputeDEDXPerVolume().
◆ CorrectionsAlongStep()
|
overridevirtual |
Reimplemented from G4VEmModel.
Definition at line 216 of file G4BraggIonModel.cc.
222{
223 // no correction for alpha
224 if (isAlpha) { return; }
225
226 // no correction at a small step at the last step
227 if(eloss >= preKinEnergy || eloss < preKinEnergy*0.05) { return; }
228
229 // corrections only for ions
231
232 // effective energy and charge at a step
237 /*
238 G4cout << "G4BraggIonModel::CorrectionsAlongStep: Epre(MeV)="
239 << preKinEnergy << " Eeff(MeV)=" << e
240 << " eloss=" << eloss << " elossnew=" << eloss*qfactor
241 << " qfactor=" << qfactor << " Qpre=" << q20
242 << p->GetParticleName() <<G4endl;
243 */
244 eloss *= qfactor;
245}
void SetParticle(const G4ParticleDefinition *p)
Definition G4BraggModel.cc:160
◆ CrossSectionPerVolume()
|
overridevirtual |
Reimplemented from G4VEmModel.
Reimplemented in G4BraggNoDeltaModel.
Definition at line 153 of file G4BraggIonModel.cc.
159{
161 ComputeCrossSectionPerElectron(p,kinEnergy,cutEnergy,maxEnergy);
162 if(isAlpha) {
164 material->GetTotNbOfAtomsPerVolume();
165 sigma *= (HeEffChargeSquare(zeff, kinEnergy/CLHEP::MeV)/chargeSquareRatio);
166 }
167 return sigma;
168}
◆ GetChargeSquareRatio()
|
overridevirtual |
Reimplemented from G4VEmModel.
Definition at line 125 of file G4BraggIonModel.cc.
128{
129 // this method is called only for ions, so no check if it is an ion
130 if(isAlpha) { return 1.0; }
132}
G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy) override
Definition G4BraggModel.cc:171
◆ Initialise()
|
overridevirtual |
Implements G4VEmModel.
Definition at line 103 of file G4BraggIonModel.cc.
105{
106 G4BraggModel::Initialise(p, ref);
108 if(pname == "alpha") { isAlpha = true; }
109 if(isAlpha && fASTAR == nullptr) {
110 G4AutoLock l(&alphaMutex);
111 if(fASTAR == nullptr) {
112 isFirstAlpha = true;
113 fASTAR = new G4ASTARStopping();
114 }
115 l.unlock();
116 }
117 if(isFirstAlpha) {
118 fASTAR->Initialise();
119 }
120}
void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
Definition G4BraggModel.cc:118
◆ operator=()
|
delete |
The documentation for this class was generated from the following files:
