G4ePolarizedBremsstrahlungModel Class Reference

#include <G4ePolarizedBremsstrahlungModel.hh>

Inheritance diagram for G4ePolarizedBremsstrahlungModel:

Public Member Functions
	G4ePolarizedBremsstrahlungModel (const G4ParticleDefinition *p=0, const G4String &nam="PolBrem")
virtual	~G4ePolarizedBremsstrahlungModel ()
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 G4ePolarizedBremsstrahlungModel.hh.

---

Constructor & Destructor Documentation

G4ePolarizedBremsstrahlungModel::G4ePolarizedBremsstrahlungModel	(	const G4ParticleDefinition *	p = 0,
		const G4String &	nam = "PolBrem"
	)

Definition at line 56 of file G4ePolarizedBremsstrahlungModel.cc.

  : G4eBremsstrahlungModel(p,nam),
    crossSectionCalculator(0)
{
}

G4ePolarizedBremsstrahlungModel::~G4ePolarizedBremsstrahlungModel

(

)

[virtual]

Definition at line 65 of file G4ePolarizedBremsstrahlungModel.cc.

References crossSectionCalculator.

{
  if (crossSectionCalculator) delete crossSectionCalculator;
}

---

Member Function Documentation

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

Reimplemented from G4eBremsstrahlungModel.

Definition at line 72 of file G4ePolarizedBremsstrahlungModel.cc.

References crossSectionCalculator, and G4eBremsstrahlungModel::Initialise().

{
  G4eBremsstrahlungModel::Initialise(p,d);
  if (!crossSectionCalculator)
    crossSectionCalculator = new G4PolarizedBremsstrahlungCrossSection();
}

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

Reimplemented from G4eBremsstrahlungModel.

Definition at line 84 of file G4ePolarizedBremsstrahlungModel.cc.

References crossSectionCalculator, G4eBremsstrahlungModel::fParticleChange, G4cout, G4VEmModel::GetCurrentElement(), G4PolarizationHelper::GetFrame(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMomentumDirection(), G4VPolarizedCrossSection::GetPol2(), G4VPolarizedCrossSection::GetPol3(), G4DynamicParticle::GetPolarization(), G4ParticleChangeForLoss::GetProposedMomentumDirection(), G4VPolarizedCrossSection::Initialize(), G4StokesVector::InvRotateAz(), G4ParticleChangeForLoss::ProposePolarization(), G4StokesVector::RotateAz(), G4eBremsstrahlungModel::SampleSecondaries(), G4VPolarizedCrossSection::SetMaterial(), G4StokesVector::SetPhoton(), and G4StokesVector::ZERO.

{
  G4eBremsstrahlungModel::SampleSecondaries(vdp,couple,dp,tmin,maxEnergy);
  G4int num = vdp->size();

  if(num > 0) {
    G4double lepEnergy0 = dp->GetKineticEnergy();
    G4double gamEnergy1 = (*vdp)[0]->GetKineticEnergy();
    G4double sintheta   = dp->GetMomentumDirection().cross((*vdp)[0]->GetMomentumDirection()).mag();
    if (sintheta>1.) sintheta=1.;


    G4StokesVector beamPol = dp->GetPolarization();

    // determine interaction plane
    G4ThreeVector  nInteractionFrame = 
      G4PolarizationHelper::GetFrame(dp->GetMomentumDirection(), 
                 fParticleChange->GetProposedMomentumDirection());

    // 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(lepEnergy0, gamEnergy1, sintheta,
                                       beamPol, G4StokesVector::ZERO);

    // deterimine final state polarization
    G4StokesVector newBeamPol = crossSectionCalculator->GetPol2();
    newBeamPol.RotateAz(nInteractionFrame,
        fParticleChange->GetProposedMomentumDirection());
    fParticleChange->ProposePolarization(newBeamPol);

    if (num!=1) G4cout<<" WARNING "<<num<<" secondaries in polarized bremsstrahlung not supported!\n"; 
    for (G4int i=0; i<num; i++) {
      G4StokesVector photonPol = crossSectionCalculator->GetPol3();
      photonPol.SetPhoton();
      photonPol.RotateAz(nInteractionFrame,(*vdp)[i]->GetMomentumDirection());
      (*vdp)[i]->SetPolarization(photonPol.p1(),
                                 photonPol.p2(),
                                 photonPol.p3());
    }
  }
  return;
}

const G4Element * G4ePolarizedBremsstrahlungModel::SelectedAtom

(

)

[inline]

Definition at line 86 of file G4ePolarizedBremsstrahlungModel.hh.

References G4VEmModel::GetCurrentElement().

{
  return GetCurrentElement();
}

---

Field Documentation

G4VPolarizedCrossSection* G4ePolarizedBremsstrahlungModel::crossSectionCalculator [protected]

Definition at line 79 of file G4ePolarizedBremsstrahlungModel.hh.

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

---

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

• G4ePolarizedBremsstrahlungModel.hh
• G4ePolarizedBremsstrahlungModel.cc

---