PATTauDiscriminationAgainstElectronMVA6 Class Reference

Inheritance diagram for PATTauDiscriminationAgainstElectronMVA6:

Public Member Functions
void	beginEvent (const edm::Event &, const edm::EventSetup &) override
reco::SingleTauDiscriminatorContainer	discriminate (const TauRef &) const override
	PATTauDiscriminationAgainstElectronMVA6 (const edm::ParameterSet &cfg)
	~PATTauDiscriminationAgainstElectronMVA6 () override
Public Member Functions inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator, TauDiscriminatorDataType, ConsumeType >
virtual void	endEvent (edm::Event &)
void	produce (edm::Event &, const edm::EventSetup &) override
	TauDiscriminationProducerBase ()
	TauDiscriminationProducerBase (const edm::ParameterSet &iConfig)
	~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, TauDiscriminatorDataType, ConsumeType >
static void	fillProducerDescriptions (edm::ParameterSetDescription &desc)
static std::string	getTauTypeString ()
	helper method to retrieve tau type name, e.g. to build correct cfi getter More...

Private Member Functions
bool	isInEcalCrack (double) const

Private Attributes
edm::Handle< pat::ElectronCollection >	Electrons
edm::EDGetTokenT< pat::ElectronCollection >	electronToken
std::string	moduleLabel_
std::unique_ptr< AntiElectronIDMVA6 >	mva_
edm::InputTag	srcElectrons
edm::Handle< TauCollection >	taus_
int	verbosity_
bool	vetoEcalCracks_

Additional Inherited Members
Public Types inherited from TauDiscriminationProducerBase< TauType, TauDiscriminator, TauDiscriminatorDataType, ConsumeType >
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, TauDiscriminatorDataType, ConsumeType >
std::string	moduleLabel_
double	prediscriminantFailValue_
edm::EDGetTokenT< TauCollection >	Tau_token
size_t	tauIndex_
edm::InputTag	TauProducer_

Detailed Description

Definition at line 27 of file PATTauDiscriminationAgainstElectronMVA6.cc.

Constructor & Destructor Documentation

PATTauDiscriminationAgainstElectronMVA6()

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

inlineexplicit

Definition at line 29 of file PATTauDiscriminationAgainstElectronMVA6.cc.

      : PATTauDiscriminationContainerProducerBase(cfg), mva_() {
    mva_ = std::make_unique<AntiElectronIDMVA6>(cfg);
 
    srcElectrons = cfg.getParameter<edm::InputTag>("srcElectrons");
    electronToken = consumes<pat::ElectronCollection>(srcElectrons);
    vetoEcalCracks_ = cfg.getParameter<bool>("vetoEcalCracks");
    verbosity_ = cfg.getParameter<int>("verbosity");
  }

References looper::cfg, and PATTauDiscriminationAgainstElectronMVA6_cfi::srcElectrons.

~PATTauDiscriminationAgainstElectronMVA6()

PATTauDiscriminationAgainstElectronMVA6::~PATTauDiscriminationAgainstElectronMVA6 ( )

inlineoverride

Definition at line 43 of file PATTauDiscriminationAgainstElectronMVA6.cc.

{}

Member Function Documentation

beginEvent()

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

overridevirtual

Reimplemented from TauDiscriminationProducerBase< TauType, TauDiscriminator, TauDiscriminatorDataType, ConsumeType >.

Definition at line 63 of file PATTauDiscriminationAgainstElectronMVA6.cc.

                                                                                                     {
  mva_->beginEvent(evt, es);
 
  evt.getByToken(Tau_token, taus_);
 
  evt.getByToken(electronToken, Electrons);
}

References HLTTauDQMOffline_cfi::Electrons, and edm::Event::getByToken().

discriminate()

reco::SingleTauDiscriminatorContainer PATTauDiscriminationAgainstElectronMVA6::discriminate ( const TauRef &  theTauRef ) const

overridevirtual

Implements TauDiscriminationProducerBase< TauType, TauDiscriminator, TauDiscriminatorDataType, ConsumeType >.

Definition at line 71 of file PATTauDiscriminationAgainstElectronMVA6.cc.

                                   {
  reco::SingleTauDiscriminatorContainer result;
  result.rawValues = {1., -1.};
  double category = -1.;
  bool isGsfElectronMatched = false;
  float deltaRDummy = 9.9;
  const float ECALBarrelEndcapEtaBorder = 1.479;
  float tauEtaAtEcalEntrance = theTauRef->etaAtEcalEntrance();
  float leadChargedPFCandEtaAtEcalEntrance = theTauRef->etaAtEcalEntranceLeadChargedCand();
 
  if ((*theTauRef).leadChargedHadrCand().isNonnull()) {
    int numSignalPFGammaCandsInSigCone = 0;
    const reco::CandidatePtrVector signalGammaCands = theTauRef->signalGammaCands();
    for (const auto& gamma : signalGammaCands) {
      double dR = deltaR(gamma->p4(), theTauRef->leadChargedHadrCand()->p4());
      double signalrad = std::max(0.05, std::min(0.10, 3.0 / std::max(1.0, theTauRef->pt())));
      // gammas inside the tau signal cone
      if (dR < signalrad) {
        numSignalPFGammaCandsInSigCone += 1;
      }
    }
    // loop over the electrons
    for (const auto& theElectron : *Electrons) {
      if (theElectron.pt() > 10.) {  // CV: only take electrons above some minimal energy/Pt into account...
        double deltaREleTau = deltaR(theElectron.p4(), theTauRef->p4());
        deltaRDummy = deltaREleTau;
        if (deltaREleTau < 0.3) {
          double mva_match = mva_->MVAValue(*theTauRef, theElectron);
          bool hasGsfTrack = false;
          pat::PackedCandidate const* packedLeadTauCand =
              dynamic_cast<pat::PackedCandidate const*>(theTauRef->leadChargedHadrCand().get());
          if (abs(packedLeadTauCand->pdgId()) == 11)
            hasGsfTrack = true;
          if (!hasGsfTrack)
            hasGsfTrack = theElectron.gsfTrack().isNonnull();
 
          // veto taus that go to Ecal crack
          if (vetoEcalCracks_ &&
              (isInEcalCrack(tauEtaAtEcalEntrance) || isInEcalCrack(leadChargedPFCandEtaAtEcalEntrance))) {
            // return MVA output value
            result.rawValues.at(0) = -99;
            return result;
          }
          // Veto taus that go to Ecal crack
          if (std::abs(tauEtaAtEcalEntrance) < ECALBarrelEndcapEtaBorder) {  // Barrel
            if (numSignalPFGammaCandsInSigCone == 0 && hasGsfTrack) {
              category = 5.;
            } else if (numSignalPFGammaCandsInSigCone >= 1 && hasGsfTrack) {
              category = 7.;
            }
          } else {  // Endcap
            if (numSignalPFGammaCandsInSigCone == 0 && hasGsfTrack) {
              category = 13.;
            } else if (numSignalPFGammaCandsInSigCone >= 1 && hasGsfTrack) {
              category = 15.;
            }
          }
          result.rawValues.at(0) = std::min(result.rawValues.at(0), float(mva_match));
          isGsfElectronMatched = true;
        }  // deltaR < 0.3
      }    // electron pt > 10
    }      // end of loop over electrons
 
    if (!isGsfElectronMatched) {
      result.rawValues.at(0) = mva_->MVAValue(*theTauRef);
      bool hasGsfTrack = false;
      pat::PackedCandidate const* packedLeadTauCand =
          dynamic_cast<pat::PackedCandidate const*>(theTauRef->leadChargedHadrCand().get());
      if (abs(packedLeadTauCand->pdgId()) == 11)
        hasGsfTrack = true;
 
      // veto taus that go to Ecal crack
      if (vetoEcalCracks_ &&
          (isInEcalCrack(tauEtaAtEcalEntrance) || isInEcalCrack(leadChargedPFCandEtaAtEcalEntrance))) {
        // add category index
        result.rawValues.at(1) = category;
        // return MVA output value
        result.rawValues.at(0) = -99;
        return result;
      }
      // veto taus that go to Ecal crack
      if (std::abs(tauEtaAtEcalEntrance) < ECALBarrelEndcapEtaBorder) {  // Barrel
        if (numSignalPFGammaCandsInSigCone == 0 && !hasGsfTrack) {
          category = 0.;
        } else if (numSignalPFGammaCandsInSigCone >= 1 && !hasGsfTrack) {
          category = 2.;
        }
      } else {  // Endcap
        if (numSignalPFGammaCandsInSigCone == 0 && !hasGsfTrack) {
          category = 8.;
        } else if (numSignalPFGammaCandsInSigCone >= 1 && !hasGsfTrack) {
          category = 10.;
        }
      }
    }
  }
  if (verbosity_) {
    edm::LogPrint("PATTauAgainstEleMVA6") << "<PATTauDiscriminationAgainstElectronMVA6::discriminate>:";
    edm::LogPrint("PATTauAgainstEleMVA6")
        << " tau: Pt = " << theTauRef->pt() << ", eta = " << theTauRef->eta() << ", phi = " << theTauRef->phi();
    edm::LogPrint("PATTauAgainstEleMVA6")
        << " deltaREleTau = " << deltaRDummy << ", isGsfElectronMatched = " << isGsfElectronMatched;
    edm::LogPrint("PATTauAgainstEleMVA6") << " #Prongs = " << theTauRef->signalChargedHadrCands().size();
    edm::LogPrint("PATTauAgainstEleMVA6") << " MVA = " << result.rawValues.at(0) << ", category = " << category;
  }
  // add category index
  result.rawValues.at(1) = category;
  // return MVA output value
  return result;
}

References funct::abs(), taus_updatedMVAIds_cff::category, PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, HLTTauDQMOffline_cfi::Electrons, dqmMemoryStats::float, CustomPhysics_cfi::gamma, edm::Ref< C, T, F >::get(), SiStripPI::max, min(), pat::PackedCandidate::pdgId(), and mps_fire::result.

fillDescriptions()

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

static

Definition at line 188 of file PATTauDiscriminationAgainstElectronMVA6.cc.

                                                                                                         {
  // patTauDiscriminationAgainstElectronMVA6
  edm::ParameterSetDescription desc;
  desc.add<double>("minMVANoEleMatchWOgWOgsfBL", 0.0);
  desc.add<double>("minMVANoEleMatchWgWOgsfBL", 0.0);
  desc.add<bool>("vetoEcalCracks", true);
  desc.add<bool>("usePhiAtEcalEntranceExtrapolation", false);
  desc.add<std::string>("mvaName_wGwGSF_EC", "gbr_wGwGSF_EC");
  desc.add<double>("minMVAWgWgsfBL", 0.0);
  desc.add<std::string>("mvaName_woGwGSF_EC", "gbr_woGwGSF_EC");
  desc.add<double>("minMVAWOgWgsfEC", 0.0);
  desc.add<std::string>("mvaName_wGwGSF_BL", "gbr_wGwGSF_BL");
  desc.add<std::string>("mvaName_woGwGSF_BL", "gbr_woGwGSF_BL");
  desc.add<bool>("returnMVA", true);
  desc.add<bool>("loadMVAfromDB", true);
  {
    edm::ParameterSetDescription psd0;
    psd0.add<std::string>("BooleanOperator", "and");
    {
      edm::ParameterSetDescription psd1;
      psd1.add<double>("cut");
      psd1.add<edm::InputTag>("Producer");
      psd0.addOptional<edm::ParameterSetDescription>("leadTrack", psd1);
    }
    desc.add<edm::ParameterSetDescription>("Prediscriminants", psd0);
  }
  desc.add<std::string>("mvaName_NoEleMatch_woGwoGSF_BL", "gbr_NoEleMatch_woGwoGSF_BL");
  desc.add<edm::InputTag>("srcElectrons", edm::InputTag("slimmedElectrons"));
  desc.add<double>("minMVANoEleMatchWOgWOgsfEC", 0.0);
  desc.add<std::string>("mvaName_NoEleMatch_wGwoGSF_BL", "gbr_NoEleMatch_wGwoGSF_BL");
  desc.add<edm::InputTag>("PATTauProducer", edm::InputTag("slimmedTaus"));
  desc.add<double>("minMVAWOgWgsfBL", 0.0);
  desc.add<double>("minMVAWgWgsfEC", 0.0);
  desc.add<int>("verbosity", 0);
  desc.add<std::string>("mvaName_NoEleMatch_wGwoGSF_EC", "gbr_NoEleMatch_wGwoGSF_EC");
  desc.add<std::string>("method", "BDTG");
  desc.add<std::string>("mvaName_NoEleMatch_woGwoGSF_EC", "gbr_NoEleMatch_woGwoGSF_EC");
  desc.add<double>("minMVANoEleMatchWgWOgsfEC", 0.0);
  descriptions.add("patTauDiscriminationAgainstElectronMVA6", desc);
}

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

isInEcalCrack()

bool PATTauDiscriminationAgainstElectronMVA6::isInEcalCrack ( double  eta ) const

private

Definition at line 183 of file PATTauDiscriminationAgainstElectronMVA6.cc.

                                                                            {
  double absEta = fabs(eta);
  return (absEta > 1.460 && absEta < 1.558);
}

References PVValHelper::eta.

Member Data Documentation

Electrons

edm::Handle<pat::ElectronCollection> PATTauDiscriminationAgainstElectronMVA6::Electrons

private

Definition at line 55 of file PATTauDiscriminationAgainstElectronMVA6.cc.

electronToken

edm::EDGetTokenT<pat::ElectronCollection> PATTauDiscriminationAgainstElectronMVA6::electronToken

private

Definition at line 54 of file PATTauDiscriminationAgainstElectronMVA6.cc.

moduleLabel_

std::string PATTauDiscriminationAgainstElectronMVA6::moduleLabel_

private

Definition at line 50 of file PATTauDiscriminationAgainstElectronMVA6.cc.

Referenced by ExternalGeneratorFilter.ExternalGeneratorFilter::insertInto(), Modules.SwitchProducer::insertInto(), Mixins._TypedParameterizable::insertInto(), Modules.ESSource::nameInProcessDesc_(), and Modules.ESProducer::nameInProcessDesc_().

mva_

std::unique_ptr<AntiElectronIDMVA6> PATTauDiscriminationAgainstElectronMVA6::mva_

private

Definition at line 51 of file PATTauDiscriminationAgainstElectronMVA6.cc.

srcElectrons

edm::InputTag PATTauDiscriminationAgainstElectronMVA6::srcElectrons

private

Definition at line 53 of file PATTauDiscriminationAgainstElectronMVA6.cc.

taus_

edm::Handle<TauCollection> PATTauDiscriminationAgainstElectronMVA6::taus_

private

Definition at line 56 of file PATTauDiscriminationAgainstElectronMVA6.cc.

verbosity_

int PATTauDiscriminationAgainstElectronMVA6::verbosity_

private

Definition at line 60 of file PATTauDiscriminationAgainstElectronMVA6.cc.

vetoEcalCracks_

bool PATTauDiscriminationAgainstElectronMVA6::vetoEcalCracks_

private

Definition at line 58 of file PATTauDiscriminationAgainstElectronMVA6.cc.