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 () override=default
virtual G4ThreeVector	IntersectLocal (const G4ThreeVector &point, const G4ThreeVector &direc) const
G4double	GetDistanceFromPoint (const G4ThreeVector &point, const G4ThreeVector &direc) const override
G4double	GetDistanceFromPoint (const G4ThreeVector &point) const override
G4Plane3D	GetTangentPlane (const G4ThreeVector &point) const override
void	Dump (const G4String &msg) const override
Public Member Functions inherited from G4ErrorSurfaceTarget
	G4ErrorSurfaceTarget ()
	~G4ErrorSurfaceTarget () override=default
Public Member Functions inherited from G4ErrorTanPlaneTarget
	G4ErrorTanPlaneTarget ()
	~G4ErrorTanPlaneTarget () override=default
Public Member Functions inherited from G4ErrorTarget
	G4ErrorTarget ()
virtual	~G4ErrorTarget ()=default
virtual G4bool	TargetReached (const G4Step *)
G4ErrorTargetType	GetType () const

Additional Inherited Members
Protected Attributes inherited from G4ErrorTarget
G4ErrorTargetType	theType {G4ErrorTarget_GeomVolume}

Detailed Description

Definition at line 44 of file G4ErrorCylSurfaceTarget.hh.

Constructor & Destructor Documentation

G4ErrorCylSurfaceTarget() [1/2]

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

Definition at line 44 of file G4ErrorCylSurfaceTarget.cc.

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

G4ErrorCylSurfaceTarget() [2/2]

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

Definition at line 63 of file G4ErrorCylSurfaceTarget.cc.

  : fradius(radius), ftransform(trans.Inverse())
{
  theType = G4ErrorTarget_CylindricalSurface;
 
#ifdef G4VERBOSE
  if(G4ErrorPropagatorData::verbose() >= 2 )
  { 
    Dump( " $$$ creating G4ErrorCylSurfaceTarget ");
  }
#endif
}

~G4ErrorCylSurfaceTarget()

G4ErrorCylSurfaceTarget::~G4ErrorCylSurfaceTarget ( )

overridedefault

Member Function Documentation

Dump()

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

overridevirtual

Implements G4ErrorSurfaceTarget.

Definition at line 215 of file G4ErrorCylSurfaceTarget.cc.

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

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

GetDistanceFromPoint() [1/2]

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

overridevirtual

Implements G4ErrorSurfaceTarget.

Definition at line 114 of file G4ErrorCylSurfaceTarget.cc.

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

GetDistanceFromPoint() [2/2]

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

overridevirtual

Implements G4ErrorSurfaceTarget.

Definition at line 80 of file G4ErrorCylSurfaceTarget.cc.

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

GetTangentPlane()

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

overridevirtual

Implements G4ErrorSurfaceTarget.

Definition at line 187 of file G4ErrorCylSurfaceTarget.cc.

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

IntersectLocal()

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

virtual

Definition at line 134 of file G4ErrorCylSurfaceTarget.cc.

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

Referenced by GetDistanceFromPoint().

---

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

• G4ErrorCylSurfaceTarget.hh
• G4ErrorCylSurfaceTarget.cc