X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=ITS%2FAliITSVertexerZTest.C;h=9c70bb291025d7a090dd1feb6c1cc3583eee3f2f;hb=3a282863bcc86e188d6f4d7697c46efef2facbbd;hp=23e1d73396f71c0c1d0a1f4c0f7cf2bb99b4a7a0;hpb=98ca44ad605c3894d2b7b5ae48a9d1f5815da3ce;p=u%2Fmrichter%2FAliRoot.git diff --git a/ITS/AliITSVertexerZTest.C b/ITS/AliITSVertexerZTest.C index 23e1d73396f..9c70bb29102 100644 --- a/ITS/AliITSVertexerZTest.C +++ b/ITS/AliITSVertexerZTest.C @@ -1,5 +1,6 @@ #if !defined(__CINT__) || defined(__MAKECINT__) #include +#include #include #include #include @@ -7,7 +8,8 @@ #include #include #include -#include +#include +#include #include #include @@ -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"<GetLoader("ITSLoader"); + + AliGeomManager::LoadGeometry("geometry.root"); - if(!ITSloader){ - cerr<<"AliITSVertexerZTest.C : ITS loader not found"<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; iTreeE()->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]);