Chi2LJets.cxx File Reference

#include "TopEventReconstructionTools/Chi2LJets.h"
#include "TMinuit.h"

Include dependency graph for Chi2LJets.cxx:

Go to the source code of this file.

Namespaces
	top
	TopConfig A simple configuration that is NOT a singleton.

Functions
void	fcn (int &, double *, double &result, double par[], int)
	this is where we write out chi2 More...

Function Documentation

fcn()

void fcn	(	int &	,
		double *	,
		double &	result,
		double	par[],
		int
	)

this is where we write out chi2

Definition at line 183 of file Chi2LJets.cxx.

                                                                 {
  const top::FitInfo* params = (top::FitInfo*) gMinuit->GetObjectFit();
 
  double mytopmass = par[0];
  double pz = par[1];
 
  TLorentzVector neutrino(params->neutrino_px, params->neutrino_py, pz, sqrt(
                            params->neutrino_px* params->neutrino_px + params->neutrino_py* params->neutrino_py +
                            pz* pz));
  double w_lep_term = ((params->lepton + neutrino).M() - top::Chi2LJets::input_wmass) / top::Chi2LJets::input_wsigma;
  double t_lep_term = ((params->lepton + neutrino + params->leptonic_b).M() - mytopmass) /
                      top::Chi2LJets::input_topsigma;
 
  double w_term = (params->hadronic_w_mass - top::Chi2LJets::input_wmass) / top::Chi2LJets::input_wsigma;
  double t_term = (params->hadronic_t_mass - mytopmass) / top::Chi2LJets::input_topsigma;
 
  result = w_lep_term * w_lep_term +
           t_lep_term * t_lep_term +
           w_term * w_term +
           t_term * t_term;
 
  //  cout << "SJH     fcn hadw " << fit->hadronicW << " " << TopWidth::input_wmass << " " << TopWidth::input_wsigma <<
  // endl;
  //  cout << "SJH     fcn hadt " << fit->hadronicTop << " " << mytopmass << " " << TopWidth::input_topsigma << endl;
}