75 crossBornPerElectron(0),
131 G4cout <<
"G4eeToHadronsModel: Cross sections per electron"
132 <<
" nbins= " << nbins
142 <<
" crossBorn(nb)= " << s2/
nanobarn
207 model->SampleSecondaries(newp, mass, dir);
208 G4int np = newp->size();
209 for(
G4int j=0; j<np; ++j) {
224 newp->push_back(gamma);
225 if(std::abs(t) >
MeV) {
226 G4cout <<
"G4eeToHadronsModel::SampleSecondaries: Ebalance(MeV)= "
227 << t/
MeV <<
" primary 4-momentum: " << inlv <<
G4endl;
248 cs += s1*(del*
G4Exp(
G4Log(x1)*bt) - bt*(x1 - 0.25*x1*
x1));
256 for(
G4int j=i-2; j>=0; --j) {
261 cs += 0.5*(x1 -
x2)*(w2*s2 + w1*s1);
325 const G4int iimax = 1000;
337 G4cout <<
"G4DynamicParticle* G4eeToHadronsModel:WARNING "
338 << f <<
" > " << grej <<
" majorant is`small!"
341 if(++ii >= iimax) {
break; }
G4double Energy(size_t index) const
void set(double x, double y, double z)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
T max(const T t1, const T t2)
brief Return the largest of the two arguments
void ComputeCMCrossSectionPerElectron()
static constexpr double MeV
static const G4int LL[nN]
G4eeToHadronsModel(G4Vee2hadrons *, G4int ver=0, const G4String &nam="eeToHadrons")
Float_t x1[n_points_granero]
const G4ThreeVector & GetMomentumDirection() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX) override
G4double Value(G4double theEnergy, size_t &lastidx) const
static constexpr double nanobarn
G4double G4Log(G4double x)
G4double LowEnergyLimit() const
G4DynamicParticle * GenerateCMPhoton(G4double)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX) override
G4ParticleDefinition * theGamma
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy) override
static constexpr double electron_mass_c2
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
void Set4Momentum(const G4LorentzVector &momentum)
G4LorentzVector Get4Momentum() const
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
G4PhysicsVector * crossBornPerElectron
Hep3Vector boostVector() const
static constexpr double pi
CLHEP::HepLorentzVector G4LorentzVector
static constexpr double fine_structure_const
Float_t x2[n_points_geant4]
virtual ~G4eeToHadronsModel()
G4PhysicsVector * crossPerElectron
G4double HighEnergyLimit() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
const XML_Char XML_Content * model
G4double GetElectronDensity() const
void PutValue(size_t index, G4double theValue)
HepLorentzVector & boost(double, double, double)
size_t GetVectorLength() const
T min(const T t1, const T t2)
brief Return the smallest of the two arguments