AliPhysics  4a7363b (4a7363b)
AliAnalysisTaskClusterQA.cxx
Go to the documentation of this file.
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  * *
4  * Authors: Nicolas Schmidt *
5  * Version 1.0 *
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 
17 //---------------------------------------------
18 // Class used to create a tree for QA of clusters
19 //---------------------------------------------
21 
22 #include "AliAnalysisTaskClusterQA.h"
23 #include "TChain.h"
24 #include "TRandom.h"
25 #include "AliAnalysisManager.h"
26 #include "TParticle.h"
27 #include "TVectorF.h"
28 #include "AliPIDResponse.h"
29 #include "TFile.h"
30 #include "AliESDtrackCuts.h"
31 #include "AliAODMCParticle.h"
32 #include "AliAODMCHeader.h"
33 #include "AliAODEvent.h"
34 #include "AliMultSelection.h"
35 
36 class iostream;
37 
38 using namespace std;
39 
40 ClassImp(AliAnalysisTaskClusterQA)
41 
42 //________________________________________________________________________
43 AliAnalysisTaskClusterQA::AliAnalysisTaskClusterQA() : AliAnalysisTaskSE(),
44  fV0Reader(NULL),
45  fV0ReaderName("V0ReaderV1"),
46  fCorrTaskSetting(""),
47  fConversionCuts(NULL),
48  fEventCuts(NULL),
49  fClusterCutsEMC(NULL),
50  fClusterCutsDMC(NULL),
51  fMesonCuts(NULL),
52  fMinNLMCut(1),
53  fMaxNLMCut(1),
54  fInputEvent(NULL),
55  fMCEvent(NULL),
56  fWeightJetJetMC(1),
57  fGeomEMCAL(NULL),
58  fClusterTree(NULL),
59  fIsHeavyIon(kFALSE),
60  ffillTree(-100),
61  ffillHistograms(kFALSE),
62  fOutputList(NULL),
63  fIsMC(kFALSE),
64  fCorrectForNonlinearity(kFALSE),
65  fSaveEventProperties(0),
66  fSaveCells(0),
67  fSaveSurroundingCells(0),
68  fSaveTracks(0),
69  fConeRadius(0),
70  fMinTrackPt(0),
71  fMinClusterEnergy(0),
72  fSaveMCInformation(0),
73  fSaveAdditionalHistos(0),
74  fExtractionPercentages(),
75  fExtractionPercentagePtBins(),
76  fBuffer_EventID(0),
77  fBuffer_ClusterE(0),
78  fBuffer_ClusterPhi(0),
79  fBuffer_ClusterEta(0),
80  fBuffer_ClusterIsEMCAL(kFALSE),
81  fBuffer_MC_Cluster_Flag(0),
82  fBuffer_ClusterNumCells(0),
83  fBuffer_LeadingCell_ID(0),
84  fBuffer_LeadingCell_E(0),
85  fBuffer_LeadingCell_Eta(0),
86  fBuffer_LeadingCell_Phi(0),
87  fBuffer_ClusterM02(0),
88  fBuffer_ClusterM20(0),
89  fBuffer_Event_Vertex_X(0),
90  fBuffer_Event_Vertex_Y(0),
91  fBuffer_Event_Vertex_Z(0),
92  fBuffer_Event_Multiplicity(0),
93  fBuffer_Event_NumActiveCells(0),
94  fBuffer_Cells_ID(0),
95  fBuffer_Cells_E(0),
96  fBuffer_Cells_RelativeEta(0),
97  fBuffer_Cells_RelativePhi(0),
98  fBuffer_Surrounding_NCells(0),
99  fBuffer_Surrounding_Cells_ID(0),
100  fBuffer_Surrounding_Cells_R(0),
101  fBuffer_Surrounding_Cells_E(0),
102  fBuffer_Surrounding_Cells_RelativeEta(0),
103  fBuffer_Surrounding_Cells_RelativePhi(0),
104  fBuffer_Surrounding_NTracks(10),
105  fBuffer_Surrounding_Tracks_R(0),
106  fBuffer_Surrounding_Tracks_Pt(0),
107  fBuffer_Surrounding_Tracks_RelativeEta(0),
108  fBuffer_Surrounding_Tracks_RelativePhi(0),
109  fBuffer_Cluster_MC_Label(-10),
110  fBuffer_Mother_MC_Label(-10),
111  hNCellsInClustersVsCentrality(NULL),
112  hNActiveCellsVsCentrality(NULL),
113  hNActiveCellsAbove50MeVVsCentrality(NULL),
114  hNActiveCellsAbove80MeVVsCentrality(NULL),
115  hNActiveCellsAbove100MeVVsCentrality(NULL),
116  hNActiveCellsAbove150MeVVsCentrality(NULL),
117  hECellsInClustersVsCentrality(NULL),
118  hEActiveCellsVsCentrality(NULL),
119  hEActiveCells50MeVVsCentrality(NULL),
120  hEActiveCells80MeVVsCentrality(NULL),
121  hEActiveCells100MeVVsCentrality(NULL),
122  hEActiveCells150MeVVsCentrality(NULL)
123 {
124  fBuffer_Cells_ID = new Int_t[kMaxActiveCells];
125  fBuffer_Cells_E = new Float_t[kMaxActiveCells];
126  fBuffer_Cells_RelativeEta = new Float_t[kMaxActiveCells];
127  fBuffer_Cells_RelativePhi = new Float_t[kMaxActiveCells];
128  fBuffer_Surrounding_Cells_ID = new Int_t[kMaxActiveCells];
129  fBuffer_Surrounding_Cells_R = new Float_t[kMaxActiveCells];
130  fBuffer_Surrounding_Cells_E = new Float_t[kMaxActiveCells];
131  fBuffer_Surrounding_Cells_RelativeEta = new Float_t[kMaxActiveCells];
132  fBuffer_Surrounding_Cells_RelativePhi = new Float_t[kMaxActiveCells];
133  fBuffer_Surrounding_Tracks_R = new Float_t[kMaxNTracks];
134  fBuffer_Surrounding_Tracks_Pt = new Float_t[kMaxNTracks];
135  fBuffer_Surrounding_Tracks_RelativeEta= new Float_t[kMaxNTracks];
136  fBuffer_Surrounding_Tracks_RelativePhi= new Float_t[kMaxNTracks];
137 }
138 
139 AliAnalysisTaskClusterQA::AliAnalysisTaskClusterQA(const char *name) : AliAnalysisTaskSE(name),
140  fV0Reader(NULL),
141  fV0ReaderName("V0ReaderV1"),
142  fCorrTaskSetting(""),
143  fConversionCuts(NULL),
144  fEventCuts(NULL),
145  fClusterCutsEMC(NULL),
146  fClusterCutsDMC(NULL),
147  fMesonCuts(NULL),
148  fMinNLMCut(1),
149  fMaxNLMCut(1),
150  fInputEvent(NULL),
151  fMCEvent(NULL),
152  fWeightJetJetMC(1),
153  fGeomEMCAL(NULL),
154  fClusterTree(NULL),
155  fIsHeavyIon(kFALSE),
156  ffillTree(-100),
157  ffillHistograms(kFALSE),
158  fOutputList(NULL),
159  fIsMC(kFALSE),
160  fCorrectForNonlinearity(kFALSE),
161  fSaveEventProperties(0),
162  fSaveCells(0),
163  fSaveSurroundingCells(0),
164  fSaveTracks(0),
165  fConeRadius(0),
166  fMinTrackPt(0),
167  fMinClusterEnergy(0),
168  fSaveMCInformation(0),
169  fSaveAdditionalHistos(0),
170  fExtractionPercentages(),
171  fExtractionPercentagePtBins(),
172  fBuffer_EventID(0),
173  fBuffer_ClusterE(0),
174  fBuffer_ClusterPhi(0),
175  fBuffer_ClusterEta(0),
176  fBuffer_ClusterIsEMCAL(kFALSE),
177  fBuffer_MC_Cluster_Flag(0),
178  fBuffer_ClusterNumCells(0),
179  fBuffer_LeadingCell_ID(0),
180  fBuffer_LeadingCell_E(0),
181  fBuffer_LeadingCell_Eta(0),
182  fBuffer_LeadingCell_Phi(0),
183  fBuffer_ClusterM02(0),
184  fBuffer_ClusterM20(0),
185  fBuffer_Event_Vertex_X(0),
186  fBuffer_Event_Vertex_Y(0),
187  fBuffer_Event_Vertex_Z(0),
188  fBuffer_Event_Multiplicity(0),
189  fBuffer_Event_NumActiveCells(0),
190  fBuffer_Cells_ID(0),
191  fBuffer_Cells_E(0),
192  fBuffer_Cells_RelativeEta(0),
193  fBuffer_Cells_RelativePhi(0),
194  fBuffer_Surrounding_NCells(0),
195  fBuffer_Surrounding_Cells_ID(0),
196  fBuffer_Surrounding_Cells_R(0),
197  fBuffer_Surrounding_Cells_E(0),
198  fBuffer_Surrounding_Cells_RelativeEta(0),
199  fBuffer_Surrounding_Cells_RelativePhi(0),
200  fBuffer_Surrounding_NTracks(10),
201  fBuffer_Surrounding_Tracks_R(0),
202  fBuffer_Surrounding_Tracks_Pt(0),
203  fBuffer_Surrounding_Tracks_RelativeEta(0),
204  fBuffer_Surrounding_Tracks_RelativePhi(0),
205  fBuffer_Cluster_MC_Label(-10),
206  fBuffer_Mother_MC_Label(-10),
207  hNCellsInClustersVsCentrality(NULL),
208  hNActiveCellsVsCentrality(NULL),
209  hNActiveCellsAbove50MeVVsCentrality(NULL),
210  hNActiveCellsAbove80MeVVsCentrality(NULL),
211  hNActiveCellsAbove100MeVVsCentrality(NULL),
212  hNActiveCellsAbove150MeVVsCentrality(NULL),
213  hECellsInClustersVsCentrality(NULL),
214  hEActiveCellsVsCentrality(NULL),
215  hEActiveCells50MeVVsCentrality(NULL),
216  hEActiveCells80MeVVsCentrality(NULL),
217  hEActiveCells100MeVVsCentrality(NULL),
218  hEActiveCells150MeVVsCentrality(NULL)
219 {
220  fBuffer_Cells_ID = new Int_t[kMaxActiveCells];
221  fBuffer_Cells_E = new Float_t[kMaxActiveCells];
222  fBuffer_Cells_RelativeEta = new Float_t[kMaxActiveCells];
223  fBuffer_Cells_RelativePhi = new Float_t[kMaxActiveCells];
224  fBuffer_Surrounding_Cells_ID = new Int_t[kMaxActiveCells];
225  fBuffer_Surrounding_Cells_R = new Float_t[kMaxActiveCells];
226  fBuffer_Surrounding_Cells_E = new Float_t[kMaxActiveCells];
227  fBuffer_Surrounding_Cells_RelativeEta = new Float_t[kMaxActiveCells];
228  fBuffer_Surrounding_Cells_RelativePhi = new Float_t[kMaxActiveCells];
229  fBuffer_Surrounding_Tracks_R = new Float_t[kMaxNTracks];
230  fBuffer_Surrounding_Tracks_Pt = new Float_t[kMaxNTracks];
231  fBuffer_Surrounding_Tracks_RelativeEta= new Float_t[kMaxNTracks];
232  fBuffer_Surrounding_Tracks_RelativePhi= new Float_t[kMaxNTracks];
233  // Default constructor
234 
235  DefineInput(0, TChain::Class());
236  DefineOutput(1, TList::Class());
237  DefineOutput(2, TTree::Class());
238 }
239 
240 //________________________________________________________________________
241 AliAnalysisTaskClusterQA::~AliAnalysisTaskClusterQA()
242 {
243  // default deconstructor
244 
245 }
246 //________________________________________________________________________
247 void AliAnalysisTaskClusterQA::UserCreateOutputObjects()
248 {
249  // Create User Output Objects
250 
251  if(fOutputList != NULL){
252  delete fOutputList;
253  fOutputList = NULL;
254  }
255  if(fOutputList == NULL){
256  fOutputList = new TList();
257  fOutputList->SetOwner(kTRUE);
258  }
259 
260  if(ffillHistograms){
261 
262  if(((AliConvEventCuts*)fEventCuts)->GetCutHistograms()){
263  fOutputList->Add(((AliConvEventCuts*)fEventCuts)->GetCutHistograms());
264  }
265 
266  if(((AliCaloPhotonCuts*)fClusterCutsEMC)->GetCutHistograms()){
267  fOutputList->Add(((AliCaloPhotonCuts*)fClusterCutsEMC)->GetCutHistograms());
268  }
269  if(((AliCaloPhotonCuts*)fClusterCutsDMC)->GetCutHistograms()){
270  fOutputList->Add(((AliCaloPhotonCuts*)fClusterCutsDMC)->GetCutHistograms());
271  }
272  if(((AliConversionMesonCuts*)fMesonCuts)->GetCutHistograms()){
273  fOutputList->Add(((AliConversionMesonCuts*)fMesonCuts)->GetCutHistograms());
274  }
275  if(fSaveAdditionalHistos){
276  hNCellsInClustersVsCentrality = new TH2F("NCellsInClusters_Centrality","NCellsInClusters_Centrality",18000, 0, 18000, 100, 0, 100);
277  fOutputList->Add(hNCellsInClustersVsCentrality);
278  hNActiveCellsVsCentrality = new TH2F("NActiveCells_Centrality","NActiveCells_Centrality",18000, 0, 18000, 100, 0, 100);
279  fOutputList->Add(hNActiveCellsVsCentrality);
280  hNActiveCellsAbove50MeVVsCentrality = new TH2F("NActiveCells50MeV_Centrality","NActiveCells50MeV_Centrality",18000, 0, 18000, 100, 0, 100);
281  fOutputList->Add(hNActiveCellsAbove50MeVVsCentrality);
282  hNActiveCellsAbove80MeVVsCentrality = new TH2F("NActiveCells80MeV_Centrality","NActiveCells80MeV_Centrality",18000, 0, 18000, 100, 0, 100);
283  fOutputList->Add(hNActiveCellsAbove80MeVVsCentrality);
284  hNActiveCellsAbove100MeVVsCentrality = new TH2F("NActiveCells100MeV_Centrality","NActiveCells100MeV_Centrality",18000, 0, 18000, 100, 0, 100);
285  fOutputList->Add(hNActiveCellsAbove100MeVVsCentrality);
286  hNActiveCellsAbove150MeVVsCentrality = new TH2F("NActiveCells150MeV_Centrality","NActiveCells150MeV_Centrality",18000, 0, 18000, 100, 0, 100);
287  fOutputList->Add(hNActiveCellsAbove150MeVVsCentrality);
288 
289  hECellsInClustersVsCentrality = new TH2F("EofCellsInClusters_Centrality","EofCellsInClusters_Centrality",2000, 0, 20, 100, 0, 100);
290  fOutputList->Add(hECellsInClustersVsCentrality);
291  hEActiveCellsVsCentrality = new TH2F("EofActiveCells_Centrality","EofActiveCells_Centrality",2000, 0, 20, 100, 0, 100);
292  fOutputList->Add(hEActiveCellsVsCentrality);
293  hEActiveCells50MeVVsCentrality = new TH2F("EofActiveCells50MeV_Centrality","EofActiveCells50MeV_Centrality",2000, 0, 20, 100, 0, 100);
294  fOutputList->Add(hEActiveCells50MeVVsCentrality);
295  hEActiveCells80MeVVsCentrality = new TH2F("EofActiveCells80MeV_Centrality","EofActiveCells80MeV_Centrality",2000, 0, 20, 100, 0, 100);
296  fOutputList->Add(hEActiveCells80MeVVsCentrality);
297  hEActiveCells100MeVVsCentrality = new TH2F("EofActiveCells100MeV_Centrality","EofActiveCells100MeV_Centrality",2000, 0, 20, 100, 0, 100);
298  fOutputList->Add(hEActiveCells100MeVVsCentrality);
299  hEActiveCells150MeVVsCentrality = new TH2F("EofActiveCells150MeV_Centrality","EofActiveCells150MeV_Centrality",2000, 0, 20, 100, 0, 100);
300  fOutputList->Add(hEActiveCells150MeVVsCentrality);
301  }
302  PostData(1, fOutputList);
303  }
304 
305  fClusterTree = new TTree(Form("ClusterQA_%s_%s_%s",(fEventCuts->GetCutNumber()).Data(),(fClusterCutsEMC->GetCutNumber()).Data(),(fClusterCutsDMC->GetCutNumber()).Data()),Form("ClusterQA_%s_%s_%s",(fEventCuts->GetCutNumber()).Data(),(fClusterCutsEMC->GetCutNumber()).Data(),(fClusterCutsDMC->GetCutNumber()).Data()));
306 
307  fClusterTree->Branch("Cluster_E", &fBuffer_ClusterE, "Cluster_E/F");
308  fClusterTree->Branch("Cluster_Eta", &fBuffer_ClusterEta, "Cluster_Eta/F");
309  fClusterTree->Branch("Cluster_Phi", &fBuffer_ClusterPhi, "Cluster_Phi/F");
310  fClusterTree->Branch("Cluster_IsEMCAL", &fBuffer_ClusterIsEMCAL, "Cluster_IsEMCAL/O");
311  fClusterTree->Branch("Cluster_NumCells", &fBuffer_ClusterNumCells, "Cluster_NumCells/I");
312  fClusterTree->Branch("Cluster_LeadingCell_ID", &fBuffer_LeadingCell_ID, "Cluster_LeadingCell_ID/I");
313  fClusterTree->Branch("Cluster_LeadingCell_E", &fBuffer_LeadingCell_E, "Cluster_LeadingCell_E/F");
314  fClusterTree->Branch("Cluster_LeadingCell_Eta", &fBuffer_LeadingCell_Eta, "Cluster_LeadingCell_Eta/F");
315  fClusterTree->Branch("Cluster_LeadingCell_Phi", &fBuffer_LeadingCell_Phi, "Cluster_LeadingCell_Phi/F");
316  fClusterTree->Branch("Cluster_M02", &fBuffer_ClusterM02, "Cluster_M02/F");
317  fClusterTree->Branch("Cluster_M20", &fBuffer_ClusterM20, "Cluster_M20/F");
318  if(fSaveEventProperties)
319  {
320  fClusterTree->Branch("Event_ID", &fBuffer_EventID, "Event_ID/l");
321  fClusterTree->Branch("Event_Vertex_X", &fBuffer_Event_Vertex_X, "Event_Vertex_X/F");
322  fClusterTree->Branch("Event_Vertex_Y", &fBuffer_Event_Vertex_Y, "Event_Vertex_Y/F");
323  fClusterTree->Branch("Event_Vertex_Z", &fBuffer_Event_Vertex_Z, "Event_Vertex_Z/F");
324  fClusterTree->Branch("Event_Multiplicity", &fBuffer_Event_Multiplicity, "Event_Multiplicity/F");
325  fClusterTree->Branch("Event_NumActiveCells", &fBuffer_Event_NumActiveCells, "Event_NumActiveCells/I");
326  }
327 
328  if(fSaveCells)
329  {
330  fClusterTree->Branch("Cluster_Cells_ID", fBuffer_Cells_ID, "Cluster_Cells_ID[Cluster_NumCells]/I");
331  fClusterTree->Branch("Cluster_Cells_E", fBuffer_Cells_E, "Cluster_Cells_E[Cluster_NumCells]/F");
332  fClusterTree->Branch("Cluster_Cells_RelativeEta", fBuffer_Cells_RelativeEta, "Cluster_Cells_RelativeEta[Cluster_NumCells]/F");
333  fClusterTree->Branch("Cluster_Cells_RelativePhi", fBuffer_Cells_RelativePhi, "Cluster_Cells_RelativePhi[Cluster_NumCells]/F");
334  }
335 
336  if(fSaveSurroundingCells)
337  {
338  fClusterTree->Branch("Surrounding_NCells", &fBuffer_Surrounding_NCells, "Surrounding_NCells/I");
339  fClusterTree->Branch("Surrounding_Cells_ID", fBuffer_Surrounding_Cells_ID, "Surrounding_Cells_ID[Surrounding_NCells]/I");
340  fClusterTree->Branch("Surrounding_Cells_R", fBuffer_Surrounding_Cells_R, "Surrounding_Cells_R[Surrounding_NCells]/F");
341  fClusterTree->Branch("Surrounding_Cells_E", fBuffer_Surrounding_Cells_E, "Surrounding_Cells_E[Surrounding_NCells]/F");
342  fClusterTree->Branch("Surrounding_Cells_RelativeEta", fBuffer_Surrounding_Cells_RelativeEta, "Surrounding_Cells_RelativeEta[Surrounding_NCells]/F");
343  fClusterTree->Branch("Surrounding_Cells_RelativePhi", fBuffer_Surrounding_Cells_RelativePhi, "Surrounding_Cells_RelativePhi[Surrounding_NCells]/F");
344  }
345  if(fSaveTracks)
346  {
347  fClusterTree->Branch("Surrounding_NTracks", &fBuffer_Surrounding_NTracks, "Surrounding_NTracks/I");
348  fClusterTree->Branch("Surrounding_Tracks_R", fBuffer_Surrounding_Tracks_R, "Surrounding_Tracks_R[Surrounding_NTracks]/F");
349  fClusterTree->Branch("Surrounding_Tracks_Pt", fBuffer_Surrounding_Tracks_Pt, "Surrounding_Tracks_Pt[Surrounding_NTracks]/F");
350  fClusterTree->Branch("Surrounding_Tracks_RelativeEta", fBuffer_Surrounding_Tracks_RelativeEta, "Surrounding_Tracks_RelativeEta[Surrounding_NTracks]/F");
351  fClusterTree->Branch("Surrounding_Tracks_RelativePhi", fBuffer_Surrounding_Tracks_RelativePhi, "Surrounding_Tracks_RelativePhi[Surrounding_NTracks]/F");
352  }
353 
354  if(fSaveMCInformation)
355  {
356  fClusterTree->Branch("Cluster_MC_Label", &fBuffer_Cluster_MC_Label, "Cluster_MC_Label/I");
357  fClusterTree->Branch("Mother_MC_Label", &fBuffer_Mother_MC_Label, "Mother_MC_Label/I");
358  }
359 
360  fV0Reader=(AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data());
361  OpenFile(2);
362  PostData(2, fClusterTree);
363 }
364 //_____________________________________________________________________________
365 Bool_t AliAnalysisTaskClusterQA::Notify()
366 {
367  if (fEventCuts->GetPeriodEnum() == AliConvEventCuts::kNoPeriod && ((AliConvEventCuts*)fV0Reader->GetEventCuts())->GetPeriodEnum() != AliConvEventCuts::kNoPeriod){
368  fEventCuts->SetPeriodEnumExplicit(((AliConvEventCuts*)fV0Reader->GetEventCuts())->GetPeriodEnum());
369  } else if (fEventCuts->GetPeriodEnum() == AliConvEventCuts::kNoPeriod ){
370  fEventCuts->SetPeriodEnum(fV0Reader->GetPeriodName());
371  }
372 
373  return kTRUE;
374 }
375 //________________________________________________________________________
376 void AliAnalysisTaskClusterQA::UserExec(Option_t *){
377 
378  Int_t eventQuality = ((AliConvEventCuts*)fV0Reader->GetEventCuts())->GetEventQuality();
379  if(eventQuality != 0){// Event Not Accepted
380  return;
381  }
382  fInputEvent = InputEvent();
383  if(fIsMC) fMCEvent = MCEvent();
384 
385  Int_t eventNotAccepted = fEventCuts->IsEventAcceptedByCut(fV0Reader->GetEventCuts(),fInputEvent,fMCEvent,fIsHeavyIon,kFALSE);
386  if(eventNotAccepted) return; // Check Centrality, PileUp, SDD and V0AND --> Not Accepted => eventQuality = 1
387 
388  fGeomEMCAL = AliEMCALGeometry::GetInstance();
389  if(!fGeomEMCAL){ AliFatal("EMCal geometry not initialized!");}
390 
391  Int_t nclus = 0;
392  TClonesArray * arrClustersProcess = NULL;
393  // fNCurrentClusterBasic = 0;
394  if(!fCorrTaskSetting.CompareTo("")){
395  nclus = fInputEvent->GetNumberOfCaloClusters();
396  } else {
397  arrClustersProcess = dynamic_cast<TClonesArray*>(fInputEvent->FindListObject(Form("%sClustersBranch",fCorrTaskSetting.Data())));
398  if(!arrClustersProcess)
399  AliFatal(Form("%sClustersBranch was not found in AliAnalysisTaskGammaCalo! Check the correction framework settings!",fCorrTaskSetting.Data()));
400  nclus = arrClustersProcess->GetEntries();
401  }
402 
403  if(nclus == 0) return;
404 
405  ((AliCaloPhotonCuts*)fClusterCutsEMC)->FillHistogramsExtendedQA(fInputEvent,fIsMC);
406  // ((AliCaloPhotonCuts*)fClusterCutsDMC)->FillHistogramsExtendedQA(fInputEvent,fIsMC);
407 
408  // match tracks to clusters
409  ((AliCaloPhotonCuts*)fClusterCutsEMC)->MatchTracksToClusters(fInputEvent,fWeightJetJetMC,kTRUE, fMCEvent);
410  // ((AliCaloPhotonCuts*)fClusterCutsDMC)->MatchTracksToClusters(fInputEvent,fWeightJetJetMC,kTRUE, fMCEvent);
411 
412  AliVCaloCells* cells = fInputEvent->GetEMCALCells();
413  Double_t centrality = 0;
414  if(fSaveAdditionalHistos){
415  centrality = GetCentrality(fInputEvent);
416 
417  Int_t numberOfActiveCells50MeV =0;
418  Int_t numberOfActiveCells80MeV =0;
419  Int_t numberOfActiveCells100MeV =0;
420  Int_t numberOfActiveCells150MeV =0;
421  for(Int_t aCell=0;aCell<cells->GetNumberOfCells();aCell++){
422  // Define necessary variables
423  Short_t cellNumber = 0;
424  Double_t cellAmplitude = 0, cellTime = 0, cellEFrac = 0;
425  Int_t cellMCLabel = 0;
426  // Get Cell ID
427  cells->GetCell(aCell,cellNumber,cellAmplitude,cellTime,cellMCLabel,cellEFrac);
428  Double_t cellEnergy = cells->GetCellAmplitude(cellNumber);
429  hEActiveCellsVsCentrality->Fill(cellEnergy,centrality);
430  if(cellEnergy >0.050 ){
431  numberOfActiveCells50MeV++;
432  hEActiveCells50MeVVsCentrality->Fill(cellEnergy,centrality);
433  }
434  if(cellEnergy >0.080 ){
435  numberOfActiveCells80MeV++;
436  hEActiveCells80MeVVsCentrality->Fill(cellEnergy,centrality);
437  }
438  if(cellEnergy >0.100 ){
439  numberOfActiveCells100MeV++;
440  hEActiveCells100MeVVsCentrality->Fill(cellEnergy,centrality);
441  }
442  if(cellEnergy >0.150 ){
443  numberOfActiveCells150MeV++;
444  hEActiveCells150MeVVsCentrality->Fill(cellEnergy,centrality);
445  }
446  }
447  hNActiveCellsVsCentrality->Fill(cells->GetNumberOfCells(),centrality);
448  hNActiveCellsAbove50MeVVsCentrality->Fill(numberOfActiveCells50MeV,centrality);
449  hNActiveCellsAbove80MeVVsCentrality->Fill(numberOfActiveCells80MeV,centrality);
450  hNActiveCellsAbove100MeVVsCentrality->Fill(numberOfActiveCells100MeV,centrality);
451  hNActiveCellsAbove150MeVVsCentrality->Fill(numberOfActiveCells150MeV,centrality);
452  }
453 
454  // vertex
455  Double_t vertex[3] = {0};
456  InputEvent()->GetPrimaryVertex()->GetXYZ(vertex);
457 
458  map<Long_t,Int_t> mapIsClusterAccepted;
459  map<Long_t,Int_t> mapIsClusterAcceptedWithoutTrackMatch;
460  Int_t nCellInCluster = 0;
461  // Loop over EMCal clusters
462  for(Long_t i = 0; i < nclus; i++){
463  Double_t tempClusterWeight = fWeightJetJetMC;
464  AliVCluster* clus = NULL;
465  if(fInputEvent->IsA()==AliESDEvent::Class()){
466  if(arrClustersProcess)
467  clus = new AliESDCaloCluster(*(AliESDCaloCluster*)arrClustersProcess->At(i));
468  else
469  clus = new AliESDCaloCluster(*(AliESDCaloCluster*)fInputEvent->GetCaloCluster(i));
470  } else if(fInputEvent->IsA()==AliAODEvent::Class()){
471  if(arrClustersProcess)
472  clus = new AliAODCaloCluster(*(AliAODCaloCluster*)arrClustersProcess->At(i));
473  else
474  clus = new AliAODCaloCluster(*(AliAODCaloCluster*)fInputEvent->GetCaloCluster(i));
475  }
476  if(!clus) continue;
477 
478 
479  if ( !clus->IsEMCAL()){ // for PHOS: cluster->GetType() == AliVCluster::kPHOSNeutral
480  delete clus;
481  continue;
482  }
483  if(fSaveAdditionalHistos){
484  nCellInCluster+=clus->GetNCells();
485  for(Int_t iCell=0;iCell<clus->GetNCells();iCell++){
486  hECellsInClustersVsCentrality->Fill(cells->GetCellAmplitude(clus->GetCellAbsId(iCell)),centrality);
487  }
488  }
489 
490  if ( clus->E() < fMinClusterEnergy){
491  delete clus;
492  continue;
493  }
494 
495  if(!((AliCaloPhotonCuts*)fClusterCutsEMC)->ClusterIsSelected(clus,fInputEvent,fMCEvent,fIsMC, tempClusterWeight,i) && !((AliCaloPhotonCuts*)fClusterCutsDMC)->ClusterIsSelected(clus,fInputEvent,fMCEvent,fIsMC, tempClusterWeight,i)){
496  delete clus;
497  continue;
498  }
499  ResetBuffer();
500  ProcessQATreeCluster(fInputEvent,clus,i);
501  }
502  if(fSaveAdditionalHistos)
503  hNCellsInClustersVsCentrality->Fill(nCellInCluster,centrality);
504 
505  PostData(1, fOutputList);
506 }
507 
509 void AliAnalysisTaskClusterQA::ProcessQATreeCluster(AliVEvent *event, AliVCluster* cluster, Long_t indexCluster){
510  Float_t clusPos[3] = { 0,0,0 };
511  cluster->GetPosition(clusPos);
512  TVector3 clusterVector(clusPos[0],clusPos[1],clusPos[2]);
513  Double_t etaCluster = clusterVector.Eta();
514  Double_t phiCluster = clusterVector.Phi();
515  if (phiCluster < 0) phiCluster += 2*TMath::Pi();
516 
517  if(fSaveEventProperties){
518  // Vertex position x, y, z
519  fBuffer_Event_Vertex_X = fInputEvent->GetPrimaryVertex()->GetX();
520  fBuffer_Event_Vertex_Y = fInputEvent->GetPrimaryVertex()->GetY();
521  fBuffer_Event_Vertex_Z = fInputEvent->GetPrimaryVertex()->GetZ();
522 
523  // Unique event ID of the current cluster
524  // AliESDEvent* esdEvent = dynamic_cast<AliESDEvent*>(event);
525  fBuffer_EventID = 0;
526  // if(esdEvent->GetHeader()){
527  // fBuffer_EventID = GetUniqueEventID(esdEvent->GetHeader());
528  // printf("event id: %lld\n",fBuffer_EventID);
529  // }
530 
531  // V0-based multiplicity of the event
532  fBuffer_Event_Multiplicity = fInputEvent->GetVZEROData()->GetMTotV0A()+fInputEvent->GetVZEROData()->GetMTotV0C();
533  }
534 
535  // Properties of the current cluster
536  fBuffer_ClusterE = cluster->E();
537 
538  if(etaCluster < 0.67 && etaCluster > -0.67){
539  if(phiCluster < 3.2 && phiCluster > 1.4){
540  fBuffer_ClusterIsEMCAL = kTRUE;
541  }
542  }
543  if(etaCluster < 0.67 && etaCluster > -0.67){
544  if(etaCluster > 0.23 || etaCluster < -0.23){
545  if(phiCluster < 5.6 && phiCluster > 4.6){
546  fBuffer_ClusterIsEMCAL = kFALSE;
547  }
548  }
549  }
550 
551  fBuffer_ClusterPhi = phiCluster;
552  fBuffer_ClusterEta = etaCluster;
553  fBuffer_ClusterM02 = cluster->GetM02();
554  fBuffer_ClusterM20 = cluster->GetM20();
555 
556  // Get all cells from the event
557  AliVCaloCells* cells = NULL;
558  if(cluster->IsEMCAL())
559  cells = event->GetEMCALCells();
560  else
561  return;
562 
563  // Determine all active cells in the event
564  if(fSaveEventProperties)
565  fBuffer_Event_NumActiveCells = cells->GetNumberOfCells();
566 
567  // Get the number of cells from the current cluster
568  Int_t nCellCluster = cluster->GetNCells();
569  fBuffer_ClusterNumCells = nCellCluster;
570 
571  // Find the leading cell in the cluster and its position
572  fBuffer_LeadingCell_ID = FindLargestCellInCluster(cluster,fInputEvent);
573  fBuffer_LeadingCell_E = cells->GetCellAmplitude(fBuffer_LeadingCell_ID);
574  Float_t leadcelleta = 0.;
575  Float_t leadcellphi = 0.;
576  fGeomEMCAL->EtaPhiFromIndex(fBuffer_LeadingCell_ID, leadcelleta, leadcellphi);
577  if ( leadcellphi < 0 ) leadcellphi+=TMath::TwoPi();
578  fBuffer_LeadingCell_Eta = leadcelleta;
579  fBuffer_LeadingCell_Phi = leadcellphi;
580  // Treat the remaining cells of the cluster and get their relative position compared to the leading cell
581  if(fSaveCells){
582  for(Int_t iCell=0;iCell<nCellCluster;iCell++){
583  if(iCell<100){ // maximum number of cells for a cluster is set to 100
584  fBuffer_Cells_ID[iCell] = cluster->GetCellAbsId(iCell);
585  fBuffer_Cells_E[iCell] = cells->GetCellAmplitude(cluster->GetCellAbsId(iCell));
586  Float_t celleta = 0.;
587  Float_t cellphi = 0.;
588  fGeomEMCAL->EtaPhiFromIndex(fBuffer_Cells_ID[iCell], celleta, cellphi);
589  if ( cellphi < 0 ) cellphi+=TMath::TwoPi();
590  fBuffer_Cells_RelativeEta[iCell] = leadcelleta-celleta;
591  fBuffer_Cells_RelativePhi[iCell] = leadcellphi-cellphi;
592  }
593  }
594  }
595  if(fSaveSurroundingCells){
596  Int_t nActiveCellsSurroundingInR = 0;
597  // fill surrounding cell buffer
598  for(Int_t aCell=0;aCell<cells->GetNumberOfCells();aCell++){
599  // Define necessary variables
600  Short_t cellNumber = 0;
601  Double_t cellAmplitude = 0, cellTime = 0, cellEFrac = 0;
602  Int_t cellMCLabel = 0;
603  Float_t surrcelleta = 0.;
604  Float_t surrcellphi = 0.;
605  // Get Cell ID
606  cells->GetCell(aCell,cellNumber,cellAmplitude,cellTime,cellMCLabel,cellEFrac);
607 
608  // Get eta and phi for the surounding cells
609  fGeomEMCAL->EtaPhiFromIndex(cellNumber, surrcelleta, surrcellphi);
610  if ( surrcellphi < 0 ) surrcellphi+=TMath::TwoPi();
611  Double_t dR2 = pow(leadcelleta-surrcelleta,2) + pow(leadcellphi-surrcellphi,2);
612  // Select those cells that are within fConeRadius of the leading cluster cell
613  if( dR2 < fConeRadius){
614  fBuffer_Surrounding_Cells_E[nActiveCellsSurroundingInR] = cells->GetCellAmplitude(cellNumber);
615  fBuffer_Surrounding_Cells_ID[nActiveCellsSurroundingInR] = cellNumber;
616  fBuffer_Surrounding_Cells_R[nActiveCellsSurroundingInR] = dR2;
617  fBuffer_Surrounding_Cells_RelativeEta[nActiveCellsSurroundingInR] = leadcelleta-surrcelleta;
618  fBuffer_Surrounding_Cells_RelativePhi[nActiveCellsSurroundingInR] = leadcellphi-surrcellphi;
619  nActiveCellsSurroundingInR+=1;
620  }
621  }
622  fBuffer_Surrounding_NCells = nActiveCellsSurroundingInR;
623  }
624 
625  // write PDG code of mother of particle that mainly contributed to cluster to tree
626  if(fIsMC> 0){
627  if (cluster->GetNLabels()>0){
628  if((cluster->GetLabelAt(0)!=-1)){
629  if(fInputEvent->IsA()==AliESDEvent::Class()){
630  fBuffer_Mother_MC_Label = (fMCEvent->Particle(cluster->GetLabelAt(0)))->GetPdgCode(); // mother of leading contribution
631  } else if(fInputEvent->IsA()==AliAODEvent::Class()){
632  TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(fInputEvent->FindListObject(AliAODMCParticle::StdBranchName()));
633  if (AODMCTrackArray == NULL) return;
634 
635  AliAODMCParticle* particle = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(cluster->GetLabelAt(0)));
636  fBuffer_Mother_MC_Label = particle->GetPdgCode();
637  }
638 
639  } else{
640  // mother is initial collision
641  fBuffer_Mother_MC_Label = 0;
642  }
643  } else{
644  // no mothers found (should not happen)
645  fBuffer_Mother_MC_Label = 0;
646  }
647  }
648  if(fIsMC) fBuffer_Cluster_MC_Label = MakePhotonCandidates(cluster, cells,indexCluster);
649  if(fSaveTracks) ProcessTracksAndMatching(cluster,indexCluster);
650  // Add everything to the tree
651  if (fClusterTree) fClusterTree->Fill();
652 }
653 
654 //________________________________________________________________________
655 void AliAnalysisTaskClusterQA::Terminate(Option_t *){
656 
657 }
658 
659 //________________________________________________________________________
660 Int_t AliAnalysisTaskClusterQA::FindLargestCellInCluster(AliVCluster* cluster, AliVEvent* event){
661 
662  const Int_t nCells = cluster->GetNCells();
663  AliVCaloCells* cells = NULL;
664 
665  if(cluster->IsEMCAL())
666  cells = event->GetEMCALCells();
667  else if(cluster->IsPHOS())
668  cells = event->GetPHOSCells();
669 
670  Float_t eMax = 0.;
671  Int_t idMax = -1;
672 
673  if (nCells < 1) return idMax;
674  for (Int_t iCell = 0;iCell < nCells;iCell++){
675  Int_t cellAbsID = cluster->GetCellsAbsId()[iCell];
676  if (cells->GetCellAmplitude(cellAbsID)> eMax){
677  eMax = cells->GetCellAmplitude(cellAbsID);
678  idMax = cellAbsID;
679  }
680  }
681  return idMax;
682 }
683 
684 //________________________________________________________________________
685 void AliAnalysisTaskClusterQA::GetRowAndColumnFromAbsCellID(Int_t cellIndex, Int_t& row, Int_t& column){
686  Int_t nSupMod=0, nModule=0, nIphi=0, nIeta=0;
687  row=0;
688  column=0;
689  // Get SM number and relative row/column for SM
690  fGeomEMCAL->GetCellIndex(cellIndex, nSupMod,nModule,nIphi,nIeta);
691  fGeomEMCAL->GetCellPhiEtaIndexInSModule(nSupMod,nModule,nIphi,nIeta, row,column);
692  row += nSupMod/2 * 24;
693  column += nSupMod%2 * 48;
694 }
695 
696 //________________________________________________________________________
697 Int_t AliAnalysisTaskClusterQA::MakePhotonCandidates(AliVCluster* clus, AliVCaloCells* cells, Long_t indexCluster){
698 
699  Double_t vertex[3] = {0};
700  InputEvent()->GetPrimaryVertex()->GetXYZ(vertex);
701 
702  TLorentzVector clusterVector;
703  clus->GetMomentum(clusterVector,vertex);
704 
705  TLorentzVector* tmpvec = new TLorentzVector();
706  tmpvec->SetPxPyPzE(clusterVector.Px(),clusterVector.Py(),clusterVector.Pz(),clusterVector.E());
707 
708  // convert to AODConversionPhoton
709  AliAODConversionPhoton *PhotonCandidate=new AliAODConversionPhoton(tmpvec);
710  if(!PhotonCandidate){
711  delete clus;
712  delete tmpvec;
713  if (PhotonCandidate) delete PhotonCandidate;
714  return -9;
715  }
716  // Flag Photon as CaloPhoton
717  PhotonCandidate->SetIsCaloPhoton();
718  PhotonCandidate->SetCaloClusterRef(indexCluster);
719 
720  // get MC label
721  if(fIsMC> 0){
722  Int_t* mclabelsCluster = clus->GetLabels();
723  Int_t nValidClusters = 0;
724  if (clus->GetNLabels()>0){
725  for (Int_t k =0; k< (Int_t)clus->GetNLabels(); k++){
726  if (mclabelsCluster[k]>0){
727  if (nValidClusters< 50)PhotonCandidate->SetCaloPhotonMCLabel(nValidClusters,mclabelsCluster[k]);
728  nValidClusters++;
729  }
730  }
731  }
732  PhotonCandidate->SetNCaloPhotonMCLabels(nValidClusters);
733  }
734 
735  AliAODCaloCluster* clusSub1 = new AliAODCaloCluster();
736  AliAODCaloCluster* clusSub2 = new AliAODCaloCluster();
737 
738 
739  // split clusters according to their shares in the cluster (NLM == 1) needs to be treated differently
740  if (fMinNLMCut == 1 && fMaxNLMCut == 1 ){
741  Int_t absCellIdFirst = ((AliCaloPhotonCuts*)fClusterCutsEMC)->FindLargestCellInCluster(clus, fInputEvent);
742  Int_t absCellIdSecond = ((AliCaloPhotonCuts*)fClusterCutsEMC)->FindSecondLargestCellInCluster(clus, fInputEvent);
743 
744  ((AliCaloPhotonCuts*)fClusterCutsEMC)->SplitEnergy(absCellIdFirst, absCellIdSecond, clus, fInputEvent, fIsMC, clusSub1, clusSub2);
745  } else if (fMinNLMCut > 1 ){
746  const Int_t nc = clus->GetNCells();
747  Int_t absCellIdList[nc];
748 
749  ((AliCaloPhotonCuts*)fClusterCutsEMC)->SplitEnergy(absCellIdList[0], absCellIdList[1], clus, fInputEvent, fIsMC, clusSub1, clusSub2);
750  }
751 
752  // TLorentzvector with sub cluster 1
753  TLorentzVector clusterVector1;
754  clusSub1->GetMomentum(clusterVector1,vertex);
755  TLorentzVector* tmpvec1 = new TLorentzVector();
756  tmpvec1->SetPxPyPzE(clusterVector1.Px(),clusterVector1.Py(),clusterVector1.Pz(),clusterVector1.E());
757  // convert to AODConversionPhoton
758  AliAODConversionPhoton *PhotonCandidate1=new AliAODConversionPhoton(tmpvec1);
759  if(!PhotonCandidate1){
760  delete clusSub1;
761  delete tmpvec1;
762  return -9;
763  }
764  // Flag Photon as CaloPhoton
765  PhotonCandidate1->SetIsCaloPhoton();
766  // TLorentzvector with sub cluster 2
767  TLorentzVector clusterVector2;
768  clusSub2->GetMomentum(clusterVector2,vertex);
769  TLorentzVector* tmpvec2 = new TLorentzVector();
770  tmpvec2->SetPxPyPzE(clusterVector2.Px(),clusterVector2.Py(),clusterVector2.Pz(),clusterVector2.E());
771  // convert to AODConversionPhoton
772  AliAODConversionPhoton *PhotonCandidate2=new AliAODConversionPhoton(tmpvec2);
773  if(!PhotonCandidate2){
774  delete clusSub2;
775  delete tmpvec2;
776  return -9;
777  }
778  // Flag Photon as CaloPhoton
779  PhotonCandidate2->SetIsCaloPhoton();
780  Int_t mclabel = -3;
781  // create pi0 candidate
782  //AliAODConversionMother *pi0cand = new AliAODConversionMother(PhotonCandidate1,PhotonCandidate2);
783 // if((((AliConversionMesonCuts*)fMesonCuts)->MesonIsSelected(pi0cand,kTRUE,fEventCuts->GetEtaShift()))){
784  if(fIsMC> 0 && PhotonCandidate && PhotonCandidate1 && PhotonCandidate2 && fSaveMCInformation){
785  if(fInputEvent->IsA()==AliESDEvent::Class())
786  mclabel = ProcessTrueClusterCandidates(PhotonCandidate, clus->GetM02(), PhotonCandidate1, PhotonCandidate2);
787  if(fInputEvent->IsA()==AliAODEvent::Class())
788  mclabel = ProcessTrueClusterCandidatesAOD(PhotonCandidate, clus->GetM02(), PhotonCandidate1, PhotonCandidate2);
789  return mclabel;
790  } else {
791  return -7;
792  }
793  return -1;
794 }
795 
796 //________________________________________________________________________
797 void AliAnalysisTaskClusterQA::ProcessTracksAndMatching(AliVCluster* clus, Long_t indexCluster){
798 
799 
800  Int_t nTracksInR = 0;
801  Int_t nModules = fGeomEMCAL->GetNumberOfSuperModules();
802 
803  AliESDEvent *esdev = dynamic_cast<AliESDEvent*>(fInputEvent);
804  AliAODEvent *aodev = 0;
805  if (!esdev) {
806  aodev = dynamic_cast<AliAODEvent*>(fInputEvent);
807  if (!aodev) {
808  AliError("Task needs AOD or ESD event, returning");
809  return;
810  }
811  }
812  AliESDtrackCuts *EsdTrackCuts = 0x0;
813  // Using standard function for setting Cuts
814  static int prevRun = -1;
815  // Using standard function for setting Cuts
816  if (esdev){
817  Int_t runNumber = fInputEvent->GetRunNumber();
818  if (prevRun!=runNumber) {
819  delete EsdTrackCuts;
820  EsdTrackCuts = 0;
821  prevRun = runNumber;
822  }
823  // if LHC11a or earlier or if LHC13g or if LHC12a-i -> use 2010 cuts
824  if( (runNumber<=146860) || (runNumber>=197470 && runNumber<=197692) || (runNumber>=172440 && runNumber<=193766) ){
825  EsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2010();
826  // else if run2 data use 2015 PbPb cuts
827  }else if (runNumber>=209122){
828  // hard coded track cuts for the moment, because AliESDtrackCuts::GetStandardITSTPCTrackCuts2015PbPb() gives spams warnings
829  EsdTrackCuts = new AliESDtrackCuts();
830  EsdTrackCuts->AliESDtrackCuts::SetMinNCrossedRowsTPC(70);
831  EsdTrackCuts->AliESDtrackCuts::SetMinRatioCrossedRowsOverFindableClustersTPC(0.8);
832  EsdTrackCuts->AliESDtrackCuts::SetMaxChi2PerClusterTPC(4);
833  EsdTrackCuts->AliESDtrackCuts::SetAcceptKinkDaughters(kFALSE);
834  EsdTrackCuts->AliESDtrackCuts::SetRequireTPCRefit(kTRUE);
835  // ITS
836  EsdTrackCuts->AliESDtrackCuts::SetRequireITSRefit(kTRUE);
837  EsdTrackCuts->AliESDtrackCuts::SetClusterRequirementITS(AliESDtrackCuts::kSPD,
838  AliESDtrackCuts::kAny);
839  EsdTrackCuts->AliESDtrackCuts::SetMaxDCAToVertexXYPtDep("0.0105+0.0350/pt^1.1");
840  EsdTrackCuts->AliESDtrackCuts::SetMaxChi2TPCConstrainedGlobal(36);
841  EsdTrackCuts->AliESDtrackCuts::SetMaxDCAToVertexZ(2);
842  EsdTrackCuts->AliESDtrackCuts::SetDCAToVertex2D(kFALSE);
843  EsdTrackCuts->AliESDtrackCuts::SetRequireSigmaToVertex(kFALSE);
844  EsdTrackCuts->AliESDtrackCuts::SetMaxChi2PerClusterITS(36);
845  // else use 2011 version of track cuts
846  }else{
847  EsdTrackCuts = AliESDtrackCuts::GetStandardITSTPCTrackCuts2011();
848  }
849  EsdTrackCuts->SetMaxDCAToVertexZ(2);
850  EsdTrackCuts->SetEtaRange(-0.8, 0.8);
851  EsdTrackCuts->SetPtRange(0.15);
852  }
853 
854  for (Int_t itr=0;itr<fInputEvent->GetNumberOfTracks();itr++){
855  AliExternalTrackParam *trackParam = 0;
856  AliVTrack *inTrack = 0x0;
857  if(esdev){
858  inTrack = esdev->GetTrack(itr);
859  if(!inTrack) continue;
860  if(inTrack->Pt()<fMinTrackPt) continue;
861  AliESDtrack *esdt = dynamic_cast<AliESDtrack*>(inTrack);
862  // check track cuts
863  if(!EsdTrackCuts->AcceptTrack(esdt)) continue;
864  const AliExternalTrackParam *in = esdt->GetInnerParam();
865  if (!in){AliError("Could not get InnerParam of Track, continue");continue;}
866  trackParam =new AliExternalTrackParam(*in);
867  } else if(aodev){
868  inTrack = dynamic_cast<AliVTrack*>(aodev->GetTrack(itr));
869  if(!inTrack) continue;
870  if(inTrack->Pt()<fMinTrackPt) continue;
871  AliAODTrack *aodt = dynamic_cast<AliAODTrack*>(inTrack);
872  // check track cuts
873  if(!aodt->IsHybridGlobalConstrainedGlobal()) continue;
874 
875  Double_t xyz[3] = {0}, pxpypz[3] = {0}, cv[21] = {0};
876  aodt->GetPxPyPz(pxpypz);
877  aodt->GetXYZ(xyz);
878  aodt->GetCovarianceXYZPxPyPz(cv);
879 
880  trackParam = new AliExternalTrackParam(xyz,pxpypz,cv,aodt->Charge());
881  }
882  AliExternalTrackParam emcParam(*trackParam);
883  Float_t eta, phi, pt;
884 
885  //propagate tracks to emc surfaces
886  if (!AliEMCALRecoUtils::ExtrapolateTrackToEMCalSurface(&emcParam, 440., 0.139, 20., eta, phi, pt)) {
887  delete trackParam;
888  continue;
889  }
890  if( TMath::Abs(eta) > 0.75 ) {
891  delete trackParam;
892  continue;
893  }
894  // Save some time and memory in case of no DCal present
895  if( nModules < 13 && ( phi < 70*TMath::DegToRad() || phi > 190*TMath::DegToRad())){
896  delete trackParam;
897  continue;
898  }
899  // Save some time and memory in case of run2
900  if( nModules > 12 ){
901  if (( phi < 70*TMath::DegToRad() || phi > 190*TMath::DegToRad()) && ( phi < 250*TMath::DegToRad() || phi > 340*TMath::DegToRad())){
902  delete trackParam;
903  continue;
904  }
905  }
906  Float_t dEta=-999, dPhi=-999;
907  Double_t trkPos[3] = {0.,0.,0.};
908  if (!emcParam.GetXYZ(trkPos)){
909  delete trackParam;
910  continue;
911  }
912 
913  AliExternalTrackParam trackParamTmp(emcParam);//Retrieve the starting point every time before the extrapolation
914  if(!AliEMCALRecoUtils::ExtrapolateTrackToCluster(&trackParamTmp, clus, 0.139, 5., dEta, dPhi)) continue;
915 
916  Float_t dR2 = dPhi*dPhi + dEta*dEta;
917  if(dR2 < fConeRadius){
918  fBuffer_Surrounding_Tracks_R[nTracksInR]=dR2;
919  fBuffer_Surrounding_Tracks_Pt[nTracksInR]=inTrack->Pt();
920  fBuffer_Surrounding_Tracks_RelativeEta[nTracksInR]=dEta;
921  fBuffer_Surrounding_Tracks_RelativePhi[nTracksInR]=dPhi;
922  nTracksInR+=1;
923  }
924  }
925  fBuffer_Surrounding_NTracks = nTracksInR;
926  if(nTracksInR==0){
927  fBuffer_Surrounding_Tracks_R[nTracksInR]=-1;
928  fBuffer_Surrounding_Tracks_Pt[nTracksInR]=-1;
929  fBuffer_Surrounding_Tracks_RelativeEta[nTracksInR]=-1;
930  fBuffer_Surrounding_Tracks_RelativePhi[nTracksInR]=-1;
931  }
932 
933  if(EsdTrackCuts){
934  delete EsdTrackCuts;
935  EsdTrackCuts=0x0;
936  }
937 }
938 
939 //________________________________________________________________________
940 Int_t AliAnalysisTaskClusterQA::ProcessTrueClusterCandidates(AliAODConversionPhoton *TrueClusterCandidate, Float_t m02,
941  AliAODConversionPhoton *TrueSubClusterCandidate1,
942  AliAODConversionPhoton *TrueSubClusterCandidate2)
943 {
944 
945  Int_t mcLabelCluster = -5;
946  const AliVVertex* primVtxMC = fMCEvent->GetPrimaryVertex();
947  Double_t mcProdVtxX = primVtxMC->GetX();
948  Double_t mcProdVtxY = primVtxMC->GetY();
949  Double_t mcProdVtxZ = primVtxMC->GetZ();
950 
951  TParticle *Photon = NULL;
952  if (!TrueClusterCandidate->GetIsCaloPhoton()) AliFatal("CaloPhotonFlag has not been set task will abort");
953  if (TrueClusterCandidate->GetCaloPhotonMCLabel(0) < 0){
954  mcLabelCluster = -10;
955  return mcLabelCluster;
956  }
957  if (TrueClusterCandidate->GetNCaloPhotonMCLabels()>0){
958  Photon = fMCEvent->Particle(TrueClusterCandidate->GetCaloPhotonMCLabel(0));
959  } else{
960  mcLabelCluster = -11;
961  return mcLabelCluster;
962  }
963  if(Photon == NULL){
964  mcLabelCluster = -12;
965  return mcLabelCluster;
966  }
967  AliAODConversionMother *mesoncand = new AliAODConversionMother(TrueSubClusterCandidate1,TrueSubClusterCandidate2);
968 // Bool_t mesonIsSelected = (((AliConversionMesonCuts*)fMesonCuts)->MesonIsSelected(mesoncand,kTRUE,fEventCuts->GetEtaShift()));
969 // if (!mesonIsSelected){
970 // delete mesoncand;
971 // mcLabelCluster = -13;
972 // return mcLabelCluster;
973 // }
974 
975  TrueClusterCandidate->SetCaloPhotonMCFlags(fMCEvent, kFALSE);
976 
977  // cluster classification:
978  // 1 - nice merged cluster (2 gamma | contributions from 2 gamma) from pi0/eta
979  // 2 - contribution from only 1 partner (1 gamma, 1 fully coverted gamma) from pi0/eta
980  // 3 - contribution from part of 1 partner (1 electron) from pi0/eta
981 
982 
983  Int_t clusterClass = 0;
984  Bool_t isPrimary = fEventCuts->IsConversionPrimaryESD( fMCEvent, TrueClusterCandidate->GetCaloPhotonMCLabel(0), mcProdVtxX, mcProdVtxY, mcProdVtxZ);
985 
986 
987  Long_t motherLab = -1;
988  if (TrueClusterCandidate->IsMerged() || TrueClusterCandidate->IsMergedPartConv()){
989  clusterClass = 1;
990  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
991  } else if (TrueClusterCandidate->GetNCaloPhotonMotherMCLabels()> 0){
992  if (TrueClusterCandidate->IsLargestComponentElectron() || TrueClusterCandidate->IsLargestComponentPhoton()){
993  if (TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0) > -1 && (fMCEvent->Particle(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0))->GetPdgCode() == 111 || fMCEvent->Particle(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0))->GetPdgCode() == 221) ){
994  if ( TrueClusterCandidate->IsConversion() && !TrueClusterCandidate->IsConversionFullyContained() ){
995  clusterClass = 3;
996  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
997  } else {
998  clusterClass = 2;
999  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1000  }
1001  }
1002  } else if (TrueClusterCandidate->IsSubLeadingEM()){
1003  if (TrueClusterCandidate->GetNCaloPhotonMotherMCLabels()> 1){
1004  if ( TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1) > -1){
1005  if (fMCEvent->Particle(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1))->GetPdgCode() == 111 || fMCEvent->Particle(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1))->GetPdgCode() == 221 ){
1006  clusterClass = 2;
1007  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1);
1008  }
1009  }
1010  }
1011  } else {
1012  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1013  }
1014  }
1015 
1016  // Get Mother particle
1017  TParticle *mother = NULL;
1018  Int_t motherPDG = -1;
1019  if (motherLab > -1)
1020  mother = fMCEvent->Particle(motherLab);
1021  if (mother)
1022  motherPDG = TMath::Abs(mother->GetPdgCode());
1023 
1024  //
1025  if (clusterClass == 1 || clusterClass == 2 || clusterClass == 3 ){
1026  // separate different components
1027  if (clusterClass == 1 && TrueClusterCandidate->IsMerged()){
1028  if (motherPDG == 111){
1029  mcLabelCluster = 10; // NOTE merged pi0
1030  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1031  mcLabelCluster = 11; // NOTE secondary merged pi0
1032  }
1033  if (motherPDG == 221)
1034  mcLabelCluster = 12; // NOTE merged eta
1035  } else if (clusterClass == 1 && TrueClusterCandidate->IsMergedPartConv()){
1036  if (motherPDG == 111){
1037  mcLabelCluster = 13; // NOTE merged pi0 with one converted gamma
1038  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1039  mcLabelCluster = 14; // NOTE merged secondary pi0 with one converted gamma
1040  }
1041  if (motherPDG == 221)
1042  mcLabelCluster = 15; // NOTE merged eta with one converted gamma
1043  } else if (clusterClass == 2){
1044  if (motherPDG == 111){
1045  mcLabelCluster = 20; // NOTE decay photon from pi0
1046  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1047  mcLabelCluster = 21; // NOTE decay photon from secondary pi0
1048  }
1049  if (motherPDG == 221)
1050  mcLabelCluster = 22; // NOTE decay photon from eta
1051  } else if (clusterClass == 3){
1052  if (motherPDG == 111) {
1053  mcLabelCluster = 30; // NOTE electron from decayed pi0
1054  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1055  mcLabelCluster = 31; // NOTE electron from decayed secondary pi0
1056  }
1057  if (motherPDG == 221)
1058  mcLabelCluster = 32; // NOTE electron from decayed eta
1059  }
1060 
1061  // leading particle is a photon or the conversion is fully contained and its not from pi0 || eta
1062  } else if (TrueClusterCandidate->IsLargestComponentPhoton() || TrueClusterCandidate->IsConversionFullyContained()
1063  || TrueClusterCandidate->IsElectronFromFragPhoton()){
1064 
1065  if(TrueClusterCandidate->IsLargestComponentPhoton()){
1066  // cluster is produced by a photon that either has no mother or the mother is neither from a pi^0 nor eta, e.g. inital collision
1067 
1068  if (motherLab == -1) mcLabelCluster = 40; // direct photon from initial collision
1069  else if ((motherLab >0) && (motherLab<9)) mcLabelCluster = 41; // photon from quark
1070  else if (motherLab == 11) mcLabelCluster = 42; // photon from electron
1071  else if (motherLab == 22){ // check for frag photon
1072 
1073  TParticle *dummyMother = fMCEvent->Particle(motherLab);
1074  Bool_t originReached = kFALSE;
1075  Bool_t isFragPhoton = kFALSE;
1076 
1077  while (dummyMother->GetPdgCode() == 22 && !originReached){ // follow photon's history, as long as the mother is a photon
1078  if (dummyMother->GetMother(0) > -1){
1079  dummyMother = fMCEvent->Particle(dummyMother->GetMother(0));
1080  if (TMath::Abs(dummyMother->GetPdgCode()) == 22){ // if mother of mother is again a photon, continue
1081  if (dummyMother->GetMother(0) > -1){
1082  dummyMother = fMCEvent->Particle(dummyMother->GetMother(0));
1083  } else {
1084  originReached = kTRUE;
1085  }
1086  } else {
1087  originReached = kTRUE;
1088  }
1089  isFragPhoton = (TMath::Abs(dummyMother->GetPdgCode()) < 6);// photon stems from quark = fragmentation photon
1090  } else {
1091  originReached = kTRUE;
1092  }
1093  }
1094 
1095  if(isFragPhoton) mcLabelCluster = 43; // Fragmentation photon
1096  else{
1097  mcLabelCluster = 47; // something like cluster <- photon <- photon <- X (not photon)
1098  AliInfo(Form("Mother of photon is photon etc. but origin is not quark. ID: %li", motherLab));
1099  }
1100  }
1101  else{
1102  mcLabelCluster = 44; // other (e.g. from meson decays that are not pi0 or eta0)
1103  AliInfo(Form("Single cluster is mainly produced by a photon and mother is %li", motherLab));
1104  }
1105  } else if(TrueClusterCandidate->IsConversionFullyContained()){
1106  // cluster is from a fully contained conversion
1107  mcLabelCluster = 45;
1108  } else if(TrueClusterCandidate->IsElectronFromFragPhoton()){
1109  mcLabelCluster = 46; // electron from frac
1110  }
1111 
1112  // leading particle is an electron and its not from pi0 || eta and no electron from fragmentation photon conversion
1113  } else if (TrueClusterCandidate->IsLargestComponentElectron() && !TrueClusterCandidate->IsElectronFromFragPhoton()){
1114  mcLabelCluster = 50; // NOTE single electron
1115  } else {
1116  // leading particle from hadron
1117  mcLabelCluster = 60; // NOTE hadron cluster
1118  AliInfo(Form("Single cluster is mainly produced by hadron with id: %li", motherLab));
1119  }
1120 
1121  delete mesoncand;
1122  return mcLabelCluster;
1123 }
1124 
1125 //________________________________________________________________________
1126 Int_t AliAnalysisTaskClusterQA::ProcessTrueClusterCandidatesAOD(AliAODConversionPhoton *TrueClusterCandidate, Float_t m02,
1127  AliAODConversionPhoton *TrueSubClusterCandidate1,
1128  AliAODConversionPhoton *TrueSubClusterCandidate2)
1129 {
1130  Int_t mcLabelCluster = -5;
1131  const AliVVertex* primVtxMC = fMCEvent->GetPrimaryVertex();
1132  Double_t mcProdVtxX = primVtxMC->GetX();
1133  Double_t mcProdVtxY = primVtxMC->GetY();
1134  Double_t mcProdVtxZ = primVtxMC->GetZ();
1135 
1136  AliAODMCParticle *Photon = NULL;
1137  TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(fInputEvent->FindListObject(AliAODMCParticle::StdBranchName()));
1138 
1139  if (AODMCTrackArray){
1140  if (!TrueClusterCandidate->GetIsCaloPhoton()) AliFatal("CaloPhotonFlag has not been set task will abort");
1141  if (TrueClusterCandidate->GetCaloPhotonMCLabel(0) < 0) {
1142  mcLabelCluster = -10;
1143  return mcLabelCluster;
1144  }
1145  if (TrueClusterCandidate->GetNCaloPhotonMCLabels()>0) Photon = (AliAODMCParticle*) AODMCTrackArray->At(TrueClusterCandidate->GetCaloPhotonMCLabel(0));
1146  else{
1147  mcLabelCluster = -11;
1148  return mcLabelCluster;
1149  }
1150  } else {
1151  AliInfo("AODMCTrackArray could not be loaded");
1152  mcLabelCluster = -90;
1153  return mcLabelCluster;
1154  }
1155  if(Photon == NULL){
1156  mcLabelCluster = -12;
1157  return mcLabelCluster;
1158  }
1159  AliAODConversionMother *mesoncand = new AliAODConversionMother(TrueSubClusterCandidate1,TrueSubClusterCandidate2);
1160 
1161  TrueClusterCandidate->SetCaloPhotonMCFlagsAOD(fInputEvent, kFALSE);
1162 
1163  Int_t clusterClass = 0;
1164  Bool_t isPrimary = ((AliConvEventCuts*)fEventCuts)->IsConversionPrimaryAOD( fInputEvent, Photon, mcProdVtxX, mcProdVtxY, mcProdVtxZ);
1165 
1166  Long_t motherLab = -1;
1167  if (TrueClusterCandidate->IsMerged() || TrueClusterCandidate->IsMergedPartConv()){
1168  clusterClass = 1;
1169  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1170  } else if (TrueClusterCandidate->GetNCaloPhotonMotherMCLabels()> 0){
1171  if (TrueClusterCandidate->IsLargestComponentElectron() || TrueClusterCandidate->IsLargestComponentPhoton()){
1172  if (TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0) > -1 && (((AliAODMCParticle*) AODMCTrackArray->At(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0)))->GetPdgCode() == 111 || ((AliAODMCParticle*) AODMCTrackArray->At(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0)))->GetPdgCode() == 221) ){
1173  if ( TrueClusterCandidate->IsConversion() && !TrueClusterCandidate->IsConversionFullyContained() ){
1174  clusterClass = 3;
1175  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1176  } else {
1177  clusterClass = 2;
1178  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1179  }
1180  }
1181  } else if (TrueClusterCandidate->IsSubLeadingEM()){
1182  if (TrueClusterCandidate->GetNCaloPhotonMotherMCLabels()> 1){
1183  if ( TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1) > -1){
1184  if (TMath::Abs(((AliAODMCParticle*) AODMCTrackArray->At(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1)))->GetPdgCode()) == 111 || TMath::Abs(((AliAODMCParticle*) AODMCTrackArray->At(TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1)))->GetPdgCode()) == 221 ){
1185  clusterClass = 2;
1186  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(1);
1187  }
1188  }
1189  }
1190  } else {
1191  motherLab = TrueClusterCandidate->GetCaloPhotonMotherMCLabel(0);
1192  }
1193  }
1194 
1195  // Get Mother particle
1196  AliAODMCParticle *mother = NULL;
1197  Int_t motherPDG = -1;
1198  if (motherLab > -1)
1199  mother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(motherLab));
1200  if (mother)
1201  motherPDG = TMath::Abs(mother->GetPdgCode());
1202 
1203  if (clusterClass == 1 || clusterClass == 2 || clusterClass == 3 ){
1204  // separate different components
1205  if (clusterClass == 1 && TrueClusterCandidate->IsMerged()){
1206  if (motherPDG == 111){
1207  mcLabelCluster = 10; // NOTE merged pi0
1208  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1209  mcLabelCluster = 11; // NOTE secondary merged pi0
1210  }
1211  if (motherPDG == 221)
1212  mcLabelCluster = 12; // NOTE merged eta
1213  } else if (clusterClass == 1 && TrueClusterCandidate->IsMergedPartConv()){
1214  if (motherPDG == 111){
1215  mcLabelCluster = 13; // NOTE merged pi0 with one converted gamma
1216  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1217  mcLabelCluster = 14; // NOTE merged secondary pi0 with one converted gamma
1218  }
1219  if (motherPDG == 221)
1220  mcLabelCluster = 15; // NOTE merged eta with one converted gamma
1221  } else if (clusterClass == 2){
1222  if (motherPDG == 111){
1223  mcLabelCluster = 20; // NOTE decay photon from pi0
1224  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1225  mcLabelCluster = 21; // NOTE decay photon from secondary pi0
1226  }
1227  if (motherPDG == 221)
1228  mcLabelCluster = 22; // NOTE decay photon from eta
1229  } else if (clusterClass == 3){
1230  if (motherPDG == 111) {
1231  mcLabelCluster = 30; // NOTE electron from decayed pi0
1232  if (!isPrimary && m02 >= 0 && m02 <= 4.8 )
1233  mcLabelCluster = 31; // NOTE electron from decayed secondary pi0
1234  }
1235  if (motherPDG == 221)
1236  mcLabelCluster = 32; // NOTE electron from decayed eta
1237  }
1238  // leading particle is a photon or the conversion is fully contained and its not from pi0 || eta
1239  } else if (TrueClusterCandidate->IsLargestComponentPhoton() || TrueClusterCandidate->IsConversionFullyContained()
1240  || TrueClusterCandidate->IsElectronFromFragPhoton()){
1241 
1242  if(TrueClusterCandidate->IsLargestComponentPhoton()){
1243  // cluster is produced by a photon that either has no mother or the mother is neither from a pi^0 nor eta, e.g. inital collision
1244 
1245  if (motherLab == -1) mcLabelCluster = 40; // direct photon from initial collision
1246  else if ((motherLab >0) && (motherLab<9)) mcLabelCluster = 41; // photon from quark
1247  else if (motherLab == 11) mcLabelCluster = 42; // photon from electron
1248  else if (motherLab == 22){ // check for frag photon
1249 
1250  AliAODMCParticle *dummyMother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(motherLab));
1251  Bool_t originReached = kFALSE;
1252  Bool_t isFragPhoton = kFALSE;
1253 
1254  while (dummyMother->GetPdgCode() == 22 && !originReached){ // follow photon's history, as long as the mother is a photon
1255  if (dummyMother->GetMother() > -1){
1256  dummyMother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(dummyMother->GetMother()));
1257 
1258  if (TMath::Abs(dummyMother->GetPdgCode()) == 22){ // if mother of mother is again a photon, continue
1259  if (dummyMother->GetMother() > -1){
1260  dummyMother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(dummyMother->GetMother()));
1261  } else {
1262  originReached = kTRUE;
1263  }
1264  } else {
1265  originReached = kTRUE;
1266  }
1267  isFragPhoton = (TMath::Abs(dummyMother->GetPdgCode()) < 6);// photon stems from quark = fragmentation photon
1268  } else {
1269  originReached = kTRUE;
1270  }
1271  }
1272 
1273  if(isFragPhoton) mcLabelCluster = 43; // Fragmentation photon
1274  else{
1275  mcLabelCluster = 47; // something like cluster <- photon <- photon <- X (not photon)
1276  AliInfo(Form("Mother of photon is photon etc. but origin is not quark. ID: %li", motherLab));
1277  }
1278  }
1279  else{
1280  mcLabelCluster = 44; // other (e.g. from meson decays that are not pi0 or eta0)
1281  AliInfo(Form("Single cluster is mainly produced by a photon and mother is %li", motherLab));
1282  }
1283  } else if(TrueClusterCandidate->IsConversionFullyContained()){
1284  // cluster is from a fully contained conversion
1285  mcLabelCluster = 45;
1286  } else if(TrueClusterCandidate->IsElectronFromFragPhoton()){
1287  mcLabelCluster = 46; // electron from frac
1288  }
1289 
1290  // leading particle is an electron and its not from pi0 || eta and no electron from fragmentation photon conversion
1291  } else if (TrueClusterCandidate->IsLargestComponentElectron() && !TrueClusterCandidate->IsElectronFromFragPhoton()){
1292  mcLabelCluster = 50; // NOTE single electron
1293  } else {
1294  // leading particle from hadron
1295  mcLabelCluster = 60; // NOTE hadron cluster
1296  AliInfo(Form("Single cluster is mainly produced by hadron with id: %li", motherLab));
1297  }
1298 
1299  delete mesoncand;
1300  return mcLabelCluster;
1301 }
1302 //_____________________________________________________________________________
1303 ULong64_t AliAnalysisTaskClusterQA::GetUniqueEventID(AliVHeader* header)
1304 {
1305  // To have a unique id for each event in a run!
1306  // Modified from AliRawReader.h
1307  return ((ULong64_t)header->GetBunchCrossNumber()+
1308  (ULong64_t)header->GetOrbitNumber()*3564+
1309  (ULong64_t)header->GetPeriodNumber()*16777215*3564);
1310 }
1311 
1312 //-------------------------------------------------------------
1313 Float_t AliAnalysisTaskClusterQA::GetCentrality(AliVEvent *event)
1314 { // Get Event Centrality
1315 
1316  AliESDEvent *esdEvent=dynamic_cast<AliESDEvent*>(event);
1317  if(esdEvent){
1318  AliMultSelection *MultSelection = (AliMultSelection*)event->FindListObject("MultSelection");
1319  if(!MultSelection){
1320  AliWarning ("AliMultSelection object not found !");
1321  return -1;
1322  }else{
1323  return MultSelection->GetMultiplicityPercentile("V0M");// default
1324  }
1325  }
1326 
1327 
1328  return -1;
1329 }
1330 void AliAnalysisTaskClusterQA::ResetBuffer(){
1331  fBuffer_EventID = 0;
1332  fBuffer_ClusterE = 0;
1333  fBuffer_ClusterPhi = 0;
1334  fBuffer_ClusterEta = 0;
1335  fBuffer_ClusterIsEMCAL = kFALSE;
1336  fBuffer_MC_Cluster_Flag = 0;
1337  fBuffer_ClusterNumCells = 0;
1338  fBuffer_LeadingCell_ID = 0;
1339  fBuffer_LeadingCell_E = 0;
1340  fBuffer_LeadingCell_Eta = 0;
1341  fBuffer_LeadingCell_Phi = 0;
1342  fBuffer_ClusterM02 = 0;
1343  fBuffer_ClusterM20 = 0;
1344  fBuffer_Event_Vertex_X = 0;
1345  fBuffer_Event_Vertex_Y = 0;
1346  fBuffer_Event_Vertex_Z = 0;
1347  fBuffer_Event_Multiplicity = 0;
1348  fBuffer_Event_NumActiveCells = 0;
1349  fBuffer_Cluster_MC_Label = -10;
1350  fBuffer_Mother_MC_Label = -10;
1351 
1352  for(Int_t cell = 0; cell < kMaxActiveCells; cell++){
1353  fBuffer_Cells_E[cell] = 0;
1354  fBuffer_Cells_RelativeEta[cell] = 0;
1355  fBuffer_Cells_RelativePhi[cell] = 0;
1356  fBuffer_Surrounding_Cells_ID[cell] = 0;
1357  fBuffer_Surrounding_Cells_R[cell] = 0;
1358  fBuffer_Surrounding_Cells_E[cell] = 0;
1359  fBuffer_Surrounding_Cells_RelativeEta[cell] = 0;
1360  fBuffer_Surrounding_Cells_RelativePhi[cell] = 0;
1361  }
1362  for(Int_t track = 0; track < kMaxNTracks; track++){
1363  fBuffer_Surrounding_Tracks_R[track] = 0;
1364  fBuffer_Surrounding_Tracks_Pt[track] = 0;
1365  fBuffer_Surrounding_Tracks_RelativeEta[track] = 0;
1366  fBuffer_Surrounding_Tracks_RelativePhi[track] = 0;
1367  }
1368 }
1369 
AliConvEventCuts::SetPeriodEnum
void SetPeriodEnum(TString periodName)
Definition: AliConvEventCuts.cxx:5680
AliAnalysisTaskClusterQA::fBuffer_Surrounding_NTracks
Int_t fBuffer_Surrounding_NTracks
! array buffer
Definition: AliAnalysisTaskClusterQA.h:182
AliAnalysisTaskClusterQA.h
AliAODConversionPhoton::SetCaloClusterRef
void SetCaloClusterRef(Long_t ref)
Definition: AliAODConversionPhoton.h:71
Double_t
double Double_t
Definition: External.C:58
AliAnalysisTaskClusterQA::fBuffer_ClusterPhi
Float_t fBuffer_ClusterPhi
! array buffer
Definition: AliAnalysisTaskClusterQA.h:153
AliAnalysisTaskClusterQA::UserExec
virtual void UserExec(Option_t *option)
Definition: AliAnalysisTaskClusterQA.cxx:376
TH2F
Definition: External.C:236
AliAODConversionPhoton::SetCaloPhotonMCFlags
void SetCaloPhotonMCFlags(AliMCEvent *mcEvent, Bool_t enableSort)
Definition: AliAODConversionPhoton.cxx:153
AliAnalysisTaskClusterQA::GetUniqueEventID
ULong64_t GetUniqueEventID(AliVHeader *header)
Definition: AliAnalysisTaskClusterQA.cxx:1303
AliAnalysisTaskClusterQA::fBuffer_Event_Vertex_Y
Float_t fBuffer_Event_Vertex_Y
! array buffer
Definition: AliAnalysisTaskClusterQA.h:166
AliAnalysisTaskClusterQA::MakePhotonCandidates
Int_t MakePhotonCandidates(AliVCluster *clus, AliVCaloCells *cells, Long_t indexCluster)
Definition: AliAnalysisTaskClusterQA.cxx:697
AliConvEventCuts::GetPeriodEnum
PeriodVar GetPeriodEnum()
Definition: AliConvEventCuts.h:494
AliAnalysisTaskClusterQA::fSaveMCInformation
Bool_t fSaveMCInformation
save MC information
Definition: AliAnalysisTaskClusterQA.h:143
AliConvEventCuts::SetPeriodEnumExplicit
void SetPeriodEnumExplicit(PeriodVar periodEnum)
Definition: AliConvEventCuts.h:352
kMaxNTracks
const Int_t kMaxNTracks
Definition: AliAnalysisTaskClusterQA.h:208
AliAnalysisTaskClusterQA::fBuffer_ClusterE
Float_t fBuffer_ClusterE
! array buffer
Definition: AliAnalysisTaskClusterQA.h:152
centrality
centrality
Definition: AverageDmesonRaa.C:72
AliAnalysisTaskClusterQA::fBuffer_Cells_RelativeEta
Float_t * fBuffer_Cells_RelativeEta
! array buffer
Definition: AliAnalysisTaskClusterQA.h:173
AliAnalysisTaskClusterQA::fSaveEventProperties
Bool_t fSaveEventProperties
save general event properties (centrality etc.)
Definition: AliAnalysisTaskClusterQA.h:136
AliAnalysisTaskClusterQA::fConeRadius
Float_t fConeRadius
save arrays of all cells in event
Definition: AliAnalysisTaskClusterQA.h:140
AliV0ReaderV1::GetPeriodName
TString GetPeriodName()
Definition: AliV0ReaderV1.h:129
AliAODConversionPhoton::SetCaloPhotonMCLabel
void SetCaloPhotonMCLabel(Int_t i, Int_t labelCaloPhoton)
Definition: AliAODConversionPhoton.h:62
AliAnalysisTaskClusterQA::fBuffer_Mother_MC_Label
Int_t fBuffer_Mother_MC_Label
! array buffer
Definition: AliAnalysisTaskClusterQA.h:189
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Tracks_R
Float_t * fBuffer_Surrounding_Tracks_R
! array buffer
Definition: AliAnalysisTaskClusterQA.h:183
AliAnalysisTaskClusterQA::fBuffer_Event_Vertex_Z
Float_t fBuffer_Event_Vertex_Z
! array buffer
Definition: AliAnalysisTaskClusterQA.h:167
AliAODConversionPhoton::GetCaloPhotonMotherMCLabel
Int_t GetCaloPhotonMotherMCLabel(Int_t i)
Definition: AliAODConversionPhoton.h:67
AliAnalysisTaskClusterQA::fBuffer_ClusterIsEMCAL
Bool_t fBuffer_ClusterIsEMCAL
! array buffer
Definition: AliAnalysisTaskClusterQA.h:155
AliAnalysisTaskClusterQA::fBuffer_Event_Vertex_X
Float_t fBuffer_Event_Vertex_X
! array buffer
Definition: AliAnalysisTaskClusterQA.h:165
AliAnalysisTaskClusterQA::fBuffer_LeadingCell_ID
Int_t fBuffer_LeadingCell_ID
! array buffer
Definition: AliAnalysisTaskClusterQA.h:158
AliAnalysisTaskClusterQA::GetRowAndColumnFromAbsCellID
void GetRowAndColumnFromAbsCellID(Int_t cellIndex, Int_t &row, Int_t &column)
Definition: AliAnalysisTaskClusterQA.cxx:685
AliAnalysisTaskClusterQA::ProcessTracksAndMatching
void ProcessTracksAndMatching(AliVCluster *clus, Long_t indexCluster)
Definition: AliAnalysisTaskClusterQA.cxx:797
AliAnalysisTaskClusterQA::fBuffer_Event_Multiplicity
Float_t fBuffer_Event_Multiplicity
! array buffer
Definition: AliAnalysisTaskClusterQA.h:168
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Cells_ID
Int_t * fBuffer_Surrounding_Cells_ID
! array buffer
Definition: AliAnalysisTaskClusterQA.h:177
AliAnalysisTaskClusterQA::FindLargestCellInCluster
Int_t FindLargestCellInCluster(AliVCluster *cluster, AliVEvent *event)
Definition: AliAnalysisTaskClusterQA.cxx:660
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Cells_E
Float_t * fBuffer_Surrounding_Cells_E
! array buffer
Definition: AliAnalysisTaskClusterQA.h:179
AliAnalysisTaskClusterQA::fBuffer_LeadingCell_Phi
Float_t fBuffer_LeadingCell_Phi
! array buffer
Definition: AliAnalysisTaskClusterQA.h:161
AliConvEventCuts::IsEventAcceptedByCut
Int_t IsEventAcceptedByCut(AliConvEventCuts *ReaderCuts, AliVEvent *event, AliMCEvent *mcEvent, Int_t isHeavyIon, Bool_t isEMCALAnalysis)
Definition: AliConvEventCuts.cxx:4692
Int_t
int Int_t
Definition: External.C:63
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Cells_R
Float_t * fBuffer_Surrounding_Cells_R
! array buffer
Definition: AliAnalysisTaskClusterQA.h:178
AliCaloPhotonCuts
Class handling all kinds of selection cuts for Gamma Calo analysis.
Definition: AliCaloPhotonCuts.h:77
AliAnalysisTaskClusterQA::ProcessQATreeCluster
void ProcessQATreeCluster(AliVEvent *event, AliVCluster *cluster, Long_t indexCluster)
Definition: AliAnalysisTaskClusterQA.cxx:509
AliAnalysisTaskClusterQA::fBuffer_Cells_RelativePhi
Float_t * fBuffer_Cells_RelativePhi
! array buffer
Definition: AliAnalysisTaskClusterQA.h:174
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Tracks_RelativeEta
Float_t * fBuffer_Surrounding_Tracks_RelativeEta
! array buffer
Definition: AliAnalysisTaskClusterQA.h:185
Float_t
float Float_t
Definition: External.C:68
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Tracks_Pt
Float_t * fBuffer_Surrounding_Tracks_Pt
! array buffer
Definition: AliAnalysisTaskClusterQA.h:184
AliConvEventCuts::IsConversionPrimaryESD
Bool_t IsConversionPrimaryESD(AliMCEvent *mcEvent, Long_t eventpos, Double_t prodVtxX, Double_t prodVtxY, Double_t prodVtxZ)
Definition: AliConvEventCuts.cxx:5347
AliAnalysisTaskClusterQA::fBuffer_Cluster_MC_Label
Int_t fBuffer_Cluster_MC_Label
! array buffer
Definition: AliAnalysisTaskClusterQA.h:188
AliAnalysisTaskClusterQA::fSaveSurroundingCells
Bool_t fSaveSurroundingCells
save arrays of all cells in event
Definition: AliAnalysisTaskClusterQA.h:138
AliAnalysisTaskClusterQA::fBuffer_ClusterM20
Float_t fBuffer_ClusterM20
! array buffer
Definition: AliAnalysisTaskClusterQA.h:163
AliAnalysisTaskClusterQA::fSaveAdditionalHistos
Bool_t fSaveAdditionalHistos
save MC information
Definition: AliAnalysisTaskClusterQA.h:144
AliAnalysisTaskClusterQA::fMinClusterEnergy
Float_t fMinClusterEnergy
save arrays of all cells in event
Definition: AliAnalysisTaskClusterQA.h:142
AliAnalysisTaskClusterQA::fBuffer_Cells_ID
Int_t * fBuffer_Cells_ID
! array buffer
Definition: AliAnalysisTaskClusterQA.h:171
AliAnalysisTaskClusterQA::fMaxNLMCut
Int_t fMaxNLMCut
save MC information
Definition: AliAnalysisTaskClusterQA.h:122
Data
Bool_t Data(TH1F *h, Double_t *rangefit, Bool_t writefit, Double_t &sgn, Double_t &errsgn, Double_t &bkg, Double_t &errbkg, Double_t &sgnf, Double_t &errsgnf, Double_t &sigmafit, Int_t &status)
Definition: charmCutsOptimization.C:371
Short_t
short Short_t
Definition: External.C:23
AliAnalysisTaskClusterQA::fBuffer_ClusterM02
Float_t fBuffer_ClusterM02
! array buffer
Definition: AliAnalysisTaskClusterQA.h:162
AliAnalysisTaskClusterQA::Terminate
virtual void Terminate(Option_t *)
Definition: AliAnalysisTaskClusterQA.cxx:655
AliAnalysisTaskClusterQA::fBuffer_Event_NumActiveCells
Int_t fBuffer_Event_NumActiveCells
! array buffer
Definition: AliAnalysisTaskClusterQA.h:169
AliAnalysisTaskClusterQA::fBuffer_ClusterEta
Float_t fBuffer_ClusterEta
! array buffer
Definition: AliAnalysisTaskClusterQA.h:154
AliAnalysisTaskClusterQA::GetCentrality
Float_t GetCentrality(AliVEvent *event)
Definition: AliAnalysisTaskClusterQA.cxx:1313
AliAnalysisTaskClusterQA::fSaveCells
Bool_t fSaveCells
save arrays of cluster cells
Definition: AliAnalysisTaskClusterQA.h:137
AliEMCALRecoUtils::ExtrapolateTrackToCluster
static Bool_t ExtrapolateTrackToCluster(AliExternalTrackParam *trkParam, const AliVCluster *cluster, Double_t mass, Double_t step, Float_t &tmpEta, Float_t &tmpPhi)
Definition: AliEMCALRecoUtils.h:315
AliAnalysisTaskClusterQA::fBuffer_EventID
ULong64_t fBuffer_EventID
! array buffer
Definition: AliAnalysisTaskClusterQA.h:151
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Cells_RelativePhi
Float_t * fBuffer_Surrounding_Cells_RelativePhi
! array buffer
Definition: AliAnalysisTaskClusterQA.h:181
AliConversionMesonCuts
Class handling all kinds of selection cuts for Gamma Conversion analysis.
Definition: AliConversionMesonCuts.h:62
AliAnalysisTaskClusterQA::fBuffer_ClusterNumCells
Int_t fBuffer_ClusterNumCells
! array buffer
Definition: AliAnalysisTaskClusterQA.h:157
AliAODConversionPhoton::SetNCaloPhotonMCLabels
void SetNCaloPhotonMCLabels(Int_t nLabels)
Definition: AliAODConversionPhoton.h:64
AliAnalysisTaskClusterQA::fBuffer_LeadingCell_Eta
Float_t fBuffer_LeadingCell_Eta
! array buffer
Definition: AliAnalysisTaskClusterQA.h:160
kMaxActiveCells
const Int_t kMaxActiveCells
Definition: AliAnalysisTaskClusterQA.h:207
AliConvEventCuts
Class handling all kinds of selection cuts for Gamma Conversion analysis.
Definition: AliConvEventCuts.h:55
AliAnalysisTaskClusterQA::fBuffer_MC_Cluster_Flag
Int_t fBuffer_MC_Cluster_Flag
! array buffer
Definition: AliAnalysisTaskClusterQA.h:156
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Cells_RelativeEta
Float_t * fBuffer_Surrounding_Cells_RelativeEta
! array buffer
Definition: AliAnalysisTaskClusterQA.h:180
AliAnalysisTaskClusterQA::ProcessTrueClusterCandidatesAOD
Int_t ProcessTrueClusterCandidatesAOD(AliAODConversionPhoton *TrueClusterCandidate, Float_t m02, AliAODConversionPhoton *TrueSubClusterCandidate1, AliAODConversionPhoton *TrueSubClusterCandidate2)
Definition: AliAnalysisTaskClusterQA.cxx:1126
AliAnalysisTaskClusterQA::fMesonCuts
AliConversionMesonCuts * fMesonCuts
Definition: AliAnalysisTaskClusterQA.h:120
AliAnalysisTaskClusterQA::fBuffer_LeadingCell_E
Float_t fBuffer_LeadingCell_E
! array buffer
Definition: AliAnalysisTaskClusterQA.h:159
AliV0ReaderV1::GetEventCuts
AliConvEventCuts * GetEventCuts()
Definition: AliV0ReaderV1.h:90
Option_t
const char Option_t
Definition: External.C:48
Bool_t
bool Bool_t
Definition: External.C:53
AliAnalysisTaskClusterQA::fBuffer_Surrounding_Tracks_RelativePhi
Float_t * fBuffer_Surrounding_Tracks_RelativePhi
! array buffer
Definition: AliAnalysisTaskClusterQA.h:186
AliAODConversionPhoton::SetCaloPhotonMCFlagsAOD
void SetCaloPhotonMCFlagsAOD(AliVEvent *event, Bool_t enableSort)
Definition: AliAODConversionPhoton.cxx:516
AliAnalysisTaskClusterQA::fBuffer_Cells_E
Float_t * fBuffer_Cells_E
! array buffer
Definition: AliAnalysisTaskClusterQA.h:172
AliAnalysisTaskClusterQA::fMinNLMCut
Int_t fMinNLMCut
save MC information
Definition: AliAnalysisTaskClusterQA.h:121
AliAnalysisTaskClusterQA::fBuffer_Surrounding_NCells
Int_t fBuffer_Surrounding_NCells
! array buffer
Definition: AliAnalysisTaskClusterQA.h:176
AliAnalysisTaskClusterQA::fSaveTracks
Bool_t fSaveTracks
save arrays of all cells in event
Definition: AliAnalysisTaskClusterQA.h:139
AliAnalysisTaskClusterQA::ProcessTrueClusterCandidates
Int_t ProcessTrueClusterCandidates(AliAODConversionPhoton *TrueClusterCandidate, Float_t m02, AliAODConversionPhoton *TrueSubClusterCandidate1, AliAODConversionPhoton *TrueSubClusterCandidate2)
Definition: AliAnalysisTaskClusterQA.cxx:940
OpenFile
TList * OpenFile(const char *fname)
Definition: DrawAnaELoss.C:65
AliAnalysisTaskClusterQA::fMinTrackPt
Float_t fMinTrackPt
save arrays of all cells in event
Definition: AliAnalysisTaskClusterQA.h:141