Fix for FMD DA

[u/mrichter/AliRoot.git] / ITS / AliITSVertexerZTest.C

diff --git a/ITS/AliITSVertexerZTest.C b/ITS/AliITSVertexerZTest.C
index f9068e43aa32fc31a1be1490146a5986c70736a7..0f96dd052653537b294f7f787c235e1bfa354329 100644 (file)
--- a/ITS/AliITSVertexerZTest.C
+++ b/ITS/AliITSVertexerZTest.C
@@ -8,6 +8,7 @@
 #include <AliRun.h>
 #include <AliHeader.h>
 #include <AliGenEventHeader.h>
+#include <AliGeomManager.h>
 #include <AliITSVertexerZ.h>
 #include <AliRunLoader.h>
 #include <AliITSLoader.h>
@@ -41,13 +42,16 @@ void AliITSVertexerZTest(Float_t delphi=0.05,Float_t window=3.,Float_t initx=0.,
   
   AliGeomManager::LoadGeometry("geometry.root");
 
-  AliITSVertexerZ *dovert = new AliITSVertexerZ("default",initx,inity);
+  AliITSVertexerZ *dovert = new AliITSVertexerZ(initx,inity);
+  dovert->Init("default");
   //dovert->SetDebug(0);
   //  dovert->SetDiffPhiMax(delphi);
   //  dovert->SetWindow(window);
   dovert->PrintStatus();
   Int_t sigmazero=0;
   AliESDVertex *vert = 0;
+  AliITSLoader* itsloader =  (AliITSLoader*) rl->GetLoader("ITSLoader");
+  itsloader->LoadRecPoints("read");
   for(Int_t i=0; i<rl->TreeE()->GetEntries(); i++){
     rl->GetEvent(i);
     // The true Z coord. is fetched for comparison
@@ -55,22 +59,23 @@ void AliITSVertexerZTest(Float_t delphi=0.05,Float_t window=3.,Float_t initx=0.,
     AliGenEventHeader* genEventHeader = header->GenEventHeader();
     TArrayF primaryVertex(3);
     genEventHeader->PrimaryVertex(primaryVertex);
-    vert = dovert->FindVertexForCurrentEvent(i);
+    TTree* cltree = itsloader->TreeR();
+    vert = dovert->FindVertexForCurrentEvent(cltree);
     if(kDebug>0){
       // Prints the results
       cout <<"========================================================\n";
       cout << "Event number: "<<i<<")  Z Vertex:"<<endl;
       if(vert){
-       cout<<"FOUND: "<<vert->GetZv()<<"; ";
+       cout<<"FOUND: "<<vert->GetZ()<<"; ";
        cout<<vert->GetZRes()<<"; "<<vert->GetNContributors()<<endl;
        cout <<" True Z position "<<primaryVertex[2]<<", diff= ";
-       cout<<(primaryVertex[2]-vert->GetZv())*10000.<<endl;
+       cout<<(primaryVertex[2]-vert->GetZ())*10000.<<endl;
       } else {
        cout<<"NOT FOUND"<<endl;
       }
     }
     if(vert){
-      Float_t found = vert->GetZv();
+      Float_t found = vert->GetZ();
       diff2->Fill(primaryVertex[2],found);
       found = 10000.*(found-primaryVertex[2]);
       if(vert->GetZRes()!=0){