TtresNeutrinoBuilder.cxx File Reference

#include <TMinuit.h>
#include "TopEventReconstructionTools/TtresNeutrinoBuilder.h"

Include dependency graph for TtresNeutrinoBuilder.cxx:

Go to the source code of this file.

Functions
void	delta2_fcn (Int_t &, Double_t *, Double_t &, Double_t *, Int_t)

Function Documentation

delta2_fcn()

void delta2_fcn	(	Int_t &	,
		Double_t *	,
		Double_t &	f,
		Double_t *	par,
		Int_t
	)

Definition at line 39 of file TtresNeutrinoBuilder.cxx.

                                                                            {
  Double_t delta2 = 0;
  Double_t alpha = par[0];
  Double_t r = par[1];
  Double_t dphi = par[2];
  Double_t l_pt = par[3];
  Double_t l_m = par[4];
  Double_t n_px = par[5];
  Double_t n_py = par[6];
 
  r /= sqrt(l_pt * l_pt + l_m * l_m) - l_pt * cos(dphi + alpha);
  TLorentzVector* neut = new TLorentzVector(n_px, n_py, 0., 0.);
  neut->SetE(neut->P());
 
  TLorentzVector* neut_new =
    new TLorentzVector(r * neut->P() * cos(neut->Phi() + alpha), r * neut->P() * sin(neut->Phi() + alpha), 0., 0.);
  neut_new->SetE(neut_new->P());
 
  delta2 = pow((neut_new->Px() - neut->Px()), 2) + pow((neut_new->Py() - neut->Py()), 2);
  r *= sqrt(l_pt * l_pt + l_m * l_m) - l_pt * cos(dphi + alpha);
  delete neut;
  delete neut_new;
  f = delta2;
}