First commit.

[u/mrichter/AliRoot.git] / ITS / AliITSComparisonV1.C

diff --git a/ITS/AliITSComparisonV1.C b/ITS/AliITSComparisonV1.C
index 9882e6e01f8124b60c4a4509e417d7afe48bb605..e7166c4b3afc9abea45c6d5c3fbdf482cc1ed9a1 100644 (file)
--- a/ITS/AliITSComparisonV1.C
+++ b/ITS/AliITSComparisonV1.C
@@ -34,7 +34,7 @@
     Bool_t flag;
   };   
     Int_t i;
-  GoodTrack gt[15000]; 
+  GoodTrack gt[50000]; 
   Int_t ngood=0;  
   cerr<<"Reading itsgood tracks...\n";
   while (in>>gt[ngood].fEventN>>gt[ngood].lab>>gt[ngood].code
@@ -44,7 +44,7 @@
          >>gt[ngood].ptg >>gt[ngood].flag) {
     ngood++;
     cerr<<ngood<<'\r';
-    if (ngood==15000) {
+    if (ngood==50000) {
       cerr<<"Too many good tracks !\n";
       break;
     }
@@ -122,8 +122,7 @@
         iotrack=(AliITSIOTrack*)tarray.UncheckedAt(i);
         if(!iotrack) continue;
      Int_t labITS=iotrack->GetLabel();
-     Int_t labTPC=iotrack->GetTPCLabel();  
-     //Int_t labTPC=labITS;  //provvisoria  
+     Int_t labTPC=iotrack->GetTPCLabel();   
          Double_t phistate=iotrack->GetStatePhi();
          Double_t tgl=iotrack->GetStateTgl();  
          Double_t Zstate=iotrack->GetStateZ();
@@ -136,9 +135,9 @@
          Int_t c = iotrack->GetCharge(); 
          Double_t x=iotrack->GetX();
          Double_t y=iotrack->GetY();
-         Double_t z= iotrack->GetZ(); 
-         Double_t Dz=z;   //non e' vero bisogna levare vertice
-         Double_t Dtot= TMath::Sqrt(Dr*Dr+Dz*Dz);
+         Double_t z= iotrack->GetZ();
+         Double_t Dz=iotrack->GetDz(); 
+
          
      // cout<<" track label = "<<label<<"\n";
      // cout<<" phi z D tanl C = "<<phistate<<" "<<Zstate<<" "<<Dr<<" "<<tgl<<" "<<C<<"\n";      
@@ -146,7 +145,8 @@
   //  Int_t ilab=TMath::Abs(iotrack->GetLabel());
     Int_t flaglab=0;   
     Int_t iii=0;
-   Double_t ptg=0.,pxg=0.,pyg=0.,pzg=0.; 
+   Double_t ptg=0.,pxg=0.,pyg=0.,pzg=0.;
+       Double_t xo=0., yo=0., zo=0.;    
 
     Int_t mingood=0, maxgood=0, jj=0;
     if(nev==evNumber1) mingood=0;
@@ -163,7 +163,10 @@
        ptg=gt[iii].ptg; 
        pxg=gt[iii].pxg;
        pyg=gt[iii].pyg;
-       pzg=gt[iii].pzg;        
+       pzg=gt[iii].pzg;
+       xo=gt[iii].x;
+       yo=gt[iii].y;
+       zo=gt[iii].z;                   
        break;
       }
     }   
@@ -191,7 +194,7 @@
      if(phi<0.) phi+=duepi;      
       Double_t signC=0.; 
       if(c>0) signC=1.; else signC=-1.;
-        
+         Double_t Dtot= TMath::Sqrt(Dr*Dr+Dz*Dz);       
       Double_t difphi = (phi - phig)*1000.;
       dataOut(kkk)=difphi; kkk++;
       dataOut(kkk)=Dtot*1.e4; kkk++;