+
+void AliExternalTrackParam::CheckCovariance() {
+
+ // This function forces the diagonal elements of the covariance matrix to be positive.
+ // In case the diagonal element is bigger than the maximal allowed value, it is set to
+ // the limit and the off-diagonal elements that correspond to it are set to zero.
+
+ fC[0] = TMath::Abs(fC[0]);
+ if (fC[0]>kC0max) {
+ fC[0] = kC0max;
+ fC[1] = 0;
+ fC[3] = 0;
+ fC[6] = 0;
+ fC[10] = 0;
+ }
+ fC[2] = TMath::Abs(fC[2]);
+ if (fC[2]>kC2max) {
+ fC[2] = kC2max;
+ fC[1] = 0;
+ fC[4] = 0;
+ fC[7] = 0;
+ fC[11] = 0;
+ }
+ fC[5] = TMath::Abs(fC[5]);
+ if (fC[5]>kC5max) {
+ fC[5] = kC5max;
+ fC[3] = 0;
+ fC[4] = 0;
+ fC[8] = 0;
+ fC[12] = 0;
+ }
+ fC[9] = TMath::Abs(fC[9]);
+ if (fC[9]>kC9max) {
+ fC[9] = kC9max;
+ fC[6] = 0;
+ fC[7] = 0;
+ fC[8] = 0;
+ fC[13] = 0;
+ }
+ fC[14] = TMath::Abs(fC[14]);
+ if (fC[14]>kC14max) {
+ fC[14] = kC14max;
+ fC[10] = 0;
+ fC[11] = 0;
+ fC[12] = 0;
+ fC[13] = 0;
+ }
+
+ // The part below is used for tests and normally is commented out
+// TMatrixDSym m(5);
+// TVectorD eig(5);
+
+// m(0,0)=fC[0];
+// m(1,0)=fC[1]; m(1,1)=fC[2];
+// m(2,0)=fC[3]; m(2,1)=fC[4]; m(2,2)=fC[5];
+// m(3,0)=fC[6]; m(3,1)=fC[7]; m(3,2)=fC[8]; m(3,3)=fC[9];
+// m(4,0)=fC[10]; m(4,1)=fC[11]; m(4,2)=fC[12]; m(4,3)=fC[13]; m(4,4)=fC[14];
+
+// m(0,1)=m(1,0);
+// m(0,2)=m(2,0); m(1,2)=m(2,1);
+// m(0,3)=m(3,0); m(1,3)=m(3,1); m(2,3)=m(3,2);
+// m(0,4)=m(4,0); m(1,4)=m(4,1); m(2,4)=m(4,2); m(3,4)=m(4,3);
+// m.EigenVectors(eig);
+
+// // assert(eig(0)>=0 && eig(1)>=0 && eig(2)>=0 && eig(3)>=0 && eig(4)>=0);
+// if (!(eig(0)>=0 && eig(1)>=0 && eig(2)>=0 && eig(3)>=0 && eig(4)>=0)) {
+// AliWarning("Negative eigenvalues of the covariance matrix!");
+// this->Print();
+// eig.Print();
+// }
+}