L3MuonCleaner.cc

Go to the documentation of this file.

// -*- C++ -*-
// Framework
#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"

class L3MuonCleaner : public edm::global::EDProducer<> {
 public:
  L3MuonCleaner(const edm::ParameterSet&);
  virtual ~L3MuonCleaner(){}
  virtual void produce(edm::StreamID, edm::Event&, const edm::EventSetup&) const override;
 private:
  edm::InputTag m_input; 
  int m_minTrkHits;
  int m_minMuonHits;
  double m_maxNormalizedChi2;
  edm::EDGetTokenT<reco::TrackCollection> inputToken_;



};

L3MuonCleaner::L3MuonCleaner(const edm::ParameterSet& parameterSet){
  m_input = parameterSet.getParameter<edm::InputTag>("input");
  m_minTrkHits = parameterSet.getParameter<int>("minTrkHits");
  m_minMuonHits = parameterSet.getParameter<int>("minMuonHits");
  m_maxNormalizedChi2 = parameterSet.getParameter<double>("maxNormalizedChi2");
  inputToken_ = consumes<reco::TrackCollection>(m_input);

  produces<reco::TrackCollection>();
}

void L3MuonCleaner::produce(edm::StreamID, edm::Event& event, const edm::EventSetup&) const{
  edm::Handle<reco::TrackCollection> tracks; 
  event.getByToken(inputToken_,tracks);
  auto outTracks = std::make_unique<reco::TrackCollection>();
  for ( reco::TrackCollection::const_iterator trk=tracks->begin(); trk!=tracks->end(); ++trk ){
    if (trk->normalizedChi2()>m_maxNormalizedChi2) continue;
    if (trk->hitPattern().numberOfValidTrackerHits()<m_minTrkHits) continue;
    if (trk->hitPattern().numberOfValidMuonHits()<m_minMuonHits) continue;
    outTracks->push_back(*trk);
  }
  event.put(std::move(outTracks));
}
DEFINE_FWK_MODULE(L3MuonCleaner);