#include <TtSemiLepJetCombMVAComputer.h>

Inheritance diagram for TtSemiLepJetCombMVAComputer:

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

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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	produce (edm::Event &evt, const edm::EventSetup &setup) override

Private Attributes
edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_

edm::EDGetTokenT< edm::View< reco::RecoCandidate > >	lepsToken_

const int	maxNComb_

const int	maxNJets_

edm::EDGetTokenT< std::vector< pat::MET > >	metsToken_

PhysicsTools::MVAComputerCache	mvaComputer

edm::ESGetToken< PhysicsTools::Calibration::MVAComputerContainer, TtSemiLepJetCombMVARcd >	mvaToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 21 of file TtSemiLepJetCombMVAComputer.h.

Constructor & Destructor Documentation

◆ TtSemiLepJetCombMVAComputer()

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

explicit

Definition at line 5 of file TtSemiLepJetCombMVAComputer.cc.

     : mvaToken_(esConsumes()),
       lepsToken_(consumes<edm::View<reco::RecoCandidate>>(cfg.getParameter<edm::InputTag>("leps"))),
       jetsToken_(consumes<std::vector<pat::Jet>>(cfg.getParameter<edm::InputTag>("jets"))),
       metsToken_(consumes<std::vector<pat::MET>>(cfg.getParameter<edm::InputTag>("mets"))),
       maxNJets_(cfg.getParameter<int>("maxNJets")),
       maxNComb_(cfg.getParameter<int>("maxNComb")) {
   produces<std::vector<std::vector<int>>>();
   produces<std::vector<double>>("Discriminators");
   produces<std::string>("Method");
   produces<int>("NumberOfConsideredJets");
 }

Member Function Documentation

◆ produce()

void TtSemiLepJetCombMVAComputer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	setup
	)

overrideprivate

Definition at line 18 of file TtSemiLepJetCombMVAComputer.cc.

References PhysicsTools::MVAComputer::eval(), evaluateTtSemiLepJetComb(), edm::Event::getByToken(), mps_fire::i, createfilelist::int, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, lepsToken_, HLT_2023v12_cff::leptons, TtSemiLepEvtPartons::LightQ, TtSemiLepEvtPartons::LightQBar, maxNComb_, maxNJets_, BTaggingMonitor_cfi::met, singleTopDQM_cfi::mets, metsToken_, eostools::move(), mvaComputer, mvaToken_, stdcomb::next_combination(), nPartons, edm::Event::put(), singleTopDQM_cfi::setup, AlCaHLTBitMon_QueryRunRegistry::string, PhysicsTools::MVAComputerCache::update(), contentValuesCheck::values, and trackerHitRTTI::vector.

                                                                                    {
   std::unique_ptr<std::vector<std::vector<int>>> pOut(new std::vector<std::vector<int>>);
   std::unique_ptr<std::vector<double>> pOutDisc(new std::vector<double>);
   std::unique_ptr<std::string> pOutMeth(new std::string);
   std::unique_ptr<int> pJetsConsidered(new int);
 
   const auto& calibContainer = setup.getData(mvaToken_);
   mvaComputer.update(&calibContainer, "ttSemiLepJetCombMVA");
 
   // read name of the processor that provides the MVA discriminator
   // (to be used as meta information)
   std::vector<const PhysicsTools::Calibration::VarProcessor*> processors =
       (calibContainer.find("ttSemiLepJetCombMVA")).getProcessors();
   *pOutMeth = (processors[processors.size() - 3])->getInstanceName();
   evt.put(std::move(pOutMeth), "Method");
 
   // get lepton, jets and mets
   edm::Handle<edm::View<reco::RecoCandidate>> leptons;
   evt.getByToken(lepsToken_, leptons);
 
   edm::Handle<std::vector<pat::Jet>> jets;
   evt.getByToken(jetsToken_, jets);
 
   edm::Handle<std::vector<pat::MET>> mets;
   evt.getByToken(metsToken_, mets);
 
   const unsigned int nPartons = 4;
 
   // skip events with no appropriate lepton candidate,
   // empty METs vector or less jets than partons
   if (leptons->empty() || mets->empty() || jets->size() < nPartons) {
     std::vector<int> invalidCombi;
     for (unsigned int i = 0; i < nPartons; ++i)
       invalidCombi.push_back(-1);
     pOut->push_back(invalidCombi);
     evt.put(std::move(pOut));
     pOutDisc->push_back(0.);
     evt.put(std::move(pOutDisc), "Discriminators");
     *pJetsConsidered = jets->size();
     evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
     return;
   }
 
   const math::XYZTLorentzVector lepton = leptons->begin()->p4();
 
   const pat::MET* met = &(*mets)[0];
 
   // analyze jet combinations
   std::vector<int> jetIndices;
   for (unsigned int i = 0; i < jets->size(); ++i) {
     if (maxNJets_ >= (int)nPartons && maxNJets_ == (int)i) {
       *pJetsConsidered = i;
       break;
     }
     jetIndices.push_back(i);
   }
 
   std::vector<int> combi;
   for (unsigned int i = 0; i < nPartons; ++i)
     combi.push_back(i);
 
   typedef std::pair<double, std::vector<int>> discCombPair;
   std::list<discCombPair> discCombList;
 
   do {
     for (int cnt = 0; cnt < TMath::Factorial(combi.size()); ++cnt) {
       // take into account indistinguishability of the two jets from the hadr. W decay,
       // reduces combinatorics by a factor of 2
       if (combi[TtSemiLepEvtPartons::LightQ] < combi[TtSemiLepEvtPartons::LightQBar]) {
         TtSemiLepJetComb jetComb(*jets, combi, lepton, *met);
 
         // feed MVA input variables into a ValueList
         PhysicsTools::Variable::ValueList values;
         evaluateTtSemiLepJetComb(values, jetComb);
 
         // get discriminator from the MVAComputer
         double discrim = mvaComputer->eval(values);
 
         discCombList.push_back(std::make_pair(discrim, combi));
       }
       next_permutation(combi.begin(), combi.end());
     }
   } while (stdcomb::next_combination(jetIndices.begin(), jetIndices.end(), combi.begin(), combi.end()));
 
   // sort results w.r.t. discriminator values
   discCombList.sort();
 
   // write result into the event
   // (starting with the JetComb having the highest discriminator value -> reverse iterator)
   unsigned int iDiscComb = 0;
   typedef std::list<discCombPair>::reverse_iterator discCombIterator;
   for (discCombIterator discCombPair = discCombList.rbegin(); discCombPair != discCombList.rend(); ++discCombPair) {
     if (maxNComb_ >= 1 && iDiscComb == (unsigned int)maxNComb_)
       break;
     pOut->push_back(discCombPair->second);
     pOutDisc->push_back(discCombPair->first);
     iDiscComb++;
   }
   evt.put(std::move(pOut));
   evt.put(std::move(pOutDisc), "Discriminators");
   evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
 }