MuonVPlusJetsIDSelectionFunctor Class Reference

#include <MuonVPlusJetsIDSelectionFunctor.h>

Inheritance diagram for MuonVPlusJetsIDSelectionFunctor:

Public Types
enum	Version_t {   SUMMER08, FIRSTDATA, SPRING10, FALL10,   N_VERSIONS, KITQCD }
Public Types inherited from Selector< pat::Muon >
typedef std::pair< index_type, size_t >	cut_flow_item
typedef std::vector< cut_flow_item >	cut_flow_map
typedef pat::Muon	data_type
typedef std::map< index_type, double >	double_map
typedef pat::strbitset::index_type	index_type
typedef std::map< index_type, int >	int_map

Public Member Functions
bool	fall10Cuts (const pat::Muon &muon, edm::EventBase const &event, pat::strbitset &ret)
bool	firstDataCuts (const pat::Muon &muon, pat::strbitset &ret)
void	initialize (Version_t version, double chi2=10.0, double d0=999.0, double ed0=999.0, double sd0=3.0, int nhits=11, int nValidMuonHits=0, double ecalveto=4.0, double hcalveto=6.0, double reliso=0.05, double maxLepZ=1.0, int minPixelHits=1, int minNMatches=1)
bool	kitQCDCuts (const pat::Muon &muon, edm::EventBase const &event, pat::strbitset &ret)
	MuonVPlusJetsIDSelectionFunctor ()
	MuonVPlusJetsIDSelectionFunctor (edm::ParameterSet const &parameters, edm::ConsumesCollector &iC)
	MuonVPlusJetsIDSelectionFunctor (edm::ParameterSet const &parameters)
	MuonVPlusJetsIDSelectionFunctor (Version_t version, double chi2=10.0, double d0=0.2, double ed0=999.0, double sd0=999.0, int nhits=11, int nValidMuonHits=0, double ecalveto=4.0, double hcalveto=6.0, double reliso=0.05, double maxLepZ=1.0, int minPixelHits=1, int minNMatches=1)
bool	operator() (const pat::Muon &muon, edm::EventBase const &event, pat::strbitset &ret) override
	This provides an alternative signature that includes extra information. More...
bool	operator() (const pat::Muon &muon, pat::strbitset &ret) override
	This provides the interface for base classes to select objects. More...
bool	spring10Cuts (const pat::Muon &muon, edm::EventBase const &event, pat::strbitset &ret)
bool	summer08Cuts (const pat::Muon &muon, pat::strbitset &ret)
Public Member Functions inherited from Selector< pat::Muon >
void	clear (std::string const &s)
	Turn off a given selection cut. More...
void	clear (index_type const &i)
bool	considerCut (std::string const &s) const
	consider the cut at index "s" More...
bool	considerCut (index_type const &i) const
int	cut (index_type const &i, int val) const
	Access the int cut values at index "s". More...
double	cut (index_type const &i, double val) const
	Access the double cut values at index "s". More...
int	cut (std::string s, int val) const
	Access the int cut values at index "s". More...
double	cut (std::string s, double val) const
	Access the double cut values at index "s". More...
pat::strbitset	getBitTemplate () const
	Get an empty bitset with the proper names. More...
double	getPasses (std::string const &s) const
	Return the number of passing cases. More...
double	getPasses (index_type const &i) const
bool	ignoreCut (std::string const &s) const
	ignore the cut at index "s" More...
bool	ignoreCut (index_type const &i) const
virtual bool	operator() (pat::Muon const &t)
	This provides an alternative signature without the second ret. More...
virtual bool	operator() (pat::Muon const &t, edm::EventBase const &e)
	This provides an alternative signature that includes extra information. More...
bool	operator[] (std::string const &s) const
bool	operator[] (index_type const &i) const
void	passCut (pat::strbitset &ret, std::string const &s)
	Passing cuts. More...
void	passCut (pat::strbitset &ret, index_type const &i)
void	print (std::ostream &out) const
	Print the cut flow. More...
void	printActiveCuts (std::ostream &out) const
	Print the cuts being considered. More...
virtual void	push_back (std::string const &s)
	This is the registration of an individual cut string. More...
virtual void	push_back (std::string const &s, int cut)
	This is the registration of an individual cut string, with an int cut value. More...
virtual void	push_back (std::string const &s, double cut)
	This is the registration of an individual cut string, with a double cut value. More...
	Selector ()
	Constructor clears the bits. More...
void	set (std::string const &s, bool val=true)
	Set a given selection cut, on or off. More...
void	set (index_type const &i, bool val=true)
void	set (std::string const &s, int cut, bool val=true)
	Set a given selection cut, on or off, and reset int cut value. More...
void	set (index_type const &i, int cut, bool val=true)
void	set (std::string const &s, double cut, bool val=true)
	Set a given selection cut, on or off, and reset int cut value. More...
void	set (index_type const &i, double cut, bool val=true)
void	setIgnored (pat::strbitset &ret)
	set ignored bits More...
void	setIgnoredCuts (std::vector< std::string > const &bitsToIgnore)
	set the bits to ignore from a vector More...
virtual	~Selector ()

Public Attributes
bool	verbose_

Private Attributes
edm::InputTag	beamLineSrc_
edm::EDGetTokenT< reco::BeamSpot >	beamLineSrcToken_
index_type	indexChi2_
index_type	indexD0_
index_type	indexECalVeto_
index_type	indexED0_
index_type	indexHCalVeto_
index_type	indexLepZ_
index_type	indexNHits_
index_type	indexNValMuHits_
index_type	indexPixHits_
index_type	indexRelIso_
index_type	indexSD0_
index_type	indexStations_
edm::InputTag	pvSrc_
edm::EDGetTokenT< std::vector< reco::Vertex > >	pvSrcToken_
bool	recalcDBFromBSp_
Version_t	version_

Additional Inherited Members
Protected Attributes inherited from Selector< pat::Muon >
pat::strbitset	bits_
	the bitset indexed by strings More...
cut_flow_map	cutFlow_
	map of cut flows in "human" order More...
double_map	doubleCuts_
	the double-value cut map More...
int_map	intCuts_
	the int-value cut map More...
pat::strbitset	retInternal_
	internal ret if users don't care about return bits More...

Detailed Description

Definition at line 17 of file MuonVPlusJetsIDSelectionFunctor.h.

Member Enumeration Documentation

Version_t

enum MuonVPlusJetsIDSelectionFunctor::Version_t

Enumerator
SUMMER08
FIRSTDATA
SPRING10
FALL10
N_VERSIONS
KITQCD

Definition at line 21 of file MuonVPlusJetsIDSelectionFunctor.h.

{ SUMMER08, FIRSTDATA, SPRING10, FALL10, N_VERSIONS, KITQCD };

Constructor & Destructor Documentation

MuonVPlusJetsIDSelectionFunctor() [1/4]

MuonVPlusJetsIDSelectionFunctor::MuonVPlusJetsIDSelectionFunctor ( )

inline

Definition at line 23 of file MuonVPlusJetsIDSelectionFunctor.h.

{}

MuonVPlusJetsIDSelectionFunctor() [2/4]

MuonVPlusJetsIDSelectionFunctor::MuonVPlusJetsIDSelectionFunctor ( edm::ParameterSet const &  parameters, edm::ConsumesCollector &  iC  )

inline

Definition at line 26 of file MuonVPlusJetsIDSelectionFunctor.h.

References beamLineSrc_, beamLineSrcToken_, edm::ConsumesCollector::consumes(), pvSrc_, and pvSrcToken_.

      : MuonVPlusJetsIDSelectionFunctor(parameters) {
    beamLineSrcToken_ = iC.consumes<reco::BeamSpot>(beamLineSrc_);
    pvSrcToken_ = iC.consumes<std::vector<reco::Vertex> >(pvSrc_);
  }

MuonVPlusJetsIDSelectionFunctor() [3/4]

MuonVPlusJetsIDSelectionFunctor::MuonVPlusJetsIDSelectionFunctor ( edm::ParameterSet const &  parameters )

inline

Definition at line 33 of file MuonVPlusJetsIDSelectionFunctor.h.

References gather_cfg::cout, Exception, FALL10, FIRSTDATA, Selector< pat::Muon >::getBitTemplate(), initialize(), KITQCD, N_VERSIONS, pvSrc_, recalcDBFromBSp_, Selector< pat::Muon >::retInternal_, Selector< pat::Muon >::setIgnoredCuts(), SPRING10, AlCaHLTBitMon_QueryRunRegistry::string, SUMMER08, verbose_, and BeamSplash_cfg::version.

                                                                     {
    verbose_ = false;

    std::string versionStr = parameters.getParameter<std::string>("version");

    Version_t version = N_VERSIONS;

    if (versionStr == "SUMMER08") {
      version = SUMMER08;
    } else if (versionStr == "FIRSTDATA") {
      version = FIRSTDATA;
    } else if (versionStr == "SPRING10") {
      version = SPRING10;
    } else if (versionStr == "FALL10") {
      version = FALL10;
      if (verbose_)
        std::cout << "\nMUON SELECTION - you are using FALL10 Selection" << std::endl;
    } else if (versionStr == "KITQCD") {
      version = KITQCD;
      if (verbose_)
        std::cout << "\nMUON SELECTION - you are using KITQCD Selection" << std::endl;
    } else {
      throw cms::Exception("InvalidInput")
          << "Expect version to be one of SUMMER08, FIRSTDATA, SPRING10, FALL10" << std::endl;
    }

    initialize(version,
               parameters.getParameter<double>("Chi2"),
               parameters.getParameter<double>("D0"),
               parameters.getParameter<double>("ED0"),
               parameters.getParameter<double>("SD0"),
               parameters.getParameter<int>("NHits"),
               parameters.getParameter<int>("NValMuHits"),
               parameters.getParameter<double>("ECalVeto"),
               parameters.getParameter<double>("HCalVeto"),
               parameters.getParameter<double>("RelIso"),
               parameters.getParameter<double>("LepZ"),
               parameters.getParameter<int>("nPixelHits"),
               parameters.getParameter<int>("nMatchedStations"));
    if (parameters.exists("cutsToIgnore"))
      setIgnoredCuts(parameters.getParameter<std::vector<std::string> >("cutsToIgnore"));

    retInternal_ = getBitTemplate();

    recalcDBFromBSp_ = parameters.getParameter<bool>("RecalcFromBeamSpot");
    pvSrc_ = parameters.getParameter<edm::InputTag>("pvSrc");
  }

MuonVPlusJetsIDSelectionFunctor() [4/4]

MuonVPlusJetsIDSelectionFunctor::MuonVPlusJetsIDSelectionFunctor ( Version_t  version, double  chi2 = 10.0, double  d0 = 0.2, double  ed0 = 999.0, double  sd0 = 999.0, int  nhits = 11, int  nValidMuonHits = 0, double  ecalveto = 4.0, double  hcalveto = 6.0, double  reliso = 0.05, double  maxLepZ = 1.0, int  minPixelHits = 1, int  minNMatches = 1  )

inline

Definition at line 81 of file MuonVPlusJetsIDSelectionFunctor.h.

References hltPixelTracks_cff::chi2, d0, Selector< pat::Muon >::getBitTemplate(), initialize(), beam_dqm_sourceclient-live_cfg::minPixelHits, nhits, Selector< pat::Muon >::retInternal_, and BeamSplash_cfg::version.

      : recalcDBFromBSp_(false) {
    initialize(version,
               chi2,
               d0,
               ed0,
               sd0,
               nhits,
               nValidMuonHits,
               ecalveto,
               hcalveto,
               reliso,
               maxLepZ,
               minPixelHits,
               minNMatches);

    retInternal_ = getBitTemplate();
  }

Member Function Documentation

fall10Cuts()

bool MuonVPlusJetsIDSelectionFunctor::fall10Cuts ( const pat::Muon &  muon, edm::EventBase const &  event, pat::strbitset &  ret  )

inline

Definition at line 388 of file MuonVPlusJetsIDSelectionFunctor.h.

References beamLineSrc_, pwdgSkimBPark_cfi::beamSpot, Selector< pat::Muon >::cut(), Exception, Selector< pat::Muon >::ignoreCut(), indexChi2_, indexD0_, indexECalVeto_, indexED0_, indexHCalVeto_, indexLepZ_, indexNHits_, indexNValMuHits_, indexPixHits_, indexRelIso_, indexSD0_, indexStations_, susybsm::HSCParticleType::innerTrack, edm::HandleBase::isValid(), nhits, Selector< pat::Muon >::passCut(), DiDispStaMuonMonitor_cfi::pt, pvSrc_, recalcDBFromBSp_, objects.autophobj::relIso, runTheMatrix::ret, Selector< pat::Muon >::setIgnored(), and mathSSE::sqrt().

Referenced by operator()().

                                                                                     {
    ret.set(false);

    double norm_chi2 = muon.normChi2();
    double corr_d0 = muon.dB();
    double corr_ed0 = muon.edB();
    double corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

    // Get the PV for the muon z requirement
    edm::Handle<std::vector<reco::Vertex> > pvtxHandle_;
    event.getByLabel(pvSrc_, pvtxHandle_);

    double zvtx = -999;
    if (!pvtxHandle_->empty()) {
      zvtx = pvtxHandle_->at(0).z();
    } else {
      throw cms::Exception("InvalidInput")
          << " There needs to be at least one primary vertex in the event." << std::endl;
    }

    //If required, recalculate the impact parameter using the beam spot
    if (recalcDBFromBSp_) {
      //Get the beam spot
      reco::TrackBase::Point beamPoint(0, 0, 0);
      reco::BeamSpot beamSpot;
      edm::Handle<reco::BeamSpot> beamSpotHandle;
      event.getByLabel(beamLineSrc_, beamSpotHandle);

      if (beamSpotHandle.isValid()) {
        beamSpot = *beamSpotHandle;
      } else {
        edm::LogError("DataNotAvailable")
            << "No beam spot available from EventSetup, not adding high level selection \n";
      }
      beamPoint = reco::TrackBase::Point(beamSpot.x0(), beamSpot.y0(), beamSpot.z0());

      //Use the beamspot to correct the impact parameter and uncertainty
      reco::TrackRef innerTrack = muon.innerTrack();
      if (innerTrack.isNonnull() && innerTrack.isAvailable()) {
        corr_d0 = -1.0 * innerTrack->dxy(beamPoint);
        corr_ed0 =
            sqrt(innerTrack->d0Error() * innerTrack->d0Error() + 0.5 * beamSpot.BeamWidthX() * beamSpot.BeamWidthX() +
                 0.5 * beamSpot.BeamWidthY() * beamSpot.BeamWidthY());
        corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

      } else {
        corr_d0 = 999.;
        corr_ed0 = 999.;
      }
    }

    int nhits = static_cast<int>(muon.numberOfValidHits());
    int nValidMuonHits = static_cast<int>(muon.globalTrack()->hitPattern().numberOfValidMuonHits());

    double ecalVeto = muon.isolationR03().emVetoEt;
    double hcalVeto = muon.isolationR03().hadVetoEt;

    double hcalIso = muon.hcalIso();
    double ecalIso = muon.ecalIso();
    double trkIso = muon.trackIso();
    double pt = muon.pt();

    double relIso = (ecalIso + hcalIso + trkIso) / pt;

    double z_mu = muon.vertex().z();

    int nPixelHits = muon.innerTrack()->hitPattern().pixelLayersWithMeasurement();

    int nMatchedStations = muon.numberOfMatches();

    if (norm_chi2 < cut(indexChi2_, double()) || ignoreCut(indexChi2_))
      passCut(ret, indexChi2_);
    if (fabs(corr_d0) < cut(indexD0_, double()) || ignoreCut(indexD0_))
      passCut(ret, indexD0_);
    if (fabs(corr_ed0) < cut(indexED0_, double()) || ignoreCut(indexED0_))
      passCut(ret, indexED0_);
    if (fabs(corr_sd0) < cut(indexSD0_, double()) || ignoreCut(indexSD0_))
      passCut(ret, indexSD0_);
    if (nhits >= cut(indexNHits_, int()) || ignoreCut(indexNHits_))
      passCut(ret, indexNHits_);
    if (nValidMuonHits > cut(indexNValMuHits_, int()) || ignoreCut(indexNValMuHits_))
      passCut(ret, indexNValMuHits_);
    if (hcalVeto < cut(indexHCalVeto_, double()) || ignoreCut(indexHCalVeto_))
      passCut(ret, indexHCalVeto_);
    if (ecalVeto < cut(indexECalVeto_, double()) || ignoreCut(indexECalVeto_))
      passCut(ret, indexECalVeto_);
    if (relIso < cut(indexRelIso_, double()) || ignoreCut(indexRelIso_))
      passCut(ret, indexRelIso_);
    if (fabs(z_mu - zvtx) < cut(indexLepZ_, double()) || ignoreCut(indexLepZ_))
      passCut(ret, indexLepZ_);
    if (nPixelHits > cut(indexPixHits_, int()) || ignoreCut(indexPixHits_))
      passCut(ret, indexPixHits_);
    if (nMatchedStations > cut(indexStations_, int()) || ignoreCut(indexStations_))
      passCut(ret, indexStations_);

    setIgnored(ret);

    return (bool)ret;
  }

firstDataCuts()

bool MuonVPlusJetsIDSelectionFunctor::firstDataCuts ( const pat::Muon &  muon, pat::strbitset &  ret  )

inline

Definition at line 269 of file MuonVPlusJetsIDSelectionFunctor.h.

References Selector< pat::Muon >::cut(), Selector< pat::Muon >::ignoreCut(), indexChi2_, indexD0_, indexECalVeto_, indexED0_, indexHCalVeto_, indexNHits_, indexRelIso_, indexSD0_, nhits, Selector< pat::Muon >::passCut(), DiDispStaMuonMonitor_cfi::pt, objects.autophobj::relIso, runTheMatrix::ret, and Selector< pat::Muon >::setIgnored().

Referenced by operator()().

                                                             {
    ret.set(false);

    double norm_chi2 = muon.normChi2();
    double corr_d0 = muon.dB();
    double corr_ed0 = muon.edB();
    double corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;
    int nhits = static_cast<int>(muon.numberOfValidHits());

    double ecalVeto = muon.isolationR03().emVetoEt;
    double hcalVeto = muon.isolationR03().hadVetoEt;

    double hcalIso = muon.hcalIso();
    double ecalIso = muon.ecalIso();
    double trkIso = muon.trackIso();
    double pt = muon.pt();

    double relIso = (ecalIso + hcalIso + trkIso) / pt;

    if (norm_chi2 < cut(indexChi2_, double()) || ignoreCut(indexChi2_))
      passCut(ret, indexChi2_);
    if (fabs(corr_d0) < cut(indexD0_, double()) || ignoreCut(indexD0_))
      passCut(ret, indexD0_);
    if (fabs(corr_ed0) < cut(indexED0_, double()) || ignoreCut(indexED0_))
      passCut(ret, indexED0_);
    if (fabs(corr_sd0) < cut(indexSD0_, double()) || ignoreCut(indexSD0_))
      passCut(ret, indexSD0_);
    if (nhits >= cut(indexNHits_, int()) || ignoreCut(indexNHits_))
      passCut(ret, indexNHits_);
    if (hcalVeto < cut(indexHCalVeto_, double()) || ignoreCut(indexHCalVeto_))
      passCut(ret, indexHCalVeto_);
    if (ecalVeto < cut(indexECalVeto_, double()) || ignoreCut(indexECalVeto_))
      passCut(ret, indexECalVeto_);
    if (relIso < cut(indexRelIso_, double()) || ignoreCut(indexRelIso_))
      passCut(ret, indexRelIso_);

    setIgnored(ret);

    return (bool)ret;
  }

initialize()

void MuonVPlusJetsIDSelectionFunctor::initialize ( Version_t  version, double  chi2 = 10.0, double  d0 = 999.0, double  ed0 = 999.0, double  sd0 = 3.0, int  nhits = 11, int  nValidMuonHits = 0, double  ecalveto = 4.0, double  hcalveto = 6.0, double  reliso = 0.05, double  maxLepZ = 1.0, int  minPixelHits = 1, int  minNMatches = 1  )

inline

Definition at line 112 of file MuonVPlusJetsIDSelectionFunctor.h.

References Selector< pat::Muon >::bits_, hltPixelTracks_cff::chi2, d0, FALL10, FIRSTDATA, indexChi2_, indexD0_, indexECalVeto_, indexED0_, indexHCalVeto_, indexLepZ_, indexNHits_, indexNValMuHits_, indexPixHits_, indexRelIso_, indexSD0_, indexStations_, beam_dqm_sourceclient-live_cfg::minPixelHits, nhits, Selector< pat::Muon >::push_back(), SPRING10, SUMMER08, BeamSplash_cfg::version, and version_.

Referenced by MuonVPlusJetsIDSelectionFunctor().

                                       {
    version_ = version;

    push_back("Chi2", chi2);
    push_back("D0", d0);
    push_back("ED0", ed0);
    push_back("SD0", sd0);
    push_back("NHits", nhits);
    push_back("NValMuHits", nValidMuonHits);
    push_back("ECalVeto", ecalveto);
    push_back("HCalVeto", hcalveto);
    push_back("RelIso", reliso);
    push_back("LepZ", maxLepZ);
    push_back("nPixelHits", minPixelHits);
    push_back("nMatchedStations", minNMatches);

    set("Chi2");
    set("D0");
    set("ED0");
    set("SD0");
    set("NHits");
    set("NValMuHits");
    set("ECalVeto");
    set("HCalVeto");
    set("RelIso");
    set("LepZ");
    set("nPixelHits");
    set("nMatchedStations");

    indexChi2_ = index_type(&bits_, "Chi2");
    indexD0_ = index_type(&bits_, "D0");
    indexED0_ = index_type(&bits_, "ED0");
    indexSD0_ = index_type(&bits_, "SD0");
    indexNHits_ = index_type(&bits_, "NHits");
    indexNValMuHits_ = index_type(&bits_, "NValMuHits");
    indexECalVeto_ = index_type(&bits_, "ECalVeto");
    indexHCalVeto_ = index_type(&bits_, "HCalVeto");
    indexRelIso_ = index_type(&bits_, "RelIso");
    indexLepZ_ = index_type(&bits_, "LepZ");
    indexPixHits_ = index_type(&bits_, "nPixelHits");
    indexStations_ = index_type(&bits_, "nMatchedStations");

    if (version == FALL10) {
      set("ED0", false);
      set("SD0", false);
      set("ECalVeto", false);
      set("HCalVeto", false);
    } else if (version == SPRING10) {
      set("ED0", false);
      set("SD0", false);
      set("ECalVeto", false);
      set("HCalVeto", false);
      set("LepZ", false);
      set("nPixelHits", false);
      set("nMatchedStations", false);
    } else if (version_ == FIRSTDATA) {
      set("D0", false);
      set("ED0", false);
      set("NValMuHits", false);
      set("LepZ", false);
      set("nPixelHits", false);
      set("nMatchedStations", false);
    } else if (version == SUMMER08) {
      set("SD0", false);
      set("NValMuHits", false);
      set("LepZ", false);
      set("nPixelHits", false);
      set("nMatchedStations", false);
    }
  }

kitQCDCuts()

bool MuonVPlusJetsIDSelectionFunctor::kitQCDCuts ( const pat::Muon &  muon, edm::EventBase const &  event, pat::strbitset &  ret  )

inline

Definition at line 492 of file MuonVPlusJetsIDSelectionFunctor.h.

References beamLineSrc_, pwdgSkimBPark_cfi::beamSpot, Selector< pat::Muon >::cut(), Exception, Selector< pat::Muon >::ignoreCut(), indexChi2_, indexD0_, indexECalVeto_, indexED0_, indexHCalVeto_, indexLepZ_, indexNHits_, indexNValMuHits_, indexPixHits_, indexRelIso_, indexSD0_, indexStations_, susybsm::HSCParticleType::innerTrack, edm::HandleBase::isValid(), nhits, Selector< pat::Muon >::passCut(), DiDispStaMuonMonitor_cfi::pt, pvSrc_, recalcDBFromBSp_, objects.autophobj::relIso, runTheMatrix::ret, Selector< pat::Muon >::setIgnored(), and mathSSE::sqrt().

Referenced by operator()().

                                                                                     {
    ret.set(false);

    double norm_chi2 = muon.normChi2();
    double corr_d0 = muon.dB();
    double corr_ed0 = muon.edB();
    double corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

    // Get the PV for the muon z requirement
    edm::Handle<std::vector<reco::Vertex> > pvtxHandle_;
    event.getByLabel(pvSrc_, pvtxHandle_);

    double zvtx = -999;
    if (!pvtxHandle_->empty()) {
      zvtx = pvtxHandle_->at(0).z();
    } else {
      throw cms::Exception("InvalidInput")
          << " There needs to be at least one primary vertex in the event." << std::endl;
    }

    //If required, recalculate the impact parameter using the beam spot
    if (recalcDBFromBSp_) {
      //Get the beam spot
      reco::TrackBase::Point beamPoint(0, 0, 0);
      reco::BeamSpot beamSpot;
      edm::Handle<reco::BeamSpot> beamSpotHandle;
      event.getByLabel(beamLineSrc_, beamSpotHandle);

      if (beamSpotHandle.isValid()) {
        beamSpot = *beamSpotHandle;
      } else {
        edm::LogError("DataNotAvailable")
            << "No beam spot available from EventSetup, not adding high level selection \n";
      }
      beamPoint = reco::TrackBase::Point(beamSpot.x0(), beamSpot.y0(), beamSpot.z0());

      //Use the beamspot to correct the impact parameter and uncertainty
      reco::TrackRef innerTrack = muon.innerTrack();
      if (innerTrack.isNonnull() && innerTrack.isAvailable()) {
        corr_d0 = -1.0 * innerTrack->dxy(beamPoint);
        corr_ed0 =
            sqrt(innerTrack->d0Error() * innerTrack->d0Error() + 0.5 * beamSpot.BeamWidthX() * beamSpot.BeamWidthX() +
                 0.5 * beamSpot.BeamWidthY() * beamSpot.BeamWidthY());
        corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

      } else {
        corr_d0 = 999.;
        corr_ed0 = 999.;
      }
    }

    int nhits = static_cast<int>(muon.numberOfValidHits());
    int nValidMuonHits = static_cast<int>(muon.globalTrack()->hitPattern().numberOfValidMuonHits());

    double ecalVeto = muon.isolationR03().emVetoEt;
    double hcalVeto = muon.isolationR03().hadVetoEt;

    double hcalIso = muon.hcalIso();
    double ecalIso = muon.ecalIso();
    double trkIso = muon.trackIso();
    double pt = muon.pt();

    double relIso = (ecalIso + hcalIso + trkIso) / pt;

    double z_mu = muon.vertex().z();

    int nPixelHits = muon.innerTrack()->hitPattern().pixelLayersWithMeasurement();

    int nMatchedStations = muon.numberOfMatches();

    if (norm_chi2 < cut(indexChi2_, double()) || ignoreCut(indexChi2_))
      passCut(ret, indexChi2_);
    if (fabs(corr_d0) < cut(indexD0_, double()) || ignoreCut(indexD0_))
      passCut(ret, indexD0_);
    if (fabs(corr_ed0) < cut(indexED0_, double()) || ignoreCut(indexED0_))
      passCut(ret, indexED0_);
    if (fabs(corr_sd0) < cut(indexSD0_, double()) || ignoreCut(indexSD0_))
      passCut(ret, indexSD0_);
    if (nhits >= cut(indexNHits_, int()) || ignoreCut(indexNHits_))
      passCut(ret, indexNHits_);
    if (nValidMuonHits > cut(indexNValMuHits_, int()) || ignoreCut(indexNValMuHits_))
      passCut(ret, indexNValMuHits_);
    if (hcalVeto < cut(indexHCalVeto_, double()) || ignoreCut(indexHCalVeto_))
      passCut(ret, indexHCalVeto_);
    if (ecalVeto < cut(indexECalVeto_, double()) || ignoreCut(indexECalVeto_))
      passCut(ret, indexECalVeto_);
    if (fabs(z_mu - zvtx) < cut(indexLepZ_, double()) || ignoreCut(indexLepZ_))
      passCut(ret, indexLepZ_);
    if (nPixelHits > cut(indexPixHits_, int()) || ignoreCut(indexPixHits_))
      passCut(ret, indexPixHits_);
    if (nMatchedStations > cut(indexStations_, int()) || ignoreCut(indexStations_))
      passCut(ret, indexStations_);

    //
    // JMS Dec 13 2010
    // HACK
    // Need double-sided relIso cut to implement data-driven QCD
    //
    //
    if (((relIso > 0.2) && (relIso < 0.75)) || ignoreCut(indexRelIso_))
      passCut(ret, indexRelIso_);

    setIgnored(ret);

    return (bool)ret;
  }

operator()() [1/2]

bool MuonVPlusJetsIDSelectionFunctor::operator() ( const pat::Muon &  t, edm::EventBase const &  e, pat::strbitset &  ret  )

inlineoverridevirtual

This provides an alternative signature that includes extra information.

Reimplemented from Selector< pat::Muon >.

Definition at line 196 of file MuonVPlusJetsIDSelectionFunctor.h.

References gather_cfg::cout, FALL10, fall10Cuts(), FIRSTDATA, firstDataCuts(), KITQCD, kitQCDCuts(), runTheMatrix::ret, SPRING10, spring10Cuts(), SUMMER08, summer08Cuts(), verbose_, and version_.

                                                                                              {
    if (version_ == FALL10)
      return fall10Cuts(muon, event, ret);
    else if (version_ == SPRING10)
      return spring10Cuts(muon, event, ret);
    else if (version_ == SUMMER08)
      return summer08Cuts(muon, ret);
    else if (version_ == FIRSTDATA)
      return firstDataCuts(muon, ret);
    else if (version_ == KITQCD) {
      if (verbose_)
        std::cout << "Calling KIT selection method" << std::endl;
      return kitQCDCuts(muon, event, ret);
    } else {
      return false;
    }
  }

operator()() [2/2]

bool MuonVPlusJetsIDSelectionFunctor::operator() ( const pat::Muon &  t, pat::strbitset &  ret  )

inlineoverridevirtual

This provides the interface for base classes to select objects.

Implements Selector< pat::Muon >.

Definition at line 215 of file MuonVPlusJetsIDSelectionFunctor.h.

References Exception, FALL10, FIRSTDATA, firstDataCuts(), runTheMatrix::ret, SPRING10, SUMMER08, summer08Cuts(), and version_.

                                                                   {
    if (version_ == SPRING10 || version_ == FALL10)
      throw cms::Exception("LogicError") << "MuonVPlusJetsSelectionFunctor SPRING10 and FALL10 versions needs the "
                                            "event! Call operator()(muon,event,ret)"
                                         << std::endl;

    else if (version_ == SUMMER08)
      return summer08Cuts(muon, ret);
    else if (version_ == FIRSTDATA)
      return firstDataCuts(muon, ret);
    else {
      return false;
    }
  }

spring10Cuts()

bool MuonVPlusJetsIDSelectionFunctor::spring10Cuts ( const pat::Muon &  muon, edm::EventBase const &  event, pat::strbitset &  ret  )

inline

Definition at line 311 of file MuonVPlusJetsIDSelectionFunctor.h.

References beamLineSrc_, pwdgSkimBPark_cfi::beamSpot, Selector< pat::Muon >::cut(), Selector< pat::Muon >::ignoreCut(), indexChi2_, indexD0_, indexECalVeto_, indexED0_, indexHCalVeto_, indexNHits_, indexNValMuHits_, indexRelIso_, indexSD0_, susybsm::HSCParticleType::innerTrack, edm::HandleBase::isValid(), nhits, Selector< pat::Muon >::passCut(), DiDispStaMuonMonitor_cfi::pt, recalcDBFromBSp_, objects.autophobj::relIso, runTheMatrix::ret, Selector< pat::Muon >::setIgnored(), and mathSSE::sqrt().

Referenced by operator()().

                                                                                       {
    ret.set(false);

    double norm_chi2 = muon.normChi2();
    double corr_d0 = muon.dB();
    double corr_ed0 = muon.edB();
    double corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

    //If required, recalculate the impact parameter using the beam spot
    if (recalcDBFromBSp_) {
      //Get the beam spot
      reco::TrackBase::Point beamPoint(0, 0, 0);
      reco::BeamSpot beamSpot;
      edm::Handle<reco::BeamSpot> beamSpotHandle;
      event.getByLabel(beamLineSrc_, beamSpotHandle);

      if (beamSpotHandle.isValid()) {
        beamSpot = *beamSpotHandle;
      } else {
        edm::LogError("DataNotAvailable")
            << "No beam spot available from EventSetup, not adding high level selection \n";
      }
      beamPoint = reco::TrackBase::Point(beamSpot.x0(), beamSpot.y0(), beamSpot.z0());

      //Use the beamspot to correct the impact parameter and uncertainty
      reco::TrackRef innerTrack = muon.innerTrack();
      if (innerTrack.isNonnull() && innerTrack.isAvailable()) {
        corr_d0 = -1.0 * innerTrack->dxy(beamPoint);
        corr_ed0 =
            sqrt(innerTrack->d0Error() * innerTrack->d0Error() + 0.5 * beamSpot.BeamWidthX() * beamSpot.BeamWidthX() +
                 0.5 * beamSpot.BeamWidthY() * beamSpot.BeamWidthY());
        corr_sd0 = (corr_ed0 > 0.000000001) ? corr_d0 / corr_ed0 : 999.0;

      } else {
        corr_d0 = 999.;
        corr_ed0 = 999.;
      }
    }

    int nhits = static_cast<int>(muon.numberOfValidHits());
    int nValidMuonHits = static_cast<int>(muon.globalTrack()->hitPattern().numberOfValidMuonHits());

    double ecalVeto = muon.isolationR03().emVetoEt;
    double hcalVeto = muon.isolationR03().hadVetoEt;

    double hcalIso = muon.hcalIso();
    double ecalIso = muon.ecalIso();
    double trkIso = muon.trackIso();
    double pt = muon.pt();

    double relIso = (ecalIso + hcalIso + trkIso) / pt;

    if (norm_chi2 < cut(indexChi2_, double()) || ignoreCut(indexChi2_))
      passCut(ret, indexChi2_);
    if (fabs(corr_d0) < cut(indexD0_, double()) || ignoreCut(indexD0_))
      passCut(ret, indexD0_);
    if (fabs(corr_ed0) < cut(indexED0_, double()) || ignoreCut(indexED0_))
      passCut(ret, indexED0_);
    if (fabs(corr_sd0) < cut(indexSD0_, double()) || ignoreCut(indexSD0_))
      passCut(ret, indexSD0_);
    if (nhits >= cut(indexNHits_, int()) || ignoreCut(indexNHits_))
      passCut(ret, indexNHits_);
    if (nValidMuonHits > cut(indexNValMuHits_, int()) || ignoreCut(indexNValMuHits_))
      passCut(ret, indexNValMuHits_);
    if (hcalVeto < cut(indexHCalVeto_, double()) || ignoreCut(indexHCalVeto_))
      passCut(ret, indexHCalVeto_);
    if (ecalVeto < cut(indexECalVeto_, double()) || ignoreCut(indexECalVeto_))
      passCut(ret, indexECalVeto_);
    if (relIso < cut(indexRelIso_, double()) || ignoreCut(indexRelIso_))
      passCut(ret, indexRelIso_);

    setIgnored(ret);

    return (bool)ret;
  }

summer08Cuts()

bool MuonVPlusJetsIDSelectionFunctor::summer08Cuts ( const pat::Muon &  muon, pat::strbitset &  ret  )

inline

Definition at line 233 of file MuonVPlusJetsIDSelectionFunctor.h.

References Selector< pat::Muon >::cut(), Selector< pat::Muon >::ignoreCut(), indexChi2_, indexD0_, indexECalVeto_, indexHCalVeto_, indexNHits_, indexRelIso_, nhits, Selector< pat::Muon >::passCut(), DiDispStaMuonMonitor_cfi::pt, objects.autophobj::relIso, runTheMatrix::ret, and Selector< pat::Muon >::setIgnored().

Referenced by operator()().

                                                            {
    ret.set(false);

    double norm_chi2 = muon.normChi2();
    double corr_d0 = muon.dB();
    int nhits = static_cast<int>(muon.numberOfValidHits());

    double ecalVeto = muon.isolationR03().emVetoEt;
    double hcalVeto = muon.isolationR03().hadVetoEt;

    double hcalIso = muon.hcalIso();
    double ecalIso = muon.ecalIso();
    double trkIso = muon.trackIso();
    double pt = muon.pt();

    double relIso = (ecalIso + hcalIso + trkIso) / pt;

    if (norm_chi2 < cut(indexChi2_, double()) || ignoreCut(indexChi2_))
      passCut(ret, indexChi2_);
    if (fabs(corr_d0) < cut(indexD0_, double()) || ignoreCut(indexD0_))
      passCut(ret, indexD0_);
    if (nhits >= cut(indexNHits_, int()) || ignoreCut(indexNHits_))
      passCut(ret, indexNHits_);
    if (hcalVeto < cut(indexHCalVeto_, double()) || ignoreCut(indexHCalVeto_))
      passCut(ret, indexHCalVeto_);
    if (ecalVeto < cut(indexECalVeto_, double()) || ignoreCut(indexECalVeto_))
      passCut(ret, indexECalVeto_);
    if (relIso < cut(indexRelIso_, double()) || ignoreCut(indexRelIso_))
      passCut(ret, indexRelIso_);

    setIgnored(ret);

    return (bool)ret;
  }

Member Data Documentation

beamLineSrc_

edm::InputTag MuonVPlusJetsIDSelectionFunctor::beamLineSrc_

private

Definition at line 604 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), kitQCDCuts(), MuonVPlusJetsIDSelectionFunctor(), and spring10Cuts().

beamLineSrcToken_

edm::EDGetTokenT<reco::BeamSpot> MuonVPlusJetsIDSelectionFunctor::beamLineSrcToken_

private

Definition at line 606 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by MuonVPlusJetsIDSelectionFunctor().

indexChi2_

index_type MuonVPlusJetsIDSelectionFunctor::indexChi2_

private

Definition at line 613 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexD0_

index_type MuonVPlusJetsIDSelectionFunctor::indexD0_

private

Definition at line 614 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexECalVeto_

index_type MuonVPlusJetsIDSelectionFunctor::indexECalVeto_

private

Definition at line 619 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexED0_

index_type MuonVPlusJetsIDSelectionFunctor::indexED0_

private

Definition at line 615 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), and spring10Cuts().

indexHCalVeto_

index_type MuonVPlusJetsIDSelectionFunctor::indexHCalVeto_

private

Definition at line 620 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexLepZ_

index_type MuonVPlusJetsIDSelectionFunctor::indexLepZ_

private

Definition at line 622 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), initialize(), and kitQCDCuts().

indexNHits_

index_type MuonVPlusJetsIDSelectionFunctor::indexNHits_

private

Definition at line 617 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexNValMuHits_

index_type MuonVPlusJetsIDSelectionFunctor::indexNValMuHits_

private

Definition at line 618 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), initialize(), kitQCDCuts(), and spring10Cuts().

indexPixHits_

index_type MuonVPlusJetsIDSelectionFunctor::indexPixHits_

private

Definition at line 623 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), initialize(), and kitQCDCuts().

indexRelIso_

index_type MuonVPlusJetsIDSelectionFunctor::indexRelIso_

private

Definition at line 621 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), spring10Cuts(), and summer08Cuts().

indexSD0_

index_type MuonVPlusJetsIDSelectionFunctor::indexSD0_

private

Definition at line 616 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), firstDataCuts(), initialize(), kitQCDCuts(), and spring10Cuts().

indexStations_

index_type MuonVPlusJetsIDSelectionFunctor::indexStations_

private

Definition at line 624 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), initialize(), and kitQCDCuts().

pvSrc_

edm::InputTag MuonVPlusJetsIDSelectionFunctor::pvSrc_

private

Definition at line 608 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), kitQCDCuts(), and MuonVPlusJetsIDSelectionFunctor().

pvSrcToken_

edm::EDGetTokenT<std::vector<reco::Vertex> > MuonVPlusJetsIDSelectionFunctor::pvSrcToken_

private

Definition at line 610 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by MuonVPlusJetsIDSelectionFunctor().

recalcDBFromBSp_

bool MuonVPlusJetsIDSelectionFunctor::recalcDBFromBSp_

private

Definition at line 603 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by fall10Cuts(), kitQCDCuts(), MuonVPlusJetsIDSelectionFunctor(), and spring10Cuts().

verbose_

bool MuonVPlusJetsIDSelectionFunctor::verbose_

Definition at line 19 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by MuonVPlusJetsIDSelectionFunctor(), and operator()().

version_

Version_t MuonVPlusJetsIDSelectionFunctor::version_

private

Definition at line 602 of file MuonVPlusJetsIDSelectionFunctor.h.

Referenced by initialize(), and operator()().