#include <G4ErrorCylSurfaceTarget.hh>

Inheritance diagram for G4ErrorCylSurfaceTarget:

Public Member Functions
virtual void	Dump (const G4String &msg) const

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

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

virtual G4double	GetDistanceFromPoint (const G4ThreeVector &point) const

virtual G4double	GetDistanceFromPoint (const G4ThreeVector &point, const G4ThreeVector &direc) const

virtual G4Plane3D	GetTangentPlane (const G4ThreeVector &point) const

G4ErrorTargetType	GetType () const

virtual G4ThreeVector	IntersectLocal (const G4ThreeVector &point, const G4ThreeVector &direc) const

virtual G4bool	TargetReached (const G4Step *)

	~G4ErrorCylSurfaceTarget ()=default

Protected Attributes
G4ErrorTargetType	theType

Private Attributes
G4double	fradius

G4AffineTransform	ftransform

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 43 of file G4ErrorCylSurfaceTarget.cc.

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

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

◆ G4ErrorCylSurfaceTarget() [2/2]

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

Definition at line 62 of file G4ErrorCylSurfaceTarget.cc.

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

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

◆ ~G4ErrorCylSurfaceTarget()

G4ErrorCylSurfaceTarget::~G4ErrorCylSurfaceTarget ( )

default

Member Function Documentation

◆ Dump()

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

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 214 of file G4ErrorCylSurfaceTarget.cc.

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

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

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

◆ GetDistanceFromPoint() [1/2]

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

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 113 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();
}

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

◆ GetDistanceFromPoint() [2/2]

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

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 79 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;
}

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

◆ GetTangentPlane()

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

virtual

Implements G4ErrorSurfaceTarget.

Definition at line 186 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 );
}

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

◆ GetType()

G4ErrorTargetType G4ErrorTarget::GetType ( ) const

inlineinherited

Definition at line 74 of file G4ErrorTarget.hh.

{
  return theType;
}

References G4ErrorTarget::theType.

Referenced by G4ErrorPropagator::GetFinalTrajState(), G4ErrorPropagationNavigator::GetGlobalExitNormal(), and G4ErrorPropagator::MakeOneStep().

◆ IntersectLocal()

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

virtual

Definition at line 133 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;
}

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

Referenced by GetDistanceFromPoint().

◆ TargetReached()

G4bool G4ErrorTarget::TargetReached ( const G4Step * )

virtualinherited

Reimplemented in G4ErrorGeomVolumeTarget.

Definition at line 46 of file G4ErrorTarget.cc.

{
  return false;
}

Field Documentation

◆ fradius

G4double G4ErrorCylSurfaceTarget::fradius

private

Definition at line 80 of file G4ErrorCylSurfaceTarget.hh.

Referenced by Dump(), GetDistanceFromPoint(), GetTangentPlane(), and IntersectLocal().

◆ ftransform

G4AffineTransform G4ErrorCylSurfaceTarget::ftransform

private

Definition at line 81 of file G4ErrorCylSurfaceTarget.hh.

Referenced by Dump(), G4ErrorCylSurfaceTarget(), GetDistanceFromPoint(), and GetTangentPlane().

◆ theType

G4ErrorTargetType G4ErrorTarget::theType

protectedinherited

Definition at line 69 of file G4ErrorTarget.hh.

Referenced by G4ErrorCylSurfaceTarget(), G4ErrorGeomVolumeTarget::G4ErrorGeomVolumeTarget(), G4ErrorPlaneSurfaceTarget::G4ErrorPlaneSurfaceTarget(), G4ErrorTrackLengthTarget::G4ErrorTrackLengthTarget(), and G4ErrorTarget::GetType().

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

source/geometry/management/include/G4ErrorCylSurfaceTarget.hh
source/geometry/management/src/G4ErrorCylSurfaceTarget.cc

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4ErrorCylSurfaceTarget() [1/2]

◆ G4ErrorCylSurfaceTarget() [2/2]

◆ ~G4ErrorCylSurfaceTarget()

Member Function Documentation

◆ Dump()

◆ GetDistanceFromPoint() [1/2]

◆ GetDistanceFromPoint() [2/2]

◆ GetTangentPlane()

◆ GetType()

◆ IntersectLocal()

◆ TargetReached()

Field Documentation

◆ fradius

◆ ftransform

◆ theType