4 void ReadITSTracks(Int_t evNumber1=0,Int_t evNumber2=0) {
6 const char *filename="itstracks.root";
8 ///////////////// Dynamically link some shared libs ////////////////////////////////
10 if (gClassTable->GetID("AliRun") < 0) {
11 gROOT->LoadMacro("loadlibs.C");
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);
26 for (int nev=evNumber1; nev<= evNumber2; nev++) {
29 // for (int nev=0; nev<= evNumber2; nev++) {
31 sprintf(tname,"TreeT%d",nev);
32 TTree *tracktree=(TTree*)file->Get(tname);
33 TBranch *tbranch=tracktree->GetBranch("ITStracks");
34 cout<<" nev = "<<nev<<"\n";
35 //cout<<" open the file \n";
37 Int_t nentr=tracktree->GetEntries();
39 TObjArray tarray(nentr);
40 // AliITSIOTrack *iotrack=0;
41 printf("nentr %d\n",nentr);
43 for (Int_t i=0; i<nentr; i++) {
44 AliITSIOTrack *iotrack=new AliITSIOTrack;
45 // tarray.AddAt(new AliITSIOTrack,i);
46 // iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
47 tbranch->SetAddress(&iotrack);
48 tracktree->GetEvent(i);
49 tarray.AddLast(iotrack);
53 AliITSIOTrack *iotrack;
54 for (Int_t i=0; i<nentr; i++) {
55 AliITSIOTrack *iotrack=new AliITSIOTrack;
56 iotrack=(AliITSIOTrack*)tarray.UncheckedAt(i);
57 if(!iotrack) continue;
58 Int_t label=iotrack->GetLabel();
59 Double_t phistate=iotrack->GetStatePhi();
60 Double_t tgl=iotrack->GetStateTgl();
61 Double_t Zstate=iotrack->GetStateZ();
62 Double_t Dr=iotrack->GetStateD();
63 Double_t C=iotrack->GetStateC();
65 // cout<<" track label = "<<label<<"\n";
66 // cout<<" phi z D tanl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n";
72 Int_t label=iotrack->GetLabel();
73 Float_t x=iotrack->GetX();
74 Float_t y=iotrack->GetY();
75 Float_t z=iotrack->GetZ();
76 printf(" i %d label x y z %d %f %f %f\n",i,label,x,y,z);