101 if(vdp && vdp->size()>0) {
103 G4double lepEnergy1 = (*vdp)[0]->GetKineticEnergy();
105 if (sintheta>1.) sintheta=1.;
113 (*vdp)[0]->GetMomentumDirection());
127 lep1Pol.
RotateAz(nInteractionFrame,(*vdp)[0]->GetMomentumDirection());
128 (*vdp)[0]->SetPolarization(lep1Pol.
p1(),
132 size_t num = vdp->size();
133 if (num!=2)
G4cout<<
" WARNING "<<num<<
" secondaries in polarized pairproduction not supported!\n";
134 for (
size_t i =1; i<num; ++i) {
136 lep2Pol.
RotateAz(nInteractionFrame,(*vdp)[i]->GetMomentumDirection());
137 (*vdp)[i]->SetPolarization(lep2Pol.
p1(),
static G4ThreeVector GetFrame(const G4ThreeVector &, const G4ThreeVector &)
virtual G4StokesVector GetPol2()
void SetMaterial(G4double A, G4double Z, G4double coul)
G4VPolarizedCrossSection * crossSectionCalculator
const G4ThreeVector & GetMomentumDirection() const
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
G4PolarizedGammaConversionModel(const G4ParticleDefinition *p=nullptr, const G4String &nam="polConv")
virtual G4StokesVector GetPol3()
void InvRotateAz(G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
void RotateAz(G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
Hep3Vector cross(const Hep3Vector &) const
static const G4StokesVector ZERO
virtual ~G4PolarizedGammaConversionModel()
const G4ThreeVector & GetPolarization() const
G4double GetKineticEnergy() const
G4GLOB_DLL std::ostream G4cout
const G4Element * GetCurrentElement() const
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &) override
virtual void Initialize(G4double, G4double, G4double, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0)