/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/RecoEgamma/EgammaHLTProducers/interface/EgammaHLTNxNClusterProducer.h

Go to the documentation of this file.

#ifndef RecoEcal_EgammaClusterProducers_EgammaHLTNxNClusterProducer_h_
#define RecoEcal_EgammaClusterProducers_EgammaHLTNxNClusterProducer_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 "RecoEcal/EgammaCoreTools/interface/PositionCalc.h"
#include "Geometry/CaloTopology/interface/CaloSubdetectorTopology.h"

//
#include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
#include "DataFormats/Math/interface/Point3D.h"
#include "DataFormats/CaloRecHit/interface/CaloID.h"

//Ecal status
#include "CondFormats/DataRecord/interface/EcalChannelStatusRcd.h"
#include "CondFormats/EcalObjects/interface/EcalChannelStatus.h"

#include "TVector3.h"
#include <vector>



// Less than operator for sorting EcalRecHits according to energy.
class ecalRecHitSort : public std::binary_function<EcalRecHit, EcalRecHit, bool> 
{
 public:
  bool operator()(EcalRecHit x, EcalRecHit y) 
    { 
      return (x.energy() > y.energy()); 
    }
};




class EgammaHLTNxNClusterProducer : public edm::EDProducer 
{
  public:

      EgammaHLTNxNClusterProducer(const edm::ParameterSet& ps);

      ~EgammaHLTNxNClusterProducer();
      
     
      
      virtual void produce(edm::Event&, const edm::EventSetup&);
      
   private:
      
      void makeNxNClusters(edm::Event &evt, const edm::EventSetup &es,const EcalRecHitCollection *hits, const reco::CaloID::Detectors detector); 
      
      bool checkStatusOfEcalRecHit(const EcalChannelStatus &channelStatus, const EcalRecHit &rh);
      
      
      std::string barrelClusterCollection_;
      std::string endcapClusterCollection_;
      

      std::string barrelHits_;
      std::string endcapHits_;
      
      
      
      PositionCalc posCalculator_; // position calculation algorithm
      
      

      std::map<std::string,double> providedParameters;
      
      
   
      edm::InputTag barrelHitProducer_;
      edm::InputTag endcapHitProducer_;
      
      
      double clusSeedThr_;
      double clusSeedThrEndCap_;
      
      bool doBarrel_;
      bool doEndcaps_;
      
      
      bool useRecoFlag_; 
      bool useDBStatus_; 
      int flagLevelRecHitsToUse_; 
      int statusLevelRecHitsToUse_;

      int clusEtaSize_ ;
      int clusPhiSize_;
      
      
      int debug_; 
      
      
      int maxNumberofSeeds_ ; 
      int maxNumberofClusters_; 
      


      
};


#endif