class that does the fitting More...

#include <TtFullHadKinFitter.h>

Public Member Functions
std::list< TtFullHadKinFitter::KinFitResult >	fit (const std::vector< pat::Jet > &jets)
	do the fitting and return fit result More...

	KinFit ()
	default constructor More...

	KinFit (bool useBTagging, unsigned int bTags, std::string bTagAlgo, double minBTagValueBJet, double maxBTagValueNonBJet, const std::vector< edm::ParameterSet > &udscResolutions, const std::vector< edm::ParameterSet > &bResolutions, const std::vector< double > &jetEnergyResolutionScaleFactors, const std::vector< double > &jetEnergyResolutionEtaBinning, std::string jetCorrectionLevel, int maxNJets, int maxNComb, unsigned int maxNrIter, double maxDeltaS, double maxF, unsigned int jetParam, const std::vector< unsigned > &constraints, double mW, double mTop)
	special constructor More...

void	setBTagging (bool useBTagging, unsigned int bTags, std::string bTagAlgo, double minBTagValueBJet, double maxBTagValueNonBJet)
	set all parameters for b-tagging More...

void	setFitter (int maxNJets, unsigned int maxNrIter, double maxDeltaS, double maxF, unsigned int jetParam, const std::vector< unsigned > &constraints, double mW, double mTop)
	set parameters for fitter More...

void	setJEC (std::string jetCorrectionLevel)
	set jec level More...

void	setMatch (const std::vector< int > &match)
	set match to be used More...

void	setMatchInvalidity (bool invalidMatch)
	set the validity of a match More...

void	setOutput (int maxNComb)
	set number of combinations of output More...

void	setResolutions (const std::vector< edm::ParameterSet > &udscResolutions, const std::vector< edm::ParameterSet > &bResolutions, const std::vector< double > &jetEnergyResolutionScaleFactors, const std::vector< double > &jetEnergyResolutionEtaBinning)
	set resolutions More...

void	setUseOnlyMatch (bool useOnlyMatch)
	set useOnlyMatch More...

	~KinFit ()
	default destructor More...

Private Member Functions
TtFullHadKinFitter::Constraint	constraint (unsigned int configParameter)

std::vector< TtFullHadKinFitter::Constraint >	constraints (const std::vector< unsigned int > &configParameters)

pat::Jet	corJet (const pat::Jet &jet, const std::string &quarkType)
	helper function to construct the proper corrected jet for its corresponding quarkType More...

bool	doBTagging (const std::vector< pat::Jet > &jets, const unsigned int &bJetCounter, std::vector< int > &combi)

TtFullHadKinFitter::Param	param (unsigned int configParameter)

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

TtFullHadKinFitter *	fitter
	kinematic fit interface More...

bool	invalidMatch_
	match is invalid More...

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

std::vector< double >	jetEnergyResolutionEtaBinning_

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

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

std::vector< int >	match_
	the combination that should be used 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_

bool	useOnlyMatch_
	fit or only a certain combination More...

Detailed Description

class that does the fitting

Definition at line 145 of file TtFullHadKinFitter.h.

Constructor & Destructor Documentation

◆ KinFit() [1/2]

TtFullHadKinFitter::KinFit::KinFit ( )

default constructor

Definition at line 343 of file TtFullHadKinFitter.cc.

     : useBTagging_(true),
       bTags_(2),
       bTagAlgo_("trackCountingHighPurBJetTags"),
       minBTagValueBJet_(3.41),
       maxBTagValueNonBJet_(3.41),
       udscResolutions_(std::vector<edm::ParameterSet>(0)),
       bResolutions_(std::vector<edm::ParameterSet>(0)),
       jetEnergyResolutionScaleFactors_(0),
       jetEnergyResolutionEtaBinning_(0),
       jetCorrectionLevel_("L3Absolute"),
       maxNJets_(-1),
       maxNComb_(1),
       maxNrIter_(500),
       maxDeltaS_(5e-5),
       maxF_(0.0001),
       jetParam_(1),
       mW_(80.4),
       mTop_(173.),
       useOnlyMatch_(false),
       match_(std::vector<int>(0)),
       invalidMatch_(false) {
   constraints_.push_back(1);
   constraints_.push_back(2);
   constraints_.push_back(5);
 }

References constraints_.

◆ KinFit() [2/2]

TtFullHadKinFitter::KinFit::KinFit	(	bool	useBTagging,
		unsigned int	bTags,
		std::string	bTagAlgo,
		double	minBTagValueBJet,
		double	maxBTagValueNonBJet,
		const std::vector< edm::ParameterSet > &	udscResolutions,
		const std::vector< edm::ParameterSet > &	bResolutions,
		const std::vector< double > &	jetEnergyResolutionScaleFactors,
		const std::vector< double > &	jetEnergyResolutionEtaBinning,
		std::string	jetCorrectionLevel,
		int	maxNJets,
		int	maxNComb,
		unsigned int	maxNrIter,
		double	maxDeltaS,
		double	maxF,
		unsigned int	jetParam,
		const std::vector< unsigned > &	constraints,
		double	mW,
		double	mTop
	)

special constructor

Definition at line 371 of file TtFullHadKinFitter.cc.

     : useBTagging_(useBTagging),
       bTags_(bTags),
       bTagAlgo_(bTagAlgo),
       minBTagValueBJet_(minBTagValueBJet),
       maxBTagValueNonBJet_(maxBTagValueNonBJet),
       udscResolutions_(udscResolutions),
       bResolutions_(bResolutions),
       jetEnergyResolutionScaleFactors_(jetEnergyResolutionScaleFactors),
       jetEnergyResolutionEtaBinning_(jetEnergyResolutionEtaBinning),
       jetCorrectionLevel_(jetCorrectionLevel),
       maxNJets_(maxNJets),
       maxNComb_(maxNComb),
       maxNrIter_(maxNrIter),
       maxDeltaS_(maxDeltaS),
       maxF_(maxF),
       jetParam_(jetParam),
       constraints_(constraints),
       mW_(mW),
       mTop_(mTop),
       useOnlyMatch_(false),
       invalidMatch_(false) {
   // define kinematic fit interface
   fitter = new TtFullHadKinFitter(param(jetParam_),
                                   maxNrIter_,
                                   maxDeltaS_,
                                   maxF_,
                                   TtFullHadKinFitter::KinFit::constraints(constraints_),
                                   mW_,
                                   mTop_,
                                   &udscResolutions_,
                                   &bResolutions_,
                                   &jetEnergyResolutionScaleFactors_,
                                   &jetEnergyResolutionEtaBinning_);
 }

References bResolutions_, constraints(), constraints_, fitter, jetEnergyResolutionEtaBinning_, jetEnergyResolutionScaleFactors_, jetParam_, maxDeltaS_, maxF_, maxNrIter_, mTop_, mW_, param(), TtFullHadKinFitter::TtFullHadKinFitter(), and udscResolutions_.

◆ ~KinFit()

TtFullHadKinFitter::KinFit::~KinFit ( )

default destructor

Definition at line 426 of file TtFullHadKinFitter.cc.

426 { delete fitter; }

Member Function Documentation

◆ constraint()

TtFullHadKinFitter::Constraint TtFullHadKinFitter::KinFit::constraint ( unsigned int configParameter )

private

Definition at line 671 of file TtFullHadKinFitter.cc.

                                                                                           {
   TtFullHadKinFitter::Constraint result;
   switch (configParameter) {
     case TtFullHadKinFitter::kWPlusMass:
       result = TtFullHadKinFitter::kWPlusMass;
       break;
     case TtFullHadKinFitter::kWMinusMass:
       result = TtFullHadKinFitter::kWMinusMass;
       break;
     case TtFullHadKinFitter::kTopMass:
       result = TtFullHadKinFitter::kTopMass;
       break;
     case TtFullHadKinFitter::kTopBarMass:
       result = TtFullHadKinFitter::kTopBarMass;
       break;
     case TtFullHadKinFitter::kEqualTopMasses:
       result = TtFullHadKinFitter::kEqualTopMasses;
       break;
     default:
       throw cms::Exception("Configuration") << "Chosen fit constraint is not supported: " << configParameter << "\n";
       break;
   }
   return result;
 }

References Exception, TtFullHadKinFitter::kEqualTopMasses, TtFullHadKinFitter::kTopBarMass, TtFullHadKinFitter::kTopMass, TtFullHadKinFitter::kWMinusMass, TtFullHadKinFitter::kWPlusMass, and mps_fire::result.

◆ constraints()

std::vector< TtFullHadKinFitter::Constraint > TtFullHadKinFitter::KinFit::constraints ( const std::vector< unsigned int > & configParameters )

private

Definition at line 696 of file TtFullHadKinFitter.cc.

                                                  {
   std::vector<TtFullHadKinFitter::Constraint> result;
   for (unsigned i = 0; i < configParameters.size(); ++i) {
     result.push_back(constraint(configParameters[i]));
   }
   return result;
 }

References HLT_FULL_cff::constraint, mps_fire::i, and mps_fire::result.

Referenced by KinFit(), and setFitter().

◆ corJet()

pat::Jet TtFullHadKinFitter::KinFit::corJet	(	const pat::Jet &	jet,
		const std::string &	quarkType
	)

private

helper function to construct the proper corrected jet for its corresponding quarkType

Definition at line 485 of file TtFullHadKinFitter.cc.

                                                                                      {
   // jetCorrectionLevel was not configured
   if (jetCorrectionLevel_.empty())
     throw cms::Exception("Configuration")
         << "Unconfigured jetCorrectionLevel. Please use an appropriate, non-empty string.\n";
  
   // quarkType is unknown
   if (!(quarkType == "wMix" || quarkType == "uds" || quarkType == "charm" || quarkType == "bottom"))
     throw cms::Exception("Configuration") << quarkType << " is unknown as a quarkType for the jetCorrectionLevel.\n";
  
   float jecFactor = 1.;
   if (quarkType == "wMix")
     jecFactor = 0.75 * jet.jecFactor(jetCorrectionLevel_, "uds") + 0.25 * jet.jecFactor(jetCorrectionLevel_, "charm");
   else
     jecFactor = jet.jecFactor(jetCorrectionLevel_, quarkType);
  
   pat::Jet ret = jet;
   ret.setP4(ret.p4() * jecFactor);
   return ret;
 }

References metsig::jet, and runTheMatrix::ret.

◆ doBTagging()

bool TtFullHadKinFitter::KinFit::doBTagging	(	const std::vector< pat::Jet > &	jets,
		const unsigned int &	bJetCounter,
		std::vector< int > &	combi
	)

private

Definition at line 428 of file TtFullHadKinFitter.cc.

                                                                    {
   if (!useBTagging_) {
     return true;
   }
   if (bTags_ == 2 && jets[combi[TtFullHadEvtPartons::B]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
       jets[combi[TtFullHadEvtPartons::BBar]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
       jets[combi[TtFullHadEvtPartons::LightQ]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
       jets[combi[TtFullHadEvtPartons::LightQBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
       jets[combi[TtFullHadEvtPartons::LightP]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
       jets[combi[TtFullHadEvtPartons::LightPBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_) {
     return true;
   } else if (bTags_ == 1) {
     if (bJetCounter == 1 &&
         (jets[combi[TtFullHadEvtPartons::B]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ ||
          jets[combi[TtFullHadEvtPartons::BBar]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_) &&
         jets[combi[TtFullHadEvtPartons::LightQ]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
         jets[combi[TtFullHadEvtPartons::LightQBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
         jets[combi[TtFullHadEvtPartons::LightP]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
         jets[combi[TtFullHadEvtPartons::LightPBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_) {
       return true;
     } else if (bJetCounter > 1 && jets[combi[TtFullHadEvtPartons::B]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
                jets[combi[TtFullHadEvtPartons::BBar]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightQ]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightQBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightP]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightPBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_) {
       return true;
     }
   } else if (bTags_ == 0) {
     if (bJetCounter == 0) {
       return true;
     } else if (bJetCounter == 1 &&
                (jets[combi[TtFullHadEvtPartons::B]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ ||
                 jets[combi[TtFullHadEvtPartons::BBar]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_) &&
                jets[combi[TtFullHadEvtPartons::LightQ]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightQBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightP]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightPBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_) {
       return true;
     } else if (bJetCounter > 1 && jets[combi[TtFullHadEvtPartons::B]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
                jets[combi[TtFullHadEvtPartons::BBar]].bDiscriminator(bTagAlgo_) >= minBTagValueBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightQ]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightQBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightP]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_ &&
                jets[combi[TtFullHadEvtPartons::LightPBar]].bDiscriminator(bTagAlgo_) < maxBTagValueNonBJet_) {
       return true;
     }
   } else if (bTags_ > 2) {
     throw cms::Exception("Configuration") << "Wrong number of bTags (" << bTags_ << " bTags not supported)!\n";
     return true;
   }
   return false;
 }

References TtFullHadEvtPartons::B, TtFullHadEvtPartons::BBar, Exception, singleTopDQM_cfi::jets, TtFullHadEvtPartons::LightP, TtFullHadEvtPartons::LightPBar, TtFullHadEvtPartons::LightQ, and TtFullHadEvtPartons::LightQBar.

◆ fit()

std::list< TtFullHadKinFitter::KinFitResult > TtFullHadKinFitter::KinFit::fit ( const std::vector< pat::Jet > & jets )

do the fitting and return fit result

skip events with less jets than partons or invalid match

analyze different jet combinations using the KinFitter (or only a given jet combination if useOnlyMatch=true)

feed out result starting with the JetComb having the smallest chi2

Definition at line 506 of file TtFullHadKinFitter.cc.

                                                                                                    {
   std::list<TtFullHadKinFitter::KinFitResult> fitResults;
  
   if (jets.size() < nPartons || invalidMatch_) {
     // indices referring to the jet combination
     std::vector<int> invalidCombi;
     for (unsigned int i = 0; i < nPartons; ++i)
       invalidCombi.push_back(-1);
  
     KinFitResult result;
     // status of the fitter
     result.Status = -1;
     // chi2
     result.Chi2 = -1.;
     // chi2 probability
     result.Prob = -1.;
     // the kinFit getters return empty objects here
     result.B = fitter->fittedB();
     result.BBar = fitter->fittedBBar();
     result.LightQ = fitter->fittedLightQ();
     result.LightQBar = fitter->fittedLightQBar();
     result.LightP = fitter->fittedLightP();
     result.LightPBar = fitter->fittedLightPBar();
     result.JetCombi = invalidCombi;
     // push back fit result
     fitResults.push_back(result);
     return fitResults;
   }
  
   std::vector<int> jetIndices;
   if (!useOnlyMatch_) {
     for (unsigned int idx = 0; idx < jets.size(); ++idx) {
       if (maxNJets_ >= (int)nPartons && maxNJets_ == (int)idx)
         break;
       jetIndices.push_back(idx);
     }
   }
  
   std::vector<int> combi;
   for (unsigned int idx = 0; idx < nPartons; ++idx) {
     useOnlyMatch_ ? combi.push_back(match_[idx]) : combi.push_back(idx);
   }
  
   unsigned int bJetCounter = 0;
   for (std::vector<pat::Jet>::const_iterator jet = jets.begin(); jet < jets.end(); ++jet) {
     if (jet->bDiscriminator(bTagAlgo_) >= minBTagValueBJet_)
       ++bJetCounter;
   }
  
   do {
     for (int cnt = 0; cnt < TMath::Factorial(combi.size()); ++cnt) {
       // take into account indistinguishability of the two jets from the two W decays,
       // and the two decay branches, this reduces the combinatorics by a factor of 2*2*2
       if (((combi[TtFullHadEvtPartons::LightQ] < combi[TtFullHadEvtPartons::LightQBar] &&
             combi[TtFullHadEvtPartons::LightP] < combi[TtFullHadEvtPartons::LightPBar] &&
             combi[TtFullHadEvtPartons::B] < combi[TtFullHadEvtPartons::BBar]) ||
            useOnlyMatch_) &&
           doBTagging(jets, bJetCounter, combi)) {
         std::vector<pat::Jet> jetCombi;
         jetCombi.resize(nPartons);
         jetCombi[TtFullHadEvtPartons::LightQ] = corJet(jets[combi[TtFullHadEvtPartons::LightQ]], "wMix");
         jetCombi[TtFullHadEvtPartons::LightQBar] = corJet(jets[combi[TtFullHadEvtPartons::LightQBar]], "wMix");
         jetCombi[TtFullHadEvtPartons::B] = corJet(jets[combi[TtFullHadEvtPartons::B]], "bottom");
         jetCombi[TtFullHadEvtPartons::BBar] = corJet(jets[combi[TtFullHadEvtPartons::BBar]], "bottom");
         jetCombi[TtFullHadEvtPartons::LightP] = corJet(jets[combi[TtFullHadEvtPartons::LightP]], "wMix");
         jetCombi[TtFullHadEvtPartons::LightPBar] = corJet(jets[combi[TtFullHadEvtPartons::LightPBar]], "wMix");
  
         // do the kinematic fit
         int status = fitter->fit(jetCombi);
  
         if (status == 0) {
           // fill struct KinFitResults if converged
           TtFullHadKinFitter::KinFitResult result;
           result.Status = status;
           result.Chi2 = fitter->fitS();
           result.Prob = fitter->fitProb();
           result.B = fitter->fittedB();
           result.BBar = fitter->fittedBBar();
           result.LightQ = fitter->fittedLightQ();
           result.LightQBar = fitter->fittedLightQBar();
           result.LightP = fitter->fittedLightP();
           result.LightPBar = fitter->fittedLightPBar();
           result.JetCombi = combi;
           // push back fit result
           fitResults.push_back(result);
         }
       }
       // don't go through combinatorics if useOnlyMatch was chosen
       if (useOnlyMatch_) {
         break;
       }
       // next permutation
       std::next_permutation(combi.begin(), combi.end());
     }
     // don't go through combinatorics if useOnlyMatch was chosen
     if (useOnlyMatch_) {
       break;
     }
   } while (stdcomb::next_combination(jetIndices.begin(), jetIndices.end(), combi.begin(), combi.end()));
  
   // sort results w.r.t. chi2 values
   fitResults.sort();
  
   if ((unsigned)fitResults.size() < 1) {
     // in case no fit results were stored in the list (i.e. when all fits were aborted)
  
     KinFitResult result;
     // status of the fitter
     result.Status = -1;
     // chi2
     result.Chi2 = -1.;
     // chi2 probability
     result.Prob = -1.;
     // the kinFit getters return empty objects here
     result.B = fitter->fittedB();
     result.BBar = fitter->fittedBBar();
     result.LightQ = fitter->fittedLightQ();
     result.LightQBar = fitter->fittedLightQBar();
     result.LightP = fitter->fittedLightP();
     result.LightPBar = fitter->fittedLightPBar();
     // indices referring to the jet combination
     std::vector<int> invalidCombi(nPartons, -1);
     result.JetCombi = invalidCombi;
     // push back fit result
     fitResults.push_back(result);
   }
   return fitResults;
 }

References TtFullHadEvtPartons::B, TtFullHadEvtPartons::BBar, mps_fire::i, heavyIonCSV_trainingSettings::idx, createfilelist::int, metsig::jet, singleTopDQM_cfi::jets, TtFullHadEvtPartons::LightP, TtFullHadEvtPartons::LightPBar, TtFullHadEvtPartons::LightQ, TtFullHadEvtPartons::LightQBar, stdcomb::next_combination(), nPartons, mps_fire::result, and mps_update::status.

Referenced by trackingPlots.Iteration::modules(), and TtFullHadKinFitProducer::produce().

◆ param()

TtFullHadKinFitter::Param TtFullHadKinFitter::KinFit::param ( unsigned int configParameter )

private

Definition at line 651 of file TtFullHadKinFitter.cc.

                                                                                     {
   TtFullHadKinFitter::Param result;
   switch (configParameter) {
     case TtFullHadKinFitter::kEMom:
       result = TtFullHadKinFitter::kEMom;
       break;
     case TtFullHadKinFitter::kEtEtaPhi:
       result = TtFullHadKinFitter::kEtEtaPhi;
       break;
     case TtFullHadKinFitter::kEtThetaPhi:
       result = TtFullHadKinFitter::kEtThetaPhi;
       break;
     default:
       throw cms::Exception("Configuration")
           << "Chosen jet parametrization is not supported: " << configParameter << "\n";
       break;
   }
   return result;
 }

References Exception, TopKinFitter::kEMom, TopKinFitter::kEtEtaPhi, TopKinFitter::kEtThetaPhi, and mps_fire::result.

Referenced by KinFit().

◆ setBTagging()

void TtFullHadKinFitter::KinFit::setBTagging	(	bool	useBTagging,
		unsigned int	bTags,
		std::string	bTagAlgo,
		double	minBTagValueBJet,
		double	maxBTagValueNonBJet
	)

inline

set all parameters for b-tagging

Definition at line 173 of file TtFullHadKinFitter.h.

                                                  {
       useBTagging_ = useBTagging;
       bTags_ = bTags;
       bTagAlgo_ = bTagAlgo;
       minBTagValueBJet_ = minBTagValueBJet;
       maxBTagValueNonBJet_ = maxBTagValueNonBJet;
     }

References TtSemiLepHitFitProducer_Electrons_cfi::bTagAlgo, bTagAlgo_, TtFullHadKinFitProducer_cfi::bTags, bTags_, TtFullHadKinFitProducer_cfi::maxBTagValueNonBJet, maxBTagValueNonBJet_, TtFullHadKinFitProducer_cfi::minBTagValueBJet, minBTagValueBJet_, TtSemiLepHitFitProducer_Electrons_cfi::useBTagging, and useBTagging_.

◆ setFitter()

void TtFullHadKinFitter::KinFit::setFitter	(	int	maxNJets,
		unsigned int	maxNrIter,
		double	maxDeltaS,
		double	maxF,
		unsigned int	jetParam,
		const std::vector< unsigned > &	constraints,
		double	mW,
		double	mTop
	)

inline

set parameters for fitter

Definition at line 195 of file TtFullHadKinFitter.h.

                                 {
       maxNJets_ = maxNJets;
       maxNrIter_ = maxNrIter;
       maxDeltaS_ = maxDeltaS;
       maxF_ = maxF;
       jetParam_ = jetParam;
       constraints_ = constraints;
       mW_ = mW;
       mTop_ = mTop;
     }

References constraints(), constraints_, jetParam_, StEvtSolProducer_cfi::maxDeltaS, maxDeltaS_, StEvtSolProducer_cfi::maxF, maxF_, SUSYDQMAnalyzer_cfi::maxNJets, maxNJets_, StEvtSolProducer_cfi::maxNrIter, maxNrIter_, TtSemiLepHitFitProducer_Electrons_cfi::mTop, mTop_, TtSemiLepHitFitProducer_Electrons_cfi::mW, and mW_.

◆ setJEC()

void TtFullHadKinFitter::KinFit::setJEC ( std::string jetCorrectionLevel )

inline

set jec level

Definition at line 213 of file TtFullHadKinFitter.h.

213 { jetCorrectionLevel_ = jetCorrectionLevel; }

References TtSemiLepHitFitProducer_Electrons_cfi::jetCorrectionLevel, and jetCorrectionLevel_.

◆ setMatch()

void TtFullHadKinFitter::KinFit::setMatch ( const std::vector< int > & match )

inline

set match to be used

Definition at line 217 of file TtFullHadKinFitter.h.

217 { match_ = match; }

References match(), and match_.

Referenced by TtFullHadKinFitProducer::produce().

◆ setMatchInvalidity()

void TtFullHadKinFitter::KinFit::setMatchInvalidity ( bool invalidMatch )

inline

set the validity of a match

Definition at line 219 of file TtFullHadKinFitter.h.

219 { invalidMatch_ = invalidMatch; }

References invalidMatch_.

Referenced by TtFullHadKinFitProducer::produce().

◆ setOutput()

void TtFullHadKinFitter::KinFit::setOutput ( int maxNComb )

inline

set number of combinations of output

Definition at line 221 of file TtFullHadKinFitter.h.

221 { maxNComb_ = maxNComb; }

References TtSemiLepHitFitProducer_Electrons_cfi::maxNComb, and maxNComb_.

◆ setResolutions()

void TtFullHadKinFitter::KinFit::setResolutions	(	const std::vector< edm::ParameterSet > &	udscResolutions,
		const std::vector< edm::ParameterSet > &	bResolutions,
		const std::vector< double > &	jetEnergyResolutionScaleFactors,
		const std::vector< double > &	jetEnergyResolutionEtaBinning
	)

inline

set resolutions

Definition at line 185 of file TtFullHadKinFitter.h.

                                                                                 {
       udscResolutions_ = udscResolutions;
       bResolutions_ = bResolutions;
       jetEnergyResolutionScaleFactors_ = jetEnergyResolutionScaleFactors;
       jetEnergyResolutionEtaBinning_ = jetEnergyResolutionEtaBinning;
     }

References TtFullHadKinFitProducer_cfi::bResolutions, bResolutions_, TtFullHadKinFitProducer_cfi::jetEnergyResolutionEtaBinning, jetEnergyResolutionEtaBinning_, TtFullHadKinFitProducer_cfi::jetEnergyResolutionScaleFactors, jetEnergyResolutionScaleFactors_, TtFullHadKinFitProducer_cfi::udscResolutions, and udscResolutions_.

◆ setUseOnlyMatch()

void TtFullHadKinFitter::KinFit::setUseOnlyMatch ( bool useOnlyMatch )

inline

set useOnlyMatch

Definition at line 215 of file TtFullHadKinFitter.h.

215 { useOnlyMatch_ = useOnlyMatch; }

References TtFullHadKinFitProducer_cfi::useOnlyMatch, and useOnlyMatch_.

Referenced by TtFullHadKinFitProducer::produce().

Member Data Documentation

◆ bResolutions_

std::vector<edm::ParameterSet> TtFullHadKinFitter::KinFit::bResolutions_

private

Definition at line 252 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setResolutions().

◆ bTagAlgo_

std::string TtFullHadKinFitter::KinFit::bTagAlgo_

private

input tag for b-tagging algorithm

Definition at line 246 of file TtFullHadKinFitter.h.

Referenced by setBTagging().

◆ bTags_

unsigned int TtFullHadKinFitter::KinFit::bTags_

private

minimal number of b-jets

Definition at line 244 of file TtFullHadKinFitter.h.

Referenced by setBTagging().

◆ constraints_

std::vector<unsigned> TtFullHadKinFitter::KinFit::constraints_

private

numbering of different possible kinematic constraints

Definition at line 271 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ fitter

TtFullHadKinFitter* TtFullHadKinFitter::KinFit::fitter

private

kinematic fit interface

Definition at line 284 of file TtFullHadKinFitter.h.

Referenced by KinFit().

◆ invalidMatch_

bool TtFullHadKinFitter::KinFit::invalidMatch_

private

match is invalid

Definition at line 281 of file TtFullHadKinFitter.h.

Referenced by setMatchInvalidity().

◆ jetCorrectionLevel_

std::string TtFullHadKinFitter::KinFit::jetCorrectionLevel_

private

correction level for jets

Definition at line 257 of file TtFullHadKinFitter.h.

Referenced by setJEC().

◆ jetEnergyResolutionEtaBinning_

std::vector<double> TtFullHadKinFitter::KinFit::jetEnergyResolutionEtaBinning_

private

Definition at line 255 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setResolutions().

◆ jetEnergyResolutionScaleFactors_

std::vector<double> TtFullHadKinFitter::KinFit::jetEnergyResolutionScaleFactors_

private

scale factors for the jet energy resolution

Definition at line 254 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setResolutions().

◆ jetParam_

unsigned int TtFullHadKinFitter::KinFit::jetParam_

private

numbering of different possible jet parametrizations

Definition at line 269 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ match_

std::vector<int> TtFullHadKinFitter::KinFit::match_

private

the combination that should be used

Definition at line 279 of file TtFullHadKinFitter.h.

Referenced by setMatch().

◆ maxBTagValueNonBJet_

double TtFullHadKinFitter::KinFit::maxBTagValueNonBJet_

private

max value of bTag for a non-b-jet

Definition at line 250 of file TtFullHadKinFitter.h.

Referenced by setBTagging().

◆ maxDeltaS_

double TtFullHadKinFitter::KinFit::maxDeltaS_

private

maximal chi2 equivalent

Definition at line 265 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ maxF_

double TtFullHadKinFitter::KinFit::maxF_

private

maximal deviation for contstraints

Definition at line 267 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ maxNComb_

int TtFullHadKinFitter::KinFit::maxNComb_

private

maximal number of combinations to be written to the event

Definition at line 261 of file TtFullHadKinFitter.h.

Referenced by setOutput().

◆ maxNJets_

int TtFullHadKinFitter::KinFit::maxNJets_

private

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

Definition at line 259 of file TtFullHadKinFitter.h.

Referenced by setFitter().

◆ maxNrIter_

unsigned int TtFullHadKinFitter::KinFit::maxNrIter_

private

maximal number of iterations to be performed for the fit

Definition at line 263 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ minBTagValueBJet_

double TtFullHadKinFitter::KinFit::minBTagValueBJet_

private

min value of bTag for a b-jet

Definition at line 248 of file TtFullHadKinFitter.h.

Referenced by setBTagging().

◆ mTop_

double TtFullHadKinFitter::KinFit::mTop_

private

top mass value used for constraints

Definition at line 275 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ mW_

double TtFullHadKinFitter::KinFit::mW_

private

W mass value used for constraints.

Definition at line 273 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setFitter().

◆ udscResolutions_

std::vector<edm::ParameterSet> TtFullHadKinFitter::KinFit::udscResolutions_

private

store the resolutions for the jets

Definition at line 252 of file TtFullHadKinFitter.h.

Referenced by KinFit(), and setResolutions().

◆ useBTagging_

bool TtFullHadKinFitter::KinFit::useBTagging_

private

switch to tell whether all possible combinations should be used for the switch to tell whether to use b-tagging or not

Definition at line 242 of file TtFullHadKinFitter.h.

Referenced by setBTagging().

◆ useOnlyMatch_

bool TtFullHadKinFitter::KinFit::useOnlyMatch_

private

fit or only a certain combination

Definition at line 277 of file TtFullHadKinFitter.h.

Referenced by setUseOnlyMatch().

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ KinFit() [1/2]

◆ KinFit() [2/2]

◆ ~KinFit()

Member Function Documentation

◆ constraint()

◆ constraints()

◆ corJet()

◆ doBTagging()

◆ fit()

skip events with less jets than partons or invalid match

◆ param()

◆ setBTagging()

◆ setFitter()

◆ setJEC()

◆ setMatch()

◆ setMatchInvalidity()

◆ setOutput()

◆ setResolutions()

◆ setUseOnlyMatch()

Member Data Documentation

◆ bResolutions_

◆ bTagAlgo_

◆ bTags_

◆ constraints_

◆ fitter

◆ invalidMatch_

◆ jetCorrectionLevel_

◆ jetEnergyResolutionEtaBinning_

◆ jetEnergyResolutionScaleFactors_

◆ jetParam_

◆ match_

◆ maxBTagValueNonBJet_

◆ maxDeltaS_

◆ maxF_

◆ maxNComb_

◆ maxNJets_

◆ maxNrIter_

◆ minBTagValueBJet_

◆ mTop_

◆ mW_

◆ udscResolutions_

◆ useBTagging_

◆ useOnlyMatch_