G4ErrorCylSurfaceTarget Class Reference

#include <G4ErrorCylSurfaceTarget.hh>

Inheritance diagram for G4ErrorCylSurfaceTarget:

Public Member Functions
	G4ErrorCylSurfaceTarget (const G4double &radius, const G4ThreeVector &trans=G4ThreeVector(), const G4RotationMatrix &rotm=G4RotationMatrix())
	G4ErrorCylSurfaceTarget (const G4double &radius, const G4AffineTransform &trans)
	~G4ErrorCylSurfaceTarget ()
virtual G4ThreeVector	IntersectLocal (const G4ThreeVector &point, const G4ThreeVector &direc) const
virtual G4double	GetDistanceFromPoint (const G4ThreeVector &point, const G4ThreeVector &direc) const
virtual G4double	GetDistanceFromPoint (const G4ThreeVector &point) const
virtual G4Plane3D	GetTangentPlane (const G4ThreeVector &point) const
virtual void	Dump (const G4String &msg) const
Public Member Functions inherited from G4ErrorSurfaceTarget
	G4ErrorSurfaceTarget ()
virtual	~G4ErrorSurfaceTarget ()
Public Member Functions inherited from G4ErrorTanPlaneTarget
	G4ErrorTanPlaneTarget ()
virtual	~G4ErrorTanPlaneTarget ()
Public Member Functions inherited from G4ErrorTarget
	G4ErrorTarget ()
virtual	~G4ErrorTarget ()
virtual G4bool	TargetReached (const G4Step *)
G4ErrorTargetType	GetType () const

Additional Inherited Members
Protected Attributes inherited from G4ErrorTarget
G4ErrorTargetType	theType

Detailed Description

Definition at line 52 of file G4ErrorCylSurfaceTarget.hh.

Constructor & Destructor Documentation

G4ErrorCylSurfaceTarget::G4ErrorCylSurfaceTarget	(	const G4double &	radius,
		const G4ThreeVector &	trans = G4ThreeVector(),
		const G4RotationMatrix &	rotm = G4RotationMatrix()
	)

Definition at line 47 of file G4ErrorCylSurfaceTarget.cc.

References Dump(), G4ErrorTarget_CylindricalSurface, CLHEP::HepRotation::inverse(), G4ErrorTarget::theType, and G4ErrorPropagatorData::verbose().

  : fradius(radius)
{
  theType = G4ErrorTarget_CylindricalSurface;

  ftransform = G4AffineTransform( rotm.inverse(), -trans );
#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 2 ) { 
    Dump( " $$$ creating G4ErrorCylSurfaceTarget ");
  }
#endif
}

G4ErrorCylSurfaceTarget::G4ErrorCylSurfaceTarget	(	const G4double &	radius,
		const G4AffineTransform &	trans
	)

Definition at line 65 of file G4ErrorCylSurfaceTarget.cc.

References Dump(), G4ErrorTarget_CylindricalSurface, G4ErrorTarget::theType, and G4ErrorPropagatorData::verbose().

  : fradius(radius), ftransform(trans.Inverse())
{
  theType = G4ErrorTarget_CylindricalSurface;

#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 2 )
  { 
    Dump( " $$$ creating G4ErrorCylSurfaceTarget ");
  }
#endif
}

G4ErrorCylSurfaceTarget::~G4ErrorCylSurfaceTarget

(

)

Definition at line 81 of file G4ErrorCylSurfaceTarget.cc.

{
}

Member Function Documentation

void G4ErrorCylSurfaceTarget::Dump ( const G4String &  msg ) const

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 222 of file G4ErrorCylSurfaceTarget.cc.

References G4cout, G4endl, G4AffineTransform::NetRotation(), and G4AffineTransform::NetTranslation().

Referenced by G4ErrorCylSurfaceTarget(), GetDistanceFromPoint(), and IntersectLocal().

{
  G4cout << msg << " radius " << fradius
                << " centre " << ftransform.NetTranslation()
                << " rotation " << ftransform.NetRotation() << G4endl;
}

G4double G4ErrorCylSurfaceTarget::GetDistanceFromPoint ( const G4ThreeVector &  point, const G4ThreeVector &  direc  ) const

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 88 of file G4ErrorCylSurfaceTarget.cc.

References Dump(), FatalException, G4cout, G4endl, G4Exception(), IntersectLocal(), CLHEP::Hep3Vector::mag(), G4AffineTransform::TransformAxis(), G4AffineTransform::TransformPoint(), and G4ErrorPropagatorData::verbose().

{
  if( dir.mag() == 0. )
  {
    G4Exception("G4ErrorCylSurfaceTarget::GetDistanceFromPoint()",
                "GeomMgt0003", FatalException, "Direction is zero !");
  }

  //----- Get intersection point
  G4ThreeVector localPoint = ftransform.TransformPoint( point );
  G4ThreeVector localDir = ftransform.TransformAxis( dir ); 
  G4ThreeVector inters = IntersectLocal(localPoint, localDir);

  G4double dist = (localPoint-inters).mag();
  
#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 3 )
  { 
    G4cout << " G4ErrorCylSurfaceTarget::GetDistanceFromPoint():" << G4endl
           << " Global point " << point << " dir " << dir << G4endl
           << " Intersection " << inters << G4endl
           << " Distance " << dist << G4endl;
    Dump( " CylSurface: " );
  }
#endif

  return dist;
}

G4double G4ErrorCylSurfaceTarget::GetDistanceFromPoint ( const G4ThreeVector &  point ) const

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 122 of file G4ErrorCylSurfaceTarget.cc.

References Dump(), G4cout, G4endl, CLHEP::Hep3Vector::perp(), G4AffineTransform::TransformPoint(), and G4ErrorPropagatorData::verbose().

{
  G4ThreeVector localPoint = ftransform.TransformPoint( point );

#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 3 )
  { 
    G4cout << " G4ErrorCylSurfaceTarget::GetDistanceFromPoint:" << G4endl
           << " Global point " << point << G4endl
           << " Distance " << fradius - localPoint.perp() << G4endl;
    Dump( " CylSurface: " );
  }
#endif

  return fradius - localPoint.perp();
}

G4Plane3D G4ErrorCylSurfaceTarget::GetTangentPlane ( const G4ThreeVector &  point ) const

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 195 of file G4ErrorCylSurfaceTarget.cc.

References G4endl, G4Exception(), G4GeometryTolerance::GetInstance(), G4GeometryTolerance::GetSurfaceTolerance(), JustWarning, G4AffineTransform::NetTranslation(), CLHEP::Hep3Vector::perp(), and G4AffineTransform::TransformPoint().

{
  G4ThreeVector localPoint = ftransform.TransformPoint( point );

  // check that point is at cylinder surface
  //
  if( std::fabs( localPoint.perp() - fradius )
      > 1000.*G4GeometryTolerance::GetInstance()->GetSurfaceTolerance() )
  {
    std::ostringstream message;
    message << "Local point not at surface !" << G4endl
            << "          Point: " << point << ", local: " << localPoint
            << G4endl
            << "          is not at surface, but far away by: "
            << localPoint.perp() - fradius << " !";
    G4Exception("G4ErrorCylSurfaceTarget::GetTangentPlane()",
                "GeomMgt1002", JustWarning, message);
  }

  G4Normal3D normal = localPoint - ftransform.NetTranslation();

  return G4Plane3D( normal, point );
}

G4ThreeVector G4ErrorCylSurfaceTarget::IntersectLocal ( const G4ThreeVector &  point, const G4ThreeVector &  direc  ) const

virtual

Definition at line 142 of file G4ErrorCylSurfaceTarget.cc.

References Dump(), G4cout, G4endl, G4Exception(), JustWarning, G4InuclParticleNames::lambda, CLHEP::Hep3Vector::mag(), CLHEP::Hep3Vector::perp(), G4ErrorPropagatorData::verbose(), CLHEP::Hep3Vector::x(), and CLHEP::Hep3Vector::y().

Referenced by GetDistanceFromPoint().

{
  G4double eqa = localDir.x()*localDir.x()+localDir.y()*localDir.y();
  G4double eqb = 2*(localPoint.x()*localDir.x()+localPoint.y()*localDir.y());
  G4double eqc = -fradius*fradius+localPoint.x()*localPoint.x()
                 +localPoint.y()*localPoint.y();
  G4int inside = (localPoint.perp() > fradius) ? -1 : 1;
  G4double lambda;

  if( eqa*inside > 0. )
  {
    lambda = (-eqb + std::sqrt(eqb*eqb-4*eqa*eqc) ) / (2.*eqa);
  }
  else if( eqa*inside < 0. )
  {
    lambda = (-eqb - std::sqrt(eqb*eqb-4*eqa*eqc) ) / (2.*eqa);
  }
  else
  {
    if( eqb != 0. )
    {
      lambda = -eqc/eqb;
    }
    else
    {
      std::ostringstream message;
      message << "Intersection not possible !" << G4endl
              << "          Point: " << localPoint << ", direction: "
              << localDir;
      Dump( " CylSurface: " ); 
      G4Exception("G4ErrorCylSurfaceTarget::IntersectLocal()",
                  "GeomMgt1002", JustWarning, message);
      lambda = kInfinity;
    }
  }

  G4ThreeVector inters = localPoint + lambda*localDir/localDir.mag();

#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 4 ) { 
    G4cout << " G4ErrorCylSurfaceTarget::IntersectLocal " << inters << " "
           << inters.perp() << " localPoint " << localPoint << " localDir "
           << localDir << G4endl;
  } 
#endif

  return inters;
}

---

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

• G4ErrorCylSurfaceTarget.hh
• G4ErrorCylSurfaceTarget.cc