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