#include <TtFullLepHypothesis.h>

Inheritance diagram for TtFullLepHypothesis:

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

	~TtFullLepHypothesis ()
	default destructor More...

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

virtual	~EDProducer ()

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

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

virtual	~ProducerBase ()

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

virtual void	produce (edm::Event &, const edm::EventSetup &)
	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::EDProducer
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Protected Attributes
reco::ShallowClonePtrCandidate *	b_

reco::ShallowClonePtrCandidate *	bBar_

edm::InputTag	elecs_

bool	getMatch_

std::string	jetCorrectionLevel_

edm::InputTag	jets_

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

reco::ShallowClonePtrCandidate *	lepton_

reco::ShallowClonePtrCandidate *	leptonBar_

edm::InputTag	match_
	input label for all necessary collections More...

edm::InputTag	mets_

edm::InputTag	mus_

reco::ShallowClonePtrCandidate *	neutrino_

reco::ShallowClonePtrCandidate *	neutrinoBar_

reco::LeafCandidate *	recNu
	candidates needed for the genmatch hypothesis More...

reco::LeafCandidate *	recNuBar

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

typedef WorkerT< EDProducer >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

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

Constructor & Destructor Documentation

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

explicit

default constructor

Definition at line 5 of file TtFullLepHypothesis.cc.

References edm::ParameterSet::exists(), getMatch_, edm::ParameterSet::getParameter(), jetCorrectionLevel_, and match_.

                                                                   :
   elecs_(cfg.getParameter<edm::InputTag>("electrons")),
   mus_  (cfg.getParameter<edm::InputTag>("muons")),
   jets_ (cfg.getParameter<edm::InputTag>("jets")),
   mets_ (cfg.getParameter<edm::InputTag>("mets")),
 
   lepton_(0), leptonBar_(0), b_(0), 
   bBar_(0), neutrino_(0), neutrinoBar_(0)
 {  
   getMatch_ = false;
   if( cfg.exists("match") ) {
     getMatch_ = true;
     match_ = 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");
 }

TtFullLepHypothesis::~TtFullLepHypothesis ( )

default destructor

Definition at line 32 of file TtFullLepHypothesis.cc.

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

 {
   if( lepton_      ) delete lepton_;
   if( leptonBar_   ) delete leptonBar_;
   if( b_           ) delete b_;
   if( bBar_        ) delete bBar_;
   if( neutrino_    ) delete neutrino_;
   if( neutrinoBar_ ) delete neutrinoBar_;
   //if( met_         ) delete met_;
 }

Member Function Documentation

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

virtual void TtFullLepHypothesis::buildKey ( )

protectedpure virtual

build the event hypothesis key

Implemented in TtFullLepHypGenMatch, and TtFullLepHypKinSolution.

Referenced by produce().

reco::CompositeCandidate TtFullLepHypothesis::hypo ( )

protected

return event hypothesis

Definition at line 111 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::LogWarning("TtFullHypothesis") << "no neutrinos for gen match" << std::endl; 
     return reco::CompositeCandidate();
   }  
   if( key()==TtFullLeptonicEvent::kKinSolution && (!neutrino_ || !neutrinoBar_) ){
     edm::LogWarning("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(*leptonBar_, TtFullLepDaughter::LepBar);
   if(key()==TtFullLeptonicEvent::kKinSolution)
     WPlus.addDaughter(*neutrino_, TtFullLepDaughter::Nu);  
   else if(key()==TtFullLeptonicEvent::kGenMatch)
     WPlus.addDaughter(*recNu, TtFullLepDaughter::Nu);     
   addFourMomenta.set(WPlus);
   Top.addDaughter(WPlus, TtFullLepDaughter::WPlus);
   Top.addDaughter(*b_,TtFullLepDaughter::B);
   addFourMomenta.set(Top);   
    
   // build up the anti top branch
   WMinus.addDaughter(*lepton_, TtFullLepDaughter::Lep); 
   if(key()==TtFullLeptonicEvent::kKinSolution)
     WMinus.addDaughter(*neutrinoBar_, TtFullLepDaughter::NuBar);  
   else if(key()==TtFullLeptonicEvent::kGenMatch)  
     WMinus.addDaughter(*recNuBar, TtFullLepDaughter::NuBar);   
   addFourMomenta.set(WMinus);
   TopBar.addDaughter(WMinus, TtFullLepDaughter::WMinus);
   TopBar.addDaughter(*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;
 }

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

References fwrapper::jets.

Referenced by TtFullLepHypGenMatch::buildHypo().

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

int TtFullLepHypothesis::key ( ) const

inlineprotected

return key

Definition at line 53 of file TtFullLepHypothesis.h.

References key_.

Referenced by hypo(), and produce().

53 { return key_; };

TtFullLepHypothesis::key_

int key_

hypothesis key (to be set by the buildKey function)

Definition: TtFullLepHypothesis.h:89

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

protectedvirtual

produce the event hypothesis as CompositeCandidate and Key

Implements edm::EDProducer.

Definition at line 45 of file TtFullLepHypothesis.cc.

References buildHypo(), buildKey(), elecs_, edm::Event::getByLabel(), getMatch_, hypo(), i, fwrapper::jets, jets_, key(), match(), match_, mets_, mus_, edm::Event::put(), and resetCandidates().

 {
   edm::Handle<std::vector<pat::Electron> > elecs;
   evt.getByLabel(elecs_, elecs);
   
   edm::Handle<std::vector<pat::Muon> > mus;
   evt.getByLabel(mus_, mus);  
 
   edm::Handle<std::vector<pat::Jet> > jets;
   evt.getByLabel(jets_, jets);
 
   edm::Handle<std::vector<pat::MET> > mets;
   evt.getByLabel(mets_, mets);
 
   std::vector<std::vector<int> > matchVec;
   if( getMatch_ ) {
     edm::Handle<std::vector<std::vector<int> > > matchHandle;
     evt.getByLabel(match_, 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 auto_ptr for products
   std::auto_ptr<std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > >
     pOut( new std::vector<std::pair<reco::CompositeCandidate, std::vector<int> > > );
   std::auto_ptr<int> pKey(new int);
   
   // build and feed out key
   buildKey();  
   *pKey=key();
   evt.put(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(pOut);
 }

void TtFullLepHypothesis::resetCandidates ( )

protected

reset candidate pointers before hypo build process

Definition at line 98 of file TtFullLepHypothesis.cc.

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

Referenced by produce().

 {
   lepton_     = 0;
   leptonBar_  = 0;
   b_          = 0;
   bBar_       = 0;
   neutrino_   = 0;
   neutrinoBar_= 0;
   //met_        = 0;
 }

template<typename C >

void TtFullLepHypothesis::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 109 of file TtFullLepHypothesis.h.

References patZpeak::handle.

Referenced by TtFullLepHypKinSolution::buildHypo(), and TtFullLepHypGenMatch::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 TtFullLepHypothesis::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 167 of file TtFullLepHypothesis.cc.

References patZpeak::handle, and jetCorrectionLevel_.

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