MaskedMeasurementTrackerEventProducer.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/stream/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "RecoTracker/MeasurementDet/interface/MeasurementTrackerEvent.h"
#include "DataFormats/DetId/interface/DetIdCollection.h"

class dso_hidden MaskedMeasurementTrackerEventProducer final : public edm::stream::EDProducer<> {
public:
      explicit MaskedMeasurementTrackerEventProducer(const edm::ParameterSet &iConfig) ;
      ~MaskedMeasurementTrackerEventProducer() override {}
private:
      void produce(edm::Event&, const edm::EventSetup&) override;

      typedef edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > StripMask;
      typedef edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > PixelMask;
      typedef edm::ContainerMask<edmNew::DetSetVector<Phase2TrackerCluster1D> > Phase2OTMask;

      edm::EDGetTokenT<MeasurementTrackerEvent> src_;

      bool skipClusters_;
      bool phase2skipClusters_;

      edm::EDGetTokenT<StripMask> maskStrips_;
      edm::EDGetTokenT<PixelMask> maskPixels_;
      edm::EDGetTokenT<Phase2OTMask> maskPhase2OTs_;
};


MaskedMeasurementTrackerEventProducer::MaskedMeasurementTrackerEventProducer(const edm::ParameterSet &iConfig) :
    src_(consumes<MeasurementTrackerEvent>(iConfig.getParameter<edm::InputTag>("src"))),
    skipClusters_(false), phase2skipClusters_(false)
{
    //FIXME:temporary solution in order to use this class for both phase0/1 and phase2
    if (iConfig.existsAs<edm::InputTag>("clustersToSkip")) {
      skipClusters_ = true;
      edm::InputTag clustersToSkip = iConfig.getParameter<edm::InputTag>("clustersToSkip");
      maskPixels_ = consumes<PixelMask>(clustersToSkip);
      maskStrips_ = consumes<StripMask>(clustersToSkip);
    }
    if (iConfig.existsAs<edm::InputTag>("phase2clustersToSkip")) {
      phase2skipClusters_ = true;
      edm::InputTag phase2clustersToSkip = iConfig.getParameter<edm::InputTag>("phase2clustersToSkip");
      maskPixels_ = consumes<PixelMask>(phase2clustersToSkip);
      maskPhase2OTs_ = consumes<Phase2OTMask>(phase2clustersToSkip);
    }
    produces<MeasurementTrackerEvent>();
}

void
MaskedMeasurementTrackerEventProducer::produce(edm::Event &iEvent, const edm::EventSetup& iSetup)
{
    edm::Handle<MeasurementTrackerEvent> mte;
    iEvent.getByToken(src_, mte);

    // prepare output
    std::unique_ptr<MeasurementTrackerEvent> out;

    if (skipClusters_) {

      edm::Handle<PixelMask> maskPixels;
      iEvent.getByToken(maskPixels_, maskPixels);
      edm::Handle<StripMask> maskStrips;
      iEvent.getByToken(maskStrips_, maskStrips);

      out = std::make_unique<MeasurementTrackerEvent>(*mte, *maskStrips, *maskPixels);

    } else if (phase2skipClusters_) {

      edm::Handle<PixelMask> maskPixels;
      iEvent.getByToken(maskPixels_, maskPixels);
      edm::Handle<Phase2OTMask> maskPhase2OTs;
      iEvent.getByToken(maskPhase2OTs_, maskPhase2OTs);

      out = std::make_unique<MeasurementTrackerEvent>(*mte, *maskPixels, *maskPhase2OTs);
    }

    // put into event
    iEvent.put(std::move(out));
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(MaskedMeasurementTrackerEventProducer);