AntiElectronDeadECAL Class Reference

#include <AntiElectronDeadECAL.h>

Classes
struct	TowerInfo

Public Member Functions
	AntiElectronDeadECAL (const edm::ParameterSet &)
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_
edm::ESWatcher< CaloGeometryRecord >	caloGeometryWatcher_
edm::ESWatcher< EcalChannelStatusRcd >	channelStatusWatcher_
const double	dR2_
const bool	extrapolateToECalEntrance_
edm::ESWatcher< IdealGeometryRecord >	idealGeometryWatcher_
bool	isFirstEvent_
const unsigned	minStatus_
PositionAtECalEntranceComputer	positionAtECalEntrance_
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 35 of file AntiElectronDeadECAL.h.

Constructor & Destructor Documentation

AntiElectronDeadECAL()

AntiElectronDeadECAL::AntiElectronDeadECAL ( const edm::ParameterSet &  cfg )

explicit

Definition at line 18 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")),
      isFirstEvent_(true) {}

~AntiElectronDeadECAL()

AntiElectronDeadECAL::~AntiElectronDeadECAL

(

)

Definition at line 25 of file AntiElectronDeadECAL.cc.

{}

Member Function Documentation

beginEvent()

void AntiElectronDeadECAL::beginEvent

(

const edm::EventSetup &

es

)

Definition at line 27 of file AntiElectronDeadECAL.cc.

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

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

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

operator()()

bool AntiElectronDeadECAL::operator()

(

const reco::Candidate *

tau

)

const

Definition at line 111 of file AntiElectronDeadECAL.cc.

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

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_.

updateBadTowers()

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

private

Definition at line 63 of file AntiElectronDeadECAL.cc.

                                                                  {
  // 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;
 
  edm::ESHandle<EcalChannelStatus> channelStatus;
  es.get<EcalChannelStatusRcd>().get(channelStatus);
 
  edm::ESHandle<CaloGeometry> caloGeometry;
  es.get<CaloGeometryRecord>().get(caloGeometry);
 
  edm::ESHandle<EcalTrigTowerConstituentsMap> ttMap;
  es.get<IdealGeometryRecord>().get(ttMap);
 
  std::map<uint32_t, unsigned> nBadCrystals, maxStatus;
  std::map<uint32_t, double> sumEta, sumPhi;
 
  loopXtals<EBDetId>(nBadCrystals,
                     maxStatus,
                     sumEta,
                     sumPhi,
                     channelStatus.product(),
                     caloGeometry.product(),
                     ttMap.product(),
                     minStatus_,
                     statusMask_);
  loopXtals<EEDetId>(nBadCrystals,
                     maxStatus,
                     sumEta,
                     sumPhi,
                     channelStatus.product(),
                     caloGeometry.product(),
                     ttMap.product(),
                     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;
}

References badTowers_, caloGeometryWatcher_, channelStatusWatcher_, edm::ESWatcher< T >::check(), edm::EventSetup::get(), get, idealGeometryWatcher_, isFirstEvent_, crabWrapper::key, minStatus_, edm::ESHandle< T >::product(), and statusMask_.

Referenced by beginEvent().

Member Data Documentation

badTowers_

std::vector<TowerInfo> AntiElectronDeadECAL::badTowers_

private

Definition at line 64 of file AntiElectronDeadECAL.h.

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

caloGeometryWatcher_

edm::ESWatcher<CaloGeometryRecord> AntiElectronDeadECAL::caloGeometryWatcher_

private

Definition at line 68 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

channelStatusWatcher_

edm::ESWatcher<EcalChannelStatusRcd> AntiElectronDeadECAL::channelStatusWatcher_

private

Definition at line 67 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

dR2_

const double AntiElectronDeadECAL::dR2_

private

Definition at line 46 of file AntiElectronDeadECAL.h.

Referenced by operator()().

extrapolateToECalEntrance_

const bool AntiElectronDeadECAL::extrapolateToECalEntrance_

private

Definition at line 47 of file AntiElectronDeadECAL.h.

Referenced by operator()().

idealGeometryWatcher_

edm::ESWatcher<IdealGeometryRecord> AntiElectronDeadECAL::idealGeometryWatcher_

private

Definition at line 69 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

isFirstEvent_

bool AntiElectronDeadECAL::isFirstEvent_

private

Definition at line 71 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

minStatus_

const unsigned AntiElectronDeadECAL::minStatus_

private

Definition at line 45 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

positionAtECalEntrance_

PositionAtECalEntranceComputer AntiElectronDeadECAL::positionAtECalEntrance_

private

Definition at line 50 of file AntiElectronDeadECAL.h.

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

statusMask_

const uint16_t AntiElectronDeadECAL::statusMask_ = 0x1F

staticprivate

Definition at line 65 of file AntiElectronDeadECAL.h.

Referenced by updateBadTowers().

verbosity_

const int AntiElectronDeadECAL::verbosity_

private

Definition at line 48 of file AntiElectronDeadECAL.h.

Referenced by operator()().