]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliRelAlignerKalman.cxx
Change the X and Y values stored by the VertexerZ from (0.,0.) which are the
[u/mrichter/AliRoot.git] / STEER / AliRelAlignerKalman.cxx
index fe59aaac8be5164d27ffed8110f3cca31bda89a2..99d35bb3b5c491edc0fc36d1290038767e8ba6d6 100644 (file)
@@ -279,7 +279,7 @@ Bool_t AliRelAlignerKalman::AddESDevent( const AliESDEvent* pEvent )
   
   Bool_t success=kFALSE;
   SetMagField( pEvent->GetMagneticField() );
-  AliESDtrack* track;
+  AliESDtrack* track=NULL;
   
   for (Int_t i=0; i<pEvent->GetNumberOfTracks(); i++)
   {
@@ -342,7 +342,7 @@ Bool_t AliRelAlignerKalman::AddCosmicEvent( const AliESDEvent* pEvent )
   TArrayI trackTArrTPC(1);
   if (!FindCosmicTrackletNumbersInEvent( trackTArrITS, trackTArrTPC, pEvent )) return kFALSE;
   SetMagField( pEvent->GetMagneticField() );
-  AliESDtrack* ptrack;
+  AliESDtrack* ptrack=NULL;
   const AliExternalTrackParam* pconstparams1;
   const AliExternalTrackParam* pconstparams2;
   AliExternalTrackParam params1;
@@ -681,9 +681,9 @@ Bool_t AliRelAlignerKalman::UpdateEstimateKalman()
 //  }
 //  else
 //  {
-    Double_t det;
-    hpht.InvertFast(&det); //since the matrix is small...
-    if (det < 2e-55) return kFALSE; //we need some sort of protection even in this case....
+    Double_t det=0.0;
+    hpht.Invert(&det); //since the matrix is small...
+    if (det < 2e-99) return kFALSE; //we need some sort of protection even in this case....
 //  }
   //printf("KalmanUpdate: det(hpht): %.4g\n",det);