NominalEfficiencyThresholds.h

Go to the documentation of this file.

#ifndef MuonIsolation_NominalEfficiencyThresholds_H
#define MuonIsolation_NominalEfficiencyThresholds_H

#include <map>
#include <string>
#include <utility>
#include <vector>

namespace muonisolation {
class NominalEfficiencyThresholds {
public:
  NominalEfficiencyThresholds() { }
  NominalEfficiencyThresholds(const std::string & infile);
  ~NominalEfficiencyThresholds() { }

  struct ThresholdLocation { float eta; int cone; };

  float thresholdValueForEfficiency(ThresholdLocation location, float eff_thr) const;

  std::vector<double> bins() const;
  void dump();

private:


  struct EfficiencyBin {
    float eff;
    float eff_previous;
    bool operator() (const EfficiencyBin & e1,
                 const EfficiencyBin & e2) const;
  };


  class EtaBounds {
  public:
    enum { NumberOfTowers = 32 };
    EtaBounds();
    int    towerFromEta(double eta) const;
    float  operator()(unsigned int i) const { return theBounds[i]; }
  private:
    float theBounds[NumberOfTowers+1]; //max eta of towers 1-32 (indx 1-32) and 0. for indx 0
  };


  struct locless {
    bool operator()(const ThresholdLocation & l1,
                const ThresholdLocation & l2) const;
    EtaBounds etabounds;
  };



  typedef std::pair<EfficiencyBin,float> ThresholdConstituent;
  typedef std::map<EfficiencyBin,float,EfficiencyBin> ThresholdConstituents;
  typedef std::map<ThresholdLocation,ThresholdConstituents,locless> MapType;

  void add(ThresholdLocation location, ThresholdConstituent threshold);
  MapType thresholds;

  EtaBounds etabounds;
};
}
#endif