RunLumiSelector.h

Go to the documentation of this file.

#ifndef Analysis_AnalysisFilters_interface_RunLumiSelector_h
#define Analysis_AnalysisFilters_interface_RunLumiSelector_h

#include "FWCore/Common/interface/EventBase.h"
#include "DataFormats/Common/interface/Handle.h"

#include "PhysicsTools/SelectorUtils/interface/EventSelector.h"
#include "DataFormats/VertexReco/interface/Vertex.h"

#include <vector>
#include <string>

// make a selector for this selection
class RunLumiSelector : public EventSelector {
public:
  RunLumiSelector() {}

  RunLumiSelector( edm::ParameterSet const & params ) {

    push_back("RunLumi");
    
    if ( params.exists("lumisToProcess") ) {
      lumis_ = params.getUntrackedParameter<std::vector<edm::LuminosityBlockRange> > ("lumisToProcess");
      set("RunLumi" );
    }
    else {
      lumis_.clear();
      set("RunLumi", false);
    }

    retInternal_ = getBitTemplate();
  }
  
  bool operator() ( edm::EventBase const & ev,  pat::strbitset & ret ) {

    if ( !ignoreCut("RunLumi") ) {
      bool goodLumi = false;
      for ( std::vector<edm::LuminosityBlockRange>::const_iterator lumisBegin = lumis_.begin(),
          lumisEnd = lumis_.end(), ilumi = lumisBegin;
        ilumi != lumisEnd; ++ilumi ) {
    if ( ev.id().run() >= ilumi->startRun() && ev.id().run() <= ilumi->endRun()  &&
         ev.id().luminosityBlock() >= ilumi->startLumi() && ev.id().luminosityBlock() <= ilumi->endLumi() )  {
      goodLumi = true;
      break;
    }
      }
      if ( goodLumi ) passCut(ret, "RunLumi" );
    } else {
      passCut(ret, "RunLumi");
    }

    setIgnored(ret);
    return (bool)ret;
  }

  using EventSelector::operator();

private:

  std::vector<edm::LuminosityBlockRange> lumis_;

};

#endif