---

Track.cxx

Go to the documentation of this file.

#include "AtlfastEvent/Track.h"

#include "GaudiKernel/MsgStream.h"

#include "GaudiKernel/ObjectList.h"
#include "GaudiKernel/ObjectVector.h"
#include "GaudiKernel/ObjectFactory.h"

#include <iomanip>
#include <iostream>

static ContainedObjectFactory < Atlfast::Track > s_AtlfastTrackFactory;
const IFactory& AtlfastTrackFactory = s_AtlfastTrackFactory;

namespace Atlfast {
  using std::setw;
  using std::setprecision;
  using std::setiosflags;
  using std::ios;

// ................................................
// Message stream operator for simple debug output

MsgStream& operator << ( MsgStream& s, const Track& t ) 
{
  if (t.truth()) {
    s << "Track created from particle " << t.pdg_id() << endreq << " with (px,py,pz) " 
    <<t.truth()->momentum().px() <<" "
    <<t.truth()->momentum().py() <<" "
    <<t.truth()->momentum().pz() <<" "
    <<endreq<<endreq;
  } 
  TrackTrajectory traj = t.trajectory();
  TrackParameters param = traj.parameters();
  TrackTrajectory truthTraj = t.truthTrajectory();
  TrackParameters truthParam = truthTraj.parameters();
  
  s << "     unsmeared Track       smeared Track"<<endreq; 

    s<<setiosflags(ios::fixed);
    s<<setprecision(7);
    s<<setw(20);

    s <<setw(20)<< truthParam.pT() 
      <<setw(20)<< param.pT()
      << "        pT"<<endreq;
    s <<setw(20)<< truthParam.eta() 
      <<setw(20)<< param.eta()
      << "        eta"<<endreq;
    s << "------------------  -------------------" << endreq;

    s <<setw(20)<< truthParam.impactParameter() 
      <<setw(20)<< param.impactParameter()
      << "        d0"<<endreq; 
    s <<setw(20)<< truthParam.zPerigee() 
      <<setw(20)<< param.zPerigee()
      << "        z0" << endreq;  
    s <<setw(20)<< truthParam.phi().val()
      <<setw(20)<< param.phi().val()
      << "       phi" <<endreq;
     s <<setw(20)<< truthParam.cotTheta()
      <<setw(20)<<param.cotTheta()
      << " cot Theta" << endreq;
    s <<setw(20)<< truthParam.invPtCharge()
      <<setw(20)<< param.invPtCharge()
      << "  1/(q.PT)" <<endreq;
    s << "-------------------  -------------------" << endreq;
    s <<setw(20)<< truthTraj.radius()
      <<setw(20)<< traj.radius()
      << "    radius" <<endreq;
    s <<setw(20)<< truthTraj.curvature()
      <<setw(20)<< traj.curvature()
      << " curvature" <<endreq<<endreq;
        
    s << "Matrix used"<<endreq ;
    s << setiosflags(2048);
    for (int a = 1; a<6; a++)
      {
        for (int b=1; b<6; b++)
          {
            s <<setw(10)<< t.smearMatrix() (a,b);
          }
        s <<endreq;
      }

  
  return s ;
}

MsgStream& operator << ( MsgStream& s, const Track* t ) 
{
  return s << *t ;
}

} // end namespace

---