SimpleCutsIsolator.h

Go to the documentation of this file.

#ifndef SimpleCutsIsolator_H
#define SimpleCutsIsolator_H

#include "RecoMuon/MuonIsolation/interface/MuIsoBaseIsolator.h"
#include "RecoMuon/MuonIsolation/interface/Cuts.h"

#include "FWCore/Framework/interface/ConsumesCollector.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

class SimpleCutsIsolator : public muonisolation::MuIsoBaseIsolator {
 public:
 SimpleCutsIsolator(const edm::ParameterSet & par, edm::ConsumesCollector && iC):
    theCuts(par.getParameter<std::vector<double> > ("EtaBounds"),
        par.getParameter<std::vector<double> > ("ConeSizes"),
        par.getParameter<std::vector<double> > ("Thresholds"))
    {
    }

  ResultType resultType() const override {return ISOL_BOOL_TYPE;}

  Result result(const DepositContainer& deposits, const edm::Event* = nullptr) const override {
    Result answer(ISOL_BOOL_TYPE);
    answer.valBool = false;
    // fail miserably...
    return answer;
  }

  Result result(const DepositContainer& deposits, const reco::Track& tk, const edm::Event* = nullptr) const override {
    Result answer(ISOL_BOOL_TYPE);

    muonisolation::Cuts::CutSpec cuts_here = theCuts(tk.eta());
    
    double conesize = cuts_here.conesize;
    double dephlt = 0;
    unsigned int nDeps = deposits.size();
    for(unsigned int iDep = 0; iDep < nDeps; ++iDep ){
      dephlt += deposits[iDep].dep->depositWithin(conesize);
    }
    answer.valFloat = dephlt;
    if (dephlt<cuts_here.threshold) {
      answer.valBool = true;
    } else {
      answer.valBool = false;
    }
    return answer;
  }
  
 private:

  // Isolation cuts
  muonisolation::Cuts theCuts;


};

#endif