Fatal message signaling the absence of calibration data

[u/mrichter/AliRoot.git] / STEER / AliTrackResidualsChi2.cxx

diff --git a/STEER/AliTrackResidualsChi2.cxx b/STEER/AliTrackResidualsChi2.cxx
index 575a8183b57eb7bc79f92a883a4f825bc96bfc44..794b3c93aed6f4d9f475965bfdec926eba9e123a 100644 (file)
--- a/STEER/AliTrackResidualsChi2.cxx
+++ b/STEER/AliTrackResidualsChi2.cxx
@@ -19,24 +19,26 @@
 //
 //-----------------------------------------------------------------
 
-#include <TMinuit.h>
+#include <TVirtualFitter.h>
 #include <TGeoMatrix.h>
 
 #include "AliAlignObj.h"
 #include "AliTrackPointArray.h"
 #include "AliTrackResidualsChi2.h"
 
-void Fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *x, Int_t iflag);
+void TrackResidualsChi2Fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *x, Int_t iflag);
+
 
 ClassImp(AliTrackResidualsChi2)
 
+
 //______________________________________________________________________________
-void Fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
+void TrackResidualsChi2Fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *par, Int_t iflag)
 {
   // This function is called by minuit
   // The corresponding member method is called
   // using SetObjectFit/GetObjectFit methods of TMinuit
-  AliTrackResidualsChi2* dummy = (AliTrackResidualsChi2 *)gMinuit->GetObjectFit();
+  AliTrackResidualsChi2* dummy = (AliTrackResidualsChi2 *)TVirtualFitter::GetFitter()->GetObjectFit();
   dummy->Chi2(npar, gin, f, par, iflag);
 }
 
@@ -45,37 +47,50 @@ Bool_t AliTrackResidualsChi2::Minimize()
 {
   // Implementation of Chi2 based minimization
   // of track residuala sum
+  // 
+  // CHOLM: Modified to use the TVirtualFitter interface only.  This
+  // makes it possible other fitters than TMinuit (say TMinuit2,
+  // TFumili) by simply changing ones configuration file, or setting
+  // it up in a script.  Much more robust and flexible. 
   Double_t arglist[10];
   Int_t ierflg = 0;
-  TMinuit *gMinuit = new TMinuit(6);  //initialize TMinuit
+  TVirtualFitter *fitter = TVirtualFitter::Fitter(this,6);  //initialize TMinuit
   arglist[0] = -1;
-  gMinuit->mnexcm("SET PRINT", arglist, 1, ierflg);
+  ierflg = fitter->ExecuteCommand("SET PRINT", arglist, 1);
 
-  gMinuit->SetObjectFit(this);
-  gMinuit->SetFCN(Fcn);
+  fitter->SetFCN(TrackResidualsChi2Fcn);
 
   arglist[0] = 1;
-  gMinuit->mnexcm("SET ERR", arglist ,1,ierflg);
+  ierflg = fitter->ExecuteCommand("SET ERR", arglist ,1);
 
   // Set starting values and step sizes for parameters
   Double_t pars[6] = {0,0,0,0,0,0};
+  for(Int_t i=0;i<6;i++) if(fBFixed[i]) pars[i]=fFixed[i];
   Double_t step[6] = {0.0001,0.0001,0.0001,0.0001,0.0001,0.0001};
-  gMinuit->mnparm(0, "dx", pars[0], step[0], 0,0,ierflg);
-  gMinuit->mnparm(1, "dy", pars[1], step[1], 0,0,ierflg);
-  gMinuit->mnparm(2, "dz", pars[2], step[2], 0,0,ierflg);
-  gMinuit->mnparm(3, "psi", pars[3], step[3], 0,0,ierflg);
-  gMinuit->mnparm(4, "theta", pars[4], step[4], 0,0,ierflg);
-  gMinuit->mnparm(5, "phi", pars[5], step[5], 0,0,ierflg);
+  ierflg = fitter->SetParameter(0, "dx", pars[0], step[0], 0,0);
+  ierflg = fitter->SetParameter(1, "dy", pars[1], step[1], 0,0);
+  ierflg = fitter->SetParameter(2, "dz", pars[2], step[2], 0,0);
+  ierflg = fitter->SetParameter(3, "psi", pars[3], step[3], 0,0);
+  ierflg = fitter->SetParameter(4, "theta", pars[4], step[4], 0,0);
+  ierflg = fitter->SetParameter(5, "phi", pars[5], step[5], 0,0);
+
+  // Fix parameters
+  if(fBFixed[0]) {printf("Fixing dx=%f\n",pars[0]); fitter->FixParameter(0);}
+  if(fBFixed[1]) {printf("Fixing dy=%f\n",pars[1]); fitter->FixParameter(1);}
+  if(fBFixed[2]) {printf("Fixing dz=%f\n",pars[2]); fitter->FixParameter(2);}
+  if(fBFixed[3]) {printf("Fixing psi=%f\n",pars[3]); fitter->FixParameter(3);}
+  if(fBFixed[4]) {printf("Fixing theta=%f\n",pars[4]); fitter->FixParameter(4);}
+  if(fBFixed[5]) {printf("Fixing phi=%f\n",pars[5]); fitter->FixParameter(5);}
 
   // Now ready for minimization step
   arglist[0] = 500;
   arglist[1] = 1.;
-  gMinuit->mnexcm("MIGRAD", arglist ,2,ierflg);
+  fitter->ExecuteCommand("MIGRAD", arglist ,2);
 
   // Print results
   Double_t amin,edm,errdef;
   Int_t nvpar,nparx,icstat;
-  gMinuit->mnstat(amin,edm,errdef,nvpar,nparx,icstat);
+  icstat = fitter->GetStats(amin,edm,errdef,nvpar,nparx);
   fChi2 = amin; fNdf -= nvpar;
 
   return kTRUE;
@@ -101,7 +116,7 @@ void AliTrackResidualsChi2::Chi2(Int_t & /* npar */, Double_t * /* gin */, Doubl
       fVolArray[itrack]->GetPoint(p1,ipoint);
       fAlignObj->Transform(p1);
       fTrackArray[itrack]->GetPoint(p2,ipoint);
-      Float_t residual = p2.GetResidual(p1,kFALSE);
+      Float_t residual = p2.GetResidual(p1,kTRUE);
       chi2 += residual;
       if (count) fNdf += 3;
     }