Geant4 11.1.1
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4AdjointIonIonisationModel Class Reference

#include <G4AdjointIonIonisationModel.hh>

+ Inheritance diagram for G4AdjointIonIonisationModel:

Public Member Functions

 G4AdjointIonIonisationModel ()
 
 ~G4AdjointIonIonisationModel () override
 
void SampleSecondaries (const G4Track &aTrack, G4bool isScatProjToProj, G4ParticleChange *fParticleChange) override
 
G4double DiffCrossSectionPerAtomPrimToSecond (G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.) override
 
void CorrectPostStepWeight (G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool isScatProjToProj) override
 
G4double GetSecondAdjEnergyMaxForScatProjToProj (G4double primAdjEnergy) override
 
G4double GetSecondAdjEnergyMinForScatProjToProj (G4double primAdjEnergy, G4double tcut=0.) override
 
G4double GetSecondAdjEnergyMaxForProdToProj (G4double primAdjEnergy) override
 
G4double GetSecondAdjEnergyMinForProdToProj (G4double primAdjEnergy) override
 
void SetUseOnlyBragg (G4bool aBool)
 
void SetIon (G4ParticleDefinition *adj_ion, G4ParticleDefinition *fwd_ion)
 
 G4AdjointIonIonisationModel (G4AdjointIonIonisationModel &)=delete
 
G4AdjointIonIonisationModel & operator= (const G4AdjointIonIonisationModel &right)=delete
 
- Public Member Functions inherited from G4VEmAdjointModel
 G4VEmAdjointModel (const G4String &nam)
 
virtual ~G4VEmAdjointModel ()
 
virtual void SampleSecondaries (const G4Track &aTrack, G4bool isScatProjToProj, G4ParticleChange *fParticleChange)=0
 
virtual G4double AdjointCrossSection (const G4MaterialCutsCouple *aCouple, G4double primEnergy, G4bool isScatProjToProj)
 
virtual G4double DiffCrossSectionPerAtomPrimToSecond (G4double kinEnergyProj, G4double kinEnergyProd, G4double Z, G4double A=0.)
 
virtual G4double DiffCrossSectionPerAtomPrimToScatPrim (G4double kinEnergyProj, G4double kinEnergyScatProj, G4double Z, G4double A=0.)
 
virtual G4double DiffCrossSectionPerVolumePrimToSecond (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyProd)
 
virtual G4double DiffCrossSectionPerVolumePrimToScatPrim (const G4Material *aMaterial, G4double kinEnergyProj, G4double kinEnergyScatProj)
 
virtual G4double GetSecondAdjEnergyMaxForScatProjToProj (G4double primAdjEnergy)
 
virtual G4double GetSecondAdjEnergyMinForScatProjToProj (G4double primAdjEnergy, G4double tcut=0.)
 
virtual G4double GetSecondAdjEnergyMaxForProdToProj (G4double primAdjEnergy)
 
virtual G4double GetSecondAdjEnergyMinForProdToProj (G4double primAdjEnergy)
 
void DefineCurrentMaterial (const G4MaterialCutsCouple *couple)
 
std::vector< std::vector< double > * > ComputeAdjointCrossSectionVectorPerAtomForSecond (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)
 
std::vector< std::vector< double > * > ComputeAdjointCrossSectionVectorPerAtomForScatProj (G4double kinEnergyProd, G4double Z, G4double A=0., G4int nbin_pro_decade=10)
 
std::vector< std::vector< double > * > ComputeAdjointCrossSectionVectorPerVolumeForSecond (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)
 
std::vector< std::vector< double > * > ComputeAdjointCrossSectionVectorPerVolumeForScatProj (G4Material *aMaterial, G4double kinEnergyProd, G4int nbin_pro_decade=10)
 
void SetCSMatrices (std::vector< G4AdjointCSMatrix * > *Vec1CSMatrix, std::vector< G4AdjointCSMatrix * > *Vec2CSMatrix)
 
G4ParticleDefinition * GetAdjointEquivalentOfDirectPrimaryParticleDefinition ()
 
G4ParticleDefinition * GetAdjointEquivalentOfDirectSecondaryParticleDefinition ()
 
G4double GetHighEnergyLimit ()
 
G4double GetLowEnergyLimit ()
 
void SetHighEnergyLimit (G4double aVal)
 
void SetLowEnergyLimit (G4double aVal)
 
void DefineDirectEMModel (G4VEmModel *aModel)
 
void SetAdjointEquivalentOfDirectPrimaryParticleDefinition (G4ParticleDefinition *aPart)
 
void SetAdjointEquivalentOfDirectSecondaryParticleDefinition (G4ParticleDefinition *aPart)
 
void SetSecondPartOfSameType (G4bool aBool)
 
G4bool GetSecondPartOfSameType ()
 
void SetUseMatrix (G4bool aBool)
 
void SetUseMatrixPerElement (G4bool aBool)
 
void SetUseOnlyOneMatrixForAllElements (G4bool aBool)
 
void SetApplyCutInRange (G4bool aBool)
 
G4bool GetUseMatrix ()
 
G4bool GetUseMatrixPerElement ()
 
G4bool GetUseOnlyOneMatrixForAllElements ()
 
G4bool GetApplyCutInRange ()
 
G4String GetName ()
 
virtual void SetCSBiasingFactor (G4double aVal)
 
void SetCorrectWeightForPostStepInModel (G4bool aBool)
 
void SetAdditionalWeightCorrectionFactorForPostStepOutsideModel (G4double factor)
 
 G4VEmAdjointModel (G4VEmAdjointModel &)=delete
 
G4VEmAdjointModel & operator= (const G4VEmAdjointModel &right)=delete
 

Additional Inherited Members

- Protected Member Functions inherited from G4VEmAdjointModel
G4double DiffCrossSectionFunction1 (G4double kinEnergyProj)
 
G4double DiffCrossSectionFunction2 (G4double kinEnergyProj)
 
G4double SampleAdjSecEnergyFromCSMatrix (std::size_t MatrixIndex, G4double prim_energy, G4bool isScatProjToProj)
 
G4double SampleAdjSecEnergyFromCSMatrix (G4double prim_energy, G4bool isScatProjToProj)
 
void SelectCSMatrix (G4bool isScatProjToProj)
 
virtual G4double SampleAdjSecEnergyFromDiffCrossSectionPerAtom (G4double prim_energy, G4bool isScatProjToProj)
 
virtual void CorrectPostStepWeight (G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool isScatProjToProj)
 
- Protected Attributes inherited from G4VEmAdjointModel
G4AdjointCSManager * fCSManager
 
G4VEmModel * fDirectModel = nullptr
 
const G4String fName
 
G4Material * fSelectedMaterial = nullptr
 
G4Material * fCurrentMaterial = nullptr
 
G4MaterialCutsCouple * fCurrentCouple = nullptr
 
G4ParticleDefinition * fAdjEquivDirectPrimPart = nullptr
 
G4ParticleDefinition * fAdjEquivDirectSecondPart = nullptr
 
G4ParticleDefinition * fDirectPrimaryPart = nullptr
 
std::vector< G4AdjointCSMatrix * > * fCSMatrixProdToProjBackScat = nullptr
 
std::vector< G4AdjointCSMatrix * > * fCSMatrixProjToProjBackScat = nullptr
 
std::vector< G4double > fElementCSScatProjToProj
 
std::vector< G4double > fElementCSProdToProj
 
G4double fKinEnergyProdForIntegration = 0.
 
G4double fKinEnergyScatProjForIntegration = 0.
 
G4double fLastCS = 0.
 
G4double fLastAdjointCSForScatProjToProj = 0.
 
G4double fLastAdjointCSForProdToProj = 0.
 
G4double fPreStepEnergy = 0.
 
G4double fTcutPrim = 0.
 
G4double fTcutSecond = 0.
 
G4double fHighEnergyLimit = 0.
 
G4double fLowEnergyLimit = 0.
 
G4double fCsBiasingFactor = 1.
 
G4double fOutsideWeightFactor = 1.
 
G4int fASelectedNucleus = 0
 
G4int fZSelectedNucleus = 0
 
std::size_t fCSMatrixUsed = 0
 
G4bool fSecondPartSameType = false
 
G4bool fInModelWeightCorr
 
G4bool fApplyCutInRange = true
 
G4bool fUseMatrix = false
 
G4bool fUseMatrixPerElement = false
 
G4bool fOneMatrixForAllElements = false
 

Detailed Description

Definition at line 46 of file G4AdjointIonIonisationModel.hh.

Constructor & Destructor Documentation

◆ G4AdjointIonIonisationModel() [1/2]

G4AdjointIonIonisationModel::G4AdjointIonIonisationModel ( )

Definition at line 41 of file G4AdjointIonIonisationModel.cc.

42 : G4VEmAdjointModel("Adjoint_IonIonisation")
43{
44 fUseMatrix = true;
45 fUseMatrixPerElement = true;
46 fApplyCutInRange = true;
47 fOneMatrixForAllElements = true;
48 fSecondPartSameType = false;
49
50 // The direct EM Model is taken as BetheBloch. It is only used for the
51 // computation of the differential cross section.
52 // The Bragg model could be used as an alternative as it offers the same
53 // differential cross section
54
55 fBetheBlochDirectEMModel = new G4BetheBlochModel(G4GenericIon::GenericIon());
56 fBraggIonDirectEMModel = new G4BraggIonModel(G4GenericIon::GenericIon());
57 fAdjEquivDirectSecondPart = G4AdjointElectron::AdjointElectron();
58 fDirectPrimaryPart = nullptr;
59}
G4AdjointElectron::AdjointElectron
static G4AdjointElectron * AdjointElectron()
Definition: G4AdjointElectron.cc:99
G4GenericIon::GenericIon
static G4GenericIon * GenericIon()
Definition: G4GenericIon.cc:92
G4VEmAdjointModel::fAdjEquivDirectSecondPart
G4ParticleDefinition * fAdjEquivDirectSecondPart
Definition: G4VEmAdjointModel.hh:260
G4VEmAdjointModel::fDirectPrimaryPart
G4ParticleDefinition * fDirectPrimaryPart
Definition: G4VEmAdjointModel.hh:261

◆ ~G4AdjointIonIonisationModel()

G4AdjointIonIonisationModel::~G4AdjointIonIonisationModel ( )
override

Definition at line 62 of file G4AdjointIonIonisationModel.cc.

62{}

◆ G4AdjointIonIonisationModel() [2/2]

G4AdjointIonIonisationModel::G4AdjointIonIonisationModel ( G4AdjointIonIonisationModel &  )
delete

Member Function Documentation

◆ CorrectPostStepWeight()

void G4AdjointIonIonisationModel::CorrectPostStepWeight ( G4ParticleChange *  fParticleChange,
G4double  old_weight,
G4double  adjointPrimKinEnergy,
G4double  projectileKinEnergy,
G4bool  isScatProjToProj 
)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 228 of file G4AdjointIonIonisationModel.cc.

231{
232 // It is needed because the direct cross section used to compute the
233 // differential cross section is not the one used in
234 // the direct model where the GenericIon stuff is considered with correction
235 // of effective charge. In the direct model the samnepl of secondaries does
236 // not reflect the integral cross section. The integral fwd cross section that
237 // we used to compute the differential CS match the sample of secondaries in
238 // the forward case despite the fact that its is not the same total CS than in
239 // the FWD case. For this reason an extra weight correction is needed at the
240 // end.
241
242 G4double new_weight = old_weight;
243
244 // the correction of CS due to the problem explained above
245 G4double kinEnergyProjScaled = fMassRatio * projectileKinEnergy;
246 fDirectModel = fBraggIonDirectEMModel;
247 if(kinEnergyProjScaled > 2. * MeV && !fUseOnlyBragg)
248 fDirectModel = fBetheBlochDirectEMModel;
249 G4double UsedFwdCS = fDirectModel->ComputeCrossSectionPerAtom(
250 fDirectPrimaryPart, projectileKinEnergy, 1, 1, fTcutSecond, 1.e20);
251 G4double chargeSqRatio = 1.;
252 if(fChargeSquare > 1.)
253 chargeSqRatio = fDirectModel->GetChargeSquareRatio(
254 fDirectPrimaryPart, fCurrentMaterial, projectileKinEnergy);
255 G4double CorrectFwdCS =
256 chargeSqRatio * fDirectModel->ComputeCrossSectionPerAtom(
257 G4GenericIon::GenericIon(), kinEnergyProjScaled, 1, 1,
258 fTcutSecond, 1.e20);
259 // May be some check is needed if UsedFwdCS ==0 probably that then we should
260 // avoid a secondary to be produced,
261 if(UsedFwdCS > 0.)
262 new_weight *= CorrectFwdCS / UsedFwdCS;
263
264 // additional CS correction needed for cross section biasing in general.
265 // May be wrong for ions. Most of the time not used.
266 new_weight *=
267 G4AdjointCSManager::GetAdjointCSManager()->GetPostStepWeightCorrection() /
268 fCsBiasingFactor;
269
270 new_weight *= projectileKinEnergy / adjointPrimKinEnergy;
271
272 fParticleChange->SetParentWeightByProcess(false);
273 fParticleChange->SetSecondaryWeightByProcess(false);
274 fParticleChange->ProposeParentWeight(new_weight);
275}
G4double
double G4double
Definition: G4Types.hh:83
G4AdjointCSManager::GetPostStepWeightCorrection
G4double GetPostStepWeightCorrection()
Definition: G4AdjointCSManager.cc:578
G4AdjointCSManager::GetAdjointCSManager
static G4AdjointCSManager * GetAdjointCSManager()
Definition: G4AdjointCSManager.cc:57
G4VEmAdjointModel::fCurrentMaterial
G4Material * fCurrentMaterial
Definition: G4VEmAdjointModel.hh:255
G4VEmAdjointModel::fDirectModel
G4VEmModel * fDirectModel
Definition: G4VEmAdjointModel.hh:250
G4VEmModel::ComputeCrossSectionPerAtom
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
Definition: G4VEmModel.cc:284
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:325
G4VParticleChange::SetSecondaryWeightByProcess
void SetSecondaryWeightByProcess(G4bool)
G4VParticleChange::SetParentWeightByProcess
void SetParentWeightByProcess(G4bool)
Definition: G4VParticleChange.cc:375
G4VParticleChange::ProposeParentWeight
void ProposeParentWeight(G4double finalWeight)

Referenced by SampleSecondaries().

◆ DiffCrossSectionPerAtomPrimToSecond()

G4double G4AdjointIonIonisationModel::DiffCrossSectionPerAtomPrimToSecond ( G4double  kinEnergyProj,
G4double  kinEnergyProd,
G4double  Z,
G4double  A = 0. 
)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 132 of file G4AdjointIonIonisationModel.cc.

134{
135 // Probably that here the Bragg Model should be also used for
136 // kinEnergyProj/nuc<2MeV
137 G4double dSigmadEprod = 0.;
138 G4double Emax_proj = GetSecondAdjEnergyMaxForProdToProj(kinEnergyProd);
139 G4double Emin_proj = GetSecondAdjEnergyMinForProdToProj(kinEnergyProd);
140
141 G4double kinEnergyProjScaled = fMassRatio * kinEnergyProj;
142
143 // the produced particle should have a kinetic energy smaller than the
144 // projectile
145 if(kinEnergyProj > Emin_proj && kinEnergyProj <= Emax_proj)
146 {
147 G4double Tmax = kinEnergyProj;
148
149 G4double E1 = kinEnergyProd;
150 G4double E2 = kinEnergyProd * 1.000001;
151 G4double dE = (E2 - E1);
152 G4double sigma1, sigma2;
153 fDirectModel = fBraggIonDirectEMModel;
154 if(kinEnergyProjScaled > 2. * MeV && !fUseOnlyBragg)
155 fDirectModel = fBetheBlochDirectEMModel;
156 sigma1 = fDirectModel->ComputeCrossSectionPerAtom(
157 fDirectPrimaryPart, kinEnergyProj, Z, A, E1, 1.e20);
158 sigma2 = fDirectModel->ComputeCrossSectionPerAtom(
159 fDirectPrimaryPart, kinEnergyProj, Z, A, E2, 1.e20);
160
161 dSigmadEprod = (sigma1 - sigma2) / dE;
162
163 if(dSigmadEprod > 1.)
164 {
165 G4cout << "sigma1 " << kinEnergyProj / MeV << '\t' << kinEnergyProd / MeV
166 << '\t' << sigma1 << G4endl;
167 G4cout << "sigma2 " << kinEnergyProj / MeV << '\t' << kinEnergyProd / MeV
168 << '\t' << sigma2 << G4endl;
169 G4cout << "dsigma " << kinEnergyProj / MeV << '\t' << kinEnergyProd / MeV
170 << '\t' << dSigmadEprod << G4endl;
171 }
172
173 if(fDirectModel == fBetheBlochDirectEMModel)
174 {
175 // correction of differential cross section at high energy to correct for
176 // the suppression of particle at secondary at high energy used in the
177 // Bethe Bloch Model. This correction consist to multiply by g the
178 // probability function used to test the rejection of a secondary Source
179 // code taken from G4BetheBlochModel::SampleSecondaries
180
181 G4double deltaKinEnergy = kinEnergyProd;
182
183 G4double x = fFormFact * deltaKinEnergy;
184 if(x > 1.e-6)
185 {
186 G4double totEnergy = kinEnergyProj + fMass;
187 G4double etot2 = totEnergy * totEnergy;
188 G4double beta2 = kinEnergyProj * (kinEnergyProj + 2.0 * fMass) / etot2;
189 G4double f1 = 0.0;
190 G4double f = 1.0 - beta2 * deltaKinEnergy / Tmax;
191 if(0.5 == fSpin)
192 {
193 f1 = 0.5 * deltaKinEnergy * deltaKinEnergy / etot2;
194 f += f1;
195 }
196 G4double x1 = 1.0 + x;
197 G4double gg = 1.0 / (x1 * x1);
198 if(0.5 == fSpin)
199 {
200 G4double x2 =
201 0.5 * electron_mass_c2 * deltaKinEnergy / (fMass * fMass);
202 gg *= (1.0 + fMagMoment2 * (x2 - f1 / f) / (1.0 + x2));
203 }
204 if(gg > 1.0)
205 {
206 G4cout << "### G4BetheBlochModel in Adjoint Sim WARNING: gg= " << gg
207 << G4endl;
208 gg = 1.;
209 }
210 dSigmadEprod *= gg;
211 }
212 }
213 }
214 return dSigmadEprod;
215}
Z
const G4int Z[17]
Definition: G4WaterStopping.cc:51
A
const G4double A[17]
Definition: G4WaterStopping.cc:53
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4AdjointIonIonisationModel::GetSecondAdjEnergyMaxForProdToProj
G4double GetSecondAdjEnergyMaxForProdToProj(G4double primAdjEnergy) override
Definition: G4AdjointIonIonisationModel.cc:325
G4AdjointIonIonisationModel::GetSecondAdjEnergyMinForProdToProj
G4double GetSecondAdjEnergyMinForProdToProj(G4double primAdjEnergy) override
Definition: G4AdjointIonIonisationModel.cc:332

◆ GetSecondAdjEnergyMaxForProdToProj()

G4double G4AdjointIonIonisationModel::GetSecondAdjEnergyMaxForProdToProj ( G4double  primAdjEnergy)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 325 of file G4AdjointIonIonisationModel.cc.

327{
328 return GetHighEnergyLimit();
329}
G4VEmAdjointModel::GetHighEnergyLimit
G4double GetHighEnergyLimit()
Definition: G4VEmAdjointModel.hh:152

Referenced by DiffCrossSectionPerAtomPrimToSecond().

◆ GetSecondAdjEnergyMaxForScatProjToProj()

G4double G4AdjointIonIonisationModel::GetSecondAdjEnergyMaxForScatProjToProj ( G4double  primAdjEnergy)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 310 of file G4AdjointIonIonisationModel.cc.

312{
313 return primAdjEnergy * fOnePlusRatio2 /
314 (fOneMinusRatio2 - 2. * fRatio * primAdjEnergy / fMass);
315}

◆ GetSecondAdjEnergyMinForProdToProj()

G4double G4AdjointIonIonisationModel::GetSecondAdjEnergyMinForProdToProj ( G4double  primAdjEnergy)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 332 of file G4AdjointIonIonisationModel.cc.

334{
335 return (2. * primAdjEnergy - 4. * fMass +
336 std::sqrt(4. * primAdjEnergy * primAdjEnergy + 16. * fMass * fMass +
337 8. * primAdjEnergy * fMass * (1. / fRatio + fRatio))) /
338 4.;
339}

Referenced by DiffCrossSectionPerAtomPrimToSecond().

◆ GetSecondAdjEnergyMinForScatProjToProj()

G4double G4AdjointIonIonisationModel::GetSecondAdjEnergyMinForScatProjToProj ( G4double  primAdjEnergy,
G4double  tcut = 0. 
)
overridevirtual

Reimplemented from G4VEmAdjointModel.

Definition at line 318 of file G4AdjointIonIonisationModel.cc.

320{
321 return primAdjEnergy + tcut;
322}

◆ operator=()

G4AdjointIonIonisationModel & G4AdjointIonIonisationModel::operator= ( const G4AdjointIonIonisationModel &  right)
delete

◆ SampleSecondaries()

void G4AdjointIonIonisationModel::SampleSecondaries ( const G4Track &  aTrack,
G4bool  isScatProjToProj,
G4ParticleChange *  fParticleChange 
)
overridevirtual

Implements G4VEmAdjointModel.

Definition at line 65 of file G4AdjointIonIonisationModel.cc.

68{
69 const G4DynamicParticle* theAdjointPrimary = aTrack.GetDynamicParticle();
70
71 // Elastic inverse scattering
72 G4double adjointPrimKinEnergy = theAdjointPrimary->GetKineticEnergy();
73 G4double adjointPrimP = theAdjointPrimary->GetTotalMomentum();
74
75 if(adjointPrimKinEnergy > GetHighEnergyLimit() * 0.999)
76 {
77 return;
78 }
79
80 // Sample secondary energy
81 G4double projectileKinEnergy =
82 SampleAdjSecEnergyFromCSMatrix(adjointPrimKinEnergy, isScatProjToProj);
83 // Caution !!!this weight correction should be always applied
84 CorrectPostStepWeight(fParticleChange, aTrack.GetWeight(),
85 adjointPrimKinEnergy, projectileKinEnergy,
86 isScatProjToProj);
87
88 // Kinematics:
89 // we consider a two body elastic scattering for the forward processes where
90 // the projectile knock on an e- at rest and gives it part of its energy
91 G4double projectileM0 = fAdjEquivDirectPrimPart->GetPDGMass();
92 G4double projectileTotalEnergy = projectileM0 + projectileKinEnergy;
93 G4double projectileP2 =
94 projectileTotalEnergy * projectileTotalEnergy - projectileM0 * projectileM0;
95
96 // Companion
97 G4double companionM0 = fAdjEquivDirectPrimPart->GetPDGMass();
98 if(isScatProjToProj)
99 {
100 companionM0 = fAdjEquivDirectSecondPart->GetPDGMass();
101 }
102 G4double companionTotalEnergy =
103 companionM0 + projectileKinEnergy - adjointPrimKinEnergy;
104 G4double companionP2 =
105 companionTotalEnergy * companionTotalEnergy - companionM0 * companionM0;
106
107 // Projectile momentum
108 G4double P_parallel =
109 (adjointPrimP * adjointPrimP + projectileP2 - companionP2) /
110 (2. * adjointPrimP);
111 G4double P_perp = std::sqrt(projectileP2 - P_parallel * P_parallel);
112 G4ThreeVector dir_parallel = theAdjointPrimary->GetMomentumDirection();
113 G4double phi = G4UniformRand() * twopi;
114 G4ThreeVector projectileMomentum =
115 G4ThreeVector(P_perp * std::cos(phi), P_perp * std::sin(phi), P_parallel);
116 projectileMomentum.rotateUz(dir_parallel);
117
118 if(!isScatProjToProj)
119 { // kill the primary and add a secondary
120 fParticleChange->ProposeTrackStatus(fStopAndKill);
121 fParticleChange->AddSecondary(
122 new G4DynamicParticle(fAdjEquivDirectPrimPart, projectileMomentum));
123 }
124 else
125 {
126 fParticleChange->ProposeEnergy(projectileKinEnergy);
127 fParticleChange->ProposeMomentumDirection(projectileMomentum.unit());
128 }
129}
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:36
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::rotateUz
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
G4AdjointIonIonisationModel::CorrectPostStepWeight
void CorrectPostStepWeight(G4ParticleChange *fParticleChange, G4double old_weight, G4double adjointPrimKinEnergy, G4double projectileKinEnergy, G4bool isScatProjToProj) override
Definition: G4AdjointIonIonisationModel.cc:228
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4DynamicParticle::GetTotalMomentum
G4double GetTotalMomentum() const
G4ParticleChange::AddSecondary
void AddSecondary(G4Track *aSecondary)
Definition: G4ParticleChange.cc:100
G4ParticleChange::ProposeEnergy
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4Track::GetWeight
G4double GetWeight() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4VEmAdjointModel::fAdjEquivDirectPrimPart
G4ParticleDefinition * fAdjEquivDirectPrimPart
Definition: G4VEmAdjointModel.hh:259
G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix
G4double SampleAdjSecEnergyFromCSMatrix(std::size_t MatrixIndex, G4double prim_energy, G4bool isScatProjToProj)
Definition: G4VEmAdjointModel.cc:413
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)

◆ SetIon()

void G4AdjointIonIonisationModel::SetIon ( G4ParticleDefinition *  adj_ion,
G4ParticleDefinition *  fwd_ion 
)

Definition at line 218 of file G4AdjointIonIonisationModel.cc.

220{
221 fDirectPrimaryPart = fwd_ion;
222 fAdjEquivDirectPrimPart = adj_ion;
223
224 DefineProjectileProperty();
225}

◆ SetUseOnlyBragg()

void G4AdjointIonIonisationModel::SetUseOnlyBragg ( G4bool  aBool)
inline

Definition at line 76 of file G4AdjointIonIonisationModel.hh.

76{ fUseOnlyBragg = aBool; }
The documentation for this class was generated from the following files: