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