#include <G4StokesVector.hh>

Inheritance diagram for G4StokesVector:

Collaboration 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.

   : G4ThreeVector(),isPhoton(false)
 {
 }

Here is the caller graph for this function:

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.

68 {

69 }

Member Function Documentation

void G4StokesVector::DiceP1 ( )

Definition at line 168 of file G4StokesVector.cc.

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

Here is the call graph for this function:

void G4StokesVector::DiceP2 ( )

Definition at line 176 of file G4StokesVector.cc.

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

Here is the call graph for this function:

void G4StokesVector::DiceP3 ( )

Definition at line 184 of file G4StokesVector.cc.

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

Here is the call graph for this function:

void G4StokesVector::DiceUniform ( )

Definition at line 158 of file G4StokesVector.cc.

 {
   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);
 }

Here is the call graph for this function:

void G4StokesVector::FlipP3 ( )

Definition at line 192 of file G4StokesVector.cc.

 {
   setZ(-z());
 }

Here is the call graph for this function:

G4double G4StokesVector::GetBeta ( )

Definition at line 151 of file G4StokesVector.cc.

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

Here is the call graph for this function:

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

Definition at line 108 of file G4StokesVector.cc.

 {
   // 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);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

G4bool G4StokesVector::IsZero ( ) const

inline

Definition at line 78 of file G4StokesVector.hh.

78 { return *this==ZERO; }

G4StokesVector::ZERO

static const G4StokesVector ZERO

Definition: G4StokesVector.hh:62

Here is the caller graph for this function:

G4double G4StokesVector::p1 ( ) const

inline

Definition at line 74 of file G4StokesVector.hh.

74 { return x(); }

CLHEP::Hep3Vector::x

double x() const

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4StokesVector::p2 ( ) const

inline

Definition at line 75 of file G4StokesVector.hh.

75 { return y(); }

CLHEP::Hep3Vector::y

double y() const

Here is the call graph for this function:

Here is the caller graph for this function:

G4double G4StokesVector::p3 ( ) const

inline

Definition at line 76 of file G4StokesVector.hh.

76 { return z(); }

CLHEP::Hep3Vector::z

double z() const

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector G4StokesVector::PolDiv ( const G4StokesVector & b )

Definition at line 204 of file G4StokesVector.cc.

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

Here is the call graph for this function:

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

Definition at line 197 of file G4StokesVector.cc.

 {
   // 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();
 }

Here is the call graph for this function:

G4ThreeVector G4StokesVector::PolSqr ( ) const

inline

Definition at line 81 of file G4StokesVector.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

G4ThreeVector G4StokesVector::PolSqrt ( ) const

inline

Definition at line 84 of file G4StokesVector.hh.

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

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 71 of file G4StokesVector.cc.

 {
   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);
 }

Here is the call graph for this function:

Here is the caller graph for this function:

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

Definition at line 132 of file G4StokesVector.cc.

 {
   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);
 }

Here is the call graph for this function:

void G4StokesVector::SetPhoton ( )

inline

Definition at line 92 of file G4StokesVector.hh.

92 { isPhoton=true; }

Here is the caller graph for this function:

G4double G4StokesVector::Transverse ( ) const

inline

Definition at line 79 of file G4StokesVector.hh.

79 { return perp(); }

CLHEP::Hep3Vector::perp

double perp() const

Here is the call graph for this function:

Member Data 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.

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

static

Definition at line 64 of file G4StokesVector.hh.

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

static

Definition at line 65 of file G4StokesVector.hh.

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

static

Definition at line 62 of file G4StokesVector.hh.

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

geant4.10.03.p01/source/processes/electromagnetic/polarisation/include/G4StokesVector.hh
geant4.10.03.p01/source/processes/electromagnetic/polarisation/src/G4StokesVector.cc

Public Member Functions

Static Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation