75 crossBornPerElectron(0),
131 G4cout <<
"G4eeToHadronsModel: Cross secsions per electron"
132 <<
" nbins= " << nbins
142 <<
" crossBorn(nb)= " << s2/
nanobarn
212 G4int np = newp->size();
213 for(
G4int j=0; j<np; ++j) {
228 newp->push_back(gamma);
230 G4cout <<
"G4eeToHadronsModel::SampleSecondaries: Ebalance(MeV)= "
231 << t/
MeV <<
" primary 4-momentum: " << inlv <<
G4endl;
246 G4double bt = 2.0*fine_structure_const*(LL - 1.0)/
pi;
248 G4double del = 1. + fine_structure_const*(1.5*LL +
pi*
pi/3. -2.)/
pi;
252 cs += s1*(del*
G4Exp(
G4Log(x1)*bt) - bt*(x1 - 0.25*x1*x1));
260 for(
G4int j=i-2; j>=0; --j) {
264 w1 = bt*(del*
G4Exp(
G4Log(x1)*btm1) - 1.0 + 0.5*x1);
265 cs += 0.5*(x1 - x2)*(w2*s2 + w1*s1);
284 G4double bt = 2.0*fine_structure_const*(LL - 1.)/
pi;
286 G4double del = 1. + fine_structure_const*(1.5*LL +
pi*
pi/3. -2.)/
pi;
329 const G4int iimax = 1000;
341 G4cout <<
"G4DynamicParticle* G4eeToHadronsModel:WARNING "
342 << f <<
" > " << grej <<
" majorant is`small!"
345 if(++ii >= iimax) {
break; }
virtual ~G4eeToHadronsModel()
G4double LowEnergy() const
G4double LowEnergyLimit() const
static const double nanobarn
G4double GetKineticEnergy() const
CLHEP::Hep3Vector G4ThreeVector
G4double HighEnergyLimit() const
G4eeToHadronsModel(G4Vee2hadrons *, G4int ver=0, const G4String &nam="eeToHadrons")
G4PhysicsVector * PhysicsVector() const
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX)
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy, G4double maxEnergy)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
size_t GetVectorLength() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kineticEnergy, G4double Z, G4double A, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
G4PhysicsVector * crossPerElectron
G4GLOB_DLL std::ostream G4cout
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, G4double, const G4ThreeVector &)=0
G4double GetElectronDensity() const
const G4ThreeVector & GetMomentumDirection() const
void PutValue(size_t index, G4double theValue)
G4PhysicsVector * crossBornPerElectron
G4double Energy(size_t index) const
G4double Value(G4double theEnergy, size_t &lastidx) const
G4LorentzVector Get4Momentum() const
virtual G4double PeakEnergy() const =0
G4double G4Log(G4double x)
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
void Set4Momentum(const G4LorentzVector &momentum)
void ComputeCMCrossSectionPerElectron()
T max(const T t1, const T t2)
brief Return the largest of the two arguments
G4double energy(const ThreeVector &p, const G4double m)
G4DynamicParticle * GenerateCMPhoton(G4double)
const G4double x[NPOINTSGL]
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static const G4int LL[nN]
G4double HighEnergy() const
G4ParticleDefinition * theGamma
virtual G4double ComputeCrossSectionPerElectron(const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
virtual G4double ComputeCrossSection(G4double) const =0
CLHEP::HepLorentzVector G4LorentzVector