TPC/AliTPCtest.C

   1 /****************************************************************************
   2  *           Origin: I.Belikov, CERN, Jouri.Belikov@cern.ch                 *
   3  ****************************************************************************/
   4
   5 Int_t AliTPCtest(Int_t n = 5) {
   6    Int_t rc=0;
   7
   8 //Test TPC simulation
   9    gROOT->LoadMacro("$(ALICE_ROOT)/macros/grun.C");
  10    grun(n);
  11
  12
  13    AliKalmanTrack::SetConvConst(1000/0.299792458/gAlice->Field()->SolenoidField());
  14
  15    Int_t ver=gAlice->GetDetector("TPC")->IsVersion();
  16
  17    AliRunLoader* rl = gAlice->GetRunLoader();
  18    if (rl == 0x0)
  19     {
  20       cerr<<"Can not get run loader from gAlice"<<endl;
  21       return 1;
  22     }
  23
  24    delete rl;//close the session left after generation (grun.C)
  25    gAlice= 0x0;
  26
  27    cout<<" \n\n\nClean -> Proceeding witg digitization \n\n\n";
  28    if ((ver!=1)&&(ver!=2)) {
  29       cerr<<"Invalid TPC version: "<<ver<<" ! (must be 1 or 2)\n";
  30       return 12345;
  31    }
  32
  33    if (ver==2) {
  34      gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCHits2Digits.C");
  35      if (rc=AliTPCHits2Digits(n)) return rc;
  36
  37      //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayDigits.C");
  38      //if (rc=AliTPCDisplayDigits(1,1)) return rc;
  39    }
  40
  41
  42 //Test TPC reconstruction
  43    gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindClusters.C");
  44    if (rc=AliTPCFindClusters(n)) return rc;
  45
  46    //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayClusters.C");
  47    //if (rc=AliTPCDisplayClusters()) return rc;
  48
  49    gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindTracks.C");
  50    if (rc=AliTPCFindTracks(n)) return rc;
  51
  52    gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS");
  53    gROOT->ProcessLine(".L $(ALICE_ROOT)/TPC/AliTPCComparison.C++");
  54    if (rc=AliTPCComparison()) return rc;
  55
  56    return rc;
  57 }
  58