From 961d500943ca2e3144476185f3f197186ad1e7fe Mon Sep 17 00:00:00 2001
From: barbera <barbera@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Thu, 26 Apr 2001 13:23:46 +0000
Subject: [PATCH] Test macro to read ITS tracks

---
 ITS/ReadITSTracks.C | 90 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 90 insertions(+)
 create mode 100644 ITS/ReadITSTracks.C

diff --git a/ITS/ReadITSTracks.C b/ITS/ReadITSTracks.C
new file mode 100644
index 00000000000..1f1bcfd0ce1
--- /dev/null
+++ b/ITS/ReadITSTracks.C
@@ -0,0 +1,90 @@
+#include "iostream.h"
+
+
+  void ReadITSTracks(Int_t evNumber1=0,Int_t evNumber2=0) {
+
+  const char *filename="itstracks.root";
+  
+  ///////////////// Dynamically link some shared libs ////////////////////////////////
+  
+  if (gClassTable->GetID("AliRun") < 0) {
+    gROOT->LoadMacro("loadlibs.C");
+    loadlibs();
+  } else {
+    delete gAlice;
+    gAlice=0;
+  }
+
+// Connect the Root Galice file containing Geometry, Kine and Hits
+   TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject(filename);
+   if (!file) file = new TFile(filename);
+
+//
+//   Loop over events 
+//
+
+   char tname[30];
+
+   for (int nev=0; nev<= evNumber2; nev++) {
+
+   sprintf(tname,"TreeT%d",nev);
+   TTree *tracktree=(TTree*)file->Get(tname);
+   TBranch *tbranch=tracktree->GetBranch("ITStracks");
+	//cout<<" open the file \n"; 
+	
+   Int_t nentr=tracktree->GetEntries();
+
+   TObjArray tarray(nentr);
+   AliITSiotrack *iotrack=0;
+   printf("nentr %d\n",nentr);
+	
+   for (Int_t i=0; i<nentr; i++) {
+      AliITSiotrack *iotrack=new AliITSiotrack;
+      // tarray.AddAt(new AliITSiotrack,i);
+      // iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
+       tbranch->SetAddress(&iotrack);
+       tracktree->GetEvent(i);
+		 tarray.AddLast(iotrack);
+   }
+	file->Close();		 
+	
+	  AliITSiotrack *iotrack;
+   for (Int_t i=0; i<nentr; i++) {	
+	 iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
+	 if(!iotrack) continue;
+     Int_t label=iotrack->GetLabel();
+	  Double_t phistate=iotrack->GetStatePhi();
+	  Double_t tgl=iotrack->GetStateTgl();	
+	  Double_t Zstate=iotrack->GetStateZ();
+	  Double_t Dr=iotrack->GetStateD();		  
+	  Double_t C=iotrack->GetStateC();
+	  
+	  Double_t C00,C10,C11,C20,C21,C22,C30,C31,C32,C33,C40,C41,C42,C43,C44;
+	 
+          iotrack->GetCovMatrix(C00,C10,C11,C20,C21,C22,C30,C31,C32,C33,C40,C41,C42,C43,C44);
+	  
+	  
+	  cout<<"\n"<<C00<<" "<<C10<<" "<<C11<<" "<<C20<<" "<<C21<<" "<<C22<<" "
+	  <<C30<<" "<<C31<<" "<<C32<<" "<<C33<<" "<<C40<<" "<<C41<<" "<<C42<<" "
+	  <<C43<<" "<<C44<<"\n\n\n";
+	  
+      cout<<" label = "<<label<<"\n";
+      cout<< " phi z D tgl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n";  getchar();	  
+	  
+	  		  		    
+	    
+
+ /*
+       Int_t label=iotrack->GetLabel();
+       Float_t x=iotrack->GetX();
+       Float_t y=iotrack->GetY();
+       Float_t z=iotrack->GetZ();
+       printf(" i %d label x y z %d %f %f %f\n",i,label,x,y,z);
+       // delete iotrack;
+		 */
+		 
+   }  
+
+   }   // event loop    
+}
+
-- 
2.43.0