AnomalousTower.cc

Go to the documentation of this file.

#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "RecoJets/JetProducers/interface/AnomalousTower.h"
#include "DataFormats/CaloTowers/interface/CaloTower.h"
 
#define init_param(type, varname) varname(ps.getParameter<type>(#varname))
 
AnomalousTower::AnomalousTower(const edm::ParameterSet& ps)
    : init_param(unsigned, maxBadEcalCells),
      init_param(unsigned, maxRecoveredEcalCells),
      init_param(unsigned, maxProblematicEcalCells),
      init_param(unsigned, maxBadHcalCells),
      init_param(unsigned, maxRecoveredHcalCells),
      init_param(unsigned, maxProblematicHcalCells) {}
 
bool AnomalousTower::operator()(const reco::Candidate& input) const {
  const CaloTower* tower = dynamic_cast<const CaloTower*>(&input);
  if (tower)
    return tower->numBadEcalCells() > maxBadEcalCells || tower->numRecoveredEcalCells() > maxRecoveredEcalCells ||
           tower->numProblematicEcalCells() > maxProblematicEcalCells || tower->numBadHcalCells() > maxBadHcalCells ||
           tower->numRecoveredHcalCells() > maxRecoveredHcalCells ||
           tower->numProblematicHcalCells() > maxProblematicHcalCells;
  else
    return false;
}