IsolatorByNominalEfficiency.h

Go to the documentation of this file.

#ifndef MuonIsolation_IsolatorByNominalEfficiency_H
#define MuonIsolation_IsolatorByNominalEfficiency_H
 
#include "RecoMuon/MuonIsolation/interface/MuIsoBaseIsolator.h"
#include "RecoMuon/MuonIsolation/interface/Cuts.h"
#include <vector>
#include <string>
 
namespace muonisolation {
  class NominalEfficiencyThresholds;
}
 
namespace muonisolation {
  class IsolatorByNominalEfficiency : public MuIsoBaseIsolator {
  public:
    typedef MuIsoBaseIsolator::DepositContainer DepositContainer;
 
    IsolatorByNominalEfficiency(const std::string& thrFile,
                                const std::vector<std::string>& ceff,
                                const std::vector<double>& weights);
 
    IsolatorByNominalEfficiency(const std::string& thrFile,
                                const std::vector<std::string>& ceff,
                                const std::vector<double>& weights,
                                const std::vector<double>& thresh);
 
    ~IsolatorByNominalEfficiency() override;
 
    Result result(const DepositContainer& deposits, const edm::Event* = nullptr) const override;
 
    Cuts cuts(float nominalEfficiency) const;
 
  private:
    class ConeSizes {
    private:
      enum IsoDim { DIM = 15 };
      static const float cone_dr[DIM];
 
    public:
      int dim() const { return DIM; }
      double size(int i) const;
      int index(float dr) const;
    };
 
    // Compute the weighted sum of deposits of different type within dRcone
    virtual double weightedSum(const DepositContainer& deposits, float dRcone) const;
 
    // Size of cone for a given nominal efficiency value.
    int bestConeForEfficiencyIndex(float eff_thr) const;
 
    typedef std::multimap<float, int> mapNomEff_Cone;
    mapNomEff_Cone cones(const std::vector<std::string>& names);
 
    std::string findPath(const std::string& fileName);
 
    ResultType resultType() const override { return ISOL_FLOAT_TYPE; }
 
  private:
    mapNomEff_Cone coneForEfficiency;
    NominalEfficiencyThresholds* thresholds;
    std::vector<double> theWeights;
    std::vector<double> theDepThresholds;
    ConeSizes theConesInfo;
  };
}  // namespace muonisolation
 
#endif