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   gROOT->LoadMacro("$ALICE_ROOT/TPC/macros/AliXRDPROOFtoolkit.cxx+")
  19   AliXRDPROOFtoolkit tool;
  20   TChain * chain = tool.MakeChain("chain.txt","esdTree",0,10000000);
  21   chain->Lookup();
  22   mgr->SetNSysInfo(20);
  23   mgr->StartAnalysis("local",chain);
  24
  25
  26 */
  27
  28
  29 AliAnalysisManager * SetupCalibTask() {
  30   //
  31   //
  32   //
  33   TStopwatch stopwatch;
  34   stopwatch.Start();
  35
  36   AliAnalysisManager *mgr=new AliAnalysisManager("TestManager");
  37
  38   AliESDInputHandler* esdH=new AliESDInputHandler;
  39   esdH->SetActiveBranches("ESDfriend");
  40   mgr->SetInputEventHandler(esdH);
  41   //
  42   //
  43   AliCDBManager::Instance()->SetRun(1) ;
  44   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
  45   AliTPCClusterParam * clusterParam = AliTPCcalibDB::Instance()->GetClusterParam();
  46
  47   AliTPCAnalysisTaskcalib *task1=new AliTPCAnalysisTaskcalib("foo bar");
  48
  49   AliTPCcalibTracksCuts *cuts = new AliTPCcalibTracksCuts(20, 0.4, 0.5, 0.13, 0.018);
  50
  51   AliTPCcalibTracks *calibTracks =  new AliTPCcalibTracks("calibTracks", "Resolution calibration object for tracks", clusterParam, cuts);
  52   calibTracks->SetDebugLevel(5);
  53   calibTracks->SetStreamLevel(5);
  54  // ---*---*-----*-*-----*----------*---
  55   // ADD CALIB JOBS HERE!!!!!!!!!!!!!!!!
  56   task1->AddJob(new AliTPCcalibAlign);//"align","The kewl alignment job"));
  57   task1->AddJob(new AliTPCcalibTracksGain("TPCGainTracks","TPCGainTracks",cuts));
  58   task1->AddJob(calibTracks);
  59   //  task1->AddJob(new AliTPCcalibBase);
  60   // task1->AddJob(new AliTPCcalibV0);
  61   // -*----*----*---*-*------*-------**--
  62   // -------*--*---------*-----*-------*-
  63
  64   mgr->AddTask(task1);
  65
  66   AliAnalysisDataContainer *cinput1
  67     =mgr->CreateContainer("cchain1",TChain::Class(),
  68                           AliAnalysisManager::kInputContainer);
  69   AliAnalysisDataContainer *coutput1
  70     =mgr->CreateContainer("asdofhaw",TObjArray::Class(),
  71                           AliAnalysisManager::kOutputContainer,
  72                           "CalibObjects.root");
  73
  74   mgr->ConnectInput(task1,0,cinput1);
  75   mgr->ConnectOutput(task1,0,coutput1);
  76
  77   if (!mgr->InitAnalysis()) return;
  78   mgr->PrintStatus();
  79
  80   stopwatch.Stop();
  81   stopwatch.Print();
  82   return mgr;
  83 }