bug fixed that was considered also under/overflow bins in vertex axis
[u/mrichter/AliRoot.git] / PWG0 / CreateESDChain.C
1 /* $Id$ */
2
3 /**************************************************************************
4  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5  *                                                                        *
6  * Author: The ALICE Off-line Project.                                    *
7  * Contributors are mentioned in the code where appropriate.              *
8  *                                                                        *
9  * Permission to use, copy, modify and distribute this software and its   *
10  * documentation strictly for non-commercial purposes is hereby granted   *
11  * without fee, provided that the above copyright notice appears in all   *
12  * copies and that both the copyright notice and this permission notice   *
13  * appear in the supporting documentation. The authors make no claims     *
14  * about the suitability of this software for any purpose. It is          *
15  * provided "as is" without express or implied warranty.                  *
16  **************************************************************************/
17
18 // This helper macros creates a chain of ESD files for you. Source can be either a text
19 // file with the file paths or a directory. In the latter case all ESD files in all subdirectories
20 // are considered.
21 //
22 // Author: Jan.Fiete.Grosse-Oetringhaus@cern.ch
23
24 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)
25 {
26   // creates chain of files in a given directory or file containing a list.
27   // In case of directory the structure is expected as:
28   // <aDataDir>/<dir0>/AliESDs.root
29   // <aDataDir>/<dir1>/AliESDs.root
30   // ...
31   //
32   // if addFileName is true the list only needs to contain the directories that contain the AliESDs.root files
33   // if addFriend is true a file AliESDfriends.root is expected in the same directory and added to the chain as friend
34   // if check is != 0 the files that work are written back into the textfile with the name check
35
36   if (!aDataDir)
37     return 0;
38
39   Long_t id, size, flags, modtime;
40   if (gSystem->GetPathInfo(aDataDir, &id, &size, &flags, &modtime))
41   {
42     printf("%s not found.\n", aDataDir);
43     return 0;
44   }
45
46   TChain* chain = new TChain("esdTree");
47   TChain* chainFriend = 0;
48   
49   if (addFriend)
50     chainFriend = new TChain("esdFriendTree");
51
52   if (flags & 2)
53   {
54     TString execDir(gSystem->pwd());
55     TSystemDirectory* baseDir = new TSystemDirectory(".", aDataDir);
56     TList* dirList            = baseDir->GetListOfFiles();
57     Int_t nDirs               = dirList->GetEntries();
58     gSystem->cd(execDir);
59
60     Int_t count = 0;
61
62     for (Int_t iDir=0; iDir<nDirs; ++iDir)
63     {
64       TSystemFile* presentDir = (TSystemFile*) dirList->At(iDir);
65       if (!presentDir || !presentDir->IsDirectory() || strcmp(presentDir->GetName(), ".") == 0 || strcmp(presentDir->GetName(), "..") == 0)
66         continue;
67
68       if (offset > 0)
69       {
70         --offset;
71         continue;
72       }
73
74       if (count++ == aRuns)
75         break;
76
77       TString presentDirName(aDataDir);
78       presentDirName += "/";
79       presentDirName += presentDir->GetName();
80
81       chain->Add(presentDirName + "/AliESDs.root/esdTree");
82     }
83   }
84   else
85   {
86     // Open the input stream
87     ifstream in;
88     in.open(aDataDir);
89
90     ofstream outfile;
91     if (check)
92       outfile.open(check);
93
94     Int_t count = 0;
95
96     // Read the input list of files and add them to the chain
97     TString line;
98     while (in.good())
99     {
100       in >> line;
101
102       if (line.Length() == 0)
103         continue;
104
105       if (offset > 0)
106       {
107         offset--;
108         continue;
109       }
110
111       if (count++ == aRuns)
112         break;
113
114       TString esdFile(line);
115
116       if (addFileName)
117         esdFile += "/AliESDs.root";
118         
119       TString esdFileFriend(esdFile);
120       esdFileFriend.ReplaceAll("AliESDs.root", "AliESDfriends.root");
121         
122       if (check)
123       {
124         TFile* file = TFile::Open(esdFile);
125         if (!file)
126           continue;
127         file->Close();
128         
129         if (chainFriend)
130         {
131           TFile* file = TFile::Open(esdFileFriend);
132           if (!file)
133             continue;
134           file->Close();
135         }
136         
137         outfile << line.Data() << endl;
138         printf("%s\n", line.Data());
139       }        
140         
141         // add esd file
142       chain->Add(esdFile);
143
144         // add file
145       if (chainFriend)
146         chainFriend->Add(esdFileFriend);
147     }
148
149     in.close();
150     
151     if (check)
152       outfile.close();
153   }
154   
155   if (chainFriend)
156     chain->AddFriend(chainFriend);
157
158   return chain;
159 }
160
161 void ChainToTextFile(TChain* chain, const char* target)
162 {
163   // write a text list of the files in the chain
164   
165   TObjArray* list = chain->GetListOfFiles();
166   TIterator* iter = list->MakeIterator();
167   TObject* obj = 0;
168
169   ofstream outfile;
170   outfile.open(target);
171
172   while ((obj = iter->Next())) {
173     TString fileName(obj->GetTitle());
174     
175     outfile << fileName.Data() << endl;
176   }
177
178   outfile.close();
179
180   delete iter;
181
182
183 TObjArray* Chain2List(TChain* chain)
184 {
185   // returns a TObjArray of TObjStrings of the file names in the chain
186
187   TObjArray* result = new TObjArray;
188
189   for (Int_t i=0; i<chain->GetListOfFiles()->GetEntries(); i++)
190     result->Add(new TObjString(chain->GetListOfFiles()->At(i)->GetTitle()));
191
192   return result;
193 }
194
195 void LookupWrite(TChain* chain, const char* target)
196 {
197   // looks up the chain and writes the remaining files to the text file target
198
199   chain->Lookup();
200
201   ChainToTextFile(chain, target);
202 }
203
204 TChain* CreateChain(const char* treeName, const char* aDataDir, Int_t aRuns, Int_t offset = 0)
205 {
206   // creates chain of files in a given directory or file containing a list.
207
208   if (!treeName || !aDataDir)
209     return 0;
210
211   TChain* chain = new TChain(treeName);
212   
213   // Open the input stream
214   ifstream in;
215   in.open(aDataDir);
216
217   Int_t count = 0;
218
219   // Read the input list of files and add them to the chain
220   TString line;
221   while(in.good()) 
222   {
223     in >> line;
224       
225     if (line.Length() == 0)
226       continue;      
227     
228     if (offset > 0)
229     {
230       --offset;
231       continue;
232     }
233
234     if (count++ == aRuns)
235       break;
236
237     chain->Add(line);
238   }
239
240   in.close();
241
242   return chain;
243 }