The present commit corresponds to an important change in the way the
[u/mrichter/AliRoot.git] / 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   AliAnalysisManager *mgr=new AliAnalysisManager("TestManager");
57
58   AliESDInputHandler* esdH=new AliESDInputHandler;
59   esdH->SetActiveBranches("ESDfriend");
60   mgr->SetInputEventHandler(esdH);  
61   //
62   //
63   AliCDBManager::Instance()->SetRun(1) ;
64   AliCDBManager::Instance()->SetDefaultStorage("local://$ALICE_ROOT");
65
66   AliTPCAnalysisTaskcalib *task1=new AliTPCAnalysisTaskcalib("TPC calibration task");
67   
68   AliTPCcalibTracksCuts *cuts = new AliTPCcalibTracksCuts(20, 0.4, 0.5, 0.13, 0.018);
69
70   //
71   AliTPCcalibV0 *calibV0 = new AliTPCcalibV0;
72   calibV0->SetName("calibV0");
73   calibV0->SetTitle("calibV0");
74   calibV0->SetDebugLevel(20);
75   calibV0->SetStreamLevel(2);
76   //
77   
78   task1->AddJob(calibV0);
79  
80   TString path=gSystem->pwd();
81   path+="/V0/";
82   gSystem->mkdir(path);
83   task1->SetDebugOuputhPath(path.Data());
84   mgr->AddTask(task1);
85
86   mgr->AddTask(task1);
87
88   AliAnalysisDataContainer *cinput1
89     =mgr->CreateContainer("cchain1",TChain::Class(),
90                           AliAnalysisManager::kInputContainer);
91   AliAnalysisDataContainer *coutput1
92     =mgr->CreateContainer("TPCCalib",TObjArray::Class(),
93                           AliAnalysisManager::kOutputContainer,
94                           "CalibObjects.root");
95
96   mgr->ConnectInput(task1,0,cinput1);
97   mgr->ConnectOutput(task1,0,coutput1);
98
99   if (!mgr->InitAnalysis()) return;
100   mgr->PrintStatus(); 
101   
102   stopwatch.Stop();
103   stopwatch.Print();
104   return mgr;
105 }