template to access pmd pbpb data: sjena
[u/mrichter/AliRoot.git] / PMD / anal / runPmdTask.C
1 //=========================================================================//
2 //                                                                         //
3 //                   A template runPmdTask for PMD analysis                //
4 //            You can copy it and add features according to your need      //
5 //                                                                         //
6 //                               Satyajit Jena                             //
7 //                               sjena@cern.ch                             //
8 //                                13/04/2012                               //
9 //                                                                         //
10 //=========================================================================//
11
12 const char * incollection = "file.txt";
13 //______________________________________________________________________________
14 void runPmdTask(Bool_t         isGrid = 0,
15                 Bool_t           isMC = 0,
16                 const char *gridmode  = "test") {
17
18  
19   if (isGrid) {
20     Printf("Strating the Grid Job ");
21     Printf("Grid Mode %s",gridmode);
22     if (isMC)  Printf("It is Data Type Run");
23     if (!isMC) Printf("It is MC Type Run");
24   }
25     
26   // Load the needed libraries most of 
27   // them already loaded by aliroot
28   //--------------------------------------
29   gSystem->Load("libTree");
30   gSystem->Load("libGeom");
31   gSystem->Load("libVMC");
32   gSystem->Load("libXMLIO");
33   gSystem->Load("libPhysics");
34
35   gSystem->Load("libSTEERBase");
36   gSystem->Load("libESD");
37   gSystem->Load("libAOD");
38   gSystem->Load("libANALYSIS");
39   gSystem->Load("libANALYSISalice");
40
41   // Use AliRoot includes to compile our task                                   
42   gROOT->ProcessLine(".include $ALICE_INSTALL/include");
43   
44   AliAnalysisManager *mgr = new AliAnalysisManager("PMDAnalysis");
45   
46   if(isGrid) {
47     gROOT->LoadMacro("CreateAlienHandler.C");
48     AliAnalysisGrid *alienHandler = CreateAlienHandler(gridmode);  
49     if (!alienHandler) return;
50     mgr->SetGridHandler(alienHandler);
51   }
52     
53   
54  
55   //  TChain *  chain;
56   //  if(!isGrid) {
57   //  gROOT->LoadMacro("CreateESDChain.C"); // use it if you know it
58   //  chain = CreateESDChain("file.txt", 10);
59   //  }
60  
61   TChain *  chain = new TChain("esdTree"); 
62   if(!isGrid) {
63     ifstream file_collect(incollection);
64     TString line;
65     while (line.ReadLine(file_collect) ) {
66       chain->Add(line.Data());
67     }
68   }
69     
70   if(isMC) {
71     AliVEventHandler* esdH = new AliESDInputHandler();
72     mgr->SetInputEventHandler(esdH);
73     AliMCEventHandler *mc = new AliMCEventHandler();
74     //mc->SetReadTR(kTRUE);
75     mgr->SetMCtruthEventHandler(mc);
76     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
77     AliPhysicsSelectionTask* physicsSelTask = AddTaskPhysicsSelection(isMC);
78     physicsSelTask->GetPhysicsSelection()->SetAnalyzeMC();
79
80     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
81     AliCentralitySelectionTask *centralityTask = AddTaskCentrality();
82   } else {
83     AliVEventHandler* esdH = new AliESDInputHandler();
84     mgr->SetInputEventHandler(esdH);
85     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskPhysicsSelection.C");
86     AliPhysicsSelectionTask* physicsSelTask = AddTaskPhysicsSelection(isMC);
87     
88     gROOT->LoadMacro("$ALICE_ROOT/ANALYSIS/macros/AddTaskCentrality.C");
89     AliCentralitySelectionTask *centralityTask = AddTaskCentrality();
90   }
91
92   gROOT->LoadMacro("AliPMDAnalysisTaskPbPb.cxx++g");   
93   gROOT->LoadMacro("AddAliPMDAnalysisTaskPbPb.C");
94   AddAliPMDAnalysisTaskPbPb("MyTask",isMC);
95   
96   // Enable debug printouts
97   mgr->SetDebugLevel(0);
98   if (!mgr->InitAnalysis())
99     return;
100   
101  mgr->PrintStatus();
102  
103  if(isGrid)
104    mgr->StartAnalysis("grid");
105  else  
106    mgr->StartAnalysis("local",chain);
107  
108 };
109
110
111 //__________________________________________________________________________________________
112 // This helper macros creates a chain of ESD files for you. Source can be either a text
113 // file with the file paths or a directory. In the latter case all ESD files in all subdirectories
114 // are considered.
115 //
116 // Inspired by: Jan.Fiete.Grosse-Oetringhaus@cern.ch
117
118 TChain* CreateESDChain(const char* aDataDir = "ESDfiles.txt", Int_t aRuns = 20, Int_t offset = 0, Bool_t addFileName = kFALSE, Bool_t addFriend = kFALSE, const char* check = 0)
119 {
120   // creates chain of files in a given directory or file containing a list.
121   // In case of directory the structure is expected as:
122   // <aDataDir>/<dir0>/AliESDs.root
123   // <aDataDir>/<dir1>/AliESDs.root
124   // ...
125   //
126   // if addFileName is true the list only needs to contain the directories that contain the AliESDs.root files
127   // if addFriend is true a file AliESDfriends.root is expected in the same directory and added to the chain as friend
128   // if check is != 0 the files that work are written back into the textfile with the name check
129
130   if (!aDataDir)
131     return 0;
132
133   Long_t id, size, flags, modtime;
134   if (gSystem->GetPathInfo(aDataDir, &id, &size, &flags, &modtime))
135   {
136     printf("%s not found.\n", aDataDir);
137     return 0;
138   }
139
140   TChain* chain = new TChain("esdTree");
141   TChain* chainFriend = 0;
142   
143   if (addFriend)
144     chainFriend = new TChain("esdFriendTree");
145
146   if (flags & 2)
147   {
148     TString execDir(gSystem->pwd());
149     TSystemDirectory* baseDir = new TSystemDirectory(".", aDataDir);
150     TList* dirList            = baseDir->GetListOfFiles();
151     Int_t nDirs               = dirList->GetEntries();
152     gSystem->cd(execDir);
153
154     Int_t count = 0;
155
156     for (Int_t iDir=0; iDir<nDirs; ++iDir)
157     {
158       TSystemFile* presentDir = (TSystemFile*) dirList->At(iDir);
159       if (!presentDir || !presentDir->IsDirectory() || strcmp(presentDir->GetName(), ".") == 0 || strcmp(presentDir->GetName(), "..") == 0)
160         continue;
161
162       if (offset > 0)
163       {
164         --offset;
165         continue;
166       }
167
168       if (count++ == aRuns)
169         break;
170
171       TString presentDirName(aDataDir);
172       presentDirName += "/";
173       presentDirName += presentDir->GetName();
174
175       chain->Add(presentDirName + "/AliESDs.root/esdTree");
176     }
177   }
178   else
179   {
180     // Open the input stream
181     ifstream in;
182     in.open(aDataDir);
183
184     ofstream outfile;
185     if (check)
186       outfile.open(check);
187
188     Int_t count = 0;
189
190     // Read the input list of files and add them to the chain
191     TString line;
192     while (in.good())
193     {
194       in >> line;
195
196       if (line.Length() == 0)
197         continue;
198
199       if (offset > 0)
200       {
201         offset--;
202         continue;
203       }
204
205       if (count++ == aRuns)
206         break;
207
208       TString esdFile(line);
209
210       if (addFileName)
211         esdFile += "/AliESDs.root";
212         
213       TString esdFileFriend(esdFile);
214       esdFileFriend.ReplaceAll("AliESDs.root", "AliESDfriends.root");
215         
216       if (check)
217       {
218         TFile* file = TFile::Open(esdFile);
219         if (!file)
220           continue;
221         file->Close();
222         
223         if (chainFriend)
224         {
225           TFile* file = TFile::Open(esdFileFriend);
226           if (!file)
227             continue;
228           file->Close();
229         }
230         
231         outfile << line.Data() << endl;
232         printf("%s\n", line.Data());
233       }        
234         
235         // add esd file
236       chain->Add(esdFile);
237
238         // add file
239       if (chainFriend)
240         chainFriend->Add(esdFileFriend);
241     }
242
243     in.close();
244     
245     if (check)
246       outfile.close();
247   }
248   
249   if (chainFriend)
250     chain->AddFriend(chainFriend);
251
252   return chain;
253 }
254
255 void ChainToTextFile(TChain* chain, const char* target)
256 {
257   // write a text list of the files in the chain
258   
259   TObjArray* list = chain->GetListOfFiles();
260   TIterator* iter = list->MakeIterator();
261   TObject* obj = 0;
262
263   ofstream outfile;
264   outfile.open(target);
265
266   while ((obj = iter->Next())) {
267     TString fileName(obj->GetTitle());
268     
269     outfile << fileName.Data() << endl;
270   }
271
272   outfile.close();
273
274   delete iter;
275
276
277 TObjArray* Chain2List(TChain* chain)
278 {
279   // returns a TObjArray of TObjStrings of the file names in the chain
280
281   TObjArray* result = new TObjArray;
282
283   for (Int_t i=0; i<chain->GetListOfFiles()->GetEntries(); i++)
284     result->Add(new TObjString(chain->GetListOfFiles()->At(i)->GetTitle()));
285
286   return result;
287 }
288
289 void LookupWrite(TChain* chain, const char* target)
290 {
291   // looks up the chain and writes the remaining files to the text file target
292
293   chain->Lookup();
294
295   ChainToTextFile(chain, target);
296 }
297
298 TChain* CreateChain(const char* treeName, const char* aDataDir, Int_t aRuns, Int_t offset = 0)
299 {
300   // creates chain of files in a given directory or file containing a list.
301
302   if (!treeName || !aDataDir)
303     return 0;
304
305   TChain* chain = new TChain(treeName);
306   
307   // Open the input stream
308   ifstream in;
309   in.open(aDataDir);
310
311   Int_t count = 0;
312
313   // Read the input list of files and add them to the chain
314   TString line;
315   while(in.good()) 
316   {
317     in >> line;
318       
319     if (line.Length() == 0)
320       continue;      
321     
322     if (offset > 0)
323     {
324       --offset;
325       continue;
326     }
327
328     if (count++ == aRuns)
329       break;
330
331     chain->Add(line);
332   }
333
334   in.close();
335
336   return chain;
337 }