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

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 16 of file AcceptJet.cc.

References reco::PFJet::chargedEmEnergy(), reco::PFJet::chargedHadronEnergy(), reco::btau::chargedHadronEnergyFraction, reco::PFJet::chargedMultiplicity(), doJetID, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), etaMax, etaMin, reco::PFJet::getPFConstituents(), edm::HandleBase::isValid(), patTestJEC_cfi::jec, metsig::jet, LogDebug, jets_cff::nConstituents, reco::PFJet::neutralEmEnergy(), reco::PFJet::neutralHadronEnergy(), reco::btau::neutralHadronEnergyFraction, reco::LeafCandidate::p(), pRecJetMax, pRecJetMin, reco::LeafCandidate::pt(), ptRecJetMax, ptRecJetMin, ratio(), ratioMax, and ratioMin.

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

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 71 of file AcceptJet.cc.

References reco::deltaR(), MillePedeFileConverter_cfg::e, reco::LeafCandidate::eta(), gen::k, and reco::LeafCandidate::phi().

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

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