TrackResiduals.h

Go to the documentation of this file.

#ifndef TrackReco_TrackResiduals_h
#define TrackReco_TrackResiduals_h

#include <iostream>
#include "DataFormats/TrackingRecHit/interface/TrackingRecHitFwd.h"
#include "DataFormats/TrackReco/interface/TrackResiduals.h"

class Trajectory;

namespace reco {
     class HitPattern;
     class TrackResiduals {
      friend class Trajectory;
     public:
      enum ResidualType { X_Y_RESIDUALS, X_Y_PULLS };
      
     protected:
      enum { numResiduals = 0x40 };
      
     public:
      TrackResiduals ();
      TrackResiduals (enum ResidualType);
      
     public:
      void setResidualXY (int idx, double residualX, double residualY);
      void setPullXY (int idx, double pullX, double pullY);
      void setResidualType (enum ResidualType);
      void print (std::ostream &stream = std::cout) const;
      void print (const HitPattern &, std::ostream &stream = std::cout) const;
      /* Example usage:
       * 
       * HitPattern p;
       * TrackResiduals r;
       * int numHits = p.numberOfValidHits() +        
       *   p.numberOfLostHits() +         
       *   p.numberOfValidMuonHits() +    
       *   p.numberOfLostMuonHits() +     
       *   p.numberOfValidTrackerHits() + 
       *   p.numberOfLostTrackerHits() +  
       *   p.numberOfValidPixelHits() +   
       *   p.numberOfLostPixelHits();
       * for (int i = 0; i < numHits; ++i) {
       *   printf("hit %d: pattern 0x%x, residual (local x) %f\n", i
       *         p.getHitPattern(i), r.residualX(p, i));
       * }
       */
      double        residualX (int i, const HitPattern &) const;
      double        residualY (int i, const HitPattern &) const;
      double        residualX (int i) const;
      double        residualY (int i) const;
       
     protected:
      static double     unpack_pull (unsigned char);
      static unsigned char      pack_pull (double);
      static double     unpack_residual (unsigned char);
      static unsigned char      pack_residual (double);
       
     private:
      unsigned char     residuals_[numResiduals];
      char          residualType;
     };
}

#endif