90 std::vector<G4VEMDataSet*>*
94 std::vector<G4VEMDataSet*>* set =
new std::vector<G4VEMDataSet*>;
102 G4int nOfBins = energyVector.size();
108 for (
size_t mLocal=0; mLocal<numOfCouples; mLocal++) {
116 G4double tcut = (*energyCuts)[mLocal];
121 for (
G4int i=0; i<nElements; i++) {
131 G4double density = nAtomsPerVolume[i];
136 energies->push_back(e);
138 log_energies->push_back(std::log10(e));
146 value *= elemCs*density;
148 cs->push_back(value);
150 if (value==0.) value=1e-300;
151 log_cs->push_back(std::log10(value));
161 set->push_back(setForMat);
174 if(energy > cutEnergy)
virtual void AddComponent(G4VEMDataSet *dataSet)=0
G4double FindValue(G4int Z, G4double e) const
G4VDataSetAlgorithm * interp
virtual G4VDataSetAlgorithm * Clone() const =0
const G4ElementVector * GetElementVector() const
const G4double * GetVecNbOfAtomsPerVolume() const
virtual G4double Probability(G4int Z, G4double minKineticEnergy, G4double maxKineticEnergy, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const =0
~G4BremsstrahlungCrossSectionHandler()
const XML_Char int const XML_Char * value
size_t GetTableSize() const
G4BremsstrahlungCrossSectionHandler(const G4VEnergySpectrum *spectrum, G4VDataSetAlgorithm *interpolation)
G4double GetCrossSectionAboveThresholdForElement(G4double energy, G4double cutEnergy, G4int Z)
static G4ProductionCutsTable * GetProductionCutsTable()
std::vector< G4Element * > G4ElementVector
std::vector< G4VEMDataSet * > * BuildCrossSectionsForMaterials(const G4DataVector &energyVector, const G4DataVector *energyCuts)
const G4MaterialCutsCouple * GetMaterialCutsCouple(G4int i) const
const G4Material * GetMaterial() const
size_t GetNumberOfElements() const
const G4VEnergySpectrum * theBR