dc740de4 |
1 | /* $Id$ */ |
2 | |
f10a1859 |
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, Bool_t check = kFALSE) |
539b6cb4 |
25 | { |
0ab29cfa |
26 | // creates chain of files in a given directory or file containing a list. |
27 | // In case of directory the structure is expected as: |
dc740de4 |
28 | // <aDataDir>/<dir0>/AliESDs.root |
dc740de4 |
29 | // <aDataDir>/<dir1>/AliESDs.root |
dc740de4 |
30 | // ... |
f10a1859 |
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 |
dc740de4 |
34 | |
539b6cb4 |
35 | if (!aDataDir) |
36 | return 0; |
37 | |
0ab29cfa |
38 | Long_t id, size, flags, modtime; |
39 | if (gSystem->GetPathInfo(aDataDir, &id, &size, &flags, &modtime)) |
40 | { |
bdfe2916 |
41 | printf("%s not found.\n", aDataDir); |
0ab29cfa |
42 | return 0; |
43 | } |
44 | |
539b6cb4 |
45 | TChain* chain = new TChain("esdTree"); |
f10a1859 |
46 | TChain* chainFriend = 0; |
47 | |
48 | if (addFriend) |
49 | chainFriend = new TChain("esdFriendTree"); |
539b6cb4 |
50 | |
0ab29cfa |
51 | if (flags & 2) |
539b6cb4 |
52 | { |
0ab29cfa |
53 | TString execDir(gSystem->pwd()); |
54 | TSystemDirectory* baseDir = new TSystemDirectory(".", aDataDir); |
55 | TList* dirList = baseDir->GetListOfFiles(); |
56 | Int_t nDirs = dirList->GetEntries(); |
57 | gSystem->cd(execDir); |
539b6cb4 |
58 | |
0ab29cfa |
59 | Int_t count = 0; |
dc740de4 |
60 | |
0ab29cfa |
61 | for (Int_t iDir=0; iDir<nDirs; ++iDir) |
62 | { |
63 | TSystemFile* presentDir = (TSystemFile*) dirList->At(iDir); |
64 | if (!presentDir || !presentDir->IsDirectory() || strcmp(presentDir->GetName(), ".") == 0 || strcmp(presentDir->GetName(), "..") == 0) |
65 | continue; |
49dc84d9 |
66 | |
0ab29cfa |
67 | if (offset > 0) |
68 | { |
69 | --offset; |
70 | continue; |
71 | } |
539b6cb4 |
72 | |
0ab29cfa |
73 | if (count++ == aRuns) |
74 | break; |
539b6cb4 |
75 | |
0ab29cfa |
76 | TString presentDirName(aDataDir); |
77 | presentDirName += "/"; |
78 | presentDirName += presentDir->GetName(); |
37dbb69e |
79 | |
0ab29cfa |
80 | chain->Add(presentDirName + "/AliESDs.root/esdTree"); |
81 | } |
82 | } |
83 | else |
84 | { |
85 | // Open the input stream |
86 | ifstream in; |
87 | in.open(aDataDir); |
dc740de4 |
88 | |
0ab29cfa |
89 | Int_t count = 0; |
37dbb69e |
90 | |
0ab29cfa |
91 | // Read the input list of files and add them to the chain |
f10a1859 |
92 | TString line; |
93 | while(in.good()) |
94 | { |
95 | in >> line; |
96 | |
97 | if (line.Length() == 0) |
98 | continue; |
99 | |
0ab29cfa |
100 | if (offset > 0) |
101 | { |
102 | --offset; |
103 | continue; |
104 | } |
37dbb69e |
105 | |
0ab29cfa |
106 | if (count++ == aRuns) |
107 | break; |
37dbb69e |
108 | |
f10a1859 |
109 | TString esdFile(line); |
110 | |
111 | if (addFileName) |
112 | esdFile += "/AliESDs.root"; |
113 | |
114 | TString esdFileFriend(esdFile); |
115 | esdFileFriend.ReplaceAll("AliESDs.root", "AliESDfriends.root"); |
116 | |
117 | if (check) |
118 | { |
119 | TFile* file = TFile::Open(esdFile); |
120 | if (!file) |
121 | continue; |
122 | file->Close(); |
123 | |
124 | if (chainFriend) |
125 | { |
126 | TFile* file = TFile::Open(esdFileFriend); |
127 | if (!file) |
128 | continue; |
129 | file->Close(); |
130 | } |
131 | |
132 | printf("%s\n", line.Data()); |
133 | } |
134 | |
0ab29cfa |
135 | // add esd file |
f10a1859 |
136 | chain->Add(esdFile); |
137 | |
138 | // add file |
139 | if (chainFriend) |
140 | chainFriend->Add(esdFileFriend); |
16e24ca3 |
141 | } |
142 | |
0ab29cfa |
143 | in.close(); |
37dbb69e |
144 | } |
f10a1859 |
145 | |
146 | if (chainFriend) |
147 | chain->AddFriend(chainFriend); |
37dbb69e |
148 | |
37dbb69e |
149 | return chain; |
150 | } |
7f8cf861 |
151 | |
6f052bf3 |
152 | void ChainToTextFile(TChain* chain, const char* target) |
7f8cf861 |
153 | { |
6f052bf3 |
154 | // write a text list of the files in the chain |
155 | |
7f8cf861 |
156 | TObjArray* list = chain->GetListOfFiles(); |
157 | TIterator* iter = list->MakeIterator(); |
158 | TObject* obj = 0; |
159 | |
160 | ofstream outfile; |
161 | outfile.open(target); |
162 | |
163 | while ((obj = iter->Next())) |
f10a1859 |
164 | outfile << obj->GetTitle() << endl; |
7f8cf861 |
165 | |
166 | outfile.close(); |
167 | |
168 | delete iter; |
6f052bf3 |
169 | } |
170 | |
171 | void LookupWrite(TChain* chain, const char* target) |
172 | { |
173 | // looks up the chain and writes the remaining files to the text file target |
174 | |
175 | chain->Lookup(); |
176 | |
177 | ChainToTextFile(chain, target); |
7f8cf861 |
178 | } |