1 #if !defined(__CINT__) || defined(__MAKECINT__)
3 #include <TGeoManager.h>
10 #include <AliGenEventHeader.h>
11 #include <AliITSVertexerZ.h>
12 #include <AliRunLoader.h>
13 #include <AliITSLoader.h>
17 void AliITSVertexerZTest(Float_t delphi=0.05,Float_t window=3.,Float_t initx=0., Float_t inity=0.,TString FileIn="galice.root"){
18 // delphi ---> azimuthal range to accept tracklets
19 // window ---> window in Z around the peak of tracklets proj. in mm
21 TH2F *diff2 = new TH2F("diff2","Zfound vs Ztrue",100,-6,6,100,-6,6);
22 TH1F *diff1 = new TH1F("diff1","Zfound - Ztrue(#mu m)",100,-500,500);
25 AliRunLoader *rl = AliRunLoader::Open(FileIn.Data());
26 Int_t retval = rl->LoadgAlice();
28 cerr<<"AliITSVertexerZTest.C: AliRun object not found"<<endl;
31 retval = rl->LoadHeader();
33 cerr<<"AliITSVertexerZTest.C : LoadHeader returned error"<<endl;
36 retval = rl->LoadKinematics();
38 cerr<<"AliITSVertexerZTest.C : LoadKinematics returned error"<<endl;
42 AliGeomManager::LoadGeometry("geometry.root");
44 AliITSVertexerZ *dovert = new AliITSVertexerZ("default",initx,inity);
45 //dovert->SetDebug(0);
46 // dovert->SetDiffPhiMax(delphi);
47 // dovert->SetWindow(window);
48 dovert->PrintStatus();
50 AliESDVertex *vert = 0;
51 for(Int_t i=0; i<rl->TreeE()->GetEntries(); i++){
53 // The true Z coord. is fetched for comparison
54 AliHeader *header = rl->GetHeader();
55 AliGenEventHeader* genEventHeader = header->GenEventHeader();
56 TArrayF primaryVertex(3);
57 genEventHeader->PrimaryVertex(primaryVertex);
58 vert = dovert->FindVertexForCurrentEvent(i);
61 cout <<"========================================================\n";
62 cout << "Event number: "<<i<<") Z Vertex:"<<endl;
64 cout<<"FOUND: "<<vert->GetZv()<<"; ";
65 cout<<vert->GetZRes()<<"; "<<vert->GetNContributors()<<endl;
66 cout <<" True Z position "<<primaryVertex[2]<<", diff= ";
67 cout<<(primaryVertex[2]-vert->GetZv())*10000.<<endl;
69 cout<<"NOT FOUND"<<endl;
73 Float_t found = vert->GetZv();
74 diff2->Fill(primaryVertex[2],found);
75 found = 10000.*(found-primaryVertex[2]);
76 if(vert->GetZRes()!=0){
81 dovert->WriteCurrentVertex();
85 cout<<"Only one tracklet (sigma = 0) "<<sigmazero<<endl;