G4StokesVector Class Reference

#include <G4StokesVector.hh>

Inheritance diagram for G4StokesVector:

Public Member Functions
	G4StokesVector ()
	G4StokesVector (const G4ThreeVector &v)
virtual	~G4StokesVector ()
G4double	p1 () const
G4double	p2 () const
G4double	p3 () const
G4bool	IsZero () const
G4double	Transverse () const
G4ThreeVector	PolSqr () const
G4ThreeVector	PolSqrt () const
G4ThreeVector	PolError (const G4StokesVector &sum2, long n)
G4ThreeVector	PolDiv (const G4StokesVector &)
void	SetPhoton ()
void	RotateAz (G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
void	InvRotateAz (G4ThreeVector nInteractionFrame, G4ThreeVector particleDirection)
void	RotateAz (G4double cosphi, G4double sinphi)
G4double	GetBeta ()
void	DiceUniform ()
void	DiceP1 ()
void	DiceP2 ()
void	DiceP3 ()
void	FlipP3 ()
Public Member Functions inherited from CLHEP::Hep3Vector
	Hep3Vector ()
	Hep3Vector (double x)
	Hep3Vector (double x, double y)
	Hep3Vector (double x, double y, double z)
	Hep3Vector (const Hep3Vector &)
	~Hep3Vector ()
double	operator() (int) const
double	operator[] (int) const
double &	operator() (int)
double &	operator[] (int)
double	x () const
double	y () const
double	z () const
void	setX (double)
void	setY (double)
void	setZ (double)
void	set (double x, double y, double z)
double	phi () const
double	theta () const
double	cosTheta () const
double	cos2Theta () const
double	mag2 () const
double	mag () const
void	setPhi (double)
void	setTheta (double)
void	setMag (double)
double	perp2 () const
double	perp () const
void	setPerp (double)
void	setCylTheta (double)
double	perp2 (const Hep3Vector &) const
double	perp (const Hep3Vector &) const
Hep3Vector &	operator= (const Hep3Vector &)
bool	operator== (const Hep3Vector &) const
bool	operator!= (const Hep3Vector &) const
bool	isNear (const Hep3Vector &, double epsilon=tolerance) const
double	howNear (const Hep3Vector &v) const
double	deltaR (const Hep3Vector &v) const
Hep3Vector &	operator+= (const Hep3Vector &)
Hep3Vector &	operator-= (const Hep3Vector &)
Hep3Vector	operator- () const
Hep3Vector &	operator*= (double)
Hep3Vector &	operator/= (double)
Hep3Vector	unit () const
Hep3Vector	orthogonal () const
double	dot (const Hep3Vector &) const
Hep3Vector	cross (const Hep3Vector &) const
double	angle (const Hep3Vector &) const
double	pseudoRapidity () const
void	setEta (double p)
void	setCylEta (double p)
Hep3Vector &	rotateX (double)
Hep3Vector &	rotateY (double)
Hep3Vector &	rotateZ (double)
Hep3Vector &	rotateUz (const Hep3Vector &)
Hep3Vector &	rotate (double, const Hep3Vector &)
Hep3Vector &	operator*= (const HepRotation &)
Hep3Vector &	transform (const HepRotation &)
void	setRThetaPhi (double r, double theta, double phi)
void	setREtaPhi (double r, double eta, double phi)
void	setRhoPhiZ (double rho, double phi, double z)
void	setRhoPhiTheta (double rho, double phi, double theta)
void	setRhoPhiEta (double rho, double phi, double eta)
double	getX () const
double	getY () const
double	getZ () const
double	getR () const
double	getTheta () const
double	getPhi () const
double	r () const
double	rho () const
double	getRho () const
double	eta () const
double	getEta () const
void	setR (double s)
void	setRho (double s)
int	compare (const Hep3Vector &v) const
bool	operator> (const Hep3Vector &v) const
bool	operator< (const Hep3Vector &v) const
bool	operator>= (const Hep3Vector &v) const
bool	operator<= (const Hep3Vector &v) const
double	diff2 (const Hep3Vector &v) const
bool	isParallel (const Hep3Vector &v, double epsilon=tolerance) const
bool	isOrthogonal (const Hep3Vector &v, double epsilon=tolerance) const
double	howParallel (const Hep3Vector &v) const
double	howOrthogonal (const Hep3Vector &v) const
double	beta () const
double	gamma () const
double	coLinearRapidity () const
double	angle () const
double	theta (const Hep3Vector &v2) const
double	cosTheta (const Hep3Vector &v2) const
double	cos2Theta (const Hep3Vector &v2) const
Hep3Vector	project () const
Hep3Vector	project (const Hep3Vector &v2) const
Hep3Vector	perpPart () const
Hep3Vector	perpPart (const Hep3Vector &v2) const
double	rapidity () const
double	rapidity (const Hep3Vector &v2) const
double	eta (const Hep3Vector &v2) const
double	polarAngle (const Hep3Vector &v2) const
double	deltaPhi (const Hep3Vector &v2) const
double	azimAngle (const Hep3Vector &v2) const
double	polarAngle (const Hep3Vector &v2, const Hep3Vector &ref) const
double	azimAngle (const Hep3Vector &v2, const Hep3Vector &ref) const
Hep3Vector &	rotate (const Hep3Vector &axis, double delta)
Hep3Vector &	rotate (const HepAxisAngle &ax)
Hep3Vector &	rotate (const HepEulerAngles &e)
Hep3Vector &	rotate (double phi, double theta, double psi)

Static Public Attributes
static const G4StokesVector	ZERO =G4ThreeVector(0.,0.,0.)
static const G4StokesVector	P1 =G4ThreeVector(1.,0.,0.)
static const G4StokesVector	P2 =G4ThreeVector(0.,1.,0.)
static const G4StokesVector	P3 =G4ThreeVector(0.,0.,1.)
static const G4StokesVector	M1 =G4ThreeVector(-1.,0.,0.)
static const G4StokesVector	M2 =G4ThreeVector(0.,-1.,0.)
static const G4StokesVector	M3 =G4ThreeVector(0.,0.,-1.)

Additional Inherited Members
Public Types inherited from CLHEP::Hep3Vector
enum	{   X =0, Y =1, Z =2, NUM_COORDINATES =3,   SIZE =NUM_COORDINATES }
enum	{ ToleranceTicks = 100 }
Static Public Member Functions inherited from CLHEP::Hep3Vector
static double	setTolerance (double tol)
static double	getTolerance ()
Protected Member Functions inherited from CLHEP::Hep3Vector
void	setSpherical (double r, double theta, double phi)
void	setCylindrical (double r, double phi, double z)
double	negativeInfinity () const
Protected Attributes inherited from CLHEP::Hep3Vector
double	dx
double	dy
double	dz
Static Protected Attributes inherited from CLHEP::Hep3Vector
static DLL_API double	tolerance = Hep3Vector::ToleranceTicks * 2.22045e-16

Detailed Description

Definition at line 58 of file G4StokesVector.hh.

Constructor & Destructor Documentation

G4StokesVector::G4StokesVector

(

)

Definition at line 57 of file G4StokesVector.cc.

Referenced by PolError().

  : G4ThreeVector(),isPhoton(false)
{
}

G4StokesVector::G4StokesVector

(

const G4ThreeVector &

v

)

Definition at line 62 of file G4StokesVector.cc.

  : G4ThreeVector(v),isPhoton(false)
{
}

G4StokesVector::~G4StokesVector ( )

virtual

Definition at line 67 of file G4StokesVector.cc.

{
}

Member Function Documentation

void G4StokesVector::DiceP1

(

)

Definition at line 168 of file G4StokesVector.cc.

References G4UniformRand, CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), and CLHEP::Hep3Vector::setZ().

{
  if (G4UniformRand()>0.5) setX(1.);
  else setX(-1.);
  setY(0.);
  setZ(0.);
}

void G4StokesVector::DiceP2

(

)

Definition at line 176 of file G4StokesVector.cc.

References G4UniformRand, CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), and CLHEP::Hep3Vector::setZ().

{
  setX(0.);
  if (G4UniformRand()>0.5) setY(1.);
  else setY(-1.);
  setZ(0.);
}

void G4StokesVector::DiceP3

(

)

Definition at line 184 of file G4StokesVector.cc.

References G4UniformRand, CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), and CLHEP::Hep3Vector::setZ().

{
  setX(0.);
  setY(0.);
  if (G4UniformRand()>0.5) setZ(1.);
  else setZ(-1.);
}

void G4StokesVector::DiceUniform

(

)

Definition at line 158 of file G4StokesVector.cc.

References G4UniformRand, python.hepunit::pi, CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), and CLHEP::Hep3Vector::setZ().

{
  G4double costheta=2.*G4UniformRand()-1.;
  G4double sintheta=std::sqrt(1.-costheta*costheta);
  G4double aphi     =2.*pi*G4UniformRand();
  setX(std::sin(aphi)*sintheta);
  setY(std::cos(aphi)*sintheta);
  setZ(costheta);
}

void G4StokesVector::FlipP3

(

)

Definition at line 192 of file G4StokesVector.cc.

References CLHEP::Hep3Vector::setZ(), and CLHEP::Hep3Vector::z().

{
  setZ(-z());
}

G4double G4StokesVector::GetBeta

(

)

Definition at line 151 of file G4StokesVector.cc.

References CLHEP::Hep3Vector::getPhi().

{
  G4double bet=getPhi();
  if (isPhoton) { bet *= 0.5; }
  return bet;
}

void G4StokesVector::InvRotateAz	(	G4ThreeVector	nInteractionFrame,
		G4ThreeVector	particleDirection
	)

Definition at line 108 of file G4StokesVector.cc.

References CLHEP::Hep3Vector::cross(), G4cout, G4PolarizationHelper::GetParticleFrameY(), and RotateAz().

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

{
  // note if incomming particle is on z-axis, 
  // we might encounter some nummerical problems, since
  // nInteratonFrame and yParticleFrame are actually (almost) the same momentum
  // and the normalization is only good to 10^-12 !

  G4ThreeVector  yParticleFrame = 
    G4PolarizationHelper::GetParticleFrameY(particleDirection);
  G4double cosphi=yParticleFrame*nInteractionFrame;

  if (cosphi>1.+1.e-8 || cosphi<-1.-1.e-8) {
        G4cout<<" warning G4StokesVector::RotateAz  cosphi>1 or cosphi<-1\n";
  }
  if (cosphi>1) cosphi=1.;
  else if (cosphi<-1)cosphi=-1.;

  // check sign once more!
  G4double hel=(yParticleFrame.cross(nInteractionFrame)*particleDirection)>0?1.:-1.;
  G4double sinphi=hel*std::sqrt(std::fabs(1.-cosphi*cosphi));
  RotateAz(cosphi,-sinphi);
}

G4bool G4StokesVector::IsZero ( ) const

inline

Definition at line 78 of file G4StokesVector.hh.

References ZERO.

Referenced by G4ePolarizedIonisation::GetMeanFreePath(), G4PolarizedBhabhaCrossSection::Initialize(), G4PolarizedMollerCrossSection::Initialize(), G4ePolarizedIonisation::PostStepGetPhysicalInteractionLength(), G4PolarizedBhabhaCrossSection::XSection(), and G4PolarizedMollerCrossSection::XSection().

{ return *this==ZERO; } 

G4double G4StokesVector::p1 ( ) const

inline

Definition at line 74 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::x().

Referenced by arorie(), attrad(), berttp(), dtwopa(), energy(), enrg(), mashel(), prod5(), pyeig4(), pytbdy(), RotateAz(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4PolarizedPEEffectModel::SampleSecondaries(), G4PolarizedComptonModel::SampleSecondaries(), and G4PolarizedAnnihilationModel::SampleSecondaries().

{ return x(); }

G4double G4StokesVector::p2 ( ) const

inline

Definition at line 75 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::y().

Referenced by attrad(), berttp(), dtwopa(), energy(), enrg(), lepto(), luindf(), mashel(), prod5(), pyeig4(), pygano(), pygbeh(), pygdir(), pyggam(), pygvmd(), pyindf(), pytbdy(), pyvacu(), RotateAz(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4PolarizedPEEffectModel::SampleSecondaries(), G4PolarizedComptonModel::SampleSecondaries(), and G4PolarizedAnnihilationModel::SampleSecondaries().

{ return y(); }

G4double G4StokesVector::p3 ( ) const

inline

Definition at line 76 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::z().

Referenced by arorie(), attrad(), G4PolarizedComptonModel::ComputeCrossSectionPerAtom(), G4PolarizedCompton::GetMeanFreePath(), pinkla(), G4PolarizedCompton::PostStepGetPhysicalInteractionLength(), prod5(), pytbdy(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4PolarizedPEEffectModel::SampleSecondaries(), G4PolarizedComptonModel::SampleSecondaries(), G4PolarizedAnnihilationModel::SampleSecondaries(), and G4PolarizedComptonCrossSection::TotalXSection().

{ return z(); }

G4ThreeVector G4StokesVector::PolDiv

(

const G4StokesVector &

b

)

Definition at line 204 of file G4StokesVector.cc.

References CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

  {return G4ThreeVector(b.x()!=0. ? x()/b.x() : 11111.,
            b.y()!=0. ? y()/b.y() : 11111., 
            b.z()!=0. ? z()/b.z() : 11111.);}

G4ThreeVector G4StokesVector::PolError	(	const G4StokesVector &	sum2,
		long	n
	)

Definition at line 197 of file G4StokesVector.cc.

References G4StokesVector(), n, PolSqr(), and PolSqrt().

{
  // delta x = sqrt[ ( <x^2> - <x>^2 )/(n-1) ]
  G4StokesVector mean=(1./n)*(*this);
  return G4StokesVector((1./(n-1.)*((1./n)*sum2 - mean.PolSqr()))).PolSqrt();
}

G4ThreeVector G4StokesVector::PolSqr ( ) const

inline

Definition at line 81 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by PolError().

                                      { 
    return G4ThreeVector(x()*x(),y()*y(),z()*z()); 
  }

G4ThreeVector G4StokesVector::PolSqrt ( ) const

inline

Definition at line 84 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

Referenced by PolError().

                                       { 
    return G4ThreeVector(std::sqrt(x()),std::sqrt(y()),std::sqrt(z())); 
  }

void G4StokesVector::RotateAz	(	G4ThreeVector	nInteractionFrame,
		G4ThreeVector	particleDirection
	)

Definition at line 71 of file G4StokesVector.cc.

References CLHEP::Hep3Vector::cross(), G4cout, G4endl, G4PolarizationHelper::GetParticleFrameY(), and CLHEP::Hep3Vector::mag().

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

{
  G4ThreeVector  yParticleFrame = 
    G4PolarizationHelper::GetParticleFrameY(particleDirection);


  G4double cosphi=yParticleFrame*nInteractionFrame;
  if (cosphi>(1.+1.e-8) || cosphi<(-1.-1.e-8)) {
    G4cout<<" warning G4StokesVector::RotateAz  cosphi>1 or cosphi<-1\n"
      <<" cosphi="<<cosphi<<"\n"
      <<" zAxis="<<particleDirection<<" ("<<particleDirection.mag()<<")\n"
      <<" yAxis="<<yParticleFrame<<" ("<<yParticleFrame.mag()<<")\n"
      <<" nAxis="<<nInteractionFrame<<" ("
      <<nInteractionFrame.mag()<<")"<<G4endl;
  }
  if (cosphi>1.) cosphi=1.;
  else if (cosphi<-1.) cosphi=-1.;

//     G4cout<<" cosphi="<<cosphi<<"\n"
//    <<" zAxis="<<particleDirection<<" ("<<particleDirection.mag()<<")\n"
//    <<" yAxis="<<yParticleFrame<<" ("<<yParticleFrame.mag()<<","<<(yParticleFrame*particleDirection)<<")\n"
//    <<" nAxis="<<nInteractionFrame<<" ("
//    <<nInteractionFrame.mag()<<")"<<G4endl;

  //  G4double hel=sgn(cross(yParticleFrame*nInteractionFrame)*zInteractionFrame);
  // Why not particleDirection instead of zInteractionFrame ???!!!
  // -> is the same, since SYSIN is called with p1, and p2 as first parameter!
  G4double hel=(yParticleFrame.cross(nInteractionFrame)*particleDirection)>0?1.:-1.;

  G4double sinphi=hel*std::sqrt(1.-cosphi*cosphi);
  //  G4cout<<" sin2 + cos2 -1 = "<<(sinphi*sinphi+cosphi*cosphi-1)<<"\n";

  RotateAz(cosphi,sinphi);
}

void G4StokesVector::RotateAz	(	G4double	cosphi,
		G4double	sinphi
	)

Definition at line 132 of file G4StokesVector.cc.

References p1(), p2(), CLHEP::Hep3Vector::setX(), and CLHEP::Hep3Vector::setY().

{
  if (!isPhoton) {
    G4double xsi1=  cosphi*p1() + sinphi*p2();
    G4double xsi2= -sinphi*p1() + cosphi*p2();
    setX(xsi1);
    setY(xsi2);
    return;
  }

  G4double sin2phi=2.*cosphi*sinphi;
  G4double cos2phi=cosphi*cosphi-sinphi*sinphi;

  G4double xsi1=  cos2phi*p1() + sin2phi*p2();
  G4double xsi2= -sin2phi*p1() + cos2phi*p2();
  setX(xsi1);
  setY(xsi2);
}

void G4StokesVector::SetPhoton ( )

inline

Definition at line 92 of file G4StokesVector.hh.

Referenced by G4PolarizedPairProductionCrossSection::Initialize(), G4PolarizedBremsstrahlungCrossSection::Initialize(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedPEEffectModel::SampleSecondaries(), G4PolarizedComptonModel::SampleSecondaries(), and G4PolarizedAnnihilationModel::SampleSecondaries().

{ isPhoton=true; }

G4double G4StokesVector::Transverse ( ) const

inline

Definition at line 79 of file G4StokesVector.hh.

References CLHEP::Hep3Vector::perp().

{ return perp(); } 

Field Documentation

const G4StokesVector G4StokesVector::M1 =G4ThreeVector(-1.,0.,0.)

static

Definition at line 66 of file G4StokesVector.hh.

const G4StokesVector G4StokesVector::M2 =G4ThreeVector(0.,-1.,0.)

static

Definition at line 67 of file G4StokesVector.hh.

const G4StokesVector G4StokesVector::M3 =G4ThreeVector(0.,0.,-1.)

static

Definition at line 68 of file G4StokesVector.hh.

const G4StokesVector G4StokesVector::P1 =G4ThreeVector(1.,0.,0.)

static

Definition at line 63 of file G4StokesVector.hh.

Referenced by G4PolarizedAnnihilationModel::ComputeAsymmetriesPerElectron(), G4VPolarizedCrossSection::GetPol2(), and G4VPolarizedCrossSection::GetPol3().

const G4StokesVector G4StokesVector::P2 =G4ThreeVector(0.,1.,0.)

static

Definition at line 64 of file G4StokesVector.hh.

Referenced by G4PolarizedAnnihilationModel::ComputeAsymmetriesPerElectron(), G4VPolarizedCrossSection::GetPol2(), and G4VPolarizedCrossSection::GetPol3().

const G4StokesVector G4StokesVector::P3 =G4ThreeVector(0.,0.,1.)

static

Definition at line 65 of file G4StokesVector.hh.

Referenced by G4PolarizedAnnihilationModel::ComputeAsymmetriesPerElectron(), G4VPolarizedCrossSection::GetPol2(), and G4VPolarizedCrossSection::GetPol3().

const G4StokesVector G4StokesVector::ZERO =G4ThreeVector(0.,0.,0.)

static

Definition at line 62 of file G4StokesVector.hh.

Referenced by G4PolarizedAnnihilationModel::ComputeAsymmetriesPerElectron(), G4PolarizedMollerBhabhaModel::ComputeCrossSectionPerElectron(), G4VPolarizedCrossSection::GetPol2(), G4VPolarizedCrossSection::GetPol3(), IsZero(), G4PolarizedGammaConversionModel::SampleSecondaries(), G4ePolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedPEEffectModel::SampleSecondaries(), G4PolarizedMollerBhabhaModel::SampleSecondaries(), and G4PolarizedComptonCrossSection::XSection().

---

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

• G4StokesVector.hh
• G4StokesVector.cc