#include <TtJetPartonMatch.h>

Inheritance diagram for TtJetPartonMatch< C >:

Public Member Functions
virtual void	produce (edm::Event &, const edm::EventSetup &)
	write jet parton match objects into the event More...

	TtJetPartonMatch (const edm::ParameterSet &)
	default conructor More...

	~TtJetPartonMatch ()
	default destructor More...

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

ModuleDescription const &	moduleDescription () const

virtual	~EDProducer ()

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

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
JetPartonMatching::algorithms	readAlgorithm (const std::string &str)
	convert string for algorithm into corresponding enumerator type More...

Private Attributes
JetPartonMatching::algorithms	algorithm_
	choice of algorithm More...

edm::EDGetTokenT< TtGenEvent >	genEvt_
	TtGenEvent collection input. More...

edm::EDGetTokenT< edm::View < reco::Jet > >	jets_
	jet collection input More...

double	maxDist_

int	maxNComb_

int	maxNJets_

C	partons_
	partons More...

bool	useDeltaR_
	switch to choose between deltaR/deltaTheta matching More...

bool	useMaxDist_

int	verbosity_
	verbosity level More...

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

Public Types inherited from edm::ProducerBase
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)

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 TtJetPartonMatch< C >

Definition at line 44 of file TtJetPartonMatch.h.

Constructor & Destructor Documentation

template<typename C >

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

explicit

default conructor

Definition at line 87 of file TtJetPartonMatch.h.

                                                                :
   partons_   (cfg.getParameter<std::vector<std::string> >("partonsToIgnore")),
   genEvt_    (consumes<TtGenEvent>(edm::InputTag("genEvt"))),
   jets_      (consumes<edm::View<reco::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
   maxNJets_  (cfg.getParameter<int>                      ("maxNJets"       )),
   maxNComb_  (cfg.getParameter<int>                      ("maxNComb"       )),
   algorithm_ (readAlgorithm(cfg.getParameter<std::string>("algorithm"      ))),
   useDeltaR_ (cfg.getParameter<bool>                     ("useDeltaR"      )),
   useMaxDist_(cfg.getParameter<bool>                     ("useMaxDist"     )),
   maxDist_   (cfg.getParameter<double>                   ("maxDist"        )),
   verbosity_ (cfg.getParameter<int>                      ("verbosity"      ))
 {
   // produces a vector of jet/lepton indices in the order of
   //  * TtSemiLepEvtPartons
   //  * TtFullHadEvtPartons
   //  * TtFullLepEvtPartons
   // and vectors of the corresponding quality parameters
   produces<std::vector<std::vector<int> > >();
   produces<std::vector<double> >("SumPt");
   produces<std::vector<double> >("SumDR");
   produces<int>("NumberOfConsideredJets");
 }

template<typename C >

TtJetPartonMatch< C >::~TtJetPartonMatch ( )

default destructor

Definition at line 111 of file TtJetPartonMatch.h.

112 {

113 }

Member Function Documentation

template<typename C >

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

virtual

write jet parton match objects into the event

Implements edm::EDProducer.

Definition at line 117 of file TtJetPartonMatch.h.

References TtGenEvtProducer_cfi::genEvt, genEvt_(), edm::Event::getByToken(), JetPartonMatching::getMatchesForPartons(), JetPartonMatching::getNumberOfAvailableCombinations(), JetPartonMatching::getSumDeltaPt(), JetPartonMatching::getSumDeltaR(), fwrapper::jets, match(), matches, JetPartonMatching::print(), and edm::Event::put().

 {
   // will write
   // * parton match
   // * sumPt
   // * sumDR
   // to the event
   std::auto_ptr<std::vector<std::vector<int> > > match(new std::vector<std::vector<int> >);
   std::auto_ptr<std::vector<double> > sumPt(new std::vector<double>);
   std::auto_ptr<std::vector<double> > sumDR(new std::vector<double>);
   std::auto_ptr<int> pJetsConsidered(new int);
 
   // get TtGenEvent and jet collection from the event
   edm::Handle<TtGenEvent> genEvt;
   evt.getByToken(genEvt_, genEvt);
 
   edm::Handle<edm::View<reco::Jet> > topJets;
   evt.getByToken(jets_, topJets);
 
   // fill vector of partons in the order of
   //  * TtFullLepEvtPartons
   //  * TtSemiLepEvtPartons
   //  * TtFullHadEvtPartons
   std::vector<const reco::Candidate*> partons = partons_.vec(*genEvt);
 
   // prepare vector of jets
   std::vector<const reco::Candidate*> jets;
   for(unsigned int ij=0; ij<topJets->size(); ++ij) {
     // take all jets if maxNJets_ == -1; otherwise use
     // maxNJets_ if maxNJets_ is big enough or use same
     // number of jets as partons if maxNJets_ < number
     // of partons
     if(maxNJets_!=-1) {
       if(maxNJets_>=(int)partons.size()) {
     if((int)ij==maxNJets_) break;
       }
       else {
     if(ij==partons.size()) break;
       }
     }
     jets.push_back( (const reco::Candidate*) &(*topJets)[ij] );
   }
   *pJetsConsidered = jets.size();
 
   // do the matching with specified parameters
   JetPartonMatching jetPartonMatch(partons, jets, algorithm_, useMaxDist_, useDeltaR_, maxDist_);
 
   // print some info for each event
   // if corresponding verbosity level set
   if(verbosity_>0)
     jetPartonMatch.print();
 
   for(unsigned int ic=0; ic<jetPartonMatch.getNumberOfAvailableCombinations(); ++ic) {
     if((int)ic>=maxNComb_ && maxNComb_>=0) break;
     std::vector<int> matches = jetPartonMatch.getMatchesForPartons(ic);
     partons_.expand(matches); // insert dummy indices for partons that were chosen to be ignored
     match->push_back( matches );
     sumPt->push_back( jetPartonMatch.getSumDeltaPt(ic) );
     sumDR->push_back( jetPartonMatch.getSumDeltaR (ic) );
   }
   evt.put(match);
   evt.put(sumPt, "SumPt");
   evt.put(sumDR, "SumDR");
   evt.put(pJetsConsidered, "NumberOfConsideredJets");
 }

template<typename C >

JetPartonMatching::algorithms TtJetPartonMatch< C >::readAlgorithm ( const std::string & str )

private

convert string for algorithm into corresponding enumerator type

Definition at line 185 of file TtJetPartonMatch.h.

References Exception, JetPartonMatching::minSumDist, JetPartonMatching::ptOrderedMinDist, JetPartonMatching::totalMinDist, and JetPartonMatching::unambiguousOnly.

 {
   if     (str == "totalMinDist"    ) return JetPartonMatching::totalMinDist;
   else if(str == "minSumDist"      ) return JetPartonMatching::minSumDist;
   else if(str == "ptOrderedMinDist") return JetPartonMatching::ptOrderedMinDist;
   else if(str == "unambiguousOnly" ) return JetPartonMatching::unambiguousOnly;
   else throw cms::Exception("Configuration")
     << "Chosen algorithm is not supported: " << str << "\n";
 }

Member Data Documentation

template<typename C >

JetPartonMatching::algorithms TtJetPartonMatch< C >::algorithm_

private

choice of algorithm

Definition at line 73 of file TtJetPartonMatch.h.

template<typename C >

edm::EDGetTokenT<TtGenEvent> TtJetPartonMatch< C >::genEvt_

private

TtGenEvent collection input.

Definition at line 63 of file TtJetPartonMatch.h.

template<typename C >

edm::EDGetTokenT<edm::View<reco::Jet> > TtJetPartonMatch< C >::jets_

private

jet collection input

Definition at line 65 of file TtJetPartonMatch.h.

template<typename C >

double TtJetPartonMatch< C >::maxDist_

private

threshold for outliers in the case that useMaxDist_ =true

Definition at line 81 of file TtJetPartonMatch.h.

template<typename C >

int TtJetPartonMatch< C >::maxNComb_

private

maximal number of combinations for which the matching should be stored

Definition at line 71 of file TtJetPartonMatch.h.

template<typename C >

int TtJetPartonMatch< C >::maxNJets_

private

maximal number of jets to be considered for the matching

Definition at line 68 of file TtJetPartonMatch.h.

template<typename C >

C TtJetPartonMatch< C >::partons_

private

partons

Definition at line 61 of file TtJetPartonMatch.h.

template<typename C >

bool TtJetPartonMatch< C >::useDeltaR_

private

switch to choose between deltaR/deltaTheta matching

Definition at line 75 of file TtJetPartonMatch.h.

template<typename C >

bool TtJetPartonMatch< C >::useMaxDist_

private

switch to choose whether an outlier rejection should be applied or not

Definition at line 78 of file TtJetPartonMatch.h.

template<typename C >

int TtJetPartonMatch< C >::verbosity_

private

verbosity level

Definition at line 83 of file TtJetPartonMatch.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

template<typename C> class TtJetPartonMatch< C >

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation

template<typename C>
class TtJetPartonMatch< C >