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	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () 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 166 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

Constructor & Destructor Documentation

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

inlineexplicit

Definition at line 41 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 51 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

51 {}

Member Function Documentation

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

inlineoverridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator >.

Definition at line 53 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

   {
     updateBadTowers(es);
   }

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

inlineoverride

Definition at line 58 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References reco::deltaR(), mvaPFMET_leptons_cfi::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 184 of file PFRecoTauDiscriminationAgainstElectronDeadECAL.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), edm::ParameterSetDescription::addOptional(), DEFINE_FWK_MODULE, 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 124 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 85 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;
   }