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(),
G4VPolarizedCrossSection * crossSectionCalculator
G4double GetKineticEnergy() const
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
void SetMaterial(G4double A, G4double Z, G4double coul)
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
virtual G4StokesVector GetPol2()
G4GLOB_DLL std::ostream G4cout
G4PolarizedGammaConversionModel(const G4ParticleDefinition *p=0, const G4String &nam="polConv")
const G4ThreeVector & GetMomentumDirection() const
static G4ThreeVector GetFrame(const G4ThreeVector &, const G4ThreeVector &)
void InvRotateAz(G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
virtual void SampleSecondaries(std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
const G4ThreeVector & GetPolarization() const
virtual G4StokesVector GetPol3()
virtual void Initialise(const G4ParticleDefinition *, const G4DataVector &)
Hep3Vector cross(const Hep3Vector &) const
virtual ~G4PolarizedGammaConversionModel()
static const G4StokesVector ZERO
void RotateAz(G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
virtual void Initialize(G4double, G4double, G4double, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0)
const G4Element * GetCurrentElement() const