AcceptJet Class Reference

#include <AcceptJet.h>

Public Member Functions
	AcceptJet (const double &etaMin_, const double &etaMax_, const double &ptMin_, const double &ptMax_, const double &pMin_, const double &pMax_, const double &ratioMin_, const double &ratioMax_, const bool &doJetID_)
bool	operator() (const reco::Jet &jet, const int &jetFlavour, const edm::Handle< reco::SoftLeptonTagInfoCollection > &infos, const double jec) const
	Returns true if jet and associated parton satisfy kinematic cuts. More...
void	setDoJetID (bool b)
void	setEtaMax (double d)
void	setEtaMin (double d)
	Set cut parameters. More...
void	setPRecJetMax (double d)
void	setPRecJetMin (double d)
void	setPtRecJetMax (double d)
void	setPtRecJetMin (double d)
void	setRatioMax (double d)
void	setRatioMin (double d)

Protected Member Functions
double	ratio (const reco::Jet &jet, const edm::Handle< reco::SoftLeptonTagInfoCollection > &infos) const
	Finds the ratio of the momentum of any leptons in the jet to jet energy. More...

Protected Attributes
bool	doJetID
double	etaMax
double	etaMin
double	pRecJetMax
double	pRecJetMin
double	ptRecJetMax
double	ptRecJetMin
double	ratioMax
double	ratioMin

Detailed Description

Decide if jet and associated parton satisfy desired kinematic cuts.

Definition at line 14 of file AcceptJet.h.

Constructor & Destructor Documentation

AcceptJet()

AcceptJet::AcceptJet	(	const double &	etaMin_,
		const double &	etaMax_,
		const double &	ptMin_,
		const double &	ptMax_,
		const double &	pMin_,
		const double &	pMax_,
		const double &	ratioMin_,
		const double &	ratioMax_,
		const bool &	doJetID_
	)

Definition at line 10 of file AcceptJet.cc.

    : etaMin(etaMin_),
      etaMax(etaMax_),
      ptRecJetMin(ptMin_),
      ptRecJetMax(ptMax_),
      pRecJetMin(pMin_),
      pRecJetMax(pMax_),
      ratioMin(ratioMin_),
      ratioMax(ratioMax_),
      doJetID(doJetID_) {}

Member Function Documentation

operator()()

bool AcceptJet::operator()	(	const reco::Jet &	jet,
		const int &	jetFlavour,
		const edm::Handle< reco::SoftLeptonTagInfoCollection > &	infos,
		const double	jec
	)		const

Returns true if jet and associated parton satisfy kinematic cuts.

Definition at line 29 of file AcceptJet.cc.

                                                   {
  // temporary fudge to correct for double loop error
  //  jetPartonMomentum /= 2.0;
 
  if (fabs(jet.eta()) < etaMin || fabs(jet.eta()) > etaMax)
    return false;
 
  //   if ( jetFlavour.underlyingParton4Vec().P() < pPartonMin  ||
  //        jetFlavour.underlyingParton4Vec().P() > pPartonMax  ) accept = false;
  //
  //   if ( jetFlavour.underlyingParton4Vec().Pt() < ptPartonMin  ||
  //        jetFlavour.underlyingParton4Vec().Pt() > ptPartonMax  ) accept = false;
 
  if (jet.pt() * jec < ptRecJetMin || jet.pt() * jec > ptRecJetMax)
    return false;
 
  if (jet.p() * jec < pRecJetMin || jet.p() * jec > pRecJetMax)
    return false;
 
  if (!infos.isValid()) {
    LogDebug("infos not valid") << "A valid SoftLeptonTagInfoCollection was not found!"
                                << " Skipping ratio check.";
  } else {
    double pToEratio = ratio(jet, infos);
    if (pToEratio / jec < ratioMin || pToEratio / jec > ratioMax)
      return false;
  }
 
  if (doJetID) {
    const reco::PFJet* pfjet = dynamic_cast<const reco::PFJet*>(&jet);
    if (pfjet) {
      double neutralHadronEnergyFraction = pfjet->neutralHadronEnergy() / (jet.energy() * jec);
      double neutralEmEnergyFraction = pfjet->neutralEmEnergy() / (jet.energy() * jec);
      int nConstituents = pfjet->getPFConstituents().size();
      double chargedHadronEnergyFraction = pfjet->chargedHadronEnergy() / (jet.energy() * jec);
      int chargedMultiplicity = pfjet->chargedMultiplicity();
      double chargedEmEnergyFraction = pfjet->chargedEmEnergy() / (jet.energy() * jec);
      if (!(neutralHadronEnergyFraction < 0.99 && neutralEmEnergyFraction < 0.99 && nConstituents > 1 &&
            chargedHadronEnergyFraction > 0.0 && chargedMultiplicity > 0.0 && chargedEmEnergyFraction < 0.99))
        return false;  //2012 values
    }
    //Looks to not work -> to be fixed
    //else std::cout<<"Jets are not PF jets, put 'doJetID' to False."<<std::endl;
  }
 
  return true;
}

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::btau::chargedHadronEnergyFraction, reco::PFJet::chargedMultiplicity(), doJetID, etaMax, etaMin, reco::PFJet::getPFConstituents(), MessageLogger_cfi::infos, jetMETDQMOfflineSource_cff::jec, metsig::jet, LogDebug, jets_cff::nConstituents, reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), reco::btau::neutralHadronEnergyFraction, pRecJetMax, pRecJetMin, ptRecJetMax, ptRecJetMin, ratio(), and ratioMax.

ratio()

double AcceptJet::ratio ( const reco::Jet &  jet, const edm::Handle< reco::SoftLeptonTagInfoCollection > &  infos  ) const

protected

Finds the ratio of the momentum of any leptons in the jet to jet energy.

Definition at line 80 of file AcceptJet.cc.

                                                                                                         {
  double jetRatio = 0.0;
  reco::SoftLeptonTagInfoCollection::const_iterator infoiter = infos->begin();
  for (; infoiter != infos->end(); ++infoiter) {
    if (reco::deltaR(jet.eta(), jet.phi(), infoiter->jet()->eta(), infoiter->jet()->phi()) > 1e-4)
      continue;
 
    if (infoiter->leptons() == 0)
      break;
 
    for (unsigned int k = 0; k != infoiter->leptons(); ++k) {
      double tempRatio = infoiter->properties(k).ratio;
      if (tempRatio > jetRatio)
        jetRatio = tempRatio;
    }
  }
 
  return jetRatio;
}

References reco::deltaR(), MillePedeFileConverter_cfg::e, MessageLogger_cfi::infos, metsig::jet, and dqmdumpme::k.

Referenced by operator()().

setDoJetID()

void AcceptJet::setDoJetID ( bool  b )

inline

Definition at line 42 of file AcceptJet.h.

{ doJetID = b; }

References b, and doJetID.

setEtaMax()

void AcceptJet::setEtaMax ( double  d )

inline

Definition at line 33 of file AcceptJet.h.

{ etaMax = d; }

References ztail::d, and etaMax.

setEtaMin()

void AcceptJet::setEtaMin ( double  d )

inline

Set cut parameters.

Definition at line 32 of file AcceptJet.h.

{ etaMin = d; }

References ztail::d, and etaMin.

setPRecJetMax()

void AcceptJet::setPRecJetMax ( double  d )

inline

Definition at line 39 of file AcceptJet.h.

{ pRecJetMax = d; }

References ztail::d, and pRecJetMax.

setPRecJetMin()

void AcceptJet::setPRecJetMin ( double  d )

inline

Definition at line 38 of file AcceptJet.h.

{ pRecJetMin = d; }

References ztail::d, and pRecJetMin.

setPtRecJetMax()

void AcceptJet::setPtRecJetMax ( double  d )

inline

Definition at line 37 of file AcceptJet.h.

{ ptRecJetMax = d; }

References ztail::d, and ptRecJetMax.

setPtRecJetMin()

void AcceptJet::setPtRecJetMin ( double  d )

inline

Definition at line 36 of file AcceptJet.h.

{ ptRecJetMin = d; }

References ztail::d, and ptRecJetMin.

setRatioMax()

void AcceptJet::setRatioMax ( double  d )

inline

Definition at line 41 of file AcceptJet.h.

{ ratioMax = d; }

References ztail::d, and ratioMax.

setRatioMin()

void AcceptJet::setRatioMin ( double  d )

inline

Definition at line 40 of file AcceptJet.h.

{ ratioMin = d; }

References ztail::d, and ratioMin.

Member Data Documentation

doJetID

bool AcceptJet::doJetID

protected

Definition at line 74 of file AcceptJet.h.

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

etaMax

double AcceptJet::etaMax

protected

Definition at line 53 of file AcceptJet.h.

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

etaMin

double AcceptJet::etaMin

protected

Definition at line 52 of file AcceptJet.h.

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

pRecJetMax

double AcceptJet::pRecJetMax

protected

Definition at line 68 of file AcceptJet.h.

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

pRecJetMin

double AcceptJet::pRecJetMin

protected

Definition at line 67 of file AcceptJet.h.

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

ptRecJetMax

double AcceptJet::ptRecJetMax

protected

Definition at line 65 of file AcceptJet.h.

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

ptRecJetMin

double AcceptJet::ptRecJetMin

protected

Definition at line 64 of file AcceptJet.h.

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

ratioMax

double AcceptJet::ratioMax

protected

Definition at line 71 of file AcceptJet.h.

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

ratioMin

double AcceptJet::ratioMin

protected

Definition at line 70 of file AcceptJet.h.

Referenced by setRatioMin().