Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
G4CrossSectionHandler Class Reference

#include <G4CrossSectionHandler.hh>

+ Inheritance diagram for G4CrossSectionHandler:

Public Member Functions

 G4CrossSectionHandler ()
 
 ~G4CrossSectionHandler ()
 
- 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 G4ElementSelectRandomElement (const G4MaterialCutsCouple *material, G4double e) const
 
G4int SelectRandomShell (G4int Z, G4double e) const
 
G4VEMDataSetBuildMeanFreePathForMaterials (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

virtual std::vector< G4VEMDataSet * > * BuildCrossSectionsForMaterials (const G4DataVector &energyVector, const G4DataVector *energyCuts=0)
 
- 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 G4VDataSetAlgorithmCreateInterpolation ()
 
const G4VDataSetAlgorithmGetInterpolation () const
 

Detailed Description

Definition at line 56 of file G4CrossSectionHandler.hh.

Constructor & Destructor Documentation

◆ G4CrossSectionHandler()

G4CrossSectionHandler::G4CrossSectionHandler ( )

Definition at line 63 of file G4CrossSectionHandler.cc.

64{ }

◆ ~G4CrossSectionHandler()

G4CrossSectionHandler::~G4CrossSectionHandler ( )

Definition at line 66 of file G4CrossSectionHandler.cc.

67{ }

Member Function Documentation

◆ BuildCrossSectionsForMaterials()

std::vector< G4VEMDataSet * > * G4CrossSectionHandler::BuildCrossSectionsForMaterials ( const G4DataVector energyVector,
const G4DataVector energyCuts = 0 
)
protectedvirtual

Implements G4VCrossSectionHandler.

Definition at line 70 of file G4CrossSectionHandler.cc.

72{
73 G4DataVector* energies;
74 G4DataVector* data;
75
76 G4DataVector* log_energies;
77 G4DataVector* log_data;
78
79 std::vector<G4VEMDataSet*>* matCrossSections = new std::vector<G4VEMDataSet*>;
80
81 const G4ProductionCutsTable* theCoupleTable=
83 size_t numOfCouples = theCoupleTable->GetTableSize();
84
85 size_t nOfBins = energyVector.size();
86 const G4VDataSetAlgorithm* interpolationAlgo = CreateInterpolation();
87
88 for (size_t mLocal=0; mLocal<numOfCouples; mLocal++)
89 {
90 const G4MaterialCutsCouple* couple = theCoupleTable->GetMaterialCutsCouple(mLocal);
91 const G4Material* material= couple->GetMaterial();
92 G4int nElements = material->GetNumberOfElements();
93 const G4ElementVector* elementVector = material->GetElementVector();
94 const G4double* nAtomsPerVolume = material->GetAtomicNumDensityVector();
95
96 G4VDataSetAlgorithm* algo = interpolationAlgo->Clone();
97
98 G4VEMDataSet* setForMat = new G4CompositeEMDataSet(algo,1.,1.);
99
100 for (G4int i=0; i<nElements; i++) {
101
102 G4int Z = (G4int) (*elementVector)[i]->GetZ();
103 G4double density = nAtomsPerVolume[i];
104
105 energies = new G4DataVector;
106 data = new G4DataVector;
107
108 log_energies = new G4DataVector;
109 log_data = new G4DataVector;
110
111
112 for (size_t bin=0; bin<nOfBins; bin++)
113 {
114 G4double e = energyVector[bin];
115 energies->push_back(e);
116 if (e==0.) e=1e-300;
117 log_energies->push_back(std::log10(e));
118 G4double cross = density*FindValue(Z,e);
119 data->push_back(cross);
120 if (cross==0.) cross=1e-300;
121 log_data->push_back(std::log10(cross));
122 }
123
124 G4VDataSetAlgorithm* algo1 = interpolationAlgo->Clone();
125
126// G4VEMDataSet* elSet = new G4EMDataSet(i,energies,data,algo1,1.,1.);
127
128 G4VEMDataSet* elSet = new G4EMDataSet(i,energies,data,log_energies,log_data,algo1,1.,1.);
129
130 setForMat->AddComponent(elSet);
131 }
132
133 matCrossSections->push_back(setForMat);
134 }
135 delete interpolationAlgo;
136 return matCrossSections;
137}
std::vector< G4Element * > G4ElementVector
double G4double
Definition: G4Types.hh:83
int G4int
Definition: G4Types.hh:85
const G4Material * GetMaterial() const
const G4ElementVector * GetElementVector() const
Definition: G4Material.hh:188
size_t GetNumberOfElements() const
Definition: G4Material.hh:184
const G4double * GetAtomicNumDensityVector() const
Definition: G4Material.hh:214
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
std::size_t GetTableSize() const
static G4ProductionCutsTable * GetProductionCutsTable()
G4double FindValue(G4int Z, G4double e) const
virtual G4VDataSetAlgorithm * CreateInterpolation()
virtual G4VDataSetAlgorithm * Clone() const =0
virtual void AddComponent(G4VEMDataSet *dataSet)=0

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