TPC/AliTPCFindClusters.C

   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;}
   6
   7    if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
   8      cerr<<"gAlice have not been found on galice.root !\n";
   9      return 3;
  10    }
  11
  12    gAlice->GetEvent(0);
  13
  14    AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
  15    Int_t ver = TPC->IsVersion();
  16    cerr<<"TPC version "<<ver<<" has been found !\n";
  17
  18    AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60");
  19    if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
  20
  21    TStopwatch timer;
  22
  23    switch (ver) {
  24    case 1:
  25       cerr<<"Making clusters...\n";
  26       {
  27        AliTPCv1 &tpc=*((AliTPCv1*)TPC);
  28        tpc.SetParam(dig); timer.Start(); tpc.Hits2Clusters(out);
  29       }
  30       break;
  31    case 2:
  32       cerr<<"Looking for clusters...\n";
  33       {
  34        delete gAlice; gAlice=0;
  35        AliTPCv2 tpc;
  36        tpc.SetParam(dig); timer.Start(); tpc.Digits2Clusters(out);
  37       }
  38       break;
  39    default:
  40       cerr<<"Invalid TPC version !\n";
  41       return 5;
  42    }
  43
  44    timer.Stop(); timer.Print();
  45
  46    delete gAlice; gAlice=0;
  47
  48    out->Close();
  49    in->Close();
  50
  51    return 0;
  52 }