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