150 if(posiKinEnergy == 0.0) {
152 G4double sint = sqrt((1. - cost)*(1. + cost));
163 pol.
set(-sinphi, cosphi, 0.0);
174 G4double sqgrate = sqrt(tau/tau2)*0.5;
180 G4double epsilqot = epsilmax/epsilmin;
189 greject = 1. - epsil + (2.*gam*epsil-1.)/(epsil*tau2*tau2);
191 }
while( greject < rndmEngine->
flat());
197 G4double cost = (epsil*tau2-1.)/(epsil*sqg2m1);
198 if(std::abs(cost) > 1.0) {
199 G4cout <<
"### G4eplusTo3GammaOKVIModel WARNING cost= " << cost
200 <<
" positron Ekin(MeV)= " << posiKinEnergy
201 <<
" gamma epsil= " << epsil
203 if(cost > 1.0) cost = 1.0;
206 G4double sint = sqrt((1.+cost)*(1.-cost));
214 G4double phot1Energy = epsil*TotalAvailableEnergy;
216 G4ThreeVector phot1Direction(sint*cos(phi), sint*sin(phi), cost);
217 phot1Direction.
rotateUz(posiDirection);
226 G4double phot2Energy =(1.-epsil)*TotalAvailableEnergy;
235 pol.
set(-sinphi, cosphi, 0.0);
237 cost = pol*phot2Direction;
238 pol -= cost*phot2Direction;
251 vdp->push_back(aGamma1);
252 vdp->push_back(aGamma2);
253 vdp->push_back(aGamma3);
virtual ~G4eplusTo3GammaOKVIModel()
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
const G4ThreeVector & GetMomentumDirection() const
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Hep3Vector & rotateUz(const Hep3Vector &)
G4ParticleChangeForGamma * fParticleChange
static constexpr double classic_electr_radius
G4double G4Log(G4double x)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) 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
G4ParticleDefinition * theGamma
G4double LowestTripletEnergy() const
virtual G4double ComputeCrossSectionPerAtom(const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0., G4double maxEnergy=DBL_MAX) final
G4eplusTo3GammaOKVIModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="eplus3ggOKVI")
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
static constexpr double pi
virtual G4double CrossSectionPerVolume(const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX) final
void ProposeTrackStatus(G4TrackStatus status)
G4double GetElectronDensity() const
static G4EmParameters * Instance()