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 z = 0) override;

  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