#include <TtFullLepHypothesis.h>

Inheritance diagram for TtFullLepHypothesis:

Public Member Functions
	TtFullLepHypothesis (const edm::ParameterSet &)
	default constructor More...

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

Protected Member Functions
virtual void	buildHypo (edm::Event &evt, const edm::Handle< std::vector< pat::Electron > > &elecs, const edm::Handle< std::vector< pat::Muon > > &mus, const edm::Handle< std::vector< pat::Jet > > &jets, const edm::Handle< std::vector< pat::MET > > &mets, std::vector< int > &match, const unsigned int iComb)=0
	build event hypothesis from the reco objects of a semi-leptonic event More...

virtual void	buildKey ()=0
	build the event hypothesis key More...

reco::CompositeCandidate	hypo ()
	return event hypothesis More...

bool	isValid (const int &idx, const edm::Handle< std::vector< pat::Jet > > &jets)
	check if index is in valid range of selected jets More...

int	key () const
	return key More...

template<typename C >
std::unique_ptr< reco::ShallowClonePtrCandidate >	makeCandidate (const edm::Handle< C > &handle, const int &idx)
	use one object in a collection to set a ShallowClonePtrCandidate More...

std::unique_ptr< reco::ShallowClonePtrCandidate >	makeCandidate (const edm::Handle< std::vector< pat::Jet > > &handle, const int &idx, const std::string &correctionLevel)
	use one object in a jet collection to set a ShallowClonePtrCandidate with proper jet corrections More...

void	produce (edm::Event &, const edm::EventSetup &) override
	produce the event hypothesis as CompositeCandidate and Key More...

void	resetCandidates ()
	reset candidate pointers before hypo build process More...

Protected Attributes
std::unique_ptr< reco::ShallowClonePtrCandidate >	b_

std::unique_ptr< reco::ShallowClonePtrCandidate >	bBar_

edm::EDGetTokenT< std::vector< pat::Electron > >	elecsToken_

bool	getMatch_

std::string	jetCorrectionLevel_

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

int	key_
	hypothesis key (to be set by the buildKey function) More...

std::unique_ptr< reco::ShallowClonePtrCandidate >	lepton_

std::unique_ptr< reco::ShallowClonePtrCandidate >	leptonBar_

edm::EDGetTokenT< std::vector< std::vector< int > > >	matchToken_
	input label for all necessary collections More...

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

edm::EDGetTokenT< std::vector< pat::Muon > >	musToken_

std::unique_ptr< reco::ShallowClonePtrCandidate >	neutrino_

std::unique_ptr< reco::ShallowClonePtrCandidate >	neutrinoBar_

std::unique_ptr< reco::LeafCandidate >	recNu
	candidates needed for the genmatch hypothesis More...

std::unique_ptr< reco::LeafCandidate >	recNuBar

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 34 of file TtFullLepHypothesis.h.

Constructor & Destructor Documentation

◆ TtFullLepHypothesis()

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

explicit

default constructor

Definition at line 5 of file TtFullLepHypothesis.cc.

References looper::cfg, getMatch_, jetCorrectionLevel_, matchToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : elecsToken_(consumes<std::vector<pat::Electron> >(cfg.getParameter<edm::InputTag>("electrons"))),
       musToken_(consumes<std::vector<pat::Muon> >(cfg.getParameter<edm::InputTag>("muons"))),
       jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
       metsToken_(consumes<std::vector<pat::MET> >(cfg.getParameter<edm::InputTag>("mets"))) {
   getMatch_ = false;
   if (cfg.exists("match")) {
     getMatch_ = true;
     matchToken_ = consumes<std::vector<std::vector<int> > >(cfg.getParameter<edm::InputTag>("match"));
   }
   // if no other correction is given apply L3 (abs) correction
   jetCorrectionLevel_ = "abs";
   if (cfg.exists("jetCorrectionLevel")) {
     jetCorrectionLevel_ = cfg.getParameter<std::string>("jetCorrectionLevel");
   } else {  // if no other correction is given apply L3 (abs) correction
     jetCorrectionLevel_ = "abs";
   }
   produces<std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > >();
   produces<int>("Key");
 }

Member Function Documentation

◆ buildHypo()

virtual void TtFullLepHypothesis::buildHypo	(	edm::Event &	evt,
		const edm::Handle< std::vector< pat::Electron > > &	elecs,
		const edm::Handle< std::vector< pat::Muon > > &	mus,
		const edm::Handle< std::vector< pat::Jet > > &	jets,
		const edm::Handle< std::vector< pat::MET > > &	mets,
		std::vector< int > &	match,
		const unsigned int	iComb
	)

protectedpure virtual

build event hypothesis from the reco objects of a semi-leptonic event

Implemented in TtFullLepHypGenMatch, and TtFullLepHypKinSolution.

Referenced by produce().

◆ buildKey()

virtual void TtFullLepHypothesis::buildKey ( )

protectedpure virtual

build the event hypothesis key

Implemented in TtFullLepHypGenMatch, and TtFullLepHypKinSolution.

Referenced by produce().

◆ hypo()

reco::CompositeCandidate TtFullLepHypothesis::hypo ( )

protected

return event hypothesis

Definition at line 89 of file TtFullLepHypothesis.cc.

Referenced by produce().

                                                  {
   // check for sanity of the hypothesis
   if (!lepton_ || !leptonBar_ || !b_ || !bBar_) {
     return reco::CompositeCandidate();
   }
 
   if (key() == TtFullLeptonicEvent::kGenMatch && (!recNu || !recNuBar)) {
     edm::LogInfo("TtFullHypothesis") << "no neutrinos for gen match" << std::endl;
     return reco::CompositeCandidate();
   }
   if (key() == TtFullLeptonicEvent::kKinSolution && (!neutrino_ || !neutrinoBar_)) {
     edm::LogInfo("TtFullHypothesis") << "no neutrinos for kin solution" << std::endl;
     return reco::CompositeCandidate();
   }
 
   // setup transient references
   reco::CompositeCandidate hyp, Top, WPlus, TopBar, WMinus;
 
   AddFourMomenta addFourMomenta;
 
   // build up the top branch
   WPlus.addDaughter(std::move(leptonBar_), TtFullLepDaughter::LepBar);
   if (key() == TtFullLeptonicEvent::kKinSolution)
     WPlus.addDaughter(std::move(neutrino_), TtFullLepDaughter::Nu);
   else if (key() == TtFullLeptonicEvent::kGenMatch)
     WPlus.addDaughter(std::move(recNu), TtFullLepDaughter::Nu);
   addFourMomenta.set(WPlus);
   Top.addDaughter(WPlus, TtFullLepDaughter::WPlus);
   Top.addDaughter(std::move(b_), TtFullLepDaughter::B);
   addFourMomenta.set(Top);
 
   // build up the anti top branch
   WMinus.addDaughter(std::move(lepton_), TtFullLepDaughter::Lep);
   if (key() == TtFullLeptonicEvent::kKinSolution)
     WMinus.addDaughter(std::move(neutrinoBar_), TtFullLepDaughter::NuBar);
   else if (key() == TtFullLeptonicEvent::kGenMatch)
     WMinus.addDaughter(std::move(recNuBar), TtFullLepDaughter::NuBar);
   addFourMomenta.set(WMinus);
   TopBar.addDaughter(WMinus, TtFullLepDaughter::WMinus);
   TopBar.addDaughter(std::move(bBar_), TtFullLepDaughter::BBar);
   addFourMomenta.set(TopBar);
 
   // build ttbar hypothesis
   hyp.addDaughter(Top, TtFullLepDaughter::Top);
   hyp.addDaughter(TopBar, TtFullLepDaughter::TopBar);
   addFourMomenta.set(hyp);
 
   // the four momentum of the met is not added to the hypothesis
   // because it is allready included through the neutrinos
   //hyp.addDaughter(*met_, TtFullLepDaughter::Met);
   return hyp;
 }

◆ isValid()

bool TtFullLepHypothesis::isValid	(	const int &	idx,
		const edm::Handle< std::vector< pat::Jet > > &	jets
	)

inlineprotected

check if index is in valid range of selected jets

Definition at line 56 of file TtFullLepHypothesis.h.

References heavyIonCSV_trainingSettings::idx, and PDWG_EXODelayedJetMET_cff::jets.

Referenced by ntupleDataFormat._Object::_checkIsValid(), TtFullLepHypGenMatch::buildHypo(), and core.AutoHandle.AutoHandle::ReallyLoad().

                                                                           {
     return (0 <= idx && idx < (int)jets->size());
   };

◆ key()

int TtFullLepHypothesis::key ( ) const

inlineprotected

return key

Definition at line 52 of file TtFullLepHypothesis.h.

References key_.

Referenced by hypo(), and produce().

52 { return key_; };

TtFullLepHypothesis::key_

int key_

hypothesis key (to be set by the buildKey function)

Definition: TtFullLepHypothesis.h:90

◆ makeCandidate() [1/2]

template<typename C >

std::unique_ptr< reco::ShallowClonePtrCandidate > TtFullLepHypothesis::makeCandidate	(	const edm::Handle< C > &	handle,
		const int &	idx
	)

protected

use one object in a collection to set a ShallowClonePtrCandidate

Definition at line 109 of file TtFullLepHypothesis.h.

References patZpeak::handle, and heavyIonCSV_trainingSettings::idx.

Referenced by TtFullLepHypKinSolution::buildHypo(), and TtFullLepHypGenMatch::buildHypo().

                                                                                                    {
   typedef typename C::value_type O;
   edm::Ptr<O> ptr = edm::Ptr<O>(handle, idx);
   return std::make_unique<reco::ShallowClonePtrCandidate>(ptr, ptr->charge(), ptr->p4(), ptr->vertex());
 }

◆ makeCandidate() [2/2]

std::unique_ptr< reco::ShallowClonePtrCandidate > TtFullLepHypothesis::makeCandidate	(	const edm::Handle< std::vector< pat::Jet > > &	handle,
		const int &	idx,
		const std::string &	correctionLevel
	)

protected

use one object in a jet collection to set a ShallowClonePtrCandidate with proper jet corrections

Definition at line 143 of file TtFullLepHypothesis.cc.

References reco::LeafCandidate::charge(), pat::Jet::correctedJet(), patZpeak::handle, heavyIonCSV_trainingSettings::idx, jetCorrectionLevel_, reco::LeafCandidate::p4(), and reco::LeafCandidate::vertex().

                                                                                                    {
   edm::Ptr<pat::Jet> ptr = edm::Ptr<pat::Jet>(handle, idx);
   return std::make_unique<reco::ShallowClonePtrCandidate>(
       ptr, ptr->charge(), ptr->correctedJet(jetCorrectionLevel_, "bottom").p4(), ptr->vertex());
 }

◆ produce()

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

overrideprotected

produce the event hypothesis as CompositeCandidate and Key

Definition at line 27 of file TtFullLepHypothesis.cc.

References buildHypo(), buildKey(), singleTopDQM_cfi::elecs, elecsToken_, edm::Event::getByToken(), getMatch_, hypo(), mps_fire::i, PDWG_EXODelayedJetMET_cff::jets, jetsToken_, key(), match(), matchToken_, singleTopDQM_cfi::mets, metsToken_, eostools::move(), musToken_, edm::Event::put(), resetCandidates(), and trackerHitRTTI::vector.

                                                                            {
   edm::Handle<std::vector<pat::Electron> > elecs;
   evt.getByToken(elecsToken_, elecs);
 
   edm::Handle<std::vector<pat::Muon> > mus;
   evt.getByToken(musToken_, mus);
 
   edm::Handle<std::vector<pat::Jet> > jets;
   evt.getByToken(jetsToken_, jets);
 
   edm::Handle<std::vector<pat::MET> > mets;
   evt.getByToken(metsToken_, mets);
 
   std::vector<std::vector<int> > matchVec;
   if (getMatch_) {
     edm::Handle<std::vector<std::vector<int> > > matchHandle;
     evt.getByToken(matchToken_, matchHandle);
     ;
     matchVec = *matchHandle;
   } else {
     std::vector<int> dummyMatch;
     for (unsigned int i = 0; i < 4; ++i)
       dummyMatch.push_back(-1);
     matchVec.push_back(dummyMatch);
   }
 
   // declare unique_ptr for products
   std::unique_ptr<std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > > pOut(
       new std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > >);
   std::unique_ptr<int> pKey(new int);
 
   // build and feed out key
   buildKey();
   *pKey = key();
   evt.put(std::move(pKey), "Key");
 
   // go through given vector of jet combinations
   unsigned int idMatch = 0;
   typedef std::vector<std::vector<int> >::iterator MatchVecIterator;
   for (MatchVecIterator match = matchVec.begin(); match != matchVec.end(); ++match) {
     // reset pointers
     resetCandidates();
     // build hypothesis
     buildHypo(evt, elecs, mus, jets, mets, *match, idMatch++);
     pOut->push_back(std::make_pair(hypo(), *match));
   }
   // feed out hyps and matches
   evt.put(std::move(pOut));
 }

◆ resetCandidates()

void TtFullLepHypothesis::resetCandidates ( )

protected

reset candidate pointers before hypo build process

Definition at line 78 of file TtFullLepHypothesis.cc.

References b_, bBar_, lepton_, leptonBar_, neutrino_, and neutrinoBar_.

Referenced by produce().

                                           {
   lepton_.reset();
   leptonBar_.reset();
   b_.reset();
   bBar_.reset();
   neutrino_.reset();
   neutrinoBar_.reset();
   //met_        = 0;
 }

Member Data Documentation

◆ b_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::b_

protected

Definition at line 95 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), TtFullLepHypGenMatch::buildHypo(), hypo(), and resetCandidates().

◆ bBar_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::bBar_

protected

Definition at line 96 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), TtFullLepHypGenMatch::buildHypo(), hypo(), and resetCandidates().

◆ elecsToken_

edm::EDGetTokenT<std::vector<pat::Electron> > TtFullLepHypothesis::elecsToken_

protected

Definition at line 82 of file TtFullLepHypothesis.h.

Referenced by produce().

◆ getMatch_

bool TtFullLepHypothesis::getMatch_

protected

internal check whether the match information exists or not, if false a blind dummy match vector will be used internally

Definition at line 79 of file TtFullLepHypothesis.h.

Referenced by produce(), and TtFullLepHypothesis().

◆ jetCorrectionLevel_

std::string TtFullLepHypothesis::jetCorrectionLevel_

protected

specify the desired jet correction level (the default should be L3Absolute-'abs')

Definition at line 88 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), TtFullLepHypGenMatch::buildHypo(), makeCandidate(), and TtFullLepHypothesis().

◆ jetsToken_

edm::EDGetTokenT<std::vector<pat::Jet> > TtFullLepHypothesis::jetsToken_

protected

Definition at line 84 of file TtFullLepHypothesis.h.

Referenced by produce().

◆ key_

int TtFullLepHypothesis::key_

protected

hypothesis key (to be set by the buildKey function)

Definition at line 90 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildKey(), TtFullLepHypGenMatch::buildKey(), and key().

◆ lepton_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::lepton_

protected

candidates for internal use for the creation of the hypothesis candidate

Definition at line 93 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), TtFullLepHypGenMatch::buildHypo(), hypo(), and resetCandidates().

◆ leptonBar_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::leptonBar_

protected

Definition at line 94 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), TtFullLepHypGenMatch::buildHypo(), hypo(), and resetCandidates().

◆ matchToken_

edm::EDGetTokenT<std::vector<std::vector<int> > > TtFullLepHypothesis::matchToken_

protected

input label for all necessary collections

Definition at line 81 of file TtFullLepHypothesis.h.

Referenced by produce(), and TtFullLepHypothesis().

◆ metsToken_

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

protected

Definition at line 85 of file TtFullLepHypothesis.h.

Referenced by produce().

◆ musToken_

edm::EDGetTokenT<std::vector<pat::Muon> > TtFullLepHypothesis::musToken_

protected

Definition at line 83 of file TtFullLepHypothesis.h.

Referenced by produce().

◆ neutrino_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::neutrino_

protected

Definition at line 97 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), hypo(), and resetCandidates().

◆ neutrinoBar_

std::unique_ptr<reco::ShallowClonePtrCandidate> TtFullLepHypothesis::neutrinoBar_

protected

Definition at line 98 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypKinSolution::buildHypo(), hypo(), and resetCandidates().

◆ recNu

std::unique_ptr<reco::LeafCandidate> TtFullLepHypothesis::recNu

protected

candidates needed for the genmatch hypothesis

Definition at line 102 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypGenMatch::buildMatchingNeutrinos(), and hypo().

◆ recNuBar

std::unique_ptr<reco::LeafCandidate> TtFullLepHypothesis::recNuBar

protected

Definition at line 103 of file TtFullLepHypothesis.h.

Referenced by TtFullLepHypGenMatch::buildMatchingNeutrinos(), and hypo().

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TtFullLepHypothesis()

Member Function Documentation

◆ buildHypo()

◆ buildKey()

◆ hypo()

◆ isValid()

◆ key()

◆ makeCandidate() [1/2]

◆ makeCandidate() [2/2]

◆ produce()

◆ resetCandidates()

Member Data Documentation

◆ b_

◆ bBar_

◆ elecsToken_

◆ getMatch_

◆ jetCorrectionLevel_

◆ jetsToken_

◆ key_

◆ lepton_

◆ leptonBar_

◆ matchToken_

◆ metsToken_

◆ musToken_

◆ neutrino_

◆ neutrinoBar_

◆ recNu

◆ recNuBar