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= "
93 out <<
" momentum direction= " << tp.
fDir <<
G4endl;
94 out.precision(oldprc);
CLHEP::Hep3Vector G4ThreeVector
const G4ThreeVector & GetPosition() const
HepGeom::Point3D< G4double > G4Point3D
HepGeom::Vector3D< G4double > G4Vector3D
void SetParameters(const G4Point3D &pos, const G4Vector3D &mom)
std::ostream & operator<<(std::ostream &out, const G4ErrorFreeTrajParam &tp)
void Update(const G4Track *aTrack)
G4ThreeVector GetMomentum() const
static const G4double pos