---

BasePilupSubtractionJetProducer.h

Go to the documentation of this file.

#ifndef JetProducers_BasePilupSubtractionJetProducer_h
#define JetProducers_BasePilupSubtractionJetProducer_h

#include "FWCore/Framework/interface/EDProducer.h"
#include "DataFormats/Common/interface/EDProductfwd.h"
#include "RecoJets/JetAlgorithms/interface/JetRecoTypes.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/RecoCandidate/interface/RecoCandidate.h"
#include "DataFormats/Candidate/interface/LeafCandidate.h"
#include "DataFormats/Candidate/interface/CandidateFwd.h"
#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include <map>
#include <vector>
#include "FWCore/Framework/interface/EventSetup.h"
#include "DataFormats/HcalDetId/interface/HcalDetId.h"

namespace edm {
  class ParameterSet;
  class Event;
  class EventSetup;
}

namespace cms
{
  class BasePilupSubtractionJetProducer : public edm::EDProducer
  {
  public:

    BasePilupSubtractionJetProducer(const edm::ParameterSet& ps);

    virtual ~BasePilupSubtractionJetProducer();
    virtual void produce(edm::Event& e, const edm::EventSetup& c);
    std::string jetType () const {return mJetType;}
    
    int ieta(const reco::Candidate*);
    int iphi(const reco::Candidate*);
    
    void beginJob( const edm::EventSetup& iSetup);

    // abstract method to be set up in actual implementations
    virtual bool runAlgorithm (const  JetReco::InputCollection& fInput,  JetReco::OutputCollection* fOutput) = 0;
    
    void calculate_pedestal(const JetReco::InputCollection&);
//    reco::CandidateCollection subtract_pedestal(const JetReco::InputCollection&);
    JetReco::InputCollection subtract_pedestal(const JetReco::InputCollection&);

  private:
    edm::InputTag mSrc;
    std::string mJetType;
    bool mVerbose;
    double mEtInputCut;
    double mEInputCut;
    double mEtJetInputCut;
    double nSigmaPU;
    double radiusPU;
    std::map<int,double> esigma;
    std::map<int,double> emean;  
    std::map<int,int> geomtowers;
    std::map<int,int> ntowers_with_jets;
    std::vector<HcalDetId> allgeomid;
    const CaloGeometry* geo;
    int ietamax;
    int ietamin;
  };
}

#endif

---