Geant4 10.7.0
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
G4eIonisationCrossSectionHandler Class Reference

#include <G4eIonisationCrossSectionHandler.hh>

+ Inheritance diagram for G4eIonisationCrossSectionHandler:

Public Member Functions

 G4eIonisationCrossSectionHandler (const G4VEnergySpectrum *spec, G4VDataSetAlgorithm *alg, G4double emin, G4double emax, G4int nbin)
 
 ~G4eIonisationCrossSectionHandler ()
 
G4double GetCrossSectionAboveThresholdForElement (G4double energy, G4double cutEnergy, G4int Z)
 
- Public Member Functions inherited from G4VCrossSectionHandler
 G4VCrossSectionHandler ()
 
 G4VCrossSectionHandler (G4VDataSetAlgorithm *interpolation, G4double minE=250 *CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int nBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)
 
virtual ~G4VCrossSectionHandler ()
 
void Initialise (G4VDataSetAlgorithm *interpolation=0, G4double minE=250 *CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int numberOfBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)
 
G4int SelectRandomAtom (const G4MaterialCutsCouple *couple, G4double e) const
 
const G4Element * SelectRandomElement (const G4MaterialCutsCouple *material, G4double e) const
 
G4int SelectRandomShell (G4int Z, G4double e) const
 
G4VEMDataSet * BuildMeanFreePathForMaterials (const G4DataVector *energyCuts=0)
 
G4double FindValue (G4int Z, G4double e) const
 
G4double FindValue (G4int Z, G4double e, G4int shellIndex) const
 
G4double ValueForMaterial (const G4Material *material, G4double e) const
 
void LoadData (const G4String &dataFile)
 
void LoadNonLogData (const G4String &dataFile)
 
void LoadShellData (const G4String &dataFile)
 
void PrintData () const
 
void Clear ()
 

Protected Member Functions

std::vector< G4VEMDataSet * > * BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts)
 
- Protected Member Functions inherited from G4VCrossSectionHandler
G4int NumberOfComponents (G4int Z) const
 
void ActiveElements ()
 
virtual std::vector< G4VEMDataSet * > * BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts=0)=0
 
virtual G4VDataSetAlgorithm * CreateInterpolation ()
 
const G4VDataSetAlgorithm * GetInterpolation () const
 

Detailed Description

Definition at line 62 of file G4eIonisationCrossSectionHandler.hh.

Constructor & Destructor Documentation

◆ G4eIonisationCrossSectionHandler()

G4eIonisationCrossSectionHandler::G4eIonisationCrossSectionHandler ( const G4VEnergySpectrum *  spec,
G4VDataSetAlgorithm *  alg,
G4double  emin,
G4double  emax,
G4int  nbin 
)

Definition at line 72 of file G4eIonisationCrossSectionHandler.cc.

75 : G4VCrossSectionHandler(),
76 theParam(spec),verbose(0)
77{
78 G4VCrossSectionHandler::Initialise(alg, emin, emax, nbin);
79 interp = new G4LinLogLogInterpolation();
80}
G4VCrossSectionHandler::Initialise
void Initialise(G4VDataSetAlgorithm *interpolation=0, G4double minE=250 *CLHEP::eV, G4double maxE=100 *CLHEP::GeV, G4int numberOfBins=200, G4double unitE=CLHEP::MeV, G4double unitData=CLHEP::barn, G4int minZ=1, G4int maxZ=99)
Definition: G4VCrossSectionHandler.cc:139

◆ ~G4eIonisationCrossSectionHandler()

G4eIonisationCrossSectionHandler::~G4eIonisationCrossSectionHandler ( )

Definition at line 83 of file G4eIonisationCrossSectionHandler.cc.

84{
85 delete interp;
86}

Member Function Documentation

◆ BuildCrossSectionsForMaterials()

std::vector< G4VEMDataSet * > * G4eIonisationCrossSectionHandler::BuildCrossSectionsForMaterials ( const G4DataVector &  energyVector,
const G4DataVector *  energyCuts 
)
protectedvirtual

Implements G4VCrossSectionHandler.

Definition at line 89 of file G4eIonisationCrossSectionHandler.cc.

92{
93 std::vector<G4VEMDataSet*>* set = new std::vector<G4VEMDataSet*>;
94
95 G4DataVector* energies;
96 G4DataVector* cs;
97
98 G4DataVector* log_energies;
99 G4DataVector* log_cs;
100
101 G4int nOfBins = energyVector.size();
102
103 const G4ProductionCutsTable* theCoupleTable=
104 G4ProductionCutsTable::GetProductionCutsTable();
105 size_t numOfCouples = theCoupleTable->GetTableSize();
106
107 for (size_t mLocal=0; mLocal<numOfCouples; mLocal++) {
108
109 const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(mLocal);
110 const G4Material* material= couple->GetMaterial();
111 const G4ElementVector* elementVector = material->GetElementVector();
112 const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
113 G4int nElements = material->GetNumberOfElements();
114
115 if(verbose > 0)
116 {
117 G4cout << "eIonisation CS for " << mLocal << "th material "
118 << material->GetName()
119 << " eEl= " << nElements << G4endl;
120 }
121
122 G4double tcut = (*energyCuts)[mLocal];
123
124 G4VDataSetAlgorithm* algo = interp->Clone();
125 G4VEMDataSet* setForMat = new G4CompositeEMDataSet(algo,1.,1.);
126
127 for (G4int i=0; i<nElements; i++) {
128
129 G4int Z = (G4int) (*elementVector)[i]->GetZ();
130 G4int nShells = NumberOfComponents(Z);
131
132 energies = new G4DataVector;
133 cs = new G4DataVector;
134
135 log_energies = new G4DataVector;
136 log_cs = new G4DataVector;
137
138 G4double density = nAtomsPerVolume[i];
139
140 for (G4int bin=0; bin<nOfBins; bin++) {
141
142 G4double e = energyVector[bin];
143 energies->push_back(e);
144 log_energies->push_back(std::log10(e));
145 G4double value = 0.0;
146 G4double log_value = -300;
147
148 if(e > tcut) {
149 for (G4int n=0; n<nShells; n++) {
150 G4double cross = FindValue(Z, e, n);
151 G4double p = theParam->Probability(Z, tcut, e, e, n);
152 value += cross * p * density;
153
154 if(verbose>0 && mLocal == 0 && e>=1. && e<=0.)
155 {
156 G4cout << "G4eIonCrossSH: e(MeV)= " << e/MeV
157 << " n= " << n
158 << " cross= " << cross
159 << " p= " << p
160 << " value= " << value
161 << " tcut(MeV)= " << tcut/MeV
162 << " rho= " << density
163 << " Z= " << Z
164 << G4endl;
165 }
166
167 }
168 if (value == 0.) value = 1e-300;
169 log_value = std::log10(value);
170 }
171 cs->push_back(value);
172 log_cs->push_back(log_value);
173 }
174 G4VDataSetAlgorithm* algoLocal = interp->Clone();
175
176 //G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,algoLocal,1.,1.);
177
178 G4VEMDataSet* elSet = new G4EMDataSet(i,energies,cs,log_energies,log_cs,algoLocal,1.,1.);
179
180 setForMat->AddComponent(elSet);
181 }
182 set->push_back(setForMat);
183 }
184
185 return set;
186}
G4ElementVector
std::vector< G4Element * > G4ElementVector
Definition: G4ElementVector.hh:43
G4double
double G4double
Definition: G4Types.hh:83
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:166
G4Material::GetElementVector
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:188
G4Material::GetNumberOfElements
size_t GetNumberOfElements() const
Definition: G4Material.hh:184
G4Material::GetAtomicNumDensityVector
const G4double * GetAtomicNumDensityVector() const
Definition: G4Material.hh:214
G4Material::GetName
const G4String & GetName() const
Definition: G4Material.hh:175
G4ProductionCutsTable::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
Definition: G4ProductionCutsTable.hh:246
G4ProductionCutsTable::GetTableSize
std::size_t GetTableSize() const
Definition: G4ProductionCutsTable.hh:239
G4ProductionCutsTable::GetProductionCutsTable
static G4ProductionCutsTable * GetProductionCutsTable()
Definition: G4ProductionCutsTable.cc:58
G4VCrossSectionHandler::FindValue
G4double FindValue(G4int Z, G4double e) const
Definition: G4VCrossSectionHandler.cc:378
G4VCrossSectionHandler::NumberOfComponents
G4int NumberOfComponents(G4int Z) const
Definition: G4VCrossSectionHandler.cc:727
G4VDataSetAlgorithm::Clone
virtual G4VDataSetAlgorithm * Clone() const =0
G4VEMDataSet::AddComponent
virtual void AddComponent(G4VEMDataSet *dataSet)=0
G4VEnergySpectrum::Probability
virtual G4double Probability(G4int Z, G4double minKineticEnergy, G4double maxKineticEnergy, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const =0

◆ GetCrossSectionAboveThresholdForElement()

G4double G4eIonisationCrossSectionHandler::GetCrossSectionAboveThresholdForElement ( G4double  energy,
G4double  cutEnergy,
G4int  Z 
)

Definition at line 188 of file G4eIonisationCrossSectionHandler.cc.

191{
192 G4int nShells = NumberOfComponents(Z);
193 G4double value = 0.;
194 if(energy > cutEnergy)
195 {
196 for (G4int n=0; n<nShells; n++) {
197 G4double cross = FindValue(Z, energy, n);
198 G4double p = theParam->Probability(Z, cutEnergy, energy, energy, n);
199 value += cross * p;
200 }
201 }
202 return value;
203}

Referenced by G4LivermoreIonisationModel::ComputeCrossSectionPerAtom().

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