158{
159
160
161
162 G4String rName = track.GetStep()->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetRegion()->GetName();
163
164 if( rName != fEnvelopeName )
165 {
167 {
168 G4cout<<
"Go out from G4NeutrinoElectronProcess::PostStepDoIt: wrong volume "<<
G4endl;
169 }
171 }
174 G4double weight = track.GetWeight();
176
177 if( track.GetTrackStatus() !=
fAlive )
178 {
180 }
181
182
183 if (track.GetTrackStatus() !=
fAlive &&
185 {
189 {
191 ed << "G4HadronicProcess: track in unusable state - "
193 ed <<
"G4HadronicProcess: returning unchanged track " <<
G4endl;
196 }
197
199 }
200
201
203
204 G4double kineticEnergy = track.GetKineticEnergy();
207
208
209
210
212
215
216
217
218 const G4StepPoint* pPostStepPoint = step.GetPostStepPoint();
222
223 if( fNuEleCcBias > 1.0 || fNuEleNcBias > 1.0)
224 {
229
232
233 G4double forward = track.GetVolume()->GetLogicalVolume()->GetSolid()->DistanceToOut(localP, localV);
234 G4double backward = track.GetVolume()->GetLogicalVolume()->GetSolid()->DistanceToOut(localP, -localV);
235
236 G4double distance = forward+backward;
237
238
239
240
241
242
244
245 newPosition =
position + range*direction;
246
248
250 }
254
255
257
258 try
259 {
261 *targNucleus);
262 }
264 {
268 ed <<
" PostStepDoIt failed on element selection" <<
G4endl;
269 G4Exception(
"G4NeutrinoElectronProcess::PostStepDoIt",
"had003",
271 }
272
274
276 {
277
279
281
283
285
287
289 }
290 else
291 {
292
294
295 size_t idx = track.GetMaterialCutsCouple()->GetIndex();
296
298
300
302 {
303 G4cout <<
"G4NeutrinoElectronProcess::PostStepDoIt for "
305 <<
" in " << material->
GetName()
308 }
309 try
310 {
312 }
314 {
318 ed <<
"Target element "<< elm->
GetName()<<
" Z= "
322 ed <<
" ApplyYourself failed" <<
G4endl;
323 G4Exception(
"G4NeutrinoElectronProcess::PostStepDoIt",
"had006",
325 }
326
327
332
334 {
338 << " dir= " << outdir
340 }
341
342
345
346 if(efinal < 0.0) { efinal = 0.0; }
347 if(edep < 0.0) { edep = 0.0; }
348
349
350
351
352 if(efinal <= lowestEnergy)
353 {
354 edep += efinal;
355 efinal = 0.0;
356 }
357
358
360
362
363 if(efinal > 0.0)
364 {
368 }
369 else
370 {
372 {
374 }
375 else
376 {
378 }
380 }
381
382
384
385
386
388 {
390
392 {
395
396
397
398
401
402
403
405
406
407
409 track.GetPosition());
413 }
414 else
415 {
417 delete p;
418 }
419 }
423 }
425}
@ fKillTrackAndSecondaries
G4GLOB_DLL std::ostream G4cout
Hep3Vector & rotateUz(const Hep3Vector &)
Hep3Vector & rotate(double, 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)
void Initialise(const G4Track &aT)
G4LorentzRotation & GetTrafoToLab()
G4DynamicParticle * GetParticle()
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()
G4ParticleChange * theTotalResult
std::vector< G4HadronicInteraction * > & GetHadronicInteractionList()
const G4String & GetName() const
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
void AddSecondary(G4Track *aSecondary)
void ProposePosition(G4double x, G4double y, G4double z)
void Initialize(const G4Track &) override
void ProposeEnergy(G4double finalEnergy)
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
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 ReLocateWithinVolume(const G4ThreeVector &pGlobalPoint)
const G4VTouchable * GetTouchable() const
const G4ThreeVector & GetPosition() const
void SetWeight(G4double aValue)
void SetTouchableHandle(const G4TouchableHandle &apValue)
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
void ProposeTrackStatus(G4TrackStatus status)
void ProposeNonIonizingEnergyDeposit(G4double anEnergyPart)
void ProposeWeight(G4double finalWeight)
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
void SetNumberOfSecondaries(G4int totSecondaries)
G4TrackStatus GetTrackStatus() const
void ClearNumberOfInteractionLengthLeft()
virtual const G4ThreeVector & GetTranslation(G4int depth=0) const =0
virtual const G4RotationMatrix * GetRotation(G4int depth=0) const =0