ParticleTowerProducer.h

Go to the documentation of this file.

#ifndef ParticleTowerProducer_h
#define ParticleTowerProducer_h

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/Framework/interface/ESHandle.h"

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

#include "DataFormats/HcalDetId/interface/HcalDetId.h"
#include "DataFormats/ParticleFlowCandidate/interface/PFCandidate.h"
#include "DataFormats/CaloTowers/interface/CaloTowerCollection.h"

#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
#include "Geometry/Records/interface/CaloGeometryRecord.h"

#include "TMath.h"
#include "TRandom.h"

class ParticleTowerProducer : public edm::EDProducer {
public:
  explicit ParticleTowerProducer(const edm::ParameterSet&);
  ~ParticleTowerProducer() override;

private:
  void beginJob() override;
  void produce(edm::Event&, const edm::EventSetup&) override;
  void endJob() override;
  void resetTowers(edm::Event& iEvent, const edm::EventSetup& iSetup);
  DetId getNearestTower(const reco::PFCandidate& in) const;
  DetId getNearestTower(double eta, double phi) const;
  //  uint32_t denseIndex(int ieta, int iphi, double eta) const;
  int eta2ieta(double eta) const;
  int phi2iphi(double phi, int ieta) const;

  // ----------member data ---------------------------

  edm::EDGetTokenT<reco::PFCandidateCollection> src_;
  bool useHF_;

  std::map<DetId, double> towers_;

  double PI;
  TRandom* random_;

  CaloGeometry const* geo_;  // geometry

  static const double etatow[];
  static const double etacent[];
  double etaedge[42];
};
#endif