/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_6_2_5/src/DataFormats/ParticleFlowReco/src/PFSuperCluster.cc

Go to the documentation of this file.

#include "DataFormats/ParticleFlowReco/interface/PFSuperCluster.h"
#include "DataFormats/GeometryVector/interface/Pi.h"
#include "DataFormats/Common/interface/PtrVector.h"


using namespace std;
using namespace reco;

PFSuperCluster::PFSuperCluster(const edm::PtrVector<reco::PFCluster>& clusters):
clusters_(clusters)
{  
}

void PFSuperCluster::reset() {
  
  PFCluster::reset();
  clusters_.clear();

}

PFSuperCluster& PFSuperCluster::operator=(const PFSuperCluster& other) {

  PFCluster::operator=((PFCluster)other); 
  clusters_ = other.clusters_;

  return *this;
}


std::ostream& reco::operator<<(std::ostream& out, 
                               const PFSuperCluster& cluster) {
  
  if(!out) return out;

  const math::XYZPoint&  pos = cluster.position();
  const PFCluster::REPPoint&  posrep = cluster.positionREP();
  const std::vector< reco::PFRecHitFraction >& fracs =
    cluster.recHitFractions();

  out<<"PFSuperCluster "
     <<", clusters: "<<cluster.clusters().size()
     <<", layer: "<<cluster.layer()
     <<"\tE = "<<cluster.energy()
     <<"\tXYZ: "
     <<pos.X()<<","<<pos.Y()<<","<<pos.Z()<<" | "
     <<"\tREP: "
     <<posrep.Rho()<<","<<posrep.Eta()<<","<<posrep.Phi()<<" | "
     <<fracs.size()<<" rechits";

  for(unsigned i=0; i<fracs.size(); i++) {
    // PFRecHit is not available, print the detID
    if( !fracs[i].recHitRef().isAvailable() )
      out<<cluster.printHitAndFraction(i)<<", ";
    else
      out<<fracs[i]<<", ";
  }

  
  return out;
}