Template class to fill the TtEvent structure. More...

#include "TopQuarkAnalysis/TopEventProducers/interface/TtEvtBuilder.h"

Inheritance diagram for TtEvtBuilder< C >:

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

	~TtEvtBuilder () 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

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)

	~ProducerBase () 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 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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

virtual	~EDConsumerBase () noexcept(false)

Private Types
typedef std::vector< edm::EDGetTokenT< int > >::const_iterator	EventHypoIntToken

typedef std::vector< edm::EDGetTokenT< std::vector< TtEvent::HypoCombPair > > >::const_iterator	EventHypoToken

Private Member Functions
virtual void	fillSpecific (C &, const edm::Event &)
	fill data members that are decay-channel specific More...

template<>
void	fillSpecific (TtFullHadronicEvent &ttEvent, const edm::Event &evt)

template<>
void	fillSpecific (TtFullLeptonicEvent &ttEvent, const edm::Event &evt)

template<>
void	fillSpecific (TtSemiLeptonicEvent &ttEvent, const edm::Event &evt)

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
int	decayChnTop1_

int	decayChnTop2_

edm::EDGetTokenT< std::vector< double > >	discToken_

edm::EDGetTokenT< std::vector< double > >	fitChi2Token_

edm::EDGetTokenT< std::vector< double > >	fitProbToken_

edm::InputTag	genEvt_
	TtGenEvent. More...

edm::EDGetTokenT< TtGenEvent >	genEvtToken_

edm::ParameterSet	genMatch_

edm::ParameterSet	hitFit_

edm::EDGetTokenT< std::vector< double > >	hitFitChi2Token_

edm::EDGetTokenT< std::vector< double > >	hitFitMTToken_

edm::EDGetTokenT< std::vector< double > >	hitFitProbToken_

edm::EDGetTokenT< std::vector< double > >	hitFitSigMTToken_

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

std::vector< edm::EDGetTokenT< int > >	hypKeyTokens_
	vector of hypothesis class names More...

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

std::vector< edm::EDGetTokenT< std::vector< TtEvent::HypoCombPair > > >	hypTokens_

edm::ParameterSet	kinFit_

edm::ParameterSet	kinSolution_

edm::EDGetTokenT< std::string >	methToken_

edm::ParameterSet	mvaDisc_

edm::EDGetTokenT< std::vector< double > >	solWeightToken_

edm::EDGetTokenT< std::vector< double > >	sumDRToken_

edm::EDGetTokenT< std::vector< double > >	sumPtToken_

int	verbosity_
	vebosity level More...

edm::EDGetTokenT< bool >	wrongChargeToken_

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	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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)

Detailed Description

template<typename C>
class TtEvtBuilder< C >

Template class to fill the TtEvent structure.

Template class to fill the TtEvent structure for:

TtSemiLeptonicEvent TtFullLeptonicEvent TtFullHadronicEvent

event hypothesis, genEvent and extra information (if available) are read from the event and contracted into the TtEvent

Definition at line 30 of file TtEvtBuilder.h.

Member Typedef Documentation

template<typename C >

typedef std::vector<edm::EDGetTokenT<int> >::const_iterator TtEvtBuilder< C >::EventHypoIntToken

private

Definition at line 56 of file TtEvtBuilder.h.

template<typename C >

typedef std::vector<edm::EDGetTokenT<std::vector<TtEvent::HypoCombPair> > >::const_iterator TtEvtBuilder< C >::EventHypoToken

private

Definition at line 57 of file TtEvtBuilder.h.

Constructor & Destructor Documentation

template<typename C >

TtEvtBuilder< C >::TtEvtBuilder ( const edm::ParameterSet & cfg )

explicit

default constructor

Definition at line 96 of file TtEvtBuilder.h.

                                                         :
   verbosity_   (cfg.getParameter<int>                        ("verbosity"    )),
   hypKeyTokens_        (edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("hypotheses"   ), [this](edm::InputTag const & tag){return consumes<int>(edm::InputTag(tag.label(), "Key"));})),
   hypTokens_        (edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("hypotheses"   ), [this](edm::InputTag const & tag){return consumes<std::vector<TtEvent::HypoCombPair> >(tag);})),
   hypNeutrTokens_        (edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("hypotheses"   ), [this](edm::InputTag const & tag){return consumes<int>(edm::InputTag(tag.label(), "NumberOfRealNeutrinoSolutions"));})),
   hypJetTokens_        (edm::vector_transform(cfg.getParameter<std::vector<edm::InputTag> >("hypotheses"   ), [this](edm::InputTag const & tag){return consumes<int>(edm::InputTag(tag.label(), "NumberOfConsideredJets"));})),
   genEvt_      (cfg.getParameter<edm::InputTag>              ("genEvent"     )),
   genEvtToken_      (mayConsume<TtGenEvent>(genEvt_)),
   decayChnTop1_(cfg.getParameter<int>                        ("decayChannel1")),
   decayChnTop2_(cfg.getParameter<int>                        ("decayChannel2"))
 {
   // parameter subsets for kKinFit
   if( cfg.exists("kinFit") ) {
     kinFit_  = cfg.getParameter<edm::ParameterSet>("kinFit");
     fitChi2Token_ = mayConsume<std::vector<double> >(kinFit_.getParameter<edm::InputTag>("chi2"));
     fitProbToken_ = mayConsume<std::vector<double> >(kinFit_.getParameter<edm::InputTag>("prob"));
   }
   // parameter subsets for kHitFit
   if( cfg.exists("hitFit") ) {
     hitFit_  = cfg.getParameter<edm::ParameterSet>("hitFit");
     hitFitChi2Token_ = mayConsume<std::vector<double> >(hitFit_.getParameter<edm::InputTag>("chi2"));
     hitFitProbToken_ = mayConsume<std::vector<double> >(hitFit_.getParameter<edm::InputTag>("prob"));
     hitFitMTToken_ = mayConsume<std::vector<double> >(hitFit_.getParameter<edm::InputTag>("mt"));
     hitFitSigMTToken_ = mayConsume<std::vector<double> >(hitFit_.getParameter<edm::InputTag>("sigmt"));
   }
   // parameter subsets for kKinSolution
   if( cfg.exists("kinSolution") ) {
     kinSolution_  = cfg.getParameter<edm::ParameterSet>("kinSolution");
     solWeightToken_    = mayConsume<std::vector<double> >(kinSolution_.getParameter<edm::InputTag>("solWeight"));
     wrongChargeToken_  = mayConsume<bool>(kinSolution_.getParameter<edm::InputTag>("wrongCharge"));
   }
   // parameter subsets for kGenMatch
   if( cfg.exists("genMatch") ) {
     genMatch_ = cfg.getParameter<edm::ParameterSet>("genMatch");
     sumPtToken_    = mayConsume<std::vector<double> >(genMatch_.getParameter<edm::InputTag>("sumPt"));
     sumDRToken_    = mayConsume<std::vector<double> >(genMatch_.getParameter<edm::InputTag>("sumDR"));
   }
   // parameter subsets for kMvaDisc
   if( cfg.exists("mvaDisc") ) {
     mvaDisc_ = cfg.getParameter<edm::ParameterSet>("mvaDisc");
     methToken_    = mayConsume<std::string>(mvaDisc_.getParameter<edm::InputTag>("meth"));
     discToken_    = mayConsume<std::vector<double> >(mvaDisc_.getParameter<edm::InputTag>("disc"));
   }
   // produces a TtEventEvent for:
   //  * TtSemiLeptonicEvent
   //  * TtFullLeptonicEvent
   //  * TtFullHadronicEvent
   // from hypotheses and associated extra information
   produces<C>();
 }

template<typename C >

TtEvtBuilder< C >::~TtEvtBuilder ( )

inlineoverride

default destructor

Definition at line 37 of file TtEvtBuilder.h.

References patCaloMETCorrections_cff::C, TtEvtBuilder< C >::fillSpecific(), and TtEvtBuilder< C >::produce().

37 {};

Member Function Documentation

template<typename C >

virtual void TtEvtBuilder< C >::fillSpecific	(	C &	,
		const edm::Event &
	)

privatevirtual

fill data members that are decay-channel specific

Referenced by TtEvtBuilder< C >::produce(), and TtEvtBuilder< C >::~TtEvtBuilder().

template<>

void TtEvtBuilder< TtFullHadronicEvent >::fillSpecific	(	TtFullHadronicEvent &	ttEvent,
		const edm::Event &	evt
	)

private

Definition at line 244 of file TtEvtBuilder.h.

245 {

246 }

template<>

void TtEvtBuilder< TtFullLeptonicEvent >::fillSpecific	(	TtFullLeptonicEvent &	ttEvent,
		const edm::Event &	evt
	)

private

Definition at line 249 of file TtEvtBuilder.h.

References edm::Event::getByToken(), TtEvent::isHypoAvailable(), TtEvent::kKinSolution, TtFullLeptonicEvent::setSolWeight(), TtFullLeptonicEvent::setWrongCharge(), TtFullLepEvtBuilder_cfi::solWeight, TtEvtBuilder< C >::solWeightToken_, TtFullLepEvtBuilder_cfi::wrongCharge, and TtEvtBuilder< C >::wrongChargeToken_.

 {
 
   // set kKinSolution extras
   if( ttEvent.isHypoAvailable(TtEvent::kKinSolution) ) {
     edm::Handle<std::vector<double> > solWeight;
     evt.getByToken(solWeightToken_, solWeight);
     ttEvent.setSolWeight( *solWeight );
 
     edm::Handle<bool> wrongCharge;
     evt.getByToken(wrongChargeToken_, wrongCharge);
     ttEvent.setWrongCharge( *wrongCharge );
   }
 
 }

template<>

void TtEvtBuilder< TtSemiLeptonicEvent >::fillSpecific	(	TtSemiLeptonicEvent &	ttEvent,
		const edm::Event &	evt
	)

private

Definition at line 266 of file TtEvtBuilder.h.

References edm::Event::getByToken(), TtEvtBuilder< C >::hypJetTokens_, TtEvtBuilder< C >::hypKeyTokens_, TtEvtBuilder< C >::hypNeutrTokens_, crabWrapper::key, TtEvent::setNumberOfConsideredJets(), and TtSemiLeptonicEvent::setNumberOfRealNeutrinoSolutions().

 {
 
   EventHypoIntToken hKey = hypKeyTokens_.begin();
   EventHypoIntToken hNeutr = hypNeutrTokens_.begin();
   EventHypoIntToken hJet = hypJetTokens_.begin();
   for( ; hKey != hypKeyTokens_.end(); ++hKey, ++hNeutr, ++hJet){
     edm::Handle<int> key;
     evt.getByToken(*hKey, key);
 
     // set number of real neutrino solutions for all hypotheses
     edm::Handle<int> numberOfRealNeutrinoSolutions;
     evt.getByToken(*hNeutr, numberOfRealNeutrinoSolutions);
     ttEvent.setNumberOfRealNeutrinoSolutions((TtEvent::HypoClassKey&)*key, *numberOfRealNeutrinoSolutions);
 
     // set number of considered jets for all hypotheses
     edm::Handle<int> numberOfConsideredJets;
     evt.getByToken(*hJet, numberOfConsideredJets);
     ttEvent.setNumberOfConsideredJets((TtEvent::HypoClassKey&)*key, *numberOfConsideredJets);
   }
 
 }

template<typename C >

void TtEvtBuilder< C >::produce	(	edm::Event &	evt,
		const edm::EventSetup &	setup
	)

overrideprivate

produce function (this one is not even accessible for derived classes)

Definition at line 149 of file TtEvtBuilder.h.

Referenced by TtEvtBuilder< C >::~TtEvtBuilder().

 {
   C ttEvent;
 
   // set leptonic decay channels
   ttEvent.setLepDecays( WDecay::LeptonType(decayChnTop1_), WDecay::LeptonType(decayChnTop2_) );
 
   // set genEvent (if available)
   edm::Handle<TtGenEvent> genEvt;
   if ( !genEvt_.label().empty() )
     if( evt.getByToken(genEvtToken_, genEvt) )
       ttEvent.setGenEvent(genEvt);
 
   // add event hypotheses for all given
   // hypothesis classes to the TtEvent
   EventHypoIntToken hKey = hypKeyTokens_.begin();
   EventHypoToken h = hypTokens_.begin();
   for( ; hKey != hypKeyTokens_.end(); ++hKey, ++h){
     edm::Handle<int> key;
     evt.getByToken(*hKey, key);
 
     edm::Handle<std::vector<TtEvent::HypoCombPair> > hypMatchVec;
     evt.getByToken(*h, hypMatchVec);
 
     typedef std::vector<TtEvent::HypoCombPair>::const_iterator HypMatch;
     for(HypMatch hm=hypMatchVec->begin(); hm != hypMatchVec->end(); ++hm){
       ttEvent.addEventHypo((TtEvent::HypoClassKey&)*key, *hm);
     }
   }
 
   // set kKinFit extras
   if( ttEvent.isHypoAvailable(TtEvent::kKinFit) ) {
     edm::Handle<std::vector<double> > fitChi2;
     evt.getByToken(fitChi2Token_, fitChi2);
     ttEvent.setFitChi2( *fitChi2 );
 
     edm::Handle<std::vector<double> > fitProb;
     evt.getByToken(fitProbToken_, fitProb);
     ttEvent.setFitProb( *fitProb );
   }
 
   // set kHitFit extras
   if( ttEvent.isHypoAvailable(TtEvent::kHitFit) ) {
     edm::Handle<std::vector<double> > hitFitChi2;
     evt.getByToken(hitFitChi2Token_, hitFitChi2);
     ttEvent.setHitFitChi2( *hitFitChi2 );
 
     edm::Handle<std::vector<double> > hitFitProb;
     evt.getByToken(hitFitProbToken_, hitFitProb);
     ttEvent.setHitFitProb( *hitFitProb );
 
     edm::Handle<std::vector<double> > hitFitMT;
     evt.getByToken(hitFitMTToken_, hitFitMT);
     ttEvent.setHitFitMT( *hitFitMT );
 
     edm::Handle<std::vector<double> > hitFitSigMT;
     evt.getByToken(hitFitSigMTToken_, hitFitSigMT);
     ttEvent.setHitFitSigMT( *hitFitSigMT );
   }
 
   // set kGenMatch extras
   if( ttEvent.isHypoAvailable(TtEvent::kGenMatch) ) {
     edm::Handle<std::vector<double> > sumPt;
     evt.getByToken(sumPtToken_, sumPt);
     ttEvent.setGenMatchSumPt( *sumPt );
 
     edm::Handle<std::vector<double> > sumDR;
     evt.getByToken(sumDRToken_, sumDR);
     ttEvent.setGenMatchSumDR( *sumDR );
   }
 
   // set kMvaDisc extras
   if( ttEvent.isHypoAvailable(TtEvent::kMVADisc) ) {
     edm::Handle<std::string> meth;
     evt.getByToken(methToken_, meth);
     ttEvent.setMvaMethod( *meth );
 
     edm::Handle<std::vector<double> > disc;
     evt.getByToken(discToken_, disc);
     ttEvent.setMvaDiscriminators( *disc );
   }
 
   // fill data members that are decay-channel specific
   fillSpecific(ttEvent, evt);
 
   // print summary via MessageLogger if verbosity_>0
   ttEvent.print(verbosity_);
 
   // write object into the edm::Event
   std::unique_ptr<C> pOut(new C);
   *pOut=ttEvent;
   evt.put(std::move(pOut));
 }