PFRecHitProducerECAL.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterProducer_PFRecHitProducerECAL_h_
#define RecoParticleFlow_PFClusterProducer_PFRecHitProducerECAL_h_

// system include files
#include <memory>
#include <vector>

// user include files
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"

#include "FWCore/ParameterSet/interface/ParameterSet.h"

#include "Geometry/CaloTopology/interface/CaloDirection.h"

#include "RecoParticleFlow/PFClusterProducer/plugins/PFRecHitProducer.h"
#include "DataFormats/ParticleFlowReco/interface/PFLayer.h"

#include "DataFormats/Math/interface/Vector3D.h"
class CaloSubdetectorTopology;
class CaloSubdetectorGeometry;
class EcalBarrelGeometry;
class EcalEndcapGeometry;
class CaloSubdetectorGeometry;
class DetId;





class PFRecHitProducerECAL : public PFRecHitProducer {

 public:
  explicit PFRecHitProducerECAL(const edm::ParameterSet&);
  ~PFRecHitProducerECAL();

  

 private:

  void createRecHits(std::vector<reco::PFRecHit>& rechits,
             std::vector<reco::PFRecHit>& rechitsCleaned,
             edm::Event&, const edm::EventSetup&);


  

  reco::PFRecHit*  createEcalRecHit( const DetId& detid,
                     double energy,
                     PFLayer::Layer layer,
                     const CaloSubdetectorGeometry* geom );



  bool findEcalRecHitGeometry ( const DetId& detid, 
                const CaloSubdetectorGeometry* geom,
                math::XYZVector& position, 
                math::XYZVector& axis );

  void 
    findRecHitNeighboursECAL( reco::PFRecHit& rh, 
                  const std::map<unsigned,unsigned >& sortedHits );

  void ecalNeighbArray( const EcalBarrelGeometry& barrelGeom,
            const CaloSubdetectorTopology& barrelTopo,
            const EcalEndcapGeometry& endcapGeom,
            const CaloSubdetectorTopology& endcapTopo );

  DetId move(DetId cell, const CaloDirection& dir ) const;

  bool stdsimplemove(DetId& cell, 
             const CaloDirection& dir,
             const CaloSubdetectorTopology& barrelTopo,
             const CaloSubdetectorTopology& endcapTopo,
             const EcalBarrelGeometry& barrelGeom,
             const EcalEndcapGeometry& endcapGeom ) const;


  bool stdmove(DetId& cell, 
           const CaloDirection& dir,
           const CaloSubdetectorTopology& barrelTopo,
           const CaloSubdetectorTopology& endcapTopo,
           const EcalBarrelGeometry& barrelGeom,
           const EcalEndcapGeometry& endcapGeom ) const;

  // ----------member data ---------------------------
  

 
  std::vector<std::vector<DetId> >  neighboursEB_;

  std::vector<std::vector<DetId> >  neighboursEE_;
  
  bool  neighbourmapcalculated_;

  bool  crossBarrelEndcapBorder_;

  // ----------access to event data
  edm::InputTag    inputTagEcalRecHitsEB_;
  edm::InputTag    inputTagEcalRecHitsEE_;

  // ---- Perform timing cleaning
  bool timingCleaning_;
  double threshCleaningEB_;
  double threshCleaningEE_;

  // ---- Perform topological cleaning
  bool topologicalCleaning_;

};

#endif