69 if( vecV.mag() < kCarTolerance )
72 vecV = Zvec.cross(plane.normal());
75 G4Vector3D vecW = plane.normal().cross( vecV );
86 if( mom.mag() > 0. ) {
101 fPV = momv.dot( vecV );
102 fPW = momv.dot( vecW );
108 fV = posv.dot( vecV );
109 fW = posv.dot( vecW );
120 out <<
" InvP= " << tp.
fInvP <<
" PV= " << tp.
fPV
121 <<
" PW= " << tp.
fPW <<
" V= " << tp.
fV <<
" W= " << tp.
fW <<
G4endl;
122 out <<
" vectorV direction= " << tp.
fVectorV
std::ostream & operator<<(std::ostream &out, const G4ErrorSurfaceTrajParam &tp)
CLHEP::Hep3Vector G4ThreeVector
G4double GetSurfaceTolerance() const
HepGeom::Point3D< G4double > G4Point3D
HepGeom::Vector3D< G4double > G4Vector3D
void SetParameters(const G4Point3D &pos, const G4Vector3D &mom, const G4Vector3D &vecV, const G4Vector3D &vecW)
G4ErrorSurfaceTrajParam()
HepGeom::Plane3D< G4double > G4Plane3D
static const G4double pos
static G4GeometryTolerance * GetInstance()