/data/refman/pasoursint/CMSSW_5_3_0/src/DataFormats/TrackReco/interface/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