1 Int_t AliTPCFindClusters() {
2 TFile *out=TFile::Open("AliTPCclusters.root","new");
3 if (!out->IsOpen()) {cerr<<"Delete old AliTPCclusters.root !\n"; return 1;}
4 TFile *in=TFile::Open("galice.root");
5 if (!in->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
7 if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
8 cerr<<"gAlice have not been found on galice.root !\n";
12 TDirectory *cwd = gDirectory;
16 AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
17 Int_t ver = TPC->IsVersion();
18 cerr<<"TPC version "<<ver<<" has been found !\n";
20 AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60");
21 if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
27 cerr<<"Making clusters...\n";
29 AliTPCv1 &tpc=*((AliTPCv1*)TPC);
30 tpc.SetParam(dig); timer.Start(); cwd->cd(); tpc.Hits2Clusters(out);
34 cerr<<"Looking for clusters...\n";
36 delete gAlice; gAlice=0;
38 tpc.SetParam(dig); timer.Start(); cwd->cd(); tpc.Digits2Clusters(out);
42 cerr<<"Invalid TPC version !\n";
46 timer.Stop(); timer.Print();
48 delete gAlice; gAlice=0;