]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - ITS/ReadITSTracks.C
Correct HitsToFastRecPoints
[u/mrichter/AliRoot.git] / ITS / ReadITSTracks.C
index 1f1bcfd0ce1df79c9e994cdaea2d768473e26de3..ba7be61584da59859d1713365aca4d0ef2dd5e20 100644 (file)
 //
 //   Loop over events 
 //
-
    char tname[30];
+   for (int nev=evNumber1; nev<= evNumber2; nev++) {
 
-   for (int nev=0; nev<= evNumber2; nev++) {
+  // for (int nev=0; nev<= evNumber2; nev++) {
 
    sprintf(tname,"TreeT%d",nev);
    TTree *tracktree=(TTree*)file->Get(tname);
    TBranch *tbranch=tracktree->GetBranch("ITStracks");
+   cout<<" nev = "<<nev<<"\n";
        //cout<<" open the file \n"; 
        
    Int_t nentr=tracktree->GetEntries();
 
    TObjArray tarray(nentr);
-   AliITSiotrack *iotrack=0;
+  // 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);
+      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();           
+   //file->Close();             
        
-         AliITSiotrack *iotrack;
-   for (Int_t i=0; i<nentr; i++) {     
-        iotrack=(AliITSiotrack*)tarray.UncheckedAt(i);
+         AliITSIOTrack *iotrack;
+   for (Int_t i=0; i<nentr; i++) {
+     AliITSIOTrack *iotrack=new AliITSIOTrack;         
+        iotrack=(AliITSIOTrack*)tarray.UncheckedAt(i);
         if(!iotrack) continue;
      Int_t label=iotrack->GetLabel();
          Double_t phistate=iotrack->GetStatePhi();
          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();    
+     // cout<<" track label = "<<label<<"\n";
+     // cout<<" phi z D tanl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n";      
          
                                            
            
        printf(" i %d label x y z %d %f %f %f\n",i,label,x,y,z);
        // delete iotrack;
                 */
-                
+    delete iotrack;             
    }  
 
-   }   // event loop    
+   }   // event loop 
+   file->Close();   
 }