TPC/macros/TestV0.C

   1 /*
   2   Simple test of the V0 finder
   3   //
   4   //0. Setup memory chcecker if you want
   5   //
   6   gSystem->Load("$ROOTSYS/lib/libGui.so");
   7   gSystem->Load("$ROOTSYS/lib/libTree.so");
   8   gSystem->Load("$MEMSTAT/libMemStat.so");
   9   TMemStat *memstat = new TMemStat(100000000,10000000,kTRUE);
  10   AliSysInfo::AddCallBack(TMemStatManager::GetInstance()->fStampCallBack);
  11   AliSysInfo::AddStamp("Start");
  12   //
  13
  14   //1. Load needed libraries
  15   gSystem->Load("libANALYSIS");
  16   gSystem->Load("libTPCcalib");
  17   //
  18   // Setup analysis manager
  19   //
  20   .L $ALICE_ROOT/TPC/macros/CalibrateTPC.C
  21   AliAnalysisManager * mgr = SetupCalibTask();
  22   //
  23   // Process data - chain
  24   //
  25   gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
  26   gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+")
  27   AliXRDPROOFtoolkit tool;
  28   TChain * chain = tool.MakeChain("esd.txt","esdTree",0,50000);
  29   chain->Lookup();
  30   // memory
  31   mgr->SetNSysInfo(100);
  32   //
  33   mgr->SetDebugLevel(1);
  34   mgr->StartAnalysis("proof",chain);
  35   //mgr->StartAnalysis("local",chain);
  36   // delete manager
  37   //
  38   delete mgr;
  39   AliSysInfo::AddStamp("End");
  40   //
  41   // analyze memstat report
  42   //
  43   delete memstat;
  44   TMemStat draw("memstat.root");
  45   draw.MakeReport(0,0,"order 0 sortstat 3 sortstamp 0 sortdeep 10 stackdeep 15 maxlength 50")
  46 */
  47
  48
  49 AliAnalysisManager * SetupV0Task() {
  50   //
  51   //
  52   //
  53   TStopwatch stopwatch;
  54   stopwatch.Start();
  55   //
  56   // set magnetic field form the cosmos - it should be provided by framework
  57   AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 2);
  58   AliTracker::SetFieldMap(field,0);
  59   //
  60   AliAnalysisManager *mgr=new AliAnalysisManager("TestManager");
  61
  62   AliESDInputHandler* esdH=new AliESDInputHandler;
  63   esdH->SetActiveBranches("ESDfriend");
  64   mgr->SetInputEventHandler(esdH);
  65   //
  66   //
  67   AliCDBManager::Instance()->SetRun(1) ;
  68   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
  69
  70   AliTPCAnalysisTaskcalib *task1=new AliTPCAnalysisTaskcalib("TPC calibration task");
  71
  72   AliTPCcalibTracksCuts *cuts = new AliTPCcalibTracksCuts(20, 0.4, 0.5, 0.13, 0.018);
  73
  74   //
  75   AliTPCcalibV0 *calibV0 = new AliTPCcalibV0;
  76   calibV0->SetName("calibV0");
  77   calibV0->SetTitle("calibV0");
  78   calibV0->SetDebugLevel(20);
  79   calibV0->SetStreamLevel(2);
  80   //
  81
  82   task1->AddJob(calibV0);
  83
  84   TString path=gSystem->pwd();
  85   path+="/V0/";
  86   gSystem->mkdir(path);
  87   task1->SetDebugOuputhPath(path.Data());
  88   mgr->AddTask(task1);
  89
  90   mgr->AddTask(task1);
  91
  92   AliAnalysisDataContainer *cinput1
  93     =mgr->CreateContainer("cchain1",TChain::Class(),
  94                           AliAnalysisManager::kInputContainer);
  95   AliAnalysisDataContainer *coutput1
  96     =mgr->CreateContainer("TPCCalib",TObjArray::Class(),
  97                           AliAnalysisManager::kOutputContainer,
  98                           "CalibObjects.root");
  99
 100   mgr->ConnectInput(task1,0,cinput1);
 101   mgr->ConnectOutput(task1,0,coutput1);
 102
 103   if (!mgr->InitAnalysis()) return;
 104   mgr->PrintStatus();
 105
 106   stopwatch.Stop();
 107   stopwatch.Print();
 108   return mgr;
 109 }