Track.cxx

Go to the documentation of this file.

#include "AtlfastEvent/Track.h"
#include "AtlfastEvent/IAOOvisitor.h"

#include "GaudiKernel/MsgStream.h"
#include <iomanip>
#include <iostream>

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(ios::scientific); 
    s << setiosflags(ios::left);
    s << setprecision(5);
    for (int a = 1; a<6; a++){
      for (int b=1; b<6; b++){
        s << setw(12) << t.smearMatrix() (a,b);
      }
      s <<endreq;
    }
    
    return s ;
  }
  void Track::accept(IAOOvisitor* iaPtr) const{
    iaPtr->process(this);
  }

  TrackTrajectory Track::trajectory() const {
    return m_trajectory;
  }
  int Track::pdg_id() const {return m_pdg_id;}
  const HepMC::GenParticle* Track::truth() const { return m_truth ; }
  const HepSymMatrix Track::smearMatrix() const {return m_smearMatrix;}
  void Track::setMatrix(const HepSymMatrix& matrix) {m_smearMatrix = matrix;}
  const TrackTrajectory Track::truthTrajectory() const {return m_truthTrajectory;}
  HepLorentzVector  Track::momentum() const {
    return m_truth->momentum();
  }
  double  Track::eta() const{ 
    return momentum().pseudoRapidity();
  }
  double  Track::phi() const{ 
    return momentum().phi();
  }
  double  Track::pT()  const{ 
    return momentum().perp();
  }
  double  Track::eT()  const{ 
    return momentum().perp();
  }
  double  Track::mT()  const{ 
    return momentum().mt();
  }
  IKinematic* Track::clone() const{
    IKinematic* ik = new Track(*this);
    return ik;
    }
  void Track::print(const std::string& coor, std::string t) const{
    AODNavigationImp::print(coor, t);
    std::cout<<s<<"\n"
             <<std::string(73,' ')
             <<"Track"<<std::endl;
  }
 
  MsgStream& operator << ( MsgStream& s, const Track* t ) {
    return s << *t ;
  }

} // end namespace

---