Inheritance diagram for PATTauIDEmbedder:

Public Member Functions
	PATTauIDEmbedder (const edm::ParameterSet &)

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

	~PATTauIDEmbedder () 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

Private Types
typedef std::pair< edm::InputTag, std::vector< NameWPIdx > >	IDContainerData

typedef std::pair< std::string, edm::InputTag >	NameTag

typedef std::pair< std::string, int >	NameWPIdx

Private Attributes
size_t	nNewPlainTauIds_

size_t	nNewTauIds_

std::vector< edm::EDGetTokenT< reco::TauDiscriminatorContainer > >	patTauIDContainerTokens_

std::vector< edm::EDGetTokenT< pat::PATTauDiscriminator > >	patTauIDTokens_

edm::EDGetTokenT< pat::TauCollection >	src_

std::vector< std::vector< NameWPIdx > >	tauIDSrcContainers_

std::vector< NameTag >	tauIDSrcs_

Additional Inherited Members
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

Detailed Description

Definition at line 11 of file PATTauIDEmbedder.cc.

Member Typedef Documentation

◆ IDContainerData

typedef std::pair<edm::InputTag, std::vector<NameWPIdx> > PATTauIDEmbedder::IDContainerData

private

Definition at line 24 of file PATTauIDEmbedder.cc.

◆ NameTag

typedef std::pair<std::string, edm::InputTag> PATTauIDEmbedder::NameTag

private

Definition at line 21 of file PATTauIDEmbedder.cc.

◆ NameWPIdx

typedef std::pair<std::string, int> PATTauIDEmbedder::NameWPIdx

private

Definition at line 22 of file PATTauIDEmbedder.cc.

Constructor & Destructor Documentation

◆ PATTauIDEmbedder()

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

explicit

Definition at line 33 of file PATTauIDEmbedder.cc.

                                                              {
   src_ = consumes<pat::TauCollection>(cfg.getParameter<edm::InputTag>("src"));
   // read the different tau ID names
   edm::ParameterSet idps = cfg.getParameter<edm::ParameterSet>("tauIDSources");
   std::vector<std::string> names = idps.getParameterNamesForType<edm::ParameterSet>();
   std::map<std::string, IDContainerData> idContainerMap;
   for (auto const& name : names) {
     edm::ParameterSet idp = idps.getParameter<edm::ParameterSet>(name);
     int wpidx = idp.getParameter<int>("workingPointIndex");
     edm::InputTag tag = idp.getParameter<edm::InputTag>("inputTag");
     if (wpidx == -99) {
       tauIDSrcs_.push_back(NameTag(name, tag));
     } else {
       std::map<std::string, IDContainerData>::iterator it;
       it = idContainerMap.insert({tag.label() + tag.instance(), {tag, std::vector<NameWPIdx>()}}).first;
       it->second.second.push_back(NameWPIdx(name, wpidx));
     }
   }
   // but in any case at least once
   if (tauIDSrcs_.empty() && idContainerMap.empty())
     throw cms::Exception("Configuration") << "PATTauProducer: id addTauID is true, you must specify:\n"
                                           << "\tPSet tauIDSources = { \n"
                                           << "\t\tInputTag <someName> = <someTag>   // as many as you want \n "
                                           << "\t}\n";
   for (auto const& mapEntry : idContainerMap) {
     tauIDSrcContainers_.push_back(mapEntry.second.second);
     patTauIDContainerTokens_.push_back(mayConsume<reco::TauDiscriminatorContainer>(mapEntry.second.first));
   }
   patTauIDTokens_ = edm::vector_transform(
       tauIDSrcs_, [this](NameTag const& tag) { return mayConsume<pat::PATTauDiscriminator>(tag.second); });
   nNewPlainTauIds_ = tauIDSrcs_.size();
   nNewTauIds_ = nNewPlainTauIds_;
   for (auto const& cont : tauIDSrcContainers_) {
     nNewTauIds_ += cont.size();
   }
  
   produces<std::vector<pat::Tau> >();
 }

References looper::cfg, generateEDF::cont, dqmdumpme::first, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNamesForType(), Skims_PA_cff::name, names, nNewPlainTauIds_, nNewTauIds_, patTauIDContainerTokens_, patTauIDTokens_, src_, GlobalPosition_Frontier_DevDB_cff::tag, tauIDSrcContainers_, tauIDSrcs_, and edm::vector_transform().

◆ ~PATTauIDEmbedder()

PATTauIDEmbedder::~PATTauIDEmbedder ( )

inlineoverride

Definition at line 14 of file PATTauIDEmbedder.cc.

14 {};

Member Function Documentation

◆ produce()

void PATTauIDEmbedder::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 72 of file PATTauIDEmbedder.cc.

                                                                      {
   edm::Handle<pat::TauCollection> inputTaus;
   evt.getByToken(src_, inputTaus);
  
   auto outputTaus = std::make_unique<std::vector<pat::Tau> >();
   outputTaus->reserve(inputTaus->size());
  
   int tau_idx = 0;
   for (pat::TauCollection::const_iterator inputTau = inputTaus->begin(); inputTau != inputTaus->end();
        ++inputTau, ++tau_idx) {
     pat::Tau outputTau(*inputTau);
     pat::TauRef inputTauRef(inputTaus, tau_idx);
     size_t nTauIds = inputTau->tauIDs().size();
     std::vector<pat::Tau::IdPair> tauIds(nTauIds + nNewTauIds_);
  
     // copy IDs that are already stored in PAT taus
     for (size_t i = 0; i < nTauIds; ++i) {
       tauIds[i] = inputTau->tauIDs().at(i);
     }
  
     // store IDs that were produced in PATTauDiscriminator format
     edm::Handle<pat::PATTauDiscriminator> tauDiscr;
     for (size_t i = 0; i < nNewPlainTauIds_; ++i) {
       evt.getByToken(patTauIDTokens_[i], tauDiscr);
       tauIds[nTauIds + i].first = tauIDSrcs_[i].first;
       tauIds[nTauIds + i].second = (*tauDiscr)[inputTauRef];
     }
  
     // store IDs that were produced in PATTauDiscriminatorContainer format
     size_t nEmbeddedIDs = nTauIds + nNewPlainTauIds_;
     for (size_t i = 0; i < tauIDSrcContainers_.size(); ++i) {
       auto const& tauDiscrCont = evt.get(patTauIDContainerTokens_[i]);
       for (size_t j = 0; j < tauIDSrcContainers_[i].size(); ++j) {
         tauIds[nEmbeddedIDs + j].first = tauIDSrcContainers_[i][j].first;
         int wpIdx = tauIDSrcContainers_[i][j].second;
         if (wpIdx < 0) {
           if (tauDiscrCont[inputTauRef].rawValues.size() == 1)
             //Only 0th component filled with default value if prediscriminor in PatTauDiscriminator failed.
             tauIds[nEmbeddedIDs + j].second = tauDiscrCont[inputTauRef].rawValues.at(0);
           else
             //uses negative indices to access rawValues. In most cases only one rawValue at WPIdx=-1 exists.
             tauIds[nEmbeddedIDs + j].second = tauDiscrCont[inputTauRef].rawValues.at(-1 - wpIdx);
         } else {
           if (tauDiscrCont[inputTauRef].workingPoints.empty())
             //WP vector not filled if prediscriminor in PatTauDiscriminator failed. Set PAT output to false in this case
             tauIds[nEmbeddedIDs + j].second = 0.0;
           else
             tauIds[nEmbeddedIDs + j].second = tauDiscrCont[inputTauRef].workingPoints.at(wpIdx);
         }
       }
       nEmbeddedIDs += tauIDSrcContainers_[i].size();
     }
  
     outputTau.setTauIDs(tauIds);
     outputTaus->push_back(outputTau);
   }
  
   evt.put(std::move(outputTaus));
 }