G4PolarizedGammaConversionModel Class Reference

#include <G4PolarizedGammaConversionModel.hh>

Inheritance diagram for G4PolarizedGammaConversionModel:

Public Member Functions
	G4PolarizedGammaConversionModel (const G4ParticleDefinition *p=0, const G4String &nam="polConv")
virtual	~G4PolarizedGammaConversionModel ()
virtual void	Initialise (const G4ParticleDefinition *, const G4DataVector &)
virtual void	SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy)
const G4Element *	SelectedAtom ()
Protected Attributes
G4VPolarizedCrossSection *	crossSectionCalculator

---

Detailed Description

Definition at line 58 of file G4PolarizedGammaConversionModel.hh.

---

Constructor & Destructor Documentation

G4PolarizedGammaConversionModel::G4PolarizedGammaConversionModel	(	const G4ParticleDefinition *	p = 0,
		const G4String &	nam = "polConv"
	)

Definition at line 68 of file G4PolarizedGammaConversionModel.cc.

  : G4BetheHeitlerModel(pd,nam), crossSectionCalculator(0)
{
}

G4PolarizedGammaConversionModel::~G4PolarizedGammaConversionModel

(

)

[virtual]

Definition at line 76 of file G4PolarizedGammaConversionModel.cc.

References crossSectionCalculator.

{
  if (crossSectionCalculator) delete crossSectionCalculator;
}

---

Member Function Documentation

void G4PolarizedGammaConversionModel::Initialise	(	const G4ParticleDefinition *	,
		const G4DataVector &
	)			[virtual]

Reimplemented from G4BetheHeitlerModel.

Definition at line 83 of file G4PolarizedGammaConversionModel.cc.

References crossSectionCalculator, and G4BetheHeitlerModel::Initialise().

{
  G4BetheHeitlerModel::Initialise(pd,dv);
  if (!crossSectionCalculator)
    crossSectionCalculator = new G4PolarizedPairProductionCrossSection();
}

void G4PolarizedGammaConversionModel::SampleSecondaries	(	std::vector< G4DynamicParticle * > *	,
		const G4MaterialCutsCouple *	,
		const G4DynamicParticle *	,
		G4double	tmin,
		G4double	maxEnergy
	)			[virtual]

Reimplemented from G4BetheHeitlerModel.

Definition at line 93 of file G4PolarizedGammaConversionModel.cc.

References crossSectionCalculator, G4cout, G4VEmModel::GetCurrentElement(), G4PolarizationHelper::GetFrame(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), G4VPolarizedCrossSection::GetPol2(), G4VPolarizedCrossSection::GetPol3(), G4DynamicParticle::GetPolarization(), G4VPolarizedCrossSection::Initialize(), G4StokesVector::InvRotateAz(), G4StokesVector::p1(), G4StokesVector::p2(), G4StokesVector::p3(), G4StokesVector::RotateAz(), G4BetheHeitlerModel::SampleSecondaries(), G4VPolarizedCrossSection::SetMaterial(), G4StokesVector::SetPhoton(), and G4StokesVector::ZERO.

{
  G4BetheHeitlerModel::SampleSecondaries(vdp, couple, dp, tmin, maxEnergy);
 
  if(vdp && vdp->size()>0) {
    G4double gamEnergy0 = dp->GetKineticEnergy();
    G4double lepEnergy1 = (*vdp)[0]->GetKineticEnergy();
    G4double sintheta   = dp->GetMomentumDirection().cross((*vdp)[0]->GetMomentumDirection()).mag();
    if (sintheta>1.) sintheta=1.;

    G4StokesVector beamPol = dp->GetPolarization();
    beamPol.SetPhoton();

    // determine interaction plane
    G4ThreeVector  nInteractionFrame = 
      G4PolarizationHelper::GetFrame(dp->GetMomentumDirection(), 
                                     (*vdp)[0]->GetMomentumDirection());

    // transform polarization into interaction frame
    beamPol.InvRotateAz(nInteractionFrame,dp->GetMomentumDirection());

    // calulcate polarization transfer
    crossSectionCalculator->SetMaterial(GetCurrentElement()->GetN(), // number of nucleons
                                        GetCurrentElement()->GetZ(), 
                                        GetCurrentElement()->GetfCoulomb());
    crossSectionCalculator->Initialize(gamEnergy0, lepEnergy1, sintheta,
                                       beamPol, G4StokesVector::ZERO);

    // deterimine final state polarization
    G4StokesVector lep1Pol = crossSectionCalculator->GetPol2();
    lep1Pol.RotateAz(nInteractionFrame,(*vdp)[0]->GetMomentumDirection());
    (*vdp)[0]->SetPolarization(lep1Pol.p1(),
                               lep1Pol.p2(),
                               lep1Pol.p3());

    size_t num = vdp->size();
    if (num!=2) G4cout<<" WARNING "<<num<<" secondaries in polarized pairproduction not supported!\n"; 
    for (size_t i =1; i<num; ++i) {
      G4StokesVector lep2Pol = crossSectionCalculator->GetPol3();
      lep2Pol.RotateAz(nInteractionFrame,(*vdp)[i]->GetMomentumDirection());
      (*vdp)[i]->SetPolarization(lep2Pol.p1(),
                                 lep2Pol.p2(),
                                 lep2Pol.p3());

    }
  }
}

const G4Element * G4PolarizedGammaConversionModel::SelectedAtom

(

)

[inline]

Definition at line 87 of file G4PolarizedGammaConversionModel.hh.

References G4VEmModel::GetCurrentElement().

{
  return GetCurrentElement();
}

---

Field Documentation

G4VPolarizedCrossSection* G4PolarizedGammaConversionModel::crossSectionCalculator [protected]

Definition at line 79 of file G4PolarizedGammaConversionModel.hh.

Referenced by Initialise(), SampleSecondaries(), and ~G4PolarizedGammaConversionModel().

---

The documentation for this class was generated from the following files:

• G4PolarizedGammaConversionModel.hh
• G4PolarizedGammaConversionModel.cc

---