TrackSelectorByRegion.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/global/EDProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"

#include "DataFormats/TrackReco/interface/TrackFwd.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/Common/interface/OwnVector.h"
#include "RecoTracker/TkTrackingRegions/interface/TrackingRegion.h"

#include <vector>
#include <memory>

class TrackSelectorByRegion final : public edm::global::EDProducer<> {
public:
  explicit TrackSelectorByRegion(const edm::ParameterSet& conf)
      : produceCollection_(conf.getParameter<bool>("produceTrackCollection")),
        produceMask_(conf.getParameter<bool>("produceMask")),
        tracksToken_(consumes<reco::TrackCollection>(conf.getParameter<edm::InputTag>("tracks"))),
        inputTrkRegionToken_(consumes<edm::OwnVector<TrackingRegion>>(conf.getParameter<edm::InputTag>("regions"))),
        outputTracksToken_(produceCollection_ ? produces<reco::TrackCollection>()
                                              : edm::EDPutTokenT<reco::TrackCollection>{}),
        outputMaskToken_(produceMask_ ? produces<std::vector<bool>>() : edm::EDPutTokenT<std::vector<bool>>{}) {}

  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag>("tracks", edm::InputTag("hltPixelTracks"));
    desc.add<edm::InputTag>("regions", edm::InputTag(""));
    desc.add<bool>("produceTrackCollection", true);
    desc.add<bool>("produceMask", true);
    descriptions.add("trackSelectorByRegion", desc);
  }

private:
  using MaskCollection = std::vector<bool>;

  void produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup&) const override {
    if (not produceCollection_ and not produceMask_)
      return;

    auto const& regions = iEvent.get(inputTrkRegionToken_);
    auto const& tracks = iEvent.get(tracksToken_);
    MaskCollection mask(tracks.size(), false);  // output mask

    for (auto const& region : regions) {
      region.checkTracks(tracks, mask);
    }

    if (produceCollection_) {
      reco::TrackCollection output_tracks;  // output collection with a (shallow) copy of the selected tracks
      size_t size = 0;
      // count the number of selected tracks
      for (size_t i = 0; i < mask.size(); i++) {
        size += mask[i];
      }
      output_tracks.reserve(size);
      for (size_t i = 0; i < mask.size(); i++) {
        if (mask[i])
          output_tracks.push_back(tracks[i]);
      }
      iEvent.emplace(outputTracksToken_, std::move(output_tracks));
    }
    if (produceMask_) {
      iEvent.emplace(outputMaskToken_, std::move(mask));
    }
  }

  const bool produceCollection_;
  const bool produceMask_;
  const edm::EDGetTokenT<reco::TrackCollection> tracksToken_;
  const edm::EDGetTokenT<edm::OwnVector<TrackingRegion>> inputTrkRegionToken_;
  const edm::EDPutTokenT<reco::TrackCollection> outputTracksToken_;
  const edm::EDPutTokenT<std::vector<bool>> outputMaskToken_;
};

#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(TrackSelectorByRegion);