pat::CaloJetSelector Class Reference

Selects good Jets. More...

#include "PhysicsTools/PatUtils/CaloJetSelector.h"

Public Member Functions
	CaloJetSelector (const JetSelection &config)
const ParticleStatus	filter (const reco::CaloJet &Jet) const
	~CaloJetSelector ()

Private Attributes
JetSelection	config_

Detailed Description

Selects good Jets.

The calo jet selector returns a flag (see pat::ParticleStatus) based on one of the possible selections. It is called by the generic JetSelector in case of "custom" selection.

Author

C. Autermann (Uni Hamburg)

Version

Id

CaloJetSelector.h,v 1.4 2008/03/05 14:51:02 fronga Exp

Definition at line 23 of file CaloJetSelector.h.

Constructor & Destructor Documentation

pat::CaloJetSelector::CaloJetSelector ( const JetSelection &  config )

inline

Definition at line 25 of file CaloJetSelector.h.

: config_(config) {}

pat::CaloJetSelector::~CaloJetSelector ( )

inline

Definition at line 26 of file CaloJetSelector.h.

References filter().

{}

Member Function Documentation

const pat::ParticleStatus CaloJetSelector::filter

(

const reco::CaloJet &

Jet

)

const

Returns 0 if Jet matches criteria, a flag otherwise. Criteria depend on the selector's configuration.

Retrieve information Pt Jet

eta region

electromagnetic fraction

(EMCalEnergyFraction + HadCalEnergyFraction) / (EMCalEnergyFraction - HadCalEnergyFraction)

n90: number of towers containing 90% of the jet's energy

Tower Number

Highest Et Tower / Et Jet

Sum(E Twr * DeltaR(Twr-Jet)) / Sum(E Twr)

Pt Jet / Towers Area

Highest Et Tower / Towers Area

Definition at line 8 of file CaloJetSelector.cc.

References pat::BAD, config_, PbPb_ZMuSkimMuonDPG_cff::deltaR, reco::CaloJet::emEnergyFraction(), pat::JetSelection::EMFmax, pat::JetSelection::EMvsHadFmax, PVValHelper::eta, pat::JetSelection::Etamax, reco::CaloJet::getCaloConstituents(), pat::GOOD, pat::JetSelection::HighestTowerOverJetmax, reco::CaloJet::n90(), pat::JetSelection::N90max, pat::JetSelection::NCaloTowersmax, pat::JetSelection::NCaloTowersmin, reco::LeafCandidate::p4(), pat::JetSelection::PtJetOverArea_max, pat::JetSelection::Ptmin, pat::JetSelection::PtTowerOverArea_max, mps_fire::result, pat::JetSelection::RWidthmax, and reco::CaloJet::towersArea().

Referenced by ~CaloJetSelector().

                                  {
  ParticleStatus result = GOOD;

  if (Jet.p4().Pt() < config_.Ptmin)
    result = BAD;

  double eta = fabs(Jet.p4().Eta());
  if (eta > config_.Etamax)
    result = BAD;

  double EMF = Jet.emEnergyFraction();
  if (EMF < config_.EMFmin || EMF > config_.EMFmax)
    result = BAD;

  double HadF = Jet.emEnergyFraction();
  double EMvsHadF = 0.;
  if (EMF - HadF != 0.)
    EMvsHadF = (EMF + HadF) / (EMF - HadF);
  if (EMvsHadF < config_.EMvsHadFmin || EMvsHadF > config_.EMvsHadFmax)
    result = BAD;

  //ratio Energy over Momentum (both from calorimeter)
  //Useful? Both come from a lorentz-vector
  //double EoverP = 0.;
  //if (Jet.p4().P()!=0.) EoverP = Jet.energy() / Jet.p4().P();
  //if (EoverP > config_.EoverPmax) result = BAD;

  double n90 = Jet.n90();
  if (n90 < config_.N90min || n90 > config_.N90max)
    result = BAD;

  std::vector<CaloTowerPtr> jetTowers = Jet.getCaloConstituents();
  if (jetTowers.size() < config_.NCaloTowersmin || jetTowers.size() > config_.NCaloTowersmax)
    result = BAD;

  //calculate tower related variables:
  double MaxEnergyTower = 0., SumTowPt = 0., SumTowPtR = 0.;
  for (std::vector<CaloTowerPtr>::const_iterator tow = jetTowers.begin(), towend = jetTowers.end(); tow != towend;
       ++tow) {
    SumTowPt += (*tow)->et();
    SumTowPtR += (*tow)->et() * deltaR(Jet.p4().Eta(), Jet.p4().Phi(), (*tow)->eta(), (*tow)->phi());
    if ((*tow)->et() > MaxEnergyTower)
      MaxEnergyTower = (*tow)->et();
  }

  double EtTowerMaxOverEtJet = 0.;
  if (Jet.p4().Et() != 0.)
    EtTowerMaxOverEtJet = MaxEnergyTower / Jet.p4().Et();
  if (EtTowerMaxOverEtJet < config_.HighestTowerOverJetmin || EtTowerMaxOverEtJet > config_.HighestTowerOverJetmax)
    result = BAD;

  double RWidth = 0.;
  if (SumTowPt != 0.)
    RWidth = SumTowPtR / SumTowPt;
  if (RWidth < config_.RWidthmin || RWidth > config_.RWidthmax)
    result = BAD;

  double PtJetoverArea = 0.;
  if (Jet.towersArea() != 0.)
    PtJetoverArea = Jet.p4().Pt() / Jet.towersArea();
  if (PtJetoverArea < config_.PtJetOverArea_min || PtJetoverArea > config_.PtJetOverArea_max)
    result = BAD;

  double PtToweroverArea = 0.;
  if (Jet.towersArea() != 0.)
    PtToweroverArea = MaxEnergyTower / Jet.towersArea();
  if (PtToweroverArea < config_.PtTowerOverArea_min || PtToweroverArea > config_.PtTowerOverArea_max)
    result = BAD;

  return result;
}

Member Data Documentation

JetSelection pat::CaloJetSelector::config_

private

Definition at line 33 of file CaloJetSelector.h.

Referenced by filter().