6 #include "AliTPCParam.h"
9 #include "TStopwatch.h"
12 Int_t AliTPCFindClusters(Int_t n=1) {
13 TFile *out=TFile::Open("AliTPCclusters.root","new");
14 if (!out->IsOpen()) {cerr<<"Delete old AliTPCclusters.root !\n"; return 1;}
15 TFile *in=TFile::Open("rfio:galice.root");
16 if (!in->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
18 if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
19 cerr<<"gAlice have not been found on galice.root !\n";
23 TDirectory *cwd = gDirectory;
25 AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
26 Int_t ver = TPC->IsVersion();
27 cerr<<"TPC version "<<ver<<" has been found !\n";
29 AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60");
30 if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
36 cerr<<"Making clusters...\n";
38 AliTPCv1 &tpc=*((AliTPCv1*)TPC);
39 tpc.SetParam(dig); timer.Start(); cwd->cd();
40 for(Int_t i=0;i<n;i++){
41 printf("Processing event %d\n",i);
43 tpc.Hits2Clusters(out,i);
48 cerr<<"Looking for clusters...\n";
50 // delete gAlice; gAlice=0;
52 tpc.SetParam(dig); timer.Start(); cwd->cd();
53 for (Int_t i=0;i<n;i++){
54 printf("Processing event %d\n",i);
55 tpc.Digits2Clusters(out,i);
56 // AliTPCclusterer::Digits2Clusters(dig, out, i);
61 cerr<<"Invalid TPC version !\n";
65 timer.Stop(); timer.Print();
67 delete gAlice; gAlice=0;