First version of the parallel TPC tracking (M.Ivanov)

[u/mrichter/AliRoot.git] / TPC / AliTPCFindClusters.C

diff --git a/TPC/AliTPCFindClusters.C b/TPC/AliTPCFindClusters.C
index cb2d4bf17e0c2449e02c244ef6642ea8f464955a..bf33e302583b7d3c488dc9d16d8f98619f16422c 100644 (file)
--- a/TPC/AliTPCFindClusters.C
+++ b/TPC/AliTPCFindClusters.C
@@ -1,3 +1,7 @@
+/****************************************************************************
+ *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
+ ****************************************************************************/
+
 #ifndef __CINT__
   #include <iostream.h>
   #include "AliRun.h"
@@ -15,6 +19,9 @@ Int_t AliTPCFindClusters(Int_t n=1) {
    TFile *in=TFile::Open("rfio:galice.root");
    if (!in->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
 
+   TFile *ind=TFile::Open("digits.root");
+   if (!ind->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
+
    if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
      cerr<<"gAlice have not been found on galice.root !\n";
      return 3;
@@ -26,16 +33,7 @@ Int_t AliTPCFindClusters(Int_t n=1) {
    Int_t ver = TPC->IsVersion(); 
    cerr<<"TPC version "<<ver<<" has been found !\n";
 
-   AliTPCParamSR *dig=(AliTPCParamSR *)in->Get("75x40_100x60");
-   if(dig){
-     cerr<<"2 pad-length geom hits with 3 pad-lengths geom digits\n";
-     delete dig;
-     dig = new AliTPCParamSR();
-   }
-   else
-   {
-     dig=(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60");
-   }
+   AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60_150x60");
    if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
 
    TStopwatch timer;
@@ -58,7 +56,7 @@ Int_t AliTPCFindClusters(Int_t n=1) {
       {
        // delete gAlice; gAlice=0;
        AliTPCv2 tpc; 
-       tpc.SetParam(dig); timer.Start(); cwd->cd();  
+       tpc.SetParam(dig); timer.Start(); ind->cd();  
        for (Int_t i=0;i<n;i++){
         printf("Processing event %d\n",i);
          tpc.Digits2Clusters(out,i);