TrackStreamer.cpp

Go to the documentation of this file.

#include "Event/TrackStreamer.h"
#include "Event/TrackKeys.h"
#include "Event/StateKeys.h"

std::ostream& str::track(const Track& track, std::ostream& os) {

  os << str::trackHeader(track,os) << std::endl;
  os << str::trackFields(track,os) << std::endl;
  os << str::states(track.states(),os) << std::endl;
  os << str::lhcbIDs(track.lhcbIDs(),os);
  return os;
}

std::ostream& str::trackHeader(const Track& track, std::ostream& os) {

  os << " chi2  \t" << track.chi2() << std::endl;
  os << " ndof  \t" << track.nDoF() << std::endl;
  os << std::endl;
  os << " flags \t " << track.flags() << std::endl;
  os << " nstates \t " << track.states().size() << std::endl;
  os << " nIDs \t " << track.lhcbIDs().size();
  return os;
}


std::ostream& str::states(const std::vector<State*>& states,  
                          std::ostream& os) {
  for (std::vector<State*>::const_iterator it = states.begin(); 
       it != states.end(); it++) {
    const State& state = *(*it);
    str::state(state,os); os << std::endl;
  }
  return os;
}


std::ostream& str::state(const State& state, std::ostream& os) {
  os << " flags \t" << state.flags() << std::endl;
  str::stateFields(state,os); os << std::endl;
  os << " z \t" << state.z() << std::endl;
  os << " state \t" << state.state();
  //  os << " cov " << state.covariance() << std::endl;
  return os;
}

std::ostream& str::lhcbIDs(const std::vector<LHCbID>& ids, std::ostream& os) {
  for (std::vector<LHCbID>::const_iterator it = ids.begin();
       it != ids.end(); it++) 
    str::lhcbID(*it,os); os << std::endl;
  return os;
}

std::ostream& str::lhcbID(const LHCbID& id, std::ostream& os) {
  os << " ID detector " << id.detectorType() 
     << " channel " << id.channelID()
     << " size " << (int) id.spareBits();
  os << std::endl;
  return os;
}


std::ostream& str::trackFields(const Track& track, std::ostream& os) {
  
  os << " type : \t" << track.type() << std::endl;
  os << " status: \t" <<track.status() << std::endl;
  os << " history: \t" << track.history() << std::endl;
  os << " history Fit: \t" << track.historyFit() << std::endl;
  // os << " flags: \t";
  os << std::endl;
  return os;
}

std::ostream& str::stateFields(const State& state, std::ostream& os) {
  
  if (state.checkType(StateKeys::HasMomentum)) os <<"T:HasMomentum, ";
  if (state.checkType(StateKeys::StraightLine)) os <<"T:StraightLine, ";
  
  if (state.checkLocation(StateKeys::LocationUnknown)) os <<"L:Unknown, ";
  if (state.checkLocation(StateKeys::ClosestToBeam)) 
    os <<"L:ClosestToBeam, ";
  if (state.checkLocation(StateKeys::FirstMeasurement)) 
    os <<"L:FirstMeasurement, ";
  if (state.checkLocation(StateKeys::EndVelo)) os <<"L:EndVelo, ";
  if (state.checkLocation(StateKeys::AtTT)) os <<"L:AtTT, ";
  if (state.checkLocation(StateKeys::AtT)) os <<"L:AtT, ";
  if (state.checkLocation(StateKeys::BegRich1)) os <<"L:BegRich1, ";
  if (state.checkLocation(StateKeys::EndRich1)) os <<"L:EndRich1, ";
  if (state.checkLocation(StateKeys::BegRich2)) os <<"L:BegRich2, ";
  if (state.checkLocation(StateKeys::EndRich2)) os <<"L:EndRich2, ";
  if (state.checkLocation(StateKeys::Calo)) os <<"L:Calo, ";
  if (state.checkLocation(StateKeys::Muon)) os <<"L:Muon, ";

  return os;  
}

---