/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_7/src/RecoParticleFlow/PFClusterTools/interface/PFResolutionMap.h

Go to the documentation of this file.

#ifndef DataFormats_ParticleFlowReco_PFResolutionMap_h
#define DataFormats_ParticleFlowReco_PFResolutionMap_h

#include <iostream>
#include <string>
#include <stdexcept>

#include <TH2.h>



class PFResolutionMap : public TH2D {

 public:

  PFResolutionMap() : TH2D() {}
  
  PFResolutionMap(const char* name, const char* mapfile);
  
  PFResolutionMap(const char* name, 
                  unsigned nbinseta, double mineta, double maxeta,
                  unsigned nbinse, double mine, double maxe, double value=-1);
  
  PFResolutionMap(const TH2D& h) : TH2D(h) {}
 

  bool ReadMapFile(const char* mapfile);

  bool WriteMapFile(const char* mapfile);

  int  FindBin(double eta, double e);

  double getRes(double eta, double phi, double e,int MapEta = -1); 

  const char* GetMapFile() const {return mapFile_.c_str();}

  friend std::ostream& operator<<(std::ostream& out, const PFResolutionMap& rm);
  
 private:
  bool IsInAPhiCrack(double phi, double eta);
  double minimum(double a,double b);
  double dCrackPhi(double phi, double eta);
  static const unsigned lineSize_;
  std::string           mapFile_;
};

#endif