//************************************************************* int fitter::GetQ(const Trajectory& traj){ //************************************************************* if (_model.compare("particle/helix")!=0) return 0; double q; /// q = traj.state(traj.last_fitted_node()).vector()[dim-1]; if (q<0) q=-1; else q=1; return (int) q; }