57 vxPerp = mom/mom.
mag();
60 vyPerp /= vyPerp.
mag();
62 vzPerp /= vzPerp.
mag();
66 if( vyPerp.
mag() != 0. ) {
70 fYPerp = posv.
dot( vyPerp );
71 fZPerp = posv.
dot( vzPerp );
89 G4int oldprc = out.precision(8);
90 out <<
" InvP= " << tp.fInvP <<
" Theta= "
91 << tp.fLambda <<
" Phi= " << tp.fPhi <<
" YPerp= " << tp.fYPerp
92 <<
" ZPerp= " << tp.fZPerp <<
G4endl;
93 out <<
" momentum direction= " << tp.fDir <<
G4endl;
94 out.precision(oldprc);
double dot(const Hep3Vector &) const
const G4ThreeVector & GetPosition() const
HepGeom::Vector3D< G4double > G4Vector3D
void SetParameters(const G4Point3D &pos, const G4Vector3D &mom)
void Update(const G4Track *aTrack)
G4ThreeVector GetMomentum() const
std::ostream & operator<<(std::ostream &, const BasicVector3D< float > &)
BasicVector3D< T > cross(const BasicVector3D< T > &v) const