64{
69
70
72
73 G4double kineticEnergy = track.GetKineticEnergy();
75 if(kineticEnergy == 0.0 || track.GetTrackStatus() !=
fAlive) {
77 }
78
83
84
87
88
92
93 if(fDiffraction)
94 {
99
101 {
102 try
103 {
105 }
107 {
111 ed <<
"Target element "<< elm->
GetName()<<
" Z= "
115 ed <<
" ApplyYourself failed" <<
G4endl;
116 G4Exception(
"G4HadronElasticProcess::PostStepDoIt",
"had006",
118 }
119
120 result =
CheckResult(theProj, *targNucleus, result);
121
124
126
129 }
131 }
132 }
133
134
135 try
136 {
138 }
140 {
143 ed <<
"Target element "<< elm->
GetName()<<
" Z= "
146 DumpState(track,
"ChooseHadronicInteraction",ed);
147 ed <<
" No HadronicInteraction found out" <<
G4endl;
148 G4Exception(
"G4HadronElasticProcess::PostStepDoIt",
"had005",
150 }
151
152 size_t idx = track.GetMaterialCutsCouple()->GetIndex();
156
158 G4cout <<
"G4HadronElasticProcess::PostStepDoIt for "
160 <<
" in " << material->
GetName()
163 <<
" Tcut(MeV)= " << tcut <<
G4endl;
164 }
165
166 try
167 {
169 }
171 {
175 ed <<
"Target element "<< elm->
GetName()<<
" Z= "
179 ed <<
" ApplyYourself failed" <<
G4endl;
180 G4Exception(
"G4HadronElasticProcess::PostStepDoIt",
"had006",
182 }
183
184
185
186
187
188
189
192
197 << " dir= " << outdir
199 }
200
201
204
205
207
208 if(efinal > 0.0) {
211 } else {
216 }
217
218
220
221
224
228
230
232
233
235 track.GetPosition());
239
240 } else {
242 delete p;
243 }
244 }
248
250}
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
std::ostringstream G4ExceptionDescription
G4GLOB_DLL std::ostream G4cout
Hep3Vector & rotateUz(const Hep3Vector &)
const G4Element * SampleZandA(const G4DynamicParticle *, const G4Material *, G4Nucleus &target)
void SetMomentumDirection(const G4ThreeVector &aDirection)
const G4ThreeVector & GetMomentumDirection() const
G4ParticleDefinition * GetDefinition() const
G4double GetKineticEnergy() const
const G4String & GetName() const
G4double GetEnergyChange() const
void SetTrafoToLab(const G4LorentzRotation &aT)
G4double GetLocalEnergyDeposit() const
const G4ThreeVector & GetMomentumChange() const
std::size_t GetNumberOfSecondaries() const
G4HadSecondary * GetSecondary(size_t i)
G4DynamicParticle * GetParticle()
void Report(std::ostream &aS) const
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
const G4String & GetModelName() const
void SetRecoilEnergyThreshold(G4double val)
void FillResult(G4HadFinalState *aR, const G4Track &aT)
G4Nucleus * GetTargetNucleusPointer()
G4HadFinalState * CheckResult(const G4HadProjectile &thePro, const G4Nucleus &targetNucleus, G4HadFinalState *result)
G4ParticleChange * theTotalResult
G4HadronicInteraction * ChooseHadronicInteraction(const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement)
void CheckEnergyMomentumConservation(const G4Track &, const G4Nucleus &)
G4CrossSectionDataStore * GetCrossSectionDataStore()
void DumpState(const G4Track &, const G4String &, G4ExceptionDescription &)
const G4String & GetName() const
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
void AddSecondary(G4Track *aSecondary)
void ProposeEnergy(G4double finalEnergy)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
virtual void Initialize(const G4Track &)
G4ProcessManager * GetProcessManager() const
const G4String & GetParticleName() const
G4ProcessVector * GetAtRestProcessVector(G4ProcessVectorTypeIndex typ=typeGPIL) const
const std::vector< G4double > * GetEnergyCutsVector(std::size_t pcIdx) const
static G4ProductionCutsTable * GetProductionCutsTable()
void SetWeight(G4double aValue)
void SetTouchableHandle(const G4TouchableHandle &apValue)
virtual G4double ComputeRatio(const G4ParticleDefinition *, G4double kinEnergy, G4int, G4int)=0
void ProposeTrackStatus(G4TrackStatus status)
void ProposeNonIonizingEnergyDeposit(G4double anEnergyPart)
void ProposeWeight(G4double finalWeight)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetNumberOfSecondaries(G4int totSecondaries)
void ClearNumberOfInteractionLengthLeft()