Go to the documentation of this file.00001 #ifndef RECOLOCALCALO_HFCLUSTERPRODUCER_HFCLUSTERALGO_H
00002 #define RECOLOCALCALO_HFCLUSTERPRODUCER_HFCLUSTERALGO_H 1
00003
00004 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
00005 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
00006 #include "DataFormats/EgammaReco/interface/HFEMClusterShape.h"
00007 #include "DataFormats/EgammaReco/interface/BasicCluster.h"
00008 #include "DataFormats/EgammaReco/interface/SuperCluster.h"
00009 #include "DataFormats/EgammaReco/interface/HFEMClusterShapeAssociation.h"
00010 #include "DataFormats/EgammaReco/interface/HFEMClusterShape.h"
00011 #include "DataFormats/EgammaReco/interface/BasicClusterFwd.h"
00012 #include "DataFormats/EgammaReco/interface/SuperClusterFwd.h"
00013 #include <map>
00014 #include <list>
00015
00020
00021
00022 class HFClusterAlgo {
00023 public:
00024 HFClusterAlgo();
00025
00026 void setup(double minTowerEnergy, double seedThreshold,double maximumSL,double m_maximumRenergy,bool usePMTflag,bool usePulseflag, bool forcePulseFlagMC, int correctionSet);
00027
00028 void isMC(bool isMC) { m_isMC=isMC; }
00029
00031 void clusterize(const HFRecHitCollection& hf,
00032 const CaloGeometry& geom,
00033 reco::HFEMClusterShapeCollection& clusters,
00034 reco::SuperClusterCollection& SuperClusters);
00035
00036
00037 private:
00038 friend class CompareHFCompleteHitET;
00039 friend class CompareHFCore;
00040
00041 double m_minTowerEnergy, m_seedThreshold,m_maximumSL,m_maximumRenergy;
00042 bool m_usePMTFlag;
00043 bool m_usePulseFlag,m_forcePulseFlagMC;
00044 bool m_isMC;
00045 int m_correctionSet;
00046 std::vector<double> m_cutByEta;
00047 std::vector<double> m_correctionByEta;
00048
00049 struct HFCompleteHit {
00050 HcalDetId id;
00051 double energy, et;
00052 };
00053 bool isPMTHit(const HFRecHit& hfr);
00054 bool makeCluster(const HcalDetId& seedid,
00055 const HFRecHitCollection& hf,
00056 const CaloGeometry& geom,
00057 reco::HFEMClusterShape& clusShp,
00058 reco::SuperCluster& SClus);
00059 };
00060
00061 #endif