#include <TtFullHadKinFitProducer.h>

Inheritance diagram for TtFullHadKinFitProducer:

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

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

Public Attributes
TtFullHadKinFitter::KinFit *	kinFitter
	kinematic fit interface More...

Private Member Functions
virtual void	produce (edm::Event &event, const edm::EventSetup &setup)
	produce fitted object collections and meta data describing fit quality More...

Private Attributes
std::vector< edm::ParameterSet >	bResolutions_

std::string	bTagAlgo_
	input tag for b-tagging algorithm More...

unsigned int	bTags_
	minimal number of b-jets More...

std::vector< unsigned >	constraints_
	numbering of different possible kinematic constraints More...

std::string	jetCorrectionLevel_
	correction level for jets More...

unsigned int	jetParam_
	numbering of different possible jet parametrizations More...

edm::InputTag	jets_
	input tag for jets More...

edm::InputTag	match_
	input tag for matches (in case the fit should be performed on certain matches) More...

double	maxBTagValueNonBJet_
	max value of bTag for a non-b-jet More...

double	maxDeltaS_
	maximal chi2 equivalent More...

double	maxF_
	maximal deviation for contstraints More...

int	maxNComb_
	maximal number of combinations to be written to the event More...

int	maxNJets_
	maximal number of jets (-1 possible to indicate 'all') More...

unsigned int	maxNrIter_
	maximal number of iterations to be performed for the fit More...

double	minBTagValueBJet_
	min value of bTag for a b-jet More...

double	mTop_
	top mass value used for constraints More...

double	mW_
	W mass value used for constraints. More...

double	resolutionSmearFactor_
	smearing factor for jet energy resolutions More...

std::vector< edm::ParameterSet >	udscResolutions_
	store the resolutions for the jets More...

bool	useBTagging_
	switch to tell whether to use b-tagging or not More...

bool	useOnlyMatch_

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)

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)

Detailed Description

Definition at line 22 of file TtFullHadKinFitProducer.h.

Constructor & Destructor Documentation

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

explicit

default constructor

Definition at line 6 of file TtFullHadKinFitProducer.cc.

References bResolutions_, bTagAlgo_, bTags_, constraints_, edm::hlt::Exception, edm::ParameterSet::exists(), edm::ParameterSet::getParameter(), jetCorrectionLevel_, jetParam_, kinFitter, maxBTagValueNonBJet_, maxDeltaS_, maxF_, maxNComb_, maxNJets_, maxNrIter_, minBTagValueBJet_, mTop_, mW_, resolutionSmearFactor_, udscResolutions_, and useBTagging_.

                                                                           :
   jets_                 (cfg.getParameter<edm::InputTag>("jets")),
   match_                (cfg.getParameter<edm::InputTag>("match")),
   useOnlyMatch_         (cfg.getParameter<bool>("useOnlyMatch")),
   bTagAlgo_             (cfg.getParameter<std::string>("bTagAlgo")),
   minBTagValueBJet_     (cfg.getParameter<double>("minBTagValueBJet")),
   maxBTagValueNonBJet_  (cfg.getParameter<double>("maxBTagValueNonBJet")),
   useBTagging_          (cfg.getParameter<bool>("useBTagging")),
   bTags_                (cfg.getParameter<unsigned int>("bTags")),
   jetCorrectionLevel_   (cfg.getParameter<std::string>("jetCorrectionLevel")),
   maxNJets_             (cfg.getParameter<int>("maxNJets")),
   maxNComb_             (cfg.getParameter<int>("maxNComb")),
   maxNrIter_            (cfg.getParameter<unsigned int>("maxNrIter")),
   maxDeltaS_            (cfg.getParameter<double>("maxDeltaS")),
   maxF_                 (cfg.getParameter<double>("maxF")),
   jetParam_             (cfg.getParameter<unsigned>("jetParametrisation")),
   constraints_          (cfg.getParameter<std::vector<unsigned> >("constraints")),
   mW_                   (cfg.getParameter<double>("mW"  )),
   mTop_                 (cfg.getParameter<double>("mTop")),
   resolutionSmearFactor_(cfg.getParameter<double>("resolutionSmearFactor"))
 {
   if(cfg.exists("udscResolutions") && cfg.exists("bResolutions")){
     udscResolutions_ = cfg.getParameter <std::vector<edm::ParameterSet> >("udscResolutions");
     bResolutions_    = cfg.getParameter <std::vector<edm::ParameterSet> >("bResolutions");
   }
   else if(cfg.exists("udscResolutions") || cfg.exists("bResolutions")){
     if(cfg.exists("udscResolutions")) throw cms::Exception("WrongConfig") << "Parameter 'bResolutions' is needed if parameter 'udscResolutions' is defined!\n";
     else                              throw cms::Exception("WrongConfig") << "Parameter 'udscResolutions' is needed if parameter 'bResolutions' is defined!\n";
   }
 
   // define kinematic fit interface
   kinFitter = new TtFullHadKinFitter::KinFit(useBTagging_, bTags_, bTagAlgo_, minBTagValueBJet_, maxBTagValueNonBJet_,
                          udscResolutions_, bResolutions_, resolutionSmearFactor_ ,
                          jetCorrectionLevel_, maxNJets_, maxNComb_,
                          maxNrIter_, maxDeltaS_, maxF_, jetParam_, constraints_, mW_, mTop_);
 
   // produces the following collections
   produces< std::vector<pat::Particle> >("PartonsB");
   produces< std::vector<pat::Particle> >("PartonsBBar");
   produces< std::vector<pat::Particle> >("PartonsLightQ");
   produces< std::vector<pat::Particle> >("PartonsLightQBar");
   produces< std::vector<pat::Particle> >("PartonsLightP");
   produces< std::vector<pat::Particle> >("PartonsLightPBar");
 
   produces< std::vector<std::vector<int> > >();
   produces< std::vector<double> >("Chi2");
   produces< std::vector<double> >("Prob");
   produces< std::vector<int> >("Status");
 }

TtFullHadKinFitProducer::~TtFullHadKinFitProducer ( )

default destructor

Definition at line 57 of file TtFullHadKinFitProducer.cc.

References kinFitter.

 {
   delete kinFitter;
 }

Member Function Documentation

void TtFullHadKinFitProducer::produce	(	edm::Event &	event,
		const edm::EventSetup &	setup
	)

privatevirtual

produce fitted object collections and meta data describing fit quality

set match to be used

set the validity of a match

Implements edm::EDProducer.

Definition at line 64 of file TtFullHadKinFitProducer.cc.

References TtFullHadKinFitter::KinFit::fit(), analyzePatCleaning_cfg::jets, jets_, kinFitter, match(), match_, maxNComb_, nPartons, TtFullHadKinFitter::KinFit::setMatch(), TtFullHadKinFitter::KinFit::setMatchInvalidity(), TtFullHadKinFitter::KinFit::setUseOnlyMatch(), and useOnlyMatch_.

 {
   // get jet collection
   edm::Handle<std::vector<pat::Jet> > jets;
   event.getByLabel(jets_, jets);
 
   // get match in case that useOnlyMatch_ is true
   std::vector<int> match;
   bool invalidMatch=false;
   if(useOnlyMatch_) {
     kinFitter->setUseOnlyMatch(true);
     // in case that only a ceratin match should be used, get match here
     edm::Handle<std::vector<std::vector<int> > > matches;
     event.getByLabel(match_, matches);
     match = *(matches->begin());
     // check if match is valid
     if( match.size()!=nPartons ){
       invalidMatch=true;
     }
     else {
       for(unsigned int idx=0; idx<match.size(); ++idx) {
     if(match[idx]<0 || match[idx]>=(int)jets->size()) {
       invalidMatch=true;
       break;
     }
       }
     }
     kinFitter->setMatch(match);
   }
 
   kinFitter->setMatchInvalidity(invalidMatch);
 
   std::list<TtFullHadKinFitter::KinFitResult> fitResults = kinFitter->fit(*jets);
 
   // pointer for output collections
   std::auto_ptr< std::vector<pat::Particle> > pPartonsB( new std::vector<pat::Particle> );
   std::auto_ptr< std::vector<pat::Particle> > pPartonsBBar( new std::vector<pat::Particle> );
   std::auto_ptr< std::vector<pat::Particle> > pPartonsLightQ   ( new std::vector<pat::Particle> );
   std::auto_ptr< std::vector<pat::Particle> > pPartonsLightQBar( new std::vector<pat::Particle> );
   std::auto_ptr< std::vector<pat::Particle> > pPartonsLightP   ( new std::vector<pat::Particle> );
   std::auto_ptr< std::vector<pat::Particle> > pPartonsLightPBar( new std::vector<pat::Particle> );
   // pointer for meta information
   std::auto_ptr< std::vector<std::vector<int> > > pCombi ( new std::vector<std::vector<int> > );
   std::auto_ptr< std::vector<double> > pChi2  ( new std::vector<double> );
   std::auto_ptr< std::vector<double> > pProb  ( new std::vector<double> );
   std::auto_ptr< std::vector<int> > pStatus( new std::vector<int> );
 
   unsigned int iComb = 0;
   for(std::list<TtFullHadKinFitter::KinFitResult>::const_iterator res = fitResults.begin(); res != fitResults.end(); ++res){
     if(maxNComb_>=1 && iComb==(unsigned int)maxNComb_){ 
       break;
     }
     ++iComb;
 
     pPartonsB        ->push_back( res->B         );
     pPartonsBBar     ->push_back( res->BBar      );
     pPartonsLightQ   ->push_back( res->LightQ    );
     pPartonsLightQBar->push_back( res->LightQBar );
     pPartonsLightP   ->push_back( res->LightP    );
     pPartonsLightPBar->push_back( res->LightPBar );
 
     pCombi ->push_back( res->JetCombi );
     pChi2  ->push_back( res->Chi2     );
     pProb  ->push_back( res->Prob     );
     pStatus->push_back( res->Status   );
 
   }
 
   event.put(pCombi);
   event.put(pPartonsB        , "PartonsB"        );
   event.put(pPartonsBBar     , "PartonsBBar"     );
   event.put(pPartonsLightQ   , "PartonsLightQ"   );
   event.put(pPartonsLightQBar, "PartonsLightQBar");
   event.put(pPartonsLightP   , "PartonsLightP"   );
   event.put(pPartonsLightPBar, "PartonsLightPBar");
   event.put(pChi2   , "Chi2"   );
   event.put(pProb   , "Prob"   );
   event.put(pStatus , "Status" );
 }