1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
17 //_________________________________________________________________________
18 // Analysis task that executes the analysis classes
19 // that depend on the PartCorr frame, frame for Particle identification and correlations.
20 // Specially designed for calorimeters but also can be used for charged tracks
21 // Input of this task is a configuration file that contains all the settings of the analyis
23 // -- Author: Gustavo Conesa (INFN-LNF)
29 #include <TInterpreter.h>
30 //#include <Riostream.h>
33 #include "AliAnalysisTaskParticleCorrelation.h"
34 #include "AliAnaPartCorrMaker.h"
35 #include "AliCaloTrackReader.h"
38 ClassImp(AliAnalysisTaskParticleCorrelation)
40 ////////////////////////////////////////////////////////////////////////
41 AliAnalysisTaskParticleCorrelation::AliAnalysisTaskParticleCorrelation():
44 fOutputContainer(0x0),
47 // Default constructor
50 //_____________________________________________________
51 AliAnalysisTaskParticleCorrelation::AliAnalysisTaskParticleCorrelation(const char* name):
52 AliAnalysisTaskSE(name),
54 fOutputContainer(0x0),
57 // Default constructor
59 DefineOutput(1, TList::Class());
63 //_____________________________________________________
64 AliAnalysisTaskParticleCorrelation::~AliAnalysisTaskParticleCorrelation()
66 // Remove all pointers
69 fOutputContainer->Clear() ;
70 delete fOutputContainer ;
75 //_____________________________________________________
76 void AliAnalysisTaskParticleCorrelation::UserCreateOutputObjects()
78 // Create the output container
79 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::UserCreateOutputObjects() - Begin\n");
81 //Get list of aod arrays, add each aod array to analysis frame
82 TClonesArray * array = 0;
83 TList * list = fAna->GetAODBranchList();
84 for(Int_t iaod = 0; iaod < list->GetEntries(); iaod++){
85 array = (TClonesArray*) list->At(iaod);
86 AddAODBranch("TClonesArray", &array);
89 //Histograms container
91 fOutputContainer = fAna->GetOutputContainer();
92 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::UserCreateOutputObjects() - End\n");
95 //_____________________________________________________
96 void AliAnalysisTaskParticleCorrelation::Init()
100 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::Init() - Begin\n");
102 // Call configuration file if specified
104 if (fConfigName.Length()) {
105 printf("AliAnalysisTaskParticleCorrelation::Init() - ### Configuration file is %s.C ###\n", fConfigName.Data());
106 gROOT->LoadMacro(fConfigName+".C");
107 fAna = (AliAnaPartCorrMaker*) gInterpreter->ProcessLine("ConfigAnalysis()");
111 printf("AliAnalysisTaskParticleCorrelation::Init() - Analysis maker pointer not initialized, no analysis specified, STOP!\n");
115 // Add different generator particles to PDG Data Base
116 // to avoid problems when reading MC generator particles
117 AliPDG::AddParticlesToPdgDataBase();
119 // Initialise analysis
122 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::Init() - End\n");
127 //_____________________________________________________
128 void AliAnalysisTaskParticleCorrelation::UserExec(Option_t */*option*/)
130 // Execute analysis for current event
132 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::UserExec() - Begin\n");
134 //Get the type of data, check if type is correct
135 Int_t datatype = fAna->GetReader()->GetDataType();
136 if(datatype != AliCaloTrackReader::kESD && datatype != AliCaloTrackReader::kAOD &&
137 datatype != AliCaloTrackReader::kMC){
138 printf("AliAnalysisTaskParticleCorrelation::UserExec() - Wrong type of data\n");
142 fAna->GetReader()->SetInputOutputMCEvent(InputEvent(), AODEvent(), MCEvent());
145 fAna->ProcessEvent((Int_t) Entry(), CurrentFileName());
146 //printf("AliAnalysisTaskParticleCorrelation::Current Event %d; Current File Name : %s\n",(Int_t) Entry(), CurrentFileName());
147 if (DebugLevel() > 1) printf("AliAnalysisTaskParticleCorrelation::UserExec() - End\n");
149 PostData(1, fOutputContainer);
153 //_____________________________________________________
154 void AliAnalysisTaskParticleCorrelation::Terminate(Option_t */*option*/)
156 // Terminate analysis
159 // Get merged histograms from the output container
160 TList *outputList = (TList*)GetOutputData(1);
161 // Propagate histagrams to maker
162 fAna->Terminate(outputList);