One more correction to allow "/" in the file name
[u/mrichter/AliRoot.git] / PHOS / AliPHOSGetter.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id:  */
17
18 /* $Log:
19    29.05.2001 Yuri Kharlov:
20               Everywhere reading the treese TTree->GetEvent(i)
21               is replaced by reading the branches TBranch->GetEntry(0)
22 */
23
24 //_________________________________________________________________________
25 //  A singleton. This class should be used in the analysis stage to get 
26 //  reconstructed objects: Digits, RecPoints, TrackSegments and RecParticles,
27 //  instead of directly reading them from galice.root file. This container 
28 //  ensures, that one reads Digits, made of these particular digits, RecPoints, 
29 //  made of these particular RecPoints, TrackSegments and RecParticles. 
30 //  This becomes non trivial if there are several identical branches, produced with
31 //  different set of parameters. 
32 //
33 //  An example of how to use (see also class AliPHOSAnalyser):
34 //  AliPHOSGetter * gime = AliPHOSGetter::GetInstance("galice.root","test") ;
35 //  for(Int_t irecp = 0; irecp < gime->NRecParticles() ; irecp++)
36 //     AliPHOSRecParticle * part = gime->RecParticle(1) ;
37 //     ................
38 //  please->GetEvent(event) ;    // reads new event from galice.root
39 //                  
40 //*-- Author: Yves Schutz (SUBATECH) & Dmitri Peressounko (RRC KI & SUBATECH)
41 //*--         Completely redesigned by Dmitri Peressounko March 2001  
42 //
43 //*-- YS June 2001 : renamed the original AliPHOSIndexToObject and make
44 //*--         systematic usage of TFolders without changing the interface        
45 //////////////////////////////////////////////////////////////////////////////
46
47
48 // --- ROOT system ---
49
50 #include "TFile.h"
51 #include "TTree.h"
52 #include "TROOT.h"
53 #include "TObjString.h"
54 #include "TFolder.h"
55 #include "TParticle.h"
56
57 // --- Standard library ---
58 #include <iostream.h>
59
60 // --- AliRoot header files ---
61
62 #include "AliRun.h"
63 #include "AliConfig.h"
64 #include "AliPHOSGetter.h"
65 #include "AliPHOS.h"
66 #include "AliPHOSDigitizer.h"
67 #include "AliPHOSSDigitizer.h"
68 #include "AliPHOSClusterizer.h"
69 #include "AliPHOSClusterizerv1.h"
70 #include "AliPHOSTrackSegmentMaker.h"
71 #include "AliPHOSTrackSegmentMakerv1.h"
72 #include "AliPHOSTrackSegment.h"
73 #include "AliPHOSPID.h" 
74 #include "AliPHOSPIDv1.h" 
75 #include "AliPHOSGeometry.h"
76
77 ClassImp(AliPHOSGetter)
78   
79   AliPHOSGetter * AliPHOSGetter::fgObjGetter = 0 ; 
80
81 //____________________________________________________________________________ 
82 AliPHOSGetter::AliPHOSGetter(const char* headerFile, const char* branchTitle )
83 {
84   //Initialize  all lists
85
86   fDebug = 0 ; 
87   
88   fHeaderFile         = headerFile ; 
89   fBranchTitle        = branchTitle ;
90   fSDigitsTitle       = branchTitle ; 
91   fDigitsTitle        = branchTitle ; 
92   fRecPointsTitle     = branchTitle ; 
93   fRecParticlesTitle  = branchTitle ; 
94   fTrackSegmentsTitle = branchTitle ; 
95
96   fPrimaries = new TObjArray(1) ;
97
98   fModuleFolder    = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/Run/Configuration/Modules")); 
99   fPrimariesFolder = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/RunMC/Event/Data")); 
100   fHitsFolder      = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/RunMC/Event/Data/Hits")); 
101   fSDigitsFolder   = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/RunMC/Event/Data/SDigits")); 
102   fDigitsFolder    = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/Run/Event/Data")); 
103   fRecoFolder      = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/Run/Event/RecData")); 
104   fQAFolder        = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/Run/Conditions/QA")); 
105   fTasksFolder     = dynamic_cast<TFolder*>(gROOT->FindObjectAny("Folders/Tasks")) ; 
106
107   if ( fHeaderFile != "aliroot"  ) { // to call the getter without a file
108
109     //open headers file
110     TFile * file = static_cast<TFile*>(gROOT->GetFile(fHeaderFile.Data() ) ) ;
111     
112     if(file == 0){    //if file was not opened yet, read gAlice
113       if(fHeaderFile.Contains("rfio")) // if we read file using HPSS
114         file =  TFile::Open(fHeaderFile.Data(),"update") ;
115       else
116         file = new TFile(fHeaderFile.Data(),"update") ;
117       
118       if (!file->IsOpen()) {
119         cerr << "ERROR : AliPHOSGetter::AliPHOSGetter -> Cannot open " << fHeaderFile.Data() << endl ; 
120         abort() ; 
121       }
122       
123       gAlice = static_cast<AliRun *>(file->Get("gAlice")) ;
124     }
125   }
126
127   if (!gAlice) {
128     cerr << "ERROR : AliPHOSGetter::AliPHOSGetter -> Cannot find gAlice in " << fHeaderFile.Data() << endl ; 
129     abort() ; 
130   }
131   if (!PHOS()) {
132     if (fDebug)
133       cout << "INFO: AliPHOSGetter -> Posting PHOS to Folders" << endl ; 
134     AliConfig * conf = AliConfig::Instance() ; 
135     conf->Add(static_cast<AliDetector*>(gAlice->GetDetector("PHOS"))) ; 
136     conf->Add(static_cast<AliModule*>(gAlice->GetDetector("PHOS"))) ; 
137   }
138   
139   fDebug=0;
140 }
141 //____________________________________________________________________________ 
142 AliPHOSGetter::~AliPHOSGetter(){
143
144 }
145
146 //____________________________________________________________________________ 
147 void AliPHOSGetter::CreateWhiteBoard() const
148 {
149
150 }
151
152 //____________________________________________________________________________ 
153 AliPHOSGetter * AliPHOSGetter::GetInstance()
154 {
155   // Returns the pointer of the unique instance already defined
156   
157   AliPHOSGetter * rv = 0 ;
158   if ( fgObjGetter )
159     rv = fgObjGetter ;
160   else
161     cout << "AliPHOSGetter::GetInstance ERROR: not yet initialized" << endl ;
162
163   return rv ;
164 }
165
166 //____________________________________________________________________________ 
167 AliPHOSGetter * AliPHOSGetter::GetInstance(const char* headerFile,
168                                            const char* branchTitle)
169 {
170   // Creates and returns the pointer of the unique instance
171   // Must be called only when the environment has changed 
172
173   if ( fgObjGetter )    
174     if((fgObjGetter->fBranchTitle.CompareTo(branchTitle) == 0) && 
175        (fgObjGetter->fHeaderFile.CompareTo(headerFile)==0))
176       return fgObjGetter ;
177     else
178       fgObjGetter->~AliPHOSGetter() ;  // delete it if already exists another version
179   
180   fgObjGetter = new AliPHOSGetter(headerFile,branchTitle) ; 
181   
182   // Posts a few item to the white board (folders)
183   // fgObjGetter->CreateWhiteBoard() ;
184     
185   return fgObjGetter ; 
186   
187 }
188
189 //____________________________________________________________________________ 
190 const AliPHOS * AliPHOSGetter::PHOS() 
191 {
192   // returns the PHOS object 
193   AliPHOS * phos = dynamic_cast<AliPHOS*>(fModuleFolder->FindObject("PHOS")) ;  
194   if (!phos) 
195     if (fDebug)
196       cout << "WARNING: AliPHOSGetter::PHOS -> PHOS module not found in Folders" << endl ; 
197   return phos ; 
198 }  
199
200 //____________________________________________________________________________ 
201 const AliPHOSGeometry * AliPHOSGetter::PHOSGeometry() 
202 {
203   AliPHOSGeometry * rv = 0 ; 
204   if (PHOS() )
205     rv =  PHOS()->GetGeometry() ;
206   return rv ; 
207
208
209 //____________________________________________________________________________ 
210 Bool_t AliPHOSGetter::PostPrimaries(void) const 
211 {  //------- Primaries ----------------------
212
213   // the hierarchy is //Folders/RunMC/Event/Data/Primaries
214   
215   TFolder * primariesFolder = dynamic_cast<TFolder*>(fPrimariesFolder->FindObject("Primaries")) ; 
216   if ( !primariesFolder ) {
217     if (fDebug) {
218       cout << "WARNING: AliPHOSGetter::Post Primaries -> Folder //" << fPrimariesFolder << "/Primaries/ not found!" << endl;
219       cout << "INFO:    AliPHOSGetter::Post Primaries -> Adding Folder //" << fPrimariesFolder << "/Primaries/"  << endl;
220     }
221     primariesFolder = fPrimariesFolder->AddFolder("Primaries", "Primaries particles from TreeK") ; 
222   }    
223   TClonesArray *primaries=  new TClonesArray("TParticle",1000) ;
224   primaries->SetName("Primaries") ;
225   primariesFolder->Add(primaries) ; 
226   
227   return kTRUE;
228
229
230 //____________________________________________________________________________ 
231 TObject** AliPHOSGetter::PrimariesRef(void) const 
232 {  //------- Primaries ----------------------
233
234   
235   // the hierarchy is //Folders/RunMC/Event/Data/Primaries
236   if ( !fPrimariesFolder ) {
237     cerr << "ERROR: AliPHOSGetter::Post PrimariesRef -> Folder //" << fPrimariesFolder << " not found!" << endl;
238     return 0;
239   }    
240  
241   TFolder * primariesFolder = dynamic_cast<TFolder *>(fPrimariesFolder->FindObject("Primaries")) ;
242   if ( !primariesFolder ) {
243     cerr << "ERROR: AliPHOSGetter::Post PrimariesRef -> Folder //" << fPrimariesFolder << "/Primaries/ not found!" << endl;  
244     return 0;
245   }
246  
247   TObject * p = primariesFolder->FindObject("Primaries") ;
248   if(!p) {
249     cerr << "ERROR: AliPHOSGetter::PrimariesRef -> " << primariesFolder->GetName() << "/Primaries not found !" << endl ; 
250     return 0 ;
251   }
252   else
253     return primariesFolder->GetListOfFolders()->GetObjectRef(p) ;
254 }
255
256 //____________________________________________________________________________ 
257 Bool_t AliPHOSGetter::PostHits(void) const 
258 {  //------- Hits ----------------------
259
260   // the hierarchy is //Folders/RunMC/Event/Data/PHOS/Hits
261   
262   TFolder * phosFolder = dynamic_cast<TFolder*>(fHitsFolder->FindObject("PHOS")) ; 
263   if ( !phosFolder ) {
264     if (fDebug) {
265       cout << "WARNING: AliPHOSGetter::Post H -> Folder //" << fHitsFolder << "/PHOS/ not found!" << endl;
266       cout << "INFO:    AliPHOSGetter::Post H -> Adding Folder //" << fHitsFolder << "/PHOS/"  << endl;
267     }
268     phosFolder = fHitsFolder->AddFolder("PHOS", "Hits from PHOS") ; 
269   }    
270   TClonesArray *hits=  new TClonesArray("AliPHOSHit",1000) ;
271   hits->SetName("Hits") ;
272   phosFolder->Add(hits) ; 
273   
274   return kTRUE;
275
276
277 //____________________________________________________________________________ 
278 TObject** AliPHOSGetter::HitsRef(void) const 
279 {  //------- Hits ----------------------
280
281   
282   // the hierarchy is //Folders/RunMC/Event/Data/PHOS/Hits
283   if ( !fHitsFolder ) {
284     cerr << "ERROR: AliPHOSGetter::Post H -> Folder //" << fHitsFolder << " not found!" << endl;
285     return 0;
286   }    
287  
288   TFolder * phosFolder = dynamic_cast<TFolder *>(fHitsFolder->FindObject("PHOS")) ;
289   if ( !phosFolder ) {
290     cerr << "ERROR: AliPHOSGetter::Post HRef -> Folder //" << fHitsFolder << "/PHOS/ not found!" << endl;  
291     return 0;
292   }
293  
294   TObject * h = phosFolder->FindObject("Hits") ;
295   if(!h) {
296     cerr << "ERROR: AliPHOSGetter::HRef -> " << phosFolder->GetName() << "/Hits not found !" << endl ; 
297     return 0 ;
298   }
299   else
300     return phosFolder->GetListOfFolders()->GetObjectRef(h) ;
301 }
302
303 //____________________________________________________________________________ 
304 Bool_t AliPHOSGetter::PostSDigits(const char * name, const char * headerFile) const 
305 {  //---------- SDigits -------------------------
306
307   
308   // the hierarchy is //Folders/RunMC/Event/Data/PHOS/SDigits/headerFile/sdigitsname
309   // because you can have sdigits from several hit files for mixing
310   
311   TFolder * phosFolder = dynamic_cast<TFolder*>(fSDigitsFolder->FindObject("PHOS")) ;
312   if ( !phosFolder ) {
313     if (fDebug) {
314       cout << "WARNING: AliPHOSGetter::Post S -> Folder //" << fSDigitsFolder << "/PHOS/ not found!" << endl;
315       cout << "INFO:    AliPHOSGetter::Post S -> Adding Folder //" << fHitsFolder << "/PHOS/" << endl;
316     }
317     phosFolder = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; 
318   }    
319   TString subdir(headerFile) ;
320   subdir.ReplaceAll("/","_") ; 
321   TFolder * phosSubFolder = dynamic_cast<TFolder*>(phosFolder->FindObject(subdir)) ; 
322   if ( !phosSubFolder ) 
323     phosSubFolder = phosFolder->AddFolder(subdir, ""); 
324   
325   TObject * sd  = phosSubFolder->FindObject(name); 
326   if ( sd ) {
327     if (fDebug)
328       cerr <<"INFO: AliPHOSGetter::Post S -> Folder " << subdir 
329            << " already exists!" << endl ;  
330   }else{
331     TClonesArray * sdigits = new TClonesArray("AliPHOSDigit",1) ;
332     sdigits->SetName(name) ;
333     phosSubFolder->Add(sdigits) ;
334   }
335   
336   return kTRUE;
337
338 //____________________________________________________________________________ 
339 TObject** AliPHOSGetter::SDigitsRef(const char * name, const char * file) const 
340 {  //------- SDigits ----------------------
341   
342   // the hierarchy is //Folders/RunMC/Event/Data/PHOS/SDigits/filename/SDigits
343
344   if ( !fSDigitsFolder ) {
345     cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //" << fSDigitsFolder << " not found!" << endl;
346     return 0;
347   }    
348  
349   TFolder * phosFolder = dynamic_cast<TFolder *>(fSDigitsFolder->FindObject("PHOS")) ;
350   if ( !phosFolder ) {
351     cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //" << fSDigitsFolder << "/PHOS/ not found!" << endl;
352     return 0;
353   }
354
355   TFolder * phosSubFolder = 0 ;
356   if(file)
357     phosSubFolder = dynamic_cast<TFolder *>(phosFolder->FindObject(file)) ;
358   else
359     phosSubFolder = dynamic_cast<TFolder *>(phosFolder->FindObject(fHeaderFile)) ;
360   
361   if(!phosSubFolder) {
362     cerr << "ERROR: AliPHOSGetter::Post SRef -> Folder //Folders/RunMC/Event/Data/PHOS/" << file << "not found!" << endl;
363     return 0;
364   }
365
366   TObject * dis = phosSubFolder->FindObject(name) ;
367   if(!dis)
368     return 0 ;
369   else
370     return phosSubFolder->GetListOfFolders()->GetObjectRef(dis) ;
371
372 }
373
374 //____________________________________________________________________________ 
375 Bool_t AliPHOSGetter::PostSDigitizer(AliPHOSSDigitizer * sdigitizer) const 
376 {  //---------- SDigitizer -------------------------
377     
378   // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname
379
380
381   TTask * sd  = dynamic_cast<TTask*>(fTasksFolder->FindObject("SDigitizer")) ; 
382
383   if ( !sd ) {
384     cerr << "ERROR: AliPHOSGetter::Post Ser -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl;
385     return kFALSE ;
386   }        
387   TTask * phos = dynamic_cast<TTask*>(sd->GetListOfTasks()->FindObject("PHOS")) ; 
388   if ( !phos )  {
389     if (fDebug) {
390       cout <<"WARNING: AliPHOSGetter::Post Ser ->//" << fTasksFolder << "/SDigitizer/PHOS/ not found!" << endl;  
391       cout <<"INFO: AliPHOSGetter::Post Ser -> Adding //" << fTasksFolder << "/SDigitizer/PHOS/" << endl;
392     }
393     phos = new TTask("PHOS", "") ; 
394     sd->Add(phos) ; 
395   } 
396   AliPHOSSDigitizer * phossd  = dynamic_cast<AliPHOSSDigitizer *>(phos->GetListOfTasks()->FindObject( sdigitizer->GetName() )); 
397   if (phossd) { 
398     if (fDebug)
399       cout << "INFO: AliPHOSGetter::Post Ser -> Task " << sdigitizer->GetName() << " already exists" << endl ; 
400     phos->GetListOfTasks()->Remove(phossd) ;
401   }
402   phos->Add(sdigitizer) ;       
403   return kTRUE; 
404   
405 }
406
407 //____________________________________________________________________________ 
408 TObject** AliPHOSGetter::SDigitizerRef(const char * name) const 
409 {  
410
411   TTask * sd  = dynamic_cast<TTask*>(fTasksFolder->FindObject("SDigitizer")) ; 
412   if ( !sd ) {
413     cerr << "ERROR: AliPHOSGetter::Post SerRef -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl;
414     abort();
415   }        
416
417   TTask * phos = dynamic_cast<TTask*>(sd->GetListOfTasks()->FindObject("PHOS")) ; 
418   if ( !phos )  {
419     cerr <<"ERROR: AliPHOSGetter::Post SerRef ->  //" << fTasksFolder << "/SDigitizer/PHOS not found!" << endl;
420     abort();
421   }        
422
423   TTask * task = dynamic_cast<TTask*>(phos->GetListOfTasks()->FindObject(name)) ; 
424
425   return phos->GetListOfTasks()->GetObjectRef(task) ;
426
427 }
428
429 //____________________________________________________________________________ 
430 Bool_t AliPHOSGetter::PostSDigitizer(const char * name, const char * file) const 
431 {  //---------- SDigitizer -------------------------
432   
433  // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname
434
435
436   TTask * sd  = dynamic_cast<TTask*>(fTasksFolder->FindObject("SDigitizer")) ; 
437   if ( !sd ) {
438     cerr << "ERROR: AliPHOSGetter::Post Ser -> Task //" << fTasksFolder << "/SDigitizer not found!" << endl;
439     return kFALSE ;
440   }        
441
442   TTask * phos = dynamic_cast<TTask*>(sd->GetListOfTasks()->FindObject("PHOS")) ; 
443   if ( !phos )  {
444     if (fDebug) {
445       cout <<"WARNING: AliPHOSGetter::Post Ser ->  //" << fTasksFolder << "/SDigitizer/PHOS/ not found!" << endl;
446       cout <<"INFO: AliPHOSGetter::Post Ser -> Adding  //" << fTasksFolder << "/SDigitizer/PHOS" << endl;
447     }
448     phos = new TTask("PHOS", "") ; 
449     sd->Add(phos) ; 
450   } 
451
452   TString sdname(name) ;
453   sdname.Append(":") ;
454   sdname.Append(file);
455   sdname.ReplaceAll("/","_") ; 
456   AliPHOSSDigitizer * phossd  = dynamic_cast<AliPHOSSDigitizer *>(phos->GetListOfTasks()->FindObject( sdname )); 
457   if (!phossd) {
458     phossd = new AliPHOSSDigitizer() ;  
459     //Note, we can not call constructor with parameters: it will call Getter and scrud up everething
460     phossd->SetName(sdname) ;
461     phossd->SetTitle(file) ;
462     phos->Add(phossd) ; 
463   }
464   return kTRUE; 
465   
466 }
467
468 //____________________________________________________________________________ 
469 Bool_t AliPHOSGetter::PostDigits(const char * name) const 
470 {  //---------- Digits -------------------------
471
472   // the hierarchy is //Folders/Run/Event/Data/PHOS/SDigits/name
473
474   TFolder * phosFolder  = dynamic_cast<TFolder*>(fDigitsFolder->FindObject("PHOS")) ;
475
476   if ( !phosFolder ) {
477     if (fDebug) {
478       cout << "WARNING: AliPHOSGetter::Post D -> Folder //" << fDigitsFolder << "/PHOS/ not found!" << endl;
479       cout << "INFO:    AliPHOSGetter::Post D -> Adding Folder //" << fDigitsFolder << "/PHOS/" << endl;
480     }
481     phosFolder = fDigitsFolder->AddFolder("PHOS", "Digits from PHOS") ;  
482   }    
483  
484   TObject*  dig = phosFolder->FindObject( name ) ;
485   if ( !dig ) {
486     TClonesArray * digits = new TClonesArray("AliPHOSDigit",1000) ;
487     digits->SetName(name) ;
488     phosFolder->Add(digits) ;  
489   }
490   return kTRUE; 
491 }
492
493 //____________________________________________________________________________ 
494 TObject** AliPHOSGetter::DigitsRef(const char * name) const 
495 { //------- Digits ----------------------
496   
497   // the hierarchy is //Folders/Run/Event/Data/PHOS/Digits/name
498
499   if ( !fDigitsFolder ) {
500     cerr << "ERROR: AliPHOSGetter::Post DRef -> Folder //" << fDigitsFolder << " not found!" << endl;
501     return 0;
502   }    
503   
504   TFolder * phosFolder  = dynamic_cast<TFolder*>(fDigitsFolder->FindObject("PHOS")) ; 
505   if ( !phosFolder ) {
506     cerr << "ERROR: AliPHOSGetter::DRef -> Folder //" << fDigitsFolder << "/PHOS/ not found!" << endl;
507     return 0;
508   }    
509
510   TObject * d = phosFolder->FindObject(name) ;
511   if(!d)
512     return 0 ;
513   else
514     return phosFolder->GetListOfFolders()->GetObjectRef(d) ;
515
516 }
517
518 //____________________________________________________________________________ 
519 Bool_t AliPHOSGetter::PostDigitizer(AliPHOSDigitizer * digitizer) const 
520 {  //---------- Digitizer -------------------------
521   
522   TTask * sd  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Digitizer")) ; 
523
524   if ( !sd ) {
525     cerr << "ERROR: AliPHOSGetter::Post Der -> Task //" << fTasksFolder << "/Digitizer not found!" << endl;
526     return kFALSE ;
527   }        
528   TTask * phos = dynamic_cast<TTask*>(sd->GetListOfTasks()->FindObject("PHOS")) ; 
529   if ( !phos )  {
530     if (fDebug) {
531       cout <<"WARNING: AliPHOSGetter::Post Der ->  //" << fTasksFolder << "/Digitizer/PHOS not found!" << endl;
532       cout <<"INFO: AliPHOSGetter::Post Der -> Adding //" << fTasksFolder << "/Digitizer/PHOS" << endl; 
533     }
534     phos = new TTask("PHOS", "") ; 
535     sd->Add(phos) ; 
536   } 
537
538     AliPHOSDigitizer * phosd = dynamic_cast<AliPHOSDigitizer*>(phos->GetListOfTasks()->FindObject(digitizer->GetName())) ; 
539     if (phosd) { 
540       phosd->Delete() ;
541       phos->GetListOfTasks()->Remove(phosd) ;
542     }
543     phos->Add(digitizer) ; 
544     return kTRUE; 
545 }  
546
547 //____________________________________________________________________________ 
548 Bool_t AliPHOSGetter::PostDigitizer(const char * name) const 
549 {  //---------- Digitizer -------------------------
550   
551  // the hierarchy is //Folders/Tasks/SDigitizer/PHOS/sdigitsname
552
553   TTask * d  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Digitizer")) ; 
554   if ( !d ) {
555     cerr << "ERROR: AliPHOSGetter::Post Der -> Task //" << fTasksFolder << "/Digitizer not found!" << endl;
556     return kFALSE ;
557   }        
558
559   TTask * phos = dynamic_cast<TTask*>(d->GetListOfTasks()->FindObject("PHOS")) ; 
560   if ( !phos )  {
561     if (fDebug) {
562       cout <<"WARNING: AliPHOSGetter::Post Der -> //" << fTasksFolder << "/Digitizer/PHOS not found!" << endl; 
563       cout <<"INFO: AliPHOSGetter::Post Der -> Adding //" << fTasksFolder << "/Digitizer/PHOS" << endl;
564     }
565     phos = new TTask("PHOS", "") ; 
566     d->Add(phos) ; 
567
568
569   AliPHOSDigitizer * phosd = dynamic_cast<AliPHOSDigitizer*>(phos->GetListOfTasks()->FindObject(name)) ; 
570   if (!phosd) { 
571     phosd = new AliPHOSDigitizer() ;
572     phosd->SetName(fDigitsTitle) ;
573     phosd->SetTitle(fHeaderFile) ;
574     phos->Add(phosd) ;
575   }
576   return kTRUE;  
577 }
578
579 //____________________________________________________________________________ 
580 TObject** AliPHOSGetter::DigitizerRef(const char * name) const 
581 {  
582   TTask * sd  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Digitizer")) ; 
583   if ( !sd ) {
584     cerr << "ERROR: AliPHOSGetter::Post DerRef -> Task //" << fTasksFolder << "/Digitizer not found!" << endl;
585     abort();
586   }        
587
588   TTask * phos = dynamic_cast<TTask*>(sd->GetListOfTasks()->FindObject("PHOS")) ; 
589   if ( !phos )  {
590     cerr <<"ERROR: AliPHOSGetter::Post DerRef ->  //" << fTasksFolder << "/Digitizer/PHOS" << endl;
591     abort();
592   }        
593
594   TTask * task = dynamic_cast<TTask*>(phos->GetListOfTasks()->FindObject(name)) ; 
595
596   return phos->GetListOfTasks()->GetObjectRef(task) ;
597
598 }
599  
600 //____________________________________________________________________________ 
601 Bool_t AliPHOSGetter::PostRecPoints(const char * name) const 
602 { // -------------- RecPoints -------------------------------------------
603   
604   // the hierarchy is //Folders/Run/Event/RecData/PHOS/EMCARecPoints/name
605   // the hierarchy is //Folders/Run/Event/RecData/PHOS/CPVRecPoints/name
606
607   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS")) ; 
608   
609   if ( !phosFolder ) {
610     if (fDebug) {
611       cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl;
612       cout << "INFO:    AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/" << endl;
613     }
614     phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ;  
615   }    
616   
617   // EMCA RecPoints 
618   TFolder * phosRPoEMCAFolder  = dynamic_cast<TFolder*>(phosFolder->FindObject("EMCARecPoints")) ;
619   if ( !phosRPoEMCAFolder ) {
620     if (fDebug) {
621       cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/EMCARecPoints/ not found!" << endl;
622       cout << "INFO:    AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/EMCARecPoints not found!" << endl;
623     }
624     phosRPoEMCAFolder = phosFolder->AddFolder("EMCARecPoints", "EMCA RecPoints from PHOS") ;  
625   }    
626   
627   TObject * erp = phosFolder->FindObject( name ) ;
628   if ( !erp )   {
629     TObjArray * emcrp = new TObjArray(100) ;
630     emcrp->SetName(name) ;
631     phosRPoEMCAFolder->Add(emcrp) ;  
632   }
633
634   // CPV RecPoints 
635   TFolder * phosRPoCPVFolder  = dynamic_cast<TFolder*>(phosFolder->FindObject("CPVRecPoints")) ;
636   if ( !phosRPoCPVFolder ) {
637     if (fDebug) {
638       cout << "WARNING: AliPHOSGetter::Post RPo -> Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/ not found!" << endl;
639       cout << "INFO:    AliPHOSGetter::Post Rpo -> Adding Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/" << endl;
640     }
641     phosRPoCPVFolder = phosFolder->AddFolder("CPVRecPoints", "CPV RecPoints from PHOS") ;  
642   }    
643   
644   TObject * crp =  phosRPoCPVFolder->FindObject( name ) ;
645   if ( !crp )   {
646     TObjArray * cpvrp = new TObjArray(100) ;
647     cpvrp->SetName(name) ;
648     phosRPoCPVFolder->Add(cpvrp) ;  
649   }
650   return kTRUE; 
651 }
652
653 //____________________________________________________________________________ 
654 TObject** AliPHOSGetter::EmcRecPointsRef(const char * name) const 
655 { // -------------- RecPoints -------------------------------------------
656   
657   // the hierarchy is //Folders/Run/Event/RecData/PHOS/EMCARecPoints/name
658    
659   if ( !fRecoFolder ) {
660     cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << " not found!" << endl;
661     return 0 ; 
662   }    
663
664   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS/EMCARecPoints")) ; 
665   if ( !phosFolder ) {
666     cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << "/PHOS/EMCARecPoints/ not found!" << endl;
667     return 0;
668   }    
669
670
671   TObject * erp = phosFolder->FindObject(name ) ;
672   if ( !erp )   {
673     return 0 ;
674   }
675   return phosFolder->GetListOfFolders()->GetObjectRef(erp) ;
676
677
678
679 //____________________________________________________________________________ 
680 TObject** AliPHOSGetter::CpvRecPointsRef(const char * name) const 
681 { // -------------- RecPoints -------------------------------------------
682   
683   // the hierarchy is //Folders/Run/Event/RecData/PHOS/CPVRecPoints/name
684    
685   if ( !fRecoFolder ) {
686     cerr << "ERROR: AliPHOSGetter::EmcRecPointsRef -> Folder //" << fRecoFolder << " not found!" << endl;
687     return 0 ; 
688   }    
689
690   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS/CPVRecPoints")) ; 
691   if ( !phosFolder ) {
692     cerr << "ERROR: AliPHOSGetter::CpvRecPointsRef -> Folder //" << fRecoFolder << "/PHOS/CPVRecPoints/" << endl;
693     return 0;
694   }    
695
696   TObject * crp = phosFolder->FindObject(name ) ;
697   if ( !crp )   {
698     return 0 ;
699   }
700   return phosFolder->GetListOfFolders()->GetObjectRef(crp) ;
701
702
703
704 //____________________________________________________________________________ 
705 Bool_t AliPHOSGetter::PostClusterizer(AliPHOSClusterizer * clu) const 
706 { // ------------------ AliPHOSClusterizer ------------------------
707   
708   // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname
709
710   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
711
712   if ( !tasks ) {
713     cerr << "ERROR: AliPHOSGetter::Post Rer -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
714     return kFALSE ;
715   }        
716         
717   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
718   if ( !phos )  {
719     if (fDebug) {
720       cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/ReconstructionerPHOS not found!" << endl; 
721       cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; 
722     }
723     phos = new TTask("PHOS", "") ; 
724     tasks->Add(phos) ; 
725   } 
726
727   AliPHOSClusterizer * phoscl = dynamic_cast<AliPHOSClusterizer*>(phos->GetListOfTasks()->FindObject(clu->GetName())) ; 
728   if (phoscl) { 
729     if (fDebug)
730       cout << "INFO: AliPHOSGetter::Post Rer -> Task " << clu->GetName() << " already exists" << endl ; 
731     phoscl->Delete() ; 
732     phos->GetListOfTasks()->Remove(phoscl) ;
733   }
734   phos->Add(clu) ;      
735   return kTRUE; 
736
737
738 //____________________________________________________________________________ 
739 TObject** AliPHOSGetter::ClusterizerRef(const char * name) const 
740 { // ------------------ AliPHOSClusterizer ------------------------
741   
742   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
743
744   if ( !tasks ) {
745     cerr << "ERROR: AliPHOSGetter::Post RerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
746     return kFALSE ;
747   }        
748         
749   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
750   if ( !phos )  {
751     cerr <<"WARNING: AliPHOSGetter::Post RerRef -> //" << fTasksFolder << "/Reconstructioner/PHOS" << endl; 
752     return 0 ; 
753   }   
754
755   TList * l = phos->GetListOfTasks() ; 
756   TIter it(l) ;
757   TTask * task ;
758   TTask * clu = 0 ;
759   TString cluname(name) ;
760   cluname+=":clu-" ;
761   while((task = static_cast<TTask *>(it.Next()) )){
762     TString taskname(task->GetName()) ;
763     if(taskname.BeginsWith(cluname)){
764       clu = task ;
765       break ;
766     }
767   }
768
769   if(clu) 
770     return l->GetObjectRef(clu) ;
771   else
772     return 0 ;
773 }
774
775 //____________________________________________________________________________ 
776 Bool_t AliPHOSGetter::PostClusterizer(const char * name) const 
777 { // ------------------ AliPHOSClusterizer ------------------------
778
779   // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname
780   
781   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
782
783   if ( !tasks ) {
784     cerr << "ERROR: AliPHOSGetter::Post Rer -> Task//" << fTasksFolder << "/Reconstructioner not found!" << endl; 
785     return kFALSE ;
786   }        
787   
788   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
789   if ( !phos )  {
790     if (fDebug) {
791       cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl;
792       cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl;
793     }
794     phos = new TTask("PHOS", "") ; 
795     tasks->Add(phos) ; 
796   } 
797
798   AliPHOSClusterizer * phoscl = new AliPHOSClusterizerv1() ;
799   TString clun(name) ;
800   clun+=":clu-v1" ;
801   phoscl->SetName(clun) ;
802   phos->Add(phoscl) ;
803   return kTRUE; 
804   
805 }
806
807 //____________________________________________________________________________ 
808 Bool_t AliPHOSGetter::PostTrackSegments(const char * name) const 
809 { // ---------------TrackSegments -----------------------------------
810   
811   // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name
812
813   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS")) ; 
814   
815   if ( !phosFolder ) {
816     if (fDebug) {
817       cout << "WARNING: AliPHOSGetter::Post TS -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl;
818       cout << "INFO:    AliPHOSGetter::Post TS -> Adding Folder //" << fRecoFolder << "/PHOS" << endl;
819     }
820     phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ;  
821   }    
822
823   TFolder * phosTSFolder  = dynamic_cast<TFolder*>(phosFolder->FindObject("TrackSegments")) ;
824   if ( !phosTSFolder ) {
825     if (fDebug) {
826       cout << "WARNING: AliPHOSGetter::Post TS -> Folder//" << fRecoFolder << "/PHOS/TrackSegments/ not found!" << endl; 
827       cout << "INFO:    AliPHOSGetter::Post TS -> Adding Folder //" << fRecoFolder << "/PHOS/TrackSegments/" << endl; 
828     }
829     phosTSFolder = phosFolder->AddFolder("TrackSegments", "TrackSegments from PHOS") ;  
830   }    
831   
832   TObject * tss =  phosTSFolder->FindObject( name ) ;
833   if (!tss) {
834     TClonesArray * ts = new TClonesArray("AliPHOSTrackSegment",100) ;
835     ts->SetName(name) ;
836     phosTSFolder->Add(ts) ;  
837   }
838   return kTRUE; 
839
840
841 //____________________________________________________________________________ 
842 TObject** AliPHOSGetter::TrackSegmentsRef(const char * name) const 
843 { // ---------------TrackSegments -----------------------------------
844   
845   // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name
846
847  if ( !fRecoFolder ) {
848     cerr << "ERROR: AliPHOSGetter::TrackSegmentsRef -> Folder //" << fRecoFolder << "not found!" << endl;
849     return 0 ; 
850   }    
851
852   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS/TrackSegments")) ; 
853   if ( !phosFolder ) {
854     cerr << "ERROR: AliPHOSGetter::TrackSegmentsRef -> Folder //" << fRecoFolder << "/PHOS/TrackSegments/ not found!" << endl;
855     return 0;
856   }    
857   
858   TObject * tss =  phosFolder->FindObject(name) ;
859   if (!tss) {
860     return 0 ;  
861   }
862   return phosFolder->GetListOfFolders()->GetObjectRef(tss) ;
863
864
865 //____________________________________________________________________________ 
866 Bool_t AliPHOSGetter::PostTrackSegmentMaker(AliPHOSTrackSegmentMaker * tsmaker) const 
867 { //------------Track Segment Maker ------------------------------
868   
869   // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname
870
871   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
872
873   if ( !tasks ) {
874     cerr << "ERROR: AliPHOSGetter::Post Ter -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
875     return kFALSE ;
876   }        
877         
878   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
879   if ( !phos )  {
880     if (fDebug) {
881       cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; 
882       cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl;
883     }
884     phos = new TTask("PHOS", "") ; 
885     tasks->Add(phos) ; 
886   } 
887
888   AliPHOSTrackSegmentMaker * phosts = 
889     dynamic_cast<AliPHOSTrackSegmentMaker*>(phos->GetListOfTasks()->FindObject(tsmaker->GetName())) ; 
890   if (phosts) { 
891     phosts->Delete() ;
892     phos->GetListOfTasks()->Remove(phosts) ;
893   }
894   phos->Add(tsmaker) ;      
895   return kTRUE; 
896   
897
898 //____________________________________________________________________________ 
899 Bool_t AliPHOSGetter::PostTrackSegmentMaker(const char * name) const 
900 { //------------Track Segment Maker ------------------------------
901   
902   // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname
903   
904   
905   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
906   
907   if ( !tasks ) {
908     cerr << "ERROR: AliPHOSGetter::Post Ter -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
909     return kFALSE ;
910   }        
911   
912   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
913   if ( !phos )  {
914     if (fDebug) {
915       cout <<"WARNING: AliPHOSGetter::Post Rer -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; 
916       cout <<"INFO: AliPHOSGetter::Post Rer -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl;
917     }
918     phos = new TTask("PHOS", "") ; 
919     tasks->Add(phos) ; 
920   } 
921
922   AliPHOSTrackSegmentMaker * phosts = 
923     dynamic_cast<AliPHOSTrackSegmentMaker*>(phos->GetListOfTasks()->FindObject(name)) ; 
924   if (!phosts) { 
925     phosts = new AliPHOSTrackSegmentMakerv1() ;
926     TString tsn(name);
927     tsn+=":tsm-v1" ;
928     phosts->SetName(tsn) ;
929     phos->Add(phosts) ;      
930   }
931   return kTRUE; 
932   
933
934
935 //____________________________________________________________________________ 
936 TObject** AliPHOSGetter::TSMakerRef(const char * name) const 
937 { //------------Track Segment Maker ------------------------------
938   
939   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
940
941   if ( !tasks ) {
942     cerr << "ERROR: AliPHOSGetter::Post TerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
943     return kFALSE ;
944   }        
945         
946   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
947   if ( !phos )  {
948     cerr <<"WARNING: AliPHOSGetter::Post TerRef -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; 
949     return 0 ; 
950   }   
951
952   TList * l = phos->GetListOfTasks() ; 
953   TIter it(l) ;
954   TTask * task ;
955   TTask * tsm = 0 ;
956   TString tsmname(name) ;
957   tsmname+=":tsm-" ;
958   while((task = static_cast<TTask *>(it.Next()) )){
959     TString taskname(task->GetName()) ;
960     if(taskname.BeginsWith(tsmname)){
961       tsm = task ;
962       break ;
963     }
964   }
965   
966   if(tsm) 
967     return l->GetObjectRef(tsm) ;
968   else
969     return 0 ;
970   
971
972
973 //____________________________________________________________________________ 
974 Bool_t AliPHOSGetter::PostRecParticles(const char * name) const 
975 {  // -------------------- RecParticles ------------------------
976   
977   // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name
978
979   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS")) ; 
980   
981   if ( !phosFolder ) {
982     if (fDebug) {
983       cout << "WARNING: AliPHOSGetter::Post RPa -> Folder //" << fRecoFolder << "/PHOS/ not found!" << endl;
984       cout << "INFO:    AliPHOSGetter::Post Rpa -> Adding Folder //" << fRecoFolder << "/PHOS/" << endl;
985     }
986     phosFolder = fRecoFolder->AddFolder("PHOS", "Reconstructed data from PHOS") ;  
987   }    
988
989  TFolder * phosRPaFolder  = dynamic_cast<TFolder*>(phosFolder->FindObject("RecParticles")) ;
990   if ( !phosRPaFolder ) {
991     if (fDebug) {
992       cout << "WARNING: AliPHOSGetter::Post RPa -> Folder //" << fRecoFolder << "/PHOS/RecParticles/ not found!" << endl;
993       cout << "INFO:    AliPHOSGetter::Post RPa -> Adding Folder //" << fRecoFolder << "/PHOS/RecParticles/" << endl;
994     }
995     phosRPaFolder = phosFolder->AddFolder("RecParticles", "RecParticles from PHOS") ;  
996   } 
997
998   TObject * rps = phosRPaFolder->FindObject( name )  ;
999   if ( !rps ) {
1000     TClonesArray * rp = new TClonesArray("AliPHOSRecParticle",100) ;
1001     rp->SetName(name) ;    
1002     phosRPaFolder->Add(rp) ;  
1003   }
1004   return kTRUE; 
1005
1006
1007 //____________________________________________________________________________ 
1008 TObject** AliPHOSGetter::RecParticlesRef(const char * name) const 
1009 { // ---------------TrackSegments -----------------------------------
1010   
1011   // the hierarchy is //Folders/Run/Event/RecData/PHOS/TrackSegments/name
1012
1013  if ( !fRecoFolder ) {
1014     cerr << "ERROR: AliPHOSGetter::RecParticlesRef -> Folder//" << fRecoFolder << " not found!" << endl; 
1015     return 0 ; 
1016   }    
1017
1018   TFolder * phosFolder  = dynamic_cast<TFolder*>(fRecoFolder->FindObject("PHOS/RecParticles")) ; 
1019   if ( !phosFolder ) {
1020     cerr << "ERROR: AliPHOSGetter::RecParticlesRef -> Folder //" << fRecoFolder << "/PHOS/RecParticles/ not found!" << endl;
1021     return 0;
1022   }    
1023
1024   TObject * tss =  phosFolder->FindObject(name  ) ;
1025   if (!tss) {
1026     return 0 ;  
1027   }
1028   return phosFolder->GetListOfFolders()->GetObjectRef(tss) ;
1029 }
1030
1031 //____________________________________________________________________________ 
1032 Bool_t AliPHOSGetter::PostPID(AliPHOSPID * pid) const 
1033 {      // ------------AliPHOS PID -----------------------------
1034
1035   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
1036
1037   if ( !tasks ) {
1038     cerr << "ERROR: AliPHOSGetter::Post Per -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
1039     return kFALSE ;
1040   }        
1041   
1042   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
1043   if ( !phos )  {
1044     if (fDebug) {
1045       cout <<"WARNING: AliPHOSGetter::Post Per -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; 
1046       cout <<"INFO: AliPHOSGetter::Post Per -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl;
1047     }
1048     phos = new TTask("PHOS", "") ; 
1049     tasks->Add(phos) ; 
1050   } 
1051
1052   AliPHOSPID * phospid = dynamic_cast<AliPHOSPID*>(phos->GetListOfTasks()->FindObject(pid->GetName())) ; 
1053   if (phospid) { 
1054     if (fDebug)
1055       cout << "INFO: AliPHOSGetter::Post Per -> Task " << pid->GetName()
1056            << " already exists" << endl ; 
1057     phos->GetListOfTasks()->Remove(phospid) ;
1058   }
1059   
1060   phos->Add(pid) ;      
1061   return kTRUE; 
1062
1063
1064 //____________________________________________________________________________ 
1065 Bool_t AliPHOSGetter::PostPID(const char * name) const 
1066 {     
1067   // the hierarchy is //Folders/Tasks/Reconstructioner/PHOS/sdigitsname
1068   
1069   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
1070
1071   if ( !tasks ) {
1072     cerr << "ERROR: AliPHOSGetter::Post Per -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
1073     return kFALSE ;
1074   }        
1075   
1076   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
1077   if ( !phos )  {
1078     if (fDebug) {
1079       cout <<"WARNING: AliPHOSGetter::Post Per -> //" << fTasksFolder << "/Reconstructioner/PHOS not found!" << endl; 
1080       cout <<"INFO: AliPHOSGetter::Post Per -> Adding //" << fTasksFolder << "/Reconstructioner/PHOS" << endl;
1081     }
1082     phos = new TTask("PHOS", "") ; 
1083     tasks->Add(phos) ; 
1084   } 
1085
1086   TList * l = phos->GetListOfTasks() ;   
1087   TIter it(l) ;
1088   TString pidname(name) ;
1089   pidname+=":pid" ; 
1090   TTask * task ;
1091   while((task = static_cast<TTask *>(it.Next()) )){
1092     TString taskname(task->GetName()) ;
1093     if(taskname.BeginsWith(pidname))
1094       return kTRUE ;
1095   }
1096  
1097   AliPHOSPIDv1 * phospid = new AliPHOSPIDv1() ;
1098   pidname+="-v1" ;
1099   phospid->SetName(pidname) ;
1100   phos->Add(phospid) ;      
1101   
1102   return kTRUE; 
1103
1104
1105 //____________________________________________________________________________ 
1106 TObject** AliPHOSGetter::PIDRef(const char * name) const 
1107 { //------------PID ------------------------------
1108
1109   TTask * tasks  = dynamic_cast<TTask*>(fTasksFolder->FindObject("Reconstructioner")) ; 
1110
1111   if ( !tasks ) {
1112     cerr << "ERROR: AliPHOSGetter::Post PerRef -> Task //" << fTasksFolder << "/Reconstructioner not found!" << endl;
1113     return kFALSE ;
1114   }        
1115         
1116   TTask * phos = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
1117   if ( !phos )  {
1118     cerr <<"WARNING: AliPHOSGetter::Post PerRef -> //" << fTasksFolder << "/ReconstructionerPHOS not found!" << endl; 
1119     return 0 ; 
1120   }   
1121   
1122   TList * l = phos->GetListOfTasks() ; 
1123   TIter it(l) ;
1124   TTask * task ;
1125   TTask * pid = 0 ;
1126   TString pidname(name) ;
1127   pidname+=":pid-" ;
1128   while((task = static_cast<TTask *>(it.Next()) )){
1129     TString taskname(task->GetName()) ;
1130     if(taskname.BeginsWith(pidname)){
1131       pid = task ;
1132       break ;
1133     }
1134   }
1135   
1136   if(pid) 
1137     return l->GetObjectRef(pid) ;
1138   else
1139     return 0 ;
1140   
1141
1142
1143 //____________________________________________________________________________ 
1144 Bool_t AliPHOSGetter::PostQA(void) const 
1145 { // ------------------ QA ---------------------------------
1146
1147   // the hierarchy is //Folders/Run/Conditions/QA/PHOS/alarmsName
1148
1149   TFolder * phosFolder = dynamic_cast<TFolder*>(fQAFolder->FindObject("PHOS")) ; 
1150   if ( !phosFolder ) {
1151     if (fDebug) {
1152       cout << "WARNING: AliPHOSGetter::Post Q -> Folder //" << fQAFolder << "/PHOS/ not found!" << endl;
1153       cout << "INFO:    AliPHOSGetter::Post Q -> Adding Folder //" << fQAFolder << "/PHOS/" << endl;
1154     }
1155     phosFolder = fQAFolder->AddFolder("PHOS", "QA from PHOS") ; 
1156   }      
1157
1158   return kTRUE;
1159 }
1160
1161 //____________________________________________________________________________ 
1162 TObject** AliPHOSGetter::AlarmsRef(void) const 
1163 {  //------- Alarms ----------------------
1164
1165   
1166   // the hierarchy is //Folders/Run/Conditions/QA/PHOS
1167   if ( !fQAFolder ) {
1168     cerr << "ERROR: AliPHOSGetter::Post QRef -> Folder //" << fQAFolder << " not found!" << endl;
1169     return 0;
1170   }    
1171  
1172   TFolder * phosFolder = dynamic_cast<TFolder *>(fQAFolder->FindObject("PHOS")) ;
1173   if ( !phosFolder ) {
1174     cerr << "ERROR: AliPHOSGetter::Post QRef -> Folder //" << fQAFolder << "/PHOS/ not found!" << endl;
1175     return 0;
1176   }
1177    
1178   return fQAFolder->GetListOfFolders()->GetObjectRef(phosFolder) ;
1179 }
1180
1181 //____________________________________________________________________________ 
1182 const TParticle * AliPHOSGetter::Primary(Int_t index) const
1183 {
1184   // Return primary particle numbered by <index>
1185
1186   if(index < 0) 
1187     return 0 ;
1188   TParticle *  p = gAlice->Particle(index) ; 
1189   if (p->GetFirstMother() != -1 ) {
1190     cout << "AliPHOSGetter::Primary : Not a primary " << endl ; 
1191   }
1192   return p ; 
1193   
1194   
1195 //   Int_t primaryIndex = index % 10000000 ; 
1196 //   Int_t primaryList = (Int_t ) ((index-primaryIndex)/10000000.)  ;
1197   
1198 //   if ( primaryList > 0  ) {
1199 //     if (fDebug) {
1200 //       cout << " Getter does not support currently Mixing of primary " << endl ;
1201 //       cout << "   can not return primary: " << index<< " (list "<< primaryList<< " primary # " << primaryIndex << " )"<<endl ;
1202 //     }
1203 //     return 0;
1204 //   }
1205   
1206 //   return gAlice->Particle(primaryIndex) ;
1207   
1208 }
1209
1210 //____________________________________________________________________________ 
1211 const TParticle * AliPHOSGetter::Secondary(TParticle* p, Int_t index) const
1212 {
1213   // Return first (index=1) or second (index=2) secondary particle of primary particle p 
1214
1215   if(index <= 0) 
1216     return 0 ;
1217   if(index > 2)
1218     return 0 ;
1219
1220   if(p) {
1221   Int_t daughterIndex = p->GetDaughter(index-1) ; 
1222   return  gAlice->Particle(daughterIndex) ; 
1223   }
1224   else
1225     return 0 ;
1226 }
1227
1228 //____________________________________________________________________________ 
1229 void AliPHOSGetter::ReadTreeD()
1230 {
1231   // Read the digit tree gAlice->TreeD()  
1232   if(gAlice->TreeD()== 0){
1233     cerr <<   "ERROR: AliPHOSGetter::ReadTreeD: can not read TreeD " << endl ;
1234   return ;
1235   }
1236   
1237   TObjArray * lob = static_cast<TObjArray*>(gAlice->TreeD()->GetListOfBranches()) ;
1238   TIter next(lob) ; 
1239   TBranch * branch = 0 ; 
1240   TBranch * digitsbranch = 0 ; 
1241   TBranch * digitizerbranch = 0 ; 
1242   Bool_t phosfound = kFALSE, digitizerfound = kFALSE ; 
1243   
1244   while ( (branch = static_cast<TBranch*>(next())) && (!phosfound || !digitizerfound) ) {
1245     if ( (strcmp(branch->GetName(), "PHOS")==0) && (strcmp(branch->GetTitle(), fDigitsTitle)==0) ) {
1246       digitsbranch = branch ; 
1247       phosfound = kTRUE ;
1248     }
1249     else if ( (strcmp(branch->GetName(), "AliPHOSDigitizer")==0) && (strcmp(branch->GetTitle(), fDigitsTitle)==0) ) {
1250       digitizerbranch = branch ; 
1251       digitizerfound = kTRUE ; 
1252     }
1253   }
1254
1255   if ( !phosfound || !digitizerfound ) {
1256     if (fDebug)
1257       cout << "WARNING: AliPHOSGetter::ReadTreeD -> Cannot find Digits and/or Digitizer with name " 
1258            << fDigitsTitle << endl ;
1259     return ; 
1260   }   
1261  
1262   //read digits
1263   if(!Digits(fDigitsTitle) ) 
1264     PostDigits(fDigitsTitle);
1265   digitsbranch->SetAddress(DigitsRef(fDigitsTitle)) ;
1266   digitsbranch->GetEntry(0) ;
1267   
1268   
1269   // read  the Digitizer
1270   if(!Digitizer(fDigitsTitle))
1271     PostDigitizer(fDigitsTitle) ;
1272   digitizerbranch->SetAddress(DigitizerRef(fDigitsTitle)) ;
1273   digitizerbranch->GetEntry(0) ;
1274  
1275   
1276 }
1277
1278 //____________________________________________________________________________ 
1279 void AliPHOSGetter::ReadTreeH()
1280 {
1281   // Read the first entry of PHOS branch in hit tree gAlice->TreeH()
1282
1283   if(gAlice->TreeH()== 0){
1284     cerr <<   "ERROR: AliPHOSGetter::ReadTreeH: -> Cannot read TreeH " << endl ;
1285     return ;
1286   }
1287   
1288   TBranch * hitsbranch = static_cast<TBranch*>(gAlice->TreeH()->GetBranch("PHOS")) ;
1289   if ( !hitsbranch ) {
1290     if (fDebug)
1291       cout << "WARNING:  AliPHOSGetter::ReadTreeH -> Cannot find branch PHOS" << endl ; 
1292     return ;
1293   }
1294   if(!Hits())
1295     PostHits() ;
1296
1297   if (hitsbranch->GetEntries() > 1 ) {
1298     TClonesArray * tempo =  new TClonesArray("AliPHOSHit",1000) ;
1299     TClonesArray * hits = dynamic_cast<TClonesArray*>(*HitsRef()) ; 
1300     hitsbranch->SetAddress(&tempo) ;
1301     Int_t index = 0 ; 
1302     Int_t i = 0 ;
1303     for (i = 0 ; i < hitsbranch->GetEntries() ; i++) {
1304       hitsbranch->GetEntry(i) ;
1305       Int_t j = 0 ; 
1306       for ( j = 0 ; j < tempo->GetEntries() ; j++) { 
1307         const AliPHOSHit * hit = static_cast<const AliPHOSHit *>(tempo->At(j)) ; 
1308         new((*hits)[index]) AliPHOSHit( *hit ) ;
1309         index++ ; 
1310       }
1311     }
1312     delete tempo ; 
1313   }
1314   else {
1315     hitsbranch->SetAddress(HitsRef()) ;
1316     hitsbranch->GetEntry(0) ;
1317   }
1318 }
1319
1320 //____________________________________________________________________________ 
1321 void AliPHOSGetter::Track(Int_t itrack)
1322 {
1323   // Read the first entry of PHOS branch in hit tree gAlice->TreeH()
1324
1325   if(gAlice->TreeH()== 0){
1326     cerr <<   "ERROR: AliPHOSGetter::ReadTreeH: -> Cannot read TreeH " << endl ;
1327     return ;
1328   }
1329   
1330   TBranch * hitsbranch = dynamic_cast<TBranch*>(gAlice->TreeH()->GetListOfBranches()->FindObject("PHOS")) ;
1331   if ( !hitsbranch ) {
1332     if (fDebug)
1333       cout << "WARNING:  AliPHOSGetter::ReadTreeH -> Cannot find branch PHOS" << endl ; 
1334     return ;
1335   }  
1336   if(!Hits())
1337     PostHits() ;
1338
1339   hitsbranch->SetAddress(HitsRef()) ;
1340   hitsbranch->GetEntry(itrack) ;
1341
1342
1343 }
1344 //____________________________________________________________________________ 
1345 void AliPHOSGetter::ReadTreeQA()
1346 {
1347   // Read the digit tree gAlice->TreeQA()
1348   // so far only PHOS knows about this Tree  
1349
1350   if(PHOS()->TreeQA()== 0){
1351     cerr <<   "ERROR: AliPHOSGetter::ReadTreeQA: can not read TreeQA " << endl ;
1352     return ;
1353   }
1354   
1355   TBranch * qabranch = PHOS()->TreeQA()->GetBranch("PHOS") ; 
1356   if (!qabranch) { 
1357     if (fDebug)
1358       cout << "WARNING: AliPHOSGetter::ReadTreeQA -> Cannot find QA Alarms for PHOS" << endl ;
1359     return ; 
1360   }   
1361   
1362   if(!Alarms())
1363     PostQA() ; 
1364
1365   qabranch->SetAddress(AlarmsRef()) ;
1366
1367   qabranch->GetEntry(0) ;
1368  
1369 //   PostQA("PHOS") ; 
1370 //   TFolder * alarmsF = Alarms() ; 
1371 //   alarmsF->Clear() ; 
1372 //   qabranch->SetAddress(&alarmsF) ;
1373 //   qabranch->GetEntry(0) ;
1374   
1375 }
1376
1377 //____________________________________________________________________________ 
1378 void AliPHOSGetter::ReadTreeR()
1379 {
1380   // Read the reconstrunction tree gAlice->TreeR()
1381
1382   if(gAlice->TreeR()== 0){
1383     cerr <<   "ERROR: AliPHOSGetter::ReadTreeR: can not read TreeR " << endl ;
1384     return ;
1385   }
1386   
1387   // RecPoints 
1388   TObjArray * lob = static_cast<TObjArray*>(gAlice->TreeR()->GetListOfBranches()) ;
1389   TIter next(lob) ; 
1390   TBranch * branch = 0 ; 
1391   TBranch * emcbranch = 0 ; 
1392   TBranch * cpvbranch = 0 ; 
1393   TBranch * clusterizerbranch = 0 ; 
1394   Bool_t phosemcrpfound = kFALSE, phoscpvrpfound = kFALSE, clusterizerfound = kFALSE ; 
1395   
1396   while ( (branch = static_cast<TBranch*>(next())) && (!phosemcrpfound || !phoscpvrpfound || !clusterizerfound) ) 
1397     if(strcmp(branch->GetTitle(), fRecPointsTitle)==0) {
1398       if ( strcmp(branch->GetName(), "PHOSEmcRP")==0) {
1399         emcbranch = branch ; 
1400         phosemcrpfound = kTRUE ;
1401       }
1402       else if ( strcmp(branch->GetName(), "PHOSCpvRP")==0) {
1403         cpvbranch = branch ; 
1404         phoscpvrpfound = kTRUE ;
1405       }
1406       else if(strcmp(branch->GetName(), "AliPHOSClusterizer")==0){
1407         clusterizerbranch = branch ; 
1408         clusterizerfound = kTRUE ; 
1409       }
1410     }
1411
1412   if ( !phosemcrpfound ) {
1413     if (fDebug)
1414       cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find EmcRecPoints with title " 
1415            << fRecPointsTitle << endl ;
1416     return ; 
1417   }   
1418   if ( !phoscpvrpfound ) {
1419     if (fDebug)
1420       cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find CpvRecPoints with title " 
1421            << fRecPointsTitle << endl ;
1422     return ; 
1423   }   
1424   if ( !clusterizerfound ) {
1425     if (fDebug)
1426       cout << "WARNING: AliPHOSGetter::ReadTreeR -> Can not find Clusterizer with title " 
1427            << fRecPointsTitle << endl ;
1428     return ; 
1429   }   
1430   
1431   // Read and Post the RecPoints
1432   if(!EmcRecPoints(fRecPointsTitle) ) 
1433     PostRecPoints(fRecPointsTitle) ;
1434
1435   emcbranch->SetAddress(EmcRecPointsRef(fRecPointsTitle)) ;
1436   emcbranch->GetEntry(0) ;
1437
1438   cpvbranch->SetAddress(CpvRecPointsRef(fRecPointsTitle)) ; 
1439   cpvbranch->GetEntry(0) ;
1440   
1441   if(!Clusterizer(fRecPointsTitle) )
1442     PostClusterizer(fRecPointsTitle) ;
1443
1444   clusterizerbranch->SetAddress(ClusterizerRef(fRecPointsTitle)) ;
1445   clusterizerbranch->GetEntry(0) ;
1446  
1447   
1448   //------------------- TrackSegments ---------------------
1449   next.Reset() ; 
1450   TBranch * tsbranch = 0 ; 
1451   TBranch * tsmakerbranch = 0 ; 
1452   Bool_t phostsfound = kFALSE, tsmakerfound = kFALSE ; 
1453     
1454   while ( (branch = static_cast<TBranch*>(next())) && (!phostsfound || !tsmakerfound) ) 
1455     if(strcmp(branch->GetTitle(), fTrackSegmentsTitle)==0)  {
1456       if ( strcmp(branch->GetName(), "PHOSTS")==0){
1457         tsbranch = branch ; 
1458         phostsfound = kTRUE ;
1459       }
1460       else if(strcmp(branch->GetName(), "AliPHOSTrackSegmentMaker")==0) {
1461         tsmakerbranch = branch ; 
1462         tsmakerfound  = kTRUE ; 
1463       }
1464     }
1465   
1466   if ( !phostsfound || !tsmakerfound ) {
1467     if (fDebug)
1468       cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find TrackSegments and/or TrackSegmentMaker with name "
1469            << fTrackSegmentsTitle << endl ;
1470     return ; 
1471   } 
1472   
1473   // Read and Post the TrackSegments
1474   if(!TrackSegments(fTrackSegmentsTitle))
1475     PostTrackSegments(fTrackSegmentsTitle) ;
1476
1477   tsbranch->SetAddress(TrackSegmentsRef(fTrackSegmentsTitle)) ;
1478   tsbranch->GetEntry(0) ;
1479   
1480   // Read and Post the TrackSegment Maker
1481   if(!TrackSegmentMaker(fTrackSegmentsTitle))
1482     PostTrackSegmentMaker(fTrackSegmentsTitle) ;
1483  
1484   tsmakerbranch->SetAddress(TSMakerRef(fTrackSegmentsTitle)) ;
1485   tsmakerbranch->GetEntry(0) ;
1486   
1487   
1488   //------------ RecParticles ----------------------------
1489   next.Reset() ; 
1490   TBranch * rpabranch = 0 ; 
1491   TBranch * pidbranch = 0 ; 
1492   Bool_t phosrpafound = kFALSE, pidfound = kFALSE ; 
1493   
1494   while ( (branch = static_cast<TBranch*>(next())) && (!phosrpafound || !pidfound) ) 
1495     if(strcmp(branch->GetTitle(), fRecParticlesTitle)==0) {   
1496       if ( strcmp(branch->GetName(), "PHOSRP")==0) {   
1497         rpabranch = branch ; 
1498         phosrpafound = kTRUE ;
1499       }
1500       else if (strcmp(branch->GetName(), "AliPHOSPID")==0) {
1501         pidbranch = branch ; 
1502         pidfound  = kTRUE ; 
1503       }
1504     }
1505   
1506   if ( !phosrpafound || !pidfound ) {
1507     if (fDebug)
1508       cout << "WARNING: AliPHOSGetter::ReadTreeR -> Cannot find RecParticles and/or PID with name " 
1509            << fRecParticlesTitle << endl ;
1510     return ; 
1511   } 
1512   
1513   // Read and Post the RecParticles
1514   if(!RecParticles(fRecParticlesTitle))
1515     PostRecParticles(fRecParticlesTitle) ;
1516
1517   rpabranch->SetAddress(RecParticlesRef(fRecParticlesTitle)) ;
1518   rpabranch->GetEntry(0) ;
1519   
1520   // Read and Post the PID
1521   if(!PID(fRecParticlesTitle))
1522     PostPID(fRecParticlesTitle) ;
1523
1524   pidbranch->SetAddress(PIDRef(fRecParticlesTitle)) ;
1525   pidbranch->GetEntry(0) ;
1526   
1527   
1528 }
1529
1530 //____________________________________________________________________________ 
1531 void AliPHOSGetter::ReadTreeS(Int_t event)
1532 {
1533   // Read the summable digits tree gAlice->TreeS()  
1534   
1535   // loop over all opened files and read their SDigits to the White Board
1536   TFolder * phosF = dynamic_cast<TFolder *>(fSDigitsFolder->FindObject("PHOS")) ;
1537   if (!phosF) 
1538     phosF = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; 
1539   TCollection * folderslist = phosF->GetListOfFolders() ; 
1540   
1541   //Add current file to list if it is not there yet
1542   
1543   TString subdir(fHeaderFile) ;
1544   subdir.ReplaceAll("/","_") ; 
1545
1546   if ( (subdir != "aliroot") && ( !folderslist->Contains(subdir) ) ){
1547     phosF->AddFolder(subdir, ""); 
1548   }
1549     
1550   TIter next(folderslist) ; 
1551   TFolder * folder = 0 ; 
1552   TFile * file; 
1553   TTree * treeS = 0;
1554   while ( (folder = static_cast<TFolder*>(next())) ) {
1555     TString fileName(folder->GetName()) ; 
1556     fileName.ReplaceAll("_","/") ; 
1557     if(fHeaderFile.CompareTo(fileName) == 0 ) 
1558       treeS=gAlice->TreeS() ;
1559     else{
1560       file = static_cast<TFile*>(gROOT->GetFile(fileName)); 
1561       file->cd() ;
1562       
1563       // Get SDigits Tree header from file
1564       TString treeName("TreeS") ;
1565       treeName += event ; 
1566       treeS = dynamic_cast<TTree*>(gDirectory->Get(treeName.Data()));
1567     }
1568     if(treeS==0){
1569       cerr << "ERROR: AliPHOSGetter::ReadTreeS There is no SDigit Tree" << endl;
1570       return ;
1571     }
1572     
1573     //set address of the SDigits and SDigitizer
1574     TBranch   * sdigitsBranch    = 0;
1575     TBranch   * sdigitizerBranch = 0;
1576     TBranch   * branch           = 0 ;  
1577     TObjArray * lob = static_cast<TObjArray*>(treeS->GetListOfBranches()) ;
1578     TIter next(lob) ; 
1579     Bool_t phosfound = kFALSE, sdigitizerfound = kFALSE ; 
1580     
1581     while ( (branch = static_cast<TBranch*>(next())) && (!phosfound || !sdigitizerfound) ) {
1582       if ( (strcmp(branch->GetName(), "PHOS")==0) && (strcmp(branch->GetTitle(), fSDigitsTitle)==0) ) {
1583         phosfound = kTRUE ;
1584         sdigitsBranch = branch ; 
1585       }
1586       
1587       else if ( (strcmp(branch->GetName(), "AliPHOSSDigitizer")==0) && (strcmp(branch->GetTitle(), fSDigitsTitle)==0) ) {
1588         sdigitizerfound = kTRUE ; 
1589         sdigitizerBranch = branch ;
1590       }
1591     }
1592     if ( !phosfound || !sdigitizerfound ) {
1593       if (fDebug)
1594         cout << "WARNING: AliPHOSDigitizer::ReadSDigits -> Digits and/or Digitizer branch with name " << GetName() 
1595              << " not found" << endl ;
1596       return ; 
1597     }   
1598     
1599     if ( !folder->FindObject(fSDigitsTitle) )  
1600       PostSDigits(fSDigitsTitle,folder->GetName()) ;
1601
1602     sdigitsBranch->SetAddress(SDigitsRef(fSDigitsTitle,folder->GetName())) ;
1603     sdigitsBranch->GetEntry(0) ;
1604     
1605     TString sdname(fSDigitsTitle) ;
1606     sdname+=":" ;
1607     sdname+=folder->GetName() ;
1608     if(!SDigitizer(sdname) ) 
1609       PostSDigitizer(fSDigitsTitle,folder->GetName()) ;
1610
1611     sdigitizerBranch->SetAddress(SDigitizerRef(sdname)) ;
1612     sdigitizerBranch->GetEntry(0) ;
1613     
1614   }    
1615   
1616   // After SDigits have been read from all files, return to the first one
1617   
1618   next.Reset();
1619   folder = static_cast<TFolder*>(next());
1620   if(folder){
1621     TString fileName(folder->GetName()) ; 
1622     fileName.ReplaceAll("_","/") ; 
1623     file   = static_cast<TFile*>(gROOT->GetFile(fileName)); 
1624     file   ->cd() ;
1625   }
1626   
1627 }
1628 //____________________________________________________________________________ 
1629 void AliPHOSGetter::ReadTreeS(TTree * treeS, Int_t input)
1630 {  // Read the summable digits fron treeS()  
1631
1632
1633   TString filename("mergefile") ;
1634   filename+= input ;
1635
1636   TFolder * phosFolder = dynamic_cast<TFolder*>(fSDigitsFolder->FindObject("PHOS")) ; 
1637   if ( !phosFolder ) { 
1638    phosFolder = fSDigitsFolder->AddFolder("PHOS", "SDigits from PHOS") ; 
1639   } 
1640   TFolder * folder=(TFolder*)phosFolder->FindObject(filename) ;
1641   //set address of the SDigits and SDigitizer
1642   TBranch   * sdigitsBranch    = 0;
1643   TBranch   * sdigitizerBranch = 0;
1644   TBranch   * branch           = 0 ;  
1645   TObjArray * lob = (TObjArray*)treeS->GetListOfBranches() ;
1646   TIter next(lob) ; 
1647   Bool_t phosfound = kFALSE, sdigitizerfound = kFALSE ; 
1648   
1649   while ( (branch = (TBranch*)next()) && (!phosfound || !sdigitizerfound) ) {
1650     if ( strcmp(branch->GetName(), "PHOS")==0) {
1651       phosfound = kTRUE ;
1652       sdigitsBranch = branch ; 
1653     }
1654     
1655     else if ( strcmp(branch->GetName(), "AliPHOSSDigitizer")==0) {
1656       sdigitizerfound = kTRUE ; 
1657       sdigitizerBranch = branch ;
1658     }
1659   }
1660   if ( !phosfound || !sdigitizerfound ) {
1661     if (fDebug)
1662       cout << "WARNING: AliPHOSGetter::ReadTreeS -> Digits and/or Digitizer branch not found" << endl ;
1663     return ; 
1664   }   
1665   
1666   if (!folder || !(folder->FindObject(sdigitsBranch->GetTitle()) ) )
1667     PostSDigits(sdigitsBranch->GetTitle(),filename) ;
1668
1669   sdigitsBranch->SetAddress(SDigitsRef(sdigitsBranch->GetTitle(),filename)) ;
1670   sdigitsBranch->GetEntry(0) ;
1671   
1672   TString sdname(sdigitsBranch->GetTitle()) ;
1673   sdname+=":" ;
1674   sdname+=filename ;
1675   if(!SDigitizer(sdigitsBranch->GetTitle()) )
1676     PostSDigitizer(sdigitsBranch->GetTitle(),filename) ;
1677
1678   sdigitizerBranch->SetAddress(SDigitizerRef(sdname)) ;
1679   sdigitizerBranch->GetEntry(0) ;
1680   
1681 }    
1682
1683
1684 //____________________________________________________________________________ 
1685 void AliPHOSGetter::ReadPrimaries()
1686 {
1687
1688   
1689   // Read the first entry of PHOS branch in hit tree gAlice->TreeH()
1690   TTree * particleTree = gAlice->TreeK()  ; 
1691   if( ! particleTree ){
1692     cerr <<   "ERROR: AliPHOSGetter::ReadTreePrimaries: -> Cannot read TreeK " << endl ;
1693     return ;
1694   }
1695   
1696   TBranch * particlesbranch = static_cast<TBranch*>(particleTree->GetBranch("Particles")) ;
1697   if ( !particlesbranch ) {
1698     if (fDebug)
1699       cout << "WARNING:  AliPHOSGetter::ReadTreePrimaries -> Cannot find branch Particles" << endl ; 
1700     return ;
1701   }
1702
1703   TParticle * particle = 0; 
1704   TClonesArray * ar ; 
1705   if(! (ar = Primaries()) ) { 
1706     PostPrimaries() ;
1707     ar = Primaries() ; 
1708   }
1709   ar->Delete() ; 
1710   particlesbranch->SetAddress(&particle) ;
1711   
1712   Int_t index = 0 ; 
1713   for (index = 0 ; index < particleTree->GetEntries(); index++) { 
1714     particlesbranch->GetEntry(index) ;
1715     new ((*ar)[index]) TParticle(*particle);
1716   }
1717   fNPrimaries= ar->GetEntries() ; 
1718   
1719   // Reads specific branches of primaries
1720   // fNPrimaries = gAlice->GetNtrack();
1721   
1722   //Check, is it necessary to open new files
1723  //  TArrayI* events = fDigitizer->GetCurrentEvents() ; 
1724 //   TClonesArray * filenames = fDigitizer->GetHeadersFiles() ;
1725 //   Int_t input ;
1726 //   for(input = 0; input < filenames->GetEntriesFast(); input++){
1727     
1728 //     TObjString * filename = (TObjString *) filenames->At(input) ;
1729
1730 //     //Test, if this file already open
1731 //     TFile *file = (TFile*) gROOT->GetFile( filename->GetString() ) ;
1732 //     if(file == 0)
1733 //       file = new TFile( filename->GetString()) ;
1734 //     file->cd() ;
1735     
1736 //     // Get Kine Tree from file
1737 //     char treeName[20];
1738 //     sprintf(treeName,"TreeK%d",events->At(input));
1739 //     TTree * treeK = (TTree*)gDirectory->Get(treeName);
1740 //     if (treeK) 
1741 //       treeK->SetBranchAddress("Particles", &fParticleBuffer);
1742 //     else    
1743 //       cout << "AliPHOSGetter: cannot find Kine Tree for event:" << events->At(input) << endl;
1744
1745 //     // Create the particle stack
1746 //     if(!fParticles) fParticles = new TClonesArray("TParticle",1000);
1747 //     // Build the pointer list
1748 //     if(fParticleMap) {     <----
1749 //       fParticleMap->Clear();
1750 //       fParticleMap->Expand(treeK->GetEntries());
1751 //     } else
1752 //       fParticleMap = new TObjArray(treeK->GetEntries());
1753     
1754 //     // From gAlice->Particle(i) 
1755
1756
1757 //   if(!(*fParticleMap)[i]) {
1758 //     Int_t nentries = fParticles->GetEntries();
1759     
1760 //     // algorithmic way of getting entry index
1761 //     // (primary particles are filled after secondaries)
1762 //     Int_t entry;
1763 //     if (i<fHeader.GetNprimary())
1764 //       entry = i+fHeader.GetNsecondary();
1765 //     else 
1766 //       entry = i-fHeader.GetNprimary();
1767       
1768 //     // only check the algorithmic way and give
1769 //     // the fatal error if it is wrong
1770 //     if (entry != fParticleFileMap[i]) {
1771 //       Fatal("Particle",
1772 //         "!!!! The algorithmic way is WRONG: !!!\n entry: %d map: %d",
1773 //      entry, fParticleFileMap[i]); 
1774 //     }  
1775       
1776 //     fTreeK->GetEntry(fParticleFileMap[i]);
1777 //     new ((*fParticles)[nentries]) TParticle(*fParticleBuffer);
1778 //     fParticleMap->AddAt((*fParticles)[nentries],i);
1779 //   }
1780 //   return (TParticle *) (*fParticleMap)[i];
1781
1782    
1783     
1784 // //   }
1785
1786
1787 // //   //scan over opened files and read corresponding TreeK##
1788
1789 //   return ;
1790 }
1791 //____________________________________________________________________________ 
1792 void AliPHOSGetter::Event(const Int_t event, const char* opt)
1793 {
1794   // Reads the content of all Tree's S, D and R
1795   
1796   if (event >= gAlice->TreeE()->GetEntries() ) {
1797     cerr << "ERROR: AliPHOSGetter::Event -> " << event << " not found in TreeE!" << endl ; 
1798     return ; 
1799   }
1800   gAlice->GetEvent(event) ; 
1801
1802   if(strstr(opt,"H") )
1803     ReadTreeH() ;
1804   
1805   if(strstr(opt,"S") )
1806     ReadTreeS(event) ;
1807
1808   if( strstr(opt,"D") )
1809     ReadTreeD() ;
1810
1811   if( strstr(opt,"R") )
1812     ReadTreeR() ;
1813
1814   if( strstr(opt,"Q") )
1815     ReadTreeQA() ;
1816
1817   //  if( strstr(opt,"P") || (strcmp(opt,"")==0) )
1818   //  ReadPrimaries() ;
1819
1820 }
1821
1822 //____________________________________________________________________________ 
1823 TObject * AliPHOSGetter::ReturnO(TString what, TString name, TString file) const 
1824 {
1825   // get the object named "what" from the folder
1826   // folders are named like //Folders
1827
1828   if ( file.IsNull() ) 
1829     file = fHeaderFile ; 
1830
1831   TFolder * folder = 0 ;
1832   TObject * phosO  = 0 ; 
1833
1834   //  if ( name.IsNull() ) {
1835   if ( what.CompareTo("Primaries") == 0 ) {
1836     folder = dynamic_cast<TFolder *>(fPrimariesFolder->FindObject("Primaries")) ; 
1837     if (folder) 
1838       phosO  = dynamic_cast<TObject *>(folder->FindObject("Primaries")) ;  
1839   }
1840   else if ( what.CompareTo("Hits") == 0 ) {
1841     folder = dynamic_cast<TFolder *>(fHitsFolder->FindObject("PHOS")) ; 
1842     if (folder) 
1843       phosO  = dynamic_cast<TObject *>(folder->FindObject("Hits")) ;  
1844   }
1845   else if ( what.CompareTo("SDigits") == 0 ) {
1846     file.ReplaceAll("/","_") ; 
1847     TString path = "PHOS/" + file  ;
1848     folder = dynamic_cast<TFolder *>(fSDigitsFolder->FindObject(path.Data())) ; 
1849     if (folder) { 
1850       if (name.IsNull())
1851         name = fSDigitsTitle ; 
1852       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1853     }
1854   }
1855   else if ( what.CompareTo("Digits") == 0 ){
1856     folder = dynamic_cast<TFolder *>(fDigitsFolder->FindObject("PHOS")) ; 
1857     if (folder) { 
1858       if (name.IsNull())
1859         name = fDigitsTitle ; 
1860       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1861     } 
1862   }
1863   else if ( what.CompareTo("EmcRecPoints") == 0 ) {
1864     folder = dynamic_cast<TFolder *>(fRecoFolder->FindObject("PHOS/EMCARecPoints")) ; 
1865     if (folder) { 
1866       if (name.IsNull())
1867         name = fRecPointsTitle ; 
1868       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1869     } 
1870   }
1871   else if ( what.CompareTo("CpvRecPoints") == 0 ) {
1872     folder = dynamic_cast<TFolder *>(fRecoFolder->FindObject("PHOS/CPVRecPoints")) ; 
1873     if (folder) { 
1874       if (name.IsNull())
1875         name = fRecPointsTitle ; 
1876       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1877     }   
1878   }
1879   else if ( what.CompareTo("TrackSegments") == 0 ) {
1880     folder = dynamic_cast<TFolder *>(fRecoFolder->FindObject("PHOS/TrackSegments")) ; 
1881     if (folder) { 
1882       if (name.IsNull())
1883         name = fTrackSegmentsTitle ; 
1884       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1885     }   
1886   }
1887   else if ( what.CompareTo("RecParticles") == 0 ) {
1888     folder = dynamic_cast<TFolder *>(fRecoFolder->FindObject("PHOS/RecParticles")) ; 
1889    if (folder) { 
1890       if (name.IsNull())
1891         name = fRecParticlesTitle ; 
1892       phosO  = dynamic_cast<TObject *>(folder->FindObject(name)) ; 
1893     }   
1894  }
1895   else if ( what.CompareTo("Alarms") == 0 ){ 
1896     if (name.IsNull() ) 
1897       phosO = dynamic_cast<TObject *>(fQAFolder->FindObject("PHOS")) ;  
1898     else {
1899       folder = dynamic_cast<TFolder *>(fQAFolder->FindObject("PHOS")) ; 
1900       if (!folder) 
1901         phosO = 0 ; 
1902       else 
1903         phosO = dynamic_cast<TObject *>(folder->FindObject(name)) ;  
1904     }
1905   }
1906   if (!phosO) {
1907     if(fDebug)
1908       cerr << "ERROR : AliPHOSGetter::ReturnO -> Object " << what << " not found in " << folder->GetName() << endl ; 
1909     return 0 ;
1910   }
1911
1912   return phosO ;
1913 }
1914   
1915 //____________________________________________________________________________ 
1916 const TTask * AliPHOSGetter::ReturnT(TString what, TString name) const 
1917 {
1918   // get the TTask named "what" from the folder
1919   // folders are named like //Folders/Tasks/what/PHOS/name
1920
1921   TString search(what) ; 
1922   if ( what.CompareTo("Clusterizer") == 0 ) 
1923     search = "Reconstructioner" ; 
1924   else if ( what.CompareTo("TrackSegmentMaker") == 0 ) 
1925     search = "Reconstructioner" ; 
1926   else if ( what.CompareTo("PID") == 0 ) 
1927     search = "Reconstructioner" ; 
1928   else if ( what.CompareTo("QATasks") == 0 ) 
1929     search = "QA" ; 
1930
1931   TTask * tasks = dynamic_cast<TTask*>(fTasksFolder->FindObject(search)) ; 
1932
1933   if (!tasks) {
1934     cerr << "ERROR: AliPHOSGetter::ReturnT -> Task " << what << " not found!" << endl ;  
1935     return 0 ; 
1936   }
1937
1938   TTask * phosT = dynamic_cast<TTask*>(tasks->GetListOfTasks()->FindObject("PHOS")) ; 
1939   if (!phosT) { 
1940     cerr << "ERROR: AliPHOSGetter::ReturnT -> Task " << what << "/PHOS not found!" << endl ;  
1941     return 0 ; 
1942   }
1943   
1944   TList * list = phosT->GetListOfTasks() ; 
1945  
1946   if (what.CompareTo("SDigitizer") == 0) {  
1947     if ( name.IsNull() )
1948       name =  fSDigitsTitle ; 
1949   } else  if (what.CompareTo("Digitizer") == 0){ 
1950     if ( name.IsNull() )
1951       name =  fDigitsTitle ;
1952   } else  if (what.CompareTo("Clusterizer") == 0){ 
1953     if ( name.IsNull() )
1954       name =  fRecPointsTitle ;
1955     name.Append(":clu") ;
1956   }
1957   else  if (what.CompareTo("TrackSegmentMaker") == 0){ 
1958     if ( name.IsNull() )
1959       name =  fTrackSegmentsTitle ;
1960     name.Append(":tsm") ;
1961   }
1962   else  if (what.CompareTo("PID") == 0){ 
1963     if ( name.IsNull() )
1964       name =  fRecParticlesTitle ;
1965     name.Append(":pid") ;
1966   }
1967   else  if (what.CompareTo("QATasks") == 0){ 
1968     if ( name.IsNull() )
1969       return phosT ;
1970   }
1971   
1972   TIter it(list) ;
1973   TTask * task = 0 ; 
1974   while((task = static_cast<TTask *>(it.Next()) )){
1975     TString taskname(task->GetName()) ;
1976     if(taskname.BeginsWith(name))
1977       return task ;
1978   }
1979   
1980   if(fDebug)
1981     cout << "WARNING: AliPHOSGetter::ReturnT -> Task " << search << "/" << name << " not found!" << endl ; 
1982   return 0 ;
1983 }