101 for(
G4int i=0; i<= nbins; ++i) {
186 if(posiKinEnergy == 0.0) {
188 G4double sint = sqrt((1. - cost)*(1. + cost));
199 pol.
set(-sinphi, cosphi, 0.0);
210 G4double sqgrate = sqrt(tau/tau2)*0.5;
216 G4double epsilqot = epsilmax/epsilmin;
225 greject = 1. - epsil + (2.*gam*epsil-1.)/(epsil*tau2*tau2);
227 }
while( greject < rndmEngine->
flat());
233 G4double cost = (epsil*tau2-1.)/(epsil*sqg2m1);
234 if(std::abs(cost) > 1.0) {
235 G4cout <<
"### G4eplusTo2GammaOKVIModel WARNING cost= " << cost
236 <<
" positron Ekin(MeV)= " << posiKinEnergy
237 <<
" gamma epsil= " << epsil
239 if(cost > 1.0) cost = 1.0;
242 G4double sint = sqrt((1.+cost)*(1.-cost));
250 G4double phot1Energy = epsil*TotalAvailableEnergy;
252 G4ThreeVector phot1Direction(sint*cos(phi), sint*sin(phi), cost);
253 phot1Direction.
rotateUz(posiDirection);
262 G4double phot2Energy =(1.-epsil)*TotalAvailableEnergy;
271 pol.
set(-sinphi, cosphi, 0.0);
273 cost = pol*phot2Direction;
274 pol -= cost*phot2Direction;
287 vdp->push_back(aGamma1);
288 vdp->push_back(aGamma2);
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 SetPolarization(G4double polX, G4double polY, G4double polZ)
static constexpr double MeV
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX) final
const G4ThreeVector & GetMomentumDirection() const
G4double ComputeCrossSectionPerElectron(G4double kinEnergy)
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Hep3Vector & rotateUz(const Hep3Vector &)
static const G4double emax
G4double Value(G4double theEnergy, size_t &lastidx) const
G4ParticleChangeForGamma * fParticleChange
static constexpr double classic_electr_radius
G4double G4Log(G4double x)
G4eplusTo2GammaOKVIModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="eplus2ggOKVI")
G4eplusTo3GammaOKVIModel * f3GModel
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) final
static constexpr double TeV
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0., G4double maxEnergy=DBL_MAX) final
static constexpr double electron_mass_c2
G4double ComputeCrossSectionPerElectron(G4double kinEnergy)
static constexpr double twopi
static constexpr double eV
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX) final
static G4PhysicsVector * f3GProbability
G4double LowestTripletEnergy() const
virtual ~G4eplusTo2GammaOKVIModel()
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
static G4PhysicsVector * fCrossSection
static constexpr double pi
G4ParticleDefinition * theGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
void ProposeTrackStatus(G4TrackStatus status)
void SetTripletModel(G4VEmModel *)
G4double GetElectronDensity() const
static G4EmParameters * Instance()
void PutValue(size_t index, G4double theValue)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double maxEnergy=DBL_MAX) final