Base class (AliMisaligner); each detector will provide its derived class,

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

diff --git a/STEER/AliTrackResidualsChi2.cxx b/STEER/AliTrackResidualsChi2.cxx
index c3d63d8a4a82a1fde286528c3836f9832073847c..0d8c70c7b57869de202a296d0903616bc60e881d 100644 (file)
--- a/STEER/AliTrackResidualsChi2.cxx
+++ b/STEER/AliTrackResidualsChi2.cxx
@@ -26,6 +26,8 @@
 #include "AliTrackPointArray.h"
 #include "AliTrackResidualsChi2.h"
 
+
+
 void TrackResidualsChi2Fcn(Int_t &npar, Double_t *gin, Double_t &f, Double_t *x, Int_t iflag);
 
 
@@ -65,6 +67,7 @@ Bool_t AliTrackResidualsChi2::Minimize()
 
   // 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};
   ierflg = fitter->SetParameter(0, "dx", pars[0], step[0], 0,0);
   ierflg = fitter->SetParameter(1, "dy", pars[1], step[1], 0,0);
@@ -73,6 +76,14 @@ Bool_t AliTrackResidualsChi2::Minimize()
   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.;
@@ -82,8 +93,47 @@ Bool_t AliTrackResidualsChi2::Minimize()
   Double_t amin,edm,errdef;
   Int_t nvpar,nparx,icstat;
   icstat = fitter->GetStats(amin,edm,errdef,nvpar,nparx);
-  fChi2 = amin; fNdf -= nvpar;
 
+   //Construct the covariance matrix for AlignObj
+  Double_t *cov=fitter->GetCovarianceMatrix(); 
+  Int_t unfixedparam=6;
+  Int_t fixedparamat[6]={0,0,0,0,0,0};
+  if(fBFixed[0]==kTRUE){
+    unfixedparam--;
+    fixedparamat[0]=1;
+  }
+  for(Int_t j=1;j<6;j++){
+    if(fBFixed[j]==kTRUE){
+      unfixedparam--;
+      fixedparamat[j]=fixedparamat[j-1]+1;
+    }
+    else fixedparamat[j]=fixedparamat[j-1];
+  }
+  
+  Double_t cov2[36];
+  for(Int_t i=0;i<6;i++){
+    for(Int_t j=0;j<6;j++){
+      if(fBFixed[i]==kTRUE||fBFixed[j]==kTRUE){
+               cov2[i+6*j]=0.;
+      }
+      else cov2[i+6*j]=cov[i-fixedparamat[i]+6*(j-fixedparamat[j])];
+    
+    }
+  }
+  
+  Double_t covmatrarray[21];
+  for(Int_t j=0;j<6;j++){
+    for(Int_t i=j;i<6;i++){
+      covmatrarray[i*(i+1)/2+j]=cov2[i+6*j];
+    }
+  }
+  
+  //  printf("covar 2:  %.10f  \n   %.10f ; %.10f ;  \n; %.10f;  %.10f;  %.10f; \n  %.10f; %.10f;  %.10f;  %.10f;  \n  %.10f ; %.10f  ;  %.10f  ;  %.10f  ;  %.10f \n %.10f ;  %.10f ;  %.10f  ;  %.10f  ;  %.10f  ;  %.10f   ;\n",covmatrarray[0],covmatrarray[1],covmatrarray[2],covmatrarray[3],covmatrarray[4],covmatrarray[5],covmatrarray[6],covmatrarray[7],covmatrarray[8],covmatrarray[9],covmatrarray[10],covmatrarray[11],covmatrarray[12],covmatrarray[13],covmatrarray[14],covmatrarray[15],covmatrarray[16],covmatrarray[17],covmatrarray[18],covmatrarray[19],covmatrarray[20]);
+  
+  
+  fAlignObj->SetCorrMatrix(covmatrarray);
+  fChi2 = amin; fNdf -= nvpar;
+  
   return kTRUE;
 }
 
@@ -114,3 +164,5 @@ void AliTrackResidualsChi2::Chi2(Int_t & /* npar */, Double_t * /* gin */, Doubl
   }
   f = chi2;
 }
+
+