Fixes for bug #52499: Field polarities inconsistiency
[u/mrichter/AliRoot.git] / STEER / AliRunLoader.cxx
CommitLineData
8364b0ef 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
504b172d 18//____________________________________________________________________
19//////////////////////////////////////////////////////////////////////
20// //
21// class AliRunLoader //
22// //
23// This class aims to be the unque interface for managing data I/O. //
24// It stores Loaders for all modules which, knows names //
25// of the files were data are to be stored. //
26// //
27// It aims to substitud AliRun in automatic data managing //
28// thus there is no necessity of loading gAlice from file in order //
29// to get access to the data. //
30// //
31// Logical place to put the specific Loader to the given //
32// detector is detector itself (i.e ITSLoader in ITS). //
33// But, to load detector object one need to load gAlice, and //
34// by the way all other detectors with their geometrieces and //
35// so on. So, if one need to open TPC clusters there is no //
36// principal need to read everything. //
37// //
38// //
39// When RunLoader is read from the file it does not connect to //
40// the folder structure automatically. It must be connected //
41// (mounted) manualy. Default event folder is defined by //
e191bb57 42// AliConfig::GetDefaultEventFolderName() //
504b172d 43// but can be mounted elsewhere. Usefull specially in merging case, //
44// when more than pone session needs to be loaded //
45// //
46//////////////////////////////////////////////////////////////////////
88cb7938 47
134d2ab3 48#include <TROOT.h>
4095d1ca 49#include <TBranch.h>
4095d1ca 50#include <TFile.h>
51#include <TFolder.h>
4095d1ca 52#include <TObjArray.h>
88cb7938 53#include <TString.h>
024a7e64 54class TTask;
88cb7938 55#include <TTree.h>
88cb7938 56
21bf7095 57#include "AliLog.h"
88cb7938 58#include "AliRun.h"
59#include "AliConfig.h"
60#include "AliLoader.h"
61#include "AliHeader.h"
62#include "AliStack.h"
88cb7938 63#include "AliDetector.h"
9e1ceb13 64#include "AliCDBManager.h"
fe913d8f 65#include "AliCDBLocal.h"
bacbe0fd 66#include "AliCentralTrigger.h"
88cb7938 67
68ClassImp(AliRunLoader)
69
70AliRunLoader* AliRunLoader::fgRunLoader = 0x0;
71
72const TString AliRunLoader::fgkRunLoaderName("RunLoader");
88cb7938 73const TString AliRunLoader::fgkHeaderBranchName("Header");
bacbe0fd 74const TString AliRunLoader::fgkTriggerBranchName("ClassMask");
88cb7938 75const TString AliRunLoader::fgkHeaderContainerName("TE");
bacbe0fd 76const TString AliRunLoader::fgkTriggerContainerName("TreeCT");
88cb7938 77const TString AliRunLoader::fgkKineContainerName("TreeK");
78const TString AliRunLoader::fgkTrackRefsContainerName("TreeTR");
79const TString AliRunLoader::fgkKineBranchName("Particles");
80const TString AliRunLoader::fgkDefaultKineFileName("Kinematics.root");
81const TString AliRunLoader::fgkDefaultTrackRefsFileName("TrackRefs.root");
82const TString AliRunLoader::fgkGAliceName("gAlice");
bacbe0fd 83const TString AliRunLoader::fgkDefaultTriggerFileName("Trigger.root");
88cb7938 84/**************************************************************************/
85
86AliRunLoader::AliRunLoader():
87 fLoaders(0x0),
88 fEventFolder(0x0),
33c3c91a 89 fRun(-1),
88cb7938 90 fCurrentEvent(0),
91 fGAFile(0x0),
92 fHeader(0x0),
93 fStack(0x0),
bacbe0fd 94 fCTrigger(0x0),
88cb7938 95 fKineDataLoader(0x0),
96 fTrackRefsDataLoader(0x0),
97 fNEventsPerFile(1),
c516e34c 98 fNEventsPerRun(0),
88cb7938 99 fUnixDirName(".")
100{
101 AliConfig::Instance();//force to build the folder structure
85a5290f 102 if (!fgRunLoader) fgRunLoader = this;
88cb7938 103}
104/**************************************************************************/
105
106AliRunLoader::AliRunLoader(const char* eventfoldername):
107 TNamed(fgkRunLoaderName,fgkRunLoaderName),
108 fLoaders(new TObjArray()),
109 fEventFolder(0x0),
33c3c91a 110 fRun(-1),
88cb7938 111 fCurrentEvent(0),
112 fGAFile(0x0),
113 fHeader(0x0),
114 fStack(0x0),
bacbe0fd 115 fCTrigger(0x0),
88cb7938 116 fKineDataLoader(new AliDataLoader(fgkDefaultKineFileName,fgkKineContainerName,"Kinematics")),
117 fTrackRefsDataLoader(new AliDataLoader(fgkDefaultTrackRefsFileName,fgkTrackRefsContainerName,"Track References")),
118 fNEventsPerFile(1),
c516e34c 119 fNEventsPerRun(0),
88cb7938 120 fUnixDirName(".")
121{
122//ctor
123 SetEventFolderName(eventfoldername);
85a5290f 124 if (!fgRunLoader) fgRunLoader = this;
88cb7938 125}
126/**************************************************************************/
127
128AliRunLoader::~AliRunLoader()
129{
f2a509af 130//dtor
1bb20a37 131 if (fgRunLoader == this) fgRunLoader = 0x0;
132
88cb7938 133 UnloadHeader();
134 UnloadgAlice();
135
136 if(fLoaders) {
137 fLoaders->SetOwner();
138 delete fLoaders;
139 }
140
141 delete fKineDataLoader;
142 delete fTrackRefsDataLoader;
143
144
145 RemoveEventFolder();
146
147 //fEventFolder is deleted by the way of removing - TopAliceFolder owns it
bacbe0fd 148 if( fCTrigger ) delete fCTrigger;
88cb7938 149 delete fHeader;
150 delete fStack;
151 delete fGAFile;
152}
153/**************************************************************************/
154
85a5290f 155AliRunLoader::AliRunLoader(TFolder* topfolder):
156 TNamed(fgkRunLoaderName,fgkRunLoaderName),
157 fLoaders(new TObjArray()),
158 fEventFolder(topfolder),
8513c72c 159 fRun(-1),
85a5290f 160 fCurrentEvent(0),
161 fGAFile(0x0),
162 fHeader(0x0),
163 fStack(0x0),
bacbe0fd 164 fCTrigger(0x0),
85a5290f 165 fKineDataLoader(new AliDataLoader(fgkDefaultKineFileName,fgkKineContainerName,"Kinematics")),
166 fTrackRefsDataLoader(new AliDataLoader(fgkDefaultTrackRefsFileName,fgkTrackRefsContainerName,"Track References")),
167 fNEventsPerFile(1),
c516e34c 168 fNEventsPerRun(0),
85a5290f 169 fUnixDirName(".")
88cb7938 170{
f2a509af 171//ctor
88cb7938 172 if(topfolder == 0x0)
173 {
131d919a 174 TString errmsg("Parameter is NULL");
21bf7095 175 AliError(errmsg.Data());
131d919a 176 throw errmsg;
88cb7938 177 return;
178 }
88cb7938 179
180 TObject* obj = fEventFolder->FindObject(fgkRunLoaderName);
181 if (obj)
182 { //if it is, then sth. is going wrong... exits aliroot session
131d919a 183 TString errmsg("In Event Folder Named ");
184 errmsg+=fEventFolder->GetName();
185 errmsg+=" object named "+fgkRunLoaderName+" already exists. I am confused ...";
186
21bf7095 187 AliError(errmsg.Data());
131d919a 188 throw errmsg;
88cb7938 189 return;//never reached
190 }
131d919a 191
192 if (!fgRunLoader) fgRunLoader = this;
88cb7938 193
88cb7938 194 fEventFolder->Add(this);//put myself to the folder to accessible for all
195
196}
88cb7938 197
024a7e64 198/**************************************************************************/
199
88cb7938 200Int_t AliRunLoader::GetEvent(Int_t evno)
201{
202//Gets event number evno
203//Reloads all data properly
04096ea4 204//PH if (fCurrentEvent == evno) return 0;
88cb7938 205
206 if (evno < 0)
207 {
21bf7095 208 AliError("Can not give the event with negative number");
88cb7938 209 return 4;
210 }
211
212 if (evno >= GetNumberOfEvents())
213 {
21bf7095 214 AliError(Form("There is no event with number %d",evno));
88cb7938 215 return 3;
216 }
217
21bf7095 218 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
219 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
220 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
221 AliDebug(1, Form(" GETTING EVENT %d",evno));
222 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
223 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
224 AliDebug(1, ">>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>");
88cb7938 225
226 fCurrentEvent = evno;
227
228 Int_t retval;
229
230 //Reload header (If header was loaded)
231 if (GetHeader())
232 {
233 retval = TreeE()->GetEvent(fCurrentEvent);
234 if ( retval == 0)
235 {
21bf7095 236 AliError(Form("Cannot find event: %d\n ",fCurrentEvent));
88cb7938 237 return 5;
238 }
239 }
240 //Reload stack (If header was loaded)
241 if (TreeE()) fStack = GetHeader()->Stack();
242 //Set event folder in stack (it does not mean that we read kinematics from file)
bacbe0fd 243 if( GetTrigger() && TreeCT() ) {
244 retval = TreeCT()->GetEvent(fCurrentEvent);
c4e6edd3 245 if ( retval < 0 ) {
bacbe0fd 246 AliError(Form("Error occured while GetEvent for Trigger. Event %d",evno));
247 return 2;
248 }
249 }
88cb7938 250
251 retval = SetEvent();
252 if (retval)
253 {
21bf7095 254 AliError(Form("Error occured while setting event %d",evno));
88cb7938 255 return 1;
256 }
257
258 //Post Track References
259 retval = fTrackRefsDataLoader->GetEvent();
260 if (retval)
261 {
21bf7095 262 AliError(Form("Error occured while GetEvent for Track References. Event %d",evno));
88cb7938 263 return 2;
264 }
265
266 //Read Kinematics if loaded
b180b606 267 retval = fKineDataLoader->GetEvent();
88cb7938 268 if (retval)
269 {
21bf7095 270 AliError(Form("Error occured while GetEvent for Kinematics. Event %d",evno));
88cb7938 271 return 2;
272 }
273
a43369dc 274 if (fStack && fKineDataLoader->GetBaseLoader(0)->IsLoaded())
275 {
9ca424c1 276 fStack->ConnectTree(TreeK());
277
a43369dc 278 if (fStack->GetEvent() == kFALSE)
279 {
280 AliError(Form("Error occured while GetEvent for Stack. Event %d",evno));
281 return 2;
282 }
283 }
284
88cb7938 285 //Trigger data reloading in all loaders
286 TIter next(fLoaders);
287 AliLoader *loader;
288 while((loader = (AliLoader*)next()))
289 {
290 retval = loader->GetEvent();
291 if (retval)
292 {
21bf7095 293 AliError(Form("Error occured while getting event for %s. Event %d.",
294 loader->GetDetectorName().Data(), evno));
88cb7938 295 return 3;
296 }
297 }
298
299 SetDetectorAddresses();
300
301 return 0;
302}
303/**************************************************************************/
304Int_t AliRunLoader::SetEvent()
305{
f2a509af 306//if kinematics was loaded Cleans folder data
307
88cb7938 308 Int_t retval;
309
310 retval = fKineDataLoader->SetEvent();
311 if (retval)
312 {
21bf7095 313 AliError("SetEvent for Kinamtics Data Loader retutned error.");
88cb7938 314 return retval;
315 }
316 retval = fTrackRefsDataLoader->SetEvent();
317 if (retval)
318 {
21bf7095 319 AliError("SetEvent for Track References Data Loader retutned error.");
88cb7938 320 return retval;
321 }
322
323 TIter next(fLoaders);
324 AliLoader *loader;
325 while((loader = (AliLoader*)next()))
326 {
327 retval = loader->SetEvent();
328 if (retval)
329 {
21bf7095 330 AliError(Form("SetEvent for %s Data Loader retutned error.",loader->GetName()));
88cb7938 331 return retval;
332 }
333 }
334
335 return 0;
336}
337/**************************************************************************/
338
339Int_t AliRunLoader::SetEventNumber(Int_t evno)
340{
341 //cleans folders and sets the root dirs in files
342 if (fCurrentEvent == evno) return 0;
343 fCurrentEvent = evno;
344 return SetEvent();
345}
346
347/**************************************************************************/
9e1ceb13 348AliCDBEntry* AliRunLoader::GetCDBEntry(const char* name) const
2c8628dd 349{
9e1ceb13 350//Get an AliCDBEntry from the run data storage
2c8628dd 351
9e1ceb13 352 if ( !(AliCDBManager::Instance()->IsDefaultStorageSet()) ) {
353 AliError("No run data storage defined!");
354 return 0x0;
2c8628dd 355 }
9e1ceb13 356 return AliCDBManager::Instance()->GetDefaultStorage()->Get(name, GetHeader()->GetRun());
357
2c8628dd 358}
359
360/**************************************************************************/
88cb7938 361AliRunLoader* AliRunLoader::Open
362 (const char* filename, const char* eventfoldername, Option_t* option)
363{
364//Opens a desired file 'filename'
365//gets the the run-Loader and mounts it desired folder
366//returns the pointer to run Loader which can be further used for accessing data
367//in case of error returns NULL
368
024a7e64 369 static const TString kwebaddress("http://alisoft.cern.ch/people/skowron/codedoc/split/index.html");
21bf7095 370 AliDebugClass(1,Form("\n\n\nNew I/O strcture: See more info:\n %s\n\n\n",kwebaddress.Data()));
88dc6d4a 371
88cb7938 372 AliRunLoader* result = 0x0;
88dc6d4a 373
88cb7938 374 /* ************************************************ */
375 /* Chceck if folder with given name already exists */
376 /* ************************************************ */
88dc6d4a 377
88cb7938 378 TObject* obj = AliConfig::Instance()->GetTopFolder()->FindObject(eventfoldername);
379 if(obj)
380 {
381 TFolder* fold = dynamic_cast<TFolder*>(obj);
382 if (fold == 0x0)
383 {
21bf7095 384 AliErrorClass("Such a obejct already exists in top alice folder and it is not a folder.");
88cb7938 385 return 0x0;
386 }
387
388 //check if we can get RL from that folder
389 result = AliRunLoader::GetRunLoader(eventfoldername);
390 if (result == 0x0)
391 {
21bf7095 392 AliErrorClass(Form("Folder %s already exists, and can not find session there. Can not mount.",eventfoldername));
88cb7938 393 return 0x0;
394 }
395
396 if (result->GetFileName().CompareTo(filename) != 0)
397 {
21bf7095 398 AliErrorClass("Other file is mounted in demanded folder. Can not mount.");
88cb7938 399 return 0x0;
400 }
401
402 //check if now is demanded (re)creation
403 if ( AliLoader::TestFileOption(option) == kFALSE)
404 {
21bf7095 405 AliErrorClass(Form("Session already exists in folder %s and this session option is %s. Unable to proceed.",
406 eventfoldername,option));
88cb7938 407 return 0x0;
408 }
409
410 //check if demanded option is update and existing one
411 TString tmpstr(option);
412 if ( (tmpstr.CompareTo("update",TString::kIgnoreCase) == 0) &&
413 (result->fGAFile->IsWritable() == kFALSE) )
414 {
21bf7095 415 AliErrorClass(Form("Session already exists in folder %s and is not writable while this session option is %s. Unable to proceed.",
416 eventfoldername,option));
88cb7938 417 return 0x0;
418 }
419
21bf7095 420 AliWarningClass("Session is already opened and mounted in demanded folder");
bacbe0fd 421 if (!fgRunLoader) fgRunLoader = result; //PH get access from any place
88cb7938 422 return result;
423 } //end of checking in case of existance of object named identically that folder session is being opened
424
425
426 TFile * gAliceFile = TFile::Open(filename,option);//open a file
427 if (!gAliceFile)
428 {//null pointer returned
21bf7095 429 AliFatalClass(Form("Can not open file %s.",filename));
88cb7938 430 return 0x0;
431 }
432
433 if (gAliceFile->IsOpen() == kFALSE)
434 {//pointer to valid object returned but file is not opened
21bf7095 435 AliErrorClass(Form("Can not open file %s.",filename));
88cb7938 436 return 0x0;
437 }
438
439 //if file is "read" or "update" than we try to find AliRunLoader there - if not found cry and exit
440 //else create new AliRunLoader
441 if ( AliLoader::TestFileOption(option) )
442 {
21bf7095 443 AliDebugClass(1, "Reading RL from file");
88cb7938 444
445 result = dynamic_cast<AliRunLoader*>(gAliceFile->Get(fgkRunLoaderName));//get the run Loader from the file
446 if (result == 0x0)
447 {//didn't get
21bf7095 448 AliErrorClass(Form("Can not find run-Loader in file %s.",filename));
88cb7938 449 delete gAliceFile;//close the file
450 return 0x0;
451 }
452 Int_t tmp = result->SetEventFolderName(eventfoldername);//mount a event folder
453 if (tmp)//if SetEvent returned error
454 {
21bf7095 455 AliErrorClass(Form("Can not mount event in folder %s.",eventfoldername));
88cb7938 456 delete result; //delete run-Loader
457 delete gAliceFile;//close the file
458 return 0x0;
459 }
460 }
461 else
462 {
21bf7095 463 AliDebugClass(1, Form("Creating new AliRunLoader. Folder name is %s",eventfoldername));
131d919a 464 try
465 {
466 result = new AliRunLoader(eventfoldername);
467 }
468 catch (TString& errmsg)
469 {
21bf7095 470 AliErrorClass(Form("AliRunLoader constrcutor has thrown exception: %s\n",errmsg.Data()));
131d919a 471 delete result;
472 delete gAliceFile;//close the file
473 return 0x0;
474 }
88cb7938 475 }
476
477//procedure for extracting dir name from the file name
478 TString fname(filename);
fca6cd9f 479 Int_t nsl = fname.Last('#');//look for hash in file name
88cb7938 480 TString dirname;
fca6cd9f 481 if (nsl < 0) {//hash not found
482 nsl = fname.Last('/');// look for slash
483 if (nsl < 0)
484 nsl = fname.Last(':');// look for colon e.g. rfio:galice.root
485 }
486
487 if (nsl < 0) dirname = "./"; // no directory path, use "."
488 else dirname = fname.Remove(nsl+1);// directory path
88cb7938 489
21bf7095 490 AliDebugClass(1, Form("Dir name is : %s",dirname.Data()));
88cb7938 491
492 result->SetDirName(dirname);
493 result->SetGAliceFile(gAliceFile);//set the pointer to gAliceFile
85a5290f 494 if (!fgRunLoader) fgRunLoader = result; //PH get access from any place
88cb7938 495 return result;
496}
497/**************************************************************************/
498Int_t AliRunLoader::GetNumberOfEvents()
499{
500 //returns number of events in Run
501 Int_t retval;
502 if( TreeE() == 0x0 )
503 {
504 retval = LoadHeader();
505 if (retval)
506 {
21bf7095 507 AliError("Error occured while loading header");
88cb7938 508 return -1;
509 }
510 }
511 return (Int_t)TreeE()->GetEntries();
512}
88cb7938 513/**************************************************************************/
f2a509af 514
88cb7938 515void AliRunLoader::MakeHeader()
516{
517 //Makes header and connects it to header tree (if it exists)
21bf7095 518 AliDebug(1, "");
88cb7938 519 if(fHeader == 0x0)
520 {
21bf7095 521 AliDebug(1, "Creating new Header Object");
88cb7938 522 fHeader= new AliHeader();
523 }
524 TTree* tree = TreeE();
525 if (tree)
526 {
21bf7095 527 AliDebug(1, "Got Tree from folder.");
88cb7938 528 TBranch* branch = tree->GetBranch(fgkHeaderBranchName);
529 if (branch == 0x0)
530 {
21bf7095 531 AliDebug(1, "Creating new branch");
88cb7938 532 branch = tree->Branch(fgkHeaderBranchName, "AliHeader", &fHeader, 4000, 0);
533 branch->SetAutoDelete(kFALSE);
534 }
535 else
536 {
21bf7095 537 AliDebug(1, "Got Branch from Tree");
88cb7938 538 branch->SetAddress(&fHeader);
539 tree->GetEvent(fCurrentEvent);
540 fStack = fHeader->Stack(); //should be safe - if we created Stack, header returns pointer to the same object
541 if (fStack)
542 {
9ca424c1 543 if (TreeK()) {
544 fStack->ConnectTree(TreeK());
545 fStack->GetEvent();
546 }
88cb7938 547 }
548 else
549 {
21bf7095 550 AliDebug(1, "Header does not have a stack.");
88cb7938 551 }
552 }
553 }
21bf7095 554 AliDebug(1, "Exiting MakeHeader method");
88cb7938 555}
556/**************************************************************************/
557
558void AliRunLoader::MakeStack()
559{
560//Creates the stack object - do not connect the tree
561 if(fStack == 0x0)
562 {
563 fStack = new AliStack(10000);
88cb7938 564 }
565}
bacbe0fd 566/**************************************************************************/
567
568void AliRunLoader::MakeTrigger()
569{
570 // Makes trigger object and connects it to trigger tree (if it exists)
571 AliDebug( 1, "" );
572 if( fCTrigger == 0x0 ) {
573 AliDebug( 1, "Creating new Trigger Object" );
574 fCTrigger = new AliCentralTrigger();
575 }
576 TTree* tree = TreeCT();
577 if( tree ) {
578 fCTrigger->MakeBranch( fgkTriggerBranchName, tree );
579 tree->GetEvent( fCurrentEvent );
580 }
88cb7938 581
bacbe0fd 582 AliDebug( 1, "Exiting MakeTrigger method" );
583}
88cb7938 584/**************************************************************************/
585
586void AliRunLoader::MakeTree(Option_t *option)
587{
588//Creates trees
bacbe0fd 589 const char *oK = strstr(option,"K"); //Kine
590 const char *oE = strstr(option,"E"); //Header
0f46f5fa 591 const char *oGG = strstr(option,"GG"); //Central TriGGer
9ca424c1 592
593 if(oK)
594 {
595 if (fKineDataLoader->GetBaseLoader(0)->IsLoaded() == kFALSE)
88cb7938 596 {
9ca424c1 597 AliError("Load Kinematics first");
88cb7938 598 }
9ca424c1 599 else
88cb7938 600 {
9ca424c1 601 if (!TreeK()) {
602 fKineDataLoader->MakeTree();
603 MakeStack();
604 }
605 fStack->ConnectTree(TreeK());
606 WriteKinematics("OVERWRITE");
607 }
608 } // TreeK
88cb7938 609
610 if(oE && !TreeE())
611 {
612 fGAFile->cd();
613 TTree* tree = new TTree(fgkHeaderContainerName,"Tree with Headers");
614 GetEventFolder()->Add(tree);
615 MakeHeader();
616 WriteHeader("OVERWRITE");
617 }
618
0f46f5fa 619 if(oGG && !TreeCT())
bacbe0fd 620 {
621 // create the CTP Trigger output file and tree
622 TFile* file = gROOT->GetFile( fgkDefaultTriggerFileName );
623 if( !file ) {
624 file = TFile::Open( gSystem->ConcatFileName( fUnixDirName.Data(), fgkDefaultTriggerFileName.Data() ), "RECREATE" ) ;
625 }
626
627 file->cd();
628 TTree* tree = new TTree( fgkTriggerContainerName, "Tree with Central Trigger Mask" );
629 GetEventFolder()->Add(tree);
630 MakeTrigger();
631 // WriteHeader("OVERWRITE");
632 }
633
88cb7938 634 TIter next(fLoaders);
635 AliLoader *loader;
636 while((loader = (AliLoader*)next()))
637 {
638 loader->MakeTree(option);
639 }
640
641}
642/**************************************************************************/
643
644Int_t AliRunLoader::LoadgAlice()
645{
646//Loads gAlice from file
647 if (GetAliRun())
648 {
21bf7095 649 AliWarning("AliRun is already in folder. Unload first.");
88cb7938 650 return 0;
651 }
652 AliRun* alirun = dynamic_cast<AliRun*>(fGAFile->Get(fgkGAliceName));
653 if (alirun == 0x0)
654 {
21bf7095 655 AliError(Form("Can not find gAlice in file %s",fGAFile->GetName()));
88cb7938 656 return 2;
657 }
658 alirun->SetRunLoader(this);
659 if (gAlice)
660 {
21bf7095 661 AliWarning(Form("gAlice already exists. Putting retrived object in folder named %s",
662 GetEventFolder()->GetName()));
88cb7938 663 }
664 else
665 {
666 gAlice = alirun;
667 }
668 SetDetectorAddresses();//calls SetTreeAddress for all detectors
669 return 0;
670}
671/**************************************************************************/
672
673Int_t AliRunLoader::LoadHeader()
674{
f2a509af 675//loads treeE and reads header object for current event
88cb7938 676 if (TreeE())
677 {
21bf7095 678 AliWarning("Header is already loaded. Use ReloadHeader to force reload. Nothing done");
88cb7938 679 return 0;
680 }
681
682 if (GetEventFolder() == 0x0)
683 {
21bf7095 684 AliError("Event folder not specified yet");
88cb7938 685 return 1;
686 }
687
688 if (fGAFile == 0x0)
689 {
21bf7095 690 AliError("Session not opened. Use AliRunLoader::Open");
88cb7938 691 return 2;
692 }
693
694 if (fGAFile->IsOpen() == kFALSE)
695 {
21bf7095 696 AliError("Session not opened. Use AliRunLoader::Open");
88cb7938 697 return 2;
698 }
699
700 TTree* tree = dynamic_cast<TTree*>(fGAFile->Get(fgkHeaderContainerName));
701 if (tree == 0x0)
702 {
21bf7095 703 AliError(Form("Can not find header tree named %s in file %s",
704 fgkHeaderContainerName.Data(),fGAFile->GetName()));
88cb7938 705 return 2;
706 }
707
708 if (tree == TreeE()) return 0;
709
710 CleanHeader();
711 GetEventFolder()->Add(tree);
712 MakeHeader();//creates header object and connects to tree
713 return 0;
714
715}
716/**************************************************************************/
717
bacbe0fd 718Int_t AliRunLoader::LoadTrigger(Option_t* option)
719{
720 //Load treeCT
721
722 if( TreeCT() ) {
723 AliWarning("Trigger is already loaded. Nothing done");
724 return 0;
725 }
726
727 if( GetEventFolder() == 0x0 ) {
728 AliError("Event folder not specified yet");
729 return 1;
730 }
731 // get the CTP Trigger output file and tree
732 TString trgfile = gSystem->ConcatFileName( fUnixDirName.Data(),
733 fgkDefaultTriggerFileName.Data() );
734 TFile* file = gROOT->GetFile( trgfile );
735 if( !file ) {
736 file = TFile::Open( trgfile, option ) ;
737 if (!file || file->IsOpen() == kFALSE ) {
738 AliError( Form( "Can not open trigger file %s", trgfile.Data() ) );
739 return 2;
740 }
741 }
742 file->cd();
743
744 TTree* tree = dynamic_cast<TTree*>(file->Get( fgkTriggerContainerName ));
745 if( !tree ) {
746 AliError( Form( "Can not find trigger tree named %s in file %s",
747 fgkTriggerContainerName.Data(), file->GetName() ) );
748 return 2;
749 }
750
751 CleanTrigger();
752
753 fCTrigger = dynamic_cast<AliCentralTrigger*>(file->Get( "AliCentralTrigger" ));
754 GetEventFolder()->Add( tree );
755 MakeTrigger();
756
757 return 0;
758}
759
760/**************************************************************************/
761
88cb7938 762Int_t AliRunLoader::LoadKinematics(Option_t* option)
763{
88dc6d4a 764//Loads the kinematics
88cb7938 765 Int_t retval = fKineDataLoader->GetBaseLoader(0)->Load(option);
766 if (retval)
767 {
21bf7095 768 AliError("Error occured while loading kinamatics tree.");
88cb7938 769 return retval;
770 }
443b7449 771 if (fStack)
772 {
9ca424c1 773 fStack->ConnectTree(TreeK());
774 retval = fStack->GetEvent();
443b7449 775 if ( retval == kFALSE)
776 {
21bf7095 777 AliError("Error occured while loading kinamatics tree.");
443b7449 778 return retval;
779 }
780
781 }
88cb7938 782 return 0;
783}
784/**************************************************************************/
785
786Int_t AliRunLoader::OpenDataFile(const TString& filename,TFile*& file,TDirectory*& dir,Option_t* opt,Int_t cl)
787{
788//Opens File with kinematics
789 if (file)
790 {
791 if (file->IsOpen() == kFALSE)
792 {//pointer is not null but file is not opened
21bf7095 793 AliWarning("Pointer to file is not null, but file is not opened");//risky any way
88cb7938 794 delete file;
795 file = 0x0; //proceed with opening procedure
796 }
797 else
798 {
21bf7095 799 AliWarning(Form("File %s already opened",filename.Data()));
88cb7938 800 return 0;
801 }
802 }
803//try to find if that file is opened somewere else
804 file = (TFile *)( gROOT->GetListOfFiles()->FindObject(filename) );
805 if (file)
806 {
807 if(file->IsOpen() == kTRUE)
808 {
21bf7095 809 AliWarning(Form("File %s already opened by sombody else.",file->GetName()));
88cb7938 810 return 0;
811 }
812 }
813
814 file = TFile::Open(filename,opt);
815 if (file == 0x0)
816 {//file is null
21bf7095 817 AliError(Form("Can not open file %s",filename.Data()));
88cb7938 818 return 1;
819 }
820 if (file->IsOpen() == kFALSE)
821 {//file is not opened
21bf7095 822 AliError(Form("Can not open file %s",filename.Data()));
88cb7938 823 return 1;
824 }
825
826 file->SetCompressionLevel(cl);
827
828 dir = AliLoader::ChangeDir(file,fCurrentEvent);
829 if (dir == 0x0)
830 {
21bf7095 831 AliError(Form("Can not change to root directory in file %s",filename.Data()));
88cb7938 832 return 3;
833 }
834 return 0;
835}
836/**************************************************************************/
837
838TTree* AliRunLoader::TreeE() const
839{
840 //returns the tree from folder; shortcut method
21bf7095 841 if (AliDebugLevel() > 10) fEventFolder->ls();
88cb7938 842 TObject *obj = fEventFolder->FindObject(fgkHeaderContainerName);
843 return (obj)?dynamic_cast<TTree*>(obj):0x0;
844}
845/**************************************************************************/
846
bacbe0fd 847TTree* AliRunLoader::TreeCT() const
848{
849 //returns the tree from folder; shortcut method
850 if (AliDebugLevel() > 10) fEventFolder->ls();
851 TObject *obj = fEventFolder->FindObject(fgkTriggerContainerName);
852 return (obj)?dynamic_cast<TTree*>(obj):0x0;
853}
854/**************************************************************************/
855
88cb7938 856AliHeader* AliRunLoader::GetHeader() const
857{
f2a509af 858//returns pointer header object
88cb7938 859 return fHeader;
860}
861/**************************************************************************/
bacbe0fd 862
863AliCentralTrigger* AliRunLoader::GetTrigger() const
864{
865//returns pointer trigger object
866 return fCTrigger;
867}
868
869/**************************************************************************/
88cb7938 870
871TTree* AliRunLoader::TreeK() const
872{
873 //returns the tree from folder; shortcut method
874 TObject *obj = GetEventFolder()->FindObject(fgkKineContainerName);
875 return (obj)?dynamic_cast<TTree*>(obj):0x0;
876}
877/**************************************************************************/
878
879TTree* AliRunLoader::TreeTR() const
880{
881 //returns the tree from folder; shortcut method
882 TObject* obj = GetEventFolder()->FindObject(fgkTrackRefsContainerName);
883 return (obj)?dynamic_cast<TTree*>(obj):0x0;
884}
885/**************************************************************************/
886
887AliRun* AliRunLoader::GetAliRun() const
888{
889//returns AliRun which sits in the folder
890 if (fEventFolder == 0x0) return 0x0;
891 TObject *obj = fEventFolder->FindObject(fgkGAliceName);
892 return (obj)?dynamic_cast<AliRun*>(obj):0x0;
893}
894/**************************************************************************/
895
88cb7938 896Int_t AliRunLoader::WriteHeader(Option_t* opt)
897{
f2a509af 898//writes treeE
21bf7095 899 AliDebug(1, "WRITING HEADER");
88cb7938 900
901 TTree* tree = TreeE();
902 if ( tree == 0x0)
903 {
21bf7095 904 AliWarning("Can not find Header Tree in Folder");
88cb7938 905 return 0;
906 }
907 if (fGAFile->IsWritable() == kFALSE)
908 {
21bf7095 909 AliError(Form("File %s is not writable",fGAFile->GetName()));
88cb7938 910 return 1;
911 }
912
913 TObject* obj = fGAFile->Get(fgkHeaderContainerName);
914 if (obj)
915 { //if they exist, see if option OVERWRITE is used
916 TString tmp(opt);
917 if(tmp.Contains("OVERWRITE",TString::kIgnoreCase) == 0)
918 {//if it is not used - give an error message and return an error code
21bf7095 919 AliError("Tree already exisists. Use option \"OVERWRITE\" to overwrite previous data");
88cb7938 920 return 3;
921 }
922 }
923 fGAFile->cd();
924 tree->SetDirectory(fGAFile);
925 tree->Write(0,TObject::kOverwrite);
926
21bf7095 927 AliDebug(1, "WRITTEN\n\n");
88cb7938 928
929 return 0;
930}
bacbe0fd 931
932/**************************************************************************/
933
934Int_t AliRunLoader::WriteTrigger(Option_t* opt)
935{
936 //writes TreeCT
937 AliDebug( 1, "WRITING TRIGGER" );
938
939 TTree* tree = TreeCT();
940 if ( tree == 0x0) {
941 AliWarning("Can not find Trigger Tree in Folder");
942 return 0;
943 }
944
945 TFile* file = gROOT->GetFile( gSystem->ConcatFileName( fUnixDirName.Data(), fgkDefaultTriggerFileName.Data() ) ) ;
946 if( !file || !file->IsOpen() ) {
947 AliError( "can't write Trigger, file is not open" );
948 return kFALSE;
949 }
950
951 TObject* obj = file->Get( fgkTriggerContainerName );
952 if( obj ) { //if they exist, see if option OVERWRITE is used
953 TString tmp(opt);
954 if( tmp.Contains( "OVERWRITE", TString::kIgnoreCase ) == 0) {
955 //if it is not used - give an error message and return an error code
956 AliError( "Tree already exisists. Use option \"OVERWRITE\" to overwrite previous data" );
957 return 3;
958 }
959 }
960 file->cd();
961 fCTrigger->Write( 0, TObject::kOverwrite );
962 tree->Write( 0, TObject::kOverwrite );
963 file->Flush();
964
965 AliDebug(1, "WRITTEN\n\n");
966
967 return 0;
968}
88cb7938 969/**************************************************************************/
970
d1898505 971Int_t AliRunLoader::WriteAliRun(Option_t* /*opt*/)
88cb7938 972{
f2a509af 973//writes AliRun object to the file
88cb7938 974 fGAFile->cd();
f2a509af 975 if (GetAliRun()) GetAliRun()->Write();
88cb7938 976 return 0;
977}
978/**************************************************************************/
979
980Int_t AliRunLoader::WriteKinematics(Option_t* opt)
981{
f2a509af 982//writes Kinematics
88cb7938 983 return fKineDataLoader->GetBaseLoader(0)->WriteData(opt);
984}
985/**************************************************************************/
986Int_t AliRunLoader::WriteTrackRefs(Option_t* opt)
987{
f2a509af 988//writes Track References tree
88cb7938 989 return fTrackRefsDataLoader->GetBaseLoader(0)->WriteData(opt);
990}
991/**************************************************************************/
992
993Int_t AliRunLoader::WriteHits(Option_t* opt)
994{
995//Calls WriteHits for all loaders
996 Int_t res;
997 Int_t result = 0;
998 TIter next(fLoaders);
999 AliLoader *loader;
1000 while((loader = (AliLoader*)next()))
1001 {
1002 res = loader->WriteHits(opt);
1003 if (res)
1004 {
21bf7095 1005 AliError(Form("Failed to write hits for %s (%d)",loader->GetDetectorName().Data(),res));
88cb7938 1006 result = 1;
1007 }
1008 }
1009 return result;
1010}
1011/**************************************************************************/
1012
1013Int_t AliRunLoader::WriteSDigits(Option_t* opt)
1014{
f2a509af 1015//Calls WriteSDigits for all loaders
88cb7938 1016 Int_t res;
1017 Int_t result = 0;
1018 TIter next(fLoaders);
1019 AliLoader *loader;
1020 while((loader = (AliLoader*)next()))
1021 {
1022 res = loader->WriteSDigits(opt);
1023 if (res)
1024 {
21bf7095 1025 AliError(Form("Failed to write summable digits for %s.",loader->GetDetectorName().Data()));
88cb7938 1026 result = 1;
1027 }
1028 }
1029 return result;
1030}
1031/**************************************************************************/
1032
1033Int_t AliRunLoader::WriteDigits(Option_t* opt)
1034{
f2a509af 1035//Calls WriteDigits for all loaders
88cb7938 1036 Int_t res;
1037 Int_t result = 0;
1038 TIter next(fLoaders);
1039 AliLoader *loader;
1040 while((loader = (AliLoader*)next()))
1041 {
1042 res = loader->WriteDigits(opt);
1043 if (res)
1044 {
21bf7095 1045 AliError(Form("Failed to write digits for %s.",loader->GetDetectorName().Data()));
88cb7938 1046 result = 1;
1047 }
1048 }
1049 return result;
1050}
1051/**************************************************************************/
1052
1053Int_t AliRunLoader::WriteRecPoints(Option_t* opt)
1054{
f2a509af 1055//Calls WriteRecPoints for all loaders
88cb7938 1056 Int_t res;
1057 Int_t result = 0;
1058 TIter next(fLoaders);
1059 AliLoader *loader;
1060 while((loader = (AliLoader*)next()))
1061 {
1062 res = loader->WriteRecPoints(opt);
1063 if (res)
1064 {
21bf7095 1065 AliError(Form("Failed to write Reconstructed Points for %s.",
1066 loader->GetDetectorName().Data()));
88cb7938 1067 result = 1;
1068 }
1069 }
1070 return result;
1071}
1072/**************************************************************************/
1073
1074Int_t AliRunLoader::WriteTracks(Option_t* opt)
1075{
f2a509af 1076//Calls WriteTracks for all loaders
88cb7938 1077 Int_t res;
1078 Int_t result = 0;
1079 TIter next(fLoaders);
1080 AliLoader *loader;
1081 while((loader = (AliLoader*)next()))
1082 {
1083 res = loader->WriteTracks(opt);
1084 if (res)
1085 {
21bf7095 1086 AliError(Form("Failed to write Tracks for %s.",
1087 loader->GetDetectorName().Data()));
88cb7938 1088 result = 1;
1089 }
1090 }
1091 return result;
1092}
1093/**************************************************************************/
1094
d1898505 1095Int_t AliRunLoader::WriteRunLoader(Option_t* /*opt*/)
88cb7938 1096{
f2a509af 1097//Writes itself to the file
88cb7938 1098 CdGAFile();
1099 this->Write(0,TObject::kOverwrite);
1100 return 0;
1101}
1102/**************************************************************************/
1103
1104Int_t AliRunLoader::SetEventFolderName(const TString& name)
f2a509af 1105{
1106//sets top folder name for this run; of alread
88cb7938 1107 if (name.IsNull())
1108 {
21bf7095 1109 AliError("Name is empty");
88cb7938 1110 return 1;
1111 }
1112
1113 //check if such a folder already exists - try to find it in alice top folder
1114 TObject* obj = AliConfig::Instance()->GetTopFolder()->FindObject(name);
1115 if(obj)
1116 {
1117 TFolder* fold = dynamic_cast<TFolder*>(obj);
1118 if (fold == 0x0)
1119 {
21bf7095 1120 AliError("Such a obejct already exists in top alice folder and it is not a folder.");
88cb7938 1121 return 2;
1122 }
1123 //folder which was found is our folder
1124 if (fEventFolder == fold)
1125 {
1126 return 0;
1127 }
1128 else
1129 {
21bf7095 1130 AliError("Such a folder already exists in top alice folder. Can not mount.");
88cb7938 1131 return 2;
1132 }
1133 }
1134
1135 //event is alredy connected, just change name of the folder
1136 if (fEventFolder)
1137 {
1138 fEventFolder->SetName(name);
1139 return 0;
1140 }
1141
1142 if (fKineDataLoader == 0x0)
1143 fKineDataLoader = new AliDataLoader(fgkDefaultKineFileName,fgkKineContainerName,"Kinematics");
1144
1145 if ( fTrackRefsDataLoader == 0x0)
1146 fTrackRefsDataLoader = new AliDataLoader(fgkDefaultTrackRefsFileName,fgkTrackRefsContainerName,"Track References");
1147
1148 //build the event folder structure
21bf7095 1149 AliDebug(1, Form("Creating new event folder named %s",name.Data()));
88cb7938 1150 fEventFolder = AliConfig::Instance()->BuildEventFolder(name,"Event Folder");
1151 fEventFolder->Add(this);//put myself to the folder to accessible for all
1152
88cb7938 1153 TIter next(fLoaders);
1154 AliLoader *loader;
1155 while((loader = (AliLoader*)next()))
1156 {
1157 loader->Register(fEventFolder);//build folder structure for this detector
1158 }
1159
1160 fKineDataLoader->SetEventFolder(GetEventFolder());
1161 fTrackRefsDataLoader->SetEventFolder(GetEventFolder());
1162 fKineDataLoader->SetFolder(GetEventFolder());
1163 fTrackRefsDataLoader->SetFolder(GetEventFolder());
1164
1165 fEventFolder->SetOwner();
1166 return 0;
1167}
1168/**************************************************************************/
1169
1170void AliRunLoader::AddLoader(AliLoader* loader)
1171 {
1172 //Adds the Loader for given detector
1173 if (loader == 0x0) //if null shout and exit
1174 {
21bf7095 1175 AliError("Parameter is NULL");
88cb7938 1176 return;
1177 }
1178 loader->SetDirName(fUnixDirName);
1179 if (fEventFolder) loader->SetEventFolder(fEventFolder); //if event folder is already defined,
1180 //pass information to the Loader
1181 fLoaders->Add(loader);//add the Loader to the array
1182 }
1183/**************************************************************************/
1184
1185void AliRunLoader::AddLoader(AliDetector* det)
1186 {
1187//Asks module (detector) ro make a Loader and stores in the array
1188 if (det == 0x0) return;
1189 AliLoader* get = det->GetLoader();//try to get loader
1190 if (get == 0x0) get = det->MakeLoader(fEventFolder->GetName());//if did not obtain, ask to make it
1191
1192 if (get)
1193 {
21bf7095 1194 AliDebug(1, Form("Detector: %s Loader : %s",det->GetName(),get->GetName()));
88cb7938 1195 AddLoader(get);
1196 }
1197 }
1198
1199/**************************************************************************/
1200
1201AliLoader* AliRunLoader::GetLoader(const char* detname) const
1202{
f2a509af 1203//returns loader for given detector
1204//note that naming convention is TPCLoader not just TPC
88cb7938 1205 return (AliLoader*)fLoaders->FindObject(detname);
1206}
1207
1208/**************************************************************************/
1209
1210AliLoader* AliRunLoader::GetLoader(AliDetector* det) const
1211{
f2a509af 1212//get loader for detector det
88cb7938 1213 if(det == 0x0) return 0x0;
1214 TString getname(det->GetName());
1215 getname+="Loader";
21bf7095 1216 AliDebug(1, Form(" Loader name is %s",getname.Data()));
88cb7938 1217 return GetLoader(getname);
1218}
1219
1220/**************************************************************************/
1221
1222void AliRunLoader::CleanFolders()
1223{
1224// fEventFolder->Add(this);//put myself to the folder to accessible for all
1225
1226 CleanDetectors();
1227 CleanHeader();
1228 CleanKinematics();
bacbe0fd 1229 CleanTrigger();
88cb7938 1230}
1231/**************************************************************************/
1232
1233void AliRunLoader::CleanDetectors()
1234{
1235//Calls CleanFolders for all detectors
1236 TIter next(fLoaders);
d0d4a6b3 1237 AliLoader *loader;
1238 while((loader = (AliLoader*)next()))
88cb7938 1239 {
d0d4a6b3 1240 loader->CleanFolders();
88cb7938 1241 }
1242}
1243/**************************************************************************/
1244
1245void AliRunLoader::RemoveEventFolder()
1246{
1247//remove all the tree of event
1248//all the stuff changing EbE stays untached (PDGDB, tasks, etc.)
1249
1250 if (fEventFolder == 0x0) return;
1251 fEventFolder->SetOwner(kFALSE);//don't we want to deleted while removing the folder that we are sitting in
1252 fEventFolder->Remove(this);//remove us drom folder
1253
1254 AliConfig::Instance()->GetTopFolder()->SetOwner(); //brings ownership back for fEventFolder since it sits in top folder
1255 AliConfig::Instance()->GetTopFolder()->Remove(fEventFolder); //remove the event tree
1256 delete fEventFolder;
1257}
1258/**************************************************************************/
1259
1260void AliRunLoader::SetGAliceFile(TFile* gafile)
1261{
f2a509af 1262//sets pointer to galice.root file
88cb7938 1263 fGAFile = gafile;
1264}
1265
1266/**************************************************************************/
1267
1268Int_t AliRunLoader::LoadHits(Option_t* detectors,Option_t* opt)
1269{
1270//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1271
21bf7095 1272 AliDebug(1, "Loading Hits");
88cb7938 1273 TObjArray* loaders;
1274 TObjArray arr;
1275
1276 const char* oAll = strstr(detectors,"all");
1277 if (oAll)
1278 {
21bf7095 1279 AliDebug(1, "Option is All");
88cb7938 1280 loaders = fLoaders;
1281 }
1282 else
1283 {
1284 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1285 loaders = &arr;//get the pointer array
1286 }
1287
21bf7095 1288 AliDebug(1, Form("For detectors. Number of detectors chosen for loading %d",loaders->GetEntries()));
88cb7938 1289
1290 TIter next(loaders);
1291 AliLoader *loader;
1292 while((loader = (AliLoader*)next()))
1293 {
21bf7095 1294 AliDebug(1, Form(" Calling LoadHits(%s) for %s",opt,loader->GetName()));
88cb7938 1295 loader->LoadHits(opt);
1296 }
21bf7095 1297 AliDebug(1, "Done");
88cb7938 1298 return 0;
1299}
1300
1301/**************************************************************************/
1302
1303Int_t AliRunLoader::LoadSDigits(Option_t* detectors,Option_t* opt)
1304{
1305//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1306
1307 TObjArray* loaders;
1308 TObjArray arr;
1309
1310 const char* oAll = strstr(detectors,"all");
1311 if (oAll)
1312 {
1313 loaders = fLoaders;
1314 }
1315 else
1316 {
1317 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1bb20a37 1318 loaders = &arr;//get the pointer to array
88cb7938 1319 }
1320
1321 TIter next(loaders);
1322 AliLoader *loader;
1323 while((loader = (AliLoader*)next()))
1324 {
1325 loader->LoadSDigits(opt);
1326 }
1327 return 0;
1328}
1329
1330/**************************************************************************/
1331
1332Int_t AliRunLoader::LoadDigits(Option_t* detectors,Option_t* opt)
1333{
1334//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1335
d0d4a6b3 1336 TObjArray* loaders;
88cb7938 1337 TObjArray arr;
1338
1339 const char* oAll = strstr(detectors,"all");
1340 if (oAll)
1341 {
d0d4a6b3 1342 loaders = fLoaders;
88cb7938 1343 }
1344 else
1345 {
1346 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
d0d4a6b3 1347 loaders = &arr;//get the pointer array
88cb7938 1348 }
1349
d0d4a6b3 1350 TIter next(loaders);
1351 AliLoader *loader;
1352 while((loader = (AliLoader*)next()))
88cb7938 1353 {
d0d4a6b3 1354 loader->LoadDigits(opt);
88cb7938 1355 }
1356 return 0;
1357}
88cb7938 1358/**************************************************************************/
1359
1360Int_t AliRunLoader::LoadRecPoints(Option_t* detectors,Option_t* opt)
1361{
1362//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1363
d0d4a6b3 1364 TObjArray* loaders;
88cb7938 1365 TObjArray arr;
1366
1367 const char* oAll = strstr(detectors,"all");
1368 if (oAll)
1369 {
d0d4a6b3 1370 loaders = fLoaders;
88cb7938 1371 }
1372 else
1373 {
1374 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
d0d4a6b3 1375 loaders = &arr;//get the pointer array
88cb7938 1376 }
1377
d0d4a6b3 1378 TIter next(loaders);
1379 AliLoader *loader;
1380 while((loader = (AliLoader*)next()))
88cb7938 1381 {
d0d4a6b3 1382 loader->LoadRecPoints(opt);
88cb7938 1383 }
1384 return 0;
1385}
1bb20a37 1386/**************************************************************************/
88cb7938 1387
1bb20a37 1388Int_t AliRunLoader::LoadRecParticles(Option_t* detectors,Option_t* opt)
1389{
1390//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1391
d0d4a6b3 1392 TObjArray* loaders;
1bb20a37 1393 TObjArray arr;
1394
1395 const char* oAll = strstr(detectors,"all");
1396 if (oAll)
1397 {
d0d4a6b3 1398 loaders = fLoaders;
1bb20a37 1399 }
1400 else
1401 {
1402 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
d0d4a6b3 1403 loaders = &arr;//get the pointer array
1bb20a37 1404 }
1405
d0d4a6b3 1406 TIter next(loaders);
1407 AliLoader *loader;
1408 while((loader = (AliLoader*)next()))
1bb20a37 1409 {
d0d4a6b3 1410 loader->LoadRecParticles(opt);
1bb20a37 1411 }
1412 return 0;
1413}
88cb7938 1414/**************************************************************************/
1415
1416Int_t AliRunLoader::LoadTracks(Option_t* detectors,Option_t* opt)
1417{
1418//LoadHits in selected detectors i.e. detectors="ITS TPC TRD" or "all"
1419
d0d4a6b3 1420 TObjArray* loaders;
88cb7938 1421 TObjArray arr;
1422
1423 const char* oAll = strstr(detectors,"all");
1424 if (oAll)
1425 {
d0d4a6b3 1426 loaders = fLoaders;
88cb7938 1427 }
1428 else
1429 {
1430 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
d0d4a6b3 1431 loaders = &arr;//get the pointer array
88cb7938 1432 }
1433
d0d4a6b3 1434 TIter next(loaders);
1435 AliLoader *loader;
1436 while((loader = (AliLoader*)next()))
88cb7938 1437 {
d0d4a6b3 1438 loader->LoadTracks(opt);
88cb7938 1439 }
1440 return 0;
1441}
1bb20a37 1442/**************************************************************************/
1443
1444void AliRunLoader::UnloadHits(Option_t* detectors)
1445{
1446 //unloads hits for detectors specified in parameter
1447 TObjArray* loaders;
1448 TObjArray arr;
1449
1450 const char* oAll = strstr(detectors,"all");
1451 if (oAll)
1452 {
1453 loaders = fLoaders;
1454 }
1455 else
1456 {
1457 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1458 loaders = &arr;//get the pointer to array
1459 }
1460
1461 TIter next(loaders);
1462 AliLoader *loader;
1463 while((loader = (AliLoader*)next()))
1464 {
1465 loader->UnloadHits();
1466 }
1467}
1468/**************************************************************************/
1469
1470void AliRunLoader::UnloadSDigits(Option_t* detectors)
1471{
1472 //unloads SDigits for detectors specified in parameter
1473 TObjArray* loaders;
1474 TObjArray arr;
1475
1476 const char* oAll = strstr(detectors,"all");
1477 if (oAll)
1478 {
1479 loaders = fLoaders;
1480 }
1481 else
1482 {
1483 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1484 loaders = &arr;//get the pointer to array
1485 }
1486
1487 TIter next(loaders);
1488 AliLoader *loader;
1489 while((loader = (AliLoader*)next()))
1490 {
1491 loader->UnloadSDigits();
1492 }
1493}
1494/**************************************************************************/
1495
1496void AliRunLoader::UnloadDigits(Option_t* detectors)
1497{
1498 //unloads Digits for detectors specified in parameter
1499 TObjArray* loaders;
1500 TObjArray arr;
1501
1502 const char* oAll = strstr(detectors,"all");
1503 if (oAll)
1504 {
1505 loaders = fLoaders;
1506 }
1507 else
1508 {
1509 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1510 loaders = &arr;//get the pointer to array
1511 }
1512
1513 TIter next(loaders);
1514 AliLoader *loader;
1515 while((loader = (AliLoader*)next()))
1516 {
1517 loader->UnloadDigits();
1518 }
1519}
1520/**************************************************************************/
1521
1522void AliRunLoader::UnloadRecPoints(Option_t* detectors)
1523{
1524 //unloads RecPoints for detectors specified in parameter
1525 TObjArray* loaders;
1526 TObjArray arr;
1527
1528 const char* oAll = strstr(detectors,"all");
1529 if (oAll)
1530 {
1531 loaders = fLoaders;
1532 }
1533 else
1534 {
1535 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1536 loaders = &arr;//get the pointer to array
1537 }
1538
1539 TIter next(loaders);
1540 AliLoader *loader;
1541 while((loader = (AliLoader*)next()))
1542 {
1543 loader->UnloadRecPoints();
1544 }
1545}
1546/**************************************************************************/
1547
1548void AliRunLoader::UnloadAll(Option_t* detectors)
1549{
1550 //calls UnloadAll for detectors names specified in parameter
1551 // option "all" passed can be passed
1552 TObjArray* loaders;
1553 TObjArray arr;
1554
1555 const char* oAll = strstr(detectors,"all");
1556 if (oAll)
1557 {
1558 loaders = fLoaders;
1559 }
1560 else
1561 {
1562 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1563 loaders = &arr;//get the pointer to array
1564 }
88cb7938 1565
1bb20a37 1566 TIter next(loaders);
1567 AliLoader *loader;
1568 while((loader = (AliLoader*)next()))
1569 {
1570 loader->UnloadAll();
1571 }
1572}
1573/**************************************************************************/
1574
1575void AliRunLoader::UnloadTracks(Option_t* detectors)
1576{
1577 //unloads Tracks for detectors specified in parameter
1578 TObjArray* loaders;
1579 TObjArray arr;
1580
1581 const char* oAll = strstr(detectors,"all");
1582 if (oAll)
1583 {
1584 loaders = fLoaders;
1585 }
1586 else
1587 {
1588 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1589 loaders = &arr;//get the pointer to array
1590 }
1591
1592 TIter next(loaders);
1593 AliLoader *loader;
1594 while((loader = (AliLoader*)next()))
1595 {
1596 loader->UnloadTracks();
1597 }
1598}
1599/**************************************************************************/
1600
1601void AliRunLoader::UnloadRecParticles(Option_t* detectors)
1602{
1603 //unloads Particles for detectors specified in parameter
1604 TObjArray* loaders;
1605 TObjArray arr;
1606
1607 const char* oAll = strstr(detectors,"all");
1608 if (oAll)
1609 {
1610 loaders = fLoaders;
1611 }
1612 else
1613 {
1614 GetListOfDetectors(detectors,arr);//this method looks for all Loaders corresponding to names (many) specified in detectors option
1615 loaders = &arr;//get the pointer to array
1616 }
1617
1618 TIter next(loaders);
1619 AliLoader *loader;
1620 while((loader = (AliLoader*)next()))
1621 {
1622 loader->UnloadRecParticles();
1623 }
1624}
88cb7938 1625/**************************************************************************/
1626
1627AliRunLoader* AliRunLoader::GetRunLoader(const char* eventfoldername)
f2a509af 1628{
1629//returns RunLoader from folder named eventfoldername
88cb7938 1630 TFolder* evfold= dynamic_cast<TFolder*>(AliConfig::Instance()->GetTopFolder()->FindObject(eventfoldername));
1631 if (evfold == 0x0)
1632 {
1633 return 0x0;
1634 }
1635 AliRunLoader* runget = dynamic_cast<AliRunLoader*>(evfold->FindObject(AliRunLoader::fgkRunLoaderName));
1636 return runget;
1637
f2a509af 1638}
88cb7938 1639/**************************************************************************/
1640
1bb20a37 1641AliLoader* AliRunLoader::GetDetectorLoader(const char* detname, const char* eventfoldername)
8de97894 1642{
1643//get the loader of the detector with the given name from the global
1644//run loader object
1bb20a37 1645 AliRunLoader* runLoader = GetRunLoader(eventfoldername);
8de97894 1646 if (!runLoader) {
21bf7095 1647 AliErrorClass("No run loader found");
8de97894 1648 return NULL;
1649 }
1bb20a37 1650 return runLoader->GetDetectorLoader(detname);
1651}
1652/**************************************************************************/
1653
1654AliLoader* AliRunLoader::GetDetectorLoader(const char* detname)
1655{
1656//get the loader of the detector with the given name from the global
1657//run loader object
1658
8de97894 1659 char loadername[256];
1660 sprintf(loadername, "%sLoader", detname);
1bb20a37 1661 AliLoader* loader = GetLoader(loadername);
8de97894 1662 if (!loader) {
21bf7095 1663 AliError(Form("No loader for %s found", detname));
8de97894 1664 return NULL;
1665 }
1666 return loader;
1667}
1bb20a37 1668/**************************************************************************/
1669
1670TTree* AliRunLoader::GetTreeH(const char* detname, Bool_t maketree, const char* eventfoldername)
1671{
1672//get the tree with hits of the detector with the given name
1673//if maketree is true and the tree does not exist, the tree is created
1674 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1675 if (!loader) return NULL;
1676 if (!loader->TreeH() && maketree) loader->MakeTree("H");
1677 return loader->TreeH();
1678}
8de97894 1679
1680/**************************************************************************/
1681
1682TTree* AliRunLoader::GetTreeH(const char* detname, Bool_t maketree)
1683{
1684//get the tree with hits of the detector with the given name
1685//if maketree is true and the tree does not exist, the tree is created
1686 AliLoader* loader = GetDetectorLoader(detname);
1687 if (!loader) return NULL;
1688 if (!loader->TreeH() && maketree) loader->MakeTree("H");
1689 return loader->TreeH();
1690}
1bb20a37 1691/**************************************************************************/
8de97894 1692
1bb20a37 1693TTree* AliRunLoader::GetTreeS(const char* detname, Bool_t maketree,const char* eventfoldername)
1694{
1695//get the tree with summable digits of the detector with the given name
1696//if maketree is true and the tree does not exist, the tree is created
1697 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1698 if (!loader) return NULL;
1699 if (!loader->TreeS() && maketree) loader->MakeTree("S");
1700 return loader->TreeS();
1701}
8de97894 1702/**************************************************************************/
1703
1704TTree* AliRunLoader::GetTreeS(const char* detname, Bool_t maketree)
1705{
1706//get the tree with summable digits of the detector with the given name
1707//if maketree is true and the tree does not exist, the tree is created
1708 AliLoader* loader = GetDetectorLoader(detname);
1709 if (!loader) return NULL;
1710 if (!loader->TreeS() && maketree) loader->MakeTree("S");
1711 return loader->TreeS();
1712}
1bb20a37 1713/**************************************************************************/
8de97894 1714
1bb20a37 1715TTree* AliRunLoader::GetTreeD(const char* detname, Bool_t maketree,const char* eventfoldername)
1716{
1717//get the tree with digits of the detector with the given name
1718//if maketree is true and the tree does not exist, the tree is created
1719 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1720 if (!loader) return NULL;
1721 if (!loader->TreeD() && maketree) loader->MakeTree("D");
1722 return loader->TreeD();
1723}
8de97894 1724/**************************************************************************/
1725
1726TTree* AliRunLoader::GetTreeD(const char* detname, Bool_t maketree)
1727{
1728//get the tree with digits of the detector with the given name
1729//if maketree is true and the tree does not exist, the tree is created
1730 AliLoader* loader = GetDetectorLoader(detname);
1731 if (!loader) return NULL;
1732 if (!loader->TreeD() && maketree) loader->MakeTree("D");
1733 return loader->TreeD();
1734}
1bb20a37 1735/**************************************************************************/
1736TTree* AliRunLoader::GetTreeR(const char* detname, Bool_t maketree,const char* eventfoldername)
1737{
1738//get the tree with clusters of the detector with the given name
1739//if maketree is true and the tree does not exist, the tree is created
1740 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1741 if (!loader) return NULL;
1742 if (!loader->TreeR() && maketree) loader->MakeTree("R");
1743 return loader->TreeR();
1744}
8de97894 1745/**************************************************************************/
1746
1747TTree* AliRunLoader::GetTreeR(const char* detname, Bool_t maketree)
1748{
1749//get the tree with clusters of the detector with the given name
1750//if maketree is true and the tree does not exist, the tree is created
1751 AliLoader* loader = GetDetectorLoader(detname);
1752 if (!loader) return NULL;
1753 if (!loader->TreeR() && maketree) loader->MakeTree("R");
1754 return loader->TreeR();
1755}
1bb20a37 1756/**************************************************************************/
8de97894 1757
1bb20a37 1758TTree* AliRunLoader::GetTreeT(const char* detname, Bool_t maketree,const char* eventfoldername)
1759{
1760//get the tree with tracks of the detector with the given name
1761//if maketree is true and the tree does not exist, the tree is created
1762 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1763 if (!loader) return NULL;
1764 if (!loader->TreeT() && maketree) loader->MakeTree("T");
1765 return loader->TreeT();
1766}
8de97894 1767/**************************************************************************/
1768
1769TTree* AliRunLoader::GetTreeT(const char* detname, Bool_t maketree)
1770{
1771//get the tree with tracks of the detector with the given name
1772//if maketree is true and the tree does not exist, the tree is created
1773 AliLoader* loader = GetDetectorLoader(detname);
1774 if (!loader) return NULL;
1775 if (!loader->TreeT() && maketree) loader->MakeTree("T");
1776 return loader->TreeT();
1777}
1bb20a37 1778/**************************************************************************/
8de97894 1779
1bb20a37 1780TTree* AliRunLoader::GetTreeP(const char* detname, Bool_t maketree,const char* eventfoldername)
1781{
1782//get the tree with particles of the detector with the given name
1783//if maketree is true and the tree does not exist, the tree is created
1784 AliLoader* loader = GetDetectorLoader(detname,eventfoldername);
1785 if (!loader) return NULL;
1786 if (!loader->TreeP() && maketree) loader->MakeTree("P");
1787 return loader->TreeP();
1788}
8de97894 1789/**************************************************************************/
1790
1791TTree* AliRunLoader::GetTreeP(const char* detname, Bool_t maketree)
1792{
1793//get the tree with particles of the detector with the given name
1794//if maketree is true and the tree does not exist, the tree is created
1795 AliLoader* loader = GetDetectorLoader(detname);
1796 if (!loader) return NULL;
1797 if (!loader->TreeP() && maketree) loader->MakeTree("P");
1798 return loader->TreeP();
1799}
1800
1801/**************************************************************************/
1802
88cb7938 1803void AliRunLoader::CdGAFile()
1804{
1805//sets gDirectory to galice file
1806//work around
1807 if(fGAFile) fGAFile->cd();
1808}
1809
1810/**************************************************************************/
1811
1812void AliRunLoader::GetListOfDetectors(const char * namelist,TObjArray& pointerarray) const
1813 {
1814//this method looks for all Loaders corresponding
1815//to names (many) specified in namelist i.e. namelist ("ITS TPC TRD")
1816
1817 char buff[10];
1818 char dets [200];
1819 strcpy(dets,namelist);//compiler cries when char* = const Option_t*;
1820 dets[strlen(dets)+1] = '\n';//set endl at the end of string
1821 char* pdet = dets;
1822 Int_t tmp;
1823 for(;;)
1824 {
1825 tmp = sscanf(pdet,"%s",buff);//read the string from the input string pdet into buff
1826 if ( (buff[0] == 0) || (tmp == 0) ) break; //if not read
1827
f2a509af 1828 pdet = strstr(pdet,buff) + strlen(buff);//move the input pointer about number of bytes (letters) read
88cb7938 1829 //I am aware that is a little complicated. I don't know the number of spaces between detector names
1830 //so I read the string, than I find where it starts (strstr) and move the pointer about length of a string
1831 //If there is a better way, please write me (Piotr.Skowronski@cern.ch)
1832 //construct the Loader name
1833 TString getname(buff);
1834 getname+="Loader";
1835 AliLoader* loader = GetLoader(getname);//get the Loader
1836 if (loader)
1837 {
1838 pointerarray.Add(loader);
1839 }
1840 else
1841 {
21bf7095 1842 AliError(Form("Can not find Loader for %s",buff));
88cb7938 1843 }
1844
1845 buff[0] = 0;
1846 }
1847 }
1848/*****************************************************************************/
1849
1850void AliRunLoader::Clean(const TString& name)
1851{
1852//removes object with given name from event folder and deletes it
1853 if (GetEventFolder() == 0x0) return;
1854 TObject* obj = GetEventFolder()->FindObject(name);
1855 if(obj)
1856 {
21bf7095 1857 AliDebug(1, Form("name=%s, cleaning %s.",GetName(),name.Data()));
88cb7938 1858 GetEventFolder()->Remove(obj);
1859 delete obj;
bacbe0fd 1860 obj = 0x0;
88cb7938 1861 }
1862}
1863
1864/*****************************************************************************/
1865
1866TTask* AliRunLoader::GetRunDigitizer()
1867{
1868//returns Run Digitizer from folder
1869
1870 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1871 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetDigitizerTaskName());
1872 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1873}
1874/*****************************************************************************/
1875
1876TTask* AliRunLoader::GetRunSDigitizer()
1877{
1878//returns SDigitizer Task from folder
1879
1880 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1881 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetSDigitizerTaskName());
1882 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1883}
1884/*****************************************************************************/
1885
1886TTask* AliRunLoader::GetRunReconstructioner()
1887{
1888//returns Reconstructioner Task from folder
1889 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1890 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetReconstructionerTaskName());
1891 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1892}
1893/*****************************************************************************/
1894
1895TTask* AliRunLoader::GetRunTracker()
1896{
1897//returns Tracker Task from folder
1898 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1899 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetTrackerTaskName());
1900 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1901}
1902/*****************************************************************************/
1903
1904TTask* AliRunLoader::GetRunPIDTask()
1905{
1906//returns Tracker Task from folder
1907 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1908 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetPIDTaskName());
1909 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1910}
1911/*****************************************************************************/
1912
1913TTask* AliRunLoader::GetRunQATask()
1914{
1915//returns Quality Assurance Task from folder
1916 TFolder* topf = AliConfig::Instance()->GetTaskFolder();
1917 if (topf == 0x0)
1918 {
21bf7095 1919 AliErrorClass("Can not get task folder from AliConfig");
88cb7938 1920 return 0x0;
1921 }
1922 TObject* obj = topf->FindObjectAny(AliConfig::Instance()->GetQATaskName());
1923 return (obj)?dynamic_cast<TTask*>(obj):0x0;
1924}
1925
1926/*****************************************************************************/
1927
1928void AliRunLoader::SetCompressionLevel(Int_t cl)
1929{
1930//Sets Compression Level in all files
1931 if (fGAFile) fGAFile->SetCompressionLevel(cl);
1932 SetKineComprLevel(cl);
1933 SetTrackRefsComprLevel(cl);
1934 TIter next(fLoaders);
1935 AliLoader *loader;
1936 while((loader = (AliLoader*)next()))
1937 {
1938 loader->SetCompressionLevel(cl);
1939 }
1940}
1941/**************************************************************************/
1942
1943void AliRunLoader::SetKineComprLevel(Int_t cl)
1944{
1945//Sets comression level in Kine File
1946 fKineDataLoader->SetCompressionLevel(cl);
1947}
1948/**************************************************************************/
1949
1950void AliRunLoader::SetTrackRefsComprLevel(Int_t cl)
1951{
1952//Sets comression level in Track Refences File
1953 fTrackRefsDataLoader->SetCompressionLevel(cl);
1954}
1955/**************************************************************************/
1956
1957void AliRunLoader::UnloadHeader()
1958{
1959 //removes TreeE from folder and deletes it
1960 // as well as fHeader object
1961 CleanHeader();
1962 delete fHeader;
1963 fHeader = 0x0;
1964}
1965/**************************************************************************/
1966
bacbe0fd 1967void AliRunLoader::UnloadTrigger()
1968{
1969 //removes TreeCT from folder and deletes it
1970 // as well as fHeader object
1971 CleanTrigger();
1972 delete fCTrigger;
1973 fCTrigger = 0x0;
1974}
1975
1976/**************************************************************************/
1977
88cb7938 1978void AliRunLoader::UnloadKinematics()
1979{
f2a509af 1980//Unloads Kinematics
88cb7938 1981 fKineDataLoader->GetBaseLoader(0)->Unload();
1982}
1983/**************************************************************************/
1984
1985void AliRunLoader::UnloadTrackRefs()
1986{
f2a509af 1987//Unloads Track Refernces
88cb7938 1988 fTrackRefsDataLoader->GetBaseLoader(0)->Unload();
1989}
1990/**************************************************************************/
1991
1992void AliRunLoader::UnloadgAlice()
1993{
f2a509af 1994//Unloads gAlice
88cb7938 1995 if (gAlice == GetAliRun())
1996 {
21bf7095 1997 AliDebug(1, "Set gAlice = 0x0");
88cb7938 1998 gAlice = 0x0;//if gAlice is the same that in folder (to be deleted by the way of folder)
1999 }
2000 AliRun* alirun = GetAliRun();
2001 if (GetEventFolder()) GetEventFolder()->Remove(alirun);
2002 delete alirun;
2003}
2004/**************************************************************************/
2005
2006void AliRunLoader::MakeTrackRefsContainer()
2007{
2008// Makes a tree for Track References
2009 fTrackRefsDataLoader->MakeTree();
2010}
2011/**************************************************************************/
2012
2013Int_t AliRunLoader::LoadTrackRefs(Option_t* option)
2014{
2015//Load track references from file (opens file and posts tree to folder)
2016
f2a509af 2017 return fTrackRefsDataLoader->GetBaseLoader(0)->Load(option);
88cb7938 2018}
2019/**************************************************************************/
2020
2021void AliRunLoader::SetDirName(TString& dirname)
2022{
2023//sets directory name
2024 if (dirname.IsNull()) return;
2025 fUnixDirName = dirname;
2026 fKineDataLoader->SetDirName(dirname);
2027 fTrackRefsDataLoader->SetDirName(dirname);
2028
2029 TIter next(fLoaders);
2030 AliLoader *loader;
2031 while((loader = (AliLoader*)next()))
2032 {
2033 loader->SetDirName(dirname);
2034 }
2035
2036}
2037/*****************************************************************************/
2038
2039Int_t AliRunLoader::GetFileOffset() const
2040{
f2a509af 2041//returns the file number that is added to the file name for current event
88cb7938 2042 return Int_t(fCurrentEvent/fNEventsPerFile);
2043}
2044
2045/*****************************************************************************/
2046const TString AliRunLoader::SetFileOffset(const TString& fname)
2047{
f2a509af 2048//adds the the number to the file name at proper place for current event
88cb7938 2049 Long_t offset = (Long_t)GetFileOffset();
2050 if (offset < 1) return fname;
2051 TString soffset;
2052 soffset += offset;//automatic conversion to string
2053 TString dotroot(".root");
2054 const TString& offfsetdotroot = offset + dotroot;
2055 TString out = fname;
2056 out = out.ReplaceAll(dotroot,offfsetdotroot);
21bf7095 2057 AliDebug(1, Form(" in=%s out=%s",fname.Data(),out.Data()));
88cb7938 2058 return out;
2059}
2060/*****************************************************************************/
2061
2062void AliRunLoader::SetDigitsFileNameSuffix(const TString& suffix)
2063{
2064//adds the suffix before ".root",
2065//e.g. TPC.Digits.root -> TPC.DigitsMerged.root
2066//made on Jiri Chudoba demand
2067
2068 TIter next(fLoaders);
d0d4a6b3 2069 AliLoader *loader;
2070 while((loader = (AliLoader*)next()))
88cb7938 2071 {
d0d4a6b3 2072 loader->SetDigitsFileNameSuffix(suffix);
88cb7938 2073 }
2074}
2075/*****************************************************************************/
2076
2077TString AliRunLoader::GetFileName() const
2078{
2079//returns name of galice file
2080 TString result;
2081 if (fGAFile == 0x0) return result;
2082 result = fGAFile->GetName();
2083 return result;
2084}
2085/*****************************************************************************/
2086
2087void AliRunLoader::SetDetectorAddresses()
2088{
2089 //calls SetTreeAddress for all detectors
2090 if (GetAliRun()==0x0) return;
2091 TIter next(GetAliRun()->Modules());
2092 AliModule* mod;
2093 while((mod = (AliModule*)next()))
2094 {
2095 AliDetector* det = dynamic_cast<AliDetector*>(mod);
2096 if (det) det->SetTreeAddress();
2097 }
2098}
2099/*****************************************************************************/
f0f6f856 2100
2101void AliRunLoader::Synchronize()
2102{
2103 //synchrinizes all writtable files
2104 TIter next(fLoaders);
2105 AliLoader *loader;
2106 while((loader = (AliLoader*)next()))
2107 {
2108 loader->Synchronize();
2109 }
2110
2111 fKineDataLoader->Synchronize();
2112 fTrackRefsDataLoader->Synchronize();
bacbe0fd 2113
2114 TFile* file = gROOT->GetFile( gSystem->ConcatFileName( fUnixDirName.Data(), fgkDefaultTriggerFileName.Data() ) ) ;
2115 if( file ) file->Flush();
f0f6f856 2116
504b172d 2117 if (fGAFile) fGAFile->Flush();
f0f6f856 2118}
88cb7938 2119/*****************************************************************************/
2120/*****************************************************************************/