TPC/AliTPCFindClusters.C

   1 /****************************************************************************
   2  *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
   3  ****************************************************************************/
   4
   5 #ifndef __CINT__
   6   #include <iostream.h>
   7   #include "AliRun.h"
   8   #include "AliTPCv1.h"
   9   #include "AliTPCv2.h"
  10   #include "AliTPCParam.h"
  11
  12   #include "TFile.h"
  13   #include "TStopwatch.h"
  14 #endif
  15
  16 Int_t AliTPCFindClusters(Int_t n=1) {
  17    TFile *out=TFile::Open("AliTPCclusters.root","new");
  18    if (!out->IsOpen()) {cerr<<"Delete old AliTPCclusters.root !\n"; return 1;}
  19    TFile *in=TFile::Open("rfio:galice.root");
  20    if (!in->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
  21
  22    TFile *ind=TFile::Open("galice.root");
  23    if (!ind->IsOpen()) {cerr<<"Can't open galice.root !\n"; return 2;}
  24
  25    if (!(gAlice=(AliRun*)in->Get("gAlice"))) {
  26      cerr<<"gAlice have not been found on galice.root !\n";
  27      return 3;
  28    }
  29
  30    TDirectory *cwd = gDirectory;
  31
  32    AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC");
  33    Int_t ver = TPC->IsVersion();
  34    cerr<<"TPC version "<<ver<<" has been found !\n";
  35
  36    AliTPCParam *dig=(AliTPCParam *)in->Get("75x40_100x60_150x60");
  37    if (!dig) {cerr<<"TPC parameters have not been found !\n"; return 4;}
  38
  39    TStopwatch timer;
  40
  41    switch (ver) {
  42    case 1:
  43       cerr<<"Making clusters...\n";
  44       {
  45        AliTPCv1 &tpc=*((AliTPCv1*)TPC);
  46        tpc.SetParam(dig); timer.Start(); cwd->cd();
  47        for(Int_t i=0;i<n;i++){
  48          printf("Processing event %d\n",i);
  49          gAlice->GetEvent(i);
  50          tpc.Hits2Clusters(out,i);
  51        }
  52       }
  53       break;
  54    case 2:
  55       cerr<<"Looking for clusters...\n";
  56       {
  57         // delete gAlice; gAlice=0;
  58        AliTPCv2 tpc;
  59        tpc.SetParam(dig); timer.Start(); ind->cd();
  60        for (Int_t i=0;i<n;i++){
  61          printf("Processing event %d\n",i);
  62          tpc.Digits2Clusters(out,i);
  63          //      AliTPCclusterer::Digits2Clusters(dig, out, i);
  64        }
  65       }
  66       break;
  67    default:
  68       cerr<<"Invalid TPC version !\n";
  69       return 5;
  70    }
  71
  72    timer.Stop(); timer.Print();
  73
  74    delete gAlice; gAlice=0;
  75
  76    out->Close();
  77
  78    in->Close();
  79
  80    return 0;
  81 }