FlavorHistory.cc

Go to the documentation of this file.

#include "DataFormats/HepMCCandidate/interface/FlavorHistory.h"

using namespace reco;
using namespace edm;
using namespace std;

FlavorHistory::FlavorHistory() { flavorSource_ = FLAVOR_NULL; }

FlavorHistory::FlavorHistory(FLAVOR_T flavorSource,
                             CandidatePtr const& parton,
                             CandidatePtr const& progenitor,
                             CandidatePtr const& sister,
                             reco::ShallowClonePtrCandidate const& matchedJet,
                             reco::ShallowClonePtrCandidate const& sisterJet)
    : flavorSource_(flavorSource),
      parton_(parton),
      progenitor_(progenitor),
      sister_(sister),
      matchedJet_(matchedJet),
      sisterJet_(sisterJet) {}

FlavorHistory::FlavorHistory(FLAVOR_T flavorSource,
                             Handle<View<Candidate> > h_partons,
                             int parton,
                             int progenitor,
                             int sister,
                             reco::ShallowClonePtrCandidate const& matchedJet,
                             reco::ShallowClonePtrCandidate const& sisterJet)
    : flavorSource_(flavorSource),
      parton_(parton >= 0 && static_cast<unsigned int>(parton) < h_partons->size() ? CandidatePtr(h_partons, parton)
                                                                                   : CandidatePtr()),
      progenitor_(progenitor >= 0 && static_cast<unsigned int>(progenitor) < h_partons->size()
                      ? CandidatePtr(h_partons, progenitor)
                      : CandidatePtr()),
      sister_(sister >= 0 && static_cast<unsigned int>(sister) < h_partons->size() ? CandidatePtr(h_partons, sister)
                                                                                   : CandidatePtr()),
      matchedJet_(matchedJet),
      sisterJet_(sisterJet) {}

FlavorHistory::FlavorHistory(FLAVOR_T flavorSource,
                             Handle<CandidateCollection> h_partons,
                             int parton,
                             int progenitor,
                             int sister,
                             reco::ShallowClonePtrCandidate const& matchedJet,
                             reco::ShallowClonePtrCandidate const& sisterJet)
    : flavorSource_(flavorSource),
      parton_(parton >= 0 && static_cast<unsigned int>(parton) < h_partons->size() ? CandidatePtr(h_partons, parton)
                                                                                   : CandidatePtr()),
      progenitor_(progenitor >= 0 && static_cast<unsigned int>(progenitor) < h_partons->size()
                      ? CandidatePtr(h_partons, progenitor)
                      : CandidatePtr()),
      sister_(sister >= 0 && static_cast<unsigned int>(sister) < h_partons->size() ? CandidatePtr(h_partons, sister)
                                                                                   : CandidatePtr()),
      matchedJet_(matchedJet),
      sisterJet_(sisterJet) {}

ostream& operator<<(ostream& out, Candidate const& cand) {
  char buff[1000];
  sprintf(buff,
          "%5d, status = %5d, nmo = %5d, nda = %5d, pt = %6.2f, eta = %6.2f, phi = %6.2f, m = %6.2f",
          cand.pdgId(),
          cand.status(),
          static_cast<int>(cand.numberOfMothers()),
          static_cast<int>(cand.numberOfDaughters()),
          cand.pt(),
          cand.eta(),
          cand.phi(),
          cand.mass());
  out << buff;
  return out;
}

ostream& operator<<(ostream& out, FlavorHistory const& cand) {
  out << "Source     = " << cand.flavorSource() << endl;
  if (cand.hasParton())
    out << "Parton     = " << cand.parton().key() << " : " << *(cand.parton()) << endl;
  if (cand.hasProgenitor())
    out << "Progenitor = " << cand.progenitor().key() << " : " << *(cand.progenitor()) << endl;
  if (cand.hasSister())
    out << "Sister     = " << cand.sister().key() << " : " << *(cand.sister()) << endl;
  if (cand.hasMatchedJet())
    out << "jet        = " << cand.matchedJet() << endl;
  if (cand.hasSisterJet())
    out << "sister jet = " << cand.sisterJet() << endl;
  if (cand.hasParton()) {
    out << "Ancestry: " << endl;
    Candidate const* ipar = cand.parton()->mother();
    while (ipar->numberOfMothers() > 0) {
      out << *ipar << endl;
      ipar = ipar->mother();
    }
  }
  return out;
}