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       cerr<<"Can not get run loader from gAlice"<<endl;
  20       return 1;
  21    }
  22
  23    delete rl;//close the session left after generation (grun.C)
  24    delete gAlice; gAlice=0x0;
  25
  26    cout<<" \n\n\nClean -> Proceeding witg digitization \n\n\n";
  27    if ((ver!=1)&&(ver!=2)) {
  28       cerr<<"Invalid TPC version: "<<ver<<" ! (must be 1 or 2)\n";
  29       return 12345;
  30    }
  31
  32    if (ver==2) {
  33      gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCHits2Digits.C");
  34      if (rc=AliTPCHits2Digits(n)) return rc;
  35
  36      //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayDigits.C");
  37      //if (rc=AliTPCDisplayDigits(1,1)) return rc;
  38    }
  39
  40
  41 //Test TPC reconstruction
  42    gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindClusters.C");
  43    if (rc=AliTPCFindClusters(n)) return rc;
  44
  45    //gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCDisplayClusters.C");
  46    //if (rc=AliTPCDisplayClusters()) return rc;
  47
  48    gROOT->LoadMacro("$(ALICE_ROOT)/TPC/AliTPCFindTracks.C");
  49    if (rc=AliTPCFindTracks(n)) return rc;
  50
  51    gSystem->SetIncludePath("-I$ROOTSYS/include -I$ALICE_ROOT/include -I$ALICE_ROOT/TPC -I$ALICE_ROOT/CONTAINERS");
  52    gROOT->ProcessLine(".L $(ALICE_ROOT)/TPC/AliTPCComparison.C++");
  53    if (rc=AliTPCComparison()) return rc;
  54
  55    return rc;
  56 }
  57