73 maxKineticEnergy = 1.2*
GeV;
132 nModels = models.size();
145 models.push_back(model);
147 ekinMin.push_back(elow);
148 if(thKineticEnergy > elow) thKineticEnergy = elow;
151 cumSum.push_back(0.0);
186 if (kinEnergy > thKineticEnergy) {
188 for(
G4int i=0; i<nModels; i++) {
207 G4cout <<
" e+ annihilation into hadrons active from "
208 << e1/
GeV <<
" GeV to " << e2/
GeV <<
" GeV"
220 G4cout <<
"### G4eeToHadronsMultiModel: The cross section for G4eeToHadronsMultiModel "
221 <<
" is increased by the Factor= " << csFactor <<
G4endl;
G4double LowEnergyLimit() const
G4double GetKineticEnergy() const
void SetLowEnergy(G4double val)
G4double HighEnergyLimit() const
void SetHighEnergy(G4double val)
void SetHighEnergyLimit(G4double)
virtual G4double ThresholdEnergy() const =0
G4GLOB_DLL std::ostream G4cout
G4double GetElectronDensity() const
const XML_Char XML_Content * model
G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
virtual G4double PeakEnergy() const =0
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4eeToHadronsMultiModel(G4int ver=0, const G4String &nam="eeToHadrons")
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
virtual ~G4eeToHadronsMultiModel()
G4double HighEnergy() const
void ProposeTrackStatus(G4TrackStatus status)
void SetLowEnergyLimit(G4double)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX)
void SetCrossSecFactor(G4double fac)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
G4ParticleChangeForGamma * GetParticleChangeForGamma()