New version including the full 3D vertex finder
[u/mrichter/AliRoot.git] / ITS / VertexMacro.C
CommitLineData
858eb2e9 1#include "iostream.h"
6b88f180 2#include "TMath.h"
858eb2e9 3
4void VertexMacro(Int_t evNumber1=0,Int_t evNumber2=0) {
5
6 const char *filename="galice.root";
7
8 ///////////////// Dynamically link some shared libs ////////////////////////////////
9
10 if (gClassTable->GetID("AliRun") < 0) {
11 gROOT->LoadMacro("loadlibs.C");
12 loadlibs();
13 } else {
14 delete gAlice;
15 gAlice=0;
16 }
17
18// Connect the Root Galice file containing Geometry, Kine and Hits
19 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(filename);
20 if (!file) file = new TFile(filename,"UPDATE");
21
22
23// Get AliRun object from file or create it if not on file
24 if (!gAlice) {
25 gAlice = (AliRun*)file->Get("gAlice");
26 if (gAlice) printf("AliRun object found on file\n");
27 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
28 }
29
30// Loop over events
31//
32 Int_t Nh=0;
33 Int_t Nh1=0;
34 for (int nev=0; nev<= evNumber2; nev++) {
35 gAlice->SetEvent(nev);
36 Int_t nparticles = gAlice->GetEvent(nev);
37 cout << "nev " << nev <<endl;
38 cout << "nparticles " << nparticles <<endl;
39 if (nev < evNumber1) continue;
40 if (nparticles <= 0) return;
41
42
43 TStopwatch timer;
44 timer.Start();
45
6b88f180 46 AliITSVertex *V = new AliITSVertex();
858eb2e9 47
48 timer.Stop();
49 timer.Print();
50
6b88f180 51 cout << endl << "Xv = " << V->GetXv() << " cm" << endl;
52 cout << "X resolution = " << V->GetXRes()*10000 << " microns" << endl;
53 cout << "Signal/Noise for X = " << V->GetXSNR() << endl;
54 cout << endl << "Yv = " << V->GetYv() << " cm" << endl;
55 cout << "Y resolution = " << V->GetYRes()*10000 << " microns" << endl;
56 cout << "Signal/Noise for Y = " << V->GetYSNR() << endl;
858eb2e9 57 cout << endl << "Zv = " << V->GetZv() << " cm" << endl;
ba2e68c4 58 cout << "Z Resolution = " << V->GetZRes()*10000 << " microns" << endl;
59 cout << "Signal/Noise for Z = " << V->GetZSNR() <<endl;
6b88f180 60
858eb2e9 61 delete V;
858eb2e9 62 }
63
6b88f180 64 file->Close();
858eb2e9 65}
66