coverity fixes
[u/mrichter/AliRoot.git] / PWGJE / UserTasks / AliAnalysisTaskJetHadronCorrelation.cxx
1 //#include <string.h>
2 //#include <TStyle.h>
3 #include <list>
4 #include <string>
5
6 #include "TTree.h"
7 #include "TCanvas.h"
8 #include "AliAnalysisTask.h"
9 #include "AliInputEventHandler.h"
10 #include "AliESDtrack.h"
11 #include "AliAODVertex.h"
12 #include "AliAODCluster.h"
13
14 #include <TROOT.h>
15 #include <TRandom.h>
16 #include <TSystem.h>
17 #include <TInterpreter.h>
18 #include <TChain.h>
19 #include <TFile.h>
20 #include <TKey.h>
21 #include <TH1F.h>
22 #include <TH2F.h>
23 #include <TH3F.h>
24 #include <TProfile.h>
25 #include <TList.h>
26 #include <TLorentzVector.h>
27 #include <TClonesArray.h>
28 #include <TRefArray.h>
29
30 #include "TDatabasePDG.h"
31 #include "AliAnalysisManager.h"
32 #include "AliJetFinder.h"
33 #include "AliJetHeader.h"
34 #include "AliJetReader.h"
35 #include "AliJetReaderHeader.h"
36 #include "AliUA1JetHeaderV1.h"
37 #include "AliSISConeJetHeader.h"
38 #include "AliESDEvent.h"
39 #include "AliAODEvent.h"
40 #include "AliAODHandler.h"
41 #include "AliAODInputHandler.h"
42 #include "AliAODTrack.h"
43 #include "AliAODJet.h"
44 #include "AliAODJetEventBackground.h"
45 #include "AliMCParticle.h"
46 #include "AliAODMCParticle.h"
47 #include "AliMCEventHandler.h"
48 #include "AliMCEvent.h"
49 #include "AliStack.h"
50
51 #include "AliAODHeader.h"
52 #include "AliAODMCHeader.h"
53 //#include "AliGenPythiaEventHeader.h"
54 #include "AliJetKineReaderHeader.h"
55 #include "AliGenCocktailEventHeader.h"
56 #include "AliInputEventHandler.h"
57 #include "AliGenEventHeader.h"
58 #include "AliGenDPMjetEventHeader.h"
59
60
61
62 #include "AliAnalysisTaskJetHadronCorrelation.h"
63 #include "AliAnalysisTaskPhiCorrelations.h"
64 //#include "AliAnalysisHelperJetTasks.h"
65 #include "AliPWG4HighPtQAMC.h"
66
67
68
69
70 using std::cout;
71 using std::endl;
72
73 ClassImp(AliAnalysisTaskJetHadronCorrelation)
74
75                                 //________________________________________________________________________
76                                 AliAnalysisTaskJetHadronCorrelation::AliAnalysisTaskJetHadronCorrelation():
77                                                                 AliAnalysisTaskSE(),
78                                                                 fUseAODInput(kFALSE),
79                                                                 fFillAOD(kFALSE),
80                                                                 fJetBranch("jets"),
81                                                                 fNonStdFile(""),
82                                                                 fAODIn(0x0),
83                                                                 fAODOut(0x0),
84                                                                 fAODExtension(0x0),
85                                                                 JFAlg("ANTIKT"),         
86                                                                 Radius(0.4),
87                                                                 Filtermask(272),
88                                                                 BackM(0),
89                                                                 TrackPtcut(0.15),
90                                                                 SkipCone(0),
91                                                                 IsMC(kTRUE),
92                 fxsec(0.),
93                 ftrial(1.),
94                 fJetRecEtaWindow(0.5),       // eta window for rec jets
95                 fMinJetPt(10), 
96                                                                 fHistList(0x0), // Output list
97                                                                 fIfiles(0),
98                                                                 fH1Events(0x0),
99                                                                 fH1Xsec(0x0),
100                                                                 fH1Trials(0x0),
101                                                                 fH1JetAKT04_pt                (0x0),
102                                                                 fH1leadJetAKT04_pt            (0x0),
103                                                                 fH1leadJetAKT04_pt_dijet      (0x0),
104                                                                 fH1subJetAKT04_pt_dijet       (0x0),
105                                                                 fH2JetsJetAKT04_dphi          (0x0),
106                                                                 fH2JetsJetAKT04_deta          (0x0),
107                                                                 fH2JetsJetAKT04_Aj            (0x0),
108                                                                 fH2JetsJetAKT04_pt            (0x0),
109                                                                 fH1JetMCAKT04_pt              (0x0),
110                                                                 fH1leadJetMCAKT04_pt          (0x0),
111                                                                 fH1leadJetMCAKT04_pt_dijet    (0x0),
112                                                                 fH1subJetMCAKT04_pt_dijet     (0x0),
113                                                                 fH2JetsJetMCAKT04_dphi        (0x0),
114                                                                 fH2JetsJetMCAKT04_deta        (0x0),
115                                                                 fH2JetsJetMCAKT04_Aj          (0x0),
116                                                                 fH2JetsJetMCAKT04_pt          (0x0)
117
118
119 {
120
121
122                                 if(IsMC){
123                                                                 for(int j=0;j<5;j++){
124                                                                                                 fH1AKT04_ndiJ_ediv[j]=0;
125                                                                                                 fH1leadJetMCAKT04_dphiResolution          [j]=0;
126                                                                                                 fH1subJetMCAKT04_dphiResolution           [j]=0;
127                                                                                                 for(int k=0;k<5;k++){
128                                                                                                                                 fH1JetHadronAKT04_dphi_ediv               [j][k]=0;
129                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv     [j][k]=0;
130                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv    [j][k]=0;
131                                                                                                 }
132                                                                 }
133                                 }else{
134                                                                 for(int j=0;j<5;j++){
135                                                                                                 fH1AKT04_ndiJ_ediv[j]=0;
136                                                                                                 for(int k=0;k<5;k++){
137                                                                                                                                 fH1JetHadronAKT04_dphi_ediv             [j][k]=0;
138                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv   [j][k]=0;
139                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv  [j][k]=0;
140                                                                                                 }
141                                                                 }
142                                 }
143
144
145                                 // Default constructor
146 }
147
148 //________________________________________________________________________
149 AliAnalysisTaskJetHadronCorrelation::AliAnalysisTaskJetHadronCorrelation(const char *name):
150                                 AliAnalysisTaskSE(name),
151                                 fUseAODInput(kFALSE),
152                                 fFillAOD(kFALSE),
153                                 fJetBranch("jets"),
154                                 fNonStdFile(""),
155                                 fAODIn(0x0), 
156                                 fAODOut(0x0), 
157                                 fAODExtension(0x0),
158                                 JFAlg("ANTIKT"),         
159                                 Radius(0.4),
160                                 Filtermask(272),
161                                 BackM(0),
162                                 TrackPtcut(0.15),
163                                 SkipCone(0),
164                                 IsMC(kTRUE),
165                                 fxsec(0.),
166                                 ftrial(1.),
167                                 fJetRecEtaWindow(0.5),       // eta window for rec jets
168                                 fMinJetPt(10), 
169                                 fHistList(0x0), // Output list
170                                 fIfiles(0),
171
172                                 fH1Events(0x0),
173                                 fH1Xsec(0x0),
174                                 fH1Trials(0x0),
175                                 fH1JetAKT04_pt                (0x0),
176                                 fH1leadJetAKT04_pt            (0x0),
177                                 fH1leadJetAKT04_pt_dijet      (0x0),
178                                 fH1subJetAKT04_pt_dijet       (0x0),
179                                 fH2JetsJetAKT04_dphi          (0x0),
180                                 fH2JetsJetAKT04_deta          (0x0),
181                                 fH2JetsJetAKT04_Aj            (0x0),
182                                 fH2JetsJetAKT04_pt            (0x0),
183                                 fH1JetMCAKT04_pt              (0x0),
184                                 fH1leadJetMCAKT04_pt          (0x0),
185                                 fH1leadJetMCAKT04_pt_dijet    (0x0),
186                                 fH1subJetMCAKT04_pt_dijet     (0x0),
187                                 fH2JetsJetMCAKT04_dphi        (0x0),
188                                 fH2JetsJetMCAKT04_deta        (0x0),
189                                 fH2JetsJetMCAKT04_Aj          (0x0),
190                                 fH2JetsJetMCAKT04_pt          (0x0)
191
192
193 {
194                                 if(IsMC){
195                                                                 for(int j=0;j<5;j++){
196                                                                                                 fH1AKT04_ndiJ_ediv[j]=0;
197                                                                                                 fH1leadJetMCAKT04_dphiResolution          [j]=0;
198                                                                                                 fH1subJetMCAKT04_dphiResolution           [j]=0;
199                                                                                                 for(int k=0;k<5;k++){
200                                                                                                                                 fH1JetHadronAKT04_dphi_ediv               [j][k]=0;
201                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv     [j][k]=0;
202                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv    [j][k]=0;
203                                                                                                 }
204                                                                 }
205                                 }else{
206                                                                 for(int j=0;j<5;j++){
207                                                                                                 fH1AKT04_ndiJ_ediv[j]=0;
208                                                                                                 for(int k=0;k<5;k++){
209                                                                                                                                 fH1JetHadronAKT04_dphi_ediv             [j][k]=0;
210                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv   [j][k]=0;
211                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv  [j][k]=0;
212                                                                                                 }
213                                                                 }
214                                 }
215
216                                 // Default constructor
217
218                                 // Constructor
219
220                                 // Define input and output slots here
221                                 // Input slot #0 works with a TChain
222                                 DefineInput(0, TChain::Class());
223                                 // Output slot #0 id reserved by the base class for AOD
224                                 // Output slot #1 writes into a TH1 container
225                                 DefineOutput(1, TList::Class());
226
227 }
228
229 //________________________________________________________________________
230 void AliAnalysisTaskJetHadronCorrelation::UserCreateOutputObjects()
231 {
232                                 // Create histograms
233                                 // Called once
234
235
236                                 fHistList = new TList();fHistList->SetOwner(kTRUE); cout<<"TList is created for output "<<endl;
237                                 //if (!fHistList){ fHistList = new TList();fHistList->SetOwner(kTRUE); cout<<"TList is created for output "<<endl;}
238
239                                 Bool_t oldStatus = TH1::AddDirectoryStatus();
240                                 TH1::AddDirectory(kFALSE);
241
242                                 Float_t pi=TMath::Pi();
243                                 //gStyle->SetPalette(1);
244
245
246                                 char *histname;
247                                 if(IsMC){
248                                                                 fH1Xsec                         = new TProfile("Xsec"               ,"Xsec"                   ,1,0,1);
249                                                                 fH1Trials                       = new TH1F    ("Trials"             ,"Trials"                 ,1,0,1);
250                                                                 fH1JetMCAKT04_pt                = new TH1F("JetMCAKT04_pt"          ,"JetMCAKT04_pt"          ,400,0,400);
251                                                                 fH1leadJetMCAKT04_pt            = new TH1F("leadJetMCAKT04_pt"      ,"leadJetMCAKT04_pt"      ,400,0,400);
252                                                                 fH1leadJetMCAKT04_pt_dijet      = new TH1F("leadJetMCAKT04_pt_dijet","leadJetMCAKT04_pt_dijet",400,0,400);
253                                                                 fH1subJetMCAKT04_pt_dijet       = new TH1F("subJetMCAKT04_pt_dijet" ,"subJetMCAKT04_pt_dijet" ,400,0,400);
254                                                                 fH1JetAKT04_pt                  = new TH1F("JetAKT04_pt"            ,"JetAKT04_pt"            ,400,0,400);
255                                                                 fH1leadJetAKT04_pt              = new TH1F("leadJetAKT04_pt"        ,"leadJetAKT04_pt"        ,400,0,400);
256                                                                 fH1leadJetAKT04_pt_dijet        = new TH1F("leadJetAKT04_pt_dijet"  ,"leadJetAKT04_pt_dijet"  ,400,0,400);
257                                                                 fH1subJetAKT04_pt_dijet         = new TH1F("subJetAKT04_pt_dijet"   ,"subJetAKT04_pt_dijet"   ,400,0,400);
258                                                                 histname = Form("JetsJetMCAKT04_dphi");
259                                                                 fH2JetsJetMCAKT04_dphi          = new TH2F(histname,histname,200,0,400,100,-2*pi,2*pi);
260                                                                 histname = Form("JetsJetMCAKT04_deta");
261                                                                 fH2JetsJetMCAKT04_deta          = new TH2F(histname,histname,200,0,400,100,-1.5,1.5);
262                                                                 histname = Form("JetsJetMCAKT04_Aj");
263                                                                 fH2JetsJetMCAKT04_Aj            = new TH2F(histname,histname,200,0,400,100,0,1.2);
264                                                                 histname = Form("JetsJetMCAKT04_pt");
265                                                                 fH2JetsJetMCAKT04_pt            = new TH2F(histname,histname,200,0,400,200,0,400);
266                                                                 histname = Form("JetsJetAKT04_dphi");
267                                                                 fH2JetsJetAKT04_dphi            = new TH2F(histname,histname,200,0,400,100,-2*pi,2*pi);
268                                                                 histname = Form("JetsJetAKT04_deta");
269                                                                 fH2JetsJetAKT04_deta            = new TH2F(histname,histname,200,0,400,100,-1.5,1.5);
270                                                                 histname = Form("JetsJetAKT04_Aj");
271                                                                 fH2JetsJetAKT04_Aj              = new TH2F(histname,histname,200,0,400,100,0,1.2);
272                                                                 histname = Form("JetsJetAKT04_pt");
273                                                                 fH2JetsJetAKT04_pt              = new TH2F(histname,histname,200,0,400,200,0,400);
274                                                                 for(int j=0;j<5;j++){
275                                                                                                 histname = Form("AKT04_ndiJ_ediv%d",j);
276                                                                                                 fH1AKT04_ndiJ_ediv[j]= new TH1F(histname,histname,1,1,2);
277                                                                                                 histname = Form("leadJetMCAKT04_dphiResolution%d",j);
278                                                                                                 fH1leadJetMCAKT04_dphiResolution[j] = new TH1F(histname,histname,200,-2*pi,2*pi);
279                                                                                                 histname = Form("subJetMCAKT04_dphiResolution%d",j);
280                                                                                                 fH1subJetMCAKT04_dphiResolution[j] = new TH1F(histname,histname,200,-2*pi,2*pi);
281                                                                                                 for(int k=0;k<5;k++){
282                                                                                                                                 histname = Form("JetHadronAKT04_dphi_ediv%d%d",j,k);
283                                                                                                                                 fH1JetHadronAKT04_dphi_ediv             [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
284                                                                                                                                 histname = Form("JetHadronAKT04_dphi_tptweight_ediv%d%d",j,k);
285                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv   [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
286                                                                                                                                 histname = Form("JetHadronAKT04_dphi_tJptweight_ediv%d%d",j,k);
287                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv  [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
288                                                                                                 }
289                                                                 }
290
291                                                                 fHistList->Add(fH1Xsec);
292                                                                 fHistList->Add(fH1Trials);
293                                                                 fHistList->Add(fH1JetMCAKT04_pt          );
294                                                                 fHistList->Add(fH1leadJetMCAKT04_pt      );
295                                                                 fHistList->Add(fH1leadJetMCAKT04_pt_dijet);
296                                                                 fHistList->Add(fH1subJetMCAKT04_pt_dijet );
297                                                                 fHistList->Add(fH1JetAKT04_pt            );
298                                                                 fHistList->Add(fH1leadJetAKT04_pt        );
299                                                                 fHistList->Add(fH1leadJetAKT04_pt_dijet  );
300                                                                 fHistList->Add(fH1subJetAKT04_pt_dijet   );
301                                                                 fHistList->Add(fH2JetsJetMCAKT04_dphi);
302                                                                 fHistList->Add(fH2JetsJetMCAKT04_deta);
303                                                                 fHistList->Add(fH2JetsJetMCAKT04_Aj  );
304                                                                 fHistList->Add(fH2JetsJetMCAKT04_pt  );
305                                                                 fHistList->Add(fH2JetsJetAKT04_dphi  );
306                                                                 fHistList->Add(fH2JetsJetAKT04_deta  );
307                                                                 fHistList->Add(fH2JetsJetAKT04_Aj    );
308                                                                 fHistList->Add(fH2JetsJetAKT04_pt    );
309                                                                 for(int j=0;j<5;j++){
310                                                                                                 fHistList->Add(fH1AKT04_ndiJ_ediv    [j]);
311                                                                                                 fHistList->Add(fH1leadJetMCAKT04_dphiResolution          [j]);
312                                                                                                 fHistList->Add(fH1subJetMCAKT04_dphiResolution           [j]);
313                                                                                                 for(int k=0;k<5;k++){
314                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_ediv               [j][k]);
315                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_tptweight_ediv     [j][k]);
316                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_tJptweight_ediv    [j][k]);
317                                                                                                 }
318                                                                 }
319                                 }
320                                 else{
321                                                                 fH1Events                     = new TH1F("Events"                  ,"Events"                  ,1,0,1);
322                                                                 fH1JetAKT04_pt                = new TH1F("JetAKT04_pt"             ,"JetAKT04_pt"             ,400,0,400);
323                                                                 fH1leadJetAKT04_pt            = new TH1F("leadJetAKT04_pt"         ,"leadJetAKT04_pt"         ,400,0,400);
324                                                                 fH1leadJetAKT04_pt_dijet      = new TH1F("leadJetAKT04_pt_dijet"   ,"leadJetAKT04_pt_dijet"   ,400,0,400);
325                                                                 fH1subJetAKT04_pt_dijet       = new TH1F("subJetAKT04_pt_dijet"    ,"subJetAKT04_pt_dijet"    ,400,0,400);
326                                                                 histname = Form("JetsJetAKT04_dphi");
327                                                                 fH2JetsJetAKT04_dphi          = new TH2F(histname,histname,200,0,400,100,-2*pi,2*pi);
328                                                                 histname = Form("JetsJetAKT04_deta");
329                                                                 fH2JetsJetAKT04_deta          = new TH2F(histname,histname,200,0,400,100,-1.5,1.5);
330                                                                 histname = Form("JetsJetAKT04_Aj");
331                                                                 fH2JetsJetAKT04_Aj            = new TH2F(histname,histname,200,0,400,100,0,1.2);
332                                                                 histname = Form("JetsJetAKT04_pt");
333                                                                 fH2JetsJetAKT04_pt            = new TH2F(histname,histname,200,0,400,200,0,400);
334                                                                 for(int j=0;j<5;j++){
335                                                                                                 histname = Form("AKT04_ndiJ_ediv%d",j);
336                                                                                                 fH1AKT04_ndiJ_ediv[j]= new TH1F(histname,histname,1,1,2);
337                                                                                                 for(int k=0;k<5;k++){
338                                                                                                                                 histname = Form("JetHadronAKT04_dphi_ediv%d%d",j,k);
339                                                                                                                                 fH1JetHadronAKT04_dphi_ediv             [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
340                                                                                                                                 histname = Form("JetHadronAKT04_dphi_tptweight_ediv%d%d",j,k);
341                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv   [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
342                                                                                                                                 histname = Form("JetHadronAKT04_dphi_tJptweight_ediv%d%d",j,k);
343                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv  [j][k]= new TH1F(histname,histname,200,-1./2.*pi,3./2.*pi);
344                                                                                                 }
345                                                                 }
346                                                                 fHistList->Add(fH1Events);
347                                                                 fHistList->Add(fH1JetAKT04_pt          );
348                                                                 fHistList->Add(fH1leadJetAKT04_pt      );
349                                                                 fHistList->Add(fH1leadJetAKT04_pt_dijet);
350                                                                 fHistList->Add(fH1subJetAKT04_pt_dijet );
351                                                                 fHistList->Add(fH2JetsJetAKT04_dphi  );
352                                                                 fHistList->Add(fH2JetsJetAKT04_deta  );
353                                                                 fHistList->Add(fH2JetsJetAKT04_Aj    );
354                                                                 fHistList->Add(fH2JetsJetAKT04_pt    );
355                                                                 for(int j=0;j<5;j++){
356                                                                                                 fHistList->Add(fH1AKT04_ndiJ_ediv    [j]);
357                                                                                                 for(int k=0;k<5;k++){
358                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_ediv             [j][k]);
359                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_tptweight_ediv   [j][k]);
360                                                                                                                                 fHistList->Add(fH1JetHadronAKT04_dphi_tJptweight_ediv  [j][k]);
361                                                                                                 }
362                                                                 }
363                                 }
364
365
366
367                                 // =========== Switch on Sumw2 for all histos ===========
368                                 for (Int_t i=0; i<fHistList->GetEntries(); ++i) 
369                                 {
370                                                                 TH1 *h1 = dynamic_cast<TH1*>(fHistList->At(i));
371                                                                 if (h1)
372                                                                 {
373                                                                                                 h1->Sumw2();
374                                                                                                 continue;
375                                                                 }
376                                                                 THnSparse *hn = dynamic_cast<THnSparse*>(fHistList->At(i));
377                                                                 if(hn)hn->Sumw2();
378                                 }
379                                 TH1::AddDirectory(oldStatus);
380
381
382
383
384                                 PostData(1,fHistList);
385
386 }
387
388 //----------------------------------------------------------------------                                                 
389 void AliAnalysisTaskJetHadronCorrelation::Init()
390 {
391                                 // Initialization                                                                                                    
392                                 if (fDebug) printf("AnalysisTaskJetHadronCorrelation::Init() \n");
393
394 }
395
396 Bool_t AliAnalysisTaskJetHadronCorrelation::Notify()
397 {
398
399
400                                 fIfiles++;
401                                 fAODIn = dynamic_cast<AliAODEvent*>(InputEvent());
402                                 fAODOut = AODEvent();
403                                 if(fNonStdFile.Length()!=0){
404                                                                 AliAODHandler *aodH = dynamic_cast<AliAODHandler*>(AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler());
405                                                                 fAODExtension = (aodH?aodH->GetExtension(fNonStdFile.Data()):0);
406                                                                 if(fAODExtension){
407                                                                                                 if(fDebug>1)Printf("AODExtension found for %s ",fNonStdFile.Data());
408                                                                 }
409                                 }
410
411                                 TTree *tree = AliAnalysisManager::GetAnalysisManager()->GetTree();
412                                 fxsec=0;
413                                 ftrial=1;
414
415                                 if(tree){
416                                                                 TFile *curfile = tree->GetCurrentFile();
417                                                                 if(!curfile){
418                                                                                                 Error("Notify","No current file");
419                                                                                                 return kFALSE;
420                                                                 }
421
422
423                                                                 if(IsMC){
424                                                                                                 AliPWG4HighPtQAMC::PythiaInfoFromFile(curfile->GetName(),fxsec,ftrial);
425                                                                                                 //cout<<" Xsec "<<fxsec<<" trial "<<ftrial<<endl;
426                                                                                                 fH1Xsec  ->Fill(0.,fxsec);
427                                                                                                 fH1Trials->Fill(0.,ftrial);
428                                                                 }
429                                                                 //else    fH1Events->Fill(0.,ftrial);
430
431                                 }
432
433                                 printf("Reading File %s ",fInputHandler->GetTree()->GetCurrentFile()->GetName());
434                                 return kTRUE;
435 }
436 void AliAnalysisTaskJetHadronCorrelation::FinishTaskOutput()
437 {
438 }
439
440
441
442 //________________________________________________________________________
443 void AliAnalysisTaskJetHadronCorrelation::UserExec(Option_t *) 
444 {
445
446
447                                 // Main loop (called each event)
448                          // Execute analysis for current event
449
450                                 AliAODEvent *fAOD;
451                                 TObject* handler = AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler();
452                                 if( handler && handler->InheritsFrom("AliAODInputHandler") ) {
453                                                                 fAOD  =  ((AliAODInputHandler*)handler)->GetEvent();
454                                                                 if(fUseAODInput) fAODIn = fAOD;
455                                                                 if (fDebug > 1)  Printf("%s:%d AOD event from input", (char*)__FILE__,__LINE__);
456                                 }
457                                 else {
458                                                                 handler = AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler();
459                                                                 if( handler && handler->InheritsFrom("AliAODHandler") ) {
460                                                                                                 fAOD = ((AliAODHandler*)handler)->GetAOD();
461                                                                                                 fAODIn = fAOD;
462                                                                                                 if (fDebug > 1)  Printf("%s:%d AOD event from output", (char*)__FILE__,__LINE__);
463                                                                 }
464                                 }
465
466                                 if(!fAODIn && !fUseAODInput){ // case we have AOD in input & output and want jets from output
467                                                                 TObject* outHandler = AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler();
468                                                                 if( outHandler && outHandler->InheritsFrom("AliAODHandler") ) {
469                                                                                                 fAODIn = ((AliAODHandler*)outHandler)->GetAOD();
470                                                                                                 if (fDebug > 1)  Printf("%s:%d jets from output AOD", (char*)__FILE__,__LINE__);
471                                                                 }
472                                 }
473
474                                 //if(fNonStdFile.Length()!=0){
475                                 //                              // case we have an AOD extension - fetch the jets from the extended output
476                                 //                              AliAODHandler *aodH = dynamic_cast<AliAODHandler*>(AliAnalysisManager::GetAnalysisManager()->GetOutputEventHandler());
477                                 //                              fAODExtension = (aodH?aodH->GetExtension(fNonStdFile.Data()):0);
478                                 //                              if(!fAODExtension){
479                                 //                                                              if(fDebug>1)Printf("AODExtension not found for %s",fNonStdFile.Data());
480                                 //                              }
481                                 //}
482                                 //fAODIn = dynamic_cast<AliAODEvent*>(InputEvent());
483                                 if (!fAODIn) {
484                                                                 Printf("ERROR %s : fAODIn not available",(char*)__FILE__);
485                                                                 return;
486                                 }
487
488                                   AliInputEventHandler* inputHandler = (AliInputEventHandler*)
489                                                                             ((AliAnalysisManager::GetAnalysisManager())->GetInputEventHandler());
490                                 if(!(inputHandler->IsEventSelected() & AliVEvent::kMB)){
491                                                                 if (fDebug > 1 ) Printf(" Trigger Selection: event REJECTED ... ");
492                                                                 return;
493                                 }
494                                 if(!IsMC)fH1Events->Fill(0);
495
496                                 // start jet analysis
497
498                                 Double_t Jet_n  [20];
499                                 Double_t Jet_pt [20][1000];
500                                 Double_t Jet_eta[20][1000];
501                                 Double_t Jet_phi[20][1000];
502                                 Double_t subJet_n  [20];
503                                 Double_t subJet_pt [20][1000];
504                                 Double_t subJet_eta[20][1000];
505                                 Double_t subJet_phi[20][1000];
506                                 Double_t Track_n  ;
507                                 Double_t Track_pt [1000];
508                                 Double_t Track_eta[1000];
509                                 Double_t Track_phi[1000];
510
511                                 Track_n=0;
512                                 for(int i=0;i<20;i++){
513                                                                 Jet_n[i]=0;
514                                                                 subJet_n[i]=0;
515                                                                 for(int j=0;j<1000;j++){
516                                                                                                 Jet_pt[i][j]=0.;
517                                                                                                 Jet_phi[i][j]=999.;
518                                                                                                 Jet_eta[i][j]=999.;
519                                                                                                 subJet_pt[i][j]=0.;
520                                                                                                 subJet_phi[i][j]=999.;
521                                                                                                 subJet_eta[i][j]=999.;
522                                                                                                 Track_pt [j]=0.;
523                                                                                                 Track_phi[j]=999.;
524                                                                                                 Track_eta[j]=999.;
525                                                                 }
526                                 }
527
528                                 //////-----------------------------------------------------------------------------------
529
530
531                                 int nLJetAOD=999; double ptLJetAOD=0;double phiLJetAOD=999.;double etaLJetAOD=999.;int nsLJetAOD=900;double ptsLJetAOD=0;double phisLJetAOD=900.;double etasLJetAOD=900.;
532                                 int nLJetMC2=999; double ptLJetMC2=0;double phiLJetMC2=999.;double etaLJetMC2=999.;int nsLJetMC2=900;double ptsLJetMC2=0;double phisLJetMC2=900.;double etasLJetMC2=900.;
533                                 int nLJetMC =999; double ptLJetMC =0;double phiLJetMC =999.;double etaLJetMC =999.;int nsLJetMC =900;double ptsLJetMC =0;double phisLJetMC =900.;double etasLJetMC =900.;
534                                 bool findLJetAOD=false;bool findsLJetAOD=false;bool findsLJetAOD_temp=false;
535                                 bool findLJetMC2=false;bool findsLJetMC2=false;bool findsLJetMC2_temp=false;
536                                 bool findLJetMC =false;bool findsLJetMC =false;bool findsLJetMC_temp =false;
537                                 bool findsLJet=false;
538                                 int nLJet = 999;
539
540
541                                 TString cAdd = "";
542                                 cAdd += Form("%02d_",(int)((Radius+0.01)*10.));
543                                 cAdd += Form("B%d",(int)BackM);
544                                 cAdd += Form("_Filter%05d",Filtermask);
545                                 cAdd += Form("_Cut%05d",(int)(1000.*TrackPtcut));
546                                 cAdd += Form("_Skip%02d",SkipCone);
547                                 TString Branchname_gen,Branchname_gen2,Branchname_rec;
548                                 //Branchname_gen  = Form("clustersMCKINE_%s%s",JFAlg.Data(),cAdd.Data());
549                                 //Branchname_gen2 = Form("clustersMCKINE2_%s%s",JFAlg.Data(),cAdd.Data());
550                                 Branchname_gen  = Form("clustersAODMC_%s%s",JFAlg.Data(),cAdd.Data());
551                                 Branchname_gen2 = Form("clustersAODMC2_%s%s",JFAlg.Data(),cAdd.Data());
552                                 Branchname_rec  = Form("clustersAOD_%s%s",JFAlg.Data(),cAdd.Data());
553
554
555                                 for(int algorithm=0;algorithm<3;algorithm++){
556                                                                 //for LHC11a1  LHC11a2
557                                                                 if(algorithm==0)fJetBranch   = Branchname_rec.Data();
558                                                                 if(algorithm==1)fJetBranch   = Branchname_gen2.Data();
559                                                                 if(algorithm==2)fJetBranch   = Branchname_gen.Data();
560
561                                                                 if((!IsMC&&(algorithm==1||algorithm==2)))continue;
562
563                                                                 TClonesArray* jets = dynamic_cast <TClonesArray*> (fAODIn->FindListObject(fJetBranch.Data()));
564                                                                 if(!jets){
565                                                                                                 printf(" Tere are no Branch named %s \n",fJetBranch.Data());
566                                                                                                 continue;
567                                                                 }
568                                                                 Int_t nj = jets->GetEntriesFast();
569                                                                 if (fDebug) printf("There are %5d jets in the event \n", nj);
570                                                                 AliAODJet* jetsAOD;
571                                                                 Jet_n[algorithm] = nj;
572                                                                 int nLjet_in05_pthdiv[20][20];
573                                                                 int nsLjet_in05_pthdiv[20][20];
574                                                                 for(int i=0;i<20;i++){
575                                                                                                 for(int j=0;j<20;j++){
576                                                                                                                                 nLjet_in05_pthdiv [i][j]=0;
577                                                                                                                                 nsLjet_in05_pthdiv[i][j]=0;
578                                                                                                 }
579                                                                 }
580                                                                 //Find Leading Jet -------------------------------------------------------
581                                                                 for(int njet =0;njet<nj;njet++){
582                                                                                                 jetsAOD = (AliAODJet*) (jets->At(njet));
583                                                                                                 Jet_pt   [algorithm][njet] = jetsAOD->Pt();
584                                                                                                 Jet_phi  [algorithm][njet] = jetsAOD->Phi();  
585                                                                                                 Jet_eta  [algorithm][njet] = jetsAOD->Eta();
586                                                                                                 //TRefArray *reftracks = jetsAOD->GetRefTracks();
587                                                                                                 double eta_cut_Jet=0.5;
588                                                                                                 if((TMath::Abs(Jet_eta[algorithm][njet])<eta_cut_Jet)&&(Jet_pt[algorithm][njet]>10.)){
589                                                                                                                                 if(algorithm==0){
590                                                                                                                                                                  fH1JetAKT04_pt->Fill(Jet_pt[algorithm][njet]);  
591                                                                                                                                                                 if(Jet_pt[algorithm][njet]>ptLJetAOD){
592                                                                                                                                                                                                 findLJetAOD=true;
593                                                                                                                                                                                                 nLJetAOD=njet;ptLJetAOD=Jet_pt[algorithm][njet];phiLJetAOD=Jet_phi[algorithm][njet];etaLJetAOD=Jet_eta[algorithm][njet];
594                                                                                                                                                                 }
595                                                                                                                                 }
596                                                                                                                                 if(algorithm==1){
597                                                                                                                                                                  fH1JetMCAKT04_pt->Fill(Jet_pt[algorithm][njet]);  
598                                                                                                                                                                 if(Jet_pt[algorithm][njet]>ptLJetMC2){
599                                                                                                                                                                                                 findLJetMC2=true;
600                                                                                                                                                                                                 nLJetMC2=njet;ptLJetMC2=Jet_pt[algorithm][njet];phiLJetMC2=Jet_phi[algorithm][njet];etaLJetMC2=Jet_eta[algorithm][njet];
601                                                                                                                                                                 }
602                                                                                                                                 }
603                                                                                                                                 if(algorithm==2){
604                                                                                                                                                                 if(Jet_pt[algorithm][njet]>ptLJetMC){
605                                                                                                                                                                                                 findLJetMC=true;
606                                                                                                                                                                                                 nLJetMC=njet;ptLJetMC=Jet_pt[algorithm][njet];phiLJetMC=Jet_phi[algorithm][njet];etaLJetMC=Jet_eta[algorithm][njet];
607                                                                                                                                                                 }
608                                                                                                                                 }
609                                                                                                 }
610                                                                 }//njet loop
611                                                                 //Leading Jet -----------------------------------------------------------
612
613                                                                 if(algorithm==0){nLJet=nLJetAOD;fH1leadJetAKT04_pt->Fill(Jet_pt[algorithm][nLJet]);}
614                                                                 if(algorithm==1){nLJet=nLJetMC2;fH1leadJetMCAKT04_pt->Fill(Jet_pt[algorithm][nLJet]);}
615                                                                 if(algorithm==2){nLJet=nLJetMC2;}
616
617                                                                 if(nj<2)continue;
618                                                                 //Find Sub leading Jet ==================================================
619                                                                 for(int njet=0;njet<nj;njet++){
620                                                                                                 if(njet==nLJet)continue;
621                                                                                                 jetsAOD = (AliAODJet *)jets->At(njet);
622                                                                                                 subJet_pt [algorithm][njet] = jetsAOD->Pt();
623                                                                                                 subJet_phi[algorithm][njet] = jetsAOD->Phi();
624                                                                                                 subJet_eta[algorithm][njet] = jetsAOD->Eta();
625                                                                                                 double eta_cut_Jet=0.5;
626                                                                                                 if((TMath::Abs(subJet_eta[algorithm][njet])<eta_cut_Jet) && (subJet_pt[algorithm][njet]>10.)){
627                                                                                                                                 if(subJet_pt[algorithm][njet]>ptsLJetAOD&&algorithm==0){
628                                                                                                                                                                 findsLJetAOD_temp=true;
629                                                                                                                                                                 nsLJetAOD=njet;ptsLJetAOD=Jet_pt[algorithm][njet];phisLJetAOD=Jet_phi[algorithm][njet];etasLJetAOD=Jet_eta[algorithm][njet];
630                                                                                                                                 }
631                                                                                                                                 if(subJet_pt[algorithm][njet]>ptsLJetMC2 &&algorithm==1){
632                                                                                                                                                                 findsLJetMC2_temp=true;
633                                                                                                                                                                 nsLJetMC2=njet;ptsLJetMC2=Jet_pt[algorithm][njet];phisLJetMC2=Jet_phi[algorithm][njet];etasLJetMC2=Jet_eta[algorithm][njet];
634                                                                                                                                 }
635                                                                                                                                 if(subJet_pt[algorithm][njet]>ptsLJetMC &&algorithm==2){
636                                                                                                                                                                 findsLJetMC_temp=true;
637                                                                                                                                                                 nsLJetMC =njet;ptsLJetMC =Jet_pt[algorithm][njet];phisLJetMC =Jet_phi[algorithm][njet];etasLJetMC =Jet_eta[algorithm][njet];
638                                                                                                                                 }
639                                                                                                 }
640                                                                 }
641                                                                 ////Sub leading Jet ======================================================
642
643                                                                 double Leading_pt=0.;double Leading_phi=999.;double Leading_eta=999.;double sLeading_pt=0.;double sLeading_phi=999.;double sLeading_eta=999.;
644                                                                 if(algorithm==0){Leading_pt=ptLJetAOD;Leading_phi=phiLJetAOD;Leading_eta=etaLJetAOD;sLeading_pt=ptsLJetAOD;sLeading_phi=phisLJetAOD;sLeading_eta=etasLJetAOD;}
645                                                                 if(algorithm==1){Leading_pt=ptLJetMC2;Leading_phi=phiLJetMC2;Leading_eta=etaLJetMC2;sLeading_pt=ptsLJetMC2;sLeading_phi=phisLJetMC2;sLeading_eta=etasLJetMC2;}
646                                                                 if(algorithm==2){Leading_pt=ptLJetMC ;Leading_phi=phiLJetMC ;Leading_eta=etaLJetMC ;sLeading_pt=ptsLJetMC ;sLeading_phi=phisLJetMC ;sLeading_eta=etasLJetMC ;}
647
648                                                                 ////Di-Jet event trigger +++++++++++++++++++++++++++++++++++++++++++++++
649                                                                 double DPhi = DeltaPhi(Leading_phi,sLeading_phi);
650                                                                 double DEta = Leading_eta-sLeading_eta;
651                                                                 if(algorithm==0){
652                                                                                                 fH2JetsJetAKT04_dphi->Fill(Leading_pt,DPhi);
653                                                                                                 fH2JetsJetAKT04_deta->Fill(Leading_pt,DEta);
654                                                                 }
655                                                                 if(algorithm==1){
656                                                                                                 fH2JetsJetMCAKT04_dphi->Fill(Leading_pt,DPhi);
657                                                                                                 fH2JetsJetMCAKT04_deta->Fill(Leading_pt,DEta);
658                                                                 }
659                                                                 if((TMath::Cos(DPhi)<-0.5)&&(Leading_pt>10.)&&(sLeading_pt>10.)){
660                                                                                                 if(algorithm==0)findsLJetAOD=true;                                                         
661                                                                                                 if(algorithm==1)findsLJetMC=true;                                                         
662
663                                                                                                 double Aj = (Leading_pt-sLeading_pt)/(Leading_pt+sLeading_pt);
664                                                                                                 if(algorithm==0){
665                                                                                                                                 fH1leadJetAKT04_pt_dijet->Fill(Leading_pt);
666                                                                                                                                 fH1subJetAKT04_pt_dijet ->Fill(sLeading_pt);
667                                                                                                                                 fH2JetsJetAKT04_Aj->Fill(Leading_pt,Aj);
668                                                                                                                                 fH2JetsJetAKT04_pt->Fill(Leading_pt,sLeading_pt);
669                                                                                                 }
670                                                                                                 if(algorithm==1){
671                                                                                                                                 fH1leadJetMCAKT04_pt_dijet->Fill(Leading_pt);
672                                                                                                                                 fH1subJetMCAKT04_pt_dijet ->Fill(sLeading_pt);
673                                                                                                                                 fH2JetsJetMCAKT04_Aj->Fill(Leading_pt,Aj);
674                                                                                                                                 fH2JetsJetMCAKT04_pt->Fill(Leading_pt,sLeading_pt);
675                                                                                                 }
676
677                                                                 }
678                                                                 if(algorithm==0)findsLJet=findsLJetAOD;
679                                                                 if(algorithm==1)findsLJet=findsLJetMC2;
680                                                                 if(algorithm==2)findsLJet=findsLJetMC;
681                                                                 ////++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
682
683                                                                 if(algorithm!=0)continue;
684
685                                                                 if((findsLJet)&&(Leading_pt>10.)&&(sLeading_pt>10.)){
686
687                                                                                                 for(int eb=0;eb<5;eb++){//count number of Di-Jet in pt bin
688                                                                                                                                 if(TMath::Abs(Leading_pt -20.*(eb+1))<10.){
689                                                                                                                                                                 if(algorithm==0)fH1AKT04_ndiJ_ediv[eb]->Fill(1);
690                                                                                                                                 }
691                                                                                                 }
692
693                                                                                                 TClonesArray* tracks = dynamic_cast <TClonesArray*> (fAODIn->GetTracks());
694                                                                                                 if(!tracks){
695                                                                                                                                 if (fDebug > 1)  Printf("%s:%d could not get AODtracks", (char*)__FILE__,__LINE__);
696                                                                                                                                 continue;
697                                                                                                 }
698                                                                                                 Int_t nt = fAODIn->GetNumberOfTracks();
699
700                                                                                                 AliAODTrack* trackAOD;
701                                                                                                 Track_n = nt;
702                                                                                                 double eta_cut_Jet=0.5;
703                                                                                                 if(TMath::Abs(Leading_eta)<eta_cut_Jet){
704                                                                                                                                 for(int eb=0;eb<5;eb++){
705                                                                                                                                                                 if(TMath::Abs(Leading_pt -20.*(eb+1))<10.){
706                                                                                                                                                                                                 for(int ntrack =0;ntrack<nt;ntrack++){
707                                                                                                                                                                                                                                 trackAOD = (AliAODTrack*) (fAODIn->GetTrack(ntrack));
708                                                                                                                                                                                                                                 Bool_t bgoodT=false;
709                                                                                                                                                                                                                                 if(Filtermask!=768){
710                                                                                                                                                                                                                                                                 if(trackAOD->TestFilterMask(Filtermask))bgoodT=true;
711                                                                                                                                                                                                                                 }
712                                                                                                                                                                                                                                 else{
713                                                                                                                                                                                                                                                                 if(trackAOD->IsHybridGlobalConstrainedGlobal())bgoodT=true; //for hybrid Track cuts
714                                                                                                                                                                                                                                 }
715                                                                                                                                                                                                                                 if(!bgoodT)continue;
716                                                                                                                                                                                                                                 Track_pt   [ntrack]      = trackAOD->Pt();
717                                                                                                                                                                                                                                 Track_phi  [ntrack]      = trackAOD->Phi();
718                                                                                                                                                                                                                                 Track_eta  [ntrack]      = trackAOD->Eta();
719                                                                                                                                                                                                                                 double DelPhi = DeltaPhi(Leading_phi,Track_phi[ntrack]);
720                                                                                                                                                                                                                                 if(TMath::Abs(Track_eta[ntrack])<0.9){
721                                                                                                                                                                                                                                                                 for(int teb=0;teb<5;teb++){
722                                                                                                                                                                                                                                                                                                 if(teb==0){if(!(Track_pt[ntrack]>0.15))continue;}
723                                                                                                                                                                                                                                                                                                 if(teb==1){if(!((Track_pt[ntrack]<1.5)&&(Track_pt[ntrack]>0.15)))continue;}
724                                                                                                                                                                                                                                                                                                 if(teb==2){if(!((Track_pt[ntrack]<3.0)&&(Track_pt[ntrack]>1.5)))continue;}
725                                                                                                                                                                                                                                                                                                 if(teb==3){if(!((Track_pt[ntrack]<4.5)&&(Track_pt[ntrack]>3.0)))continue;}
726                                                                                                                                                                                                                                                                                                 if(teb==4){if(!(Track_pt[ntrack]>4.5))continue;}
727                                                                                                                                                                                                                                                                                                 if(algorithm==0){
728                                                                                                                                                                                                                                                                                                                                 fH1JetHadronAKT04_dphi_ediv                [eb][teb]->Fill(DelPhi); 
729                                                                                                                                                                                                                                                                                                                                 fH1JetHadronAKT04_dphi_tptweight_ediv      [eb][teb]->Fill(DelPhi,Track_pt[ntrack]);
730                                                                                                                                                                                                                                                                                                                                 fH1JetHadronAKT04_dphi_tJptweight_ediv     [eb][teb]->Fill(DelPhi,Track_pt[ntrack]/Leading_pt);
731                                                                                                                                                                                                                                                                                                 }
732                                                                                                                                                                                                                                                                 }
733                                                                                                                                                                                                                                 }
734                                                                                                                                                                                                 }//Track Loop
735                                                                                                                                                                 }
736                                                                                                                                 }
737                                                                                                 }//eta cut
738                                                                 }// Di-Jet
739                                 }// algorithm LOOP
740                                 if(IsMC){
741                                                                 for(int eb=0;eb<5;eb++){
742                                                                                                 double DPhi;
743                                                                                                 if(TMath::Abs(ptLJetAOD -20.*(eb+1))<10.){
744                                                                                                                                 DPhi = DeltaPhi(phiLJetMC,phiLJetAOD);
745                                                                                                                                 fH1leadJetMCAKT04_dphiResolution[eb]->Fill(DPhi);
746                                                                                                                                 DPhi = DeltaPhi(phisLJetMC,phisLJetAOD);
747                                                                                                                                 fH1subJetMCAKT04_dphiResolution[eb]->Fill(DPhi);
748                                                                                                 }
749                                                                 }
750                                 }
751
752                                 PostData(1, fHistList);
753                                 return;
754 }      
755
756 //________________________________________________________________________
757 void AliAnalysisTaskJetHadronCorrelation::Terminate(Option_t *){
758                                 // Terminate analysis
759                                 if (fDebug) printf("AnalysisTaskPt: Terminate() \n");
760 }
761
762
763 Int_t  AliAnalysisTaskJetHadronCorrelation::GetListOfJets(TList *list,TClonesArray* jarray,Int_t type){
764
765                                 if(fDebug>2)Printf("%s:%d Selecting jets with cuts %d",(char*)__FILE__,__LINE__,type);
766                                 Int_t iCount = 0;
767
768                                 if(!jarray){
769                                                                 Printf("%s:%d no Jet array",(char*)__FILE__,__LINE__);
770                                                                 return iCount;
771                                 }
772
773
774                                 for(int ij=0;ij<jarray->GetEntries();ij++){
775                                                                 AliAODJet* jet = (AliAODJet*)jarray->At(ij);
776                                                                 if(!jet)continue;
777                                                                 if(type==0){
778                                                                                                 // No cut at all, main purpose here is sorting      
779                                                                                                 list->Add(jet);
780                                                                                                 iCount++;
781                                                                 }
782                                                                 else if(type == 1){
783                                                                                                 // eta cut
784                                                                                                 if(JetSelected(jet)){
785                                                                                                                                 list->Add(jet);
786                                                                                                                                 iCount++;
787                                                                                                 }
788                                                                 }
789                                 }
790
791                                 list->Sort();
792                                 return iCount;
793
794 }
795
796
797
798 Bool_t  AliAnalysisTaskJetHadronCorrelation::JetSelected(AliAODJet *jet){
799                                 Bool_t selected = false;
800
801                                 if(!jet)return selected;
802
803                                 if(fabs(jet->Eta())<fJetRecEtaWindow&&jet->Pt()>fMinJetPt){
804                                                                 selected = kTRUE;
805                                 }
806                                 return selected;
807
808 }
809
810
811 Double_t AliAnalysisTaskJetHadronCorrelation::DeltaPhi(Double_t phi1,Double_t phi2){
812                                 Float_t pi=TMath::Pi();
813                                 Double_t dphi = phi1-phi2;
814                                 if     (dphi<(-1./2*pi))dphi = dphi +2*pi;
815                                 else if(dphi>( 3./2*pi))dphi = dphi -2*pi;
816                                 return dphi;
817 }
818
819 //Float_t AliAnalysisTaskJetHadronCorrelation::GetTotalEvents(const char* currFile){
820 //                              Float_t totalevent;
821 //                              //TString file_es(currFile);
822 //                              //if(file_es.Contains("root_archive.zip#")){
823 //                              //                              Ssiz_t pos1 = file_es.Index("root_archive",12,TString::kExact);
824 //                              //                              Ssiz_t pos = file_es.Index("#",1,pos1,TString::kExact);
825 //                              //                              file_es.Replace(pos+1,20,"");
826 //                              //}
827 //                              //else {
828 //                              //                              // not an archive take the basename....
829 //                              //                              file_es.ReplaceAll(gSystem->BaseName(file_es.Data()),"");
830 //                              //}
831 //
832 //                              //TString cAdd = "";
833 //                              //cAdd += Form("%02d_",(int)((Radius+0.01)*10.));
834 //                              //cAdd += Form("B%d",(int)BackM);
835 //                              //cAdd += Form("_Filter%05d",Filtermask);
836 //                              //cAdd += Form("_Cut%05d",(int)(1000.*TrackPtcut));
837 //                              //cAdd += Form("_Skip%02d",SkipCone);
838 //                              //TString Dirname,Listname;
839 //                              //Dirname  = Form("PWG4_cluster_AOD__%s%s",JFAlg.Data(),cAdd.Data());
840 //                              //Listname = Form("pwg4cluster_AOD__%s%s",JFAlg.Data(),cAdd.Data());
841 //
842 //                              //TFile *feventstat = TFile::Open(Form("%s%s",file_es.Data(),"JetTasksOutput.root"));
843 //                              //gROOT->Cd(Dirname.Data());
844 //                              //TList *templist     = (TList*)gROOT->FindObject(Listname.Data());
845 //                              //TH1F* temphist = (TH1F*)templist->FindObject("fh1Trials");
846 //                              //totalevent = temphist->Integral();
847 //                              ////cout<<temphist->Integral()<<endl;
848 //                              //delete feventstat;
849 //                              return totalevent;
850 //
851 //}
852
853