HiFJGridEmptyAreaCalculator.h

Go to the documentation of this file.

#ifndef HiJetBackground_HiFJGridEmptyAreaCalculator_h
#define HiJetBackground_HiFJGridEmptyAreaCalculator_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/Common/interface/Handle.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "DataFormats/PatCandidates/interface/Jet.h"

class HiFJGridEmptyAreaCalculator : public edm::stream::EDProducer<> {
public:
  explicit HiFJGridEmptyAreaCalculator(const edm::ParameterSet&);
  ~HiFJGridEmptyAreaCalculator() override;

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

private:
  void beginStream(edm::StreamID) override;
  void produce(edm::Event&, const edm::EventSetup&) override;
  void endStream() override;

  //\{
  //----------------------------------------------------------------


private:
  void setupGrid(double eta_min, double eta_max);
  void setupGridJet(const reco::Jet* jet);

  int tileIndexJet(const reco::Candidate* pfCand);
  int tileIndexEta(const reco::Candidate* pfCand);
  int tileIndexEtaJet(const reco::Candidate* pfCand);
  int tileIndexPhi(const reco::Candidate* pfCand);

  int numJetGridCells(std::vector<std::pair<int, int>>& indices);

  void calculateAreaFractionOfJets(const edm::Event& iEvent, const edm::EventSetup& iSetup);
  void calculateGridRho(const edm::Event& iEvent, const edm::EventSetup& iSetup);

  const double twopi_ = 2 * M_PI;


  //parameters of grid covering the full acceptance
  double ymin_;
  double ymax_;
  double dy_;
  double dphi_;
  double tileArea_;

  //parameters of grid around jets
  double dyJet_;
  double yminJet_;
  double ymaxJet_;
  double totalInboundArea_;

  //leave bands at boundaries
  double etaminJet_;
  double etamaxJet_;

  //for the grid calculation covering the full acceptance
  int ny_;
  int nphi_;
  int ntotal_;

  //for the grid calculation around each jet
  int ntotalJet_;
  int nyJet_;

  double gridWidth_;
  double band_;
  int hiBinCut_;
  bool doCentrality_;
  bool keepGridInfo_;

  std::vector<double> rhoVsEta_;
  std::vector<double> meanRhoVsEta_;
  std::vector<double> etaMaxGrid_;
  std::vector<double> etaMinGrid_;

  int n_tiles() { return ntotal_; }

  edm::EDGetTokenT<edm::View<reco::Jet>> jetsToken_;
  edm::EDGetTokenT<reco::CandidateView> pfCandsToken_;
  edm::EDGetTokenT<std::vector<double>> mapEtaToken_;
  edm::EDGetTokenT<std::vector<double>> mapRhoToken_;
  edm::EDGetTokenT<std::vector<double>> mapRhoMToken_;

  edm::EDGetTokenT<int> centralityBinToken_;
};

#endif