CMS 3D CMS Logo
List of all members | Public Member Functions | Private Member Functions | Private Attributes
StKinFitter Class Reference

#include <StKinFitter.h>

Inheritance diagram for StKinFitter:
TopKinFitter

Public Member Functions

StEvtSolution addKinFitInfo (StEvtSolution *asol)
 
 StKinFitter ()
 
 StKinFitter (int jetParam, int lepParam, int metParam, int maxNrIter, double maxDeltaS, double maxF, const std::vector< int > &constraints)
 
 StKinFitter (Param jetParam, Param lepParam, Param metParam, int maxNrIter, double maxDeltaS, double maxF, const std::vector< int > &constraints)
 
 ~StKinFitter ()
 
- Public Member Functions inherited from TopKinFitter
int fitNrIter () const
 return number of used iterations More...
 
double fitProb () const
 return fit probability More...
 
double fitS () const
 return chi2 of fit (not normalized to degrees of freedom) More...
 
void setVerbosity (const int verbosityLevel)
 allows to change the verbosity of the TKinFitter More...
 
 TopKinFitter (const int maxNrIter=200, const double maxDeltaS=5e-5, const double maxF=1e-4, const double mW=80.4, const double mTop=173.)
 default constructor More...
 
 ~TopKinFitter ()
 default destructor More...
 

Private Member Functions

void setupFitter ()
 

Private Attributes

TFitConstraintM * cons1_
 
TFitConstraintM * cons2_
 
TFitConstraintM * cons3_
 
std::vector< int > constraints_
 
TAbsFitParticle * fitBottom_
 
TAbsFitParticle * fitLepton_
 
TAbsFitParticle * fitLight_
 
TAbsFitParticle * fitNeutrino_
 
Param jetParam_
 
Param lepParam_
 
Param metParam_
 

Additional Inherited Members

- Public Types inherited from TopKinFitter
enum  Param { kEMom, kEtEtaPhi, kEtThetaPhi }
 supported parameterizations More...
 
- Protected Member Functions inherited from TopKinFitter
std::string param (const Param &param) const
 convert Param to human readable form More...
 
- Protected Attributes inherited from TopKinFitter
TKinFitter * fitter_
 kinematic fitter More...
 
double maxDeltaS_
 maximal allowed chi2 (not normalized to degrees of freedom) More...
 
double maxF_
 maximal allowed distance from constraints More...
 
int maxNrIter_
 maximal allowed number of iterations to be used for the fit More...
 
double mTop_
 top mass value used for constraints More...
 
double mW_
 W mass value used for constraints. More...
 

Detailed Description

Definition at line 19 of file StKinFitter.h.

Constructor & Destructor Documentation

◆ StKinFitter() [1/3]

StKinFitter::StKinFitter ( )

Definition at line 27 of file StKinFitter.cc.

27 : TopKinFitter(), jetParam_(kEMom), lepParam_(kEMom), metParam_(kEMom) { setupFitter(); }

References setupFitter().

◆ StKinFitter() [2/3]

StKinFitter::StKinFitter ( int  jetParam,
int  lepParam,
int  metParam,
int  maxNrIter,
double  maxDeltaS,
double  maxF,
const std::vector< int > &  constraints 
)

Definition at line 29 of file StKinFitter.cc.

36  : TopKinFitter(maxNrIter, maxDeltaS, maxF),
37  jetParam_((Param)jetParam),
38  lepParam_((Param)lepParam),
39  metParam_((Param)metParam),
40  constraints_(constraints) {
41  setupFitter();
42 }

References setupFitter().

◆ StKinFitter() [3/3]

StKinFitter::StKinFitter ( Param  jetParam,
Param  lepParam,
Param  metParam,
int  maxNrIter,
double  maxDeltaS,
double  maxF,
const std::vector< int > &  constraints 
)

◆ ~StKinFitter()

StKinFitter::~StKinFitter ( )

Definition at line 59 of file StKinFitter.cc.

59  {
60  delete cons1_;
61  delete cons2_;
62  delete cons3_;
63  delete fitBottom_;
64  delete fitLight_;
65  delete fitLepton_;
66  delete fitNeutrino_;
67 }

References cons1_, cons2_, cons3_, fitBottom_, fitLepton_, fitLight_, and fitNeutrino_.

Member Function Documentation

◆ addKinFitInfo()

StEvtSolution StKinFitter::addKinFitInfo ( StEvtSolution *  asol)

Definition at line 69 of file StKinFitter.cc.

69  {
70  StEvtSolution fitsol(*asol);
71 
72  TMatrixD m1(3, 3), m2(3, 3);
73  TMatrixD m1b(4, 4), m2b(4, 4);
74  TMatrixD m3(3, 3), m4(3, 3);
75  m1.Zero();
76  m2.Zero();
77  m1b.Zero();
78  m2b.Zero();
79  m3.Zero();
80  m4.Zero();
81 
82  TLorentzVector bottomVec(
83  fitsol.getBottom().px(), fitsol.getBottom().py(), fitsol.getBottom().pz(), fitsol.getBottom().energy());
84  TLorentzVector lightVec(
85  fitsol.getLight().px(), fitsol.getLight().py(), fitsol.getLight().pz(), fitsol.getLight().energy());
86  TLorentzVector leplVec;
87  if (fitsol.getDecay() == "electron")
88  leplVec = TLorentzVector(
89  fitsol.getElectron().px(), fitsol.getElectron().py(), fitsol.getElectron().pz(), fitsol.getElectron().energy());
90  if (fitsol.getDecay() == "muon")
91  leplVec =
92  TLorentzVector(fitsol.getMuon().px(), fitsol.getMuon().py(), fitsol.getMuon().pz(), fitsol.getMuon().energy());
93  TLorentzVector lepnVec(fitsol.getNeutrino().px(), fitsol.getNeutrino().py(), 0, fitsol.getNeutrino().et());
94 
95  // jet resolutions
96  {
97  //FIXME this dirty hack needs a clean solution soon!
98  double pt = fitsol.getBottom().pt();
99  double eta = fitsol.getBottom().eta();
100  res::HelperJet jetRes;
101  if (jetParam_ == kEMom) {
102  m1b(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
103  m1b(1, 1) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
104  m1b(2, 2) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
105  m1b(3, 3) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
106  m2b(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
107  m2b(1, 1) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
108  m2b(2, 2) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
109  m2b(3, 3) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
110  } else if (jetParam_ == kEtEtaPhi) {
111  m1(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
112  m1(1, 1) = pow(jetRes.eta(pt, eta, res::HelperJet::kB), 2);
113  m1(2, 2) = pow(jetRes.phi(pt, eta, res::HelperJet::kB), 2);
114  m2(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
115  m2(1, 1) = pow(jetRes.eta(pt, eta, res::HelperJet::kUds), 2);
116  m2(2, 2) = pow(jetRes.phi(pt, eta, res::HelperJet::kUds), 2);
117  } else if (jetParam_ == kEtThetaPhi) {
118  m1(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kB), 2);
119  m1(1, 1) = pow(jetRes.eta(pt, eta, res::HelperJet::kB), 2);
120  m1(2, 2) = pow(jetRes.phi(pt, eta, res::HelperJet::kB), 2);
121  m2(0, 0) = pow(jetRes.pt(pt, eta, res::HelperJet::kUds), 2);
122  m2(1, 1) = pow(jetRes.eta(pt, eta, res::HelperJet::kUds), 2);
123  m2(2, 2) = pow(jetRes.phi(pt, eta, res::HelperJet::kUds), 2);
124  }
125  }
126  // lepton resolutions
127  {
128  //FIXME this dirty hack needs a clean solution soon!
129  double pt = fitsol.getElectron().pt();
130  double eta = fitsol.getElectron().eta();
131  res::HelperMuon muonRes;
132  res::HelperElectron elecRes;
133  if (lepParam_ == kEMom) {
134  if (fitsol.getDecay() == "electron") {
135  m3(0, 0) = pow(elecRes.pt(pt, eta), 2);
136  m3(1, 1) = pow(elecRes.pt(pt, eta), 2);
137  m3(2, 2) = pow(elecRes.pt(pt, eta), 2);
138  }
139  if (fitsol.getDecay() == "muon") {
140  m3(0, 0) = pow(muonRes.pt(pt, eta), 2);
141  m3(1, 1) = pow(muonRes.pt(pt, eta), 2);
142  m3(2, 2) = pow(muonRes.pt(pt, eta), 2);
143  }
144  } else if (lepParam_ == kEtEtaPhi) {
145  if (fitsol.getDecay() == "electron") {
146  m3(0, 0) = pow(elecRes.pt(pt, eta), 2);
147  m3(1, 1) = pow(elecRes.eta(pt, eta), 2);
148  m3(2, 2) = pow(elecRes.phi(pt, eta), 2);
149  }
150  if (fitsol.getDecay() == "muon") {
151  m3(0, 0) = pow(muonRes.pt(pt, eta), 2);
152  m3(1, 1) = pow(muonRes.eta(pt, eta), 2);
153  m3(2, 2) = pow(muonRes.phi(pt, eta), 2);
154  }
155  } else if (lepParam_ == kEtThetaPhi) {
156  if (fitsol.getDecay() == "electron") {
157  m3(0, 0) = pow(elecRes.pt(pt, eta), 2);
158  m3(1, 1) = pow(elecRes.eta(pt, eta), 2);
159  m3(2, 2) = pow(elecRes.phi(pt, eta), 2);
160  }
161  if (fitsol.getDecay() == "muon") {
162  m3(0, 0) = pow(muonRes.pt(pt, eta), 2);
163  m3(1, 1) = pow(muonRes.eta(pt, eta), 2);
164  m3(2, 2) = pow(muonRes.phi(pt, eta), 2);
165  }
166  }
167  }
168  // neutrino resolutions
169  {
170  //FIXME this dirty hack needs a clean solution soon!
171  double met = fitsol.getNeutrino().pt();
172  res::HelperMET metRes;
173  if (metParam_ == kEMom) {
174  m4(0, 0) = pow(metRes.met(met), 2);
175  m4(1, 1) = pow(9999., 2);
176  m4(2, 2) = pow(metRes.met(met), 2);
177  } else if (metParam_ == kEtEtaPhi) {
178  m4(0, 0) = pow(metRes.met(met), 2);
179  m4(1, 1) = pow(9999., 2);
180  m4(2, 2) = pow(metRes.phi(met), 2);
181  } else if (metParam_ == kEtThetaPhi) {
182  m4(0, 0) = pow(metRes.met(met), 2);
183  m4(1, 1) = pow(9999., 2);
184  m4(2, 2) = pow(metRes.phi(met), 2);
185  }
186  }
187  // set the kinematics of the objects to be fitted
188  fitBottom_->setIni4Vec(&bottomVec);
189  fitLight_->setIni4Vec(&lightVec);
190  fitLepton_->setIni4Vec(&leplVec);
191  fitNeutrino_->setIni4Vec(&lepnVec);
192  if (jetParam_ == kEMom) {
193  fitBottom_->setCovMatrix(&m1b);
194  fitLight_->setCovMatrix(&m2b);
195  } else {
196  fitBottom_->setCovMatrix(&m1);
197  fitLight_->setCovMatrix(&m2);
198  }
199  fitLepton_->setCovMatrix(&m3);
200  fitNeutrino_->setCovMatrix(&m4);
201 
202  // perform the fit!
203  fitter_->fit();
204 
205  // add fitted information to the solution
206  if (fitter_->getStatus() == 0) {
207  // read back the jet kinematics and resolutions
208  pat::Particle aFitBottom(reco::LeafCandidate(0,
209  math::XYZTLorentzVector(fitBottom_->getCurr4Vec()->X(),
210  fitBottom_->getCurr4Vec()->Y(),
211  fitBottom_->getCurr4Vec()->Z(),
212  fitBottom_->getCurr4Vec()->E()),
213  math::XYZPoint()));
214  pat::Particle aFitLight(reco::LeafCandidate(0,
215  math::XYZTLorentzVector(fitLight_->getCurr4Vec()->X(),
216  fitLight_->getCurr4Vec()->Y(),
217  fitLight_->getCurr4Vec()->Z(),
218  fitLight_->getCurr4Vec()->E()),
219  math::XYZPoint()));
220 
221  // read back the lepton kinematics and resolutions
222  pat::Particle aFitLepton(reco::LeafCandidate(0,
223  math::XYZTLorentzVector(fitLepton_->getCurr4Vec()->X(),
224  fitLepton_->getCurr4Vec()->Y(),
225  fitLepton_->getCurr4Vec()->Z(),
226  fitLepton_->getCurr4Vec()->E()),
227  math::XYZPoint()));
228 
229  // read back the MET kinematics and resolutions
230  pat::Particle aFitNeutrino(reco::LeafCandidate(0,
231  math::XYZTLorentzVector(fitNeutrino_->getCurr4Vec()->X(),
232  fitNeutrino_->getCurr4Vec()->Y(),
233  fitNeutrino_->getCurr4Vec()->Z(),
234  fitNeutrino_->getCurr4Vec()->E()),
235  math::XYZPoint()));
236 
237  // finally fill the fitted particles
238  fitsol.setFitBottom(aFitBottom);
239  fitsol.setFitLight(aFitLight);
240  fitsol.setFitLepton(aFitLepton);
241  fitsol.setFitNeutrino(aFitNeutrino);
242 
243  // store the fit's chi2 probability
244  fitsol.setChi2Prob(fitProb());
245  }
246 
247  return fitsol;
248 }

References reco::LeafCandidate::energy(), reco::LeafCandidate::et(), res::HelperElectron::eta(), res::HelperMuon::eta(), res::HelperJet::eta(), PVValHelper::eta, reco::LeafCandidate::eta(), TKinFitter::fit(), fitBottom_, fitLepton_, fitLight_, fitNeutrino_, TopKinFitter::fitProb(), TopKinFitter::fitter_, StEvtSolution::getBottom(), TAbsFitParticle::getCurr4Vec(), StEvtSolution::getDecay(), StEvtSolution::getElectron(), StEvtSolution::getLight(), StEvtSolution::getMuon(), StEvtSolution::getNeutrino(), TKinFitter::getStatus(), jetParam_, res::HelperJet::kB, TopKinFitter::kEMom, TopKinFitter::kEtEtaPhi, TopKinFitter::kEtThetaPhi, res::HelperJet::kUds, lepParam_, res::HelperMET::met(), BTaggingMonitor_cfi::met, metParam_, res::HelperMuon::phi(), res::HelperElectron::phi(), res::HelperJet::phi(), res::HelperMET::phi(), funct::pow(), res::HelperElectron::pt(), res::HelperMuon::pt(), res::HelperJet::pt(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), StEvtSolution::setChi2Prob(), TAbsFitParticle::setCovMatrix(), StEvtSolution::setFitBottom(), StEvtSolution::setFitLepton(), StEvtSolution::setFitLight(), StEvtSolution::setFitNeutrino(), and TAbsFitParticle::setIni4Vec().

Referenced by StEvtSolutionMaker::produce().

◆ setupFitter()

void StKinFitter::setupFitter ( )
private

Definition at line 253 of file StKinFitter.cc.

253  {
254  // FIXME: replace by messagelogger!!!
255 
256  std::cout << std::endl << std::endl << "+++++++++++ KINFIT SETUP ++++++++++++" << std::endl;
257  std::cout << " jet parametrisation: " << param(jetParam_) << std::endl;
258  std::cout << " lepton parametrisation: " << param(lepParam_) << std::endl;
259  std::cout << " met parametrisation: " << param(metParam_) << std::endl;
260  std::cout << " constraints: " << std::endl;
261  for (unsigned int i = 0; i < constraints_.size(); i++) {
262  if (constraints_[i] == 1)
263  std::cout << " - hadronic W-mass" << std::endl;
264  if (constraints_[i] == 2)
265  std::cout << " - leptonic W-mass" << std::endl;
266  if (constraints_[i] == 3)
267  std::cout << " - hadronic top mass" << std::endl;
268  if (constraints_[i] == 4)
269  std::cout << " - leptonic top mass" << std::endl;
270  if (constraints_[i] == 5)
271  std::cout << " - neutrino mass" << std::endl;
272  }
273  std::cout << "Max. number of iterations: " << maxNrIter_ << std::endl;
274  std::cout << "Max. deltaS: " << maxDeltaS_ << std::endl;
275  std::cout << "Max. F: " << maxF_ << std::endl;
276  std::cout << "++++++++++++++++++++++++++++++++++++++++++++" << std::endl << std::endl << std::endl;
277 
278  TMatrixD empty3(3, 3);
279  TMatrixD empty4(4, 4);
280  if (jetParam_ == kEMom) {
281  fitBottom_ = new TFitParticleEMomDev("Jet1", "Jet1", nullptr, &empty4);
282  fitLight_ = new TFitParticleEMomDev("Jet2", "Jet2", nullptr, &empty4);
283  } else if (jetParam_ == kEtEtaPhi) {
284  fitBottom_ = new TFitParticleEtEtaPhi("Jet1", "Jet1", nullptr, &empty3);
285  fitLight_ = new TFitParticleEtEtaPhi("Jet2", "Jet2", nullptr, &empty3);
286  } else if (jetParam_ == kEtThetaPhi) {
287  fitBottom_ = new TFitParticleEtThetaPhi("Jet1", "Jet1", nullptr, &empty3);
288  fitLight_ = new TFitParticleEtThetaPhi("Jet2", "Jet2", nullptr, &empty3);
289  }
290  if (lepParam_ == kEMom) {
291  fitLepton_ = new TFitParticleEScaledMomDev("Lepton", "Lepton", nullptr, &empty3);
292  } else if (lepParam_ == kEtEtaPhi) {
293  fitLepton_ = new TFitParticleEtEtaPhi("Lepton", "Lepton", nullptr, &empty3);
294  } else if (lepParam_ == kEtThetaPhi) {
295  fitLepton_ = new TFitParticleEtThetaPhi("Lepton", "Lepton", nullptr, &empty3);
296  }
297  if (metParam_ == kEMom) {
298  fitNeutrino_ = new TFitParticleEScaledMomDev("Neutrino", "Neutrino", nullptr, &empty3);
299  } else if (metParam_ == kEtEtaPhi) {
300  fitNeutrino_ = new TFitParticleEtEtaPhi("Neutrino", "Neutrino", nullptr, &empty3);
301  } else if (metParam_ == kEtThetaPhi) {
302  fitNeutrino_ = new TFitParticleEtThetaPhi("Neutrino", "Neutrino", nullptr, &empty3);
303  }
304 
305  cons1_ = new TFitConstraintM("MassConstraint", "Mass-Constraint", nullptr, nullptr, mW_);
306  cons1_->addParticles1(fitLepton_, fitNeutrino_);
307  cons2_ = new TFitConstraintM("MassConstraint", "Mass-Constraint", nullptr, nullptr, mTop_);
308  cons2_->addParticles1(fitLepton_, fitNeutrino_, fitBottom_);
309  cons3_ = new TFitConstraintM("MassConstraint", "Mass-Constraint", nullptr, nullptr, 0.);
310  cons3_->addParticle1(fitNeutrino_);
311 
312  for (unsigned int i = 0; i < constraints_.size(); i++) {
313  if (constraints_[i] == 1)
314  fitter_->addConstraint(cons1_);
315  if (constraints_[i] == 2)
316  fitter_->addConstraint(cons2_);
317  if (constraints_[i] == 3)
318  fitter_->addConstraint(cons3_);
319  }
320  fitter_->addMeasParticle(fitBottom_);
321  fitter_->addMeasParticle(fitLight_);
322  fitter_->addMeasParticle(fitLepton_);
323  fitter_->addMeasParticle(fitNeutrino_);
324 }

References TKinFitter::addConstraint(), TKinFitter::addMeasParticle(), TFitConstraintM::addParticle1(), TFitConstraintM::addParticles1(), cons1_, cons2_, cons3_, constraints_, gather_cfg::cout, fitBottom_, fitLepton_, fitLight_, fitNeutrino_, TopKinFitter::fitter_, mps_fire::i, jetParam_, TopKinFitter::kEMom, TopKinFitter::kEtEtaPhi, TopKinFitter::kEtThetaPhi, lepParam_, TopKinFitter::maxDeltaS_, TopKinFitter::maxF_, TopKinFitter::maxNrIter_, metParam_, TopKinFitter::mTop_, TopKinFitter::mW_, and TopKinFitter::param().

Referenced by StKinFitter().

Member Data Documentation

◆ cons1_

TFitConstraintM* StKinFitter::cons1_
private

Definition at line 50 of file StKinFitter.h.

Referenced by setupFitter(), and ~StKinFitter().

◆ cons2_

TFitConstraintM* StKinFitter::cons2_
private

Definition at line 51 of file StKinFitter.h.

Referenced by setupFitter(), and ~StKinFitter().

◆ cons3_

TFitConstraintM* StKinFitter::cons3_
private

Definition at line 52 of file StKinFitter.h.

Referenced by setupFitter(), and ~StKinFitter().

◆ constraints_

std::vector<int> StKinFitter::constraints_
private

Definition at line 55 of file StKinFitter.h.

Referenced by setupFitter().

◆ fitBottom_

TAbsFitParticle* StKinFitter::fitBottom_
private

Definition at line 45 of file StKinFitter.h.

Referenced by addKinFitInfo(), setupFitter(), and ~StKinFitter().

◆ fitLepton_

TAbsFitParticle* StKinFitter::fitLepton_
private

Definition at line 47 of file StKinFitter.h.

Referenced by addKinFitInfo(), setupFitter(), and ~StKinFitter().

◆ fitLight_

TAbsFitParticle* StKinFitter::fitLight_
private

Definition at line 46 of file StKinFitter.h.

Referenced by addKinFitInfo(), setupFitter(), and ~StKinFitter().

◆ fitNeutrino_

TAbsFitParticle* StKinFitter::fitNeutrino_
private

Definition at line 48 of file StKinFitter.h.

Referenced by addKinFitInfo(), setupFitter(), and ~StKinFitter().

◆ jetParam_

Param StKinFitter::jetParam_
private

Definition at line 54 of file StKinFitter.h.

Referenced by addKinFitInfo(), and setupFitter().

◆ lepParam_

Param StKinFitter::lepParam_
private

Definition at line 54 of file StKinFitter.h.

Referenced by addKinFitInfo(), and setupFitter().

◆ metParam_

Param StKinFitter::metParam_
private

Definition at line 54 of file StKinFitter.h.

Referenced by addKinFitInfo(), and setupFitter().

TopKinFitter::kEtEtaPhi
Definition: TopKinFitter.h:20
TFitParticleEMomDev
Definition: TFitParticleEMomDev.h:10
StKinFitter::cons2_
TFitConstraintM * cons2_
Definition: StKinFitter.h:51
res::HelperMuon
Definition: Muon.h:7
res::HelperMET
Definition: MET.h:7
mps_fire.i
i
Definition: mps_fire.py:428
TopKinFitter::fitter_
TKinFitter * fitter_
kinematic fitter
Definition: TopKinFitter.h:47
TopKinFitter::maxNrIter_
int maxNrIter_
maximal allowed number of iterations to be used for the fit
Definition: TopKinFitter.h:49
DiDispStaMuonMonitor_cfi.pt
pt
Definition: DiDispStaMuonMonitor_cfi.py:39
StEvtSolProducer_cfi.maxNrIter
maxNrIter
Definition: StEvtSolProducer_cfi.py:24
res::HelperJet
Definition: Jet.h:7
TFitConstraintM
Definition: TFitConstraintM.h:11
gather_cfg.cout
cout
Definition: gather_cfg.py:144
StKinFitter::setupFitter
void setupFitter()
Definition: StKinFitter.cc:253
TopKinFitter::param
std::string param(const Param &param) const
convert Param to human readable form
Definition: TopKinFitter.cc:18
res::HelperMET::phi
double phi(double pt)
Definition: MET.h:51
StKinFitter::fitLight_
TAbsFitParticle * fitLight_
Definition: StKinFitter.h:46
TFitConstraintM::addParticle1
void addParticle1(TAbsFitParticle *particle)
Definition: TFitConstraintM.cc:73
res::HelperJet::phi
double phi(double pt, double eta, Flavor flav)
Definition: Jet.h:282
res::HelperJet::pt
double pt(double pt, double eta, Flavor flav)
Definition: Jet.h:26
TopKinFitter::kEMom
Definition: TopKinFitter.h:20
res::HelperElectron::phi
double phi(double pt, double eta)
Definition: Electron.h:160
TopKinFitter::fitProb
double fitProb() const
return fit probability
Definition: TopKinFitter.h:37
StKinFitter::cons3_
TFitConstraintM * cons3_
Definition: StKinFitter.h:52
BTaggingMonitor_cfi.met
met
Definition: BTaggingMonitor_cfi.py:84
StKinFitter::fitNeutrino_
TAbsFitParticle * fitNeutrino_
Definition: StKinFitter.h:48
TopKinFitter::mW_
double mW_
W mass value used for constraints.
Definition: TopKinFitter.h:55
res::HelperElectron::eta
double eta(double pt, double eta)
Definition: Electron.h:210
TKinFitter::fit
Int_t fit()
Definition: TKinFitter.cc:318
res::HelperJet::kB
Definition: Jet.h:9
PVValHelper::eta
Definition: PVValidationHelpers.h:70
TopKinFitter::TopKinFitter
TopKinFitter(const int maxNrIter=200, const double maxDeltaS=5e-5, const double maxF=1e-4, const double mW=80.4, const double mTop=173.)
default constructor
Definition: TopKinFitter.cc:4
TopKinFitter::maxDeltaS_
double maxDeltaS_
maximal allowed chi2 (not normalized to degrees of freedom)
Definition: TopKinFitter.h:51
TFitConstraintM::addParticles1
void addParticles1(TAbsFitParticle *p1, TAbsFitParticle *p2=nullptr, TAbsFitParticle *p3=nullptr, TAbsFitParticle *p4=nullptr, TAbsFitParticle *p5=nullptr, TAbsFitParticle *p6=nullptr, TAbsFitParticle *p7=nullptr, TAbsFitParticle *p8=nullptr, TAbsFitParticle *p9=nullptr, TAbsFitParticle *p10=nullptr)
Definition: TFitConstraintM.cc:87
StKinFitter::metParam_
Param metParam_
Definition: StKinFitter.h:54
StKinFitter::constraints_
std::vector< int > constraints_
Definition: StKinFitter.h:55
StKinFitter::fitBottom_
TAbsFitParticle * fitBottom_
Definition: StKinFitter.h:45
res::HelperElectron::pt
double pt(double pt, double eta)
Definition: Electron.h:24
TAbsFitParticle::getCurr4Vec
const TLorentzVector * getCurr4Vec()
Definition: TAbsFitParticle.h:31
TKinFitter::getStatus
Int_t getStatus()
Definition: TKinFitter.h:51
math::XYZPoint
XYZPointD XYZPoint
point in space with cartesian internal representation
Definition: Point3D.h:12
res::HelperElectron
Definition: Electron.h:7
TAbsFitParticle::setCovMatrix
virtual void setCovMatrix(const TMatrixD *theCovMatrix)
Definition: TAbsFitParticle.cc:98
TopKinFitter::kEtThetaPhi
Definition: TopKinFitter.h:20
res::HelperMET::met
double met(double met)
Definition: MET.h:24
createBeamHaloJobs.constraints
string constraints
Definition: createBeamHaloJobs.py:229
StKinFitter::cons1_
TFitConstraintM * cons1_
Definition: StKinFitter.h:50
StEvtSolProducer_cfi.maxDeltaS
maxDeltaS
Definition: StEvtSolProducer_cfi.py:25
TopKinFitter::mTop_
double mTop_
top mass value used for constraints
Definition: TopKinFitter.h:57
TKinFitter::addConstraint
void addConstraint(TAbsFitConstraint *constraint)
Definition: TKinFitter.cc:292
TFitParticleEScaledMomDev
Definition: TFitParticleEScaledMomDev.h:8
StKinFitter::jetParam_
Param jetParam_
Definition: StKinFitter.h:54
res::HelperMuon::eta
double eta(double pt, double eta)
Definition: Muon.h:204
TopKinFitter::maxF_
double maxF_
maximal allowed distance from constraints
Definition: TopKinFitter.h:53
math::XYZTLorentzVector
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
Definition: LorentzVector.h:29
TFitParticleEtEtaPhi
Definition: TFitParticleEtEtaPhi.h:8
TKinFitter::addMeasParticle
void addMeasParticle(TAbsFitParticle *particle)
Definition: TKinFitter.cc:194
res::HelperJet::eta
double eta(double pt, double eta, Flavor flav)
Definition: Jet.h:378
res::HelperMuon::phi
double phi(double pt, double eta)
Definition: Muon.h:154
StKinFitter::lepParam_
Param lepParam_
Definition: StKinFitter.h:54
funct::pow
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:29
reco::LeafCandidate
Definition: LeafCandidate.h:16
TFitParticleEtThetaPhi
Definition: TFitParticleEtThetaPhi.h:8
StKinFitter::fitLepton_
TAbsFitParticle * fitLepton_
Definition: StKinFitter.h:47
TopKinFitter::Param
Param
supported parameterizations
Definition: TopKinFitter.h:20
TAbsFitParticle::setIni4Vec
virtual void setIni4Vec(const TLorentzVector *pini)=0
StEvtSolProducer_cfi.maxF
maxF
Definition: StEvtSolProducer_cfi.py:26
res::HelperJet::kUds
Definition: Jet.h:9
pat::Particle
Analysis-level particle class.
Definition: Particle.h:30
res::HelperMuon::pt
double pt(double pt, double eta)
Definition: Muon.h:24
StEvtSolution
Definition: StEvtSolution.h:21
Generated for CMSSW Reference Manual by doxygen 1.8.16