G4ErrorSurfaceTrajState Class Reference

#include <G4ErrorSurfaceTrajState.hh>

Inheritance diagram for G4ErrorSurfaceTrajState:

Public Member Functions
	G4ErrorSurfaceTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4Plane3D &plane, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
	G4ErrorSurfaceTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
	G4ErrorSurfaceTrajState (G4ErrorFreeTrajState &tpSC, const G4Plane3D &plane)
	G4ErrorSurfaceTrajState (G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW, G4ErrorMatrix &transfM)
	~G4ErrorSurfaceTrajState ()
G4ErrorMatrix	BuildErrorMatrix (G4ErrorFreeTrajState &tpSC, const G4Vector3D &vecV, const G4Vector3D &vecW)
virtual void	Dump (std::ostream &out=G4cout) const
G4ErrorSurfaceTrajParam	GetParameters () const
void	SetParameters (const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW)
void	SetParameters (const G4Point3D &pos, const G4Vector3D &mom, const G4Plane3D &plane)
G4Vector3D	GetVectorV () const
G4Vector3D	GetVectorW () const
virtual void	SetPosition (const G4Point3D pos)
virtual void	SetMomentum (const G4Vector3D &mom)
Public Member Functions inherited from G4ErrorTrajState
	G4ErrorTrajState ()
	G4ErrorTrajState (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom, const G4ErrorTrajErr &errmat=G4ErrorTrajErr(5, 0))
	G4ErrorTrajState (const G4ErrorTrajState &)
	G4ErrorTrajState (G4ErrorTrajState &&)
virtual	~G4ErrorTrajState ()
G4ErrorTrajState &	operator= (const G4ErrorTrajState &)
G4ErrorTrajState &	operator= (G4ErrorTrajState &&)
void	SetData (const G4String &partType, const G4Point3D &pos, const G4Vector3D &mom)
void	BuildCharge ()
virtual G4int	PropagateError (const G4Track *)
virtual G4int	Update (const G4Track *)
void	UpdatePosMom (const G4Point3D &pos, const G4Vector3D &mom)
void	DumpPosMomError (std::ostream &out=G4cout) const
const G4String &	GetParticleType () const
void	SetParticleType (const G4String &partType)
G4Point3D	GetPosition () const
G4Vector3D	GetMomentum () const
G4ErrorTrajErr	GetError () const
virtual void	SetError (G4ErrorTrajErr em)
G4Track *	GetG4Track () const
void	SetG4Track (G4Track *trk)
G4double	GetCharge () const
void	SetCharge (G4double ch)
virtual G4eTSType	GetTSType () const

Friends
std::ostream &	operator<< (std::ostream &out, const G4ErrorSurfaceTrajState &ts)

Additional Inherited Members
Protected Attributes inherited from G4ErrorTrajState
G4String	fParticleType
G4Point3D	fPosition
G4Vector3D	fMomentum
G4double	fCharge = 0.
G4ErrorTrajErr	fError
G4eTSType	theTSType
G4Track *	theG4Track = nullptr
G4int	iverbose = 0

Detailed Description

Definition at line 55 of file G4ErrorSurfaceTrajState.hh.

Constructor & Destructor Documentation

G4ErrorSurfaceTrajState() [1/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState	(	const G4String &	partType,
		const G4Point3D &	pos,
		const G4Vector3D &	mom,
		const G4Plane3D &	plane,
		const G4ErrorTrajErr &	errmat = G4ErrorTrajErr(5, 0) )

Definition at line 56 of file G4ErrorSurfaceTrajState.cc.

  : G4ErrorTrajState(partType, pos, mom, errmat)
{
  Init();
  fTrajParam = G4ErrorSurfaceTrajParam(pos, mom, plane);
}

G4ErrorSurfaceTrajState() [2/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState	(	const G4String &	partType,
		const G4Point3D &	pos,
		const G4Vector3D &	mom,
		const G4Vector3D &	vecV,
		const G4Vector3D &	vecW,
		const G4ErrorTrajErr &	errmat = G4ErrorTrajErr(5, 0) )

Definition at line 46 of file G4ErrorSurfaceTrajState.cc.

  : G4ErrorTrajState(partType, pos, mom, errmat)
{
  Init();
  fTrajParam = G4ErrorSurfaceTrajParam(pos, mom, vecU, vecV);
}

G4ErrorSurfaceTrajState() [3/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState	(	G4ErrorFreeTrajState &	tpSC,
		const G4Plane3D &	plane )

Definition at line 68 of file G4ErrorSurfaceTrajState.cc.

  : G4ErrorTrajState(tpSC.GetParticleType(), tpSC.GetPosition(),
                     tpSC.GetMomentum())
{
  //  fParticleType = tpSC.GetParticleType();
  //  fPosition = tpSC.GetPosition();
  //  fMomentum = tpSC.GetMomentum();
  fTrajParam = G4ErrorSurfaceTrajParam(fPosition, fMomentum, plane);
  Init();
 
  //----- Get the error matrix in SC coordinates
  BuildErrorMatrix(tpSC, GetVectorV(), GetVectorW());
}

G4ErrorSurfaceTrajState() [4/4]

G4ErrorSurfaceTrajState::G4ErrorSurfaceTrajState	(	G4ErrorFreeTrajState &	tpSC,
		const G4Vector3D &	vecV,
		const G4Vector3D &	vecW,
		G4ErrorMatrix &	transfM )

Definition at line 84 of file G4ErrorSurfaceTrajState.cc.

  : G4ErrorTrajState(tpSC.GetParticleType(), tpSC.GetPosition(),
                     tpSC.GetMomentum())
{
  Init();  // needed to define charge sign
  fTrajParam = G4ErrorSurfaceTrajParam(fPosition, fMomentum, vecU, vecV);
  //----- Get the error matrix in SC coordinates
  transfM = BuildErrorMatrix(tpSC, vecU, vecV);
}

~G4ErrorSurfaceTrajState()

G4ErrorSurfaceTrajState::~G4ErrorSurfaceTrajState ( )

inline

Definition at line 78 of file G4ErrorSurfaceTrajState.hh.

{}

Member Function Documentation

BuildErrorMatrix()

G4ErrorMatrix G4ErrorSurfaceTrajState::BuildErrorMatrix	(	G4ErrorFreeTrajState &	tpSC,
		const G4Vector3D &	vecV,
		const G4Vector3D &	vecW )

Definition at line 98 of file G4ErrorSurfaceTrajState.cc.

{
  G4double sclambda = tpSC.GetParameters().GetLambda();
  G4double scphi    = tpSC.GetParameters().GetPhi();
  if(G4ErrorPropagatorData::GetErrorPropagatorData()->GetMode() ==
     G4ErrorMode_PropBackwards)
  {
    sclambda *= -1;
    scphi += CLHEP::pi;
  }
  G4double cosLambda = std::cos(sclambda);
  G4double sinLambda = std::sin(sclambda);
  G4double sinPhi    = std::sin(scphi);
  G4double cosPhi    = std::cos(scphi);
 
#ifdef G4EVERBOSE
  if(iverbose >= 4)
    G4cout << " PM " << fMomentum.mag() << " pLambda " << sclambda << " pPhi "
           << scphi << G4endl;
#endif
 
  G4ThreeVector vTN(cosLambda * cosPhi, cosLambda * sinPhi, sinLambda);
  G4ThreeVector vUN(-sinPhi, cosPhi, 0.);
  G4ThreeVector vVN(-vTN.z() * vUN.y(), vTN.z() * vUN.x(), cosLambda);
 
#ifdef G4EVERBOSE
  if(iverbose >= 4)
    G4cout << " SC2SD: vTN " << vTN << " vUN " << vUN << " vVN " << vVN
           << G4endl;
#endif
  G4double UJ = vUN * GetVectorV();
  G4double UK = vUN * GetVectorW();
  G4double VJ = vVN * GetVectorV();
  G4double VK = vVN * GetVectorW();
 
  //--- Get transformation first
  G4ErrorMatrix transfM(5, 5, 0);
  //--- Get magnetic field
  const G4Field* field = G4TransportationManager::GetTransportationManager()
                           ->GetFieldManager()
                           ->GetDetectorField();
 
  G4Vector3D vectorU = GetVectorV().cross(GetVectorW());
  G4double T1R       = 1. / (vTN * vectorU);
 
#ifdef G4EVERBOSE
  if(iverbose >= 4)
    G4cout << "surf vectors:U,V,W " << vectorU << " " << GetVectorV() << " "
           << GetVectorW() << "  T1R:" << T1R << G4endl;
#endif
 
  if(fCharge != 0 && field)
  {
    G4double pos[3];
    pos[0] = fPosition.x() * cm;
    pos[1] = fPosition.y() * cm;
    pos[2] = fPosition.z() * cm;
    G4double Hd[3];
    field->GetFieldValue(pos, Hd);
    G4ThreeVector H =
      G4ThreeVector(Hd[0], Hd[1], Hd[2]) / tesla * 10.;  // in kilogauus
    G4double magH  = H.mag();
    G4double invP  = 1. / (fMomentum.mag() / GeV);
    G4double magHM = magH * invP;
    if(magH != 0.)
    {
      G4double magHM2 = fCharge / magH;
      G4double Q      = -magHM * c_light / (km / ns);
#ifdef G4EVERBOSE
      if(iverbose >= 4)
        G4cout << GeV << " Q " << Q << " magHM " << magHM << " c_light/(km/ns) "
               << c_light / (km / ns) << G4endl;
#endif
 
      G4double sinz = -H * vUN * magHM2;
      G4double cosz = H * vVN * magHM2;
      G4double T3R  = Q * std::pow(T1R, 3);
      G4double UI   = vUN * vectorU;
      G4double VI   = vVN * vectorU;
#ifdef G4EVERBOSE
      if(iverbose >= 4)
      {
        G4cout << " T1R " << T1R << " T3R " << T3R << G4endl;
        G4cout << " UI " << UI << " VI " << VI << " vectorU " << vectorU
               << G4endl;
        G4cout << " UJ " << UJ << " VJ " << VJ << G4endl;
        G4cout << " UK " << UK << " VK " << VK << G4endl;
      }
#endif
 
      transfM[1][3] = -UI * (VK * cosz - UK * sinz) * T3R;
      transfM[1][4] = -VI * (VK * cosz - UK * sinz) * T3R;
      transfM[2][3] = UI * (VJ * cosz - UJ * sinz) * T3R;
      transfM[2][4] = VI * (VJ * cosz - UJ * sinz) * T3R;
    }
  }
 
  G4double T2R  = T1R * T1R;
  transfM[0][0] = 1.;
  transfM[1][1] = -UK * T2R;
  transfM[1][2] = VK * cosLambda * T2R;
  transfM[2][1] = UJ * T2R;
  transfM[2][2] = -VJ * cosLambda * T2R;
  transfM[3][3] = VK * T1R;
  transfM[3][4] = -UK * T1R;
  transfM[4][3] = -VJ * T1R;
  transfM[4][4] = UJ * T1R;
 
#ifdef G4EVERBOSE
  if(iverbose >= 4)
    G4cout << " SC2SD transf matrix A= " << transfM << G4endl;
#endif
  fError = G4ErrorTrajErr(tpSC.GetError().similarity(transfM));
 
#ifdef G4EVERBOSE
  if(iverbose >= 1)
    G4cout << "G4EP: error matrix SC2SD " << fError << G4endl;
  if(iverbose >= 4)
    G4cout << "G4ErrorSurfaceTrajState from SC " << *this << G4endl;
#endif
 
  return transfM;  // want to use trasnfM for the reverse transform
}

Referenced by G4ErrorSurfaceTrajState(), and G4ErrorSurfaceTrajState().

Dump()

void G4ErrorSurfaceTrajState::Dump ( std::ostream & out = G4cout ) const

virtual

Implements G4ErrorTrajState.

Definition at line 231 of file G4ErrorSurfaceTrajState.cc.

{ out << *this; }

GetParameters()

G4ErrorSurfaceTrajParam G4ErrorSurfaceTrajState::GetParameters ( ) const

inline

Definition at line 91 of file G4ErrorSurfaceTrajState.hh.

{ return fTrajParam; }

Referenced by G4ErrorFreeTrajState::G4ErrorFreeTrajState().

GetVectorV()

G4Vector3D G4ErrorSurfaceTrajState::GetVectorV ( ) const

inline

Definition at line 108 of file G4ErrorSurfaceTrajState.hh.

{ return fTrajParam.GetVectorV(); }

Referenced by BuildErrorMatrix(), G4ErrorFreeTrajState::G4ErrorFreeTrajState(), G4ErrorSurfaceTrajState(), SetMomentum(), and SetPosition().

GetVectorW()

G4Vector3D G4ErrorSurfaceTrajState::GetVectorW ( ) const

inline

Definition at line 110 of file G4ErrorSurfaceTrajState.hh.

{ return fTrajParam.GetVectorW(); }

Referenced by BuildErrorMatrix(), G4ErrorFreeTrajState::G4ErrorFreeTrajState(), G4ErrorSurfaceTrajState(), SetMomentum(), and SetPosition().

SetMomentum()

virtual void G4ErrorSurfaceTrajState::SetMomentum ( const G4Vector3D & mom )

inlinevirtual

Reimplemented from G4ErrorTrajState.

Definition at line 117 of file G4ErrorSurfaceTrajState.hh.

  {
    SetParameters(fPosition, mom, GetVectorV(), GetVectorW());
  }

SetParameters() [1/2]

void G4ErrorSurfaceTrajState::SetParameters ( const G4Point3D & pos, const G4Vector3D & mom, const G4Plane3D & plane )

inline

Definition at line 100 of file G4ErrorSurfaceTrajState.hh.

  {
    fPosition = pos;
    fMomentum = mom;
    fTrajParam.SetParameters(pos, mom, plane);
  }

SetParameters() [2/2]

void G4ErrorSurfaceTrajState::SetParameters ( const G4Point3D & pos, const G4Vector3D & mom, const G4Vector3D & vecV, const G4Vector3D & vecW )

inline

Definition at line 92 of file G4ErrorSurfaceTrajState.hh.

  {
    fPosition = pos;
    fMomentum = mom;
    fTrajParam.SetParameters(pos, mom, vecV, vecW);
  }

Referenced by SetMomentum(), and SetPosition().

SetPosition()

virtual void G4ErrorSurfaceTrajState::SetPosition ( const G4Point3D pos )

inlinevirtual

Reimplemented from G4ErrorTrajState.

Definition at line 112 of file G4ErrorSurfaceTrajState.hh.

  {
    SetParameters(pos, fMomentum, GetVectorV(), GetVectorW());
  }

Friends And Related Symbol Documentation

operator<<

std::ostream & operator<< ( std::ostream & out, const G4ErrorSurfaceTrajState & ts )

friend

Definition at line 234 of file G4ErrorSurfaceTrajState.cc.

{
  std::ios::fmtflags oldFlags = out.flags();
  out.setf(std::ios::fixed, std::ios::floatfield);
 
  ts.DumpPosMomError(out);
 
  out << " G4ErrorSurfaceTrajState: Params: " << ts.fTrajParam << G4endl;
  out.flags(oldFlags);
  return out;
}

---

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

• G4ErrorSurfaceTrajState.hh
• G4ErrorSurfaceTrajState.cc