81bb131c453686c193d629732f158d9f39304fc0
[u/mrichter/AliRoot.git] / 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 }