/data/refman/pasoursint/CMSSW_5_3_3/src/Validation/RecoMuon/src/HTrack.h

Go to the documentation of this file.

class HTrackVariables;  
class HResolution;
class TFile;
class SimTrack;
class TrajectoryStateOnSurface;
class FreeTrajectoryState;
class DQMStore;

#include "TString.h" 
#include <string>

class HTrack{

 public: 
  HTrack(std::string, std::string name, std::string whereIs ="");
  
  
  double pull(double rec,double sim, double sigmarec);
  
  double resolution(double rec,double sim);
  
  void computeResolution(FreeTrajectoryState& fts,
                         SimTrack &simTracks,
                         HResolution* hReso);
  
  void computeTDRResolution(FreeTrajectoryState& fts,
                            SimTrack &simTracks,
                            HResolution* hReso);
  
  void computePull(FreeTrajectoryState& fts,
                   SimTrack &simTracks,
                   HResolution* hReso); 


  void computeResolutionAndPull(TrajectoryStateOnSurface &vtx,
                                SimTrack &simTrack);

  void computeResolutionAndPull(FreeTrajectoryState& fts, 
                                SimTrack& simTrack);
  
  void Fill(TrajectoryStateOnSurface &);
  void Fill(FreeTrajectoryState &);
  void FillDeltaR(double);

  double computeEfficiency(HTrackVariables *sim);

 private:
  DQMStore * dbe_;

  HTrackVariables *hVariables;  
  
  // Global Resolution
  HResolution *hResolution;  
  HResolution *hPull;

  // TDR Resolution
  HResolution *hTDRResolution; 
  HResolution *hTDRPull;

  // [5-10] GeV range
  HResolution *hResolution_5_10;
  HResolution *hTDRResolution_5_10;
  HResolution *hPull_5_10;

  // [10-40] GeV range
  HResolution *hResolution_10_40; 
  HResolution *hTDRResolution_10_40;
  HResolution *hPull_10_40; 

  // [40-70] GeV range
  HResolution *hResolution_40_70; 
  HResolution *hTDRResolution_40_70;
  HResolution *hPull_40_70; 
  
  // [70-100] GeV range
  HResolution *hResolution_70_100;
  HResolution *hTDRResolution_70_100;
  HResolution *hPull_70_100;

  // eta range |eta|<0.8
  HResolution *hResolution_08;
  HResolution *hTDRResolution_08;
  HResolution *hPull_08;
  
  // eta range 0.8<|eta|<1.2
  HResolution *hResolution_08_12;
  HResolution *hTDRResolution_08_12;
  HResolution *hPull_08_12;
  
  // eta range 1.2<|eta|<2.1
  HResolution *hResolution_12_21;
  HResolution *hTDRResolution_12_21;
  HResolution *hPull_12_21;

  // eta range 1.2<|eta|<2.4
  HResolution *hResolution_12_24;
  HResolution *hTDRResolution_12_24;
  HResolution *hPull_12_24;

  // eta range 1.2<eta<2.1
  HResolution *hResolution_12_21_plus;
  HResolution *hTDRResolution_12_21_plus;
  HResolution *hPull_12_21_plus;

  // eta range 1.2<eta<2.4
  HResolution *hResolution_12_24_plus;
  HResolution *hTDRResolution_12_24_plus;
  HResolution *hPull_12_24_plus;

  // eta range -2.1<eta<-1.2
  HResolution *hResolution_12_21_minus;
  HResolution *hTDRResolution_12_21_minus;
  HResolution *hPull_12_21_minus;

  // eta range -2.4<eta<-1.2
  HResolution *hResolution_12_24_minus;
  HResolution *hTDRResolution_12_24_minus;
  HResolution *hPull_12_24_minus;

  TString theName;
  TString where;

  bool doSubHisto;
};