#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)

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

	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 (index_type const &i)

void	clear (std::string const &s)
	Turn off a given selection cut. More...

bool	considerCut (index_type const &i) const

bool	considerCut (std::string const &s) const
	consider the cut at index "s" More...

double	cut (index_type const &i, double val) const
	Access the double cut values at index "s". More...

int	cut (index_type const &i, 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...

int	cut (std::string s, int val) const
	Access the int cut values at index "s". More...

pat::strbitset	getBitTemplate () const
	Get an empty bitset with the proper names. More...

double	getPasses (index_type const &i) const

double	getPasses (std::string const &s) const
	Return the number of passing cases. More...

bool	ignoreCut (index_type const &i) const

bool	ignoreCut (std::string const &s) const
	ignore the cut at index "s" More...

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[] (index_type const &i) const

bool	operator[] (std::string const &s) const

void	passCut (pat::strbitset &ret, index_type const &i)

void	passCut (pat::strbitset &ret, std::string const &s)
	Passing cuts. More...

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, double cut)
	This is the registration of an individual cut string, with a double cut value. 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...

	Selector ()
	Constructor clears the bits. More...

void	set (index_type const &i, bool val=true)

void	set (index_type const &i, double cut, bool val=true)

void	set (index_type const &i, int cut, bool val=true)

void	set (std::string const &s, bool val=true)
	Set a given selection cut, on or off. More...

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 (std::string const &s, int cut, bool val=true)
	Set a given selection cut, on or off, and reset int cut value. More...

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_

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.

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

Constructor & Destructor Documentation

◆ MuonVPlusJetsIDSelectionFunctor() [1/4]

MuonVPlusJetsIDSelectionFunctor::MuonVPlusJetsIDSelectionFunctor ( )

inline

Definition at line 23 of file MuonVPlusJetsIDSelectionFunctor.h.

23 {}

◆ MuonVPlusJetsIDSelectionFunctor() [2/4]

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

inline

Definition at line 26 of file MuonVPlusJetsIDSelectionFunctor.h.

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

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

◆ MuonVPlusJetsIDSelectionFunctor() [3/4]

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

inline

Definition at line 33 of file MuonVPlusJetsIDSelectionFunctor.h.

                                                                      {
     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");
   }

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.

◆ 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.

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

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.

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.

                                                                                      {
     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;
   }

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()().

◆ firstDataCuts()

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

inline

Definition at line 269 of file MuonVPlusJetsIDSelectionFunctor.h.

                                                              {
     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;
   }

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()().

◆ 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.

                                        {
     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);
     }
   }

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(), Selector< pat::Muon >::set(), SPRING10, SUMMER08, BeamSplash_cfg::version, and version_.

Referenced by MuonVPlusJetsIDSelectionFunctor().

◆ kitQCDCuts()

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

inline

Definition at line 492 of file MuonVPlusJetsIDSelectionFunctor.h.

                                                                                      {
     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;
   }

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()().

◆ 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.

                                                                                               {
     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;
     }
   }

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

◆ 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.

                                                                    {
     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;
     }
   }

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

◆ spring10Cuts()

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

inline

Definition at line 311 of file MuonVPlusJetsIDSelectionFunctor.h.

                                                                                        {
     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;
   }

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()().

◆ summer08Cuts()

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

inline

Definition at line 233 of file MuonVPlusJetsIDSelectionFunctor.h.

                                                             {
     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;
   }

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()().

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()().

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_

Public Types

Public Member Functions

Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Member Enumeration Documentation

◆ Version_t

Constructor & Destructor Documentation

◆ MuonVPlusJetsIDSelectionFunctor() [1/4]

◆ MuonVPlusJetsIDSelectionFunctor() [2/4]

◆ MuonVPlusJetsIDSelectionFunctor() [3/4]

◆ MuonVPlusJetsIDSelectionFunctor() [4/4]

Member Function Documentation

◆ fall10Cuts()

◆ firstDataCuts()

◆ initialize()

◆ kitQCDCuts()

◆ operator()() [1/2]

◆ operator()() [2/2]

◆ spring10Cuts()

◆ summer08Cuts()

Member Data Documentation

◆ beamLineSrc_

◆ beamLineSrcToken_

◆ indexChi2_

◆ indexD0_

◆ indexECalVeto_

◆ indexED0_

◆ indexHCalVeto_

◆ indexLepZ_

◆ indexNHits_

◆ indexNValMuHits_

◆ indexPixHits_

◆ indexRelIso_

◆ indexSD0_

◆ indexStations_

◆ pvSrc_

◆ pvSrcToken_

◆ recalcDBFromBSp_

◆ verbose_

◆ version_