PFHcalSuperClusterAlgo.h

Go to the documentation of this file.

#ifndef RecoParticleFlow_PFClusterProducer_PFHcalSuperClusterAlgo_h
#define RecoParticleFlow_PFClusterProducer_PFHcalSuperClusterAlgo_h

#include "DataFormats/ParticleFlowReco/interface/PFCluster.h"
#include "DataFormats/ParticleFlowReco/interface/PFClusterFwd.h"
#include "DataFormats/ParticleFlowReco/interface/PFSuperCluster.h"
#include "DataFormats/ParticleFlowReco/interface/PFSuperClusterFwd.h"
#include "DataFormats/ParticleFlowReco/interface/PFLayer.h"
#include "DataFormats/ParticleFlowReco/interface/PFRecHit.h"
#include "DataFormats/Common/interface/PtrVector.h"
#include "DataFormats/Common/interface/OrphanHandle.h"

#include <vector>

#include <memory>


class PFHcalSuperClusterAlgo {

 public:

  PFHcalSuperClusterAlgo();

  virtual ~PFHcalSuperClusterAlgo() {;}

  void enableDebugging(bool debug) { debug_ = debug;}
 

  typedef edm::Handle< reco::PFClusterCollection > PFClusterHandle;
  typedef edm::Ptr< reco::PFCluster> PFClusterPtr;
  const edm::PtrVector<reco::PFCluster>& clusters() const
  {return clusters_; }
  
  void doClustering( const reco::PFClusterCollection& clusters, const reco::PFClusterCollection& clustersHO );

  std::pair<double, double> calculateWidths(const reco::PFCluster& cluster);

  std::pair<double, double> calculatePosition(const reco::PFCluster& cluster);

  void doClustering( const PFClusterHandle& clustersHandle, const PFClusterHandle& clustersHOHandle );

  // set histogram file pointer
//  void setHistos(TFile* file);

  // write histos
  void write();
  
  
 
  std::auto_ptr< std::vector< reco::PFCluster > >& clusters()  
    {return pfClusters_;}

  std::auto_ptr< std::vector< reco::PFSuperCluster > >& superClusters()  
    {return pfSuperClusters_;}


  friend std::ostream& operator<<(std::ostream& out,const PFHcalSuperClusterAlgo& algo);


 private:
  void doClusteringWorker( const reco::PFClusterCollection& clusters, const reco::PFClusterCollection& clustersHO );

  std::auto_ptr< std::vector<reco::PFCluster> > pfClusters_;
  
  std::auto_ptr< std::vector<reco::PFSuperCluster> > pfSuperClusters_;

  reco::PFSuperCluster            SuperCluster_; 
  PFClusterHandle           clustersHandle_;
  PFClusterHandle           clustersHOHandle_;
  edm::PtrVector< reco::PFCluster >  clusters_;


  bool   debug_;

  static  unsigned  prodNum_;

};

#endif