TPC/macros/CalibrateTPC.C

   1 /*
   2   Simple calibration analysis
   3
   4
   5   //1. Load needed libraries
   6   gSystem->Load("libANALYSIS");
   7   gSystem->Load("libTPCcalib");
   8   //
   9   // Setup analysis manager
  10   //
  11   gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+")
  12   .L $ALICE_ROOT/TPC/macros/CalibrateTPC.C
  13   AliAnalysisManager * mgr = SetupCalibTask();
  14   //
  15   // Process data - chain
  16   //
  17   gSystem->AddIncludePath("-I$ALICE_ROOT/TPC/macros");
  18   AliXRDPROOFtoolkit tool;
  19   TChain * chain = tool.MakeChain("cahin.txt","esdTree",0,10)
  20   mgr->StartAnalysis("local",chain);
  21
  22
  23 */
  24
  25
  26 AliAnalysisManager * SetupCalibTask() {
  27   //
  28   //
  29   //
  30   TStopwatch stopwatch;
  31   stopwatch.Start();
  32
  33   AliAnalysisManager *mgr=new AliAnalysisManager("TestManager");
  34
  35   AliESDInputHandler* esdH=new AliESDInputHandler;
  36   esdH->SetActiveBranches("ESDfriend");
  37   mgr->SetInputEventHandler(esdH);
  38
  39   AliTPCAnalysisTaskcalib *task1=new AliTPCAnalysisTaskcalib("foo bar");
  40
  41   AliTPCcalibTracksCuts *cuts = new AliTPCcalibTracksCuts(20, 0.4, 0.5, 0.13, 0.018);
  42
  43   // ---*---*-----*-*-----*----------*---
  44   // ADD CALIB JOBS HERE!!!!!!!!!!!!!!!!
  45   task1->AddJob(new AliTPCcalibAlign);//"align","The kewl alignment job"));
  46   //  task1->AddJob(new AliTPCcalibTracks("resolution","I would have been called AliTPCcalibResolution in a bit more perfect world.",0,cuts));
  47   task1->AddJob(new AliTPCcalibTracksGain("resolution","I would have been called AliTPCcalibGain in a bit more perfect world.",cuts));
  48   //  task1->AddJob(new AliTPCcalibBase);
  49   // task1->AddJob(new AliTPCcalibV0);
  50   // -*----*----*---*-*------*-------**--
  51   // -------*--*---------*-----*-------*-
  52
  53   mgr->AddTask(task1);
  54
  55   AliAnalysisDataContainer *cinput1
  56     =mgr->CreateContainer("cchain1",TChain::Class(),
  57                           AliAnalysisManager::kInputContainer);
  58   AliAnalysisDataContainer *coutput1
  59     =mgr->CreateContainer("asdofhaw",TObjArray::Class(),
  60                           AliAnalysisManager::kOutputContainer,
  61                           "CalibObjects.root");
  62
  63   mgr->ConnectInput(task1,0,cinput1);
  64   mgr->ConnectOutput(task1,0,coutput1);
  65
  66   if (!mgr->InitAnalysis()) return;
  67   mgr->PrintStatus();
  68
  69   stopwatch.Stop();
  70   stopwatch.Print();
  71   return mgr;
  72 }