G4PolarizationHelper Class Reference

#include <G4PolarizationHelper.hh>

Static Public Member Functions
static G4ThreeVector	GetFrame (const G4ThreeVector &, const G4ThreeVector &)
static G4ThreeVector	GetParticleFrameX (const G4ThreeVector &)
static G4ThreeVector	GetParticleFrameY (const G4ThreeVector &)
static G4ThreeVector	GetRandomFrame (const G4ThreeVector &)
static G4ThreeVector	GetSpinInPRF (const G4ThreeVector &uZ, const G4ThreeVector &spin)
static void	TestPolarizationTransformations ()
static void	TestInteractionFrame ()

Detailed Description

Definition at line 48 of file G4PolarizationHelper.hh.

Member Function Documentation

GetFrame()

G4ThreeVector G4PolarizationHelper::GetFrame ( const G4ThreeVector &  mom1, const G4ThreeVector &  mom2  )

static

Definition at line 47 of file G4PolarizationHelper.cc.

{
  G4ThreeVector normal = (mom1.cross(mom2)).unit();
  return normal;
  //  return 1./normal.mag()*normal;
}

Referenced by G4PolarizedAnnihilationModel::SampleSecondaries(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedComptonModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4PolarizedMollerBhabhaModel::SampleSecondaries(), and G4PolarizedPEEffectModel::SampleSecondaries().

GetParticleFrameX()

G4ThreeVector G4PolarizationHelper::GetParticleFrameX ( const G4ThreeVector &  uZ )

static

Definition at line 66 of file G4PolarizationHelper.cc.

{
  // compare also G4ThreeVector::rotateUz()
 
  if (uZ.x()==0. && uZ.y()==0.) {
    if (uZ.z()>=0.) return G4ThreeVector(1.,0.,0.);
    return G4ThreeVector(-1.,0.,0.);
  }
 
  G4double perp    = std::sqrt(sqr(uZ.x())+sqr(uZ.y()));
  G4double invPerp = uZ.z()/perp;
  return G4ThreeVector(uZ.x()*invPerp,uZ.y()*invPerp,-perp);
}

Referenced by GetRandomFrame().

GetParticleFrameY()

G4ThreeVector G4PolarizationHelper::GetParticleFrameY ( const G4ThreeVector &  uZ )

static

Definition at line 54 of file G4PolarizationHelper.cc.

{
  // compare also G4ThreeVector::rotateUz()
 
  if (uZ.x()==0. && uZ.y()==0.) {
    return G4ThreeVector(0.,1.,0.);
  }
 
  G4double invPerp = 1./std::sqrt(sqr(uZ.x())+sqr(uZ.y()));
  return G4ThreeVector(-uZ.y()*invPerp,uZ.x()*invPerp,0);
}

Referenced by GetRandomFrame(), G4StokesVector::InvRotateAz(), G4StokesVector::RotateAz(), and TestPolarizationTransformations().

GetRandomFrame()

G4ThreeVector G4PolarizationHelper::GetRandomFrame ( const G4ThreeVector &  mom1 )

static

Definition at line 80 of file G4PolarizationHelper.cc.

{
  G4double phi     =2.*pi*G4UniformRand();
  G4ThreeVector normal = std::cos(phi)*GetParticleFrameX(mom1) 
    + std::sin(phi)*G4PolarizationHelper::GetParticleFrameY(mom1);
  return normal;
}

Referenced by G4PolarizedPEEffectModel::SampleSecondaries().

GetSpinInPRF()

G4ThreeVector G4PolarizationHelper::GetSpinInPRF ( const G4ThreeVector &  uZ, const G4ThreeVector &  spin  )

static

Definition at line 89 of file G4PolarizationHelper.cc.

{
  // compare also G4ThreeVector::rotateUz()
 
  if (uZ.x()==0. && uZ.y()==0.) {
    if (uZ.z()>=0.) return spin;
    return G4ThreeVector(-spin.x(),spin.y(),-spin.z());
  }
 
  G4double perp    = std::sqrt(sqr(uZ.x())+sqr(uZ.y()));
  G4double invPerp = 1./perp;
 
  G4ThreeVector uX(uZ.x()*uZ.z()*invPerp,uZ.y()*uZ.z()*invPerp,-perp);
  G4ThreeVector uY(-uZ.y()*invPerp,uZ.x()*invPerp,0); 
  
  return G4ThreeVector(spin*uX,spin*uY,spin*uZ);
}

TestInteractionFrame()

void G4PolarizationHelper::TestInteractionFrame ( )

static

Definition at line 143 of file G4PolarizationHelper.cc.

{
  // check transformation procedure for polarisation transfer 
  // calculation in scattering processes
  //  a) transfer target polarisation in beam particle reference frame (PRF)
  //  b) calc correct asymmetry w.r.t. scattering plane
  //  c) determine incomming polarisation in interaction frame (IF)
  //  d) transfer outgoing polarisation from IF to PRF
  G4cout<<"========================================\n\n";
 
  G4double theta=0.;
 
  G4ThreeVector dir0=G4ThreeVector(0.,0.,1.);
  G4ThreeVector dir2=G4ThreeVector(std::sin(theta),0.,std::cos(theta));
 
  G4StokesVector pol0=G4ThreeVector(0.,0.,1.);
  G4StokesVector pol1=G4ThreeVector(0.,0.,1.); 
 
  pol1.rotateUz(dir0);
 
  G4cout<<"========================================\n\n";
 
 
}

Referenced by G4PolarizationMessenger::SetNewValue().

TestPolarizationTransformations()

void G4PolarizationHelper::TestPolarizationTransformations ( )

static

Definition at line 107 of file G4PolarizationHelper.cc.

{
  G4double theta=0.;
  G4cout<<"========================================\n\n";
  for (G4int i=0; i<=10; ++i) {
    theta=pi*i/10.;
    G4ThreeVector zAxis = G4ThreeVector(std::sin(theta),0.,std::cos(theta));
    if (i==5) zAxis = G4ThreeVector(1.,0.,0.);
    if (i==10) zAxis = G4ThreeVector(0.,0.,-1.);
    G4ThreeVector yAxis = GetParticleFrameY(zAxis);
 
    G4cout<<zAxis<<" "<<zAxis.mag()<<"\n";
    G4cout<<yAxis<<" "<<yAxis.mag()<<"\n";
    G4ThreeVector xAxis = yAxis.cross(zAxis);
    G4cout<<xAxis<<" "<<xAxis.mag()<<"\n\n";
  }
 
  G4cout<<"========================================\n\n";
 
  for (G4int i=0; i<=10; ++i) {
    theta=pi*i/10.;
    G4ThreeVector zAxis = G4ThreeVector(0.,std::sin(theta),std::cos(theta));
    if (i==5) zAxis = G4ThreeVector(0.,1.,0.);
    if (i==10) zAxis = G4ThreeVector(0.,0.,-1.);
    G4ThreeVector yAxis = GetParticleFrameY(zAxis);
 
    G4cout<<zAxis<<" "<<zAxis.mag()<<"\n";
    G4cout<<yAxis<<" "<<yAxis.mag()<<"\n";
    G4ThreeVector xAxis = yAxis.cross(zAxis);
    G4cout<<xAxis<<" "<<xAxis.mag()<<"\n\n";
 
    G4cout<<"spat : "<<xAxis*yAxis.cross(zAxis)<<"\n\n";
  }
  G4cout<<"========================================\n\n";
}

Referenced by G4PolarizationMessenger::SetNewValue().

---

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

• G4PolarizationHelper.hh
• G4PolarizationHelper.cc