---

Track.h

Go to the documentation of this file.

// ================================================
// Track class description
// ================================================
//
// THIS TEXT TO BE REPLACED BY ATLAS STANDARD FORMAT
//
// Namespace Atlfast::
//
// class: Track
//
// Authors: P.Clarke, H.Phillips, E.Richter-Was, P.Sherwood, R.Steward
//
// Description: 
//
//  Class to represent Tracks 
//
//  NOT included yet are:
//
//    1. Up to 6 mothers ? as I didnt understand how there could be more
//       than 1 ?
//    
//    2. Correlation coefficients. The thought is that if these go in they
//       go into ITrackTrajectory. Need to understand their use more.
//
//..................................................


#ifndef ATLFAST_TRACK_H
#define ATLFAST_TRACK_H


// CLHEP, HepMC
#include "CLHEP/Vector/ThreeVector.h"
#include "HepMC/GenParticle.h"
#include "CLHEP/Matrix/Matrix.h"

// Atlfast classes.
#include "AtlfastEvent/TrackTrajectory.h"
#include "AtlfastEvent/AssociationManager.h"

// Gaudi
#include "GaudiKernel/ContainedObject.h"
#include "GaudiKernel/MsgStream.h"


//Class ID for Track class
static const CLID CLID_ATLFAST_Track=2306 ;

using ::HepMatrix;

namespace Atlfast {


  class Track : 
    virtual public ContainedObject, 
    public AssociationManager
    {
      
    private:

      TrackTrajectory m_trajectory;
      TrackTrajectory m_truthTrajectory;
      const HepMC::GenParticle* m_truth ;
      HepMatrix m_smearMatrix;
      int m_pdg_id;

    public:
      
      // -----------------------------------
      // Constructors

      Track( const TrackTrajectory& trajectory, 
             const HepMC::GenParticle* progenitor  ) 
        : ContainedObject(),
        AssociationManager(),
        m_trajectory(trajectory),
        m_truthTrajectory(trajectory),
        m_truth( progenitor ),
        m_smearMatrix(HepMatrix(5,5,0.)),
        m_pdg_id(progenitor->pdg_id())
        { }
      
      Track( const TrackTrajectory& truthTrajectory, 
             const TrackTrajectory& smearTrajectory,
             const HepMC::GenParticle* progenitor,
             const HepMatrix& matrix  ) 
        : ContainedObject(),
        AssociationManager(),
        m_trajectory(smearTrajectory),
        m_truthTrajectory(truthTrajectory),
        m_truth( progenitor ),
        m_smearMatrix(matrix),
        m_pdg_id(progenitor->pdg_id())
        { }

      Track( const TrackTrajectory& truthTrajectory, 
             const TrackTrajectory& smearTrajectory,
             const HepMatrix& matrix, int pdg  ) 
        : ContainedObject(),
        AssociationManager(),
        m_trajectory(smearTrajectory),
        m_truthTrajectory(truthTrajectory),
        m_truth( NULL ),
        m_smearMatrix(matrix),
        m_pdg_id(pdg)
        { }
      
      // Copy
      Track( const Track& other  ) : ContainedObject( ),
        AssociationManager(other),
        m_trajectory(other.m_trajectory),
        m_truthTrajectory(other.m_truthTrajectory),
        m_truth( other.m_truth ),
        m_smearMatrix(other.m_smearMatrix),
        m_pdg_id(other.m_pdg_id)
        { }
      
      // Default which seems to be required by Athena, although doesnt have any 
      // meaning for the Track class. Hence a bogus trajectory has to be invented.
      Track() 
        : ContainedObject(), 
        AssociationManager(),
        m_truth(NULL)
        { }
      
      ~Track() {}
      
      virtual TrackTrajectory trajectory() const {return m_trajectory;}
      virtual int pdg_id() const {return m_pdg_id;}
      const HepMC::GenParticle* truth() const { return m_truth ; }
      const HepMatrix smearMatrix() const {return m_smearMatrix;}
      void setMatrix(const HepMatrix& matrix) {m_smearMatrix = matrix;}
      const TrackTrajectory truthTrajectory() const {return m_truthTrajectory;}
      // --------------------------------
      // Athena mandatory methods
      static const CLID& classID() { return CLID_ATLFAST_Track; }
      virtual const CLID& clID() const {  return CLID_ATLFAST_Track; }
      
    };
  
  
  // ................................................
  // Message stream operator for simple debug output
  
  MsgStream& operator << ( MsgStream& s, const Track& t ) ;
  MsgStream& operator << ( MsgStream& s, const Track* t ) ;
  
  
} // end of namespace 


#endif

---