ANALYSIS/AliAnalysisTask.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 // Author: Andrei Gheata, 31/05/2006
  18
  19 //==============================================================================
  20 //   AliAnalysysTask - Class representing a basic analysis task. Any
  21 // user-defined task should derive from it and implement the Exec() virtual
  22 // method.
  23 //==============================================================================
  24 //
  25 // A specific user analysis task have to derive from this class. The list of
  26 // specific input and output slots have to be defined in the derived class ctor:
  27 //
  28 //   UserTask::UserTask(name, title)
  29 //   {
  30 //      DefineInput(0, TTree::Class());
  31 //      DefineInput(1, TH1::Class());
  32 //      ...
  33 //      DefineOutput(0, TTree::Class());
  34 //      DefineOutput(1, MyObject::Class());
  35 //      ...
  36 //   }
  37 //
  38 // An existing data contaner (AliAnalysisDataContainer) can be connected to the
  39 // input/output slots of an analysis task. Containers should not be defined and
  40 // connected by the derived analysis task, but from the level of AliAnalysisManager:
  41 //
  42 //   AliAnalysisManager::ConnectInput(AliAnalysisTask *task, Int_t islot,
  43 //                                   AliAnalysisDataContainer *cont)
  44 //   AliAnalysisManager::ConnectOutput(AliAnalysisTask *task, Int_t islot,
  45 //                                    AliAnalysisDataContainer *cont)
  46 // To connect a slot to a data container, the data types declared by both must
  47 // match.
  48 //
  49 // The method ConnectInputData() has to be overloaded by the derived class in order to
  50 // set the branch address or connect to a branch address in case the input
  51 // slots are connected to trees.
  52 // Example:
  53 // MyAnalysisTask::ConnectInputData(Option_t *)
  54 // {
  55 //  // One should first check if the branch address was taken by some other task
  56 //    char ** address = (char **)GetBranchAddress(0, "ESD");
  57 //    if (address) {
  58 //      fESD = (AliESD*)(*address);
  59 //    } else {
  60 //      fESD = new AliESD();
  61 //      SetBranchAddress(0, "ESD", &fESD);
  62 //    }
  63 // }
  64 //
  65 // The method LocalInit() may be implemented to call locally (on the client)
  66 // all initialization methods of the class. It is not mandatory and was created
  67 // in order to minimize the complexity and readability of the analysis macro.
  68 // DO NOT create in this method the histigrams or task output objects that will
  69 // go in the task output containers. Use CreateOutputObjects for that.
  70 //
  71 // The method CreateOutputObjects() has to be implemented an will contain the
  72 // objects that should be created only once per session (e.g. output
  73 // histograms)
  74 //
  75 // void MyAnalysisTask::CreateOutputObjects()
  76 //{
  77   // create histograms
  78 //  fhPt = new TH1F("fhPt","This is the Pt distribution",15,0.1,3.1);
  79 //  fhPt->SetStats(kTRUE);
  80 //  fhPt->GetXaxis()->SetTitle("P_{T} [GeV]");
  81 //  fhPt->GetYaxis()->SetTitle("#frac{dN}{dP_{T}}");
  82 //  fhPt->GetXaxis()->SetTitleColor(1);
  83 //  fhPt->SetMarkerStyle(kFullCircle);
  84 // }
  85 //
  86 // The method Terminate() will be called by the framework once at the end of
  87 // data processing. Overload this if needed. DO NOT ASSUME that the pointers
  88 // to histograms defined in  CreateOutputObjects() are valid, since this is
  89 // not true in case of PROOF. Restore the pointer values like:
  90 //
  91 //void MyAnalysisTask::Terminate(Option_t *)
  92 //{
  93 //  fhPt = (TH1F*)GetOutputData(0);
  94 // ...
  95 //}
  96
  97 //
  98 //==============================================================================
  99
 100 #include <Riostream.h>
 101 #include <TFile.h>
 102 #include <TClass.h>
 103 #include <TCollection.h>
 104 #include <TTree.h>
 105 #include <TROOT.h>
 106
 107 #include "AliAnalysisTask.h"
 108 #include "AliAnalysisDataSlot.h"
 109 #include "AliAnalysisDataContainer.h"
 110 #include "AliAnalysisManager.h"
 111
 112 ClassImp(AliAnalysisTask)
 113
 114 //______________________________________________________________________________
 115 AliAnalysisTask::AliAnalysisTask()
 116                 :fReady(kFALSE),
 117                  fInitialized(kFALSE),
 118                  fNinputs(0),
 119                  fNoutputs(0),
 120                  fOutputReady(NULL),
 121                  fPublishedData(NULL),
 122                  fInputs(NULL),
 123                  fOutputs(NULL),
 124                  fBranchNames()
 125 {
 126 // Default constructor.
 127 }
 128
 129 //______________________________________________________________________________
 130 AliAnalysisTask::AliAnalysisTask(const char *name, const char *title)
 131                 :TTask(name,title),
 132                  fReady(kFALSE),
 133                  fInitialized(kFALSE),
 134                  fNinputs(0),
 135                  fNoutputs(0),
 136                  fOutputReady(NULL),
 137                  fPublishedData(NULL),
 138                  fInputs(NULL),
 139                  fOutputs(NULL),
 140                  fBranchNames()
 141 {
 142 // Constructor.
 143    fInputs      = new TObjArray(2);
 144    fOutputs     = new TObjArray(2);
 145 }
 146
 147 //______________________________________________________________________________
 148 AliAnalysisTask::AliAnalysisTask(const AliAnalysisTask &task)
 149                 :TTask(task),
 150                  fReady(task.fReady),
 151                  fInitialized(task.fInitialized),
 152                  fNinputs(task.fNinputs),
 153                  fNoutputs(task.fNoutputs),
 154                  fOutputReady(NULL),
 155                  fPublishedData(NULL),
 156                  fInputs(NULL),
 157                  fOutputs(NULL),
 158                  fBranchNames(task.fBranchNames)
 159 {
 160 // Copy ctor.
 161    fInputs      = new TObjArray((fNinputs)?fNinputs:2);
 162    fOutputs     = new TObjArray((fNoutputs)?fNoutputs:2);
 163    fPublishedData = 0;
 164    Int_t i;
 165    for (i=0; i<fNinputs; i++) fInputs->AddAt(task.GetInputSlot(i),i);
 166    fOutputReady = new Bool_t[(fNoutputs)?fNoutputs:2];
 167    for (i=0; i<fNoutputs; i++) {
 168       fOutputReady[i] = IsOutputReady(i);
 169       fOutputs->AddAt(task.GetOutputSlot(i),i);
 170    }
 171 }
 172
 173 //______________________________________________________________________________
 174 AliAnalysisTask::~AliAnalysisTask()
 175 {
 176 // Dtor.
 177    if (fTasks) fTasks->Clear();
 178    if (fInputs)  {fInputs->Delete(); delete fInputs;}
 179    if (fOutputs) {fOutputs->Delete(); delete fOutputs;}
 180 }
 181
 182 //______________________________________________________________________________
 183 AliAnalysisTask& AliAnalysisTask::operator=(const AliAnalysisTask& task)
 184 {
 185 // Assignment
 186    if (&task == this) return *this;
 187    TTask::operator=(task);
 188    fReady       = task.IsReady();
 189    fInitialized = task.IsInitialized();
 190    fNinputs     = task.GetNinputs();
 191    fNoutputs    = task.GetNoutputs();
 192    fInputs      = new TObjArray((fNinputs)?fNinputs:2);
 193    fOutputs     = new TObjArray((fNoutputs)?fNoutputs:2);
 194    fPublishedData = 0;
 195    Int_t i;
 196    for (i=0; i<fNinputs; i++) fInputs->AddAt(new AliAnalysisDataSlot(*task.GetInputSlot(i)),i);
 197    fOutputReady = new Bool_t[(fNoutputs)?fNoutputs:2];
 198    for (i=0; i<fNoutputs; i++) {
 199       fOutputReady[i] = IsOutputReady(i);
 200       fOutputs->AddAt(new AliAnalysisDataSlot(*task.GetOutputSlot(i)),i);
 201    }
 202    fBranchNames = task.fBranchNames;
 203    return *this;
 204 }
 205
 206 //______________________________________________________________________________
 207 Bool_t AliAnalysisTask::AreSlotsConnected()
 208 {
 209 // Check if all input/output slots are connected. If this is the case fReady=true
 210    fReady = kFALSE;
 211    if (!fNinputs || !fNoutputs) return kFALSE;
 212    Int_t i;
 213    AliAnalysisDataSlot *slot;
 214    for (i=0; i<fNinputs; i++) {
 215       slot = (AliAnalysisDataSlot*)fInputs->At(i);
 216       if (!slot) {
 217               Error("AreSlotsConnected", "Input slot %d of task %s not defined !",i,GetName());
 218          return kFALSE;
 219       }
 220       if (!slot->IsConnected()) return kFALSE;
 221    }
 222    for (i=0; i<fNoutputs; i++) {
 223       slot = (AliAnalysisDataSlot*)fOutputs->At(i);
 224       if (!slot) {
 225          Error("AreSlotsConnected", "Output slot %d of task %s not defined !",i,GetName());
 226          return kFALSE;
 227       }
 228       if (!slot->IsConnected()) return kFALSE;
 229    }
 230    fReady = kTRUE;
 231    return kTRUE;
 232 }
 233
 234 //______________________________________________________________________________
 235 void AliAnalysisTask::CheckNotify(Bool_t init)
 236 {
 237 // Check if data is available from all inputs. Change the status of the task
 238 // accordingly. This method is called automatically for all tasks connected
 239 // to a container where the data was published.
 240    if (init) fInitialized = kFALSE;
 241    Bool_t single_shot = IsPostEventLoop();
 242    AliAnalysisDataContainer *cinput;
 243    for (Int_t islot=0; islot<fNinputs; islot++) {
 244       cinput = GetInputSlot(islot)->GetContainer();
 245       if (!cinput->GetData() || (single_shot && !cinput->IsPostEventLoop())) {
 246          SetActive(kFALSE);
 247          return;
 248       }
 249    }
 250    SetActive(kTRUE);
 251    if (fInitialized) return;
 252    TDirectory *cursav = gDirectory;
 253    ConnectInputData();
 254    if (cursav) cursav->cd();
 255    fInitialized = kTRUE;
 256 }
 257
 258 //______________________________________________________________________________
 259 Bool_t AliAnalysisTask::CheckPostData() const
 260 {
 261 // Checks if data was posted to all outputs defined by the task. If task does
 262 // not have output slots this returns always kTRUE.
 263    Bool_t dataPosted = kTRUE;
 264    AliAnalysisDataContainer *coutput;
 265    AliAnalysisManager *mgr = AliAnalysisManager::GetAnalysisManager();
 266    for (Int_t islot=0; islot<fNoutputs; islot++) {
 267       coutput = GetOutputSlot(islot)->GetContainer();
 268       if (!coutput) continue;
 269       if (coutput->IsExchange() || !mgr->GetOutputs()->FindObject(coutput) || coutput==mgr->GetCommonOutputContainer()) continue;
 270       if (!coutput->GetData()) {
 271          Error("CheckPostData", "Data not posted for slot #%d of task %s (%s)",
 272                islot, GetName(), ClassName());
 273          dataPosted = kFALSE;
 274       }
 275    }
 276    CheckOwnership();
 277    return dataPosted;
 278 }
 279
 280 //______________________________________________________________________________
 281 Bool_t AliAnalysisTask::CheckOwnership() const
 282 {
 283 // Check ownership of containers posted on output slots (1 level only)
 284    TObject *outdata;
 285    for (Int_t islot=0; islot<fNoutputs; islot++) {
 286       outdata = GetOutputData(islot);
 287       if (outdata && outdata->InheritsFrom(TCollection::Class())) {
 288          TCollection *coll = (TCollection*)outdata;
 289          if (!coll->IsOwner()) {
 290             Error("CheckOwnership","####### IMPORTANT! ####### \n\n\n\
 291                 Task %s (%s) posts a container that is not owner at output #%d. This may apply for other embedded containers. \n\n\
 292                 ####### FIX YOUR CODE, THIS WILL PRODUCE A FATAL ERROR IN FUTURE! ##########", GetName(), ClassName(), islot);
 293              return kFALSE;
 294          }
 295       }
 296    }
 297    return kTRUE;
 298 }
 299
 300 //______________________________________________________________________________
 301 Bool_t AliAnalysisTask::ConnectInput(Int_t islot, AliAnalysisDataContainer *cont)
 302 {
 303 // Connect an input slot to a data container.
 304    AliAnalysisDataSlot *input = GetInputSlot(islot);
 305    if (!input) {
 306       Error("ConnectInput","Input slot %i not defined for analysis task %s", islot, GetName());
 307       return kFALSE;
 308    }
 309    // Check type matching
 310    if (!input->GetType()->InheritsFrom(cont->GetType())) {
 311       Error("ConnectInput","Data type %s for input %i of task %s not matching container %s of type %s",input->GetType()->GetName(), islot, GetName(), cont->GetName(), cont->GetType()->GetName());
 312       return kFALSE;
 313    }
 314    // Connect the slot to the container as input
 315    if (!input->ConnectContainer(cont)) return kFALSE;
 316    // Add this to the list of container consumers
 317    cont->AddConsumer(this, islot);
 318    AreSlotsConnected();
 319    return kTRUE;
 320 }
 321
 322 //______________________________________________________________________________
 323 Bool_t AliAnalysisTask::ConnectOutput(Int_t islot, AliAnalysisDataContainer *cont)
 324 {
 325 // Connect an output slot to a data container.
 326    AliAnalysisDataSlot *output = GetOutputSlot(islot);
 327    if (!output) {
 328       Error("ConnectOutput","Output slot %i not defined for analysis task %s", islot, GetName());
 329       return kFALSE;
 330    }
 331    // Check type matching
 332    if (!output->GetType()->InheritsFrom(cont->GetType())) {
 333       Error("ConnectOutput","Data type %s for output %i of task %s not matching container %s of type %s",output->GetType()->GetName(), islot, GetName(), cont->GetName(), cont->GetType()->GetName());
 334       return kFALSE;
 335    }
 336    // Connect the slot to the container as output
 337    if (!output->ConnectContainer(cont)) return kFALSE;
 338    // Set event loop type the same as for the task
 339    cont->SetPostEventLoop(IsPostEventLoop());
 340    // Declare this as the data producer
 341    cont->SetProducer(this, islot);
 342    AreSlotsConnected();
 343    return kTRUE;
 344 }
 345
 346 //______________________________________________________________________________
 347 void AliAnalysisTask::DefineInput(Int_t islot, TClass *type)
 348 {
 349 // Define an input slot and its type.
 350    AliAnalysisDataSlot *input = new AliAnalysisDataSlot(type, this);
 351    if (fNinputs<islot+1) fNinputs = islot+1;
 352    fInputs->AddAtAndExpand(input, islot);
 353 }
 354
 355 //______________________________________________________________________________
 356 void AliAnalysisTask::DefineOutput(Int_t islot, TClass *type)
 357 {
 358 // Define an output slot and its type.
 359    AliAnalysisDataSlot *output = new AliAnalysisDataSlot(type, this);
 360    if (fNoutputs<islot+1) {
 361       fNoutputs = islot+1;
 362       if (fOutputReady) delete [] fOutputReady;
 363       fOutputReady = new Bool_t[fNoutputs];
 364       memset(fOutputReady, 0, fNoutputs*sizeof(Bool_t));
 365    }
 366    fOutputs->AddAtAndExpand(output, islot);
 367 }
 368
 369 //______________________________________________________________________________
 370 TClass *AliAnalysisTask::GetInputType(Int_t islot) const
 371 {
 372 // Retreive type of a given input slot.
 373    AliAnalysisDataSlot *input = GetInputSlot(islot);
 374    if (!input) {
 375       Error("GetInputType","Input slot %d not defined for analysis task %s", islot, GetName());
 376       return NULL;
 377    }
 378    return (input->GetType());
 379 }
 380
 381 //______________________________________________________________________________
 382 TClass *AliAnalysisTask::GetOutputType(Int_t islot) const
 383 {
 384 // Retreive type of a given output slot.
 385    AliAnalysisDataSlot *output = GetOutputSlot(islot);
 386    if (!output) {
 387       Error("GetOutputType","Output slot %d not defined for analysis task %s", islot, GetName());
 388       return NULL;
 389    }
 390    return (output->GetType());
 391 }
 392
 393 //______________________________________________________________________________
 394 TObject *AliAnalysisTask::GetInputData(Int_t islot) const
 395 {
 396 // Retreive input data for a slot if ready. Normally called by Exec() and
 397 // the object has to be statically cast to the appropriate type.
 398    AliAnalysisDataSlot *input = GetInputSlot(islot);
 399    if (!input) {
 400       Error("GetInputData","Input slot %d not defined for analysis task %s", islot, GetName());
 401       return NULL;
 402    }
 403    return (input->GetData());
 404 }
 405
 406 //______________________________________________________________________________
 407 TObject *AliAnalysisTask::GetOutputData(Int_t islot) const
 408 {
 409 // Retreive output data for a slot. Normally called in UserTask::Terminate to
 410 // get a valid pointer to data even in case of Proof.
 411    AliAnalysisDataSlot *output = GetOutputSlot(islot);
 412    if (!output) {
 413       Error("GetOutputData","Input slot %d not defined for analysis task %s", islot, GetName());
 414       return NULL;
 415    }
 416    return (output->GetData());
 417 }
 418
 419 //______________________________________________________________________________
 420 char *AliAnalysisTask::GetBranchAddress(Int_t islot, const char *branch) const
 421 {
 422 // Check if a branch with a given name from the specified input is connected
 423 // to some address. Call this in Init() before trying to call SetBranchAddress()
 424 // since the adress may be set by other task.
 425    return (char *)GetInputSlot(islot)->GetBranchAddress(branch);
 426 }
 427
 428 //______________________________________________________________________________
 429 Bool_t AliAnalysisTask::SetBranchAddress(Int_t islot, const char *branch, void *address) const
 430 {
 431 // Connect an object address to a branch of the specified input.
 432    return GetInputSlot(islot)->SetBranchAddress(branch, address);
 433 }
 434
 435 //______________________________________________________________________________
 436 void AliAnalysisTask::EnableBranch(Int_t islot, const char *bname) const
 437 {
 438 // Call this in ConnectInputData() to enable only the branches needed by this
 439 // task. "*" will enable everything.
 440    AliAnalysisDataSlot *input = GetInputSlot(islot);
 441    if (!input || !input->GetType()->InheritsFrom(TTree::Class())) {
 442       Error("EnableBranch", "Wrong slot type #%d for task %s: not TTree-derived type", islot, GetName());
 443       return;
 444    }
 445    TTree *tree = (TTree*)input->GetData();
 446    if (!strcmp(bname, "*")) {
 447       tree->SetBranchStatus("*",1);
 448       return;
 449    }
 450    AliAnalysisDataSlot::EnableBranch(bname, tree);
 451 }
 452
 453 //______________________________________________________________________________
 454 void AliAnalysisTask::FinishTaskOutput()
 455 {
 456 // Optional method that is called in SlaveTerminate phase.
 457 // Used for calling aditional methods just after the last event was processed ON
 458 // THE WORKING NODE. The call is made also in local case.
 459 // Do NOT delete output objects here since they will have to be sent for
 460 // merging in PROOF mode - use class destructor for cleanup.
 461 }
 462
 463 //______________________________________________________________________________
 464 void AliAnalysisTask::ConnectInputData(Option_t *)
 465 {
 466 // Overload and connect your branches here.
 467 }
 468
 469 //______________________________________________________________________________
 470 void AliAnalysisTask::LocalInit()
 471 {
 472 // The method LocalInit() may be implemented to call locally (on the client)
 473 // all initialization methods of the class. It is not mandatory and was created
 474 // in order to minimize the complexity and readability of the analysis macro.
 475 // DO NOT create in this method the histigrams or task output objects that will
 476 // go in the task output containers. Use CreateOutputObjects for that.
 477 }
 478
 479 //______________________________________________________________________________
 480 void AliAnalysisTask::CreateOutputObjects()
 481 {
 482 // Called once per task either in PROOF or local mode. Overload to put some
 483 // task initialization and/or create your output objects here.
 484 }
 485
 486 //______________________________________________________________________________
 487 TFile *AliAnalysisTask::OpenFile(Int_t iout, Option_t *option) const
 488 {
 489 // This method has to be called INSIDE the user redefined CreateOutputObjects
 490 // method, before creating each object corresponding to the output containers
 491 // that are to be written to a file. This need to be done in general for the big output
 492 // objects that may not fit memory during processing.
 493 // - 'option' is the file opening option.
 494 //=========================================================================
 495 // NOTE !: The method call will be ignored in PROOF mode, in which case the
 496 // results have to be streamed back to the client and written just before Terminate()
 497 //=========================================================================
 498 //
 499 // Example:
 500 // void MyAnaTask::CreateOutputObjects() {
 501 //    OpenFile(0);   // Will open the file for the object to be written at output #0
 502 //    fAOD = new TTree("AOD for D0toKPi");
 503 //    OpenFile(1);
 504 // now some histos that should go in the file of the second output container
 505 //    fHist1 = new TH1F("my quality check hist1",...);
 506 //    fHist2 = new TH2F("my quality check hist2",...);
 507 // }
 508
 509    if (iout<0 || iout>=fNoutputs) {
 510       Error("OpenFile", "No output slot for task %s with index %d", GetName(), iout);
 511       return NULL;
 512    }
 513    // Method delegated to the analysis manager (A.G. 02/11/09)
 514    AliAnalysisDataContainer *cont = GetOutputSlot(iout)->GetContainer();
 515    return AliAnalysisManager::OpenFile(cont, option);
 516 }
 517
 518 //______________________________________________________________________________
 519 Bool_t AliAnalysisTask::Notify()
 520 {
 521 // Overload this IF you need to treat input file change.
 522    return kTRUE;
 523 }
 524
 525 //______________________________________________________________________________
 526 Bool_t AliAnalysisTask::NotifyBinChange()
 527 {
 528 // Overload this IF you need to treat bin change in event mixing.
 529    return kTRUE;
 530 }
 531
 532 //______________________________________________________________________________
 533 void AliAnalysisTask::Terminate(Option_t *)
 534 {
 535 // Method called by the framework at the end of data processing.
 536 }
 537
 538 //______________________________________________________________________________
 539 Bool_t AliAnalysisTask::PostData(Int_t iout, TObject *data, Option_t *option)
 540 {
 541 // Post output data for a given ouput slot in the corresponding data container.
 542 // Published data becomes owned by the data container.
 543 // If option is specified, the container connected to the output slot must have
 544 // an associated file name defined. The option represents the method to open the file.
 545    fPublishedData = 0;
 546    AliAnalysisDataSlot *output = GetOutputSlot(iout);
 547    if (!output) {
 548       Error("PostData","Output slot %i not defined for analysis task %s", iout, GetName());
 549       return kFALSE;
 550    }
 551    if (!output->IsConnected()) {
 552       Error("PostData","Output slot %i of analysis task %s not connected to any data container", iout, GetName());
 553       return kFALSE;
 554    }
 555    if (!fOutputReady) {
 556       fOutputReady = new Bool_t[fNoutputs];
 557       memset(fOutputReady, 0, fNoutputs*sizeof(Bool_t));
 558    }
 559    fOutputReady[iout] = kTRUE;
 560    fPublishedData = data;
 561    return (output->GetContainer()->SetData(data, option));
 562 }
 563
 564 //______________________________________________________________________________
 565 void AliAnalysisTask::SetUsed(Bool_t flag)
 566 {
 567 // Set 'used' flag recursively to task and all daughter tasks.
 568    if (TestBit(kTaskUsed)==flag) return;
 569    TObject::SetBit(kTaskUsed,flag);
 570    Int_t nd = fTasks->GetSize();
 571    AliAnalysisTask *task;
 572    for (Int_t i=0; i<nd; i++) {
 573       task = (AliAnalysisTask*)fTasks->At(i);
 574       task->SetUsed(flag);
 575    }
 576 }
 577
 578 //______________________________________________________________________________
 579 void AliAnalysisTask::Reset()
 580 {
 581 // Clear activity flag. Reset data for exchange containers.
 582    fActive = kFALSE;
 583    fHasExecuted = kFALSE;
 584    // Call PostData(islot, 0) for all slots connected to exchange containers
 585    Int_t islot;
 586    AliAnalysisDataContainer *cont;
 587    for (islot=0; islot<fNinputs; islot++) {
 588       cont = GetInputSlot(islot)->GetContainer();
 589       if (cont && cont->IsExchange()) cont->Reset();
 590    }
 591 }
 592
 593 //______________________________________________________________________________
 594 Bool_t AliAnalysisTask::CheckCircularDeps()
 595 {
 596 // Check for illegal circular dependencies, e.g. a daughter task should not have
 597 // a hierarchical parent as subtask.
 598    if (IsChecked()) return kTRUE;
 599    SetChecked();
 600    TList *tasks = GetListOfTasks();
 601    Int_t ntasks = tasks->GetSize();
 602    AliAnalysisTask *task;
 603    for (Int_t i=0; i<ntasks; i++) {
 604       task = (AliAnalysisTask*)tasks->At(i);
 605       if (task->CheckCircularDeps()) return kTRUE;
 606    }
 607    SetChecked(kFALSE);
 608    return kFALSE;
 609 }
 610
 611 //______________________________________________________________________________
 612 Bool_t AliAnalysisTask::ProducersTouched() const
 613 {
 614 // Check if all producer containers are in the "touched" state.
 615    Int_t islot;
 616    AliAnalysisDataContainer *cont;
 617    for (islot=0; islot<fNinputs; islot++) {
 618       cont = GetInputSlot(islot)->GetContainer();
 619       // Simulate the data flow so the tasks are printed only when all inputs
 620       // are touched
 621       if (cont && !cont->IsTouched()) return kFALSE;
 622    }
 623    return kTRUE;
 624 }
 625
 626 //______________________________________________________________________________
 627 void AliAnalysisTask::PrintTask(Option_t *option, Int_t indent) const
 628 {
 629 // Print task info.
 630    Int_t islot;
 631    AliAnalysisDataContainer *cont;
 632    if (fActive) return;
 633    if (!ProducersTouched()) return;
 634    TString opt(option);
 635    opt.ToLower();
 636    Bool_t dep = (opt.Contains("dep"))?kTRUE:kFALSE;
 637    TString ind;
 638    for (Int_t i=0; i<indent; i++) ind += " ";
 639    if (!dep || (dep && IsChecked())) {
 640       printf("______________________________________________________________________________\n");
 641       printf("%s\n", Form("%stask: %s  ACTIVE=%i POST_LOOP=%i", ind.Data(), GetName(),IsActive(),IsPostEventLoop()));
 642       if (dep) const_cast<AliAnalysisTask*>(this)->SetChecked(kFALSE);
 643       else {
 644          for (islot=0; islot<fNinputs; islot++) {
 645             printf("%s", Form("%s   INPUT #%i: %s <- ",ind.Data(),islot, GetInputType(islot)->GetName()));
 646             cont = GetInputSlot(islot)->GetContainer();
 647             if (cont) printf(" [%s]\n", cont->GetName());
 648             else printf(" [NO CONTAINER]\n");
 649          }
 650          for (islot=0; islot<fNoutputs; islot++) {
 651             printf("%s", Form("%s   OUTPUT #%i: %s -> ",ind.Data(),islot, GetOutputType(islot)->GetName()));
 652             cont = GetOutputSlot(islot)->GetContainer();
 653             if (cont) printf(" [%s]\n", cont->GetName());
 654             else printf(" [NO CONTAINER]\n");
 655          }
 656       }
 657    }
 658    ((AliAnalysisTask*)this)->SetActive(kTRUE);
 659    PrintContainers(option, indent+3);
 660    if (!fBranchNames.IsNull()) printf("%sRequested branches:   %s\n", ind.Data(), fBranchNames.Data());
 661 }
 662
 663 //______________________________________________________________________________
 664 void AliAnalysisTask::PrintContainers(Option_t *option, Int_t indent) const
 665 {
 666 // Print containers info.
 667    AliAnalysisDataContainer *cont;
 668    TString ind;
 669    for (Int_t i=0; i<indent; i++) ind += " ";
 670    Int_t islot;
 671    for (islot=0; islot<fNoutputs; islot++) {
 672       cont = GetOutputSlot(islot)->GetContainer();
 673       if (cont) cont->PrintContainer(option, indent);
 674    }
 675 }
 676
 677 //______________________________________________________________________________
 678 void AliAnalysisTask::SetPostEventLoop(Bool_t flag)
 679 {
 680 // Set the task execution mode - run after event loop or not. All output
 681 // containers of this task will get the same type.
 682    TObject::SetBit(kTaskPostEventLoop,flag);
 683    AliAnalysisDataContainer *cont;
 684    Int_t islot;
 685    for (islot=0; islot<fNoutputs; islot++) {
 686       cont = GetOutputSlot(islot)->GetContainer();
 687       if (cont) cont->SetPostEventLoop(flag);
 688    }
 689 }
 690
 691 //______________________________________________________________________________
 692 void AliAnalysisTask::GetBranches(const char *type, TString &result) const
 693 {
 694 // Get the list of branches for a given type (ESD, AOD). The list of branches
 695 // requested by a task has to ve declared in the form:
 696 //   SetBranches("ESD:branch1,branch2,...,branchN AOD:branch1,branch2,...,branchM")
 697    result = "";
 698    if (fBranchNames.IsNull()) return;
 699    Int_t index1 = fBranchNames.Index(type);
 700    if (index1<0) return;
 701    index1 += 1+strlen(type);
 702    Int_t index2 = fBranchNames.Index(" ", index1);
 703    if (index2<0) index2 = fBranchNames.Length();
 704    result = fBranchNames(index1, index2-index1);
 705 }