GenericParticle.cc

Go to the documentation of this file.

//
//

#include "DataFormats/PatCandidates/interface/GenericParticle.h"

using pat::GenericParticle;

GenericParticle::GenericParticle() : PATObject<reco::RecoCandidate>() {}

GenericParticle::GenericParticle(const Candidate &cand) : PATObject<reco::RecoCandidate>() { fillInFrom(cand); }

GenericParticle::GenericParticle(const edm::RefToBase<Candidate> &cand) : PATObject<reco::RecoCandidate>() {
  fillInFrom(*cand);
  refToOrig_ = edm::Ptr<reco::Candidate>(cand.id(), cand.get(), cand.key());  // correct RefToBase=>Ptr conversion
}

GenericParticle::GenericParticle(const edm::Ptr<Candidate> &cand) : PATObject<reco::RecoCandidate>() {
  fillInFrom(*cand);
  refToOrig_ = cand;
}

GenericParticle::~GenericParticle() {}

// ====== SETTERS =====
void GenericParticle::setTrack(const reco::TrackRef &ref, bool embed) {
  trackRef_ = ref;
  if (embed)
    embedTrack();
  else
    track_.clear();
}

void GenericParticle::setTracks(const reco::TrackRefVector &refs, bool embed) {
  trackRefs_ = refs;
  if (embed)
    embedTracks();
  else
    tracks_.clear();
}

void GenericParticle::setStandAloneMuon(const reco::TrackRef &ref, bool embed) {
  standaloneTrackRef_ = ref;
  if (embed)
    embedStandalone();
  else
    standaloneTrack_.clear();
}

void GenericParticle::setCombinedMuon(const reco::TrackRef &ref, bool embed) {
  combinedTrackRef_ = ref;
  if (embed)
    embedCombined();
  else
    combinedTrack_.clear();
}

void GenericParticle::setGsfTrack(const reco::GsfTrackRef &ref, bool embed) {
  gsfTrackRef_ = ref;
  if (embed)
    embedGsfTrack();
  else
    gsfTrack_.clear();
}

void GenericParticle::setSuperCluster(const reco::SuperClusterRef &ref, bool embed) {
  superClusterRef_ = ref;
  if (embed)
    embedSuperCluster();
  else
    superCluster_.clear();
}

void GenericParticle::setCaloTower(const CaloTowerRef &ref, bool embed) {
  caloTowerRef_ = ref;
  if (embed) {
    embedCaloTower();
  } else if (!caloTower_.empty()) {
    CaloTowerCollection().swap(caloTower_);
  }
}

// ========== EMBEDDER METHODS
void GenericParticle::embedTrack() {
  track_.clear();
  if (trackRef_.isNonnull())
    track_.push_back(*trackRef_);  // import
  trackRef_ = reco::TrackRef();    // clear, to save space (zeroes compress better)
}
void GenericParticle::embedTracks() {
  tracks_.clear();
  tracks_.reserve(trackRefs_.size());
  for (reco::TrackRefVector::const_iterator it = trackRefs_.begin(); it != trackRefs_.end(); ++it) {
    if (it->isNonnull())
      tracks_.push_back(**it);  // embed track
  }
  trackRefs_ = reco::TrackRefVector();  // clear, to save space
}
void GenericParticle::embedStandalone() {
  standaloneTrack_.clear();
  if (standaloneTrackRef_.isNonnull())
    standaloneTrack_.push_back(*standaloneTrackRef_);  // import
  standaloneTrackRef_ = reco::TrackRef();              // clear, to save space (zeroes compress better)
}
void GenericParticle::embedCombined() {
  combinedTrack_.clear();
  if (combinedTrackRef_.isNonnull())
    combinedTrack_.push_back(*combinedTrackRef_);  // import
  combinedTrackRef_ = reco::TrackRef();            // clear, to save space (zeroes compress better)
}
void GenericParticle::embedGsfTrack() {
  gsfTrack_.clear();
  if (gsfTrackRef_.isNonnull())
    gsfTrack_.push_back(*gsfTrackRef_);  // import
  gsfTrackRef_ = reco::GsfTrackRef();    // clear, to save space (zeroes compress better)
}

void GenericParticle::embedSuperCluster() {
  superCluster_.clear();
  if (superClusterRef_.isNonnull())
    superCluster_.push_back(*superClusterRef_);  // import
  superClusterRef_ = reco::SuperClusterRef();    // clear, to save space (zeroes compress better)
}
void GenericParticle::embedCaloTower() {
  if (!caloTower_.empty())
    CaloTowerCollection().swap(caloTower_);
  if (caloTowerRef_.isNonnull())
    caloTower_.push_back(*caloTowerRef_);  // import
  caloTowerRef_ = CaloTowerRef();          // clear, to save space (zeroes compress better)
}

void GenericParticle::fillInFrom(const reco::Candidate &cand) {
  // first, kinematics & status
  setCharge(cand.charge());
  setP4(cand.polarP4());
  setVertex(cand.vertex());
  setPdgId(cand.pdgId());
  setStatus(cand.status());
  // then RECO part, if available
  const reco::RecoCandidate *rc = dynamic_cast<const reco::RecoCandidate *>(&cand);
  if (rc != nullptr) {
    setTrack(rc->track());
    setGsfTrack(rc->gsfTrack());
    setStandAloneMuon(rc->standAloneMuon());
    setCombinedMuon(rc->combinedMuon());
    setSuperCluster(rc->superCluster());
    setCaloTower(rc->caloTower());
    size_t ntracks = rc->numberOfTracks();
    if (ntracks > 0) {
      reco::TrackRefVector tracks;
      for (size_t i = 0; i < ntracks; ++i) {
        tracks.push_back(rc->track(i));
      }
      setTracks(tracks);
    }
  }
}

bool GenericParticle::overlap(const reco::Candidate &cand) const {
  const reco::RecoCandidate *rc = dynamic_cast<const reco::RecoCandidate *>(&cand);
  if (rc != nullptr) {
    if (rc->track().isNonnull() && (track() == rc->track()))
      return true;
    if (rc->gsfTrack().isNonnull() && (gsfTrack() == rc->gsfTrack()))
      return true;
    if (rc->standAloneMuon().isNonnull() && (standAloneMuon() == rc->standAloneMuon()))
      return true;
    if (rc->combinedMuon().isNonnull() && (combinedMuon() == rc->combinedMuon()))
      return true;
    if (rc->superCluster().isNonnull() && (superCluster() == rc->superCluster()))
      return true;
    if (rc->caloTower().isNonnull() && (caloTower() == rc->caloTower()))
      return true;
  }
  const GenericParticle *rc2 = dynamic_cast<const GenericParticle *>(&cand);
  if (rc2 != nullptr) {
    if (rc2->track().isNonnull() && (track() == rc2->track()))
      return true;
    if (rc2->gsfTrack().isNonnull() && (gsfTrack() == rc2->gsfTrack()))
      return true;
    if (rc2->standAloneMuon().isNonnull() && (standAloneMuon() == rc2->standAloneMuon()))
      return true;
    if (rc2->combinedMuon().isNonnull() && (combinedMuon() == rc2->combinedMuon()))
      return true;
    if (rc2->superCluster().isNonnull() && (superCluster() == rc2->superCluster()))
      return true;
    if (rc2->caloTower().isNonnull() && (caloTower() == rc2->caloTower()))
      return true;
  }
  return false;
}