Inheritance diagram for PFRecoTauDiscriminationAgainstElectronDeadECAL:

Classes
struct	towerInfo

Public Member Functions
void	beginEvent (const edm::Event &evt, const edm::EventSetup &es) override

double	discriminate (const PFTauRef &pfTau) const override

	PFRecoTauDiscriminationAgainstElectronDeadECAL (const edm::ParameterSet &cfg)

	~PFRecoTauDiscriminationAgainstElectronDeadECAL () override

Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
virtual double	discriminate (const TauRef &tau) const =0

virtual void	endEvent (edm::Event &)

void	produce (edm::Event &, const edm::EventSetup &) override

	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)

	TauDiscriminationProducerBase ()

	~TauDiscriminationProducerBase () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
static void	fillProducerDescriptions (edm::ParameterSetDescription &desc)

Private Types
typedef ROOT::Math::LorentzVector< ROOT::Math::PtEtaPhiE4D< double > >	PolarLorentzVector

Private Member Functions
template<class Id >
void	loopXtals (std::map< uint32_t, unsigned > &nBadCrystals, std::map< uint32_t, unsigned > &maxStatus, std::map< uint32_t, double > &sumEta, std::map< uint32_t, double > &sumPhi, const EcalChannelStatus channelStatus, const CaloGeometry caloGeometry, const EcalTrigTowerConstituentsMap *ttMap) const

void	updateBadTowers (const edm::EventSetup &es)

Private Attributes
std::vector< towerInfo >	badTowers_

uint32_t	caloGeometryId_cache_

uint32_t	channelStatusId_cache_

double	dR_

uint32_t	idealGeometryId_cache_

bool	isFirstEvent_

unsigned	minStatus_

std::string	moduleLabel_

int	verbosity_

Static Private Attributes
static const uint16_t	statusMask_ = 0x1F

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
typedef std::vector< TauType >	TauCollection

typedef edm::Ref< TauCollection >	TauRef

typedef edm::RefProd< TauCollection >	TauRefProd

Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Protected Attributes inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator >
std::string	moduleLabel_

double	prediscriminantFailValue_

edm::EDGetTokenT< TauCollection >	Tau_token

size_t	tauIndex_

edm::InputTag	TauProducer_

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

Authors: Lauri Andreas Wendland, Christian Veelken

Definition at line 38 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

Member Typedef Documentation

typedef ROOT::Math::LorentzVector<ROOT::Math::PtEtaPhiE4D<double> > PFRecoTauDiscriminationAgainstElectronDeadECAL::PolarLorentzVector

private

Definition at line 158 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 40 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References edm::ParameterSet::getParameter().

       : PFTauDiscriminationProducerBase(cfg),
         moduleLabel_(cfg.getParameter<std::string>("@module_label")),
         isFirstEvent_(true) {
     minStatus_ = cfg.getParameter<uint32_t>("minStatus");
     dR_ = cfg.getParameter<double>("dR");
 
     verbosity_ = cfg.getParameter<int>("verbosity");
   }

PFRecoTauDiscriminationAgainstElectronDeadECAL::~PFRecoTauDiscriminationAgainstElectronDeadECAL ( )

inlineoverride

Definition at line 49 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

49 {}

Member Function Documentation

void PFRecoTauDiscriminationAgainstElectronDeadECAL::beginEvent	(	const edm::Event &	evt,
		const edm::EventSetup &	es
	)

inlineoverridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 51 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

51 { updateBadTowers(es); }

PFRecoTauDiscriminationAgainstElectronDeadECAL::updateBadTowers

void updateBadTowers(const edm::EventSetup &es)

Definition: PFRecoTauDiscriminationAgainstElectronDeadECAL.cc:81

double PFRecoTauDiscriminationAgainstElectronDeadECAL::discriminate ( const PFTauRef & pfTau ) const

inlineoverride

Definition at line 53 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References reco::deltaR(), PDWG_TauSkim_cff::discriminator, and fillDescriptions().

                                                             {
     if (verbosity_) {
       edm::LogPrint("PFTauAgainstEleDeadECAL") << "<PFRecoTauDiscriminationAgainstElectronDeadECAL::discriminate>:";
       edm::LogPrint("PFTauAgainstEleDeadECAL") << " moduleLabel = " << moduleLabel_;
       edm::LogPrint("PFTauAgainstEleDeadECAL") << "#badTowers = " << badTowers_.size();
       edm::LogPrint("PFTauAgainstEleDeadECAL")
           << "tau: Pt = " << pfTau->pt() << ", eta = " << pfTau->eta() << ", phi = " << pfTau->phi();
     }
     double discriminator = 1.;
     for (std::vector<towerInfo>::const_iterator badTower = badTowers_.begin(); badTower != badTowers_.end();
          ++badTower) {
       if (deltaR(badTower->eta_, badTower->phi_, pfTau->eta(), pfTau->phi()) < dR_) {
         if (verbosity_) {
           edm::LogPrint("PFTauAgainstEleDeadECAL")
               << " matches badTower: eta = " << badTower->eta_ << ", phi = " << badTower->phi_;
         }
         discriminator = 0.;
       }
     }
     if (verbosity_) {
       edm::LogPrint("PFTauAgainstEleDeadECAL") << "--> discriminator = " << discriminator;
     }
     return discriminator;
   }

void PFRecoTauDiscriminationAgainstElectronDeadECAL::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 175 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptional(), DEFINE_FWK_MODULE, HLT_2018_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                                 {
   // pfRecoTauDiscriminationAgainstElectronDeadECAL
   edm::ParameterSetDescription desc;
   desc.add<int>("verbosity", 0);
   {
     edm::ParameterSetDescription pset_Prediscriminants;
     pset_Prediscriminants.add<std::string>("BooleanOperator", "and");
     {
       edm::ParameterSetDescription psd1;
       psd1.add<double>("cut");
       psd1.add<edm::InputTag>("Producer");
       pset_Prediscriminants.addOptional<edm::ParameterSetDescription>("leadTrack", psd1);
     }
     {
       // encountered this at
       // RecoTauTag/Configuration/python/HPSPFTaus_cff.py
       edm::ParameterSetDescription psd1;
       psd1.add<double>("cut");
       psd1.add<edm::InputTag>("Producer");
       pset_Prediscriminants.addOptional<edm::ParameterSetDescription>("decayMode", psd1);
     }
     desc.add<edm::ParameterSetDescription>("Prediscriminants", pset_Prediscriminants);
   }
   desc.add<double>("dR", 0.08);
   desc.add<edm::InputTag>("PFTauProducer", edm::InputTag("pfTauProducer"));
   desc.add<unsigned int>("minStatus", 12);
   descriptions.add("pfRecoTauDiscriminationAgainstElectronDeadECAL", desc);
 }

template<class Id >

void PFRecoTauDiscriminationAgainstElectronDeadECAL::loopXtals	(	std::map< uint32_t, unsigned > &	nBadCrystals,
		std::map< uint32_t, unsigned > &	maxStatus,
		std::map< uint32_t, double > &	sumEta,
		std::map< uint32_t, double > &	sumPhi,
		const EcalChannelStatus *	channelStatus,
		const CaloGeometry *	caloGeometry,
		const EcalTrigTowerConstituentsMap *	ttMap
	)		const

inlineprivate

Definition at line 123 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References EcalCondObjectContainer< T >::end(), PV3DBase< T, PVType, FrameType >::eta(), EcalCondObjectContainer< T >::find(), EcalCondObjectContainer< T >::getMap(), CaloGeometry::getPosition(), mps_fire::i, crabWrapper::key, Max(), PV3DBase< T, PVType, FrameType >::phi(), point, mps_update::status, and EcalTrigTowerConstituentsMap::towerOf().

                                                                   {
     // NOTE: modified version of SUSY CAF code
     //         UserCode/SusyCAF/plugins/SusyCAF_EcalDeadChannels.cc
     for (int i = 0; i < Id::kSizeForDenseIndexing; ++i) {
       Id id = Id::unhashIndex(i);
       if (id == Id(0))
         continue;
       EcalChannelStatusMap::const_iterator it = channelStatus->getMap().find(id.rawId());
       unsigned status = (it == channelStatus->end()) ? 0 : (it->getStatusCode() & statusMask_);
       if (status >= minStatus_) {
         const GlobalPoint& point = caloGeometry->getPosition(id);
         uint32_t key = ttMap->towerOf(id);
         maxStatus[key] = TMath::Max(status, maxStatus[key]);
         ++nBadCrystals[key];
         sumEta[key] += point.eta();
         sumPhi[key] += point.phi();
       }
     }
   }

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

inlineprivate

Definition at line 81 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References edm::EventSetup::get(), crabWrapper::key, and edm::ESHandle< T >::product().

                                                 {
     // NOTE: modified version of SUSY CAF code
     //         UserCode/SusyCAF/plugins/SusyCAF_EcalDeadChannels.cc
     const uint32_t channelStatusId = es.get<EcalChannelStatusRcd>().cacheIdentifier();
     const uint32_t caloGeometryId = es.get<CaloGeometryRecord>().cacheIdentifier();
     const uint32_t idealGeometryId = es.get<IdealGeometryRecord>().cacheIdentifier();
 
     if (!isFirstEvent_ && channelStatusId == channelStatusId_cache_ && caloGeometryId == caloGeometryId_cache_ &&
         idealGeometryId == idealGeometryId_cache_)
       return;
 
     edm::ESHandle<EcalChannelStatus> channelStatus;
     es.get<EcalChannelStatusRcd>().get(channelStatus);
     channelStatusId_cache_ = channelStatusId;
 
     edm::ESHandle<CaloGeometry> caloGeometry;
     es.get<CaloGeometryRecord>().get(caloGeometry);
     caloGeometryId_cache_ = caloGeometryId;
 
     edm::ESHandle<EcalTrigTowerConstituentsMap> ttMap;
     es.get<IdealGeometryRecord>().get(ttMap);
     idealGeometryId_cache_ = idealGeometryId;
 
     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());
     loopXtals<EEDetId>(
         nBadCrystals, maxStatus, sumEta, sumPhi, channelStatus.product(), caloGeometry.product(), ttMap.product());
 
     badTowers_.clear();
     for (std::map<uint32_t, unsigned>::const_iterator it = nBadCrystals.begin(); it != nBadCrystals.end(); ++it) {
       uint32_t key = it->first;
       badTowers_.push_back(
           towerInfo(key, it->second, maxStatus[key], sumEta[key] / it->second, sumPhi[key] / it->second));
     }
 
     isFirstEvent_ = false;
   }