Correct initialization of symmetric matrix. Now we have to put both upper and lower...
authorcvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 May 2006 14:17:03 +0000 (14:17 +0000)
committercvetan <cvetan@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 4 May 2006 14:17:03 +0000 (14:17 +0000)
STEER/AliTrackFitterStraight.cxx
STEER/AliTrackResidualsFast.cxx

index a92b014a6f552f68be3722c7eb8b10cd67c187a5..36f207e82976b6bbe88e132c39b71f6b6a366529 100644 (file)
@@ -224,7 +224,7 @@ void AliTrackFitterStraight::Update(){
   TMatrixDSym     smatrix(2);
   TMatrixD        sums(1,2);
   smatrix(0,0) = fSumXY[0]; smatrix(1,1)=fSumXY[2];
-  smatrix(0,1) = fSumXY[1];
+  smatrix(0,1) = fSumXY[1]; smatrix(1,0)=fSumXY[1];
   sums(0,0)    = fSumXY[3]; sums(0,1)   =fSumXY[4];
   smatrix.Invert();
   if (smatrix.IsValid()){
@@ -246,7 +246,7 @@ void AliTrackFitterStraight::Update(){
   TMatrixDSym     smatrixz(2);
   TMatrixD        sumsxz(1,2);
   smatrixz(0,0) = fSumXZ[0]; smatrixz(1,1) = fSumXZ[2];
-  smatrixz(0,1) = fSumXZ[1];
+  smatrixz(0,1) = fSumXZ[1]; smatrixz(1,0) = fSumXZ[1];
   sumsxz(0,0)   = fSumXZ[3]; sumsxz(0,1)   = fSumXZ[4];
   smatrixz.Invert();
   if (smatrixz.IsValid()){
index 7fc0c55031ba95959f5fcd4d4f0edd0ea5decd92..6d5e0c87620997ca5a02b2cf5b683e7d459b11d6 100644 (file)
@@ -184,14 +184,26 @@ Bool_t AliTrackResidualsFast::Update()
   TMatrixDSym     smatrix(6);
   TMatrixD        sums(1,6);
 
-  smatrix(0,0) = fSum[0]; smatrix(0,1) = fSum[1]; smatrix(0,2)=fSum[2];
-  smatrix(0,3) = fSum[3]; smatrix(0,4) = fSum[4]; smatrix(0,5)=fSum[5];
-  smatrix(1,1) = fSum[6]; smatrix(1,2) = fSum[7]; smatrix(1,3)=fSum[8];
-  smatrix(1,4) = fSum[9]; smatrix(1,5) = fSum[10];
-  smatrix(2,2) = fSum[11];smatrix(2,3) = fSum[12];smatrix(2,4)=fSum[13];
-  smatrix(2,5) = fSum[14];
-  smatrix(3,3) = fSum[15];smatrix(3,4) = fSum[16];smatrix(3,5)=fSum[17];
-  smatrix(4,4) = fSum[18];smatrix(4,5) = fSum[19];
+  smatrix(0,0) = fSum[0];
+  smatrix(0,1) = smatrix(1,0) = fSum[1];
+  smatrix(0,2) = smatrix(2,0) = fSum[2];
+  smatrix(0,3) = smatrix(3,0) = fSum[3];
+  smatrix(0,4) = smatrix(4,0) = fSum[4];
+  smatrix(0,5) = smatrix(5,0) = fSum[5];
+  smatrix(1,1) = fSum[6];
+  smatrix(1,2) = smatrix(2,1) = fSum[7];
+  smatrix(1,3) = smatrix(3,1) = fSum[8];
+  smatrix(1,4) = smatrix(4,1) = fSum[9];
+  smatrix(1,5) = smatrix(5,1) = fSum[10];
+  smatrix(2,2) = fSum[11];
+  smatrix(2,3) = smatrix(3,2) = fSum[12];
+  smatrix(2,4) = smatrix(4,2) = fSum[13];
+  smatrix(2,5) = smatrix(5,2) = fSum[14];
+  smatrix(3,3) = fSum[15];
+  smatrix(3,4) = smatrix(4,3) = fSum[16];
+  smatrix(3,5) = smatrix(5,3) = fSum[17];
+  smatrix(4,4) = fSum[18];
+  smatrix(4,5) = smatrix(5,4) = fSum[19];
   smatrix(5,5) = fSum[20];
 
   sums(0,0)    = fSum[21]; sums(0,1) = fSum[22]; sums(0,2) = fSum[23];