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

#include <G4CoulombScattering.hh>

+ Inheritance diagram for G4CoulombScattering:

Public Member Functions

 G4CoulombScattering (const G4String &name="CoulombScat")
 
 ~G4CoulombScattering () override
 
G4bool IsApplicable (const G4ParticleDefinition &p) final
 
void ProcessDescription (std::ostream &) const override
 
G4CoulombScattering & operator= (const G4CoulombScattering &right)=delete
 
 G4CoulombScattering (const G4CoulombScattering &)=delete
 
- Public Member Functions inherited from G4VEmProcess
 G4VEmProcess (const G4String &name, G4ProcessType type=fElectromagnetic)
 
 ~G4VEmProcess () override
 
virtual G4bool IsApplicable (const G4ParticleDefinition &p) override=0
 
void ProcessDescription (std::ostream &outFile) const override
 
void PreparePhysicsTable (const G4ParticleDefinition &) override
 
void BuildPhysicsTable (const G4ParticleDefinition &) override
 
void StartTracking (G4Track *) override
 
G4double PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
 
G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &) override
 
G4bool StorePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii=false) override
 
G4bool RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &directory, G4bool ascii) override
 
virtual G4VEmProcess * GetEmProcess (const G4String &name)
 
G4double GetLambda (G4double kinEnergy, const G4MaterialCutsCouple *couple, G4double logKinEnergy)
 
G4double GetCrossSection (const G4double kinEnergy, const G4MaterialCutsCouple *couple) override
 
G4double ComputeCrossSectionPerAtom (G4double kineticEnergy, G4double Z, G4double A=0., G4double cut=0.0)
 
G4double MeanFreePath (const G4Track &track)
 
void SetLambdaBinning (G4int nbins)
 
void SetMinKinEnergy (G4double e)
 
void SetMinKinEnergyPrim (G4double e)
 
void SetMaxKinEnergy (G4double e)
 
G4PhysicsTable * LambdaTable () const
 
G4PhysicsTable * LambdaTablePrim () const
 
void SetLambdaTable (G4PhysicsTable *)
 
void SetLambdaTablePrim (G4PhysicsTable *)
 
std::vector< G4double > * EnergyOfCrossSectionMax () const
 
void SetEnergyOfCrossSectionMax (std::vector< G4double > *)
 
G4CrossSectionType CrossSectionType () const
 
void SetCrossSectionType (G4CrossSectionType val)
 
const G4ParticleDefinition * Particle () const
 
const G4ParticleDefinition * SecondaryParticle () const
 
G4VEmModel * SelectModelForMaterial (G4double kinEnergy, size_t idxCouple) const
 
void AddEmModel (G4int, G4VEmModel *, const G4Region *region=nullptr)
 
void SetEmModel (G4VEmModel *, G4int index=0)
 
G4int NumberOfModels () const
 
G4VEmModel * EmModel (size_t index=0) const
 
const G4VEmModel * GetCurrentModel () const
 
G4VEmModel * GetModelByIndex (G4int idx=0, G4bool ver=false) 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 SetEmMasterProcess (const G4VEmProcess *)
 
void SetBuildTableFlag (G4bool val)
 
void CurrentSetup (const G4MaterialCutsCouple *, G4double energy)
 
G4bool UseBaseMaterial () const
 
void BuildLambdaTable ()
 
void StreamInfo (std::ostream &outFile, const G4ParticleDefinition &, G4bool rst=false) const
 
 G4VEmProcess (G4VEmProcess &)=delete
 
G4VEmProcess & operator= (const G4VEmProcess &right)=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 &)
 

Protected Member Functions

void StreamProcessInfo (std::ostream &outFile) const override
 
void InitialiseProcess (const G4ParticleDefinition *) override
 
G4double MinPrimaryEnergy (const G4ParticleDefinition *, const G4Material *) final
 
- Protected Member Functions inherited from G4VEmProcess
virtual void StreamProcessInfo (std::ostream &) const
 
virtual void InitialiseProcess (const G4ParticleDefinition *)=0
 
G4VEmModel * SelectModel (G4double kinEnergy, size_t)
 
G4double GetMeanFreePath (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override
 
G4PhysicsVector * LambdaPhysicsVector (const G4MaterialCutsCouple *)
 
void DefineMaterial (const G4MaterialCutsCouple *couple)
 
G4int LambdaBinning () const
 
G4double MinKinEnergy () const
 
G4double MaxKinEnergy () const
 
G4double PolarAngleLimit () const
 
G4ParticleChangeForGamma * GetParticleChange ()
 
void SetParticle (const G4ParticleDefinition *p)
 
void SetSecondaryParticle (const G4ParticleDefinition *p)
 
size_t CurrentMaterialCutsCoupleIndex () const
 
const G4MaterialCutsCouple * MaterialCutsCouple () const
 
G4bool ApplyCuts () const
 
G4double GetGammaEnergyCut ()
 
G4double GetElectronEnergyCut ()
 
void SetStartFromNullFlag (G4bool val)
 
void SetSplineFlag (G4bool val)
 
const G4Element * GetTargetElement () const
 
const G4Isotope * GetTargetIsotope () const
 
G4int DensityIndex (G4int idx) const
 
G4double DensityFactor (G4int idx) const
 
virtual G4double GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Attributes inherited from G4VEmProcess
const G4MaterialCutsCouple * currentCouple = nullptr
 
const G4Material * currentMaterial = nullptr
 
G4EmBiasingManager * biasManager = nullptr
 
std::vector< G4double > * theEnergyOfCrossSectionMax = nullptr
 
G4double mfpKinEnergy = DBL_MAX
 
G4double preStepKinEnergy = 0.0
 
G4double preStepLogKinEnergy = LOG_EKIN_MIN
 
G4double preStepLambda = 0.0
 
G4int mainSecondaries = 1
 
G4int secID = _EM
 
G4int fluoID = _Fluorescence
 
G4int augerID = _AugerElectron
 
G4int biasID = _EM
 
G4int tripletID = _TripletElectron
 
size_t currentCoupleIndex = 0
 
size_t basedCoupleIndex = 0
 
size_t coupleIdxLambda = 0
 
size_t idxLambda = 0
 
G4bool isTheMaster = true
 
G4bool baseMat = false
 
std::vector< G4DynamicParticle * > secParticles
 
G4ParticleChangeForGamma fParticleChange
 
- 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 54 of file G4CoulombScattering.hh.

Constructor & Destructor Documentation

◆ G4CoulombScattering() [1/2]

G4CoulombScattering::G4CoulombScattering ( const G4String &  name = "CoulombScat")
explicit

Definition at line 60 of file G4CoulombScattering.cc.

61 : G4VEmProcess(name),q2Max(CLHEP::TeV*CLHEP::TeV),isInitialised(false)
62{
63 // G4cout << "G4CoulombScattering constructor "<< G4endl;
64 SetBuildTableFlag(true);
65 SetStartFromNullFlag(false);
66 SetSplineFlag(false);
67 SetCrossSectionType(fEmOnePeak);
68 SetSecondaryParticle(G4Proton::Proton());
69 SetProcessSubType(fCoulombScattering);
70}
fCoulombScattering
@ fCoulombScattering
Definition: G4EmProcessSubType.hh:45
fEmOnePeak
@ fEmOnePeak
Definition: G4EmTableType.hh:59
G4Proton::Proton
static G4Proton * Proton()
Definition: G4Proton.cc:92
G4VEmProcess::SetBuildTableFlag
void SetBuildTableFlag(G4bool val)
Definition: G4VEmProcess.hh:713
G4VEmProcess::SetSecondaryParticle
void SetSecondaryParticle(const G4ParticleDefinition *p)
Definition: G4VEmProcess.hh:735
G4VEmProcess::SetSplineFlag
void SetSplineFlag(G4bool val)
Definition: G4VEmProcess.hh:749
G4VEmProcess::SetCrossSectionType
void SetCrossSectionType(G4CrossSectionType val)
Definition: G4VEmProcess.hh:699
G4VEmProcess::SetStartFromNullFlag
void SetStartFromNullFlag(G4bool val)
Definition: G4VEmProcess.hh:742
G4VProcess::SetProcessSubType
void SetProcessSubType(G4int)
Definition: G4VProcess.hh:410

◆ ~G4CoulombScattering()

G4CoulombScattering::~G4CoulombScattering ( )
overridedefault

◆ G4CoulombScattering() [2/2]

G4CoulombScattering::G4CoulombScattering ( const G4CoulombScattering &  )
delete

Member Function Documentation

◆ InitialiseProcess()

void G4CoulombScattering::InitialiseProcess ( const G4ParticleDefinition *  p)
overrideprotectedvirtual

Implements G4VEmProcess.

Definition at line 85 of file G4CoulombScattering.cc.

86{
87 // second initialisation not allowed for the time being
88 // this means that polar angle limit change will not be appled
89 // after first initialisation
90 if(isInitialised) { return; }
91
92 G4EmParameters* param = G4EmParameters::Instance();
93 G4double a = param->FactorForAngleLimit()*CLHEP::hbarc/CLHEP::fermi;
94 q2Max = 0.5*a*a;
95 G4double theta = param->MscThetaLimit();
96
97 // restricted or non-restricted cross section table
98 G4bool yes = false;
99 if(theta == CLHEP::pi) {
100 yes = true;
101 // for restriced single scattering change cross section shape
102 SetCrossSectionType(fEmIncreasing);
103 }
104 SetStartFromNullFlag(yes);
105 /*
106 G4cout << "### G4CoulombScattering::InitialiseProcess: "
107 << p->GetParticleName()
108 << " Emin(MeV)= " << MinKinEnergy()/MeV
109 << " Emax(TeV)= " << MaxKinEnergy()/TeV
110 << " nbins= " << LambdaBinning()
111 << " theta= " << theta
112 << G4endl;
113 */
114
115 isInitialised = true;
116 G4double mass = p->GetPDGMass();
117 G4String name = p->GetParticleName();
118 //G4cout << name << " type: " << p->GetParticleType()
119 //<< " mass= " << mass << G4endl;
120 yes = true;
121 if (mass > CLHEP::GeV || p->GetParticleType() == "nucleus") {
122 SetBuildTableFlag(false);
123 yes = false;
124 if(name != "GenericIon") { SetVerboseLevel(0); }
125 } else {
126 if(name != "e-" && name != "e+" &&
127 name != "mu+" && name != "mu-" && name != "pi+" &&
128 name != "kaon+" && name != "proton" ) { SetVerboseLevel(0); }
129 }
130
131 if(nullptr == EmModel(0)) {
132 if(yes) { SetEmModel(new G4eCoulombScatteringModel()); }
133 else { SetEmModel(new G4IonCoulombScatteringModel()); }
134 }
135 G4VEmModel* model = EmModel(0);
136 G4double emin = std::max(param->MinKinEnergy(),model->LowEnergyLimit());
137 G4double emax = std::min(param->MaxKinEnergy(),model->HighEnergyLimit());
138 model->SetPolarAngleLimit(theta);
139 model->SetLowEnergyLimit(emin);
140 model->SetHighEnergyLimit(emax);
141 AddEmModel(1, model);
142}
fEmIncreasing
@ fEmIncreasing
Definition: G4EmTableType.hh:57
G4double
double G4double
Definition: G4Types.hh:83
G4bool
bool G4bool
Definition: G4Types.hh:86
G4EmParameters::Instance
static G4EmParameters * Instance()
Definition: G4EmParameters.cc:70
G4EmParameters::MinKinEnergy
G4double MinKinEnergy() const
Definition: G4EmParameters.cc:574
G4EmParameters::MscThetaLimit
G4double MscThetaLimit() const
Definition: G4EmParameters.cc:775
G4EmParameters::MaxKinEnergy
G4double MaxKinEnergy() const
Definition: G4EmParameters.cc:593
G4EmParameters::FactorForAngleLimit
G4double FactorForAngleLimit() const
Definition: G4EmParameters.cc:757
G4ParticleDefinition::GetParticleType
const G4String & GetParticleType() const
Definition: G4ParticleDefinition.hh:130
G4ParticleDefinition::GetParticleName
const G4String & GetParticleName() const
Definition: G4ParticleDefinition.hh:109
G4VEmModel::SetPolarAngleLimit
void SetPolarAngleLimit(G4double)
Definition: G4VEmModel.hh:781
G4VEmModel::SetHighEnergyLimit
void SetHighEnergyLimit(G4double)
Definition: G4VEmModel.hh:746
G4VEmModel::LowEnergyLimit
G4double LowEnergyLimit() const
Definition: G4VEmModel.hh:641
G4VEmModel::HighEnergyLimit
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:634
G4VEmModel::SetLowEnergyLimit
void SetLowEnergyLimit(G4double)
Definition: G4VEmModel.hh:753
G4VEmProcess::EmModel
G4VEmModel * EmModel(size_t index=0) const
Definition: G4VEmProcess.hh:798
G4VEmProcess::AddEmModel
void AddEmModel(G4int, G4VEmModel *, const G4Region *region=nullptr)
Definition: G4VEmProcess.cc:125
G4VEmProcess::SetEmModel
void SetEmModel(G4VEmModel *, G4int index=0)
Definition: G4VEmProcess.cc:136
G4VProcess::SetVerboseLevel
void SetVerboseLevel(G4int value)
Definition: G4VProcess.hh:416
G4InuclParticleNames::name
const char * name(G4int ptype)
Definition: G4InuclParticleNames.hh:76

◆ IsApplicable()

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

Implements G4VEmProcess.

Definition at line 78 of file G4CoulombScattering.cc.

79{
80 return (p.GetPDGCharge() != 0.0);
81}
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:113

◆ MinPrimaryEnergy()

G4double G4CoulombScattering::MinPrimaryEnergy ( const G4ParticleDefinition *  part,
const G4Material *  mat 
)
finalprotectedvirtual

Reimplemented from G4VDiscreteProcess.

Definition at line 146 of file G4CoulombScattering.cc.

148{
149 // Pure Coulomb scattering
150 G4double emin = 0.0;
151
152 // Coulomb scattering combined with multiple or hadronic scattering
153 G4double theta = G4EmParameters::Instance()->MscThetaLimit();
154
155 if(0.0 < theta) {
156 G4double p2 = q2Max*mat->GetIonisation()->GetInvA23()/(1.0 - cos(theta));
157 G4double mass = part->GetPDGMass();
158 emin = sqrt(p2 + mass*mass) - mass;
159 }
160
161 return emin;
162}
G4IonisParamMat::GetInvA23
G4double GetInvA23() const
Definition: G4IonisParamMat.hh:145
G4Material::GetIonisation
G4IonisParamMat * GetIonisation() const
Definition: G4Material.hh:221

◆ operator=()

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

◆ ProcessDescription()

void G4CoulombScattering::ProcessDescription ( std::ostream &  out) const
overridevirtual

Reimplemented from G4VProcess.

Definition at line 180 of file G4CoulombScattering.cc.

181{
182 out <<
183 " Coulomb scattering. Simulation of elastic scattering\n" <<
184 " events individually. May be used in combination with multiple\n" <<
185 " scattering, where Coulomb scattering is used for hard (large angle)\n" <<
186 " collisions and multiple scattering for soft collisions.";
187 G4VEmProcess::ProcessDescription(out);
188}
G4VEmProcess::ProcessDescription
void ProcessDescription(std::ostream &outFile) const override
Definition: G4VEmProcess.cc:881

◆ StreamProcessInfo()

void G4CoulombScattering::StreamProcessInfo ( std::ostream &  outFile) const
overrideprotectedvirtual

Reimplemented from G4VEmProcess.

Definition at line 166 of file G4CoulombScattering.cc.

167{
168 G4double tetmin = G4EmParameters::Instance()->MscThetaLimit()/CLHEP::degree;
169 outFile << " ";
170 if(tetmin > 179.) { outFile << "ThetaMin(p)"; }
171 else { outFile << tetmin; }
172 outFile << " < Theta(degree) < 180";
173
174 if(q2Max < DBL_MAX) { outFile << "; pLimit(GeV^1)= " << sqrt(q2Max)/GeV; }
175 outFile << G4endl;
176}
G4endl
#define G4endl
Definition: G4ios.hh:57
DBL_MAX
#define DBL_MAX
Definition: templates.hh:62
The documentation for this class was generated from the following files: