From 81ef6d2ecc53d1035208124e43361dd318785c21 Mon Sep 17 00:00:00 2001 From: LucasNL Date: Mon, 19 Jul 2021 18:58:11 -0300 Subject: [PATCH] Bug fix on minuit fitter --- TrackFitter.cxx | 12 +++++------- 1 file changed, 5 insertions(+), 7 deletions(-) diff --git a/TrackFitter.cxx b/TrackFitter.cxx index 2af1a0e..aa549af 100644 --- a/TrackFitter.cxx +++ b/TrackFitter.cxx @@ -273,12 +273,10 @@ void TrackFitter::MinuitFit(FT3Track& track) //} double arglist[100]; - arglist[0] = 0; + arglist[0] = 1; // set print level minuit->ExecuteCommand("SET PRINT", arglist, 2); - //eu posso simplesmente buscar os parametors e imprimir de novo - // minimize arglist[0] = 10000; // number of function calls arglist[1] = 0.001; // tolerance @@ -334,10 +332,10 @@ void myFitFcn(Int_t&, Double_t*, Double_t& fval, Double_t* p, Int_t) for (auto z : zPositionsMFT) { // Propagate to Z tempTrack.propagateParamToZhelix(z, fieldZ); - tmp = (TrackFitter::PosX[i] - tempTrack.getX()) / TrackFitter::ErrorsX[i]; - chi2 += tmp * tmp; - tmp = (TrackFitter::PosY[i] - tempTrack.getY()) / TrackFitter::ErrorsY[i]; - chi2 += tmp * tmp; + tmp = TrackFitter::PosX[i] - tempTrack.getX(); + chi2 += tmp*tmp / TrackFitter::ErrorsX[i]; + tmp = TrackFitter::PosY[i] - tempTrack.getY(); + chi2 += tmp*tmp / TrackFitter::ErrorsY[i]; i++; } fval = chi2;