Inheritance diagram for TtFullHadKinFitProducer:

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

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Public Attributes
std::unique_ptr< TtFullHadKinFitter::KinFit >	kinFitter
	kinematic fit interface More...

Private Member Functions
void	produce (edm::Event &event, const edm::EventSetup &setup) override
	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...

std::vector< double >	jetEnergyResolutionEtaBinning_

std::vector< double >	jetEnergyResolutionScaleFactors_
	scale factors for jet energy resolution More...

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

edm::EDGetTokenT< std::vector< pat::Jet > >	jetsToken_
	input tag for jets More...

edm::EDGetTokenT< std::vector< std::vector< int > > >	matchToken_
	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...

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::stream::EDProducer<>
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 19 of file TtFullHadKinFitProducer.cc.

Constructor & Destructor Documentation

◆ TtFullHadKinFitProducer()

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

explicit

default constructor

Definition at line 80 of file TtFullHadKinFitProducer.cc.

References bResolutions_, bTagAlgo_, bTags_, looper::cfg, constraints_, Exception, jetCorrectionLevel_, jetEnergyResolutionEtaBinning_, jetEnergyResolutionScaleFactors_, jetParam_, kinFitter, maxBTagValueNonBJet_, maxDeltaS_, maxF_, maxNComb_, maxNJets_, maxNrIter_, minBTagValueBJet_, mTop_, mW_, udscResolutions_, and useBTagging_.

     : jetsToken_(consumes<std::vector<pat::Jet> >(cfg.getParameter<edm::InputTag>("jets"))),
       matchToken_(mayConsume<std::vector<std::vector<int> > >(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")),
       jetEnergyResolutionScaleFactors_(cfg.getParameter<std::vector<double> >("jetEnergyResolutionScaleFactors")),
       jetEnergyResolutionEtaBinning_(cfg.getParameter<std::vector<double> >("jetEnergyResolutionEtaBinning")) {
   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("Configuration")
           << "Parameter 'bResolutions' is needed if parameter 'udscResolutions' is defined!\n";
     else
       throw cms::Exception("Configuration")
           << "Parameter 'udscResolutions' is needed if parameter 'bResolutions' is defined!\n";
   }
 
   // define kinematic fit interface
   kinFitter = std::make_unique<TtFullHadKinFitter::KinFit>(useBTagging_,
                                                            bTags_,
                                                            bTagAlgo_,
                                                            minBTagValueBJet_,
                                                            maxBTagValueNonBJet_,
                                                            udscResolutions_,
                                                            bResolutions_,
                                                            jetEnergyResolutionScaleFactors_,
                                                            jetEnergyResolutionEtaBinning_,
                                                            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");
 }

Member Function Documentation

◆ produce()

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

overrideprivate

produce fitted object collections and meta data describing fit quality

set match to be used

set the validity of a match

Definition at line 149 of file TtFullHadKinFitProducer.cc.

References heavyIonCSV_trainingSettings::idx, singleTopDQM_cfi::jets, jetsToken_, kinFitter, match(), oniaPATMuonsWithTrigger_cff::matches, matchToken_, maxNComb_, eostools::move(), nPartons, useOnlyMatch_, and trackerHitRTTI::vector.

                                                                                  {
   // get jet collection
   edm::Handle<std::vector<pat::Jet> > jets;
   event.getByToken(jetsToken_, 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.getByToken(matchToken_, 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::unique_ptr<std::vector<pat::Particle> > pPartonsB(new std::vector<pat::Particle>);
   std::unique_ptr<std::vector<pat::Particle> > pPartonsBBar(new std::vector<pat::Particle>);
   std::unique_ptr<std::vector<pat::Particle> > pPartonsLightQ(new std::vector<pat::Particle>);
   std::unique_ptr<std::vector<pat::Particle> > pPartonsLightQBar(new std::vector<pat::Particle>);
   std::unique_ptr<std::vector<pat::Particle> > pPartonsLightP(new std::vector<pat::Particle>);
   std::unique_ptr<std::vector<pat::Particle> > pPartonsLightPBar(new std::vector<pat::Particle>);
   // pointer for meta information
   std::unique_ptr<std::vector<std::vector<int> > > pCombi(new std::vector<std::vector<int> >);
   std::unique_ptr<std::vector<double> > pChi2(new std::vector<double>);
   std::unique_ptr<std::vector<double> > pProb(new std::vector<double>);
   std::unique_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(std::move(pCombi));
   event.put(std::move(pPartonsB), "PartonsB");
   event.put(std::move(pPartonsBBar), "PartonsBBar");
   event.put(std::move(pPartonsLightQ), "PartonsLightQ");
   event.put(std::move(pPartonsLightQBar), "PartonsLightQBar");
   event.put(std::move(pPartonsLightP), "PartonsLightP");
   event.put(std::move(pPartonsLightPBar), "PartonsLightPBar");
   event.put(std::move(pChi2), "Chi2");
   event.put(std::move(pProb), "Prob");
   event.put(std::move(pStatus), "Status");
 }

Member Data Documentation

◆ bResolutions_

std::vector<edm::ParameterSet> TtFullHadKinFitProducer::bResolutions_

private

Definition at line 67 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ bTagAlgo_

std::string TtFullHadKinFitProducer::bTagAlgo_

private

input tag for b-tagging algorithm

Definition at line 37 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ bTags_

unsigned int TtFullHadKinFitProducer::bTags_

private

minimal number of b-jets

Definition at line 45 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ constraints_

std::vector<unsigned> TtFullHadKinFitProducer::constraints_

private

numbering of different possible kinematic constraints

Definition at line 61 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ jetCorrectionLevel_

std::string TtFullHadKinFitProducer::jetCorrectionLevel_

private

correction level for jets

Definition at line 47 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ jetEnergyResolutionEtaBinning_

std::vector<double> TtFullHadKinFitProducer::jetEnergyResolutionEtaBinning_

private

Definition at line 70 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ jetEnergyResolutionScaleFactors_

std::vector<double> TtFullHadKinFitProducer::jetEnergyResolutionScaleFactors_

private

scale factors for jet energy resolution

Definition at line 69 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ jetParam_

unsigned int TtFullHadKinFitProducer::jetParam_

private

numbering of different possible jet parametrizations

Definition at line 59 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ jetsToken_

edm::EDGetTokenT<std::vector<pat::Jet> > TtFullHadKinFitProducer::jetsToken_

private

input tag for jets

Definition at line 30 of file TtFullHadKinFitProducer.cc.

Referenced by produce().

◆ kinFitter

std::unique_ptr<TtFullHadKinFitter::KinFit> TtFullHadKinFitProducer::kinFitter

kinematic fit interface

Definition at line 74 of file TtFullHadKinFitProducer.cc.

Referenced by produce(), and TtFullHadKinFitProducer().

◆ matchToken_

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

private

input tag for matches (in case the fit should be performed on certain matches)

Definition at line 32 of file TtFullHadKinFitProducer.cc.

Referenced by produce().

◆ maxBTagValueNonBJet_

double TtFullHadKinFitProducer::maxBTagValueNonBJet_

private

max value of bTag for a non-b-jet

Definition at line 41 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ maxDeltaS_

double TtFullHadKinFitProducer::maxDeltaS_

private

maximal chi2 equivalent

Definition at line 55 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ maxF_

double TtFullHadKinFitProducer::maxF_

private

maximal deviation for contstraints

Definition at line 57 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ maxNComb_

int TtFullHadKinFitProducer::maxNComb_

private

maximal number of combinations to be written to the event

Definition at line 51 of file TtFullHadKinFitProducer.cc.

Referenced by produce(), and TtFullHadKinFitProducer().

◆ maxNJets_

int TtFullHadKinFitProducer::maxNJets_

private

maximal number of jets (-1 possible to indicate 'all')

Definition at line 49 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ maxNrIter_

unsigned int TtFullHadKinFitProducer::maxNrIter_

private

maximal number of iterations to be performed for the fit

Definition at line 53 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ minBTagValueBJet_

double TtFullHadKinFitProducer::minBTagValueBJet_

private

min value of bTag for a b-jet

Definition at line 39 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ mTop_

double TtFullHadKinFitProducer::mTop_

private

top mass value used for constraints

Definition at line 65 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ mW_

double TtFullHadKinFitProducer::mW_

private

W mass value used for constraints.

Definition at line 63 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ udscResolutions_

std::vector<edm::ParameterSet> TtFullHadKinFitProducer::udscResolutions_

private

store the resolutions for the jets

Definition at line 67 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ useBTagging_

bool TtFullHadKinFitProducer::useBTagging_

private

switch to tell whether to use b-tagging or not

Definition at line 43 of file TtFullHadKinFitProducer.cc.

Referenced by TtFullHadKinFitProducer().

◆ useOnlyMatch_

bool TtFullHadKinFitProducer::useOnlyMatch_

private

switch to tell whether all possible combinations should be used for the fit or only a certain combination

Definition at line 35 of file TtFullHadKinFitProducer.cc.

Referenced by produce().

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TtFullHadKinFitProducer()

Member Function Documentation

◆ produce()

Member Data Documentation

◆ bResolutions_

◆ bTagAlgo_

◆ bTags_

◆ constraints_

◆ jetCorrectionLevel_

◆ jetEnergyResolutionEtaBinning_

◆ jetEnergyResolutionScaleFactors_

◆ jetParam_

◆ jetsToken_

◆ kinFitter

◆ matchToken_

◆ maxBTagValueNonBJet_

◆ maxDeltaS_

◆ maxF_

◆ maxNComb_

◆ maxNJets_

◆ maxNrIter_

◆ minBTagValueBJet_

◆ mTop_

◆ mW_

◆ udscResolutions_

◆ useBTagging_

◆ useOnlyMatch_