/data/doxygen/doxygen-1.7.3/gen/CMSSW_4_2_8/src/DPGAnalysis/Skims/src/MuonPtFilter.cc

Go to the documentation of this file.

/* This Class Header */
#include "DPGAnalysis/Skims/interface/MuonPtFilter.h"

/* Collaborating Class Header */
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/ESHandle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include "DataFormats/TrackReco/interface/Track.h"
#include "DataFormats/TrackReco/interface/TrackFwd.h"




/* C++ Headers */
using namespace std;
using namespace edm;

/* ====================================================================== */

MuonPtFilter::MuonPtFilter(const edm::ParameterSet& pset) {

  // the name of the STA rec hits collection
  theSTAMuonLabel = pset.getParameter<std::string>("SALabel");

  theMinPt = pset.getParameter<double>("minPt"); // pt min (GeV)

  LogDebug("MuonPt") << " SALabel : " << theSTAMuonLabel 
    << " Min Pt : " << theMinPt;
}

MuonPtFilter::~MuonPtFilter() {
}

/* Operations */ 
bool MuonPtFilter::filter(edm::Event& event, const edm::EventSetup& eventSetup) {
  bool accept = false;

  // Get the RecTrack collection from the event
  Handle<reco::TrackCollection> staTracks;
  event.getByLabel(theSTAMuonLabel, staTracks);
  
  reco::TrackCollection::const_iterator staTrack;
  
  for (staTrack = staTracks->begin(); staTrack != staTracks->end(); ++staTrack){
    
    if(staTrack->pt()>theMinPt){
      accept=true;
      return accept;
    }

  }

  return accept;


}

// define this as a plug-in
DEFINE_FWK_MODULE(MuonPtFilter);