#include <TtSemiLepHypothesis.h>

Inheritance diagram for TtSemiLepHypothesis:

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

	~TtSemiLepHypothesis () override
	default destructor More...

Public Member Functions inherited from edm::EDProducer
	EDProducer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

std::vector< bool > const &	recordProvenanceList () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
void	buildHypo (const edm::Handle< edm::View< reco::RecoCandidate > > &leps, const edm::Handle< std::vector< pat::MET > > &mets, const edm::Handle< std::vector< pat::Jet > > &jets, std::vector< int > &jetPartonAssociation)
	minimalistic build function for simple hypotheses More...

virtual void	buildHypo (edm::Event &event, const edm::Handle< edm::View< reco::RecoCandidate > > &lepton, const edm::Handle< std::vector< pat::MET > > &neutrino, const edm::Handle< std::vector< pat::Jet > > &jets, std::vector< int > &jetPartonAssociation, 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...

std::string	jetCorrectionLevel (const std::string &quarkType)
	helper function to construct the proper correction level string for corresponding quarkType More...

int	key () const

WDecay::LeptonType	leptonType (const reco::RecoCandidate *cand)
	determine lepton type of reco candidate and return a corresponding WDecay::LeptonType; the type is kNone if it is whether a muon nor an electron 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...

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

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

void	setNeutrino (const edm::Handle< std::vector< pat::MET > > &met, const edm::Handle< edm::View< reco::RecoCandidate > > &leps, const int &idx, const int &type)
	set neutrino, using mW = 80.4 to calculate the neutrino pz More...

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Protected Attributes
bool	getMatch_

reco::ShallowClonePtrCandidate *	hadronicB_

std::string	jetCorrectionLevel_

edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_
	input label for all necessary collections More...

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

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

reco::ShallowClonePtrCandidate *	lepton_

reco::ShallowClonePtrCandidate *	leptonicB_

reco::ShallowClonePtrCandidate *	lightQ_

reco::ShallowClonePtrCandidate *	lightQBar_

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

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

reco::ShallowClonePtrCandidate *	neutrino_

int	neutrinoSolutionType_
	algorithm used to calculate neutrino solutions (see cfi for further details) More...

edm::EDGetTokenT< int >	nJetsConsideredToken_

int	numberOfRealNeutrinoSolutions_

Additional Inherited Members
Public Types inherited from edm::EDProducer
typedef EDProducer	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex > >

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDProducer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Detailed Description

Definition at line 31 of file TtSemiLepHypothesis.h.

Constructor & Destructor Documentation

◆ TtSemiLepHypothesis()

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

explicit

default constructor

Definition at line 10 of file TtSemiLepHypothesis.cc.

     : jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
       lepsToken_(consumes<edm::View<reco::RecoCandidate> >(cfg.getParameter<edm::InputTag>("leps"))),
       metsToken_(consumes<std::vector<pat::MET> >(cfg.getParameter<edm::InputTag>("mets"))),
       nJetsConsideredToken_(consumes<int>(cfg.getParameter<edm::InputTag>("nJetsConsidered"))),
       numberOfRealNeutrinoSolutions_(-1),
       lightQ_(nullptr),
       lightQBar_(nullptr),
       hadronicB_(nullptr),
       leptonicB_(nullptr),
       neutrino_(nullptr),
       lepton_(nullptr) {
   getMatch_ = false;
   if (cfg.exists("match")) {
     getMatch_ = true;
     matchToken_ = consumes<std::vector<std::vector<int> > >(cfg.getParameter<edm::InputTag>("match"));
   }
   if (cfg.exists("neutrinoSolutionType"))
     neutrinoSolutionType_ = cfg.getParameter<int>("neutrinoSolutionType");
   else
     neutrinoSolutionType_ = -1;
   if (cfg.exists("jetCorrectionLevel")) {
     jetCorrectionLevel_ = cfg.getParameter<std::string>("jetCorrectionLevel");
   }
   produces<std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > >();
   produces<int>("Key");
   produces<int>("NumberOfRealNeutrinoSolutions");
   produces<int>("NumberOfConsideredJets");
 }

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

◆ ~TtSemiLepHypothesis()

TtSemiLepHypothesis::~TtSemiLepHypothesis ( )

override

default destructor

Definition at line 41 of file TtSemiLepHypothesis.cc.

                                           {
   if (lightQ_)
     delete lightQ_;
   if (lightQBar_)
     delete lightQBar_;
   if (hadronicB_)
     delete hadronicB_;
   if (leptonicB_)
     delete leptonicB_;
   if (neutrino_)
     delete neutrino_;
   if (lepton_)
     delete lepton_;
 }

References hadronicB_, lepton_, leptonicB_, lightQ_, lightQBar_, and neutrino_.

Member Function Documentation

◆ buildHypo() [1/2]

void TtSemiLepHypothesis::buildHypo	(	const edm::Handle< edm::View< reco::RecoCandidate > > &	leps,
		const edm::Handle< std::vector< pat::MET > > &	mets,
		const edm::Handle< std::vector< pat::Jet > > &	jets,
		std::vector< int > &	jetPartonAssociation
	)

protected

minimalistic build function for simple hypotheses

return key minimalistic build function for simple hypotheses

Definition at line 249 of file TtSemiLepHypothesis.cc.

                                                            {
   // -----------------------------------------------------
   // add jets
   // -----------------------------------------------------
   for (unsigned idx = 0; idx < match.size(); ++idx) {
     if (isValid(match[idx], jets)) {
       switch (idx) {
         case TtSemiLepEvtPartons::LightQ:
           setCandidate(jets, match[idx], lightQ_, jetCorrectionLevel("wQuarkMix"));
           break;
         case TtSemiLepEvtPartons::LightQBar:
           setCandidate(jets, match[idx], lightQBar_, jetCorrectionLevel("wQuarkMix"));
           break;
         case TtSemiLepEvtPartons::HadB:
           setCandidate(jets, match[idx], hadronicB_, jetCorrectionLevel("bQuark"));
           break;
         case TtSemiLepEvtPartons::LepB:
           setCandidate(jets, match[idx], leptonicB_, jetCorrectionLevel("bQuark"));
           break;
       }
     }
   }
  
   // -----------------------------------------------------
   // add lepton
   // -----------------------------------------------------
   if (leps->empty())
     return;
   setCandidate(leps, 0, lepton_);
   match.push_back(0);
  
   // -----------------------------------------------------
   // add neutrino
   // -----------------------------------------------------
   if (mets->empty())
     return;
   if (neutrinoSolutionType_ == -1)
     setCandidate(mets, 0, neutrino_);
   else
     setNeutrino(mets, leps, 0, neutrinoSolutionType_);
 }

References TtSemiLepEvtPartons::HadB, hadronicB_, heavyIonCSV_trainingSettings::idx, isValid(), jetCorrectionLevel(), singleTopDQM_cfi::jets, TtSemiLepEvtPartons::LepB, TtSemiLepHitFitProducer_Electrons_cfi::leps, lepton_, leptonicB_, TtSemiLepEvtPartons::LightQ, lightQ_, TtSemiLepEvtPartons::LightQBar, lightQBar_, match(), singleTopDQM_cfi::mets, neutrino_, neutrinoSolutionType_, setCandidate(), and setNeutrino().

Referenced by TtSemiLepHypGeom::buildHypo(), TtSemiLepHypMaxSumPtWMass::buildHypo(), TtSemiLepHypWMassMaxSumPt::buildHypo(), TtSemiLepHypWMassDeltaTopMass::buildHypo(), TtSemiLepHypMVADisc::buildHypo(), and produce().

◆ buildHypo() [2/2]

virtual void TtSemiLepHypothesis::buildHypo	(	edm::Event &	event,
		const edm::Handle< edm::View< reco::RecoCandidate > > &	lepton,
		const edm::Handle< std::vector< pat::MET > > &	neutrino,
		const edm::Handle< std::vector< pat::Jet > > &	jets,
		std::vector< int > &	jetPartonAssociation,
		const unsigned int	iComb
	)

protectedpure virtual

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

Implemented in TtSemiLepHypGeom, TtSemiLepHypMaxSumPtWMass, TtSemiLepHypMVADisc, TtSemiLepHypWMassDeltaTopMass, TtSemiLepHypWMassMaxSumPt, TtSemiLepHypHitFit, TtSemiLepHypKinFit, and TtSemiLepHypGenMatch.

◆ buildKey()

virtual void TtSemiLepHypothesis::buildKey ( )

protectedpure virtual

build the event hypothesis key

Implemented in TtSemiLepHypHitFit, TtSemiLepHypKinFit, TtSemiLepHypGenMatch, TtSemiLepHypGeom, TtSemiLepHypMaxSumPtWMass, TtSemiLepHypMVADisc, TtSemiLepHypWMassDeltaTopMass, and TtSemiLepHypWMassMaxSumPt.

Referenced by produce().

◆ hypo()

reco::CompositeCandidate TtSemiLepHypothesis::hypo ( )

protected

return event hypothesis

Definition at line 128 of file TtSemiLepHypothesis.cc.

                                                  {
   // check for sanity of the hypothesis
   if (!lightQ_ || !lightQBar_ || !hadronicB_ || !leptonicB_ || !neutrino_ || !lepton_)
     return reco::CompositeCandidate();
  
   // setup transient references
   reco::CompositeCandidate hyp, hadTop, hadW, lepTop, lepW;
  
   AddFourMomenta addFourMomenta;
   // build up the top branch that decays leptonically
   lepW.addDaughter(*lepton_, TtSemiLepDaughter::Lep);
   lepW.addDaughter(*neutrino_, TtSemiLepDaughter::Nu);
   addFourMomenta.set(lepW);
   lepTop.addDaughter(lepW, TtSemiLepDaughter::LepW);
   lepTop.addDaughter(*leptonicB_, TtSemiLepDaughter::LepB);
   addFourMomenta.set(lepTop);
  
   // build up the top branch that decays hadronically
   hadW.addDaughter(*lightQ_, TtSemiLepDaughter::HadP);
   hadW.addDaughter(*lightQBar_, TtSemiLepDaughter::HadQ);
   addFourMomenta.set(hadW);
   hadTop.addDaughter(hadW, TtSemiLepDaughter::HadW);
   hadTop.addDaughter(*hadronicB_, TtSemiLepDaughter::HadB);
   addFourMomenta.set(hadTop);
  
   // build ttbar hypotheses
   hyp.addDaughter(lepTop, TtSemiLepDaughter::LepTop);
   hyp.addDaughter(hadTop, TtSemiLepDaughter::HadTop);
   addFourMomenta.set(hyp);
  
   return hyp;
 }

References reco::CompositeCandidate::addDaughter(), TtSemiLepDaughter::HadB, TtSemiLepDaughter::HadP, TtSemiLepDaughter::HadQ, hadronicB_, TtSemiLepDaughter::HadTop, TtSemiLepDaughter::HadW, TtSemiLepDaughter::Lep, TtSemiLepDaughter::LepB, lepton_, leptonicB_, TtSemiLepDaughter::LepTop, TtSemiLepDaughter::LepW, lightQ_, lightQBar_, neutrino_, TtSemiLepDaughter::Nu, and AddFourMomenta::set().

Referenced by produce().

◆ isValid()

bool TtSemiLepHypothesis::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 69 of file TtSemiLepHypothesis.h.

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

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

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

◆ jetCorrectionLevel()

std::string TtSemiLepHypothesis::jetCorrectionLevel ( const std::string & quarkType )

protected

helper function to construct the proper correction level string for corresponding quarkType

helper function to construct the proper correction level string for corresponding quarkType, for unknown quarkTypes an empty string is returned

Definition at line 174 of file TtSemiLepHypothesis.cc.

                                                                             {
   // jetCorrectionLevel was not configured
   if (jetCorrectionLevel_.empty())
     throw cms::Exception("Configuration")
         << "Unconfigured jetCorrectionLevel. Please use an appropriate, non-empty string.\n";
  
   // quarkType is unknown
   if (!(quarkType == "wQuarkMix" || quarkType == "udsQuark" || quarkType == "cQuark" || quarkType == "bQuark"))
     throw cms::Exception("Configuration") << quarkType << " is unknown as a quarkType for the jetCorrectionLevel.\n";
  
   // combine correction level; start with a ':' even if
   // there is no flavor tag to be added, as it is needed
   // by setCandidate to disentangle the correction tag
   // from a potential flavor tag, which can be empty
   std::string level = jetCorrectionLevel_ + ":";
   if (level == "L5Flavor:" || level == "L6UE:" || level == "L7Parton:") {
     if (quarkType == "wQuarkMix") {
       level += "wMix";
     }
     if (quarkType == "udsQuark") {
       level += "uds";
     }
     if (quarkType == "cQuark") {
       level += "charm";
     }
     if (quarkType == "bQuark") {
       level += "bottom";
     }
   } else {
     level += "none";
   }
   return level;
 }

References jetCorrectionLevel_, personalPlayback::level, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by TtSemiLepHypGenMatch::buildHypo(), and buildHypo().

◆ key()

int TtSemiLepHypothesis::key ( ) const

inlineprotected

Definition at line 65 of file TtSemiLepHypothesis.h.

65 { return key_; };

References key_.

Referenced by produce().

◆ leptonType()

WDecay::LeptonType TtSemiLepHypothesis::leptonType ( const reco::RecoCandidate * cand )

protected

determine lepton type of reco candidate and return a corresponding WDecay::LeptonType; the type is kNone if it is whether a muon nor an electron

Definition at line 162 of file TtSemiLepHypothesis.cc.

                                                                               {
   // check whetherwe are dealing with a reco muon or a reco electron
   WDecay::LeptonType type = WDecay::kNone;
   if (dynamic_cast<const reco::Muon*>(cand)) {
     type = WDecay::kMuon;
   } else if (dynamic_cast<const reco::GsfElectron*>(cand)) {
     type = WDecay::kElec;
   }
   return type;
 }

References WDecay::kElec, WDecay::kMuon, and WDecay::kNone.

Referenced by TtSemiLepHypGenMatch::findMatchingLepton(), and setNeutrino().

◆ produce()

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

overrideprotectedvirtual

produce the event hypothesis as CompositeCandidate and Key

Implements edm::EDProducer.

Definition at line 57 of file TtSemiLepHypothesis.cc.

                                                                            {
   edm::Handle<std::vector<pat::Jet> > jets;
   evt.getByToken(jetsToken_, jets);
  
   edm::Handle<edm::View<reco::RecoCandidate> > leps;
   evt.getByToken(lepsToken_, leps);
  
   edm::Handle<std::vector<pat::MET> > mets;
   evt.getByToken(metsToken_, mets);
  
   edm::Handle<int> nJetsConsidered;
   evt.getByToken(nJetsConsideredToken_, nJetsConsidered);
  
   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);
   std::unique_ptr<int> pNeutrinoSolutions(new int);
   std::unique_ptr<int> pJetsConsidered(new int);
  
   // 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, leps, mets, jets, *match, idMatch++);
     pOut->push_back(std::make_pair(hypo(), *match));
   }
   // feed out hyps and matches
   evt.put(std::move(pOut));
  
   // build and feed out key
   buildKey();
   *pKey = key();
   evt.put(std::move(pKey), "Key");
  
   // feed out number of real neutrino solutions
   *pNeutrinoSolutions = numberOfRealNeutrinoSolutions_;
   evt.put(std::move(pNeutrinoSolutions), "NumberOfRealNeutrinoSolutions");
  
   // feed out number of considered jets
   *pJetsConsidered = *nJetsConsidered;
   evt.put(std::move(pJetsConsidered), "NumberOfConsideredJets");
 }

References buildHypo(), buildKey(), edm::Event::getByToken(), getMatch_, hypo(), mps_fire::i, singleTopDQM_cfi::jets, jetsToken_, key(), TtSemiLepHitFitProducer_Electrons_cfi::leps, lepsToken_, match(), matchToken_, singleTopDQM_cfi::mets, metsToken_, eostools::move(), TtSemiLepHypGenMatch_cfi::nJetsConsidered, nJetsConsideredToken_, numberOfRealNeutrinoSolutions_, L1TrackObjectNtupleMaker_cfg::pOut, edm::Event::put(), resetCandidates(), and trackerHitRTTI::vector.

◆ resetCandidates()

void TtSemiLepHypothesis::resetCandidates ( )

protected

reset candidate pointers before hypo build process

Definition at line 117 of file TtSemiLepHypothesis.cc.

                                           {
   numberOfRealNeutrinoSolutions_ = -1;
   lightQ_ = nullptr;
   lightQBar_ = nullptr;
   hadronicB_ = nullptr;
   leptonicB_ = nullptr;
   neutrino_ = nullptr;
   lepton_ = nullptr;
 }

References hadronicB_, lepton_, leptonicB_, lightQ_, lightQBar_, neutrino_, and numberOfRealNeutrinoSolutions_.

Referenced by produce().

◆ setCandidate() [1/2]

template<typename C >

void TtSemiLepHypothesis::setCandidate	(	const edm::Handle< C > &	handle,
		const int &	idx,
		reco::ShallowClonePtrCandidate *&	clone
	)

protected

use one object in a collection to set a ShallowClonePtrCandidate

Definition at line 123 of file TtSemiLepHypothesis.h.

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

References clone(), patZpeak::handle, and heavyIonCSV_trainingSettings::idx.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), and buildHypo().

◆ setCandidate() [2/2]

void TtSemiLepHypothesis::setCandidate	(	const edm::Handle< std::vector< pat::Jet > > &	handle,
		const int &	idx,
		reco::ShallowClonePtrCandidate *&	clone,
		const std::string &	correctionLevel
	)

protected

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

Definition at line 209 of file TtSemiLepHypothesis.cc.

                                                                          {
   edm::Ptr<pat::Jet> ptr = edm::Ptr<pat::Jet>(handle, idx);
   // disentangle the correction from the potential flavor tag
   // by the separating ':'; the flavor tag can be empty though
   std::string step = correctionLevel.substr(0, correctionLevel.find(':'));
   std::string flavor = correctionLevel.substr(1 + correctionLevel.find(':'));
   float corrFactor = 1.;
   if (flavor == "wMix")
     corrFactor = 0.75 * ptr->jecFactor(step, "uds") + 0.25 * ptr->jecFactor(step, "charm");
   else
     corrFactor = ptr->jecFactor(step, flavor);
   clone = new reco::ShallowClonePtrCandidate(ptr, ptr->charge(), ptr->p4() * corrFactor, ptr->vertex());
 }

References reco::LeafCandidate::charge(), clone(), patZpeak::handle, heavyIonCSV_trainingSettings::idx, pat::Jet::jecFactor(), reco::LeafCandidate::p4(), AlCaHLTBitMon_QueryRunRegistry::string, and reco::LeafCandidate::vertex().

◆ setNeutrino()

void TtSemiLepHypothesis::setNeutrino	(	const edm::Handle< std::vector< pat::MET > > &	met,
		const edm::Handle< edm::View< reco::RecoCandidate > > &	leps,
		const int &	idx,
		const int &	type
	)

protected

set neutrino, using mW = 80.4 to calculate the neutrino pz

Definition at line 227 of file TtSemiLepHypothesis.cc.

                                                        {
   edm::Ptr<pat::MET> ptr = edm::Ptr<pat::MET>(met, 0);
   MEzCalculator mez;
   mez.SetMET(*(met->begin()));
   if (leptonType(&(leps->front())) == WDecay::kMuon)
     mez.SetLepton((*leps)[idx], true);
   else if (leptonType(&(leps->front())) == WDecay::kElec)
     mez.SetLepton((*leps)[idx], false);
   else
     throw cms::Exception("UnimplementedFeature")
         << "Type of lepton given together with MET for solving neutrino kinematics is neither muon nor electron.\n";
   double pz = mez.Calculate(type);
   numberOfRealNeutrinoSolutions_ = mez.IsComplex() ? 0 : 2;
   const math::XYZTLorentzVector p4(
       ptr->px(), ptr->py(), pz, sqrt(ptr->px() * ptr->px() + ptr->py() * ptr->py() + pz * pz));
   neutrino_ = new reco::ShallowClonePtrCandidate(ptr, ptr->charge(), p4, ptr->vertex());
 }

References MEzCalculator::Calculate(), reco::LeafCandidate::charge(), Exception, heavyIonCSV_trainingSettings::idx, MEzCalculator::IsComplex(), WDecay::kElec, WDecay::kMuon, TtSemiLepHitFitProducer_Electrons_cfi::leps, leptonType(), BTaggingMonitor_cfi::met, neutrino_, numberOfRealNeutrinoSolutions_, p4, reco::LeafCandidate::px(), reco::LeafCandidate::py(), MEzCalculator::SetLepton(), MEzCalculator::SetMET(), mathSSE::sqrt(), and reco::LeafCandidate::vertex().

Referenced by TtSemiLepHypGenMatch::buildHypo(), and buildHypo().

Member Data Documentation

◆ getMatch_

bool TtSemiLepHypothesis::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 93 of file TtSemiLepHypothesis.h.

Referenced by produce(), and TtSemiLepHypothesis().

◆ hadronicB_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::hadronicB_

protected

Definition at line 114 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), and ~TtSemiLepHypothesis().

◆ jetCorrectionLevel_

std::string TtSemiLepHypothesis::jetCorrectionLevel_

protected

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

Definition at line 102 of file TtSemiLepHypothesis.h.

Referenced by jetCorrectionLevel(), and TtSemiLepHypothesis().

◆ jetsToken_

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

protected

input label for all necessary collections

Definition at line 95 of file TtSemiLepHypothesis.h.

Referenced by produce().

◆ key_

int TtSemiLepHypothesis::key_

protected

hypothesis key (to be set by the buildKey function)

Definition at line 104 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildKey(), TtSemiLepHypWMassDeltaTopMass::buildKey(), TtSemiLepHypWMassMaxSumPt::buildKey(), TtSemiLepHypMVADisc::buildKey(), TtSemiLepHypMaxSumPtWMass::buildKey(), TtSemiLepHypGeom::buildKey(), TtSemiLepHypHitFit::buildKey(), TtSemiLepHypKinFit::buildKey(), and key().

◆ lepsToken_

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

protected

Definition at line 96 of file TtSemiLepHypothesis.h.

Referenced by produce().

◆ lepton_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::lepton_

protected

Definition at line 117 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), and ~TtSemiLepHypothesis().

◆ leptonicB_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::leptonicB_

protected

Definition at line 115 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), and ~TtSemiLepHypothesis().

◆ lightQ_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::lightQ_

protected

candidates for internal use for the creation of the hypothesis candidate

Definition at line 112 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), and ~TtSemiLepHypothesis().

◆ lightQBar_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::lightQBar_

protected

Definition at line 113 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), and ~TtSemiLepHypothesis().

◆ matchToken_

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

protected

Definition at line 98 of file TtSemiLepHypothesis.h.

Referenced by produce(), and TtSemiLepHypothesis().

◆ metsToken_

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

protected

Definition at line 97 of file TtSemiLepHypothesis.h.

Referenced by produce().

◆ neutrino_

reco::ShallowClonePtrCandidate* TtSemiLepHypothesis::neutrino_

protected

Definition at line 116 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), TtSemiLepHypHitFit::buildHypo(), TtSemiLepHypKinFit::buildHypo(), buildHypo(), hypo(), resetCandidates(), setNeutrino(), and ~TtSemiLepHypothesis().

◆ neutrinoSolutionType_

int TtSemiLepHypothesis::neutrinoSolutionType_

protected

algorithm used to calculate neutrino solutions (see cfi for further details)

Definition at line 106 of file TtSemiLepHypothesis.h.

Referenced by TtSemiLepHypGenMatch::buildHypo(), buildHypo(), and TtSemiLepHypothesis().

◆ nJetsConsideredToken_

edm::EDGetTokenT<int> TtSemiLepHypothesis::nJetsConsideredToken_

protected

Definition at line 99 of file TtSemiLepHypothesis.h.

Referenced by produce().

◆ numberOfRealNeutrinoSolutions_

int TtSemiLepHypothesis::numberOfRealNeutrinoSolutions_

protected

number of real neutrino solutions: -1 if not determined, 0 if only complex, 2 if two real solutions

Definition at line 109 of file TtSemiLepHypothesis.h.

Referenced by produce(), resetCandidates(), and setNeutrino().

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TtSemiLepHypothesis()

◆ ~TtSemiLepHypothesis()

Member Function Documentation

◆ buildHypo() [1/2]

◆ buildHypo() [2/2]

◆ buildKey()

◆ hypo()

◆ isValid()

◆ jetCorrectionLevel()

◆ key()

◆ leptonType()

◆ produce()

◆ resetCandidates()

◆ setCandidate() [1/2]

◆ setCandidate() [2/2]

◆ setNeutrino()

Member Data Documentation

◆ getMatch_

◆ hadronicB_

◆ jetCorrectionLevel_

◆ jetsToken_

◆ key_

◆ lepsToken_

◆ lepton_

◆ leptonicB_

◆ lightQ_

◆ lightQBar_

◆ matchToken_

◆ metsToken_

◆ neutrino_

◆ neutrinoSolutionType_

◆ nJetsConsideredToken_

◆ numberOfRealNeutrinoSolutions_