Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90
#include "DataFormats/ParticleFlowReco/interface/PFBlockElement.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementTrack.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementGsfTrack.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementBrem.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementCluster.h"
#include "DataFormats/ParticleFlowReco/interface/PFBlockElementSuperCluster.h"

const reco::TrackRef reco::PFBlockElement::nullTrack_ = reco::TrackRef();
const reco::PFRecTrackRef reco::PFBlockElement::nullPFRecTrack_ = reco::PFRecTrackRef();
const reco::PFClusterRef reco::PFBlockElement::nullPFCluster_ = reco::PFClusterRef();
const reco::PFDisplacedTrackerVertexRef reco::PFBlockElement::nullPFDispVertex_ = reco::PFDisplacedTrackerVertexRef();
const reco::ConversionRefVector reco::PFBlockElement::nullConv_ = reco::ConversionRefVector();
const reco::MuonRef reco::PFBlockElement::nullMuon_ = reco::MuonRef();
const reco::VertexCompositeCandidateRef reco::PFBlockElement::nullVertex_ = reco::VertexCompositeCandidateRef();

using namespace reco;

// int PFBlockElement::instanceCounter_ = 0;

// int PFBlockElement::instanceCounter() {
//   return instanceCounter_;
// }

void PFBlockElement::Dump(std::ostream& out, const char* pad) const {
  if (!out)
    return;
  out << pad << "base element";
}

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

  out << "element " << element.index() << "- type " << element.type() << " ";

  try {
    switch (element.type()) {
      case PFBlockElement::TRACK: {
        const reco::PFBlockElementTrack& et = dynamic_cast<const reco::PFBlockElementTrack&>(element);
        et.Dump(out);
        if (et.trackType(PFBlockElement::T_FROM_DISP))
          out << " from displaced;";
        if (et.trackType(PFBlockElement::T_TO_DISP))
          out << " to displaced;";
        if (et.trackType(PFBlockElement::T_FROM_GAMMACONV))
          out << " from gammaconv;";
        if (et.trackType(PFBlockElement::T_FROM_V0))
          out << " from v0 decay;";
        break;
      }
      case PFBlockElement::ECAL:
      case PFBlockElement::HCAL:
      case PFBlockElement::HGCAL:
      case PFBlockElement::HO:
      case PFBlockElement::HFEM:
      case PFBlockElement::HFHAD:
      case PFBlockElement::PS1:
      case PFBlockElement::PS2: {
        const reco::PFBlockElementCluster& ec = dynamic_cast<const reco::PFBlockElementCluster&>(element);
        ec.Dump(out);
        break;
      }
      case PFBlockElement::GSF: {
        const reco::PFBlockElementGsfTrack& eg = dynamic_cast<const reco::PFBlockElementGsfTrack&>(element);
        eg.Dump(out);
        out << " from gsf;";
        break;
      }
      case PFBlockElement::BREM: {
        const reco::PFBlockElementBrem& em = dynamic_cast<const reco::PFBlockElementBrem&>(element);
        em.Dump(out);
        out << " from brem;";
        break;
      }
      case PFBlockElement::SC: {
        const reco::PFBlockElementSuperCluster& sc = dynamic_cast<const reco::PFBlockElementSuperCluster&>(element);
        sc.Dump(out);
        out << " from SuperCluster;";
        break;
      }
      default:
        out << " unknown type" << std::endl;
        break;
    }
  } catch (std::exception& err) {
    out << err.what() << std::endl;
  }

  return out;
}