From: marian <marian@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Sat, 27 Jun 2009 12:00:22 +0000 (+0000)
Subject: Test of the Propagation and Update
X-Git-Url: http://git.uio.no/git/?a=commitdiff_plain;h=c45d414a54a37b000c7c2c3ea4e6798e1d2309ee;p=u%2Fmrichter%2FAliRoot.git

Test of the Propagation and Update
Change of the mag filed orientation.

(Marian)
---

diff --git a/PWG1/AliMCTrackingTestTask.cxx b/PWG1/AliMCTrackingTestTask.cxx
index a648bdd880e..a2a27a2ea22 100644
--- a/PWG1/AliMCTrackingTestTask.cxx
+++ b/PWG1/AliMCTrackingTestTask.cxx
@@ -441,6 +441,8 @@ void  AliMCTrackingTestTask::FitTrackRefs(TParticle * part, TClonesArray * trefs
   Float_t alphaOut= TMath::ATan2(refOut->Y(),refOut->X());
   Float_t radiusOut= refOut->R();
 */
+  Bool_t isOKP=kTRUE;
+  Bool_t isOKU=kTRUE;
   AliMagF * field = (AliMagF*) TGeoGlobalMagField::Instance()->GetField();
   for (Int_t iref = iref0; iref<=iref1; iref++){
     AliTrackReference * ref = (AliTrackReference*)trefs->At(iref);
@@ -448,17 +450,17 @@ void  AliMCTrackingTestTask::FitTrackRefs(TParticle * part, TClonesArray * trefs
     Double_t pos[3] = {ref->X(), ref->Y(), ref->Z()};
     Double_t mag[3];
     field->Field(pos,mag);
-    paramPropagate->Rotate(alphaC);
-    paramUpdate->Rotate(alphaC);
+    isOKP&=paramPropagate->Rotate(alphaC);
+    isOKU&=paramUpdate->Rotate(alphaC);
     for (Float_t xref= paramPropagate->GetX(); xref<ref->R(); xref++){
-      paramPropagate->PropagateTo(xref, -mag[2]);
-      paramUpdate->PropagateTo(xref, -mag[2]);
+      isOKP&=paramPropagate->PropagateTo(xref, mag[2]);
+      isOKU&=paramUpdate->PropagateTo(xref, mag[2]);
     }
-    paramPropagate->PropagateTo(ref->R(), -mag[2]);
-    paramUpdate->PropagateTo(ref->R(), -mag[2]);
+    isOKP&=paramPropagate->PropagateTo(ref->R(), mag[2]);
+    isOKU&=paramUpdate->PropagateTo(ref->R(), mag[2]);
     Double_t clpos[2] = {0, ref->Z()};
     Double_t clcov[3] = { 0.005,0,0.005};
-    paramUpdate->Update(clpos, clcov);  
+    isOKU&= paramUpdate->Update(clpos, clcov);  
   }
   TTreeSRedirector *pcstream = GetDebugStreamer();
   if (pcstream){
@@ -472,6 +474,8 @@ void  AliMCTrackingTestTask::FitTrackRefs(TParticle * part, TClonesArray * trefs
     paramPropagate->GetPxPyPz(gmomP.GetMatrixArray());
 
      (*pcstream)<<"MCupdate"<<
+       "isOKU="<<isOKU<<
+       "isOKP="<<isOKP<<
        "m="<<mass<<
        "q="<<charge<<
        "part.="<<part<<