-
- for (Int_t i=0;i<n;i++){
- AliTPCclustererMI clusterer;
- char dname[100];
- char cname[100];
- sprintf(dname,"TreeD_75x40_100x60_150x60_%d",i);
- sprintf(cname,"TreeC_TPC_%d",i);
- TTree * input = (TTree*)ind->Get(dname);
- out->cd();
- TTree * output = new TTree(cname,cname);
-
- printf("Processing event %d\n",i);
- clusterer.SetInput(input);
- clusterer.SetOutput(output);
- clusterer.Digits2Clusters(dig, i);
- //tpc.Digits2Clusters(out,i);
- // AliTPCclusterer::Digits2Clusters(dig, out, i);
+ n = rl->GetNumberOfEvents();
+ for (Int_t i=0;i<n;i++)
+ {
+ rl->GetEvent(i);
+ AliTPCclustererMI 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");