+ tpc.SetParam(dig); timer.Start(); cwd->cd();
+
+ n = rl->GetNumberOfEvents();
+ for (Int_t i=0;i<n;i++)
+ {
+ rl->GetEvent(i);
+ AliTPCclusterer clusterer(dig);
+
+ TTree * input = tpcl->TreeD();
+ if (input == 0x0)
+ {
+ cerr << "Can not get TreeD for event " << i <<endl;
+ continue;
+ }
+
+ TTree * output = tpcl->TreeR();
+ if (output == 0x0)
+ {
+ tpcl->MakeTree("R");
+ output = tpcl->TreeR();
+ if (output == 0x0)
+ {
+ cerr << "Problems with output tree (TreeR) for event " << i <<endl;
+ continue;
+ }
+ }
+
+ printf("Processing event %d\n",i);
+ clusterer.SetInput(input);
+ clusterer.SetOutput(output);
+ clusterer.Digits2Clusters();
+
+ tpcl->WriteRecPoints("OVERWRITE");
+ }
+ }
+ break;