AntiElectronDeadECAL Class Reference

#include <AntiElectronDeadECAL.h>

Classes
struct	TowerInfo

Public Member Functions
	AntiElectronDeadECAL (const edm::ParameterSet &, edm::ConsumesCollector &&)
void	beginEvent (const edm::EventSetup &)
bool	operator() (const reco::Candidate *tau) const
	~AntiElectronDeadECAL ()

Private Member Functions
void	updateBadTowers (const edm::EventSetup &)

Private Attributes
std::vector< TowerInfo >	badTowers_
const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeometryToken_
edm::ESWatcher< CaloGeometryRecord >	caloGeometryWatcher_
const edm::ESGetToken< EcalChannelStatus, EcalChannelStatusRcd >	channelStatusToken_
edm::ESWatcher< EcalChannelStatusRcd >	channelStatusWatcher_
const double	dR2_
const bool	extrapolateToECalEntrance_
edm::ESWatcher< IdealGeometryRecord >	idealGeometryWatcher_
bool	isFirstEvent_
const unsigned	minStatus_
PositionAtECalEntranceComputer	positionAtECalEntrance_
const edm::ESGetToken< EcalTrigTowerConstituentsMap, IdealGeometryRecord >	ttMapToken_
const int	verbosity_

Static Private Attributes
static const uint16_t	statusMask_ = 0x1F

Detailed Description

Flag tau candidates reconstructed near dead ECAL channels, in order to reduce e -> tau fakes not rejected by anti-e MVA discriminator

The motivation for this flag is this presentation: https://indico.cern.ch/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=177223

Code adapted from: RecoTauTag/RecoTau/plugins/PFRecoTauDiscriminationAgainstElectronDeadECAL.cc

Authors

Lauri Andreas Wendland, Christian Veelken

Definition at line 37 of file AntiElectronDeadECAL.h.

Constructor & Destructor Documentation

AntiElectronDeadECAL()

AntiElectronDeadECAL::AntiElectronDeadECAL ( const edm::ParameterSet &  cfg, edm::ConsumesCollector &&  cc  )

explicit

Definition at line 16 of file AntiElectronDeadECAL.cc.

    : minStatus_(cfg.getParameter<uint32_t>("minStatus")),
      dR2_(std::pow(cfg.getParameter<double>("dR"), 2)),
      extrapolateToECalEntrance_(cfg.getParameter<bool>("extrapolateToECalEntrance")),
      verbosity_(cfg.getParameter<int>("verbosity")),
      channelStatusToken_(cc.esConsumes()),
      caloGeometryToken_(cc.esConsumes()),
      ttMapToken_(cc.esConsumes()),
      positionAtECalEntrance_(PositionAtECalEntranceComputer(cc)),
      isFirstEvent_(true) {}

~AntiElectronDeadECAL()

AntiElectronDeadECAL::~AntiElectronDeadECAL

(

)

Definition at line 27 of file AntiElectronDeadECAL.cc.

{}

Member Function Documentation

beginEvent()

void AntiElectronDeadECAL::beginEvent

(

const edm::EventSetup &

es

)

Definition at line 29 of file AntiElectronDeadECAL.cc.

References PositionAtECalEntranceComputer::beginEvent(), positionAtECalEntrance_, and updateBadTowers().

Referenced by TauDiscriminationAgainstElectronDeadECAL< TauType, TauDiscriminator >::beginEvent().

                                                             {
  updateBadTowers(es);
  positionAtECalEntrance_.beginEvent(es);
}

operator()()

bool AntiElectronDeadECAL::operator()

(

const reco::Candidate *

tau

)

const

Definition at line 94 of file AntiElectronDeadECAL.cc.

References badTowers_, HLTMuonOfflineAnalyzer_cfi::deltaR2, dR2_, reco::Candidate::eta(), extrapolateToECalEntrance_, edm::Ptr< T >::get(), edm::Ptr< T >::isNonnull(), reco::PFTau::leadChargedHadrCand(), pat::Tau::leadChargedHadrCand(), reco::Candidate::pdgId(), reco::Candidate::phi(), positionAtECalEntrance_, reco::Candidate::pt(), summarizeEdmComparisonLogfiles::success, metsig::tau, runTauDisplay::tau_eta, runTauDisplay::tau_phi, and verbosity_.

                                                                    {
  bool isNearBadTower = false;
  double tau_eta = tau->eta();
  double tau_phi = tau->phi();
  const reco::Candidate* leadChargedHadron = nullptr;
  if (extrapolateToECalEntrance_) {
    const reco::PFTau* pfTau = dynamic_cast<const reco::PFTau*>(tau);
    if (pfTau != nullptr) {
      leadChargedHadron = pfTau->leadChargedHadrCand().isNonnull() ? pfTau->leadChargedHadrCand().get() : nullptr;
    } else {
      const pat::Tau* patTau = dynamic_cast<const pat::Tau*>(tau);
      if (patTau != nullptr) {
        leadChargedHadron = patTau->leadChargedHadrCand().isNonnull() ? patTau->leadChargedHadrCand().get() : nullptr;
      }
    }
  }
  if (leadChargedHadron != nullptr) {
    bool success = false;
    reco::Candidate::Point positionAtECalEntrance = positionAtECalEntrance_(leadChargedHadron, success);
    if (success) {
      tau_eta = positionAtECalEntrance.eta();
      tau_phi = positionAtECalEntrance.phi();
    }
  }
  if (verbosity_) {
    edm::LogPrint("TauAgainstEleDeadECAL") << "<AntiElectronDeadECal::operator()>:";
    edm::LogPrint("TauAgainstEleDeadECAL") << " #badTowers = " << badTowers_.size();
    if (leadChargedHadron != nullptr) {
      edm::LogPrint("TauAgainstEleDeadECAL")
          << " leadChargedHadron (" << leadChargedHadron->pdgId() << "): Pt = " << leadChargedHadron->pt()
          << ", eta at ECAL (vtx) = " << tau_eta << " (" << leadChargedHadron->eta() << ")"
          << ", phi at ECAL (vtx) = " << tau_phi << " (" << leadChargedHadron->phi() << ")";
    } else {
      edm::LogPrint("TauAgainstEleDeadECAL")
          << " tau: Pt = " << tau->pt() << ", eta at vtx = " << tau_eta << ", phi at vtx = " << tau_phi;
    }
  }
  for (auto const& badTower : badTowers_) {
    if (deltaR2(badTower.eta_, badTower.phi_, tau_eta, tau_phi) < dR2_) {
      if (verbosity_) {
        edm::LogPrint("TauAgainstEleDeadECAL")
            << " matches badTower: eta = " << badTower.eta_ << ", phi = " << badTower.phi_;
      }
      isNearBadTower = true;
      break;
    }
  }
  return isNearBadTower;
}

updateBadTowers()

void AntiElectronDeadECAL::updateBadTowers ( const edm::EventSetup &  es )

private

Definition at line 65 of file AntiElectronDeadECAL.cc.

References badTowers_, caloGeometryToken_, caloGeometryWatcher_, channelStatusToken_, channelStatusWatcher_, edm::ESWatcher< T >::check(), edm::EventSetup::getData(), idealGeometryWatcher_, isFirstEvent_, submitPVResolutionJobs::key, minStatus_, statusMask_, and ttMapToken_.

Referenced by beginEvent().

                                                                  {
  // NOTE: modified version of SUSY CAF code
  //         UserCode/SusyCAF/plugins/SusyCAF_EcalDeadChannels.cc

  if (!isFirstEvent_ && !channelStatusWatcher_.check(es) && !caloGeometryWatcher_.check(es) &&
      !idealGeometryWatcher_.check(es))
    return;

  auto const& channelStatus = es.getData(channelStatusToken_);
  auto const& caloGeometry = es.getData(caloGeometryToken_);
  auto const& ttMap = es.getData(ttMapToken_);

  std::map<uint32_t, unsigned> nBadCrystals, maxStatus;
  std::map<uint32_t, double> sumEta, sumPhi;

  loopXtals<EBDetId>(
      nBadCrystals, maxStatus, sumEta, sumPhi, &channelStatus, &caloGeometry, &ttMap, minStatus_, statusMask_);
  loopXtals<EEDetId>(
      nBadCrystals, maxStatus, sumEta, sumPhi, &channelStatus, &caloGeometry, &ttMap, minStatus_, statusMask_);

  badTowers_.clear();
  for (auto it : nBadCrystals) {
    uint32_t key = it.first;
    badTowers_.push_back(TowerInfo(key, it.second, maxStatus[key], sumEta[key] / it.second, sumPhi[key] / it.second));
  }

  isFirstEvent_ = false;
}

Member Data Documentation

badTowers_

std::vector<TowerInfo> AntiElectronDeadECAL::badTowers_

private

Definition at line 69 of file AntiElectronDeadECAL.h.

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

caloGeometryToken_

const edm::ESGetToken<CaloGeometry, CaloGeometryRecord> AntiElectronDeadECAL::caloGeometryToken_

private

Definition at line 52 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

caloGeometryWatcher_

edm::ESWatcher<CaloGeometryRecord> AntiElectronDeadECAL::caloGeometryWatcher_

private

Definition at line 73 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

channelStatusToken_

const edm::ESGetToken<EcalChannelStatus, EcalChannelStatusRcd> AntiElectronDeadECAL::channelStatusToken_

private

Definition at line 51 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

channelStatusWatcher_

edm::ESWatcher<EcalChannelStatusRcd> AntiElectronDeadECAL::channelStatusWatcher_

private

Definition at line 72 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

dR2_

const double AntiElectronDeadECAL::dR2_

private

Definition at line 48 of file AntiElectronDeadECAL.h.

Referenced by operator()().

extrapolateToECalEntrance_

const bool AntiElectronDeadECAL::extrapolateToECalEntrance_

private

Definition at line 49 of file AntiElectronDeadECAL.h.

Referenced by operator()().

idealGeometryWatcher_

edm::ESWatcher<IdealGeometryRecord> AntiElectronDeadECAL::idealGeometryWatcher_

private

Definition at line 74 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

isFirstEvent_

bool AntiElectronDeadECAL::isFirstEvent_

private

Definition at line 76 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

minStatus_

const unsigned AntiElectronDeadECAL::minStatus_

private

Definition at line 47 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

positionAtECalEntrance_

PositionAtECalEntranceComputer AntiElectronDeadECAL::positionAtECalEntrance_

private

Definition at line 55 of file AntiElectronDeadECAL.h.

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

statusMask_

const uint16_t AntiElectronDeadECAL::statusMask_ = 0x1F

staticprivate

Definition at line 70 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

ttMapToken_

const edm::ESGetToken<EcalTrigTowerConstituentsMap, IdealGeometryRecord> AntiElectronDeadECAL::ttMapToken_

private

Definition at line 53 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

verbosity_

const int AntiElectronDeadECAL::verbosity_

private

Definition at line 50 of file AntiElectronDeadECAL.h.

Referenced by operator()().