CosmicClusterProducer.h

Go to the documentation of this file.

#ifndef RecoEcal_EgammaClusterProducers_CosmicClusterProducer_h_
#define RecoEcal_EgammaClusterProducers_CosmicClusterProducer_h_
 
#include <memory>
#include <ctime>
#include <vector>  //TEMP JHAUPT 4-27
 
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/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/CosmicClusterAlgo.h"
#include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"
#include "RecoEcal/EgammaCoreTools/interface/ClusterShapeAlgo.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"
 
#include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
 
//
 
class CosmicClusterProducer : public edm::stream::EDProducer<> {
public:
  CosmicClusterProducer(const edm::ParameterSet& ps);
 
  ~CosmicClusterProducer() override;
 
  void produce(edm::Event&, const edm::EventSetup&) override;
 
private:
  int nMaxPrintout_;  // max # of printouts
  int nEvt_;          // internal counter of events
 
  CosmicClusterAlgo::VerbosityLevel verbosity;
 
  edm::EDGetTokenT<EcalRecHitCollection> ebHitsToken_;
  edm::EDGetTokenT<EcalRecHitCollection> eeHitsToken_;
 
  edm::EDGetTokenT<EcalUncalibratedRecHitCollection> ebUHitsToken_;
  edm::EDGetTokenT<EcalUncalibratedRecHitCollection> eeUHitsToken_;
 
  std::string barrelClusterCollection_;
  std::string endcapClusterCollection_;
 
  std::string clustershapecollectionEB_;
  std::string clustershapecollectionEE_;
 
  //BasicClusterShape AssociationMap
  std::string barrelClusterShapeAssociation_;
  std::string endcapClusterShapeAssociation_;
 
  PositionCalc posCalculator_;  // position calculation algorithm
  ClusterShapeAlgo shapeAlgo_;  // cluster shape algorithm
  CosmicClusterAlgo* island_p;
 
  bool counterExceeded() const { return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0)); }
 
  void clusterizeECALPart(edm::Event& evt,
                          const edm::EventSetup& es,
                          const edm::EDGetTokenT<EcalRecHitCollection>& hitsToken,
                          const edm::EDGetTokenT<EcalUncalibratedRecHitCollection>& uhitsToken,
                          const std::string& clusterCollection,
                          const std::string& clusterShapeAssociation,
                          const CosmicClusterAlgo::EcalPart& ecalPart);
 
  void outputValidationInfo(reco::CaloClusterPtrVector& clusterPtrVector);
};
 
#endif