FastTrackerRecHitMaskProducer.cc

Go to the documentation of this file.

// system includes
#include <memory>

// framework stuff
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/ESHandle.h"

// data formats
#include "DataFormats/Common/interface/ValueMap.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/TrackerRecHit2D/interface/FastTrackerRecHit.h"
#include "DataFormats/TrackerRecHit2D/interface/FastTrackerRecHitCollection.h"

// other
#include "TrackingTools/PatternTools/interface/TrackCollectionTokens.h"

class FastTrackerRecHitMaskProducer : public edm::stream::EDProducer<> {
public:
  explicit FastTrackerRecHitMaskProducer(const edm::ParameterSet&);

  ~FastTrackerRecHitMaskProducer() override {}

  void produce(edm::Event&, const edm::EventSetup&) override;

private:
  // an alias
  using QualityMaskCollection = std::vector<unsigned char>;

  // tokens
  const int minNumberOfLayersWithMeasBeforeFiltering_;
  const reco::TrackBase::TrackQuality trackQuality_;

  const TrackCollectionTokens trajectories_;
  edm::EDGetTokenT<QualityMaskCollection> srcQuals;

  edm::EDGetTokenT<FastTrackerRecHitCollection> recHits_;

  edm::EDGetTokenT<std::vector<bool> > oldHitMaskToken_;
};

FastTrackerRecHitMaskProducer::FastTrackerRecHitMaskProducer(const edm::ParameterSet& iConfig)
    : minNumberOfLayersWithMeasBeforeFiltering_(iConfig.getParameter<int>("minNumberOfLayersWithMeasBeforeFiltering")),
      trackQuality_(reco::TrackBase::qualityByName(iConfig.getParameter<std::string>("TrackQuality"))),
      trajectories_(iConfig.getParameter<edm::InputTag>("trajectories"), consumesCollector()),
      recHits_(consumes<FastTrackerRecHitCollection>(iConfig.getParameter<edm::InputTag>("recHits"))) {
  produces<std::vector<bool> >();

  auto const& classifier = iConfig.getParameter<edm::InputTag>("trackClassifier");
  if (!classifier.label().empty())
    srcQuals = consumes<QualityMaskCollection>(classifier);

  auto const& oldHitRemovalInfo = iConfig.getParameter<edm::InputTag>("oldHitRemovalInfo");
  if (!oldHitRemovalInfo.label().empty()) {
    oldHitMaskToken_ = consumes<std::vector<bool> >(oldHitRemovalInfo);
  }
}

void FastTrackerRecHitMaskProducer::produce(edm::Event& iEvent, const edm::EventSetup& es) {
  // the product

  std::unique_ptr<std::vector<bool> > collectedHits(new std::vector<bool>());

  // input

  edm::Handle<FastTrackerRecHitCollection> recHits;
  iEvent.getByToken(recHits_, recHits);

  if (!oldHitMaskToken_.isUninitialized()) {
    edm::Handle<std::vector<bool> > oldHitMasks;
    iEvent.getByToken(oldHitMaskToken_, oldHitMasks);
    collectedHits->insert(collectedHits->begin(), oldHitMasks->begin(), oldHitMasks->end());
  }
  collectedHits->resize(recHits->size(), false);

  auto const& tracks = trajectories_.tracks(iEvent);

  QualityMaskCollection const* pquals = nullptr;
  if (!srcQuals.isUninitialized()) {
    edm::Handle<QualityMaskCollection> hqual;
    iEvent.getByToken(srcQuals, hqual);
    pquals = hqual.product();
  }

  // required quality
  unsigned char qualMask = ~0;
  if (trackQuality_ != reco::TrackBase::undefQuality)
    qualMask = 1 << trackQuality_;

  // loop over tracks and mask hits of selected tracks
  for (auto i = 0U; i < tracks.size(); ++i) {
    const reco::Track& track = tracks[i];
    bool goodTk = (pquals) ? (*pquals)[i] & qualMask : track.quality(trackQuality_);
    if (!goodTk)
      continue;
    if (track.hitPattern().trackerLayersWithMeasurement() < minNumberOfLayersWithMeasBeforeFiltering_)
      continue;
    for (auto hitIt = track.recHitsBegin(); hitIt != track.recHitsEnd(); ++hitIt) {
      if (!(*hitIt)->isValid())
        continue;
      const FastTrackerRecHit& hit = static_cast<const FastTrackerRecHit&>(*(*hitIt));
      // note: for matched hits nIds() returns 2, otherwise 1
      for (unsigned id_index = 0; id_index < hit.nIds(); id_index++) {
        (*collectedHits)[unsigned(hit.id(id_index))] = true;
      }
    }
  }

  iEvent.put(std::move(collectedHits));
}

DEFINE_FWK_MODULE(FastTrackerRecHitMaskProducer);