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
G4TauNeutrinoNucleusProcess Class Reference

#include <G4TauNeutrinoNucleusProcess.hh>

+ Inheritance diagram for G4TauNeutrinoNucleusProcess:

Public Member Functions

 G4TauNeutrinoNucleusProcess (G4String anEnvelopeName, const G4String &procName="tau-neutrino-nucleus")
 
virtual ~G4TauNeutrinoNucleusProcess ()
 
G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep) override
 
void PreparePhysicsTable (const G4ParticleDefinition &) override
 
virtual void SetLowestEnergy (G4double)
 
void ProcessDescription (std::ostream &outFile) const override
 
void SetBiasingFactors (G4double bfCc, G4double bfNc)
 
void SetBiasingFactor (G4double bf)
 
G4double GetMeanFreePath (const G4Track &aTrack, G4double, G4ForceCondition *) override
 
- Public Member Functions inherited from G4HadronicProcess
 G4HadronicProcess (const G4String &processName="Hadronic", G4ProcessType procType=fHadronic)
 
 G4HadronicProcess (const G4String &processName, G4HadronicProcessType subType)
 
 ~G4HadronicProcess () override
 
void RegisterMe (G4HadronicInteraction *a)
 
G4double GetElementCrossSection (const G4DynamicParticle *part, const G4Element *elm, const G4Material *mat=nullptr)
 
G4double GetMicroscopicCrossSection (const G4DynamicParticle *part, const G4Element *elm, const G4Material *mat=nullptr)
 
void StartTracking (G4Track *track) override
 
G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double, G4ForceCondition *) override
 
G4VParticleChange * PostStepDoIt (const G4Track &aTrack, const G4Step &aStep) override
 
void PreparePhysicsTable (const G4ParticleDefinition &) override
 
void BuildPhysicsTable (const G4ParticleDefinition &) override
 
void DumpPhysicsTable (const G4ParticleDefinition &p)
 
void AddDataSet (G4VCrossSectionDataSet *aDataSet)
 
std::vector< G4HadronicInteraction * > & GetHadronicInteractionList ()
 
G4HadronicInteraction * GetHadronicModel (const G4String &)
 
G4HadronicInteraction * GetHadronicInteraction () const
 
G4double GetMeanFreePath (const G4Track &aTrack, G4double, G4ForceCondition *) override
 
const G4Nucleus * GetTargetNucleus () const
 
G4Nucleus * GetTargetNucleusPointer ()
 
const G4Isotope * GetTargetIsotope ()
 
G4double ComputeCrossSection (const G4ParticleDefinition *, const G4Material *, const G4double kinEnergy)
 
G4HadXSType CrossSectionType () const
 
void SetCrossSectionType (G4HadXSType val)
 
void ProcessDescription (std::ostream &outFile) const override
 
void BiasCrossSectionByFactor (G4double aScale)
 
void MultiplyCrossSectionBy (G4double factor)
 
G4double CrossSectionFactor () const
 
void SetIntegral (G4bool val)
 
void SetEpReportLevel (G4int level)
 
void SetEnergyMomentumCheckLevels (G4double relativeLevel, G4double absoluteLevel)
 
std::pair< G4double, G4double > GetEnergyMomentumCheckLevels () const
 
G4CrossSectionDataStore * GetCrossSectionDataStore ()
 
std::vector< G4TwoPeaksHadXS * > * TwoPeaksXS () const
 
std::vector< G4double > * EnergyOfCrossSectionMax () const
 
G4HadronicProcess & operator= (const G4HadronicProcess &right)=delete
 
 G4HadronicProcess (const G4HadronicProcess &)=delete
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &aName, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
G4VDiscreteProcess & operator= (const G4VDiscreteProcess &)=delete
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
virtual G4double GetCrossSection (const G4double, const G4MaterialCutsCouple *)
 
virtual G4double MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4VProcess & operator= (const G4VProcess &)=delete
 
G4bool operator== (const G4VProcess &right) const
 
G4bool operator!= (const G4VProcess &right) const
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChange * AlongStepDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &track, const G4Step &stepData)=0
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0
 
G4double GetCurrentInteractionLength () const
 
void SetPILfactor (G4double value)
 
G4double GetPILfactor () const
 
G4double AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
G4double AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
 
G4double PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
 
virtual G4bool IsApplicable (const G4ParticleDefinition &)
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &)
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &)
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
 
const G4String & GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
 
const G4String & GetProcessName () const
 
G4ProcessType GetProcessType () const
 
void SetProcessType (G4ProcessType)
 
G4int GetProcessSubType () const
 
void SetProcessSubType (G4int)
 
virtual const G4VProcess * GetCreatorProcess () const
 
virtual void StartTracking (G4Track *)
 
virtual void EndTracking ()
 
virtual void SetProcessManager (const G4ProcessManager *)
 
virtual const G4ProcessManager * GetProcessManager ()
 
virtual void ResetNumberOfInteractionLengthLeft ()
 
G4double GetNumberOfInteractionLengthLeft () const
 
G4double GetTotalNumberOfInteractionLengthTraversed () const
 
G4bool isAtRestDoItIsEnabled () const
 
G4bool isAlongStepDoItIsEnabled () const
 
G4bool isPostStepDoItIsEnabled () const
 
virtual void DumpInfo () const
 
virtual void ProcessDescription (std::ostream &outfile) const
 
void SetVerboseLevel (G4int value)
 
G4int GetVerboseLevel () const
 
virtual void SetMasterProcess (G4VProcess *masterP)
 
const G4VProcess * GetMasterProcess () const
 
virtual void BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
 
virtual void PrepareWorkerPhysicsTable (const G4ParticleDefinition &)
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Member Functions inherited from G4HadronicProcess
G4HadronicInteraction * ChooseHadronicInteraction (const G4HadProjectile &aHadProjectile, G4Nucleus &aTargetNucleus, const G4Material *aMaterial, const G4Element *anElement)
 
G4double GetLastCrossSection ()
 
void FillResult (G4HadFinalState *aR, const G4Track &aT)
 
void DumpState (const G4Track &, const G4String &, G4ExceptionDescription &)
 
G4HadFinalState * CheckResult (const G4HadProjectile &thePro, const G4Nucleus &targetNucleus, G4HadFinalState *result)
 
void CheckEnergyMomentumConservation (const G4Track &, const G4Nucleus &)
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 
- Protected Attributes inherited from G4HadronicProcess
G4HadProjectile thePro
 
G4ParticleChange * theTotalResult
 
G4CrossSectionDataStore * theCrossSectionDataStore
 
G4double fWeight = 1.0
 
G4double aScaleFactor = 1.0
 
G4double theLastCrossSection = 0.0
 
G4double mfpKinEnergy = DBL_MAX
 
G4long epReportLevel = 0
 
G4HadXSType fXSType = fHadNoIntegral
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager = nullptr
 
G4VParticleChange * pParticleChange = nullptr
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft = -1.0
 
G4double currentInteractionLength = -1.0
 
G4double theInitialNumberOfInteractionLength = -1.0
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType = fNotDefined
 
G4int theProcessSubType = -1
 
G4double thePILfactor = 1.0
 
G4int verboseLevel = 0
 
G4bool enableAtRestDoIt = true
 
G4bool enableAlongStepDoIt = true
 
G4bool enablePostStepDoIt = true
 

Detailed Description

Definition at line 49 of file G4TauNeutrinoNucleusProcess.hh.

Constructor & Destructor Documentation

◆ G4TauNeutrinoNucleusProcess()

G4TauNeutrinoNucleusProcess::G4TauNeutrinoNucleusProcess ( G4String  anEnvelopeName,
const G4String &  procName = "tau-neutrino-nucleus" 
)

Definition at line 64 of file G4TauNeutrinoNucleusProcess.cc.

65 : G4HadronicProcess( pName, fHadronInelastic ), isInitialised(false), fBiased(true) // fHadronElastic???
66{
67 lowestEnergy = 1.*keV;
68 fEnvelope = nullptr;
69 fEnvelopeName = anEnvelopeName;
70 fTotXsc = nullptr; // new G4TauNeutrinoNucleusTotXsc();
71 fNuNuclCcBias=1.;
72 fNuNuclNcBias=1.;
73 fNuNuclTotXscBias=1.;
74 safetyHelper = G4TransportationManager::GetTransportationManager()->GetSafetyHelper();
75 safetyHelper->InitialiseHelper();
76}
fHadronInelastic
@ fHadronInelastic
Definition: G4HadronicProcessType.hh:45
G4SafetyHelper::InitialiseHelper
void InitialiseHelper()
Definition: G4SafetyHelper.cc:63
G4TransportationManager::GetTransportationManager
static G4TransportationManager * GetTransportationManager()
Definition: G4TransportationManager.cc:104
G4TransportationManager::GetSafetyHelper
G4SafetyHelper * GetSafetyHelper() const

◆ ~G4TauNeutrinoNucleusProcess()

G4TauNeutrinoNucleusProcess::~G4TauNeutrinoNucleusProcess ( )
virtual

Definition at line 78 of file G4TauNeutrinoNucleusProcess.cc.

79{
80 if( fTotXsc ) delete fTotXsc;
81}

Member Function Documentation

◆ GetMeanFreePath()

G4double G4TauNeutrinoNucleusProcess::GetMeanFreePath ( const G4Track &  aTrack,
G4double  ,
G4ForceCondition *   
)
overridevirtual

Reimplemented from G4HadronicProcess.

Definition at line 106 of file G4TauNeutrinoNucleusProcess.cc.

108{
109 //G4cout << "GetMeanFreePath " << aTrack.GetDefinition()->GetParticleName()
110 // << " Ekin= " << aTrack.GetKineticEnergy() << G4endl;
111 G4String rName = aTrack.GetStep()->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetRegion()->GetName();
112 G4double totxsc(0.);
113
114 if( rName == fEnvelopeName && fNuNuclTotXscBias > 1.)
115 {
116 totxsc = fNuNuclTotXscBias*
117 GetCrossSectionDataStore()->ComputeCrossSection(aTrack.GetDynamicParticle(),
118 aTrack.GetMaterial());
119 }
120 else
121 {
122 totxsc = GetCrossSectionDataStore()->ComputeCrossSection(aTrack.GetDynamicParticle(),
123 aTrack.GetMaterial());
124 }
125 G4double res = (totxsc>0.0) ? 1.0/totxsc : DBL_MAX;
126 //G4cout << " xsection= " << totxsc << G4endl;
127 return res;
128}
G4double
double G4double
Definition: G4Types.hh:83
G4CrossSectionDataStore::ComputeCrossSection
G4double ComputeCrossSection(const G4DynamicParticle *, const G4Material *)
Definition: G4CrossSectionDataStore.cc:68
G4HadronicProcess::GetCrossSectionDataStore
G4CrossSectionDataStore * GetCrossSectionDataStore()
Definition: G4HadronicProcess.hh:337
G4LogicalVolume::GetRegion
G4Region * GetRegion() const
G4Region::GetName
const G4String & GetName() const
G4StepPoint::GetPhysicalVolume
G4VPhysicalVolume * GetPhysicalVolume() const
G4Step::GetPreStepPoint
G4StepPoint * GetPreStepPoint() const
G4Track::GetMaterial
G4Material * GetMaterial() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4Track::GetStep
const G4Step * GetStep() const
G4VPhysicalVolume::GetLogicalVolume
G4LogicalVolume * GetLogicalVolume() const
DBL_MAX
#define DBL_MAX
Definition: templates.hh:62

◆ PostStepDoIt()

G4VParticleChange * G4TauNeutrinoNucleusProcess::PostStepDoIt ( const G4Track &  aTrack,
const G4Step &  aStep 
)
overridevirtual

! is not needed for models inheriting G4TauNeutrinoNucleus

Reimplemented from G4HadronicProcess.

Definition at line 144 of file G4TauNeutrinoNucleusProcess.cc.

145{
146 // track.GetVolume()->GetLogicalVolume()->GetName()
147 // if( track.GetVolume()->GetLogicalVolume() != fEnvelope )
148
149 G4String rName = track.GetStep()->GetPreStepPoint()->GetPhysicalVolume()->GetLogicalVolume()->GetRegion()->GetName();
150
151 if( rName != fEnvelopeName )
152 {
153 if( verboseLevel > 0 )
154 {
155 G4cout<<"Go out from G4TauNeutrinoNucleusProcess::PostStepDoIt: wrong volume "<<G4endl;
156 }
157 return G4VDiscreteProcess::PostStepDoIt( track, step );
158 }
159 theTotalResult->Clear();
160 theTotalResult->Initialize(track);
161 G4double weight = track.GetWeight();
162 theTotalResult->ProposeWeight(weight);
163
164 if( track.GetTrackStatus() != fAlive )
165 {
166 return theTotalResult;
167 }
168 // Next check for illegal track status
169 //
170 if (track.GetTrackStatus() != fAlive &&
171 track.GetTrackStatus() != fSuspend)
172 {
173 if (track.GetTrackStatus() == fStopAndKill ||
174 track.GetTrackStatus() == fKillTrackAndSecondaries ||
175 track.GetTrackStatus() == fPostponeToNextEvent)
176 {
177 G4ExceptionDescription ed;
178 ed << "G4TauNeutrinoNucleusProcess: track in unusable state - "
179 << track.GetTrackStatus() << G4endl;
180 ed << "G4TauNeutrinoNucleusProcess: returning unchanged track " << G4endl;
181 DumpState(track,"PostStepDoIt",ed);
182 G4Exception("G4TauNeutrinoNucleusProcess::PostStepDoIt", "had004", JustWarning, ed);
183 }
184 // No warning for fStopButAlive which is a legal status here
185 return theTotalResult;
186 }
187
188 // For elastic scattering, _any_ result is considered an interaction
189 ClearNumberOfInteractionLengthLeft();
190
191 G4double kineticEnergy = track.GetKineticEnergy();
192 const G4DynamicParticle* dynParticle = track.GetDynamicParticle();
193 const G4ParticleDefinition* part = dynParticle->GetDefinition();
194 const G4String pName = part->GetParticleName();
195
196 // NOTE: Very low energy scatters were causing numerical (FPE) errors
197 // in earlier releases; these limits have not been changed since.
198
199 if ( kineticEnergy <= lowestEnergy ) return theTotalResult;
200
201 const G4Material* material = track.GetMaterial();
202 G4Nucleus* targNucleus = GetTargetNucleusPointer();
203
204 //////////////// uniform random spread of the neutrino interaction point ////////////
205
206 const G4StepPoint* pPostStepPoint = step.GetPostStepPoint();
207 const G4DynamicParticle* aParticle = track.GetDynamicParticle();
208 G4ThreeVector position = pPostStepPoint->GetPosition(), newPosition=position;
209 G4ParticleMomentum direction = aParticle->GetMomentumDirection();
210
211 if( fNuNuclCcBias > 1.0 || fNuNuclNcBias > 1.0) // = true, if fBiasingfactor != 1., i.e. xsc is biased
212 {
213 const G4RotationMatrix* rotM = pPostStepPoint->GetTouchable()->GetRotation();
214 G4ThreeVector transl = pPostStepPoint->GetTouchable()->GetTranslation();
215 G4AffineTransform transform = G4AffineTransform(rotM,transl);
216 transform.Invert();
217
218 G4ThreeVector localP = transform.TransformPoint(position);
219 G4ThreeVector localV = transform.TransformAxis(direction);
220
221 G4double forward = track.GetVolume()->GetLogicalVolume()->GetSolid()->DistanceToOut(localP, localV);
222 G4double backward = track.GetVolume()->GetLogicalVolume()->GetSolid()->DistanceToOut(localP, -localV);
223
224 G4double distance = forward+backward;
225
226 // G4cout<<distance/cm<<", ";
227
228 // uniform sampling of nu-e interaction point
229 // along neutrino direction in current volume
230
231 G4double range = -backward+G4UniformRand()*distance;
232
233 newPosition = position + range*direction;
234
235 safetyHelper->ReLocateWithinVolume(newPosition);
236
237 theTotalResult->ProposePosition(newPosition); // G4Exception : GeomNav1002
238 }
239 G4HadProjectile theProj( track );
240 G4HadronicInteraction* hadi = nullptr;
241 G4HadFinalState* result = nullptr;
242
243 G4double ccTotRatio = fTotXsc->GetCcTotRatio();
244
245 if( G4UniformRand() < ccTotRatio ) // Cc-model
246 {
247 // Initialize the hadronic projectile from the track
248 thePro.Initialise(track);
249
250 if (pName == "nu_tau" ) hadi = (GetHadronicInteractionList())[0];
251 else hadi = (GetHadronicInteractionList())[2];
252
253 result = hadi->ApplyYourself( thePro, *targNucleus);
254
255 result->SetTrafoToLab(thePro.GetTrafoToLab());
256
257 ClearNumberOfInteractionLengthLeft();
258
259 FillResult(result, track);
260 }
261 else // Nc-model
262 {
263
264 if (pName == "nu_tau" ) hadi = (GetHadronicInteractionList())[1];
265 else hadi = (GetHadronicInteractionList())[3];
266
267 size_t idx = track.GetMaterialCutsCouple()->GetIndex();
268
269 G4double tcut = (*(G4ProductionCutsTable::GetProductionCutsTable()->GetEnergyCutsVector(3)))[idx];
270
271 hadi->SetRecoilEnergyThreshold(tcut);
272
273 if( verboseLevel > 1 )
274 {
275 G4cout << "G4TauNeutrinoNucleusProcess::PostStepDoIt for "
276 << part->GetParticleName()
277 << " in " << material->GetName()
278 << " Target Z= " << targNucleus->GetZ_asInt()
279 << " A= " << targNucleus->GetA_asInt() << G4endl;
280 }
281 try
282 {
283 result = hadi->ApplyYourself( theProj, *targNucleus);
284 }
285 catch(G4HadronicException & aR)
286 {
287 G4ExceptionDescription ed;
288 aR.Report(ed);
289 ed << "Call for " << hadi->GetModelName() << G4endl;
290 ed << " Z= "
291 << targNucleus->GetZ_asInt()
292 << " A= " << targNucleus->GetA_asInt() << G4endl;
293 DumpState(track,"ApplyYourself",ed);
294 ed << " ApplyYourself failed" << G4endl;
295 G4Exception("G4TauNeutrinoNucleusProcess::PostStepDoIt", "had006",
296 FatalException, ed);
297 }
298 // directions
299
300 G4ThreeVector indir = track.GetMomentumDirection();
301 G4double phi = CLHEP::twopi*G4UniformRand();
302 G4ThreeVector it(0., 0., 1.);
303 G4ThreeVector outdir = result->GetMomentumChange();
304
305 if(verboseLevel>1)
306 {
307 G4cout << "Efin= " << result->GetEnergyChange()
308 << " de= " << result->GetLocalEnergyDeposit()
309 << " nsec= " << result->GetNumberOfSecondaries()
310 << " dir= " << outdir
311 << G4endl;
312 }
313 // energies
314
315 G4double edep = result->GetLocalEnergyDeposit();
316 G4double efinal = result->GetEnergyChange();
317
318 if(efinal < 0.0) { efinal = 0.0; }
319 if(edep < 0.0) { edep = 0.0; }
320
321 // NOTE: Very low energy scatters were causing numerical (FPE) errors
322 // in earlier releases; these limits have not been changed since.
323
324 if(efinal <= lowestEnergy)
325 {
326 edep += efinal;
327 efinal = 0.0;
328 }
329 // primary change
330
331 theTotalResult->ProposeEnergy(efinal);
332
333 G4TrackStatus status = track.GetTrackStatus();
334
335 if(efinal > 0.0)
336 {
337 outdir.rotate(phi, it);
338 outdir.rotateUz(indir);
339 theTotalResult->ProposeMomentumDirection(outdir);
340 }
341 else
342 {
343 if( part->GetProcessManager()->GetAtRestProcessVector()->size() > 0)
344 {
345 status = fStopButAlive;
346 }
347 else
348 {
349 status = fStopAndKill;
350 }
351 theTotalResult->ProposeTrackStatus(status);
352 }
353 //G4cout << "Efinal= " << efinal << " TrackStatus= " << status << G4endl;
354
355 theTotalResult->SetNumberOfSecondaries(0);
356
357 // recoil
358
359 if( result->GetNumberOfSecondaries() > 0 )
360 {
361 G4DynamicParticle* p = result->GetSecondary(0)->GetParticle();
362
363 if(p->GetKineticEnergy() > tcut)
364 {
365 theTotalResult->SetNumberOfSecondaries(1);
366 G4ThreeVector pdir = p->GetMomentumDirection();
367
368 // G4cout << "recoil " << pdir << G4endl;
369 //!! is not needed for models inheriting G4TauNeutrinoNucleus
370
371 pdir.rotate(phi, it);
372 pdir.rotateUz(indir);
373
374 // G4cout << "recoil rotated " << pdir << G4endl;
375
376 p->SetMomentumDirection(pdir);
377
378 // in elastic scattering time and weight are not changed
379
380 G4Track* t = new G4Track(p, track.GetGlobalTime(),
381 track.GetPosition());
382 t->SetWeight(weight);
383 t->SetTouchableHandle(track.GetTouchableHandle());
384 theTotalResult->AddSecondary(t);
385 }
386 else
387 {
388 edep += p->GetKineticEnergy();
389 delete p;
390 }
391 }
392 theTotalResult->ProposeLocalEnergyDeposit(edep);
393 theTotalResult->ProposeNonIonizingEnergyDeposit(edep);
394 result->Clear();
395 }
396 return theTotalResult;
397}
JustWarning
@ JustWarning
Definition: G4ExceptionSeverity.hh:65
FatalException
@ FatalException
Definition: G4ExceptionSeverity.hh:61
G4Exception
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *description)
Definition: G4Exception.cc:59
G4ExceptionDescription
std::ostringstream G4ExceptionDescription
Definition: G4Exception.hh:40
G4TrackStatus
G4TrackStatus
Definition: G4TrackStatus.hh:41
fKillTrackAndSecondaries
@ fKillTrackAndSecondaries
Definition: G4TrackStatus.hh:48
fSuspend
@ fSuspend
Definition: G4TrackStatus.hh:50
fAlive
@ fAlive
Definition: G4TrackStatus.hh:42
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
fStopButAlive
@ fStopButAlive
Definition: G4TrackStatus.hh:44
fPostponeToNextEvent
@ fPostponeToNextEvent
Definition: G4TrackStatus.hh:52
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4UniformRand
#define G4UniformRand()
Definition: Randomize.hh:52
CLHEP::Hep3Vector::rotateUz
Hep3Vector & rotateUz(const Hep3Vector &)
Definition: ThreeVector.cc:33
CLHEP::Hep3Vector::rotate
Hep3Vector & rotate(double, const Hep3Vector &)
Definition: ThreeVectorR.cc:24
G4AffineTransform::Invert
G4AffineTransform & Invert()
G4AffineTransform::TransformPoint
G4ThreeVector TransformPoint(const G4ThreeVector &vec) const
G4AffineTransform::TransformAxis
G4ThreeVector TransformAxis(const G4ThreeVector &axis) const
G4DynamicParticle::SetMomentumDirection
void SetMomentumDirection(const G4ThreeVector &aDirection)
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::GetDefinition
G4ParticleDefinition * GetDefinition() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4HadFinalState::GetEnergyChange
G4double GetEnergyChange() const
Definition: G4HadFinalState.hh:54
G4HadFinalState::SetTrafoToLab
void SetTrafoToLab(const G4LorentzRotation &aT)
Definition: G4HadFinalState.hh:74
G4HadFinalState::GetLocalEnergyDeposit
G4double GetLocalEnergyDeposit() const
Definition: G4HadFinalState.hh:84
G4HadFinalState::GetMomentumChange
const G4ThreeVector & GetMomentumChange() const
Definition: G4HadFinalState.hh:59
G4HadFinalState::GetNumberOfSecondaries
std::size_t GetNumberOfSecondaries() const
Definition: G4HadFinalState.hh:51
G4HadFinalState::GetSecondary
G4HadSecondary * GetSecondary(size_t i)
Definition: G4HadFinalState.cc:79
G4HadProjectile::Initialise
void Initialise(const G4Track &aT)
Definition: G4HadProjectile.cc:54
G4HadProjectile::GetTrafoToLab
G4LorentzRotation & GetTrafoToLab()
Definition: G4HadProjectile.hh:96
G4HadSecondary::GetParticle
G4DynamicParticle * GetParticle()
Definition: G4HadSecondary.hh:43
G4HadronicException::Report
void Report(std::ostream &aS) const
Definition: G4HadronicException.cc:52
G4HadronicInteraction::ApplyYourself
virtual G4HadFinalState * ApplyYourself(const G4HadProjectile &aTrack, G4Nucleus &targetNucleus)
Definition: G4HadronicInteraction.cc:63
G4HadronicInteraction::GetModelName
const G4String & GetModelName() const
Definition: G4HadronicInteraction.hh:115
G4HadronicInteraction::SetRecoilEnergyThreshold
void SetRecoilEnergyThreshold(G4double val)
Definition: G4HadronicInteraction.hh:138
G4HadronicProcess::FillResult
void FillResult(G4HadFinalState *aR, const G4Track &aT)
Definition: G4HadronicProcess.cc:507
G4HadronicProcess::thePro
G4HadProjectile thePro
Definition: G4HadronicProcess.hh:218
G4HadronicProcess::GetTargetNucleusPointer
G4Nucleus * GetTargetNucleusPointer()
Definition: G4HadronicProcess.hh:365
G4HadronicProcess::theTotalResult
G4ParticleChange * theTotalResult
Definition: G4HadronicProcess.hh:220
G4HadronicProcess::GetHadronicInteractionList
std::vector< G4HadronicInteraction * > & GetHadronicInteractionList()
Definition: G4HadronicProcess.cc:896
G4HadronicProcess::DumpState
void DumpState(const G4Track &, const G4String &, G4ExceptionDescription &)
Definition: G4HadronicProcess.cc:860
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:172
G4Material::GetMaterial
static G4Material * GetMaterial(const G4String &name, G4bool warning=true)
Definition: G4Material.cc:691
G4Nucleus::GetA_asInt
G4int GetA_asInt() const
Definition: G4Nucleus.hh:99
G4Nucleus::GetZ_asInt
G4int GetZ_asInt() const
Definition: G4Nucleus.hh:105
G4ParticleChange::AddSecondary
void AddSecondary(G4Track *aSecondary)
Definition: G4ParticleChange.cc:100
G4ParticleChange::ProposePosition
void ProposePosition(G4double x, G4double y, G4double z)
G4ParticleChange::Initialize
void Initialize(const G4Track &) override
Definition: G4ParticleChange.cc:107
G4ParticleChange::ProposeEnergy
void ProposeEnergy(G4double finalEnergy)
G4ParticleChange::ProposeMomentumDirection
void ProposeMomentumDirection(G4double Px, G4double Py, G4double Pz)
G4ParticleDefinition::GetProcessManager
G4ProcessManager * GetProcessManager() const
Definition: G4ParticleDefinition.cc:238
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:109
G4ProcessManager::GetAtRestProcessVector
G4ProcessVector * GetAtRestProcessVector(G4ProcessVectorTypeIndex typ=typeGPIL) const
G4ProcessVector::size
std::size_t size() const
G4ProductionCutsTable::GetEnergyCutsVector
const std::vector< G4double > * GetEnergyCutsVector(std::size_t pcIdx) const
Definition: G4ProductionCutsTable.hh:233
G4ProductionCutsTable::GetProductionCutsTable
static G4ProductionCutsTable * GetProductionCutsTable()
Definition: G4ProductionCutsTable.cc:58
G4SafetyHelper::ReLocateWithinVolume
void ReLocateWithinVolume(const G4ThreeVector &pGlobalPoint)
Definition: G4SafetyHelper.cc:141
G4StepPoint::GetTouchable
const G4VTouchable * GetTouchable() const
G4StepPoint::GetPosition
const G4ThreeVector & GetPosition() const
G4Track::SetWeight
void SetWeight(G4double aValue)
G4Track::SetTouchableHandle
void SetTouchableHandle(const G4TouchableHandle &apValue)
G4VDiscreteProcess::PostStepDoIt
virtual G4VParticleChange * PostStepDoIt(const G4Track &, const G4Step &)
Definition: G4VDiscreteProcess.cc:120
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4VParticleChange::ProposeNonIonizingEnergyDeposit
void ProposeNonIonizingEnergyDeposit(G4double anEnergyPart)
G4VParticleChange::ProposeWeight
void ProposeWeight(G4double finalWeight)
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4VParticleChange::SetNumberOfSecondaries
void SetNumberOfSecondaries(G4int totSecondaries)
G4VParticleChange::GetTrackStatus
G4TrackStatus GetTrackStatus() const
G4VProcess::ClearNumberOfInteractionLengthLeft
void ClearNumberOfInteractionLengthLeft()
Definition: G4VProcess.hh:428
G4VProcess::verboseLevel
G4int verboseLevel
Definition: G4VProcess.hh:360
G4VTouchable::GetTranslation
virtual const G4ThreeVector & GetTranslation(G4int depth=0) const =0
G4VTouchable::GetRotation
virtual const G4RotationMatrix * GetRotation(G4int depth=0) const =0

◆ PreparePhysicsTable()

void G4TauNeutrinoNucleusProcess::PreparePhysicsTable ( const G4ParticleDefinition &  part)
overridevirtual

Reimplemented from G4HadronicProcess.

Definition at line 400 of file G4TauNeutrinoNucleusProcess.cc.

401{
402 if(!isInitialised) {
403 isInitialised = true;
404 // if(G4Neutron::Neutron() == &part) { lowestEnergy = 1.e-6*eV; }
405 }
406 G4HadronicProcess::PreparePhysicsTable(part);
407}
G4HadronicProcess::PreparePhysicsTable
void PreparePhysicsTable(const G4ParticleDefinition &) override
Definition: G4HadronicProcess.cc:169

◆ ProcessDescription()

void G4TauNeutrinoNucleusProcess::ProcessDescription ( std::ostream &  outFile) const
overridevirtual

Reimplemented from G4HadronicProcess.

Definition at line 132 of file G4TauNeutrinoNucleusProcess.cc.

133{
134
135 outFile << "G4TauNeutrinoNucleusProcess handles the inelastic scattering of \n"
136 << "tau-neutrino on nucleus by invoking the following model(s) and \n"
137 << "cross section(s).\n";
138
139}

◆ SetBiasingFactor()

void G4TauNeutrinoNucleusProcess::SetBiasingFactor ( G4double  bf)

Definition at line 85 of file G4TauNeutrinoNucleusProcess.cc.

86{
87 fNuNuclTotXscBias = bf;
88
89 fTotXsc = new G4TauNeutrinoNucleusTotXsc();
90 fTotXsc->SetBiasingFactor(bf);
91}

Referenced by G4EmExtraPhysics::ConstructProcess().

◆ SetBiasingFactors()

void G4TauNeutrinoNucleusProcess::SetBiasingFactors ( G4double  bfCc,
G4double  bfNc 
)

Definition at line 95 of file G4TauNeutrinoNucleusProcess.cc.

96{
97 fNuNuclCcBias = bfCc;
98 fNuNuclNcBias = bfNc;
99
100 fTotXsc = new G4TauNeutrinoNucleusTotXsc();
101 // fTotXsc->SetBiasingFactors(bfCc, bfNc);
102}

◆ SetLowestEnergy()

void G4TauNeutrinoNucleusProcess::SetLowestEnergy ( G4double  val)
virtual

Definition at line 410 of file G4TauNeutrinoNucleusProcess.cc.

411{
412 lowestEnergy = val;
413}
The documentation for this class was generated from the following files: