EventDependentAbsVetos.h

Go to the documentation of this file.

#ifndef PhysicsTools_IsolationAlgos_EventDependentAbsVetos_h
#define PhysicsTools_IsolationAlgos_EventDependentAbsVetos_h

#include "PhysicsTools/IsolationAlgos/interface/EventDependentAbsVeto.h"
#include "FWCore/Utilities/interface/InputTag.h"

namespace reco {
 namespace isodeposit {
    class OtherCandidatesDeltaRVeto : public EventDependentAbsVeto {
      public:
          OtherCandidatesDeltaRVeto(const edm::InputTag& candidates, double deltaR) :
            src_(candidates), deltaR2_(deltaR*deltaR) { }

          // Virtual destructor (should always be there)
          virtual ~OtherCandidatesDeltaRVeto() {}

          virtual bool veto(double eta, double phi, float value) const ;

          virtual void centerOn(double eta, double phi) { }

          virtual void setEvent(const edm::Event &iEvent, const edm::EventSetup &iSetup) ;

      private:
          edm::InputTag src_;
          float         deltaR2_;
          std::vector<Direction> items_;
    };

    class OtherCandVeto : public EventDependentAbsVeto {
      public:
          OtherCandVeto(const edm::InputTag& candidates, AbsVeto *veto) :
            src_(candidates), veto_(veto) { }

          // Virtual destructor (should always be there)
          virtual ~OtherCandVeto() {}

          virtual bool veto(double eta, double phi, float value) const ;

          virtual void centerOn(double eta, double phi) { }

          virtual void setEvent(const edm::Event &iEvent, const edm::EventSetup &iSetup) ;

      private:
          edm::InputTag src_;
          std::vector<Direction> items_;
          std::auto_ptr<AbsVeto> veto_;
    };

    class OtherJetConstituentsDeltaRVeto : public EventDependentAbsVeto {
      public:
          OtherJetConstituentsDeltaRVeto(Direction dir, const edm::InputTag& jets, double dRjet, const edm::InputTag& pfCandAssocMap, double dRconstituent)
        : evt_(0),
          vetoDir_(dir),
          srcJets_(jets),
          dR2jet_(dRjet*dRjet),
          srcPFCandAssocMap_(pfCandAssocMap),
          dR2constituent_(dRconstituent*dRconstituent)
      {
        //std::cout << "<OtherJetConstituentsDeltaRVeto::OtherJetConstituentsDeltaRVeto>:" << std::endl;
        //std::cout << " vetoDir: eta = " << vetoDir_.eta() << ", phi = " << vetoDir_.phi() << std::endl;
        //std::cout << " srcJets = " << srcJets_.label() << ":" << srcJets_.instance() << std::endl;
        //std::cout << " dRjet = " << sqrt(dR2jet_) << std::endl;
        //std::cout << " srcPFCandAssocMap = " << srcPFCandAssocMap_.label() << ":" << srcPFCandAssocMap_.instance() << std::endl;
        //std::cout << " dRconstituent = " << sqrt(dR2constituent_) << std::endl;
      }

          // Virtual destructor (should always be there)
          virtual ~OtherJetConstituentsDeltaRVeto() {}

          virtual bool veto(double eta, double phi, float value) const;

          virtual void centerOn(double eta, double phi);

          virtual void setEvent(const edm::Event& evt, const edm::EventSetup& es);

      private:
      void initialize();

      const edm::Event* evt_;

      Direction vetoDir_;
          edm::InputTag srcJets_;
      double dR2jet_;
      edm::InputTag srcPFCandAssocMap_;
      double dR2constituent_;
          std::vector<Direction> items_;
    };
 }
}
#endif