57 Dump(
" $$$ creating G4ErrorCylSurfaceTarget ");
67 : fradius(radius), ftransform(trans.Inverse())
74 Dump(
" $$$ creating G4ErrorCylSurfaceTarget ");
93 G4Exception(
"G4ErrorCylSurfaceTarget::GetDistanceFromPoint()",
102 G4double dist = (localPoint-inters).mag();
107 G4cout <<
" G4ErrorCylSurfaceTarget::GetDistanceFromPoint():" <<
G4endl
108 <<
" Global point " << point <<
" dir " << dir <<
G4endl
109 <<
" Intersection " << inters <<
G4endl
110 <<
" Distance " << dist <<
G4endl;
111 Dump(
" CylSurface: " );
129 G4cout <<
" G4ErrorCylSurfaceTarget::GetDistanceFromPoint:" <<
G4endl
130 <<
" Global point " << point <<
G4endl
131 <<
" Distance " << fradius - localPoint.
perp() <<
G4endl;
132 Dump(
" CylSurface: " );
136 return fradius - localPoint.
perp();
145 G4double eqa = localDir.
x()*localDir.
x()+localDir.
y()*localDir.
y();
146 G4double eqb = 2*(localPoint.
x()*localDir.
x()+localPoint.
y()*localDir.
y());
147 G4double eqc = -fradius*fradius+localPoint.
x()*localPoint.
x()
148 +localPoint.
y()*localPoint.
y();
149 G4int inside = (localPoint.
perp() > fradius) ? -1 : 1;
152 if( eqa*inside > 0. )
154 lambda = (-eqb + std::sqrt(eqb*eqb-4*eqa*eqc) ) / (2.*eqa);
156 else if( eqa*inside < 0. )
158 lambda = (-eqb - std::sqrt(eqb*eqb-4*eqa*eqc) ) / (2.*eqa);
168 std::ostringstream message;
169 message <<
"Intersection not possible !" <<
G4endl
170 <<
" Point: " << localPoint <<
", direction: "
172 Dump(
" CylSurface: " );
173 G4Exception(
"G4ErrorCylSurfaceTarget::IntersectLocal()",
183 G4cout <<
" G4ErrorCylSurfaceTarget::IntersectLocal " << inters <<
" "
184 << inters.
perp() <<
" localPoint " << localPoint <<
" localDir "
201 if( std::fabs( localPoint.
perp() - fradius )
204 std::ostringstream message;
205 message <<
"Local point not at surface !" <<
G4endl
206 <<
" Point: " << point <<
", local: " << localPoint
208 <<
" is not at surface, but far away by: "
209 << localPoint.
perp() - fradius <<
" !";
210 G4Exception(
"G4ErrorCylSurfaceTarget::GetTangentPlane()",
224 G4cout << msg <<
" radius " << fradius
G4double GetSurfaceTolerance() const
HepRotation inverse() const
G4ErrorCylSurfaceTarget(const G4double &radius, const G4ThreeVector &trans=G4ThreeVector(), const G4RotationMatrix &rotm=G4RotationMatrix())
G4GLOB_DLL std::ostream G4cout
~G4ErrorCylSurfaceTarget()
void G4Exception(const char *originOfException, const char *exceptionCode, G4ExceptionSeverity severity, const char *comments)
G4ErrorTargetType theType
virtual G4Plane3D GetTangentPlane(const G4ThreeVector &point) const
virtual G4ThreeVector IntersectLocal(const G4ThreeVector &point, const G4ThreeVector &direc) const
HepGeom::Plane3D< G4double > G4Plane3D
virtual void Dump(const G4String &msg) const
virtual G4double GetDistanceFromPoint(const G4ThreeVector &point, const G4ThreeVector &direc) const
static G4GeometryTolerance * GetInstance()