Add option for building ideal (no dead channels) calibration objects

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

diff --git a/ITS/AliITSVertexerZTest.C b/ITS/AliITSVertexerZTest.C
index 23e1d73396f71c0c1d0a1f4c0f7cf2bb99b4a7a0..9c70bb291025d7a090dd1feb6c1cc3583eee3f2f 100644 (file)
--- a/ITS/AliITSVertexerZTest.C
+++ b/ITS/AliITSVertexerZTest.C
@@ -1,5 +1,6 @@
 #if !defined(__CINT__) || defined(__MAKECINT__)
 #include <TFile.h>
+#include <TGeoManager.h>
 #include <TH1F.h>
 #include <TH2F.h>
 #include <TTree.h>
@@ -7,7 +8,8 @@
 #include <AliRun.h>
 #include <AliHeader.h>
 #include <AliGenEventHeader.h>
-#include <AliITSVertexerPPZ.h>
+#include <AliGeomManager.h>
+#include <AliITSVertexerZ.h>
 #include <AliRunLoader.h>
 #include <AliITSLoader.h>
 
@@ -37,22 +39,19 @@ void AliITSVertexerZTest(Float_t delphi=0.05,Float_t window=3.,Float_t initx=0.,
     cerr<<"AliITSVertexerZTest.C : LoadKinematics returned error"<<endl;
     return;
   }
-  AliITSLoader* ITSloader =  (AliITSLoader*) rl->GetLoader("ITSLoader");
+  
+  AliGeomManager::LoadGeometry("geometry.root");
 
-  if(!ITSloader){
-    cerr<<"AliITSVertexerZTest.C :  ITS loader not found"<<endl;
-    return;
-  }
-  //  ITSloader->LoadRecPoints("read");
-  //  TFile *fo = new TFile("vertici.root","recreate");
-  //  AliITSVertexerPPZ *dovert = new AliITSVertexerPPZ("default",initx,inity);
-  AliITSVertexerZ *dovert = new AliITSVertexerZ("default",initx,inity);
-  dovert->SetDebug(0);
+  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
@@ -60,7 +59,8 @@ 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";
@@ -75,9 +75,6 @@ void AliITSVertexerZTest(Float_t delphi=0.05,Float_t window=3.,Float_t initx=0.,
       }
     }
     if(vert){
-      Double_t pos[3];
-      for(Int_t kk=0;kk<3;kk++)pos[kk]=(Double_t)primaryVertex[kk];
-      vert->SetTruePos(pos);
       Float_t found = vert->GetZv();
       diff2->Fill(primaryVertex[2],found);
       found = 10000.*(found-primaryVertex[2]);