/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/RecoParticleFlow/PFClusterProducer/plugins/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