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

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () 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, 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));
 }