HiFJRhoProducer.h

Go to the documentation of this file.

#ifndef HiJetBackground_HiFJRhoProducer_h
#define HiJetBackground_HiFJRhoProducer_h
 
// system include files
#include <memory>
#include <sstream>
#include <string>
#include <vector>
 
// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
 
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/StreamID.h"
 
#include "DataFormats/JetReco/interface/Jet.h"
 
//
// class declaration
//
 
class HiFJRhoProducer : public edm::stream::EDProducer<> {
public:
  explicit HiFJRhoProducer(const edm::ParameterSet&);
  ~HiFJRhoProducer() override;
 
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
 
private:
  void beginStream(edm::StreamID) override;
  void produce(edm::Event&, const edm::EventSetup&) override;
  void endStream() override;
 
  double calcMedian(std::vector<double>& v);
  double calcMd(const reco::Jet* jet);
  bool isPackedCandidate(const reco::Candidate* candidate);
 
  // ----------member data ---------------------------
  //input
  edm::EDGetTokenT<edm::View<reco::Jet>> jetsToken_;
 
  //members
  edm::InputTag src_;             // input kt jet source
  unsigned int nExcl_;            //Number of leading jets to exclude
  double etaMaxExcl_;             //max eta for jets to exclude
  double ptMinExcl_;              //min pt for excluded jets
  unsigned int nExcl2_;           //Number of leading jets to exclude in 2nd eta region
  double etaMaxExcl2_;            //max eta for jets to exclude in 2nd eta region
  double ptMinExcl2_;             //min pt for excluded jets in 2nd eta region
  std::vector<double> etaRanges;  //eta boundaries for rho calculation regions
  bool checkJetCand, usingPackedCand;
};
 
#endif