00001 #ifndef RecoEcal_EgammaClusterProducers_EgammaHLTMulti5x5ClusterProducer_h_ 00002 #define RecoEcal_EgammaClusterProducers_EgammaHLTMulti5x5ClusterProducer_h_ 00003 00004 #include <memory> 00005 #include <time.h> 00006 00007 #include "FWCore/Framework/interface/Frameworkfwd.h" 00008 #include "FWCore/Framework/interface/EDProducer.h" 00009 #include "FWCore/Framework/interface/Event.h" 00010 #include "FWCore/Framework/interface/EventSetup.h" 00011 00012 #include "FWCore/ParameterSet/interface/ParameterSet.h" 00013 00014 #include "DataFormats/EcalRecHit/interface/EcalRecHitCollections.h" 00015 #include "RecoEcal/EgammaClusterAlgos/interface/Multi5x5ClusterAlgo.h" 00016 #include "RecoEcal/EgammaCoreTools/interface/PositionCalc.h" 00017 #include "RecoEcal/EgammaCoreTools/interface/ClusterShapeAlgo.h" 00018 00019 #include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h" 00020 // 00021 00022 00023 class EgammaHLTMulti5x5ClusterProducer : public edm::EDProducer 00024 { 00025 public: 00026 00027 EgammaHLTMulti5x5ClusterProducer(const edm::ParameterSet& ps); 00028 00029 ~EgammaHLTMulti5x5ClusterProducer(); 00030 00031 virtual void produce(edm::Event&, const edm::EventSetup&); 00032 00033 private: 00034 00035 int nMaxPrintout_; // max # of printouts 00036 int nEvt_; // internal counter of events 00037 00038 Multi5x5ClusterAlgo::VerbosityLevel verbosity; 00039 00040 bool doBarrel_; 00041 bool doEndcaps_; 00042 bool doIsolated_; 00043 00044 edm::InputTag barrelHitProducer_;//it was string 00045 edm::InputTag endcapHitProducer_;//it was string 00046 std::string barrelHitCollection_; 00047 std::string endcapHitCollection_; 00048 00049 std::string barrelClusterCollection_; 00050 std::string endcapClusterCollection_; 00051 00052 /* 00053 std::string clustershapecollectionEB_;//new 00054 std::string clustershapecollectionEE_;//new 00055 //BasicClusterShape AssociationMap 00056 std::string barrelClusterShapeAssociation_;//new 00057 std::string endcapClusterShapeAssociation_;//new 00058 */ 00059 00060 edm::InputTag l1TagIsolated_; 00061 edm::InputTag l1TagNonIsolated_; 00062 double l1LowerThr_; 00063 double l1UpperThr_; 00064 double l1LowerThrIgnoreIsolation_; 00065 00066 double regionEtaMargin_; 00067 double regionPhiMargin_; 00068 00069 PositionCalc posCalculator_; // position calculation algorithm 00070 Multi5x5ClusterAlgo * Multi5x5_p; 00071 //ClusterShapeAlgo shapeAlgo_;//new 00072 00073 bool counterExceeded() const { return ((nEvt_ > nMaxPrintout_) || (nMaxPrintout_ < 0)); } 00074 00075 const EcalRecHitCollection * getCollection(edm::Event& evt, 00076 const std::string& hitProducer_, 00077 const std::string& hitCollection_); 00078 00079 00080 void clusterizeECALPart(edm::Event &evt, const edm::EventSetup &es, 00081 const std::string& hitProducer, 00082 const std::string& hitCollection, 00083 const std::string& clusterCollection, 00084 const std::vector<EcalEtaPhiRegion>& regions, 00085 const Multi5x5ClusterAlgo::EcalPart& ecalPart); 00086 /* 00087 void clusterizeECALPart(edm::Event &evt, const edm::EventSetup &es, 00088 const std::string& hitProducer, 00089 const std::string& hitCollection, 00090 const std::string& clusterCollection, 00091 const std::string& clusterShapeAssociation, 00092 const Multi5x5ClusterAlgo::EcalPart& ecalPart); 00093 */ 00094 void outputValidationInfo(reco::BasicClusterRefVector &clusterRefVector); 00095 00096 }; 00097 00098 00099 #endif