Example #1
0
//*************************************************************
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;
  
}