Updated version (Yu.Belikov)

diff --git a/TPC/AliTPCComparison.C b/TPC/AliTPCComparison.C
index 2e34701611ac2bdd9b58986dfd633e834d91397c..632b5752ebde285b4d966b9b59bfe21ec3c60773 100644 (file)
--- a/TPC/AliTPCComparison.C
+++ b/TPC/AliTPCComparison.C
@@ -54,7 +54,7 @@ Int_t AliTPCComparison(Bool_t thcut=1.0) {
    }
 
    rl->LoadgAlice();
    }
 
    rl->LoadgAlice();

-   if (rl->GetAliRun()) 
+   if (rl->GetAliRun())

    AliKalmanTrack::SetConvConst(
      1000/0.299792458/rl->GetAliRun()->Field()->SolenoidField()
    );
    AliKalmanTrack::SetConvConst(
      1000/0.299792458/rl->GetAliRun()->Field()->SolenoidField()
    );


@@ -62,7 +62,7 @@ Int_t AliTPCComparison(Bool_t thcut=1.0) {
        cerr<<"AliTPCComparison.C :Can't get AliRun !\n";
        return 1;
    }
        cerr<<"AliTPCComparison.C :Can't get AliRun !\n";
        return 1;
    }

-   rl->UnloadgAlice();
+   //rl->UnloadgAlice();

   
    AliTPCLoader * tpcl = (AliTPCLoader *)rl->GetLoader("TPCLoader");
    if (tpcl == 0x0) {
   
    AliTPCLoader * tpcl = (AliTPCLoader *)rl->GetLoader("TPCLoader");
    if (tpcl == 0x0) {


@@ -71,26 +71,6 @@ Int_t AliTPCComparison(Bool_t thcut=1.0) {
        return 3;
    }
 
        return 3;
    }
 

-   Int_t nentr=0,i=0; TObjArray tarray(MAX);
-   { /*Load tracks*/
-   
-     tpcl->LoadTracks();
-     
-     TTree *tracktree=tpcl->TreeT();
-     if (!tracktree) {cerr<<"Can't get a tree with TPC tracks !\n"; return 4;}
-
-     TBranch *tbranch=tracktree->GetBranch("tracks");
-     nentr=(Int_t)tracktree->GetEntries();
-     AliTPCtrack *iotrack=0;
-     for (i=0; i<nentr; i++) {
-       iotrack=new AliTPCtrack;
-       tbranch->SetAddress(&iotrack);
-       tracktree->GetEvent(i);
-       tarray.AddLast(iotrack);
-     }   
-     tpcl->UnloadTracks();
-   }
-

    /* Generate a list of "good" tracks */
 
    GoodTrackTPC gt[MAX];
    /* Generate a list of "good" tracks */
 
    GoodTrackTPC gt[MAX];


@@ -133,6 +113,25 @@ Int_t AliTPCComparison(Bool_t thcut=1.0) {
       out.close();
    }
 
       out.close();
    }
 

+   Int_t nentr=0,i=0; TObjArray tarray(MAX);
+   { /*Load tracks*/
+   
+     tpcl->LoadTracks();
+     
+     TTree *tracktree=tpcl->TreeT();
+     if (!tracktree) {cerr<<"Can't get a tree with TPC tracks !\n"; return 4;}
+
+     TBranch *tbranch=tracktree->GetBranch("tracks");
+     nentr=(Int_t)tracktree->GetEntries();
+     AliTPCtrack *iotrack=0;
+     for (i=0; i<nentr; i++) {
+       iotrack=new AliTPCtrack;
+       tbranch->SetAddress(&iotrack);
+       tracktree->GetEvent(i);
+       tarray.AddLast(iotrack);
+     }   
+     tpcl->UnloadTracks();
+   }

 
    TH1F *hp=new TH1F("hp","PHI resolution",50,-20.,20.); hp->SetFillColor(4);
    TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-20,20);hl->SetFillColor(4);
 
    TH1F *hp=new TH1F("hp","PHI resolution",50,-20.,20.); hp->SetFillColor(4);
    TH1F *hl=new TH1F("hl","LAMBDA resolution",50,-20,20);hl->SetFillColor(4);