EgammaHLTMulti5x5ClusterProducer.h

Go to the documentation of this file.

#ifndef RecoEcal_EgammaClusterProducers_EgammaHLTMulti5x5ClusterProducer_h_
#define RecoEcal_EgammaClusterProducers_EgammaHLTMulti5x5ClusterProducer_h_

#include <memory>
#include <time.h>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"
#include "DataFormats/CaloRecHit/interface/CaloClusterFwd.h"
#include "RecoEcal/EgammaClusterAlgos/interface/Multi5x5ClusterAlgo.h"

#include "DataFormats/L1Trigger/interface/L1EmParticle.h"
#include "DataFormats/L1Trigger/interface/L1EmParticleFwd.h"

#include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h"

namespace edm {
  class ConfigurationDescriptions;
}

class EgammaHLTMulti5x5ClusterProducer : public edm::EDProducer {
 public:
  EgammaHLTMulti5x5ClusterProducer(const edm::ParameterSet& ps);
  ~EgammaHLTMulti5x5ClusterProducer();
  virtual void produce(edm::Event&, const edm::EventSetup&);
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);

 private:
  int nMaxPrintout_; // max # of printouts
  int nEvt_;         // internal counter of events

  bool doBarrel_;
  bool doEndcaps_;
  bool doIsolated_;

  edm::EDGetTokenT<l1extra::L1EmParticleCollection> l1TagIsolated_;
  edm::EDGetTokenT<l1extra::L1EmParticleCollection> l1TagNonIsolated_;
  edm::InputTag barrelHitCollection_;
  edm::InputTag endcapHitCollection_;
  edm::EDGetTokenT<EcalRecHitCollection> barrelHitToken_;
  edm::EDGetTokenT<EcalRecHitCollection> endcapHitToken_;
  
  std::string barrelClusterCollection_;
  std::string endcapClusterCollection_;

  double l1LowerThr_;
  double l1UpperThr_;
  double l1LowerThrIgnoreIsolation_;
  
  double regionEtaMargin_;
  double regionPhiMargin_;
  
  PositionCalc posCalculator_; // position calculation algorithm
  Multi5x5ClusterAlgo * Multi5x5_p;
    
  bool counterExceeded() const { return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0)); }

  const EcalRecHitCollection * getCollection(edm::Event& evt,
                         edm::EDGetTokenT<EcalRecHitCollection>& hitToken);
  
  
  void clusterizeECALPart(edm::Event &evt, const edm::EventSetup &es,
              edm::EDGetTokenT<EcalRecHitCollection>& hitToken,
              const std::string& clusterCollection,
              const std::vector<EcalEtaPhiRegion>& regions,
              const reco::CaloID::Detectors detector);

  void outputValidationInfo(reco::CaloClusterPtrVector &clusterPtrVector);

};
#endif