TrgTrackStreamer.cpp

Go to the documentation of this file.

#include "TrgConverter/TrgTrackStreamer.h"

std::ostream& str::trgTrack(const TrgTrack& gtrack, std::ostream& os) {
  
  os << str::trgTrackHeader(gtrack,os) << std::endl;
  os << str::trgTrackFields(gtrack,os) << std::endl;
  os << str::trgStates(gtrack.otherStates(),os) <<std::endl;
  os << str::trgTrackReferences(gtrack,os);
  return os;
}

std::ostream& str::trgTrackHeader(const TrgTrack& gtrack, std::ostream& os) 
{
  os << " chi2  \t" << gtrack.chiSquared() << std::endl;
  os << " ndof  \t" << gtrack.degreeOfFreedom() << std::endl;
  os << " first state \t" 
     << str::trgState(gtrack.firstState(),os) << std::endl;
  os << " flags \t " << gtrack.flags() << std::endl;
  os << " nStates \t" << gtrack.otherStates().size();
  return os;
}

std::ostream& str::trgStates(const std::vector<TrgState>& states,
                             std::ostream& os)  
{
  os << " states \t " << states.size() << std::endl;
  for (std::vector<TrgState>::const_iterator it = states.begin();
       it != states.end(); it++) {
    const TrgState& gstate = *it;
    os << str::trgState(gstate,os) << std::endl;
  }
  return os;
}

std::ostream& str::trgTrackReferences(const TrgTrack& gtrack,
                                      std::ostream& os) 
{
  os << " references :" << std::endl;
  const std::vector< std::pair<VeloChannelID,int> >& vref 
    = gtrack.veloReference();
  for (std::vector< std::pair<VeloChannelID,int> >::const_iterator
         it = vref.begin(); it != vref.end(); it++) {
    const std::pair<VeloChannelID,int>& velo = *it;
    os << " \t velo id "<<velo.first.channelID() 
       << " size " << velo.second << std::endl;
  }  
  const std::vector< std::pair<ITChannelID,int> >& itref 
    = gtrack.itReference();
  for (std::vector< std::pair<ITChannelID,int> >::const_iterator
         it = itref.begin(); it != itref.end(); it++) {
    const std::pair<ITChannelID,int>& si = *it;
    os << " \t si id "<< si.first.channelID() 
       << " size " << si.second << std::endl;
  }
  const std::vector< std::pair<OTChannelID,int> >& otref 
    = gtrack.otReference();
  for (std::vector< std::pair<OTChannelID,int> >::const_iterator
         it = otref.begin(); it != otref.end(); it++) {
    const std::pair<VeloChannelID,int>& ot = *it;
    os << " \t ot id "<< ot.first.channelID() 
       << " size " << ot.second << std::endl;
  }
  return os; 
}


std::ostream& str::trgTrackFields(const TrgTrack& gtrack, std::ostream& os) {
  
  
  TrgTrack::TrgTrackType gtype = (TrgTrack::TrgTrackType) gtrack.type();
  
  if (gtype == TrgTrack::TypeVeloR) os << "T:VeloR, ";
  if (gtype == TrgTrack::TypeVelo3D) os << "T:Velo3D, ";
  if (gtype == TrgTrack::TypeVeloTT) os << "T:VeloTT, ";
  if (gtype == TrgTrack::TypeLong) os << "T:Long, ";
  if (gtype == TrgTrack::TypeDownstream) os << "T:Downstream, ";
  if (gtype == TrgTrack::TypeKShort) os << "T:KShort, ";
  
  if (gtrack.badTrack()) os << "F:BadTrack, ";
  if (gtrack.alreadyUsed()) os << "F:AlreadyUsed, ";
  if (gtrack.backward()) os << "F:Backward, ";
  if (gtrack.selectedIP2D()) os << "F:SelectedIP2D, ";
  if (gtrack.selectedMuon2D()) os << "F:SelectedMuon2D, ";
  
  return os;
 }

std::ostream& str::trgState(const TrgState& state, std::ostream& os) {
  os << " z " << state.z() << std::endl;
  os << " state " << state.stateVector();
  //  os << " cov " << state.covariance() << std::endl;
  return os;
}
 

---