#include <TtFullHadHypothesis.h>

Inheritance diagram for TtFullHadHypothesis:

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

	~TtFullHadHypothesis () override
	default destructor More...

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

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

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)

virtual	~ProducerBase () noexcept(false)

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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

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

virtual	~EDConsumerBase () noexcept(false)

Protected Member Functions
virtual void	buildHypo (edm::Event &event, 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 full-hadronic 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
	return key 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...

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

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

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 ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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

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

Protected Attributes
reco::ShallowClonePtrCandidate *	b_

reco::ShallowClonePtrCandidate *	bBar_

bool	getMatch_

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...

reco::ShallowClonePtrCandidate *	lightP_

reco::ShallowClonePtrCandidate *	lightPBar_

reco::ShallowClonePtrCandidate *	lightQ_

reco::ShallowClonePtrCandidate *	lightQBar_

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

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)

Detailed Description

Definition at line 29 of file TtFullHadHypothesis.h.

Constructor & Destructor Documentation

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

explicit

default constructor

Definition at line 5 of file TtFullHadHypothesis.cc.

References edm::ParameterSet::exists(), getMatch_, edm::ParameterSet::getParameter(), jetCorrectionLevel_, matchToken_, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                   :
   jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
   lightQ_(nullptr), lightQBar_(nullptr), b_(nullptr), bBar_(nullptr), lightP_(nullptr), lightPBar_(nullptr)
 {
   getMatch_ = false;
   if( cfg.exists("match") ) {
     getMatch_ = true;
     matchToken_ = consumes<std::vector<std::vector<int> > >(cfg.getParameter<edm::InputTag>("match"));
   }
   if( cfg.exists("jetCorrectionLevel") ) {
     jetCorrectionLevel_ = cfg.getParameter<std::string>("jetCorrectionLevel");
   }
   produces<std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > >();
   produces<int>("Key");
 }

TtFullHadHypothesis::~TtFullHadHypothesis ( )

override

default destructor

Definition at line 22 of file TtFullHadHypothesis.cc.

References b_, bBar_, lightP_, lightPBar_, lightQ_, and lightQBar_.

 {
   if( lightQ_    ) delete lightQ_;
   if( lightQBar_ ) delete lightQBar_;
   if( b_         ) delete b_;
   if( bBar_      ) delete bBar_;
   if( lightP_    ) delete lightP_;
   if( lightPBar_ ) delete lightPBar_;
 }

Member Function Documentation

virtual void TtFullHadHypothesis::buildHypo	(	edm::Event &	event,
		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 full-hadronic event

Implemented in TtFullHadHypKinFit, and TtFullHadHypGenMatch.

Referenced by isValid(), and produce().

virtual void TtFullHadHypothesis::buildKey ( )

protectedpure virtual

build the event hypothesis key

Implemented in TtFullHadHypKinFit, and TtFullHadHypGenMatch.

Referenced by isValid(), and produce().

reco::CompositeCandidate TtFullHadHypothesis::hypo ( )

protected

return event hypothesis

Definition at line 90 of file TtFullHadHypothesis.cc.

References reco::CompositeCandidate::addDaughter(), TtFullHadDaughter::B, b_, TtFullHadDaughter::BBar, bBar_, TtFullHadDaughter::LightP, lightP_, TtFullHadDaughter::LightPBar, lightPBar_, TtFullHadDaughter::LightQ, lightQ_, TtFullHadDaughter::LightQBar, lightQBar_, AddFourMomenta::set(), AlCaHLTBitMon_QueryRunRegistry::string, TtFullHadDaughter::Top, TtFullHadDaughter::TopBar, w, TtFullHadDaughter::WMinus, and TtFullHadDaughter::WPlus.

Referenced by key(), and produce().

 {
   // check for sanity of the hypothesis
   if( !lightQ_ || !lightQBar_ || !b_ ||
       !bBar_ || !lightP_ || !lightPBar_ )
     return reco::CompositeCandidate();
 
   // setup transient references
   reco::CompositeCandidate hyp, top, w, topBar, wBar;
 
   AddFourMomenta addFourMomenta;
   // build up the top bar branch
   wBar  .addDaughter(*lightP_,    TtFullHadDaughter::LightP    );
   wBar  .addDaughter(*lightPBar_, TtFullHadDaughter::LightPBar );
   addFourMomenta.set( wBar );
   topBar.addDaughter( wBar,  TtFullHadDaughter::WMinus );
   topBar.addDaughter(*bBar_, TtFullHadDaughter::BBar   );
   addFourMomenta.set( topBar );
 
   // build up the top branch that decays hadronically
   w  .addDaughter(*lightQ_,    TtFullHadDaughter::LightQ    );
   w  .addDaughter(*lightQBar_, TtFullHadDaughter::LightQBar );
   addFourMomenta.set( w );
   top.addDaughter( w,  TtFullHadDaughter::WPlus );
   top.addDaughter(*b_, TtFullHadDaughter::B     );
   addFourMomenta.set( top );
 
   // build ttbar hypotheses
   hyp.addDaughter( topBar, TtFullHadDaughter::TopBar );
   hyp.addDaughter( top,    TtFullHadDaughter::Top    );
   addFourMomenta.set( hyp );
 
   return hyp;
 }

bool TtFullHadHypothesis::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 55 of file TtFullHadHypothesis.h.

References buildHypo(), buildKey(), and fwrapper::jets.

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

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

training_settings.idx

idx

Definition: training_settings.py:16

std::string TtFullHadHypothesis::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 127 of file TtFullHadHypothesis.cc.

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

Referenced by TtFullHadHypGenMatch::buildHypo().

 {
   // 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;
 }

int TtFullHadHypothesis::key ( ) const

inlineprotected

return key

Definition at line 51 of file TtFullHadHypothesis.h.

References hypo(), and key_.

Referenced by produce().

51 { return key_; };

TtFullHadHypothesis::key_

int key_

hypothesis key (to be set by the buildKey function)

Definition: TtFullHadHypothesis.h:81

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

overrideprotected

produce the event hypothesis as CompositeCandidate and Key

Definition at line 34 of file TtFullHadHypothesis.cc.

References buildHypo(), buildKey(), edm::Event::getByToken(), getMatch_, hypo(), mps_fire::i, fwrapper::jets, jetsToken_, key(), match(), matchToken_, eostools::move(), edm::Event::put(), and resetCandidates().

 {
   edm::Handle<std::vector<pat::Jet> > jets;
   evt.getByToken(jetsToken_, jets);
 
   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);
 
   // 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, 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");
 }

void TtFullHadHypothesis::resetCandidates ( )

protected

reset candidate pointers before hypo build process

Definition at line 78 of file TtFullHadHypothesis.cc.

References b_, bBar_, lightP_, lightPBar_, lightQ_, and lightQBar_.

Referenced by produce().

 {
   lightQ_    = nullptr;
   lightQBar_ = nullptr;
   b_         = nullptr;
   bBar_      = nullptr;
   lightP_    = nullptr;
   lightPBar_ = nullptr;
 }

template<typename C >

void TtFullHadHypothesis::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 96 of file TtFullHadHypothesis.h.

References cmsBatch::handle, and training_settings::idx.

Referenced by TtFullHadHypGenMatch::buildHypo(), and TtFullHadHypKinFit::buildHypo().

                                                                                                                   {
   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() );
 }

void TtFullHadHypothesis::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 161 of file TtFullHadHypothesis.cc.

References reco::LeafCandidate::charge(), cmsBatch::handle, training_settings::idx, pat::Jet::jecFactor(), reco::LeafCandidate::p4(), AlCaHLTBitMon_QueryRunRegistry::string, and reco::LeafCandidate::vertex().

 {
   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() );
 }