Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Protected Member Functions | List of all members
G4NuclearStopping Class Reference

#include <G4NuclearStopping.hh>

Inheritance diagram for G4NuclearStopping:
Inheritance graph
[legend]
Collaboration diagram for G4NuclearStopping:
Collaboration graph
[legend]

Public Member Functions

 G4NuclearStopping (const G4String &processName="nuclearStopping")
 
virtual ~G4NuclearStopping ()
 
G4bool IsApplicable (const G4ParticleDefinition &p) final
 
G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection) final
 
G4VParticleChange * AlongStepDoIt (const G4Track &track, const G4Step &step) final
 
virtual void PrintInfo () final
 
- Public Member Functions inherited from G4VEmProcess
 G4VEmProcess (const G4String &name, G4ProcessType type=fElectromagnetic)
 
virtual ~G4VEmProcess ()
 
virtual void ProcessDescription (std::ostream &outFile) const
 
virtual void PreparePhysicsTable (const G4ParticleDefinition &) override
 
virtual void BuildPhysicsTable (const G4ParticleDefinition &) override
 
virtual void StartTracking (G4Track *) override
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &) override
 
virtual G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false) override
 
virtual G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii) override
 
G4double CrossSectionPerVolume (G4double kineticEnergy, const G4MaterialCutsCouple *couple)
 
G4double ComputeCrossSectionPerAtom (G4double kineticEnergy, G4double Z, G4double A=0., G4double cut=0.0)
 
G4double MeanFreePath (const G4Track &track)
 
G4double GetLambda (G4double &kinEnergy, const G4MaterialCutsCouple *couple)
 
void SetLambdaBinning (G4int nbins)
 
void SetMinKinEnergy (G4double e)
 
void SetMinKinEnergyPrim (G4double e)
 
void SetMaxKinEnergy (G4double e)
 
G4PhysicsTable * LambdaTable () const
 
G4PhysicsTable * LambdaTablePrim () const
 
const G4ParticleDefinition * Particle () const
 
const G4ParticleDefinition * SecondaryParticle () const
 
G4VEmModel * SelectModelForMaterial (G4double kinEnergy, size_t &idxRegion) const
 
void AddEmModel (G4int, G4VEmModel *, const G4Region *region=nullptr)
 
G4VEmModel * EmModel (G4int index=1) const
 
void SetEmModel (G4VEmModel *, G4int index=1)
 
void UpdateEmModel (const G4String &, G4double, G4double)
 
G4int GetNumberOfModels () const
 
G4int GetNumberOfRegionModels (size_t couple_index) const
 
G4VEmModel * GetRegionModel (G4int idx, size_t couple_index) const
 
G4VEmModel * GetModelByIndex (G4int idx=0, G4bool ver=false) const
 
const G4VEmModel * GetCurrentModel () const
 
const G4Element * GetCurrentElement () const
 
void SetCrossSectionBiasingFactor (G4double f, G4bool flag=true)
 
G4double CrossSectionBiasingFactor () const
 
void ActivateForcedInteraction (G4double length=0.0, const G4String &r="", G4bool flag=true)
 
void ActivateSecondaryBiasing (const G4String &region, G4double factor, G4double energyLimit)
 
void SetIntegral (G4bool val)
 
void SetBuildTableFlag (G4bool val)
 
- Public Member Functions inherited from G4VDiscreteProcess
 G4VDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VDiscreteProcess (G4VDiscreteProcess &)
 
virtual ~G4VDiscreteProcess ()
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4int operator== (const G4VProcess &right) const
 
G4int operator!= (const G4VProcess &right) const
 
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)
 
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 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
 
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 &)
 

Protected Member Functions

void InitialiseProcess (const G4ParticleDefinition *) final
 
- Protected Member Functions inherited from G4VEmProcess
virtual G4double MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *)
 
G4VEmModel * SelectModel (G4double &kinEnergy, size_t index)
 
virtual G4double GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
 
G4PhysicsVector * LambdaPhysicsVector (const G4MaterialCutsCouple *)
 
G4int LambdaBinning () const
 
G4double MinKinEnergy () const
 
G4double MaxKinEnergy () const
 
G4double PolarAngleLimit () const
 
G4bool IsIntegral () const
 
G4double RecalculateLambda (G4double kinEnergy, const G4MaterialCutsCouple *couple)
 
G4ParticleChangeForGamma * GetParticleChange ()
 
void SetParticle (const G4ParticleDefinition *p)
 
void SetSecondaryParticle (const G4ParticleDefinition *p)
 
size_t CurrentMaterialCutsCoupleIndex () const
 
G4double GetGammaEnergyCut ()
 
G4double GetElectronEnergyCut ()
 
void SetStartFromNullFlag (G4bool val)
 
void SetSplineFlag (G4bool val)
 
const G4Element * GetTargetElement () const
 
const G4Isotope * GetTargetIsotope () const
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Attributes inherited from G4VEmProcess
G4ParticleChangeForGamma fParticleChange
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager
 
G4VParticleChange * pParticleChange
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft
 
G4double currentInteractionLength
 
G4double theInitialNumberOfInteractionLength
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType
 
G4int theProcessSubType
 
G4double thePILfactor
 
G4bool enableAtRestDoIt
 
G4bool enableAlongStepDoIt
 
G4bool enablePostStepDoIt
 
G4int verboseLevel
 

Detailed Description

Definition at line 66 of file G4NuclearStopping.hh.

Constructor & Destructor Documentation

G4NuclearStopping::G4NuclearStopping ( const G4String &  processName = "nuclearStopping")
explicit

Definition at line 56 of file G4NuclearStopping.cc.

57  : G4VEmProcess(processName)
58 {
59  isInitialized = false;
60  SetProcessSubType(fNuclearStopping);
61  SetBuildTableFlag(false);
62  enableAlongStepDoIt = true;
63  enablePostStepDoIt = false;
64  modelICRU49 = nullptr;
65 }
G4VEmProcess::G4VEmProcess
G4VEmProcess(const G4String &name, G4ProcessType type=fElectromagnetic)
Definition: G4VEmProcess.cc:91
G4VEmProcess::SetBuildTableFlag
void SetBuildTableFlag(G4bool val)
Definition: G4VEmProcess.hh:641
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:432
G4VProcess::enablePostStepDoIt
G4bool enablePostStepDoIt
Definition: G4VProcess.hh:352
G4VProcess::enableAlongStepDoIt
G4bool enableAlongStepDoIt
Definition: G4VProcess.hh:351
fNuclearStopping
Definition: G4EmProcessSubType.hh:53

Here is the call graph for this function:

G4NuclearStopping::~G4NuclearStopping ( )
virtual

Definition at line 69 of file G4NuclearStopping.cc.

70 {}

Member Function Documentation

G4VParticleChange * G4NuclearStopping::AlongStepDoIt ( const G4Track &  track,
const G4Step &  step 
)
finalvirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 107 of file G4NuclearStopping.cc.

109 {
110  nParticleChange.InitializeForAlongStep(track);
111 
112  // this line only valid if nuclear stopping
113  // is computed after G4ionIonisation process
114  nParticleChange.SetProposedCharge(step.GetPostStepPoint()->GetCharge());
115 
116  G4double T2 = step.GetPostStepPoint()->GetKineticEnergy();
117 
118  const G4ParticleDefinition* part = track.GetParticleDefinition();
119  G4double Z = std::fabs(part->GetPDGCharge()/eplus);
120 
121  if(T2 > 0.0 && T2*proton_mass_c2 < Z*Z*MeV*part->GetPDGMass()) {
122 
123  G4double length = step.GetStepLength();
124  if(length > 0.0) {
125 
126  // primary
127  G4double T1= step.GetPreStepPoint()->GetKineticEnergy();
128  G4double T = 0.5*(T1 + T2);
129  const G4MaterialCutsCouple* couple = track.GetMaterialCutsCouple();
130  G4VEmModel* mod = SelectModel(T, couple->GetIndex());
131 
132  // sample stopping
133  if(modelICRU49) { modelICRU49->SetFluctuationFlag(true); }
134  G4double nloss =
135  length*mod->ComputeDEDXPerVolume(couple->GetMaterial(), part, T);
136  if(nloss > T1) { nloss = T1; }
137  nParticleChange.SetProposedKineticEnergy(T1 - nloss);
138  nParticleChange.ProposeLocalEnergyDeposit(nloss);
139  nParticleChange.ProposeNonIonizingEnergyDeposit(nloss);
140  if(modelICRU49) { modelICRU49->SetFluctuationFlag(false); }
141  }
142  }
143  return &nParticleChange;
144 }
G4VEmProcess::SelectModel
G4VEmModel * SelectModel(G4double &kinEnergy, size_t index)
Definition: G4VEmProcess.hh:465
G4Step::GetStepLength
G4double GetStepLength() const
G4Track::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4ParticleChangeForLoss::InitializeForAlongStep
void InitializeForAlongStep(const G4Track &)
Definition: G4ParticleChangeForLoss.hh:232
G4VParticleChange::ProposeLocalEnergyDeposit
void ProposeLocalEnergyDeposit(G4double anEnergyPart)
G4Step::GetPreStepPoint
G4StepPoint * GetPreStepPoint() const
G4VParticleChange::ProposeNonIonizingEnergyDeposit
void ProposeNonIonizingEnergyDeposit(G4double anEnergyPart)
eplus
static constexpr double eplus
Definition: G4SIunits.hh:199
G4StepPoint::GetCharge
G4double GetCharge() const
G4ParticleChangeForLoss::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForLoss.hh:150
G4Track::GetParticleDefinition
const G4ParticleDefinition * GetParticleDefinition() const
G4VEmModel::ComputeDEDXPerVolume
virtual G4double ComputeDEDXPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
Definition: G4VEmModel.cc:248
G4Step::GetPostStepPoint
G4StepPoint * GetPostStepPoint() const
G4ParticleChangeForLoss::SetProposedCharge
void SetProposedCharge(G4double theCharge)
Definition: G4ParticleChangeForLoss.hh:165
G4StepPoint::GetKineticEnergy
G4double GetKineticEnergy() const
G4double
double G4double
Definition: G4Types.hh:76
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:124
Z
G4int Z
Definition: G4ParticleHPJENDLHEData.cc:262
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:150

Here is the call graph for this function:

G4double G4NuclearStopping::AlongStepGetPhysicalInteractionLength ( const G4Track &  track,
G4double  previousStepSize,
G4double  currentMinimumStep,
G4double &  proposedSafety,
G4GPILSelection *  selection 
)
finalvirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 98 of file G4NuclearStopping.cc.

100 {
101  *selection = NotCandidateForSelection;
102  return DBL_MAX;
103 }
NotCandidateForSelection
Definition: G4GPILSelection.hh:57
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
void G4NuclearStopping::InitialiseProcess ( const G4ParticleDefinition *  )
finalprotectedvirtual

Implements G4VEmProcess.

Definition at line 81 of file G4NuclearStopping.cc.

82 {
83  if(!isInitialized) {
84  isInitialized = true;
85 
86  if(!EmModel(1)) {
87  modelICRU49 = new G4ICRU49NuclearStoppingModel();
88  SetEmModel(modelICRU49);
89  }
90  AddEmModel(1, EmModel());
91  EmModel()->SetParticleChange(&nParticleChange);
92  }
93 }
G4VEmModel::SetParticleChange
void SetParticleChange(G4VParticleChange *, G4VEmFluctuationModel *f=nullptr)
Definition: G4VEmModel.cc:418
G4VEmProcess::EmModel
G4VEmModel * EmModel(G4int index=1) const
Definition: G4VEmProcess.cc:224
G4VEmProcess::SetEmModel
void SetEmModel(G4VEmModel *, G4int index=1)
Definition: G4VEmProcess.cc:213
G4VEmProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
Definition: G4VEmProcess.cc:203

Here is the call graph for this function:

G4bool G4NuclearStopping::IsApplicable ( const G4ParticleDefinition &  p)
finalvirtual

Implements G4VEmProcess.

Definition at line 74 of file G4NuclearStopping.cc.

75 {
76  return (p.GetPDGCharge() != 0.0 && !p.IsShortLived());
77 }
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:124

Here is the call graph for this function:

void G4NuclearStopping::PrintInfo ( )
finalvirtual

Implements G4VEmProcess.

Definition at line 148 of file G4NuclearStopping.cc.

149 {}
The documentation for this class was generated from the following files: