AliPhysics  58f3d52 (58f3d52)
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
AliConvEventCuts.cxx
Go to the documentation of this file.
1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Authors: Friederike Bock, Daniel Muehlheim *
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 handling all kinds of selection cuts for
19 // Gamma Conversion analysis
20 //---------------------------------------------
22 
23 #include "AliConvEventCuts.h"
24 
25 #include "AliAODTrack.h"
26 #include "AliESDtrack.h"
27 #include "AliAnalysisManager.h"
28 #include "AliInputEventHandler.h"
29 #include "AliMCEventHandler.h"
30 #include "AliAODHandler.h"
31 #include "TH1.h"
32 #include "TH2.h"
33 #include "TF1.h"
34 #include "TObjString.h"
35 #include "AliMCEvent.h"
36 #include "AliAODEvent.h"
37 #include "AliESDEvent.h"
38 #include "AliCentrality.h"
39 #include "AliMultSelection.h"
40 #include "TList.h"
41 #include "TFile.h"
42 #include "AliLog.h"
43 #include "AliGenCocktailEventHeader.h"
44 #include "AliGenDPMjetEventHeader.h"
45 #include "AliGenPythiaEventHeader.h"
46 #include "AliGenHijingEventHeader.h"
47 #include "AliTriggerAnalysis.h"
48 #include "AliV0ReaderV1.h"
49 #include "AliVCaloCells.h"
50 #include "AliAODMCParticle.h"
51 #include "AliAODMCHeader.h"
52 #include "AliEMCALTriggerPatchInfo.h"
53 
54 class iostream;
55 
56 using namespace std;
57 
59 ClassImp(AliConvEventCuts)
61 
62 
63 const char* AliConvEventCuts::fgkCutNames[AliConvEventCuts::kNCuts] = {
64  "HeavyIon", //0
65  "CentralityMin", //1
66  "CentralityMax", //2
67  "SelectSpecialTrigger", //3
68  "SelectSpecialSubTriggerClass", //4
69  "RemovePileUp", //5
70  "RejectExtraSignals", //6
71  "VertexCut", //7
72 };
73 
74 
75 //________________________________________________________________________
76 AliConvEventCuts::AliConvEventCuts(const char *name,const char *title) :
77  AliAnalysisCuts(name,title),
78  fHistograms(NULL),
79  fHeaderList(NULL),
80  fDoLightOutput(kFALSE),
81  fEventQuality(-1),
82  fIsHeavyIon(0),
83  fDetectorCentrality(0),
84  fModCentralityClass(0),
85  fEnableVertexCut(kTRUE),
86  fMaxVertexZ(10),
87  fCentralityMin(0),
88  fCentralityMax(0),
89  fMultiplicityMethod(0),
90  fSpecialTrigger(0),
91  fSpecialSubTrigger(0),
92  fRemovePileUp(kFALSE),
93  fPastFutureRejectionLow(0),
94  fPastFutureRejectionHigh(0),
95  fDoPileUpRejectV0MTPCout(0),
96  fFPileUpRejectV0MTPCout(0),
97  fRejectExtraSignals(0),
98  fOfflineTriggerMask(0),
99  fHasV0AND(kTRUE),
100  fIsSDDFired(kTRUE),
101  fRandom(0),
102  fnHeaders(0),
103  fNotRejectedStart(NULL),
104  fNotRejectedEnd(NULL),
105  fGeneratorNames(NULL),
106  fPeriodEnum(kNoPeriod),
107  fEnergyEnum(kUnset),
108  fCutString(NULL),
109  fCutStringRead(""),
110  fUtils(NULL),
111  fEtaShift(0.0),
112  fDoEtaShift(kFALSE),
113  fDoCentralityFlat(0),
114  fPathWeightsFlatCent(""),
115  fNameHistoNotFlatCentrality(""),
116  fDoReweightHistoMCPi0(kFALSE),
117  fDoReweightHistoMCEta(kFALSE),
118  fDoReweightHistoMCK0s(kFALSE),
119  fPathTrFReweighting(""),
120  fNameHistoReweightingPi0(""),
121  fNameHistoReweightingEta(""),
122  fNameHistoReweightingK0s(""),
123  fNameFitDataPi0(""),
124  fNameFitDataEta(""),
125  fNameFitDataK0s(""),
126  fHistoEventCuts(NULL),
127  fHistoPastFutureBits(NULL),
128  hCentrality(NULL),
129  hCentralityNotFlat(NULL),
130  //hCentralityVsNumberOfPrimaryTracks(NULL),
131  hVertexZ(NULL),
132  hNPileupVertices(NULL),
133  hPileupVertexToPrimZ(NULL),
134  hPileupVertexToPrimZSPDPileup(NULL),
135  hPileupVertexToPrimZTrackletvsHits(NULL),
136  hEventPlaneAngle(NULL),
137  fEventPlaneAngle(0),
138  hTriggerClass(NULL),
139  hTriggerClassSelected(NULL),
140  hTriggerClassesCorrelated(NULL),
141  hReweightMCHistPi0(NULL),
142  hReweightMCHistEta(NULL),
143  hReweightMCHistK0s(NULL),
144  fFitDataPi0(NULL),
145  fFitDataEta(NULL),
146  fFitDataK0s(NULL),
147  fAddedSignalPDGCode(0),
148  fPreSelCut(kFALSE),
149  fTriggerSelectedManually(kFALSE),
150  fSpecialTriggerName(""),
151  fSpecialSubTriggerName(""),
152  fNSpecialSubTriggerOptions(0),
153  hSPDClusterTrackletBackgroundBefore(NULL),
154  hSPDClusterTrackletBackground(NULL),
155  fV0ReaderName(""),
156  fCaloTriggers(NULL),
157  fTriggerPatchInfo(NULL),
158  fMainTriggerPatchEMCAL(NULL),
159  fCaloTriggersName(""),
160  fCaloTriggerPatchInfoName(""),
161  fTriggersEMCAL(0),
162  fTriggersEMCALSelected(-1),
163  fEMCALTrigInitialized(kFALSE),
164  fSecProdBoundary(1.0),
165  fMaxPtJetMC(0),
166  fMaxFacPtHard(2.5),
167  fMaxFacPtHardSingleParticle(1.5),
168  fMimicTrigger(kFALSE),
169  fRejectTriggerOverlap(kFALSE),
170  fDoMultiplicityWeighting(kFALSE),
171  fPathReweightingMult(""),
172  fNameHistoReweightingMultData(""),
173  fNameHistoReweightingMultMC(""),
174  hReweightMultData(NULL),
175  hReweightMultMC(NULL),
176  fDebugLevel(0)
177 {
178  for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=0;}
179  fCutString=new TObjString((GetCutNumber()).Data());
180 
181  fUtils = new AliAnalysisUtils();
182  //if you do not want to apply the cut on the distance between the SPD and TRK vertex:
183  //fUtils->SetCutOnZVertexSPD(kFALSE);
184 
185 
186 }
187 
188 //________________________________________________________________________
190  AliAnalysisCuts(ref),
191  fHistograms(NULL),
192  fHeaderList(ref.fHeaderList),
193  fDoLightOutput(ref.fDoLightOutput),
194  fEventQuality(ref.fEventQuality),
195  fIsHeavyIon(ref.fIsHeavyIon),
196  fDetectorCentrality(ref.fDetectorCentrality),
197  fModCentralityClass(ref.fModCentralityClass),
198  fEnableVertexCut(ref.fEnableVertexCut),
199  fMaxVertexZ(ref.fMaxVertexZ),
200  fCentralityMin(ref.fCentralityMin),
201  fCentralityMax(ref.fCentralityMax),
202  fMultiplicityMethod(ref.fMultiplicityMethod),
203  fSpecialTrigger(ref.fSpecialTrigger),
204  fSpecialSubTrigger(ref.fSpecialSubTrigger),
205  fRemovePileUp(ref.fRemovePileUp),
206  fPastFutureRejectionLow(ref.fPastFutureRejectionLow),
207  fPastFutureRejectionHigh(ref.fPastFutureRejectionHigh),
208  fDoPileUpRejectV0MTPCout(ref.fDoPileUpRejectV0MTPCout),
209  fFPileUpRejectV0MTPCout(ref.fFPileUpRejectV0MTPCout),
210  fRejectExtraSignals(ref.fRejectExtraSignals),
211  fOfflineTriggerMask(ref.fOfflineTriggerMask),
212  fHasV0AND(ref.fHasV0AND),
213  fIsSDDFired(ref.fIsSDDFired),
214  fRandom(ref.fRandom),
215  fnHeaders(ref.fnHeaders),
216  fNotRejectedStart(NULL),
217  fNotRejectedEnd(NULL),
218  fGeneratorNames(ref.fGeneratorNames),
219  fPeriodEnum(ref.fPeriodEnum),
220  fEnergyEnum(kUnset),
221  fCutString(NULL),
222  fCutStringRead(""),
223  fUtils(NULL),
224  fEtaShift(ref.fEtaShift),
225  fDoEtaShift(ref.fDoEtaShift),
226  fDoCentralityFlat(ref.fDoCentralityFlat),
227  fPathWeightsFlatCent(ref.fPathWeightsFlatCent),
228  fNameHistoNotFlatCentrality(ref.fNameHistoNotFlatCentrality),
229  fDoReweightHistoMCPi0(ref.fDoReweightHistoMCPi0),
230  fDoReweightHistoMCEta(ref.fDoReweightHistoMCEta),
231  fDoReweightHistoMCK0s(ref.fDoReweightHistoMCK0s),
232  fPathTrFReweighting(ref.fPathTrFReweighting),
233  fNameHistoReweightingPi0(ref.fNameHistoReweightingPi0),
234  fNameHistoReweightingEta(ref.fNameHistoReweightingEta),
235  fNameHistoReweightingK0s(ref.fNameHistoReweightingK0s),
236  fNameFitDataPi0(ref.fNameFitDataPi0),
237  fNameFitDataEta(ref.fNameFitDataEta),
238  fNameFitDataK0s(ref.fNameFitDataK0s),
239  fHistoEventCuts(NULL),
240  fHistoPastFutureBits(NULL),
241  hCentrality(ref.hCentrality),
242  hCentralityNotFlat(ref.hCentralityNotFlat),
243  //hCentralityVsNumberOfPrimaryTracks(ref.hCentralityVsNumberOfPrimaryTracks),
244  hVertexZ(ref.hVertexZ),
245  hNPileupVertices(ref.hNPileupVertices),
246  hPileupVertexToPrimZ(ref.hPileupVertexToPrimZ),
247  hPileupVertexToPrimZSPDPileup(ref.hPileupVertexToPrimZSPDPileup),
248  hPileupVertexToPrimZTrackletvsHits(ref.hPileupVertexToPrimZTrackletvsHits),
249  hEventPlaneAngle(ref.hEventPlaneAngle),
250  fEventPlaneAngle(ref.fEventPlaneAngle),
251  hTriggerClass(NULL),
252  hTriggerClassSelected(NULL),
253  hTriggerClassesCorrelated(NULL),
254  hReweightMCHistPi0(ref.hReweightMCHistPi0),
255  hReweightMCHistEta(ref.hReweightMCHistEta),
256  hReweightMCHistK0s(ref.hReweightMCHistK0s),
257  fFitDataPi0(ref.fFitDataPi0),
258  fFitDataEta(ref.fFitDataEta),
259  fFitDataK0s(ref.fFitDataK0s),
260  fAddedSignalPDGCode(ref.fAddedSignalPDGCode),
261  fPreSelCut(ref.fPreSelCut),
262  fTriggerSelectedManually(ref.fTriggerSelectedManually),
263  fSpecialTriggerName(ref.fSpecialTriggerName),
264  fSpecialSubTriggerName(ref.fSpecialSubTriggerName),
265  fNSpecialSubTriggerOptions(ref.fNSpecialSubTriggerOptions),
266  hSPDClusterTrackletBackgroundBefore(NULL),
267  hSPDClusterTrackletBackground(NULL),
268  fV0ReaderName(ref.fV0ReaderName),
269  fCaloTriggers(NULL),
270  fTriggerPatchInfo(NULL),
271  fMainTriggerPatchEMCAL(NULL),
272  fCaloTriggersName(ref.fCaloTriggersName),
273  fCaloTriggerPatchInfoName(ref.fCaloTriggerPatchInfoName),
274  fTriggersEMCAL(ref.fTriggersEMCAL),
275  fTriggersEMCALSelected(ref.fTriggersEMCALSelected),
276  fEMCALTrigInitialized(kFALSE),
277  fSecProdBoundary(ref.fSecProdBoundary),
278  fMaxPtJetMC(ref.fMaxPtJetMC),
279  fMaxFacPtHard(ref.fMaxFacPtHard),
280  fMaxFacPtHardSingleParticle(ref.fMaxFacPtHardSingleParticle),
281  fMimicTrigger(ref.fMimicTrigger),
282  fRejectTriggerOverlap(ref.fRejectTriggerOverlap),
283  fDoMultiplicityWeighting(ref.fDoMultiplicityWeighting),
284  fPathReweightingMult(ref.fPathReweightingMult),
285  fNameHistoReweightingMultData(ref.fNameHistoReweightingMultData),
286  fNameHistoReweightingMultMC(ref.fNameHistoReweightingMultMC),
287  hReweightMultData(ref.hReweightMultData),
288  hReweightMultMC(ref.hReweightMultMC),
290 {
291  // Copy Constructor
292  for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=ref.fCuts[jj];}
293  fCutString=new TObjString((GetCutNumber()).Data());
294  fUtils = new AliAnalysisUtils();
295  // dont copy histograms (if you like histograms, call InitCutHistograms())
296 
297 }
298 
299 
300 //________________________________________________________________________
302  // Destructor
303  //Deleting fHistograms leads to seg fault it it's added to output collection of a task
304  // if(fHistograms)
305  // delete fHistograms;
306  // fHistograms = NULL;
307  if(fCutString != NULL){
308  delete fCutString;
309  fCutString = NULL;
310  }
311  if(fNotRejectedStart){
312  delete[] fNotRejectedStart;
313  fNotRejectedStart = NULL;
314  }
315  if(fNotRejectedEnd){
316  delete[] fNotRejectedEnd;
317  fNotRejectedEnd = NULL;
318  }
319  if(fGeneratorNames){
320  delete[] fGeneratorNames;
321  fGeneratorNames = NULL;
322  }
323  if(fUtils){
324  delete fUtils;
325  fUtils = NULL;
326  }
327 
328 }
329 
330 //________________________________________________________________________
332 
333  // Initialize Cut Histograms for QA (only initialized and filled if function is called)
334  TH1::AddDirectory(kFALSE);
335 
336  if(fHistograms != NULL){
337  delete fHistograms;
338  fHistograms=NULL;
339  }
340  if(fHistograms==NULL){
341  fHistograms=new TList();
342  fHistograms->SetOwner(kTRUE);
343  if(name=="")fHistograms->SetName(Form("ConvEventCuts_%s",GetCutNumber().Data()));
344  else fHistograms->SetName(Form("%s_%s",name.Data(),GetCutNumber().Data()));
345  }
346 
347  if (hReweightMCHistPi0){
348  hReweightMCHistPi0->SetName("MCInputForWeightingPi0");
350  }
351  if (hReweightMCHistEta){
352  hReweightMCHistEta->SetName("MCInputForWeightingEta");
354  }
355  if (hReweightMCHistK0s){
356  hReweightMCHistK0s->SetName("MCInputForWeightingK0s");
358  }
359 
360  if (hReweightMultData){
361  hReweightMultData->SetName(Form("hReweightMultData_%s",GetCutNumber().Data()));
363  }
364  if (hReweightMultMC){
365  hReweightMultMC->SetName(Form("hReweightMultMC_%s",GetCutNumber().Data()));
367  }
368 
369  if(!fDoLightOutput){
370  hSPDClusterTrackletBackgroundBefore = new TH2F(Form("SPD tracklets vs SPD clusters %s before Pileup Cut",GetCutNumber().Data()),"SPD tracklets vs SPD clusters",100,0,200,250,0,1000);
372 
373  hSPDClusterTrackletBackground = new TH2F(Form("SPD tracklets vs SPD clusters %s",GetCutNumber().Data()),"SPD tracklets vs SPD clusters",100,0,200,250,0,1000);
375  }
376 
377  if(fIsHeavyIon > 0){
378  hCentrality=new TH1F(Form("Centrality %s",GetCutNumber().Data()),"Centrality",400,0,100);
379  fHistograms->Add(hCentrality);
380  }
381 
382  //hCentralityVsNumberOfPrimaryTracks=new TH2F(Form("Centrality vs Primary Tracks %s",GetCutNumber().Data()),"Centrality vs Primary Tracks ",400,0,100,4000,0,4000);
383  //fHistograms->Add(hCentralityVsNumberOfPrimaryTracks); commented on 3.3.2015 because it's in the main Task
384 
385  hVertexZ = new TH1F(Form("VertexZ %s",GetCutNumber().Data()),"VertexZ",1000,-50,50);
386  fHistograms->Add(hVertexZ);
387 
388  hNPileupVertices = new TH1F(Form("NPileupVertices %s",GetCutNumber().Data()),"NPileupVertices",30,-0.5,29.5);
390 
391  hPileupVertexToPrimZ = new TH1F(Form("PileupVertexDistance %s",GetCutNumber().Data()),"PileupVertexDistance",600,-15,15);
393  hPileupVertexToPrimZSPDPileup = new TH1F(Form("PileupVertexDistance_SPDPileup %s",GetCutNumber().Data()),"PileupVertexDistance_SPDPileup",600,-15,15);
395  hPileupVertexToPrimZTrackletvsHits = new TH1F(Form("PileupVertexDistance_TrackletvsHits %s",GetCutNumber().Data()),"PileupVertexDistance_TrackletvsHits",600,-15,15);
397 
398  if(fIsHeavyIon == 1){
399  hEventPlaneAngle = new TH1F(Form("EventPlaneAngle %s",GetCutNumber().Data()),"EventPlaneAngle",60, 0, TMath::Pi());
401  }
402  fHistoPastFutureBits=new TH1F(Form("PastFutureBits %s",GetCutNumber().Data()),"Past Future Bits",180,-90*25,90*25);
404 
405  // Event Cuts and Info
406  if(preCut){
407  fHistoEventCuts=new TH1F(Form("ESD_EventCuts %s",GetCutNumber().Data()),"Event Cuts",8,-0.5,7.5);
408  fHistoEventCuts->GetXaxis()->SetBinLabel(1,"in");
409  fHistoEventCuts->GetXaxis()->SetBinLabel(2,"OfflineTrigger");
410  fHistoEventCuts->GetXaxis()->SetBinLabel(3,"nvtxcontr");
411  fHistoEventCuts->GetXaxis()->SetBinLabel(4,"VertexZ");
412  fHistoEventCuts->GetXaxis()->SetBinLabel(5,"pileup");
413  fHistoEventCuts->GetXaxis()->SetBinLabel(6,"centrsel");
414  fHistoEventCuts->GetXaxis()->SetBinLabel(7,"OOB-pileup");
415  fHistoEventCuts->GetXaxis()->SetBinLabel(8,"out");
417 
418  hTriggerClass= new TH1F(Form("OfflineTrigger %s",GetCutNumber().Data()),"OfflineTrigger",36,-0.5,35.5);
419  hTriggerClass->GetXaxis()->SetBinLabel( 1,"kMB");
420  hTriggerClass->GetXaxis()->SetBinLabel( 2,"kINT7");
421  hTriggerClass->GetXaxis()->SetBinLabel( 3,"kMUON");
422  hTriggerClass->GetXaxis()->SetBinLabel( 4,"kHighMult");
423  hTriggerClass->GetXaxis()->SetBinLabel( 5,"kEMC1");
424  hTriggerClass->GetXaxis()->SetBinLabel( 6,"kCINT5");
425  hTriggerClass->GetXaxis()->SetBinLabel( 7,"kCMUS5/kMUSPB");
426  hTriggerClass->GetXaxis()->SetBinLabel( 8,"kMUSH7/kMUSHPB");
427  hTriggerClass->GetXaxis()->SetBinLabel( 9,"kMUL7/kMuonLikePB");
428  hTriggerClass->GetXaxis()->SetBinLabel(10,"kMUU7/kMuonUnlikePB");
429  hTriggerClass->GetXaxis()->SetBinLabel(11,"kEMC7/kEMC8");
430  hTriggerClass->GetXaxis()->SetBinLabel(12,"kMUS7");
431  hTriggerClass->GetXaxis()->SetBinLabel(13,"kPHI1");
432  hTriggerClass->GetXaxis()->SetBinLabel(14,"kPHI7/kPHI8/kPHOSPb");
433  hTriggerClass->GetXaxis()->SetBinLabel(15,"kEMCEJE");
434  hTriggerClass->GetXaxis()->SetBinLabel(16,"kEMCEGA");
435  hTriggerClass->GetXaxis()->SetBinLabel(17,"kCentral");
436  hTriggerClass->GetXaxis()->SetBinLabel(18,"kSemiCentral");
437  hTriggerClass->GetXaxis()->SetBinLabel(19,"kDG5");
438  hTriggerClass->GetXaxis()->SetBinLabel(20,"kZED");
439  hTriggerClass->GetXaxis()->SetBinLabel(21,"kSPI7/kSPI");
440  hTriggerClass->GetXaxis()->SetBinLabel(22,"kINT8");
441  hTriggerClass->GetXaxis()->SetBinLabel(23,"kMuonSingleLowPt8");
442  hTriggerClass->GetXaxis()->SetBinLabel(24,"kMuonSingleHighPt8");
443  hTriggerClass->GetXaxis()->SetBinLabel(25,"kMuonLikeLowPt8");
444  hTriggerClass->GetXaxis()->SetBinLabel(26,"kMuonUnlikeLowPt8");
445  hTriggerClass->GetXaxis()->SetBinLabel(27,"kMuonUnlikeLowPt0");
446  hTriggerClass->GetXaxis()->SetBinLabel(28,"kUserDefined");
447  hTriggerClass->GetXaxis()->SetBinLabel(29,"kTRD");
448  hTriggerClass->GetXaxis()->SetBinLabel(30,"kFastOnly");
449  hTriggerClass->GetXaxis()->SetBinLabel(31,"kAnyINT");
450  hTriggerClass->GetXaxis()->SetBinLabel(32,"kAny");
451  hTriggerClass->GetXaxis()->SetBinLabel(33,"V0AND");
452  hTriggerClass->GetXaxis()->SetBinLabel(34,"NOT kFastOnly");
453  hTriggerClass->GetXaxis()->SetBinLabel(35,"failed Physics Selection");
454  hTriggerClass->GetXaxis()->SetBinLabel(36,"mimickedTrigger");
456  }
457  if(!preCut){
458  hTriggerClassSelected= new TH1F(Form("OfflineTriggerSelected %s",GetCutNumber().Data()),"OfflineTriggerSelected",35,-0.5,34.5);
459  hTriggerClassSelected->GetXaxis()->SetBinLabel( 1,"kMB");
460  hTriggerClassSelected->GetXaxis()->SetBinLabel( 2,"kINT7");
461  hTriggerClassSelected->GetXaxis()->SetBinLabel( 3,"kMUON");
462  hTriggerClassSelected->GetXaxis()->SetBinLabel( 4,"kHighMult");
463  hTriggerClassSelected->GetXaxis()->SetBinLabel( 5,"kEMC1");
464  hTriggerClassSelected->GetXaxis()->SetBinLabel( 6,"kCINT5");
465  hTriggerClassSelected->GetXaxis()->SetBinLabel( 7,"kCMUS5/kMUSPB");
466  hTriggerClassSelected->GetXaxis()->SetBinLabel( 8,"kMUSH7/kMUSHPB");
467  hTriggerClassSelected->GetXaxis()->SetBinLabel( 9,"kMUL7/kMuonLikePB");
468  hTriggerClassSelected->GetXaxis()->SetBinLabel(10,"kMUU7/kMuonUnlikePB");
469  hTriggerClassSelected->GetXaxis()->SetBinLabel(11,"kEMC7/kEMC8");
470  hTriggerClassSelected->GetXaxis()->SetBinLabel(12,"kMUS7");
471  hTriggerClassSelected->GetXaxis()->SetBinLabel(13,"kPHI1");
472  hTriggerClassSelected->GetXaxis()->SetBinLabel(14,"kPHI7/kPHI8/kPHOSPb");
473  hTriggerClassSelected->GetXaxis()->SetBinLabel(15,"kEMCEJE");
474  hTriggerClassSelected->GetXaxis()->SetBinLabel(16,"kEMCEGA");
475  hTriggerClassSelected->GetXaxis()->SetBinLabel(17,"kCentral");
476  hTriggerClassSelected->GetXaxis()->SetBinLabel(18,"kSemiCentral");
477  hTriggerClassSelected->GetXaxis()->SetBinLabel(19,"kDG5");
478  hTriggerClassSelected->GetXaxis()->SetBinLabel(20,"kZED");
479  hTriggerClassSelected->GetXaxis()->SetBinLabel(21,"kSPI7/kSPI");
480  hTriggerClassSelected->GetXaxis()->SetBinLabel(22,"kINT8");
481  hTriggerClassSelected->GetXaxis()->SetBinLabel(23,"kMuonSingleLowPt8");
482  hTriggerClassSelected->GetXaxis()->SetBinLabel(24,"kMuonSingleHighPt8");
483  hTriggerClassSelected->GetXaxis()->SetBinLabel(25,"kMuonLikeLowPt8");
484  hTriggerClassSelected->GetXaxis()->SetBinLabel(26,"kMuonUnlikeLowPt8");
485  hTriggerClassSelected->GetXaxis()->SetBinLabel(27,"kMuonUnlikeLowPt0");
486  hTriggerClassSelected->GetXaxis()->SetBinLabel(28,"kUserDefined");
487  hTriggerClassSelected->GetXaxis()->SetBinLabel(29,"kTRD");
488  hTriggerClassSelected->GetXaxis()->SetBinLabel(30,"kFastOnly");
489  hTriggerClassSelected->GetXaxis()->SetBinLabel(31,"kAnyINT");
490  hTriggerClassSelected->GetXaxis()->SetBinLabel(32,"kAny");
491  hTriggerClassSelected->GetXaxis()->SetBinLabel(33,"V0AND");
492  hTriggerClassSelected->GetXaxis()->SetBinLabel(34,"NOT kFastOnly");
493  hTriggerClassSelected->GetXaxis()->SetBinLabel(35,"mimickedTrigger");
495 
496  if (fSpecialTrigger == 5 || fSpecialTrigger == 8 || fSpecialTrigger == 9){
497  hTriggerClassesCorrelated= new TH1F(Form("TriggerCorrelations %s",GetCutNumber().Data()),"Triggers Correlated with EMCal triggers",10,-0.5,9.5);
498  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 1,"kMB");
499  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 2,"kINT7");
500  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 3,"kEMC1");
501  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 4,"kEMC7");
502  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 5,"kEMCEJE");
503  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 6,"kEMCEJ1");
504  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 7,"kEMCEJ2");
505  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 8,"kEMCEGA");
506  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 9,"kEMCEG1");
507  hTriggerClassesCorrelated->GetXaxis()->SetBinLabel( 10,"kEMCEG2");
509  }
510 
511  }
512  TH1::AddDirectory(kTRUE);
513 }
514 
516 Bool_t AliConvEventCuts::EventIsSelected(AliVEvent *event, AliMCEvent *mcEvent){
517  // Process Event Selection
518 
519  Int_t cutindex=0;
520  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
521  cutindex++;
522 
523  // Check for MC event
524  Bool_t isMC = kFALSE;
525  if(mcEvent && event->IsA()==AliESDEvent::Class()){
526  // Check if MC event is correctly loaded
527  AliMCEventHandler* mcHandler = dynamic_cast<AliMCEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
528  if (!mcHandler){
529  fEventQuality = 2;
530  return kFALSE;
531  }
532  if (!mcHandler->InitOk() ){
533  fEventQuality = 2;
534  return kFALSE;
535  }
536  if (!mcHandler->TreeK() ){
537  fEventQuality = 2;
538  return kFALSE;
539  }
540  if (!mcHandler->TreeTR() ) {
541  fEventQuality = 2;
542  return kFALSE;
543  }
544  isMC = kTRUE;
545  }
546 
547 
548 
549  // Event Trigger
550  // cout << "before event trigger" << endl;
551  if(!IsTriggerSelected(event, isMC )){
552  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
553  fEventQuality = 3;
554  return kFALSE;
555  }
556  cutindex++;
557 
558  if(event->IsA()==AliESDEvent::Class()){
559  AliTriggerAnalysis fTriggerAnalysis;// = new AliTriggerAnalysis;
560  fHasV0AND = fTriggerAnalysis.IsOfflineTriggerFired((AliESDEvent*)event, AliTriggerAnalysis::kV0AND);
561  if(fHasV0AND&&hTriggerClass)hTriggerClass->Fill(32);
562  }
563  // cout << "event number " << ((AliESDEvent*)event)->GetEventNumberInFile() << " entered"<< endl;
564 
565 
566  // Number of Contributors Cut
567  if (fEnableVertexCut){
568  if(GetNumberOfContributorsVtx(event)<=0) {
569  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
570  fEventQuality = 5;
571  return kFALSE;
572  }
573  }
574  cutindex++;
575 
576  // Z Vertex Position Cut
577  if (fEnableVertexCut){
578  if(!VertexZCut(event)){
579  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
580  fEventQuality = 4;
581  return kFALSE;
582  }
583  }
584  cutindex++;
585 
586  // SPD clusters vs tracklets to check for pileup/background
587  Int_t nClustersLayer0 = event->GetNumberOfITSClusters(0);
588  Int_t nClustersLayer1 = event->GetNumberOfITSClusters(1);
589  Int_t nTracklets = event->GetMultiplicity()->GetNumberOfTracklets();
590  if(hSPDClusterTrackletBackgroundBefore) hSPDClusterTrackletBackgroundBefore->Fill(nTracklets, (nClustersLayer0 + nClustersLayer1));
591 
592 
593  Double_t distZMax = 0;
594  if(event->IsA()==AliESDEvent::Class()){
595  Int_t nPileVert = ((AliESDEvent*)event)->GetNumberOfPileupVerticesSPD();
596  if (hNPileupVertices) hNPileupVertices->Fill(nPileVert);
597  if (nPileVert > 0){
598  for(Int_t i=0; i<nPileVert;i++){
599  const AliESDVertex* pv = ((AliESDEvent*)event)->GetPileupVertexSPD(i);
600  Int_t nc2 = pv->GetNContributors();
601  if(nc2>=3){
602  Double_t z1 = ((AliESDEvent*)event)->GetPrimaryVertexSPD()->GetZ();
603  Double_t z2 = pv->GetZ();
604  Double_t distZ = z2-z1;
605  if (TMath::Abs(distZMax) < TMath::Abs(distZ) ){
606  distZMax = distZ;
607  }
608  }
609  }
610  if (hPileupVertexToPrimZ) hPileupVertexToPrimZ->Fill(distZMax);
611  }
612  }
613 
614  // Pile Up Rejection
615  if (fIsHeavyIon == 2){
616  if(GetUseNewMultiplicityFramework()){// for Run2 pPb
617  if(fUtils->IsPileUpMV(event)){
618  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
619  fEventQuality = 6;
620  return kFALSE;
621  }
622  } else{
623  if(fUtils->IsFirstEventInChunk(event)){
624  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
625  fEventQuality = 6;
626  return kFALSE;
627  }
628  if(fRemovePileUp){
629  if(fUtils->IsPileUpEvent(event)){
630  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
632  fEventQuality = 6;
633  return kFALSE;
634  }
635  if (fUtils->IsSPDClusterVsTrackletBG(event)){
636  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
638  fEventQuality = 11;
639  return kFALSE;
640  }
641  }
642  }
643  } else if(fRemovePileUp){
644  if(event->IsPileupFromSPD(3,0.8,3.,2.,5.) ){
645  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
647  fEventQuality = 6;
648  return kFALSE;
649  }
650  if (fUtils->IsSPDClusterVsTrackletBG(event)){
651  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
653  fEventQuality = 11;
654  return kFALSE;
655  }
656  }
657  cutindex++;
658 
659  // Centrality Selection
660  if(!IsCentralitySelected(event,mcEvent)){
661  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
662  fEventQuality = 1;
663  return kFALSE;
664  }
665  cutindex++;
666 
668  if(fHistoEventCuts) fHistoEventCuts->Fill(cutindex);
669  fEventQuality = 12;
670  return kFALSE;
671  }
672  cutindex++;
673  // Fill Event Histograms
674  if(fHistoEventCuts)fHistoEventCuts->Fill(cutindex);
675  if(hCentrality)hCentrality->Fill(GetCentrality(event));
676  if(hVertexZ)hVertexZ->Fill(event->GetPrimaryVertex()->GetZ());
677  // if(hCentralityVsNumberOfPrimaryTracks)
678  // hCentralityVsNumberOfPrimaryTracks->Fill(GetCentrality(event),
679  // ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()
680  // ->GetTask(fV0ReaderName.Data()))->GetNumberOfPrimaryTracks());
681 
682  if(fIsHeavyIon == 1){
683  AliEventplane *EventPlane = event->GetEventplane();
684  fEventPlaneAngle = EventPlane->GetEventplane("V0",event,2);
685  if(hEventPlaneAngle)hEventPlaneAngle->Fill(TMath::Abs(fEventPlaneAngle));
686  }
687  if(hSPDClusterTrackletBackground) hSPDClusterTrackletBackground->Fill(nTracklets, (nClustersLayer0 + nClustersLayer1));
688 
689  fEventQuality = 0;
690  return kTRUE;
691 }
692 
696 
697  if(fCutString && fCutString->GetString().Length() == kNCuts) {
698  fCutString->SetString(GetCutNumber());
699  } else {
700  return kFALSE;
701  }
702  return kTRUE;
703 }
704 
707 
708  AliInfo("Entering loading of weights for centrality flattening");
709  TFile *w = TFile::Open(fPathWeightsFlatCent.Data());
710  if(!w){
711  AliError(Form("file for centrality flattening %s not found",fPathWeightsFlatCent.Data()));
712  return;
713  }
714 
715  if (fNameHistoNotFlatCentrality.CompareTo("") != 0 && (fDoCentralityFlat > 0)){
716  cout << "I have to find: " << fNameHistoNotFlatCentrality.Data() << endl;
717  TH1D *hCentralityNotFlattemp = (TH1D*)w->Get(fNameHistoNotFlatCentrality.Data());
718  hCentralityNotFlat = new TH1D(*hCentralityNotFlattemp);
719  if (hCentralityNotFlat) AliInfo(Form("%s has been loaded from %s", fNameHistoNotFlatCentrality.Data(),fPathWeightsFlatCent.Data() ));
720  else AliWarning(Form("%s not found in %s", fNameHistoNotFlatCentrality.Data() ,fPathWeightsFlatCent.Data()));
721  hCentralityNotFlat->SetDirectory(0);
722  }
723 
724  w->Close();
725  delete w;
726 }
727 
730 
731  AliInfo("Entering loading of weights for multiplicity weighting");
732  TFile *w = TFile::Open(fPathReweightingMult.Data());
733  if(!w){
734  AliError(Form("file for multiplicity reweighting %s not found",fPathReweightingMult.Data()));
735  return;
736  }
737 
738  if (fNameHistoReweightingMultData.CompareTo("") != 0 && (fDoMultiplicityWeighting > 0)){
739  cout << "I have to find: " << fNameHistoReweightingMultData.Data() << endl;
740  TH1D *hReweightMultDatatemp = (TH1D*)w->Get(fNameHistoReweightingMultData.Data());
741  if (hReweightMultDatatemp == NULL) AliError(Form("%s was not contained in %s", fNameHistoReweightingMultData.Data(),fPathReweightingMult.Data() ));
742  hReweightMultData = new TH1D(*hReweightMultDatatemp);
743  if (hReweightMultData) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingMultData.Data(),fPathReweightingMult.Data() ));
744  else AliWarning(Form("%s not found in %s", fNameHistoReweightingMultData.Data() ,fPathReweightingMult.Data()));
745  hReweightMultData->SetDirectory(0);
746  }
747  if (fNameHistoReweightingMultMC.CompareTo("") != 0 && (fDoMultiplicityWeighting > 0)){
748  cout << "I have to find: " << fNameHistoReweightingMultMC.Data() << endl;
749  TH1D *hReweightMultMCtemp = (TH1D*)w->Get(fNameHistoReweightingMultMC.Data());
750  if (hReweightMultMCtemp == NULL) AliError(Form("%s was not contained in %s", fNameHistoReweightingMultMC.Data(),fPathReweightingMult.Data() ));
751  hReweightMultMC = new TH1D(*hReweightMultMCtemp);
752  if (hReweightMultData) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingMultMC.Data(),fPathReweightingMult.Data() ));
753  else AliWarning(Form("%s not found in %s", fNameHistoReweightingMultMC.Data() ,fPathReweightingMult.Data()));
754  hReweightMultMC->SetDirectory(0);
755  }
756 
757  w->Close();
758  delete w;
759 }
760 
761 
764 
765  AliInfo("Entering loading of histograms for weighting");
766  TFile *f = TFile::Open(fPathTrFReweighting.Data());
767  if(!f){
768  AliError(Form("file for weighting %s not found",fPathTrFReweighting.Data()));
769  return;
770  }
771  if (fNameHistoReweightingPi0.CompareTo("") != 0 && fDoReweightHistoMCPi0 ){
772  cout << "I have to find: " << fNameHistoReweightingPi0.Data() << endl;
773  TH1D *hReweightMCHistPi0temp = (TH1D*)f->Get(fNameHistoReweightingPi0.Data());
774  hReweightMCHistPi0 = new TH1D(*hReweightMCHistPi0temp);
775  if (hReweightMCHistPi0) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingPi0.Data(),fPathTrFReweighting.Data() ));
776  else AliWarning(Form("%s not found in %s", fNameHistoReweightingPi0.Data() ,fPathTrFReweighting.Data()));
777  hReweightMCHistPi0->SetDirectory(0);
778  }
779  if (fNameFitDataPi0.CompareTo("") != 0 && fDoReweightHistoMCPi0 ){
780  cout << "I have to find: " << fNameFitDataPi0.Data() << endl;
781  TF1 *fFitDataPi0temp = (TF1*)f->Get(fNameFitDataPi0.Data());
782  fFitDataPi0 = new TF1(*fFitDataPi0temp);
783  if (fFitDataPi0) AliInfo(Form("%s has been loaded from %s", fNameFitDataPi0.Data(),fPathTrFReweighting.Data() ));
784  else AliWarning(Form("%s not found in %s",fPathTrFReweighting.Data(), fNameFitDataPi0.Data() ));
785  }
786 
787  if (fNameHistoReweightingEta.CompareTo("") != 0 && fDoReweightHistoMCEta){
788  cout << "I have to find: " << fNameHistoReweightingEta.Data() << endl;
789  TH1D *hReweightMCHistEtatemp = (TH1D*)f->Get(fNameHistoReweightingEta.Data());
790  hReweightMCHistEta = new TH1D(*hReweightMCHistEtatemp);
791  if (hReweightMCHistEta) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingEta.Data(),fPathTrFReweighting.Data() ));
792  else AliWarning(Form("%s not found in %s", fNameHistoReweightingEta.Data(),fPathTrFReweighting.Data() ));
793  hReweightMCHistEta->SetDirectory(0);
794  }
795 
796  if (fNameFitDataEta.CompareTo("") != 0 && fDoReweightHistoMCEta){
797  cout << "I have to find: " << fNameFitDataEta.Data() << endl;
798  TF1 *fFitDataEtatemp = (TF1*)f->Get(fNameFitDataEta.Data());
799  fFitDataEta = new TF1(*fFitDataEtatemp);
800  if (fFitDataEta) AliInfo(Form("%s has been loaded from %s", fNameFitDataEta.Data(),fPathTrFReweighting.Data() ));
801  else AliWarning(Form("%s not found in %s", fNameFitDataEta.Data(),fPathTrFReweighting.Data() ));
802 
803  }
804  if (fNameHistoReweightingK0s.CompareTo("") != 0 && fDoReweightHistoMCK0s){
805  cout << "I have to find: " << fNameHistoReweightingK0s.Data() << endl;
806  TH1D *hReweightMCHistK0stemp = (TH1D*)f->Get(fNameHistoReweightingK0s.Data());
807  hReweightMCHistK0s = new TH1D(*hReweightMCHistK0stemp);
808  if (hReweightMCHistK0s) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingK0s.Data(),fPathTrFReweighting.Data() ));
809  else AliWarning(Form("%s not found in %s", fNameHistoReweightingK0s.Data(),fPathTrFReweighting.Data() ));
810  hReweightMCHistK0s->SetDirectory(0);
811  }
812 
813  if (fNameFitDataK0s.CompareTo("") != 0 && fDoReweightHistoMCK0s){
814  cout << "I have to find: " << fNameFitDataK0s.Data() << endl;
815  TF1 *fFitDataK0stemp = (TF1*)f->Get(fNameFitDataK0s.Data());
816  fFitDataK0s = new TF1(*fFitDataK0stemp);
817  if (fFitDataK0s) AliInfo(Form("%s has been loaded from %s", fNameFitDataK0s.Data(),fPathTrFReweighting.Data() ));
818  else AliWarning(Form("%s not found in %s", fNameFitDataK0s.Data(),fPathTrFReweighting.Data() ));
819  }
820  f->Close();
821  delete f;
822 }
823 
824 
827  fCutStringRead = Form("%s",analysisCutSelection.Data());
828 
829  // Initialize Cuts from a given Cut string
830  if(fDoCentralityFlat > 0){
831  AliInfo("Centrality flattening was enabled");
833  }
834 
836  AliInfo("Multiplicity weighting was enabled");
838  }
840  AliInfo("Particle Weighting was enabled");
842  }
843 
844 
845  AliInfo(Form("Set Event Cut Number: %s",analysisCutSelection.Data()));
846  if(analysisCutSelection.Length()!=kNCuts) {
847  AliError(Form("Cut selection has the wrong length! size is %d, number of cuts is %d", analysisCutSelection.Length(), kNCuts));
848  return kFALSE;
849  }
850  if(!analysisCutSelection.IsAlnum()){
851  AliError("Cut selection is not alphanumeric");
852  return kFALSE;
853  }
854 
855  if (fV0ReaderName.CompareTo("") == 0){
856  fV0ReaderName = "V0ReaderV1";
857  }
858  TString analysisCutSelectionLowerCase = Form("%s",analysisCutSelection.Data());
859  analysisCutSelectionLowerCase.ToLower();
860  const char *cutSelection = analysisCutSelectionLowerCase.Data();
861  #define ASSIGNARRAY(i) fCuts[i] = ((int)cutSelection[i]>=(int)'a') ? cutSelection[i]-'a'+10 : cutSelection[i]-'0'
862  for(Int_t ii=0;ii<kNCuts;ii++){
863  ASSIGNARRAY(ii);
864  }
865 
866  // Set Individual Cuts
867  for(Int_t ii=0;ii<kNCuts;ii++){
868  if(!SetCut(cutIds(ii),fCuts[ii]))return kFALSE;
869  }
870 
872 
873  return kTRUE;
874 }
875 
879 // "HeavyIon", //0
880 // "CentralityMin", //1
881 // "CentralityMax", //2
882 // "SelectSpecialTrigger", //3
883 // "SelectSpecialSubTriggerClass", //4
884 // "RemovePileUp", //5
885 // "RejectExtraSignals", //6
886 // "VertexCut", //7
887 
888  switch (cutID) {
889  case kisHeavyIon:
890  if( SetIsHeavyIon(value)) {
891  fCuts[kisHeavyIon] = value;
892  UpdateCutString();
893  return kTRUE;
894  } else return kFALSE;
895  case kCentralityMin:
896  if( SetCentralityMin(value)) {
897  fCuts[kCentralityMin] = value;
898  UpdateCutString();
899  return kTRUE;
900  } else return kFALSE;
901  case kCentralityMax:
902  if( SetCentralityMax(value)) {
903  fCuts[kCentralityMax] = value;
904  UpdateCutString();
905  return kTRUE;
906  } else return kFALSE;
908  if( SetSelectSpecialTrigger(value)) {
910  UpdateCutString();
911  return kTRUE;
912  } else return kFALSE;
914  if( SetSelectSubTriggerClass(value)) {
915  fCuts[kSelectSubTriggerClass] = value;
916  UpdateCutString();
917  return kTRUE;
918  } else return kFALSE;
919  case kremovePileUp:
920  if( SetRemovePileUp(value)) {
921  fCuts[kremovePileUp] = value;
922  UpdateCutString();
923  return kTRUE;
924  } else return kFALSE;
925  case kExtraSignals:
926  if( SetRejectExtraSignalsCut(value)) {
927  fCuts[kExtraSignals] = value;
928  UpdateCutString();
929  return kTRUE;
930  } else return kFALSE;
931  case kVertex:
932  if( SetVertexCut(value)) {
933  fCuts[kVertex] = value;
934  UpdateCutString();
935  return kTRUE;
936  } else return kFALSE;
937 
938  case kNCuts:
939  AliError("Cut id out of range");
940  return kFALSE;
941  }
942 
943  AliError("Cut id %d not recognized");
944  return kFALSE;
945 }
946 
949  // Print out current Cut Selection
950  for(Int_t ic = 0; ic < kNCuts; ic++) {
951  printf("%-30s : %d \n", fgkCutNames[ic], fCuts[ic]);
952  }
953 }
954 
956  // Print out current Cut Selection with value
957  printf("\nEvent cutnumber \n");
958  for(Int_t ic = 0; ic < kNCuts; ic++) {
959  printf("%d",fCuts[ic]);
960  }
961  printf("\n\n");
962  if (fIsHeavyIon == 0) {
963  printf("Running in pp mode \n");
964  if (fSpecialTrigger == 0){
965  if (fSpecialSubTrigger == 0){
966  printf("\t only events triggered by V0OR will be analysed \n");
967  } else if (fSpecialSubTrigger == 1){
968  printf("\t only events where SDD was present will be analysed \n");
969  }
970  } else if (fSpecialTrigger == 1){
971  if (fSpecialSubTrigger == 0){
972  printf("\t only events triggered by V0AND will be analysed \n");
973  } else if(fSpecialSubTrigger == 1){
974  printf("\t only events where SDD was present will be analysed and triggered by VOAND\n");
975  }
976  if (fRejectTriggerOverlap) printf("\t reject trigger overlaps");
977  } else if (fSpecialTrigger > 1){
978  printf("\t only events triggered by %s %s\n", fSpecialTriggerName.Data(), fSpecialSubTriggerName.Data());
979  if (fRejectTriggerOverlap) printf("\t reject trigger overlaps\n\n");
980  }
981  if ( !(fCentralityMin == 0 && fCentralityMax == 0) && !(fCentralityMax < fCentralityMin) ){
982  printf("\t Multiplicity cut %d - %d \n", fCentralityMin, fCentralityMax);
983  }
984  } else if (fIsHeavyIon == 1){
985  printf("Running in PbPb mode \n");
986  if (fDetectorCentrality == 0){
987  printf("\t centrality selection based on V0M \n");
988  } else if (fDetectorCentrality == 1){
989  printf("\t centrality selection based on Cl1 \n");
990  }
991  if (fModCentralityClass == 0){
992  printf("\t %d - %d \n", fCentralityMin*10, fCentralityMax*10);
993  } else if ( fModCentralityClass == 1){
994  printf("\t %d - %d \n", fCentralityMin*5, fCentralityMax*5);
995  } else if ( fModCentralityClass == 2){
996  printf("\t %d - %d \n", fCentralityMin*5+45, fCentralityMax*5+45);
997  } else if (fModCentralityClass == 3){
998  printf("\t %d - %d, with Track mult in MC as data \n", fCentralityMin*10, fCentralityMax*10);
999  } else if ( fModCentralityClass == 4){
1000  printf("\t %d - %d, with Track mult in MC as data \n", fCentralityMin*5, fCentralityMax*5);
1001  } else if ( fModCentralityClass == 5){
1002  printf("\t %d - %d, with Track mult in MC as data \n", fCentralityMin*5+45, fCentralityMax*5+45);
1003  }
1004  if (fSpecialTrigger == 0){
1005  printf("\t only events triggered by kMB, kCentral, kSemiCentral will be analysed \n");
1006  } else if (fSpecialTrigger > 1){
1007  printf("\t only events triggered by %s %s\n", fSpecialTriggerName.Data(), fSpecialSubTriggerName.Data());
1008  printf("\n\t SpecialTrigger is: %s\n", fSpecialTriggerName.Data());
1009  printf("\t SpecialSubTrigger is: %s\n", fSpecialSubTriggerName.Data());
1010  if (fRejectTriggerOverlap) printf("\t reject trigger overlaps\n\n");
1011  }
1012  } else if (fIsHeavyIon == 2){
1013  printf("Running in pPb mode \n");
1014  if (fDetectorCentrality == 0){
1015  printf("\t centrality selection based on V0A \n");
1016  } else if (fDetectorCentrality == 1){
1017  printf("\t centrality selection based on Cl1 \n");
1018  }
1019  if (fModCentralityClass == 0){
1020  printf("\t %d - %d \n", fCentralityMin*10, fCentralityMax*10);
1021  }
1022  if (fSpecialTrigger == 0){
1023  printf("\t only events triggered by kINT7 will be analysed \n");
1024  } else if (fSpecialTrigger > 1){
1025  printf("\t only events triggered by %s %s\n", fSpecialTriggerName.Data(), fSpecialSubTriggerName.Data());
1026  if (fRejectTriggerOverlap) printf("\t reject trigger overlaps\n\n");
1027  }
1028  }
1029  if (fEnableVertexCut) printf("\t Vertex cut with |Z_{vtx}| <%2.2f \n",fMaxVertexZ);
1030  else printf("\t No vertex cut \n");
1031 
1032  if (fRemovePileUp ==1 ) {
1033  printf("\t Doing pile up removal \n");
1034  if (fDoPileUpRejectV0MTPCout ==1 ){
1035  printf("\t Doing extra pile up removal V0M vs TPCout \n");
1036  }
1037  }
1038 
1039  printf("MC event cuts: \n");
1040  if (fRejectExtraSignals == 0) printf("\t no rejection was applied \n");
1041  else if (fRejectExtraSignals == 1) printf("\t only MB header will be inspected \n");
1042  else if (fRejectExtraSignals > 1) printf("\t special header have been selected \n");
1043  printf("\t maximum factor between jet and pt hard = %2.2f \n", fMaxFacPtHard);
1044 }
1045 
1048 { // Set Cut
1049  switch(isHeavyIon){
1050  case 0:
1051  fIsHeavyIon=0;
1052  break;
1053  case 1:
1054  fIsHeavyIon=1;
1056  break;
1057  case 2:
1058  fIsHeavyIon=1;
1060  break;
1061  case 3: //allows to select centrality 0-45% in steps of 5% for V0 Multiplicity
1062  fIsHeavyIon=1;
1065  break;
1066  case 4: //allows to select centrality 45-90% in steps of 5% for V0 Multiplicity
1067  fIsHeavyIon=1;
1070  break;
1071  case 5: //strict cut on v0 tracks for MC
1072  fIsHeavyIon=1;
1075  break;
1076  case 6: //allows to select centrality 0-45% in steps of 5% for track mult
1077  //strict cut on v0 tracks for MC
1078  fIsHeavyIon=1;
1081  break;
1082  case 7: //allows to select centrality 45-90% in steps of 5% for V0 Multiplicity
1083  //strict cut on v0 tracks for MC
1084  fIsHeavyIon=1;
1087  break;
1088  case 8:
1089  fIsHeavyIon=2;
1091  break;
1092  case 9:
1093  fIsHeavyIon=2;
1095  break;
1096  default:
1097  AliError(Form("SetHeavyIon not defined %d",isHeavyIon));
1098  return kFALSE;
1099  }
1100  return kTRUE;
1101 }
1102 
1103 //___________________________________________________________________
1105 {
1106  // Set Cut
1107  if(minCentrality<0||minCentrality>9){
1108  AliError(Form("minCentrality not defined %d",minCentrality));
1109  return kFALSE;
1110  }
1111 
1112  fCentralityMin=minCentrality;
1113  return kTRUE;
1114 }
1115 
1116 //___________________________________________________________________
1118 {
1119  // Set Cut
1120  if(maxCentrality<0||maxCentrality>9){
1121  AliError(Form("maxCentrality not defined %d",maxCentrality));
1122  return kFALSE;
1123  }
1124  fCentralityMax=maxCentrality;
1125  return kTRUE;
1126 }
1127 
1130 {
1131  // Set Cut
1132  switch(selectSpecialTrigger){
1133  case 0:
1134  fSpecialTrigger=0; // V0OR
1135  break;
1136  case 1:
1137  fSpecialTrigger=1; // V0AND
1138  break;
1139 // case 2:
1140 // fSpecialTrigger=2; //
1141 // break;
1142  case 3:
1143  fSpecialTrigger=3; //specific centrality trigger selection
1144  fSpecialTriggerName="AliVEvent::kCentral/kSemiCentral/kMB";
1145  break;
1146  case 4:
1147  fSpecialTrigger=4; // trigger alias kTRD
1148  fOfflineTriggerMask=AliVEvent::kTRD;
1149  fTriggerSelectedManually = kTRUE;
1150  fSpecialTriggerName="AliVEvent::kTRD";
1151  break;
1152  case 5:
1153  fSpecialTrigger=5; // trigger alias kEMC
1154  fOfflineTriggerMask=AliVEvent::kEMC7 | AliVEvent::kEMC8 | AliVEvent::kEMC1 ;
1155  fTriggerSelectedManually = kTRUE;
1157  SETBIT(fTriggersEMCALSelected, kL0);
1158  fSpecialTriggerName="AliVEvent::kEMC7/kEMC8/kEMC1";
1159  break;
1160  case 6:
1161  fSpecialTrigger=6; // trigger alias kPHI
1162  fOfflineTriggerMask=AliVEvent::kPHI7 | AliVEvent::kPHI1 | AliVEvent::kPHI8 | AliVEvent::kPHOSPb;
1163  fTriggerSelectedManually = kTRUE;
1164  fSpecialTriggerName="AliVEvent::kPHI7/kPHI1/kPHI8/kPHOSPb";
1165  break;
1166  case 7:
1167  fSpecialTrigger=7; // trigger alias kHighMult
1168  fOfflineTriggerMask=AliVEvent::kHighMult;
1169  fTriggerSelectedManually = kTRUE;
1170  fSpecialTriggerName="AliVEvent::kHighMult";
1171  break;
1172  case 8:
1173  fSpecialTrigger=8; // trigger alias kEMCEGA
1174  fOfflineTriggerMask=AliVEvent::kEMCEGA;
1175  fTriggerSelectedManually = kTRUE;
1177  SETBIT(fTriggersEMCALSelected, kG2);
1178  fSpecialTriggerName="AliVEvent::kEMCEGA";
1179  break;
1180  case 9:
1181  fSpecialTrigger=9; // trigger alias kEMCEJE
1182  fOfflineTriggerMask=AliVEvent::kEMCEJE;
1183  fTriggerSelectedManually = kTRUE;
1185  SETBIT(fTriggersEMCALSelected, kJ2);
1186  fSpecialTriggerName="AliVEvent::kEMCEJE";
1187  break;
1188  default:
1189  AliError("Warning: Special Trigger Not known");
1190  return 0;
1191  }
1192  return 1;
1193 }
1194 
1197 {
1198  // Set Cut
1199  if (fSpecialTrigger == 0){ //OR
1200  switch(selectSpecialSubTriggerClass){
1201  case 0://with VZERO
1202  fSpecialTrigger=0;
1204 // AliInfo("Info: Nothing to be done");
1205  break;
1206  case 3: //V0OR with SDD requested (will only work with LHC11a dataset)
1208 // cout << "V0OR with SDD requested" << endl;
1209  break;
1210  default:
1211  AliError("Warning: Special Subtrigger Class Not known");
1212  return 0;
1213  }
1214  } else if (fSpecialTrigger == 1){ //AND with different detectors
1215  switch(selectSpecialSubTriggerClass){
1216  case 0: //with VZERO general implementation of V0AND (periods LHC11c onwards)
1217  fSpecialTrigger=0;
1219  fOfflineTriggerMask=AliVEvent::kINT7;
1220  fTriggerSelectedManually = kTRUE;
1221  fSpecialTriggerName="AliVEvent::kINT7";
1222  break;
1223  case 1: //with TZERO
1224  fSpecialTrigger=0;
1226  fOfflineTriggerMask=AliVEvent::kINT8;
1227  fTriggerSelectedManually = kTRUE;
1228  fSpecialTriggerName="AliVEvent::kINT8";
1229  break;
1230  case 2: //with VZERO (will only work with LHC11a dataset)
1231  fSpecialTrigger=1;
1233 // AliInfo("Info: Nothing to be done");
1234  break;
1235  case 3: //V0AND with SDD requested (will only work with LHC11a dataset)
1236  fSpecialTrigger=1;
1238  break;
1239  default:
1240  AliError("Warning: Special Subtrigger Class Not known");
1241  return 0;
1242  }
1243  } else if (fSpecialTrigger == 3){ // Selecting kCentral and kSemiCentral from trigger classes, not aliases
1244  switch(selectSpecialSubTriggerClass){
1245  case 0: // all together
1248 // AliInfo("Info: Nothing to be done");
1249  break;
1250  case 1: // kCentral - no vertex restriction
1253  fSpecialSubTriggerName="CVHN";
1254 // cout << "kCentralOpen" << endl;
1255  break;
1256  case 2: // kCentral - T00 +- 10 cm
1259  fSpecialSubTriggerName="CCENT";
1260 // cout << "kCentralVertex" << endl;
1261  break;
1262  case 3: // kCentral - both
1265  fSpecialSubTriggerName="CVHN|CCENT|CSEMI|CVLN";
1266 // cout << "kCentral both" << endl;
1267  break;
1268  case 4: // kSemiCentral - no vertex restriction
1271  fSpecialSubTriggerName="CVLN";
1272 // cout << "kSemiCentralOpen" << endl;
1273  break;
1274  case 5: // kSemiCentral - T00 +- 10 cm
1277  fSpecialSubTriggerName="CSEMI";
1278 // cout << "kSemiCentralVertex" << endl;
1279  break;
1280  case 6: // kSemiCentral - both
1283  fSpecialSubTriggerName="CSEMI%CVLN";
1284 // cout << "kSemiCentral both" << endl;
1285  break;
1286  case 7: // kMB
1289  fSpecialSubTriggerName="CPBI1_|CPBI1-";
1290 // cout << "kMB 1" << endl;
1291  break;
1292  case 8: // kMB
1295  fSpecialSubTriggerName="CPBI2_|CPBI2-";
1296 // cout << "kMB 2" << endl;
1297  break;
1298  case 9: // kMB
1301  fSpecialSubTriggerName="CPBI2_@CPBI2-@CPBI2_@CPBI2-";
1302 // cout << "kMB both" << endl;
1303  break;
1304  default:
1305  AliError("Warning: Special Subtrigger Class Not known");
1306  return 0;
1307  }
1308  } else if (fSpecialTrigger == 4){ // Subdivision of TRD trigger classes
1309  switch(selectSpecialSubTriggerClass){
1310  case 0: // all together
1313 // AliInfo("Info: Nothing to be done");
1314  break;
1315  case 1: // 7WUHSH - V0AND with single electron in TRD & EMCAL
1318  fSpecialSubTriggerName="7WUHEE";
1319  break;
1320  case 2: // 8WUHSH - T0AND with single electron in TRD & EMCAL
1323  fSpecialSubTriggerName="8WUHEE";
1324  break;
1325  case 3: // 7WUHSE - V0AND with single high pt electron in TRD
1328  fSpecialSubTriggerName="7WUHSE";
1329  break;
1330  case 4: // 8WUHSE - T0AND with single high pt electron in TRD
1333  fSpecialSubTriggerName="8WUHSE";
1334  break;
1335  case 5: // 7WUHJE - V0AND with jet in TRD
1338  fSpecialSubTriggerName="7WUHJT";
1339  break;
1340  case 6: // 8WUHJE - T0AND with jet in TRD
1343  fSpecialSubTriggerName="8WUHJT";
1344  break;
1345  case 7: // 7WUHQU - V0AND with dielectron pair in TRD
1348  fSpecialSubTriggerName="7WUHQU";
1349  break;
1350  case 8: // 8WUHQU - T0AND with dielectron pair in TRD
1353  fSpecialSubTriggerName="8WUHQU";
1354  break;
1355  default:
1356  AliError("Warning: Special Subtrigger Class Not known");
1357  return 0;
1358  }
1359  } else if (fSpecialTrigger == 5){ // Subdivision of kEMC trigger classes
1360  switch(selectSpecialSubTriggerClass){
1361  case 0: // all together
1364 // AliInfo("Info: Nothing to be done");
1365  break;
1366  case 1: // CEMC1 - V0OR and EMCAL fired
1367  fOfflineTriggerMask=AliVEvent::kEMC1;
1368  fSpecialTriggerName="AliVEvent::kEMC1";
1371  fSpecialSubTriggerName="CEMC1";
1372  break;
1373  case 2: // CEMC7 - V0AND and EMCAL fired
1375  fOfflineTriggerMask=AliVEvent::kEMC7;
1376  fSpecialTriggerName="AliVEvent::kEMC7";
1378  fSpecialSubTriggerName="CEMC7";
1379  break;
1380  case 3: // CEMC8 - T0OR and EMCAL fired
1381  fOfflineTriggerMask=AliVEvent::kEMC8;
1382  fSpecialTriggerName="AliVEvent::kEMC8";
1385  fSpecialSubTriggerName="CEMC8";
1386  break;
1387  case 4: // CDMC1 - V0OR and DCAL fired
1388  fOfflineTriggerMask=AliVEvent::kEMC1;
1389  fSpecialTriggerName="AliVEvent::kEMC1";
1392  fSpecialSubTriggerName="CDMC1";
1393  break;
1394  case 5: // CDMC7 - V0AND and DCAL fired
1396  fOfflineTriggerMask=AliVEvent::kEMC7;
1397  fSpecialTriggerName="AliVEvent::kEMC7";
1399  fSpecialSubTriggerName="CDMC7";
1400  break;
1401  case 6: // CDMC8 - T0OR and DCAL fired
1402  fOfflineTriggerMask=AliVEvent::kEMC8;
1403  fSpecialTriggerName="AliVEvent::kEMC8";
1406  fSpecialSubTriggerName="CDMC8";
1407  break;
1408  default:
1409  AliError("Warning: Special Subtrigger Class Not known");
1410  return 0;
1411  }
1412  }else if (fSpecialTrigger == 6){ // Subdivision of kPHI trigger classes
1413  switch(selectSpecialSubTriggerClass){
1414  case 0: // all together
1417 // AliInfo("Info: Nothing to be done");
1418  break;
1419  case 1: // CEMC1 - V0OR and EMCAL fired
1420  fOfflineTriggerMask=AliVEvent::kPHI1;
1421  fSpecialTriggerName="AliVEvent::kPHI1";
1424  fSpecialSubTriggerName="CPHI1";
1425  break;
1426  case 2: // CEMC7 - V0AND and EMCAL fired
1428  fOfflineTriggerMask=AliVEvent::kPHI7;
1429  fSpecialTriggerName="AliVEvent::kPHI7";
1431  fSpecialSubTriggerName="CPHI7";
1432  break;
1433  case 3: // CEMC8 - T0OR and EMCAL fired
1434  fOfflineTriggerMask=AliVEvent::kPHI8;
1435  fSpecialTriggerName="AliVEvent::kPHI8";
1438  fSpecialSubTriggerName="CPHI8";
1439  break;
1440  default:
1441  AliError("Warning: Special Subtrigger Class Not known");
1442  return 0;
1443  }
1444  } else if (fSpecialTrigger == 7){ // Subdivision of kHighMult trigger classes
1445  switch(selectSpecialSubTriggerClass){
1446  case 0: // all together
1449 // AliInfo("Info: Nothing to be done");
1450  break;
1451  case 1: // CSHM1 - V0OR and high mult fired
1454  fSpecialSubTriggerName="CSHM1";
1455  break;
1456  case 2: // CSHM7 - V0AND and high mult fired
1459  fSpecialSubTriggerName="CSHM7";
1460  break;
1461  case 3: // CSHM8 - T0OR and high mult fired
1464  fSpecialSubTriggerName="CSHM8";
1465  break;
1466  case 4: // V0 high mult trigger
1468  fOfflineTriggerMask=AliVEvent::kAny;
1469  fSpecialTriggerName="V0Mult";
1471  fSpecialSubTriggerName="CVHMV0M-B-";
1472  break;
1473  case 5: // SPD high mult trigger
1475  fOfflineTriggerMask=AliVEvent::kAny;
1476  fSpecialTriggerName="SPMult";
1478  fSpecialSubTriggerName="CVHMSH2-B-";
1479  break;
1480  case 6: // V0 high mult trigger with pileup condition on
1482  fOfflineTriggerMask=AliVEvent::kAny;
1483  fSpecialTriggerName="V0Mult";
1485  fSpecialSubTriggerName="CVHMV0M-B-SPD2";
1486  break;
1487 
1488  default:
1489  AliError("Warning: Special Subtrigger Class Not known");
1490  return 0;
1491  }
1492  }else if (fSpecialTrigger == 8){ // Subdivision of kEMCEGA trigger classes
1493  switch(selectSpecialSubTriggerClass){
1494  case 0: // all together
1497 // AliInfo("Info: Nothing to be done");
1498  break;
1499  case 1: // 7EGA - CINT7 EGA
1502  fSpecialSubTriggerName="7EGA";
1504  SETBIT(fTriggersEMCALSelected, kG2);
1505  break;
1506  case 2: // 8EGA - CINT8 EGA
1509  fSpecialSubTriggerName="8EGA";
1511  SETBIT(fTriggersEMCALSelected, kG2);
1512  break;
1513  case 3: // 7EG1 - CINT7 EG1
1516  fSpecialSubTriggerName="7EG1";
1518  SETBIT(fTriggersEMCALSelected, kG1);
1519  break;
1520  case 4: // 8EG1 - CINT8 EG1
1523  fSpecialSubTriggerName="8EG1";
1525  SETBIT(fTriggersEMCALSelected, kG1);
1526  break;
1527  case 5: // 7EG2 - CINT7 EG2
1530  fSpecialSubTriggerName="7EG2";
1532  SETBIT(fTriggersEMCALSelected, kG2);
1533  break;
1534  case 6: // 8EG2 - CINT8 EG2
1537  fSpecialSubTriggerName="8EG2";
1539  SETBIT(fTriggersEMCALSelected, kG2);
1540  break;
1541  case 7: // 7DGA - CINT7 DGA
1544  fSpecialSubTriggerName="7DGA";
1546  SETBIT(fTriggersEMCALSelected, kG2);
1547  break;
1548  case 8: // 8DGA - CINT8 DGA
1551  fSpecialSubTriggerName="8DGA";
1553  SETBIT(fTriggersEMCALSelected, kG2);
1554  break;
1555  case 9: // 7DG1 - CINT7 DG1
1558  fSpecialSubTriggerName="7DG1";
1560  SETBIT(fTriggersEMCALSelected, kG1);
1561  break;
1562  case 10: // 8DG1 - CINT8 DG1
1565  fSpecialSubTriggerName="8DG1";
1567  SETBIT(fTriggersEMCALSelected, kG1);
1568  break;
1569  case 11: // 7DG2 - CINT7 DG2
1572  fSpecialSubTriggerName="7DG2";
1574  SETBIT(fTriggersEMCALSelected, kG2);
1575  break;
1576  case 12: // 8DG2 - CINT8 DG2
1579  fSpecialSubTriggerName="8DG2";
1581  SETBIT(fTriggersEMCALSelected, kG2);
1582  break;
1583 
1584  default:
1585  AliError("Warning: Special Subtrigger Class Not known");
1586  return 0;
1587  }
1588  } else if (fSpecialTrigger == 9){ // Subdivision of kEMCEGA trigger classes
1589  switch(selectSpecialSubTriggerClass){
1590  case 0: // all together
1593 // AliInfo("Info: Nothing to be done");
1594  break;
1595  case 1: // 7EJE - CINT7 EJE
1598  fSpecialSubTriggerName="7EJE";
1600  SETBIT(fTriggersEMCALSelected, kJ2);
1601  break;
1602  case 2: // 8EJE - CINT8 EJE
1605  fSpecialSubTriggerName="8EJE";
1607  SETBIT(fTriggersEMCALSelected, kJ2);
1608  break;
1609  case 3: // 7EJ1 - CINT7 EJ1
1612  fSpecialSubTriggerName="7EJ1";
1614  SETBIT(fTriggersEMCALSelected, kJ1);
1615  break;
1616  case 4: // 8EJ1 - CINT8 EJ1
1619  fSpecialSubTriggerName="8EJ1";
1621  SETBIT(fTriggersEMCALSelected, kJ1);
1622  break;
1623  case 5: // 7EJ2 - CINT7 EJ2
1626  fSpecialSubTriggerName="7EJ2";
1628  SETBIT(fTriggersEMCALSelected, kJ2);
1629  break;
1630  case 6: // 8EJ2 - CINT8 EJ2
1633  fSpecialSubTriggerName="8EJ2";
1635  SETBIT(fTriggersEMCALSelected, kJ2);
1636  break;
1637  case 7: // 7DJ1 - CINT7 DJ1
1640  fSpecialSubTriggerName="7DJ1";
1642  SETBIT(fTriggersEMCALSelected, kJ1);
1643  break;
1644  case 8: // 8DJ1 - CINT8 DJ1
1647  fSpecialSubTriggerName="8DJ1";
1649  SETBIT(fTriggersEMCALSelected, kJ1);
1650  break;
1651  case 9: // 7DJ2 - CINT7 DJ2
1654  fSpecialSubTriggerName="7DJ2";
1656  SETBIT(fTriggersEMCALSelected, kJ2);
1657  break;
1658  case 10: // 8DJ2 - CINT8 DJ2
1661  fSpecialSubTriggerName="8DJ2";
1663  SETBIT(fTriggersEMCALSelected, kJ2);
1664  break;
1665 
1666  default:
1667  AliError("Warning: Special Subtrigger Class Not known");
1668  return 0;
1669  }
1670  }
1671  return 1;
1672 }
1673 
1674 //________________________________________________________________________
1676 {
1677  // Set Cut
1678  fMultiplicityMethod=multiplicityMethod;
1679 
1680  // 0 Photon Multiplicity
1681  // 1 TPC Track multiplicity
1682  // 2 V0 Mult
1683  // 3 SPD Mult
1684 
1685  return kTRUE;
1686 }
1687 
1688 //________________________________________________________________________
1690 {// Set Cut
1691  switch(removePileUp){
1692  case 0:
1693  fRemovePileUp = kFALSE;
1694  break;
1695  case 1:
1696  fRemovePileUp = kTRUE;
1697  break;
1698  case 2:
1699  fRemovePileUp = kTRUE;
1702  break;
1703  case 3:
1704  fRemovePileUp = kTRUE;
1707  break;
1708  case 4:
1709  fRemovePileUp = kTRUE;
1712  break;
1713  case 5:
1714  fRemovePileUp = kTRUE;
1717  break;
1718  case 6:
1719  fRemovePileUp = kTRUE;
1720  fDoPileUpRejectV0MTPCout = kTRUE;
1721  fFPileUpRejectV0MTPCout = new TF1("fFPileUpRejectV0MTPCout","[0] + [1]*x",0.,10000.);
1722  fFPileUpRejectV0MTPCout->SetParameter(0,0.);
1723  fFPileUpRejectV0MTPCout->SetParameter(1,0.);
1724  if (fIsHeavyIon==1){
1725  if(fPeriodEnum == kLHC15o){
1726  fFPileUpRejectV0MTPCout->SetParameter(0,-2500.);
1727  fFPileUpRejectV0MTPCout->SetParameter(1,5.0);
1728  break;
1729  }else{
1730  fFPileUpRejectV0MTPCout->SetParameter(0,-1500.);
1731  fFPileUpRejectV0MTPCout->SetParameter(1,3.0);
1732  break;
1733  }
1734  } else if(fIsHeavyIon == 2){
1735  fFPileUpRejectV0MTPCout->SetParameter(0,-200.);
1736  fFPileUpRejectV0MTPCout->SetParameter(1,2.0);
1737  break;
1738  }else{
1739  fFPileUpRejectV0MTPCout->SetParameter(0,-300.);
1740  fFPileUpRejectV0MTPCout->SetParameter(1,4.0);
1741  break;
1742  }
1743  break;
1744  case 7:
1745  fRemovePileUp = kTRUE;
1746  fDoPileUpRejectV0MTPCout = kTRUE;
1747  fFPileUpRejectV0MTPCout = new TF1("fFPileUpRejectV0MTPCout","[0] + [1]*x",0.,10000.);
1748  fFPileUpRejectV0MTPCout->SetParameter(0,0.);
1749  fFPileUpRejectV0MTPCout->SetParameter(1,0.);
1750  if (fIsHeavyIon==1){
1751  if(fPeriodEnum == kLHC15o){
1752  fFPileUpRejectV0MTPCout->SetParameter(0,-2500.);
1753  fFPileUpRejectV0MTPCout->SetParameter(1,5.0);
1754  break;
1755  }else{
1756  fFPileUpRejectV0MTPCout->SetParameter(0,-2500.);
1757  fFPileUpRejectV0MTPCout->SetParameter(1,3.0);
1758  break;
1759  }
1760  } else if(fIsHeavyIon == 2){
1761  fFPileUpRejectV0MTPCout->SetParameter(0,-300.);
1762  fFPileUpRejectV0MTPCout->SetParameter(1,1.5);
1763  break;
1764  }else{
1765  fFPileUpRejectV0MTPCout->SetParameter(0,-300.);
1766  fFPileUpRejectV0MTPCout->SetParameter(1,3.0);
1767  break;
1768  }
1769  break;
1770  case 8:
1771  fRemovePileUp = kTRUE;
1772  fDoPileUpRejectV0MTPCout = kTRUE;
1773  fFPileUpRejectV0MTPCout = new TF1("fFPileUpRejectV0MTPCout","[0] + [1]*x",0.,10000.);
1774  fFPileUpRejectV0MTPCout->SetParameter(0,0.);
1775  fFPileUpRejectV0MTPCout->SetParameter(1,0.);
1776  if (fIsHeavyIon==1){
1777  if(fPeriodEnum == kLHC15o){
1778  fFPileUpRejectV0MTPCout->SetParameter(0,-2500.);
1779  fFPileUpRejectV0MTPCout->SetParameter(1,5.0);
1780  break;
1781  }else{
1782  fFPileUpRejectV0MTPCout->SetParameter(0,-1500.);
1783  fFPileUpRejectV0MTPCout->SetParameter(1,3.0);
1784  break;
1785  }
1786  } else if(fIsHeavyIon == 2){
1787  fFPileUpRejectV0MTPCout->SetParameter(0,-200.);
1788  fFPileUpRejectV0MTPCout->SetParameter(1,1.5);
1789  break;
1790  }else{
1791  fFPileUpRejectV0MTPCout->SetParameter(0,-300.);
1792  fFPileUpRejectV0MTPCout->SetParameter(1,4.0);
1793  break;
1794  }
1795  break;
1796  case 9:
1797  fRemovePileUp = kTRUE;
1800  fDoPileUpRejectV0MTPCout = kTRUE;
1801  fFPileUpRejectV0MTPCout = new TF1("fFPileUpRejectV0MTPCout","[0] + [1]*x",0.,10000.);
1802  fFPileUpRejectV0MTPCout->SetParameter(0,0.);
1803  fFPileUpRejectV0MTPCout->SetParameter(1,0.);
1804  if (fIsHeavyIon==1){
1805  if(fPeriodEnum == kLHC15o){
1806  fFPileUpRejectV0MTPCout->SetParameter(0,-2500.);
1807  fFPileUpRejectV0MTPCout->SetParameter(1,5.0);
1808  break;
1809  }else{
1810  fFPileUpRejectV0MTPCout->SetParameter(0,-1500.);
1811  fFPileUpRejectV0MTPCout->SetParameter(1,3.0);
1812  break;
1813  }
1814  } else if(fIsHeavyIon == 2){
1815  fFPileUpRejectV0MTPCout->SetParameter(0,-200.);
1816  fFPileUpRejectV0MTPCout->SetParameter(1,1.5);
1817  break;
1818  }else{
1819  fFPileUpRejectV0MTPCout->SetParameter(0,-300.);
1820  fFPileUpRejectV0MTPCout->SetParameter(1,4.0);
1821  break;
1822  }
1823  break;
1824  default:
1825  AliError("RemovePileUpCut not defined");
1826  return kFALSE;
1827  }
1828  return kTRUE;
1829 }
1830 
1831 //________________________________________________________________________
1833 
1834  switch(extraSignal){
1835  case 0:
1836  fRejectExtraSignals = 0;
1837  break; // No Rejection
1838  case 1:
1839  fRejectExtraSignals = 1;
1840  break; // MinBias Header
1841  case 2:
1842  fRejectExtraSignals = 2;
1843  break; // User String Array
1844  case 3:
1845  fRejectExtraSignals = 3;
1846  break; // Rejection for Gamma Correction only
1847  default:
1848  AliError(Form("Extra Signal Rejection not defined %d",extraSignal));
1849  return kFALSE;
1850  }
1851  return kTRUE;
1852 }
1853 
1854 //________________________________________________________________________
1856 
1857  switch(vertexCut){
1858  case 0: // no Vertex required // NOT fully working yet
1859  fEnableVertexCut = kFALSE;
1860  fMaxVertexZ = 1000;
1861  break;
1862  case 1: // vertex within +-15 cm
1863  fEnableVertexCut = kTRUE;
1864  fMaxVertexZ = 15;
1865  break;
1866  case 2: // vertex within +-12.5 cm
1867  fEnableVertexCut = kTRUE;
1868  fMaxVertexZ = 12.5;
1869  break;
1870  case 3: // vertex within +-10 cm
1871  fEnableVertexCut = kTRUE;
1872  fMaxVertexZ = 10.0;
1873  break;
1874  case 4: // vertex within +-7.5 cm
1875  fEnableVertexCut = kTRUE;
1876  fMaxVertexZ = 7.5;
1877  break;
1878  case 5: // vertex within +-5 cm
1879  fEnableVertexCut = kTRUE;
1880  fMaxVertexZ = 5.;
1881  break;
1882  default:
1883  AliError(Form("Vertex Cut not defined %d",vertexCut));
1884  return kFALSE;
1885  }
1886  return kTRUE;
1887 }
1888 
1889 //-------------------------------------------------------------
1891  if (fPeriodEnum == kLHC15n || // pp 5TeV
1892  fPeriodEnum == kLHC15o || // PbPb 5TeV
1893  fPeriodEnum == kLHC15k1a1 || fPeriodEnum == kLHC15k1a2 || fPeriodEnum == kLHC15k1a3 || fPeriodEnum == kLHC16j7 || // MC PbPb 5TeV LowIR
1894  fPeriodEnum == kLHC16h4 || // MC PbPb 5TeV added signals
1895  fPeriodEnum == kLHC16g1 || fPeriodEnum == kLHC16g1a || fPeriodEnum == kLHC16g1b || fPeriodEnum == kLHC16g1c || // MC PbPb 5TeV general purpose
1897  fPeriodEnum == kLHC16h2a || fPeriodEnum == kLHC16h2b || fPeriodEnum == kLHC16h2c || // MC PbPb 5TeV jet-jet
1899  fPeriodEnum == kLHC15g3a3 || fPeriodEnum == kLHC15g3c3 || // MC pp 13TeV
1900  fPeriodEnum == kLHC16q || fPeriodEnum == kLHC16t || // pPb 5TeV LHC16qt
1901  fPeriodEnum == kLHC16r || fPeriodEnum == kLHC16s || // pPb 8TeV LHC16rs
1905  fPeriodEnum == kLHC17g8a_fast || fPeriodEnum == kLHC17g8a_cent_woSDD || // MC pPb 5TeV LHC16qt
1906  fPeriodEnum == kLHC17g8b || fPeriodEnum == kLHC17g8c || // MC pPb 8TeV LHC16sr
1907  fPeriodEnum == kLHC17n || // Xe-Xe LHC17n
1908  fPeriodEnum == kLHC17j7 // MC Xe-Xe LHC17n
1909  ){
1910  return kTRUE;
1911  } else {
1912  return kFALSE;
1913  }
1914 }
1915 
1916 //-------------------------------------------------------------
1918 { // Get Event Centrality
1919 
1920  AliESDEvent *esdEvent=dynamic_cast<AliESDEvent*>(event);
1921  if(esdEvent){
1923  AliMultSelection *MultSelection = (AliMultSelection*)event->FindListObject("MultSelection");
1924  if(!MultSelection){
1925  AliWarning ("AliMultSelection object not found !");
1926  return -1;
1927  }else{
1928  if(fDetectorCentrality==0){
1929  if(fIsHeavyIon==2) return MultSelection->GetMultiplicityPercentile("V0A");// default for pPb
1930  else return MultSelection->GetMultiplicityPercentile("V0M");// default
1931  }else if(fDetectorCentrality==1) return MultSelection->GetMultiplicityPercentile("CL1",kTRUE);
1932  }
1933  }else{
1934  AliCentrality *fESDCentrality = (AliCentrality*)esdEvent->GetCentrality();
1935  if(fDetectorCentrality==0){
1936  if(fIsHeavyIon==2) return fESDCentrality->GetCentralityPercentile("V0A"); // default for pPb
1937  else return fESDCentrality->GetCentralityPercentile("V0M"); // default
1938  }else if(fDetectorCentrality==1) return fESDCentrality->GetCentralityPercentile("CL1");
1939  }
1940  }
1941 
1942  AliAODEvent *aodEvent=dynamic_cast<AliAODEvent*>(event);
1943  if(aodEvent){
1945  AliMultSelection *MultSelection = (AliMultSelection*)aodEvent->FindListObject("MultSelection");
1946  if(!MultSelection){
1947  AliWarning ("AliMultSelection object not found !");
1948  return -1;
1949  } else{
1950  if(fDetectorCentrality==0){
1951  if(fIsHeavyIon==2) return MultSelection->GetMultiplicityPercentile("V0A");// default for pPb
1952  else return MultSelection->GetMultiplicityPercentile("V0M",kTRUE);
1953  }else if(fDetectorCentrality==1) return MultSelection->GetMultiplicityPercentile("CL1",kTRUE);
1954  }
1955  }else{
1956  if(aodEvent->GetHeader()){return ((AliVAODHeader*)aodEvent->GetHeader())->GetCentrality();}
1957  }
1958  }
1959 
1960  return -1;
1961 }
1962 
1963 //_____________________________________________________________________________________
1964 Bool_t AliConvEventCuts::IsCentralitySelected(AliVEvent *event, AliMCEvent *mcEvent)
1965 {
1966  // Centrality Selection
1967  if(!fIsHeavyIon){
1968  if ((fCentralityMin == 0 && fCentralityMax == 0) || (fCentralityMin > fCentralityMax) ){
1969  return kTRUE;
1970  } else {
1971  Int_t primaryTracksPP[9] = { 0, 2, 5, 10, 15,
1972  30, 50, 100, 1000
1973  };
1974  Int_t nprimaryTracks = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetNumberOfPrimaryTracks();
1975  if ( nprimaryTracks >= primaryTracksPP[fCentralityMin] && nprimaryTracks < primaryTracksPP[fCentralityMax]){
1976  return kTRUE;
1977  } else {
1978  return kFALSE;
1979  }
1980  return kFALSE;
1981  }
1982  }
1983  if(fCentralityMin == fCentralityMax ) return kTRUE;//0-100%
1984  else if ( fCentralityMax==0) fCentralityMax=10; //CentralityRange = fCentralityMin-10*multfactor
1986  if(centrality<0 && !mcEvent)return kFALSE;
1987 
1988  Int_t centralityC=0;
1989  if (fModCentralityClass == 0){
1990  centralityC= Int_t(centrality/10);
1991  if(centralityC >= fCentralityMin && centralityC < fCentralityMax)
1992  return kTRUE;
1993  else return kFALSE;
1994  }
1995  else if (fModCentralityClass ==1){
1996  centralityC= Int_t(centrality);
1997  if(centralityC >= fCentralityMin*5 && centralityC < fCentralityMax*5){
1998  return kTRUE;
1999  } else return kFALSE;
2000  }
2001  else if (fModCentralityClass ==2){
2002  centralityC= Int_t(centrality);
2003  if(centralityC >= ((fCentralityMin*5)+45) && centralityC < ((fCentralityMax*5)+45))
2004  return kTRUE;
2005  else return kFALSE;
2006  }
2007 
2008  Int_t nprimaryTracks = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetNumberOfPrimaryTracks();
2009  Int_t PrimaryTracks10[11][2] =
2010  {
2011  {9999,9999}, // 0 //1550 changed to 9999 on 9 Dec
2012  {1210, 928}, // 10
2013  { 817, 658}, // 20
2014  { 536, 435}, // 30
2015  { 337, 276}, // 40
2016  { 197, 162}, // 50
2017  { 106, 100}, // 60
2018  { 51, 44}, // 70
2019  { 21, 18}, // 80
2020  { 0, 0}, // 90
2021  { 0, 0}// 100 // only max accessible
2022  };
2023  Int_t PrimaryTracksLHC11h10[11][2] =
2024  {
2025  {9999,9999}, // 0 //1550 changed to 9999 on 9 Dec
2026  { 985, 928}, // 10
2027  { 661, 658}, // 20
2028  { 434, 435}, // 30
2029  { 275, 276}, // 40
2030  { 173, 162}, // 50
2031  { 100, 100}, // 60
2032  { 42, 44}, // 70
2033  { 19, 18}, // 80
2034  { 0, 0}, // 90
2035  { 0, 0}// 100 // only max accessible
2036  };
2037  Int_t PrimaryTracksLHC15o10[11][2] =
2038  { // values calculated from ESD analysis so far
2039  {2700, 2700}, // 0-10% cent class max # of tracks: max value of the data distribution
2040  {1498, 1498}, // 0-10% cent class min # of tracks
2041  {1012, 1012}, // 10-20
2042  { 669, 669}, // 20-30
2043  { 423, 423}, // 30-40
2044  { 251, 251}, // 40-50
2045  { 136, 136}, // 50-60
2046  { 67, 67}, // 60-70
2047  { 28, 28}, // 70-80
2048  { 0, 0}, // 80-90% cent class min # of tracks
2049  { 0, 0} // not used
2050  };
2051  Int_t PrimaryTracks5a[11][2] =
2052  {
2053  {9999,9999}, // 0 ///1550 changed to 9999 on 9 Dec
2054  {1485,1168}, // 5
2055  {1210, 928}, // 10
2056  { 995, 795}, // 15
2057  { 817, 658}, // 20
2058  { 666, 538}, // 25
2059  { 536, 435}, // 30
2060  { 428, 350}, // 35
2061  { 337, 276}, // 40
2062  { 260, 214}, // 45
2063  { 0, 162}// 50 only max accessible
2064  };
2065  Int_t PrimaryTracksLHC11h5a[11][2] =
2066  {
2067  {9999,9999}, // 0 ///1550 changed to 9999 on 9 Dec
2068  {1166,1168}, // 5
2069  { 953, 928}, // 10
2070  { 805, 795}, // 15
2071  { 655, 658}, // 20
2072  { 535, 538}, // 25
2073  { 435, 435}, // 30
2074  { 349, 350}, // 35
2075  { 275, 276}, // 40
2076  { 214, 214}, // 45
2077  { 165, 162}// 50 only max accessible
2078  };
2079  Int_t PrimaryTracksLHC15o5a[11][2] =
2080  {
2081  { 2700, 2700}, // 0-5% cent class max # of tracks: max value of the data distribution
2082  { 1827, 1827}, // 0-5% cent class min # of tracks
2083  { 1498, 1498}, // 5-10
2084  { 1234, 1234}, // 10-15
2085  { 1012, 1012}, // 15-20
2086  { 827, 827}, // 20-25
2087  { 669, 669}, // 25-30
2088  { 536, 536}, // 30-35
2089  { 423, 423}, // 35-40
2090  { 329, 329}, // 40-45 cent class min # of tracks (cut digit 9)
2091  { 0, 0} // not used
2092  };
2093  Int_t PrimaryTracks5b[11][2] =
2094  {
2095  { 260, 214}, // 45
2096  { 197, 162}, // 50
2097  { 147, 125}, // 55
2098  { 106, 100}, // 60
2099  { 75, 63}, // 65
2100  { 51, 44}, // 70
2101  { 34, 29}, // 75
2102  { 21, 18}, // 80
2103  { 13, 11}, // 85
2104  { 0, 0}, // 90
2105  { 0, 0}// 100 only max accessible
2106  };
2107  Int_t PrimaryTracksLHC11h5b[11][2] =
2108  {
2109  { 214, 214}, // 45
2110  { 165, 162}, // 50
2111  { 127, 125}, // 55
2112  { 93, 100}, // 60
2113  { 64, 63}, // 65
2114  { 44, 44}, // 70
2115  { 30, 29}, // 75
2116  { 18, 18}, // 80
2117  { 11, 11}, // 85
2118  { 0, 0}, // 90
2119  { 0, 0}// 100 only max accessible
2120  };
2121  Int_t PrimaryTracksLHC15o5b[11][2] =
2122  {
2123  { 329, 329}, // 45-50% cent class max # of tracks
2124  { 251, 251}, // 45-50% cent class min # of tracks
2125  { 188, 188}, // 50-55
2126  { 136, 136}, // 55-60
2127  { 97, 97}, // 60-65
2128  { 67, 67}, // 65-70
2129  { 44, 44}, // 70-75
2130  { 28, 28}, // 75-80
2131  { 17, 17}, // 80-85
2132  { 0, 0}, // 85-90 cent class minimum # of tracks
2133  { 0, 0} // not used
2134  };
2135  Int_t PrimaryTracksLHC17n10[11][2] =
2136  {
2137  {9999,9999}, // 0 // 1500 max in hist but set to real max
2138  { 800, 800}, // 10 // guess
2139  { 628, 628}, // 20
2140  { 350, 350}, // 30 // guess
2141  { 268, 268}, // 40
2142  { 200, 200}, // 50 // guess
2143  { 100, 100}, // 60 // guess
2144  { 51, 44}, // 70 // guess
2145  { 21, 18}, // 80 // guess
2146  { 0, 0}, // 90 // guess
2147  { 0, 0} // 100 // only max accessible
2148  };
2149 
2150  Int_t column = 0;
2151  if(event->IsA()==AliESDEvent::Class()) column = 0;
2152  if(event->IsA()==AliAODEvent::Class()) column = 1;
2153 
2154  if (fModCentralityClass == 3){
2155  if(mcEvent){
2156  // setting specific arry for LHC11h for MC track mult
2158  if(nprimaryTracks > PrimaryTracksLHC11h10[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC11h10[fCentralityMin][column])
2159  return kTRUE;
2160  else return kFALSE;
2161  // setting specific arry for LHC17n for MC track mult
2162  } else if(fPeriodEnum == kLHC17j7 ){
2163  if(nprimaryTracks > PrimaryTracksLHC17n10[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC17n10[fCentralityMin][column])
2164  return kTRUE;
2165  else return kFALSE;
2166  // setting specific arry for LHC15o for MC track mult
2168  if(nprimaryTracks > PrimaryTracksLHC15o10[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC15o10[fCentralityMin][column])
2169  return kTRUE;
2170  else return kFALSE;
2171  // setting specific arry for LHC10h for MC track mult
2172  } else {
2173  if(nprimaryTracks > PrimaryTracks10[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks10[fCentralityMin][column])
2174  return kTRUE;
2175  else return kFALSE;
2176  }
2177  }
2178  else{
2179  centralityC= Int_t(centrality/10);
2180  if(centralityC >= fCentralityMin && centralityC < fCentralityMax)
2181  return kTRUE;
2182  else return kFALSE;
2183  }
2184  }
2185  else if (fModCentralityClass ==4){
2186  if(mcEvent){
2187  // setting specific arry for LHC11h for MC track mult
2189  if(nprimaryTracks > PrimaryTracksLHC11h5a[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC11h5a[fCentralityMin][column])
2190  return kTRUE;
2191  else return kFALSE;
2192  // setting specific arry for LHC15o for MC track mult
2194  if(nprimaryTracks > PrimaryTracksLHC15o5a[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC15o5a[fCentralityMin][column])
2195  return kTRUE;
2196  else return kFALSE;
2197  // setting specific arry for LHC10h for MC track mult
2198  } else {
2199  if(nprimaryTracks > PrimaryTracks5a[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks5a[fCentralityMin][column])
2200  return kTRUE;
2201  else return kFALSE;
2202  }
2203  }
2204  else{
2205  centralityC= Int_t(centrality);
2206  if(centralityC >= fCentralityMin*5 && centralityC < fCentralityMax*5){
2207  return kTRUE;
2208  } else return kFALSE;
2209  }
2210  }
2211  else if (fModCentralityClass ==5){
2212  if(mcEvent){
2213  // setting specific arry for LHC11h for MC track mult
2215  if(nprimaryTracks > PrimaryTracksLHC11h5b[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC11h5b[fCentralityMin][column])
2216  return kTRUE;
2217  else return kFALSE;
2218  // setting specific arry for LHC15o for MC track mult
2220  if(nprimaryTracks > PrimaryTracksLHC15o5b[fCentralityMax][column] && nprimaryTracks <= PrimaryTracksLHC15o5b[fCentralityMin][column])
2221  return kTRUE;
2222  else return kFALSE;
2223  // setting specific arry for LHC10h for MC track mult
2224  } else {
2225  if(nprimaryTracks > PrimaryTracks5b[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks5b[fCentralityMin][column])
2226  return kTRUE;
2227  else return kFALSE;
2228  }
2229  }
2230  else{
2231  centralityC= Int_t(centrality);
2232  if(centralityC >= ((fCentralityMin*5)+45) && centralityC < ((fCentralityMax*5)+45))
2233  return kTRUE;
2234  else return kFALSE;
2235  }
2236  }
2237 
2238  return kFALSE;
2239 }
2240 
2241 //________________________________________________________________________
2243  // Cut on z position of primary vertex
2244  Double_t fVertexZ=event->GetPrimaryVertex()->GetZ();
2245  Double_t fVertexZSPD = 0;
2246  AliESDEvent *fESDEvent=dynamic_cast<AliESDEvent*>(event);
2247  if(fESDEvent){
2248  fVertexZSPD = fESDEvent->GetPrimaryVertexSPD()->GetZ();
2249  }
2250  AliAODEvent *fAODEvent=dynamic_cast<AliAODEvent*>(event);
2251  if(fAODEvent){
2252  fVertexZSPD = fAODEvent->GetPrimaryVertexSPD()->GetZ();
2253  }
2254 
2255  if(TMath::Abs(fVertexZ)>fMaxVertexZ)return kFALSE;
2256 
2257 
2258  if (fPeriodEnum == kLHC11h){
2259  if (TMath::Abs(fVertexZ-fVertexZSPD) > 0.1) return kFALSE;
2260  }
2261  if (fIsHeavyIon == 2){
2262  if(!fUtils->IsVertexSelected2013pA(event)) return kFALSE;
2263  }
2264 
2265  return kTRUE;
2266 }
2267 
2268 //________________________________________________________________________
2270 {
2272  return kFALSE;
2273  TBits fIR1 = event->GetHeader()->GetIRInt1InteractionMap(); // IR1 contains V0 information (VIR)
2274  TBits fIR2 = event->GetHeader()->GetIRInt2InteractionMap(); // IR2 contains T0 information
2275  UShort_t bunchCrossings = event->GetBunchCrossNumber();
2277  for(Int_t i = 0; i<180;i++){
2278  if(fIR1.TestBitNumber(i))
2279  fHistoPastFutureBits->Fill((i*25)-90*25);
2280  }
2281  }
2282 
2283  Bool_t isOutOfBunchPileup = 0;
2284  Int_t pf1 = fPastFutureRejectionLow +bunchCrossings%4;
2285  Int_t pf2 = fPastFutureRejectionHigh+bunchCrossings%4;
2286  if(pf1 < -89) pf1 = -89;
2287  if(pf2 > 89) pf2 = 89;
2288  Int_t pf2maxForT0 = pf2;
2289  Int_t ir1skip = 0;
2290  for (Int_t i=pf1;i<=pf2;i++) {
2291  if (i==0) continue;
2292  if (i<=pf2maxForT0) isOutOfBunchPileup|=fIR2.TestBitNumber(90+i); // T0-based clean-up
2293  if (i>0 && i<=ir1skip) continue; // skip next 2 for old IR definitions
2294  isOutOfBunchPileup|=fIR1.TestBitNumber(90+i); // V0-based clean-up
2295  }
2296  return isOutOfBunchPileup;
2297 }
2298 //________________________________________________________________________
2299 
2301 {
2302  Bool_t isPileUpV0MTPCout=0;
2303 
2304  Double_t multV0M;
2305  Double_t valFunc;
2306  if (fIsHeavyIon==2){
2307  multV0M = event->GetVZEROData()->GetMTotV0A();
2308  }else{
2309  multV0M = event->GetVZEROData()->GetMTotV0A() + event->GetVZEROData()->GetMTotV0C() ;
2310  }
2311 
2312  if ( fFPileUpRejectV0MTPCout != 0x0 ){
2313  valFunc= fFPileUpRejectV0MTPCout->Eval(((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetNumberOfTPCoutTracks());
2314  if (multV0M < valFunc ) isPileUpV0MTPCout=1;
2315  }
2316 
2317  return isPileUpV0MTPCout;
2318 
2319 }
2320 //________________________________________________________________________
2322  // returns number of contributors to the vertex
2323 
2324  AliESDEvent *fESDEvent=dynamic_cast<AliESDEvent*>(event);
2325  if(fESDEvent){
2326  if (fESDEvent->GetPrimaryVertex() != NULL){
2327  if(fESDEvent->GetPrimaryVertex()->GetNContributors()>0) {
2328  // cout << "accepted global" << fESDEvent->GetEventNumberInFile() << " with NCont: " << fESDEvent->GetPrimaryVertex()->GetNContributors() << endl;
2329  return fESDEvent->GetPrimaryVertex()->GetNContributors();
2330  }
2331  }
2332 
2333  if(fESDEvent->GetPrimaryVertexSPD() !=NULL){
2334  if(fESDEvent->GetPrimaryVertexSPD()->GetNContributors()>0) {
2335  // cout << "accepted SPD" << fESDEvent->GetEventNumberInFile() << " with NCont: " << fESDEvent->GetPrimaryVertexSPD()->GetNContributors() << endl;
2336  return fESDEvent->GetPrimaryVertexSPD()->GetNContributors();
2337  }else {
2338  AliWarning(Form("Number of contributors from bad vertex type:: %s",fESDEvent->GetPrimaryVertex()->GetName()));
2339  // cout << "rejected " << fESDEvent->GetEventNumberInFile() << endl;
2340  return 0;
2341  }
2342  }
2343  }
2344 
2345  AliAODEvent *fAODEvent=dynamic_cast<AliAODEvent*>(event);
2346  if(fAODEvent){
2347  if (fAODEvent->GetPrimaryVertex() != NULL){
2348  if(fAODEvent->GetPrimaryVertex()->GetNContributors()>0) {
2349  return fAODEvent->GetPrimaryVertex()->GetNContributors();
2350  }
2351  }
2352  if(fAODEvent->GetPrimaryVertexSPD() !=NULL){
2353  if(fAODEvent->GetPrimaryVertexSPD()->GetNContributors()>0) {
2354  return fAODEvent->GetPrimaryVertexSPD()->GetNContributors();
2355  } else {
2356  AliWarning(Form("Number of contributors from bad vertex type:: %s",fAODEvent->GetPrimaryVertex()->GetName()));
2357  return 0;
2358  }
2359  }
2360  }
2361  // cout << "rejected " << fESDEvent->GetEventNumberInFile() << endl;
2362  return 0;
2363 }
2364 
2365 //________________________________________________________________________
2366 // Analysing Jet-Jet MC's
2367 //________________________________________________________________________
2369  AliGenCocktailEventHeader *cHeader = 0x0;
2370  Bool_t headerFound = kFALSE;
2371  weight = -1;
2372  fMaxPtJetMC = 0;
2373 
2374  if ( fPeriodEnum != kLHC17g8a_fast && fPeriodEnum != kLHC17g8a_cent_woSDD && // LHC16qt pPb 5TeV JetJet MC's
2375  fPeriodEnum != kLHC17g8b && fPeriodEnum != kLHC17g8c && // LHC16sr pPb 8TeV JetJet MC's
2376  fPeriodEnum != kLHC16P1JJ && fPeriodEnum != kLHC16P1JJLowB && // LHC16X Jet Jet MC's
2377  fPeriodEnum != kLHC16h3 && // LHC15n Jet Jet MC's
2378  fPeriodEnum != kLHC15a3a && fPeriodEnum != kLHC15a3a_plus && fPeriodEnum != kLHC15a3b && // LHC13g Jet Jet MC's
2379  fPeriodEnum != kLHC15g1a && fPeriodEnum != kLHC15g1b && // LHC11a Jet Jet MC's
2380  fPeriodEnum != kLHC13b4_fix && fPeriodEnum != kLHC13b4_plus && // LHC13 pPb Jet Jet MC's
2381  fPeriodEnum != kLHC16c3a && fPeriodEnum != kLHC16c3b && fPeriodEnum != kLHC16c3c && // LHC13 pPb Jet Jet MC's
2382  fPeriodEnum != kLHC16c2 && fPeriodEnum != kLHC16c2_plus // LHC12 JetJet MC
2383  ){
2384 
2385  weight = 1;
2386  return kTRUE;
2387  }
2388 
2389  if(mcEvent){
2390  cHeader = dynamic_cast<AliGenCocktailEventHeader*>(mcEvent->GenEventHeader());
2391  if(cHeader) headerFound = kTRUE;
2392  }else{
2393  //no mcEvent available -> not running on MC
2394  weight = 1;
2395  return kTRUE;
2396  }
2397 
2398  if(headerFound){
2399  TList *genHeaders = 0x0;
2400  if(cHeader) genHeaders = cHeader->GetHeaders();
2401  AliGenEventHeader* gh = 0;
2402  for(Int_t i = 0; i<genHeaders->GetEntries();i++){
2403  gh = (AliGenEventHeader*)genHeaders->At(i);
2404  TString GeneratorName = gh->GetName();
2405  if (GeneratorName.CompareTo("AliGenPythiaEventHeader") == 0){
2406  Bool_t eventAccepted = kTRUE;
2407  TParticle * jet = 0;
2408  Int_t nTriggerJets = dynamic_cast<AliGenPythiaEventHeader*>(gh)->NTriggerJets();
2409  Float_t ptHard = dynamic_cast<AliGenPythiaEventHeader*>(gh)->GetPtHard();
2410  Float_t tmpjet[]={0,0,0,0};
2411  for(Int_t ijet = 0; ijet< nTriggerJets; ijet++){
2412  dynamic_cast<AliGenPythiaEventHeader*>(gh)->TriggerJet(ijet, tmpjet);
2413  jet = new TParticle(94, 21, -1, -1, -1, -1, tmpjet[0],tmpjet[1],tmpjet[2],tmpjet[3], 0,0,0,0);
2414  //Compare jet pT and pt Hard
2415  if(jet->Pt() > fMaxFacPtHard * ptHard){
2416  eventAccepted= kFALSE;
2417  }
2418  if (jet->Pt() > fMaxPtJetMC) fMaxPtJetMC = jet->Pt();
2419  }
2420  if (jet) delete jet;
2421  if (mcEvent){
2422  for(Long_t i = 0; i < mcEvent->GetNumberOfPrimaries(); i++) {
2423  TParticle* particle = (TParticle *)mcEvent->Particle(i);
2424  if (!particle) continue;
2425  if (TMath::Abs(particle->GetPdgCode()) == 111 || TMath::Abs(particle->GetPdgCode()) == 221){
2426  if (particle->Pt() > fMaxFacPtHardSingleParticle*ptHard){
2427  eventAccepted= kFALSE;
2428  }
2429  }
2430 
2431  }
2432  }
2433 
2434  if ( fPeriodEnum == kLHC16P1JJLowB) {
2435  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2436  21, 28, 36, 45, 57,
2437  70, 85, 99, 115, 132,
2438  150, 169, 190, 212, 235,
2439  1000000};
2440  Double_t weightsBins[20] = { 43.7553, 13.5848, 6.788, 2.67826, 0.975255,
2441  0.39069, 0.127342, 0.0465597, 0.0206539, 0.00750243,
2442  0.00319118, 0.00122291, 0.000641232, 0.000321437, 0.000168273,
2443  9.17033e-05, 5.34755e-05, 3.01354e-05, 1.74518e-05, 2.8004e-05};
2444  Int_t bin = 0;
2445  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2446  if (bin < 20) weight = weightsBins[bin];
2447 
2448  } else if ( fPeriodEnum == kLHC16P1JJ ){
2449  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2450  21, 28, 36, 45, 57,
2451  70, 85, 99, 115, 132,
2452  150, 169, 190, 212, 235,
2453  1000000};
2454  Double_t weightsBins[20] = { 43.8654, 13.6215, 6.79856, 2.67526, 0.978794,
2455  0.390797, 0.127769, 0.0465714, 0.0206173, 0.00750282,
2456  0.00318773, 0.00122533, 0.000644385, 0.000321225, 0.00016846,
2457  9.18305e-05, 5.33507e-05, 3.00677e-05, 1.74608e-05, 2.80823e-05};
2458 
2459  Int_t bin = 0;
2460  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2461  if (bin < 20) weight = weightsBins[bin];
2462 
2463  } else if ( fPeriodEnum == kLHC16h3 ){
2464  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2465  21, 28, 36, 45, 57,
2466  70, 85, 99, 115, 132,
2467  150, 169, 190, 212, 235,
2468  1000000};
2469  Double_t weightsBins[20] = { 16.0869, 4.61169, 2.14976, 0.782544, 0.264854,
2470  9.7619E-02, 2.92747E-02, 9.89515E-03, 4.05152E-03, 1.35393E-03,
2471  5.29864E-04, 1.88317E-04, 9.23E-05, 4.29E-05, 2.09E-05,
2472  1.06E-05, 5.76E-06, 3.00E-06, 1.62E-06, 2.10E-06 };
2473  Int_t bin = 0;
2474  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2475  if (bin < 20) weight = weightsBins[bin];
2476 
2477 
2478  } else if ( fPeriodEnum == kLHC15a3b || fPeriodEnum == kLHC15g1b ){
2479  Double_t ptHardBinRanges[13] = { 5, 7, 9, 12, 16,
2480  21, 28, 36, 45, 57,
2481  70, 85, 1000};
2482  Double_t weightsBins[12] = { 7.858393e-03, 4.718691e-03, 4.077575e-03, 2.814527e-03, 1.669625e-03,
2483  1.007535e-03, 4.536554e-04, 2.111041e-04, 1.094840e-04, 4.404973e-05,
2484  1.933238e-05, 1.562895e-05};
2485  Int_t bin = 0;
2486  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2487  if (bin < 12) weight = weightsBins[bin];
2488 
2489  } else if ( fPeriodEnum == kLHC15g1a ){
2490  Double_t ptHardBinRanges[20] = { 5, 11, 21, 36, 57,
2491  84, 117, 152, 191, 1000000,
2492  5, 7, 9, 12, 16,
2493  21, 28, 36, 45, 57 };
2494  Double_t weightsBins[19] = { 4.43629 , 0.49523, 0.0394921, 0.00383174, 0.000446559,
2495  6.37374e-05, 1.03134e-05, 2.27012e-06, 7.59281e-07, 0,
2496  2.62906, 1.12884, 0.656873, 0.262822, 0.0876732,
2497  0.0307759, 0.0087083, 0.0027664, 0.00106203};
2498 
2499  Int_t bin = 0;
2500  Int_t binFromFile = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetPtHardFromFile();
2501  if (binFromFile != -1 && binFromFile >9 && ptHard < 57) bin = 9;
2502  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2503  if (bin < 19) weight = weightsBins[bin];
2504 
2505  } else if ( fPeriodEnum == kLHC15a3a || fPeriodEnum == kLHC15a3a_plus ) {
2506  Double_t ptHardBinRanges[20] = { 5, 11, 21, 36, 57,
2507  84, 117, 152, 191, 1000000,
2508  5, 7, 9, 12, 16,
2509  21, 28, 36, 45, 57 };
2510  // LHC15a3a
2511  Double_t weightsBins[19] = { 4.43897 , 0.495766, 0.039486, 0.00383011, 0.000447104,
2512  6.37277e-05, 1.03166e-05, 2.26971e-06, 7.59023e-07, 0,
2513  2.63331, 1.12815, 0.657034, 0.262756, 0.0877227,
2514  0.0307638, 0.00870635, 0.00276658, 0.00106229};
2515  Int_t bin = 0;
2516  Int_t binFromFile = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetPtHardFromFile();
2517  if (binFromFile != -1 && binFromFile >9 && ptHard < 57) bin = 9;
2518  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2519  if (bin < 19) weight = weightsBins[bin];
2520 
2521  } else if ( fPeriodEnum == kLHC16c2 || fPeriodEnum == kLHC16c2_plus ){
2522  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2523  21, 28, 36, 45, 57,
2524  70, 85, 99, 115, 132,
2525  150, 169, 190, 212, 235,
2526  1000000};
2527  Double_t weightsBins[20] = { 28.3084, 8.43277, 4.07753, 1.54359, 0.543318,
2528  0.208394, 0.0652349, 0.0186904, 0.00834528, 0.00301414,
2529  0.00125939, 0.000474403, 0.000244052, 0.00011924, 6.09838e-05,
2530  3.24148e-05, 1.84314e-05, 1.00926e-05, 5.68632e-06, 8.38092e-06};
2531  Int_t bin = 0;
2532  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2533  if (bin < 20) weight = weightsBins[bin];
2534 
2535  } else if ( fPeriodEnum == kLHC16c3a ){
2536  Double_t ptHardBinRanges[6] = { 7, 9, 12, 16, 21, 1000};
2537  Double_t weightsBins[5] = { 0.00672445, 0.00799158, 0.00678934, 0.00463908, 0.00600068}; //preliminary estimates
2538  Int_t bin = 0;
2539  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2540  if (bin < 5) weight = weightsBins[bin];
2541 
2542  } else if (fPeriodEnum == kLHC16c3b ){
2543  Double_t ptHardBinRanges[7] = { 14, 19, 26, 35, 48, 66, 1000};
2544  Double_t weightsBins[6] = { 0.00608281, 0.00393646, 0.00200138, 0.000986267, 0.000389051, 0.0001863}; //preliminary estimates
2545  Int_t bin = 0;
2546  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2547  if (bin < 6) weight = weightsBins[bin];
2548 
2549  } else if (fPeriodEnum == kLHC16c3c ){
2550  Double_t ptHardBinRanges[8] = { 0, 5, 11, 21, 36, 57, 84, 1000};
2551  Double_t weightsBins[7] = { 0.00151999, 0.000100346, 1.27688e-05, 1.82388e-06, 3.08506e-07, 6.00308e-08, 1.88414e-08}; //preliminary estimates
2552  Int_t bin = 0;
2553  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2554  if (bin < 7) weight = weightsBins[bin];
2555 
2556  } else if ( fPeriodEnum == kLHC13b4_fix || fPeriodEnum == kLHC13b4_plus ){
2557  Double_t ptHardBinRanges[11] = { 5, 11, 21, 36, 57,
2558  84, 117, 152, 191, 234,
2559  1000};
2560  Double_t weightsBins[10] = { 2.24185e-6 , 2.48463e-7, 2.23171e-8, 2.43667e-9, 3.29934e-10,
2561  5.34592e-11, 1.00937e-11, 2.6493e-12, 8.53912e-13, 5.43077e-13};
2562  Int_t bin = 0;
2563  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2564  if (bin < 10) weight = weightsBins[bin];
2565 
2566  } else {
2567  weight = 1;
2568  }
2569 
2570  if (weight == -1) return kFALSE;
2571  else return eventAccepted;
2572 
2573  }
2574  }
2575  } else {
2576  AliGenEventHeader * eventHeader = mcEvent->GenEventHeader();
2577  TString eventHeaderName = eventHeader->ClassName();
2578  if (eventHeaderName.CompareTo("AliGenPythiaEventHeader") == 0){
2579  Bool_t eventAccepted = kTRUE;
2580  TParticle * jet = 0;
2581  Int_t nTriggerJets = dynamic_cast<AliGenPythiaEventHeader*>(eventHeader)->NTriggerJets();
2582  Float_t ptHard = dynamic_cast<AliGenPythiaEventHeader*>(eventHeader)->GetPtHard();
2583  Float_t tmpjet[]={0,0,0,0};
2584  for(Int_t ijet = 0; ijet< nTriggerJets; ijet++){
2585  dynamic_cast<AliGenPythiaEventHeader*>(eventHeader)->TriggerJet(ijet, tmpjet);
2586  jet = new TParticle(94, 21, -1, -1, -1, -1, tmpjet[0],tmpjet[1],tmpjet[2],tmpjet[3], 0,0,0,0);
2587  //Compare jet pT and pt Hard
2588  if(jet->Pt() > fMaxFacPtHard * ptHard){
2589  eventAccepted= kFALSE;
2590  }
2591  if (jet->Pt() > fMaxPtJetMC) fMaxPtJetMC = jet->Pt();
2592  }
2593  if (mcEvent){
2594  for(Long_t i = 0; i < mcEvent->GetNumberOfPrimaries(); i++) {
2595  TParticle* particle = (TParticle *)mcEvent->Particle(i);
2596  if (!particle) continue;
2597  if (TMath::Abs(particle->GetPdgCode()) == 111 || TMath::Abs(particle->GetPdgCode()) == 221){
2598  if (particle->Pt() > fMaxFacPtHardSingleParticle*ptHard){
2599  eventAccepted= kFALSE;
2600  }
2601  }
2602 
2603  }
2604  }
2605 
2606  if ( fPeriodEnum == kLHC16P1JJLowB) {
2607  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2608  21, 28, 36, 45, 57,
2609  70, 85, 99, 115, 132,
2610  150, 169, 190, 212, 235,
2611  1000000};
2612  Double_t weightsBins[20] = { 43.7553, 13.5848, 6.788, 2.67826, 0.975255,
2613  0.39069, 0.127342, 0.0465597, 0.0206539, 0.00750243,
2614  0.00319118, 0.00122291, 0.000641232, 0.000321437, 0.000168273,
2615  9.17033e-05, 5.34755e-05, 3.01354e-05, 1.74518e-05, 2.8004e-05};
2616  Int_t bin = 0;
2617  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2618  if (bin < 20) weight = weightsBins[bin];
2619 
2620  } else if ( fPeriodEnum == kLHC16P1JJ ){
2621  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2622  21, 28, 36, 45, 57,
2623  70, 85, 99, 115, 132,
2624  150, 169, 190, 212, 235,
2625  1000000};
2626  Double_t weightsBins[20] = { 43.8654, 13.6215, 6.79856, 2.67526, 0.978794,
2627  0.390797, 0.127769, 0.0465714, 0.0206173, 0.00750282,
2628  0.00318773, 0.00122533, 0.000644385, 0.000321225, 0.00016846,
2629  9.18305e-05, 5.33507e-05, 3.00677e-05, 1.74608e-05, 2.80823e-05};
2630 
2631  Int_t bin = 0;
2632  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2633  if (bin < 20) weight = weightsBins[bin];
2634 
2635  } else if ( fPeriodEnum == kLHC16h3 ){
2636  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2637  21, 28, 36, 45, 57,
2638  70, 85, 99, 115, 132,
2639  150, 169, 190, 212, 235,
2640  1000000};
2641  Double_t weightsBins[20] = { 16.0869, 4.61169, 2.14976, 0.782544, 0.264854,
2642  9.7619E-02, 2.92747E-02, 9.89515E-03, 4.05152E-03, 1.35393E-03,
2643  5.29864E-04, 1.88317E-04, 9.23E-05, 4.29E-05, 2.09E-05,
2644  1.06E-05, 5.76E-06, 3.00E-06, 1.62E-06, 2.10E-06 };
2645  Int_t bin = 0;
2646  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2647  if (bin < 20) weight = weightsBins[bin];
2648 
2649 
2650  } else if ( fPeriodEnum == kLHC15a3b || fPeriodEnum == kLHC15g1b ){
2651  Double_t ptHardBinRanges[13] = { 5, 7, 9, 12, 16,
2652  21, 28, 36, 45, 57,
2653  70, 85, 1000};
2654  Double_t weightsBins[12] = { 7.858393e-03, 4.718691e-03, 4.077575e-03, 2.814527e-03, 1.669625e-03,
2655  1.007535e-03, 4.536554e-04, 2.111041e-04, 1.094840e-04, 4.404973e-05,
2656  1.933238e-05, 1.562895e-05};
2657  Int_t bin = 0;
2658  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2659  if (bin < 12) weight = weightsBins[bin];
2660 
2661  } else if ( fPeriodEnum == kLHC15g1a ){
2662  Double_t ptHardBinRanges[20] = { 5, 11, 21, 36, 57,
2663  84, 117, 152, 191, 1000000,
2664  5, 7, 9, 12, 16,
2665  21, 28, 36, 45, 57 };
2666  Double_t weightsBins[19] = { 4.43629 , 0.49523, 0.0394921, 0.00383174, 0.000446559,
2667  6.37374e-05, 1.03134e-05, 2.27012e-06, 7.59281e-07, 0,
2668  2.62906, 1.12884, 0.656873, 0.262822, 0.0876732,
2669  0.0307759, 0.0087083, 0.0027664, 0.00106203};
2670 
2671  Int_t bin = 0;
2672  Int_t binFromFile = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetPtHardFromFile();
2673  if (binFromFile != -1 && binFromFile >9 && ptHard < 57) bin = 9;
2674  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2675  if (bin < 19) weight = weightsBins[bin];
2676 
2677  } else if ( fPeriodEnum == kLHC15a3a || fPeriodEnum == kLHC15a3a_plus ) {
2678  Double_t ptHardBinRanges[20] = { 5, 11, 21, 36, 57,
2679  84, 117, 152, 191, 1000000,
2680  5, 7, 9, 12, 16,
2681  21, 28, 36, 45, 57 };
2682  // LHC15a3a
2683  Double_t weightsBins[19] = { 4.43897 , 0.495766, 0.039486, 0.00383011, 0.000447104,
2684  6.37277e-05, 1.03166e-05, 2.26971e-06, 7.59023e-07, 0,
2685  2.63331, 1.12815, 0.657034, 0.262756, 0.0877227,
2686  0.0307638, 0.00870635, 0.00276658, 0.00106229};
2687  Int_t bin = 0;
2688  Int_t binFromFile = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetPtHardFromFile();
2689  if (binFromFile != -1 && binFromFile >9 && ptHard < 57) bin = 9;
2690  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2691  if (bin < 19) weight = weightsBins[bin];
2692 
2693  } else if ( fPeriodEnum == kLHC16c2 || fPeriodEnum == kLHC16c2_plus ){
2694  Double_t ptHardBinRanges[21] = { 5, 7, 9, 12, 16,
2695  21, 28, 36, 45, 57,
2696  70, 85, 99, 115, 132,
2697  150, 169, 190, 212, 235,
2698  1000000};
2699  Double_t weightsBins[20] = { 28.3084, 8.43277, 4.07753, 1.54359, 0.543318,
2700  0.208394, 0.0652349, 0.0186904, 0.00834528, 0.00301414,
2701  0.00125939, 0.000474403, 0.000244052, 0.00011924, 6.09838e-05,
2702  3.24148e-05, 1.84314e-05, 1.00926e-05, 5.68632e-06, 8.38092e-06};
2703  Int_t bin = 0;
2704  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2705  if (bin < 20) weight = weightsBins[bin];
2706 
2707  } else if ( fPeriodEnum == kLHC16c3a ){
2708  Double_t ptHardBinRanges[6] = { 7, 9, 12, 16, 21, 1000};
2709  Double_t weightsBins[5] = { 0.00672445, 0.00799158, 0.00678934, 0.00463908, 0.00600068}; //preliminary estimates
2710  Int_t bin = 0;
2711  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2712  if (bin < 5) weight = weightsBins[bin];
2713 
2714  } else if (fPeriodEnum == kLHC16c3b ){
2715  Double_t ptHardBinRanges[7] = { 14, 19, 26, 35, 48, 66, 1000};
2716  Double_t weightsBins[6] = { 0.00608281, 0.00393646, 0.00200138, 0.000986267, 0.000389051, 0.0001863}; //preliminary estimates
2717  Int_t bin = 0;
2718  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2719  if (bin < 6) weight = weightsBins[bin];
2720 
2721  } else if (fPeriodEnum == kLHC16c3c ){
2722  Double_t ptHardBinRanges[8] = { 0, 5, 11, 21, 36, 57, 84, 1000};
2723  Double_t weightsBins[7] = { 0.00151999, 0.000100346, 1.27688e-05, 1.82388e-06, 3.08506e-07, 6.00308e-08, 1.88414e-08}; //preliminary estimates
2724  Int_t bin = 0;
2725  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2726  if (bin < 7) weight = weightsBins[bin];
2727 
2728  } else if ( fPeriodEnum == kLHC13b4_fix || fPeriodEnum == kLHC13b4_plus ){
2729  Double_t ptHardBinRanges[11] = { 5, 11, 21, 36, 57,
2730  84, 117, 152, 191, 234,
2731  1000};
2732  Double_t weightsBins[10] = { 2.24185e-6 , 2.48463e-7, 2.23171e-8, 2.43667e-9, 3.29934e-10,
2733  5.34592e-11, 1.00937e-11, 2.6493e-12, 8.53912e-13, 5.43077e-13};
2734  Int_t bin = 0;
2735  while (!((ptHard< ptHardBinRanges[bin+1] && ptHard > ptHardBinRanges[bin]) || (ptHard == ptHardBinRanges[bin]) ) )bin++;
2736  if (bin < 10) weight = weightsBins[bin];
2737  } else {
2738  weight = 1;
2739  }
2740 
2741  if (weight == -1) return kFALSE;
2742  else return eventAccepted;
2743 
2744  } else {
2745  return kFALSE;
2746  }
2747  }
2748 
2749  return kFALSE;
2750 }
2751 
2752 //________________________________________________________________________
2753 // Analysing Jet-Jet MC's
2754 //________________________________________________________________________
2755 void AliConvEventCuts::GetXSectionAndNTrials(AliMCEvent *mcEvent, Float_t &XSection, Float_t &NTrials){
2756 
2757  AliGenCocktailEventHeader *cHeader = 0x0;
2758  Bool_t headerFound = kFALSE;
2759 
2760  if ( fPeriodEnum != kLHC17g8a_fast && fPeriodEnum != kLHC17g8a_cent_woSDD && // LHC16qt pPb 5TeV JetJet MC's
2761  fPeriodEnum != kLHC17g8b && fPeriodEnum != kLHC17g8c && // LHC16sr pPb 8TeV JetJet MC's
2762  fPeriodEnum != kLHC16P1JJ && fPeriodEnum != kLHC16P1JJLowB && // LHC16X Jet Jet MC's
2763  fPeriodEnum != kLHC16h3 && // LHC15n Jet Jet MC's
2764  fPeriodEnum != kLHC15a3a && fPeriodEnum != kLHC15a3a_plus && fPeriodEnum != kLHC15a3b && // LHC13g Jet Jet MC's
2765  fPeriodEnum != kLHC15g1a && fPeriodEnum != kLHC15g1b && // LHC11a Jet Jet MC's
2766  fPeriodEnum != kLHC13b4_fix && fPeriodEnum != kLHC13b4_plus && // LHC13 pPb Jet Jet MC's
2767  fPeriodEnum != kLHC16c3a && fPeriodEnum != kLHC16c3b && fPeriodEnum != kLHC16c3c && // LHC13 pPb Jet Jet MC's
2768  fPeriodEnum != kLHC16c2 && fPeriodEnum != kLHC16c2_plus // LHC12 JetJet MC
2769  ){
2770  NTrials = -1;
2771  XSection = -1;
2772  return;
2773  }
2774 
2775  if(mcEvent){
2776  cHeader = dynamic_cast<AliGenCocktailEventHeader*>(mcEvent->GenEventHeader());
2777  if(cHeader) headerFound = kTRUE;
2778  }else{
2779  //no mcEvent available -> not running on MC
2780  NTrials = -1;
2781  XSection = -1;
2782  return;
2783  }
2784 
2785  if(headerFound){
2786  TList *genHeaders = 0x0;
2787  if(cHeader) genHeaders = cHeader->GetHeaders();
2788  AliGenEventHeader* gh = 0;
2789  for(Int_t i = 0; i<genHeaders->GetEntries();i++){
2790  gh = (AliGenEventHeader*)genHeaders->At(i);
2791  TString GeneratorName = gh->GetName();
2792  if (GeneratorName.CompareTo("AliGenPythiaEventHeader") == 0){
2793  AliGenPythiaEventHeader* gPythia = dynamic_cast<AliGenPythiaEventHeader*>(gh);
2794  NTrials = gPythia->Trials();
2795  XSection = gPythia->GetXsection();
2796  return;
2797  }
2798  }
2799  } else {
2800  AliGenEventHeader * eventHeader = mcEvent->GenEventHeader();
2801  if(eventHeader){
2802  TString eventHeaderName = eventHeader->ClassName();
2803  if (eventHeaderName.CompareTo("AliGenPythiaEventHeader") == 0){
2804  AliGenPythiaEventHeader* gPythia = dynamic_cast<AliGenPythiaEventHeader*>(eventHeader);
2805  NTrials = gPythia->Trials();
2806  XSection = gPythia->GetXsection();
2807  return;
2808  }
2809  }
2810  }
2811 
2812  NTrials = -1;
2813  XSection = -1;
2814  return;
2815 }
2816 
2817 
2818 //________________________________________________________________________
2819 // Analysing Jet-Jet MC's
2820 //________________________________________________________________________
2822  AliGenCocktailEventHeader *cHeader = 0x0;
2823  Bool_t headerFound = kFALSE;
2824 
2825  if ( fPeriodEnum != kLHC17g8a_fast && fPeriodEnum != kLHC17g8a_cent_woSDD && // LHC16qt pPb 5TeV JetJet MC's
2826  fPeriodEnum != kLHC17g8b && fPeriodEnum != kLHC17g8c && // LHC16sr pPb 8TeV JetJet MC's
2827  fPeriodEnum != kLHC16P1JJ && fPeriodEnum != kLHC16P1JJLowB && // LHC16X Jet Jet MC's
2828  fPeriodEnum != kLHC16h3 && // LHC15n Jet Jet MC's
2829  fPeriodEnum != kLHC15a3a && fPeriodEnum != kLHC15a3a_plus && fPeriodEnum != kLHC15a3b && // LHC13g Jet Jet MC's
2830  fPeriodEnum != kLHC15g1a && fPeriodEnum != kLHC15g1b && // LHC11a Jet Jet MC's
2831  fPeriodEnum != kLHC13b4_fix && fPeriodEnum != kLHC13b4_plus && // LHC13 pPb Jet Jet MC's
2832  fPeriodEnum != kLHC16c3a && fPeriodEnum != kLHC16c3b && fPeriodEnum != kLHC16c3c && // LHC13 pPb Jet Jet MC's
2833  fPeriodEnum != kLHC16c2 && fPeriodEnum != kLHC16c2_plus // LHC12 JetJet MC
2834  ) return -1;
2835 
2836  if(mcEvent){
2837  cHeader = dynamic_cast<AliGenCocktailEventHeader*>(mcEvent->GenEventHeader());
2838  if(cHeader) headerFound = kTRUE;
2839  }else{
2840  //no mcEvent available -> not running on MC
2841  return -1;
2842  }
2843 
2844  if(headerFound){
2845  TList *genHeaders = 0x0;
2846  if(cHeader) genHeaders = cHeader->GetHeaders();
2847  AliGenEventHeader* gh = 0;
2848  for(Int_t i = 0; i<genHeaders->GetEntries();i++){
2849  gh = (AliGenEventHeader*)genHeaders->At(i);
2850  TString GeneratorName = gh->GetName();
2851  if (GeneratorName.CompareTo("AliGenPythiaEventHeader") == 0){
2852  return dynamic_cast<AliGenPythiaEventHeader*>(gh)->GetPtHard();
2853  }
2854  }
2855  } else {
2856  AliGenEventHeader * eventHeader = mcEvent->GenEventHeader();
2857  if(eventHeader){
2858  TString eventHeaderName = eventHeader->ClassName();
2859  if (eventHeaderName.CompareTo("AliGenPythiaEventHeader") == 0){
2860  return dynamic_cast<AliGenPythiaEventHeader*>(eventHeader)->GetPtHard();
2861  }
2862  }
2863  }
2864 
2865  return -1;
2866 }
2867 
2868 
2869 //________________________________________________________________________
2871  // abort if mimicing not enabled
2872 
2873  if (!fMimicTrigger) return kTRUE;
2874 
2875  Int_t runRangesEMCalL0 [35] = { 144871, 145288, 146375, 146382, // LHC11a
2876  146502, 148522, // LHC11a
2877  150209, 153056, 153911, 153915, // LHC11b,c,d
2878  158135, 158136, 158178, 158182, 160683,
2879  160764, 161139, 161256, 161379, 161457,
2880  161525, 161556, 161558, 161609, 161630,
2881  161724, // LHC11d,e
2882  173731, 177144, 177147, 177653, 177724, 178327,
2883  195180, // LHC13b-f
2884  197469, 197692 // LHC13g
2885  };
2886 
2887  Double_t thresholdEMCalL0[34] = { 2.11, 3.43, 1.71, 2.05, // LHC11a 7 TeV
2888  3.43, // LHC11a 2.76TeV
2889  1.94, 3.39, 4.01, 5.25, 5.5, // LHC11b, LHC11c, LHC11d
2890  2.05, 5.50, 2.05, 5.50, 2.05, 1.71, 5.50, 1.71, 5.50, 1.71, 5.50, 1.71, 5.50, 1.71, 5.50, 1.71,
2891  2.01, 1.75, 1.52, 2.01, 1.52, 1.85,
2892  3.2,
2893  /*2.01*/1.8
2894  };
2895  Double_t spreadEMCalL0[34] = { 0., 0., 0, 0, // LHC11a 7TeV
2896  /*0.7*/0.65, // LHC11a 2.76TeV
2897  0., 0., 0., 0., 0., // LHC11b, LHC11c, LHC11d
2898  0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
2899  0., 0., 0., 0., 0.2, 0.2,/*0.,0.,*/
2900  0.1,
2901  /*0.1*/0.12
2902  };
2903 
2904  Int_t runRangesEMCalL1[4] = { 179796, // LHC12c-i
2905  195180, // LHC13b-f
2906  197469, 197692 // LHC13g
2907  };
2908 
2909  Double_t thresholdEMCalL1[3] = { 9.5/*8.398*/, 11.5, /*6.*/5.5};
2910  Double_t spreadEMCalL1[3] = { 1.0/*0.*/, 0.5, /*0.4*/0.6};
2911 
2912  Int_t runRangesEMCalL1G2[3] = { 195180, // LHC13b-f
2913  197469, 197692 // LHC13g
2914  };
2915 
2916  Double_t thresholdEMCalL1G2[2] = { 7.2, /*3.9*/3.75};
2917  Double_t spreadEMCalL1G2[2] = { 0.3, /*0.2*/0.25};
2918 
2919  Int_t runnumber = event->GetRunNumber();
2920 
2921  if (fSpecialTrigger == 5 ){
2922  if (runnumber < runRangesEMCalL0[0]) return kTRUE;
2923  Int_t binRun = 0;
2924  while (!(runnumber >= runRangesEMCalL0[binRun] && runnumber < runRangesEMCalL0[binRun+1] ) && binRun < 34 ){
2925 // cout << runnumber << "\t" << binRun << "\t" << runRangesEMCalL0[binRun] << "\t" << runRangesEMCalL0[binRun+1] << endl;
2926  binRun++;
2927  }
2928  if (binRun==34) return kFALSE;
2929  Double_t threshold = thresholdEMCalL0[binRun];
2930 
2931  if (isMC && spreadEMCalL0[binRun] != 0.){
2932  TF1* triggerSmearing = new TF1("triggerSmearing","[0]*exp(-0.5*((x-[1])/[2])**2)",0,15);
2933  triggerSmearing->SetParameter(0, 1/(spreadEMCalL0[binRun]*TMath::Sqrt(TMath::Pi()*2)));
2934  triggerSmearing->SetParameter(1, thresholdEMCalL0[binRun]);
2935  triggerSmearing->SetParameter(2, spreadEMCalL0[binRun]);
2936  threshold = triggerSmearing->GetRandom();
2937  delete triggerSmearing;
2938  }
2939 
2940 // cout << runnumber << "\t"<< binRun << "\t"<< threshold << endl;
2941 
2942  Int_t nclus = 0;
2943  nclus = event->GetNumberOfCaloClusters();
2944 
2945  if(nclus == 0) return kFALSE;
2946 
2947  // Loop over EMCal clusters
2948  Bool_t eventIsAccepted = kFALSE;
2949  for(Int_t i = 0; i < nclus; i++){
2950  AliVCluster* clus = NULL;
2951  clus = event->GetCaloCluster(i);
2952  if (!clus) continue;
2953  if (!clus->IsEMCAL()) continue;
2954  if (clus->GetM02()<0.1) continue;
2955  if (clus->GetNCells()<2) continue;
2956  if (clus->E() > threshold ){
2957 // cout << "found L0" << endl;
2958  eventIsAccepted = kTRUE;
2959  }
2960  }
2961  return eventIsAccepted;
2962 
2963  } else if (fSpecialTrigger == 6 ) {
2964 
2965  return kTRUE;
2966  } else if (fSpecialTrigger == 8 ) {
2967  if (fSpecialSubTriggerName.CompareTo("7EGA")==0 || fSpecialSubTriggerName.CompareTo("8EGA")==0 || fSpecialSubTriggerName.CompareTo("7EG1")==0 ||fSpecialSubTriggerName.CompareTo("8EG1")==0 ){
2968  if (runnumber < runRangesEMCalL1[0]) return kTRUE;
2969  Int_t binRun = 0;
2970  while (!(runnumber >= runRangesEMCalL1[binRun] && runnumber < runRangesEMCalL1[binRun+1] ) && binRun < 3 ){
2971  // cout << runnumber << "\t" << binRun << "\t" << runRangesEMCalL0[binRun] << "\t" << runRangesEMCalL0[binRun+1] << endl;
2972  binRun++;
2973  }
2974  if (binRun==3) return kFALSE;
2975  Double_t threshold = thresholdEMCalL1[binRun];
2976 
2977  if (isMC && spreadEMCalL1[binRun] != 0.){
2978  TF1* triggerSmearing = new TF1("triggerSmearing","[0]*exp(-0.5*((x-[1])/[2])**2)",0,15);
2979  triggerSmearing->SetParameter(0, 1/(spreadEMCalL1[binRun]*TMath::Sqrt(TMath::Pi()*2)));
2980  triggerSmearing->SetParameter(1, thresholdEMCalL1[binRun]);
2981  triggerSmearing->SetParameter(2, spreadEMCalL1[binRun]);
2982  threshold = triggerSmearing->GetRandom();
2983  delete triggerSmearing;
2984  }
2985 
2986 // cout << runnumber << "\t"<< binRun << "\t L1 \t"<< threshold << endl;
2987 
2988  Int_t nclus = 0;
2989  nclus = event->GetNumberOfCaloClusters();
2990 
2991  if(nclus == 0) return kFALSE;
2992 
2993  // Loop over EMCal clusters
2994  Bool_t eventIsAccepted = kFALSE;
2995  for(Int_t i = 0; i < nclus; i++){
2996  AliVCluster* clus = NULL;
2997  clus = event->GetCaloCluster(i);
2998  if (!clus) continue;
2999  if (!clus->IsEMCAL()) continue;
3000  if (clus->GetM02()<0.1) continue;
3001  if (clus->GetNCells()<2) continue;
3002  if (clus->E() > threshold ){
3003 // cout << "found L1G1" << endl;
3004  eventIsAccepted = kTRUE;
3005  }
3006  }
3007  return eventIsAccepted;
3008  } else if ( fSpecialSubTriggerName.CompareTo("7EG2")==0 ||fSpecialSubTriggerName.CompareTo("8EG2")==0 ){
3009  if (runnumber < runRangesEMCalL1G2[0]) return kTRUE;
3010  Int_t binRun = 0;
3011  while (!(runnumber >= runRangesEMCalL1G2[binRun] && runnumber < runRangesEMCalL1G2[binRun+1] ) && binRun < 2 ){
3012  // cout << runnumber << "\t" << binRun << "\t" << runRangesEMCalL0[binRun] << "\t" << runRangesEMCalL0[binRun+1] << endl;
3013  binRun++;
3014  }
3015  if (binRun==2) return kFALSE;
3016  Double_t threshold = thresholdEMCalL1G2[binRun];
3017  if (isMC && spreadEMCalL1G2[binRun] != 0.){
3018  TF1* triggerSmearing = new TF1("triggerSmearing","[0]*exp(-0.5*((x-[1])/[2])**2)",0,15);
3019  triggerSmearing->SetParameter(0, 1/(spreadEMCalL1G2[binRun]*TMath::Sqrt(TMath::Pi()*2)));
3020  triggerSmearing->SetParameter(1, thresholdEMCalL1G2[binRun]);
3021  triggerSmearing->SetParameter(2, spreadEMCalL1G2[binRun]);
3022  threshold = triggerSmearing->GetRandom();
3023  delete triggerSmearing;
3024  }
3025 // cout << runnumber << "\t"<< binRun << "\t L2 \t"<< threshold << endl;
3026 
3027  Int_t nclus = 0;
3028  nclus = event->GetNumberOfCaloClusters();
3029 
3030  if(nclus == 0) return kFALSE;
3031 
3032  // Loop over EMCal clusters
3033  Bool_t eventIsAccepted = kFALSE;
3034  for(Int_t i = 0; i < nclus; i++){
3035  AliVCluster* clus = NULL;
3036  clus = event->GetCaloCluster(i);
3037  if (!clus) continue;
3038  if (!clus->IsEMCAL()) continue;
3039  if (clus->GetM02()<0.1) continue;
3040  if (clus->GetNCells()<2) continue;
3041  if (clus->E() > threshold ){
3042 // cout << "found L1G2" << endl;
3043  eventIsAccepted = kTRUE;
3044  }
3045  }
3046  return eventIsAccepted;
3047  }
3048  return kTRUE;
3049  } else if (fSpecialTrigger == 9 ) {
3050  return kTRUE;
3051  } else {
3052  return kTRUE;
3053  }
3054 
3055  return kTRUE;
3056 }
3057 
3058 
3059 //________________________________________________________________________
3061 {
3062 
3063  AliInputEventHandler *fInputHandler=(AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
3064 
3065 
3066  UInt_t isSelected = AliVEvent::kAny;
3067 
3068  if (fInputHandler==NULL) return kFALSE;
3069  if( fInputHandler->GetEventSelection() || event->IsA()==AliAODEvent::Class()) {
3070 
3071  TString firedTrigClass = event->GetFiredTriggerClasses();
3072  // if no trigger has been selected manually, select kAny in case of presel (also important for AOD filtering!)
3073  // in other cases select standards depending on system
3075  if (fPreSelCut) fOfflineTriggerMask = AliVEvent::kAny;
3076  else {
3077  if (fIsHeavyIon == 1){
3078  fOfflineTriggerMask = AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral;
3079  } else if (fIsHeavyIon == 2){
3080  fOfflineTriggerMask = AliVEvent::kINT7;
3081  } else {
3082  fOfflineTriggerMask = AliVEvent::kMB;
3083  }
3084  }
3085  }
3086 
3087  // in case of MC switch to kAny if no MB/INT7/INT8 has been selected
3088  if(isMC){
3089  if( fIsHeavyIon == 0){
3090  if( fOfflineTriggerMask != AliVEvent::kMB && fOfflineTriggerMask != AliVEvent::kINT7 && fOfflineTriggerMask != AliVEvent::kINT8 ){
3091  fOfflineTriggerMask = AliVEvent::kAny;
3092  }
3093  }else{
3094  fOfflineTriggerMask = AliVEvent::kAny;
3095  }
3096  }
3097 
3098  if (fOfflineTriggerMask){
3099  isSelected = fOfflineTriggerMask & fInputHandler->IsEventSelected();
3100  if (isSelected && !fPreSelCut){
3101 // cout << firedTrigClass.Data() << endl;
3102 // cout << "Special trigger: "<< fSpecialTrigger << " initialized " << fEMCALTrigInitialized << endl;
3103 // if (fSpecialTrigger == 5 || fSpecialTrigger == 8 || fSpecialTrigger == 9){ // EMCAL triggers
3104 // if (!fEMCALTrigInitialized ) InitializeEMCALTrigger(event);
3105 // fTriggersEMCAL= GetTriggerList();
3106 // }
3107  if (fSpecialSubTrigger>0 && !isMC){
3108  if (!firedTrigClass.Contains(fSpecialSubTriggerName.Data())) isSelected = 0;
3109  if (fRejectTriggerOverlap){
3110  // trigger rejection EMC1,7,8
3111  if (fSpecialTrigger == 5 && fSpecialSubTriggerName.CompareTo("CEMC7") == 0){
3112  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3113  } else if (fSpecialTrigger == 5 && fSpecialSubTriggerName.CompareTo("CEMC1") == 0){
3114  if (fInputHandler->IsEventSelected() & AliVEvent::kMB) isSelected = 0;
3115  } else if (fSpecialTrigger == 5 && fSpecialSubTriggerName.CompareTo("CEMC8") == 0){
3116  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8) isSelected = 0;
3117  }
3118  // trigger rejection EGA
3119  if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("7EGA") == 0){
3120  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3121  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3122  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("8EGA") == 0){
3123  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8) isSelected = 0;
3124  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3125  }
3126  // trigger rejection EG1 & EG2
3127  if (fPeriodEnum == kLHC13g){
3128  // EG1 is the trigger with the highest threshold
3129  if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("7EG1") == 0){
3130 // cout << firedTrigClass.Data() << endl;
3131  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3132 // cout << "INT7? " << isSelected << endl;
3133  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3134 // cout << "CEM7? " << isSelected << endl;
3135  if (firedTrigClass.Contains("7EG2")) isSelected = 0;
3136 // cout << "7EG2? " << isSelected << endl;
3137  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("8EG1") == 0){
3138  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8) isSelected = 0;
3139  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3140  if (firedTrigClass.Contains("8EG2")) isSelected = 0;
3141  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("7EG2") == 0){
3142 // cout << firedTrigClass.Data() << endl;
3143  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3144 // cout << "INT7? " << isSelected << endl;
3145  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3146 // cout << "CEM7? " << isSelected << endl;
3147  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("8EG2") == 0){
3148  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3149  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3150  }
3151  } else {
3152  // EG2 is the trigger with the highest threshold
3153  if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("7EG2") == 0){
3154  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3155  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3156  if (firedTrigClass.Contains("7EG1")) isSelected = 0;
3157  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("8EG2") == 0){
3158  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8) isSelected = 0;
3159  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3160  if (firedTrigClass.Contains("8EG1")) isSelected = 0;
3161  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("7EG1") == 0){
3162  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3163  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3164  } else if (fSpecialTrigger == 8 && fSpecialSubTriggerName.CompareTo("8EG1") == 0){
3165  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7) isSelected = 0;
3166  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7) isSelected = 0;
3167  }
3168  }
3169  }
3170  if (isSelected != 0 ){
3171 // cout << "I am here" << " :" << fSpecialSubTriggerName.Data() <<endl;
3172  if (fSpecialTrigger == 5 || fSpecialTrigger == 8 || fSpecialTrigger == 9 ){
3174  if (fInputHandler->IsEventSelected() & AliVEvent::kMB)hTriggerClassesCorrelated->Fill(0);
3175  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7)hTriggerClassesCorrelated->Fill(1);
3176  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC1)hTriggerClassesCorrelated->Fill(2);
3177  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7)hTriggerClassesCorrelated->Fill(3);
3178  if (firedTrigClass.Contains("7EJE") || firedTrigClass.Contains("8EJE")) hTriggerClassesCorrelated->Fill(4);
3179  if (firedTrigClass.Contains("7EJ1") || firedTrigClass.Contains("8EJ1")) hTriggerClassesCorrelated->Fill(5);
3180  if (firedTrigClass.Contains("7EJ2") || firedTrigClass.Contains("8EJ2")) hTriggerClassesCorrelated->Fill(6);
3181  if (firedTrigClass.Contains("7EGA") || firedTrigClass.Contains("8EGA")) hTriggerClassesCorrelated->Fill(7);
3182  if (firedTrigClass.Contains("7EG1") || firedTrigClass.Contains("8EG1")) hTriggerClassesCorrelated->Fill(8);
3183  if (firedTrigClass.Contains("7EG2") || firedTrigClass.Contains("8EG2")) hTriggerClassesCorrelated->Fill(9);
3184  }
3185  }
3186  }
3187 
3188  } else if (isMC){
3189  if (fSpecialTrigger == 5 || fSpecialTrigger == 8 || fSpecialTrigger == 9){ // EMCAL triggers
3190 // isSelected = 0;
3191 // if (fTriggersEMCAL > 0)cout << "Special Trigger " << fSpecialTrigger << " triggers: " << fTriggersEMCAL << " selected triggers: " << fTriggersEMCALSelected << " run number: " <<event->GetRunNumber()<<endl;
3192 // if (fTriggersEMCAL&fTriggersEMCALSelected){
3193 // cout << "accepted ++++++++++++++++++++" << endl;
3194  isSelected = 1;
3195 // }
3196  }
3197  }
3198  //if for specific centrality trigger selection
3199  if(fSpecialSubTrigger == 1){
3200  if(fSpecialSubTriggerName.Contains("|") && GetCentrality(event) <= 10.){
3201  TObjArray *ClassesList = fSpecialSubTriggerName.Tokenize("|");
3202  for (Int_t i=0; i<ClassesList->GetEntriesFast();++i){
3203  TObjString *NameClass = (TObjString*)ClassesList->At(i);
3204  if (firedTrigClass.Contains(NameClass->GetString())) isSelected = 1;
3205  }
3206  } else if(fSpecialSubTriggerName.Contains("%")){
3207  TObjArray *ClassesList = fSpecialSubTriggerName.Tokenize("%");
3208  for (Int_t i=0; i<ClassesList->GetEntriesFast();++i){
3209  TObjString *NameClass = (TObjString*)ClassesList->At(i);
3210  if (firedTrigClass.Contains(NameClass->GetString())) isSelected = 1;
3211  }
3212  } else if(fSpecialSubTriggerName.Contains("@")){
3213  TObjArray *ClassesList = fSpecialSubTriggerName.Tokenize("@");
3214  for (Int_t i=0; i<ClassesList->GetEntriesFast();++i){
3215  TObjString *NameClass = (TObjString*)ClassesList->At(i);
3216  if (firedTrigClass.Contains(NameClass->GetString())) isSelected = 1;
3217  }
3218  } else if(fSpecialSubTriggerName.Contains("&")){ //logic AND of two classes
3219  TObjArray *ClassesList = fSpecialSubTriggerName.Tokenize("&");
3220  TString CheckClass = "";
3221  for (Int_t i=0; i<ClassesList->GetEntriesFast(); i++){
3222  TObjString *NameClass = (TObjString*)ClassesList->At(i);
3223  if (firedTrigClass.Contains(NameClass->GetString())) CheckClass+="1";
3224  else CheckClass+="0";
3225  }
3226  if(CheckClass.Contains("0")) isSelected = 0;
3227  }
3228  else if(firedTrigClass.Contains(fSpecialSubTriggerName.Data())) isSelected = 1;
3229  }
3230  }
3231  }
3232  }
3233  fIsSDDFired = !(fInputHandler->IsEventSelected() & AliVEvent::kFastOnly);
3234 
3235  Bool_t mimickedTrigger = kTRUE;
3236  if (fMimicTrigger) mimickedTrigger = MimicTrigger(event, isMC);
3237 // cout << "mimicked decision \t" << mimickedTrigger << "expect decision? "<< fMimicTrigger<< endl;
3238 
3239  // Fill Histogram
3240  if(hTriggerClass){
3241  if (fIsSDDFired) hTriggerClass->Fill(33);
3242  if (mimickedTrigger){
3243  if (fInputHandler->IsEventSelected() & AliVEvent::kMB)hTriggerClass->Fill(0);
3244  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7)hTriggerClass->Fill(1);
3245  if (fInputHandler->IsEventSelected() & AliVEvent::kMUON)hTriggerClass->Fill(2);
3246  if (fInputHandler->IsEventSelected() & AliVEvent::kHighMult)hTriggerClass->Fill(3);
3247  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC1)hTriggerClass->Fill(4);
3248  if (fInputHandler->IsEventSelected() & AliVEvent::kCINT5)hTriggerClass->Fill(5);
3249  if (fInputHandler->IsEventSelected() & AliVEvent::kCMUS5)hTriggerClass->Fill(6);
3250  // if (fInputHandler->IsEventSelected() & AliVEvent::kMUSPB)hTriggerClass->Fill(6);
3251  if (fInputHandler->IsEventSelected() & AliVEvent::kMUSH7)hTriggerClass->Fill(7);
3252  // if (fInputHandler->IsEventSelected() & AliVEvent::kMUSHPB)hTriggerClass->Fill(7);
3253  if (fInputHandler->IsEventSelected() & AliVEvent::kMUL7)hTriggerClass->Fill(8);
3254  // if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikePB)hTriggerClass->Fill(8);
3255  if (fInputHandler->IsEventSelected() & AliVEvent::kMUU7)hTriggerClass->Fill(9);
3256  // if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikePB)hTriggerClass->Fill(9);
3257  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7)hTriggerClass->Fill(10);
3258  // if (fInputHandler->IsEventSelected() & AliVEvent::kEMC8)hTriggerClass->Fill(10);
3259  if (fInputHandler->IsEventSelected() & AliVEvent::kMUS7)hTriggerClass->Fill(11);
3260  if (fInputHandler->IsEventSelected() & AliVEvent::kPHI1)hTriggerClass->Fill(12);
3261  if (fInputHandler->IsEventSelected() & AliVEvent::kPHI7)hTriggerClass->Fill(13);
3262  // if (fInputHandler->IsEventSelected() & AliVEvent::kPHI8)hTriggerClass->Fill(13);
3263  // if (fInputHandler->IsEventSelected() & AliVEvent::kPHOSPb)hTriggerClass->Fill(13);
3264  if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEJE)hTriggerClass->Fill(14);
3265  if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEGA)hTriggerClass->Fill(15);
3266  if (fInputHandler->IsEventSelected() & AliVEvent::kCentral)hTriggerClass->Fill(16);
3267  if (fInputHandler->IsEventSelected() & AliVEvent::kSemiCentral)hTriggerClass->Fill(17);
3268  if (fInputHandler->IsEventSelected() & AliVEvent::kDG5)hTriggerClass->Fill(18);
3269  if (fInputHandler->IsEventSelected() & AliVEvent::kZED)hTriggerClass->Fill(19);
3270  if (fInputHandler->IsEventSelected() & AliVEvent::kSPI7)hTriggerClass->Fill(20);
3271  // if (fInputHandler->IsEventSelected() & AliVEvent::kSPI)hTriggerClass->Fill(20);
3272  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8)hTriggerClass->Fill(21);
3273  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleLowPt8)hTriggerClass->Fill(22);
3274  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleHighPt8)hTriggerClass->Fill(23);
3275  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikeLowPt8)hTriggerClass->Fill(24);
3276  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt8)hTriggerClass->Fill(25);
3277  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt0)hTriggerClass->Fill(26);
3278  if (fInputHandler->IsEventSelected() & AliVEvent::kUserDefined)hTriggerClass->Fill(27);
3279  if (fInputHandler->IsEventSelected() & AliVEvent::kTRD)hTriggerClass->Fill(28);
3280  if (fInputHandler->IsEventSelected() & AliVEvent::kFastOnly)hTriggerClass->Fill(29);
3281  if (fInputHandler->IsEventSelected() & AliVEvent::kAnyINT)hTriggerClass->Fill(30);
3282  if (fInputHandler->IsEventSelected() & AliVEvent::kAny)hTriggerClass->Fill(31);
3283  if (!fInputHandler->IsEventSelected()) hTriggerClass->Fill(34);
3284  }
3285  if (mimickedTrigger && fMimicTrigger) hTriggerClass->Fill(35);
3286  }
3287 
3288  if(hTriggerClassSelected && isSelected){
3289  if (mimickedTrigger){
3290  if (!fIsSDDFired) hTriggerClassSelected->Fill(33);
3291  if (fInputHandler->IsEventSelected() & AliVEvent::kMB)hTriggerClassSelected->Fill(0);
3292  if (fInputHandler->IsEventSelected() & AliVEvent::kINT7)hTriggerClassSelected->Fill(1);
3293  if (fInputHandler->IsEventSelected() & AliVEvent::kMUON)hTriggerClassSelected->Fill(2);
3294  if (fInputHandler->IsEventSelected() & AliVEvent::kHighMult)hTriggerClassSelected->Fill(3);
3295  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC1)hTriggerClassSelected->Fill(4);
3296  if (fInputHandler->IsEventSelected() & AliVEvent::kCINT5)hTriggerClassSelected->Fill(5);
3297  if (fInputHandler->IsEventSelected() & AliVEvent::kCMUS5)hTriggerClassSelected->Fill(6);
3298  // if (fInputHandler->IsEventSelected() & AliVEvent::kMUSPB)hTriggerClassSelected->Fill(6);
3299  if (fInputHandler->IsEventSelected() & AliVEvent::kMUSH7)hTriggerClassSelected->Fill(7);
3300  // if (fInputHandler->IsEventSelected() & AliVEvent::kMUSHPB)hTriggerClassSelected->Fill(7);
3301  if (fInputHandler->IsEventSelected() & AliVEvent::kMUL7)hTriggerClassSelected->Fill(8);
3302  // if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikePB)hTriggerClassSelected->Fill(8);
3303  if (fInputHandler->IsEventSelected() & AliVEvent::kMUU7)hTriggerClassSelected->Fill(9);
3304  // if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikePB)hTriggerClassSelected->Fill(9);
3305  if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7)hTriggerClassSelected->Fill(10);
3306  // if (fInputHandler->IsEventSelected() & AliVEvent::kEMC8)hTriggerClassSelected->Fill(10);
3307  if (fInputHandler->IsEventSelected() & AliVEvent::kMUS7)hTriggerClassSelected->Fill(11);
3308  if (fInputHandler->IsEventSelected() & AliVEvent::kPHI1)hTriggerClassSelected->Fill(12);
3309  if (fInputHandler->IsEventSelected() & AliVEvent::kPHI7)hTriggerClassSelected->Fill(13);
3310  // if (fInputHandler->IsEventSelected() & AliVEvent::kPHI8)hTriggerClassSelected->Fill(13);
3311  // if (fInputHandler->IsEventSelected() & AliVEvent::kPHOSPb)hTriggerClassSelected->Fill(13);
3312  if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEJE)hTriggerClassSelected->Fill(14);
3313  if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEGA)hTriggerClassSelected->Fill(15);
3314  if (fInputHandler->IsEventSelected() & AliVEvent::kCentral)hTriggerClassSelected->Fill(16);
3315  if (fInputHandler->IsEventSelected() & AliVEvent::kSemiCentral)hTriggerClassSelected->Fill(17);
3316  if (fInputHandler->IsEventSelected() & AliVEvent::kDG5)hTriggerClassSelected->Fill(18);
3317  if (fInputHandler->IsEventSelected() & AliVEvent::kZED)hTriggerClassSelected->Fill(19);
3318  if (fInputHandler->IsEventSelected() & AliVEvent::kSPI7)hTriggerClassSelected->Fill(20);
3319  // if (fInputHandler->IsEventSelected() & AliVEvent::kSPI)hTriggerClassSelected->Fill(20);
3320  if (fInputHandler->IsEventSelected() & AliVEvent::kINT8)hTriggerClassSelected->Fill(21);
3321  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleLowPt8)hTriggerClassSelected->Fill(22);
3322  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleHighPt8)hTriggerClassSelected->Fill(23);
3323  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikeLowPt8)hTriggerClassSelected->Fill(24);
3324  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt8)hTriggerClassSelected->Fill(25);
3325  if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt0)hTriggerClassSelected->Fill(26);
3326  if (fInputHandler->IsEventSelected() & AliVEvent::kUserDefined)hTriggerClassSelected->Fill(27);
3327  if (fInputHandler->IsEventSelected() & AliVEvent::kTRD)hTriggerClassSelected->Fill(28);
3328  if (fInputHandler->IsEventSelected() & AliVEvent::kFastOnly)hTriggerClassSelected->Fill(29);
3329  if (fInputHandler->IsEventSelected() & AliVEvent::kAnyINT)hTriggerClassSelected->Fill(30);
3330  if (fInputHandler->IsEventSelected() & AliVEvent::kAny)hTriggerClassSelected->Fill(31);
3331  }
3332  if (mimickedTrigger && fMimicTrigger) hTriggerClassSelected->Fill(34);
3333  }
3334 
3335  if(!isSelected)return kFALSE;
3336  if (fMimicTrigger)
3337  if (!mimickedTrigger ) return kFALSE;
3338  return kTRUE;
3339 
3340 }
3341 
3342 //________________________________________________________________________
3344  // returns TString with current cut number
3345  return fCutStringRead;
3346 }
3347 
3348 //________________________________________________________________________
3349 void AliConvEventCuts::GetNotRejectedParticles(Int_t rejection, TList *HeaderList, AliVEvent *event){
3350 
3351  if(fNotRejectedStart){
3352  delete[] fNotRejectedStart;
3353  fNotRejectedStart = NULL;
3354  }
3355  if(fNotRejectedEnd){
3356  delete[] fNotRejectedEnd;
3357  fNotRejectedEnd = NULL;
3358  }
3359  if(fGeneratorNames){
3360  delete[] fGeneratorNames;
3361  fGeneratorNames = NULL;
3362  }
3363 
3364  if(rejection == 0) return; // No Rejection
3365 
3366  AliGenCocktailEventHeader *cHeader = 0x0;
3367  AliAODMCHeader *cHeaderAOD = 0x0;
3368  Bool_t headerFound = kFALSE;
3369  AliMCEvent *fMCEvent = 0x0;
3370  TClonesArray *fMCEventAOD = 0x0;
3371  if(event->IsA()==AliMCEvent::Class()){
3372  if(dynamic_cast<AliMCEvent*>(event)){
3373  cHeader = dynamic_cast<AliGenCocktailEventHeader*>(dynamic_cast<AliMCEvent*>(event)->GenEventHeader());
3374  fMCEvent = dynamic_cast<AliMCEvent*>(event);
3375  if(cHeader) headerFound = kTRUE;
3376  }
3377  }
3378  if(event->IsA()==AliAODEvent::Class()){ // event is a AODEvent in case of AOD
3379  cHeaderAOD = dynamic_cast<AliAODMCHeader*>(event->FindListObject(AliAODMCHeader::StdBranchName()));
3380  fMCEventAOD = dynamic_cast<TClonesArray*>(event->FindListObject(AliAODMCParticle::StdBranchName()));
3381  if(cHeaderAOD) headerFound = kTRUE;
3382  }
3383 
3384  if (fDebugLevel > 0 ) cout << "event starts here" << endl;
3385  if(headerFound){
3386  TList *genHeaders = 0x0;
3387  if(cHeader) genHeaders = cHeader->GetHeaders();
3388  if(cHeaderAOD){
3389  genHeaders = cHeaderAOD->GetCocktailHeaders();
3390  if(genHeaders->GetEntries()==1){
3392  return;
3393  }
3394  }
3395  AliGenEventHeader* gh = 0;
3396  fnHeaders = 0;
3397  Int_t firstindexA = 0;
3398  Int_t lastindexA = -1;
3399  if(rejection == 1 || rejection == 3) fnHeaders = 1; // MinBiasHeader
3400  if(rejection == 2){ // TList of Headers Names
3401  for(Int_t i = 0; i<genHeaders->GetEntries();i++){
3402  gh = (AliGenEventHeader*)genHeaders->At(i);
3403  TString GeneratorName = gh->GetName();
3404  lastindexA = lastindexA + gh->NProduced();
3405  if (fDebugLevel > 0 ) cout << i << "\t" << GeneratorName.Data() << endl;
3406  for(Int_t j = 0; j<HeaderList->GetEntries();j++){
3407  TString GeneratorInList = ((TObjString*)HeaderList->At(j))->GetString();
3408  if (fDebugLevel > 0 ) cout << GeneratorInList.Data() << endl;
3409  if(GeneratorName.CompareTo(GeneratorInList) == 0){
3410  if (fDebugLevel > 0 ) cout << "accepted" << endl;
3411  if (GeneratorInList.CompareTo("PARAM") == 0 || GeneratorInList.CompareTo("BOX") == 0 ){
3412  if(fMCEvent){
3413  if (fPeriodEnum == kLHC14a1b || fPeriodEnum == kLHC14a1c ){
3414  if (fMCEvent->Particle(firstindexA)->GetPdgCode() == fAddedSignalPDGCode ) {
3415  if (gh->NProduced() > 10 && fMCEvent->Particle(firstindexA+10)->GetPdgCode() == fAddedSignalPDGCode ){
3416  if (fDebugLevel > 0 ) cout << "cond 1: "<< fnHeaders << endl;
3417  fnHeaders++;
3418  continue;
3419  }
3420  continue;
3421  }
3422  } else {
3423  if (fDebugLevel > 0 ) cout << "cond 2: " << fnHeaders << endl;
3424  fnHeaders++;
3425  continue;
3426 
3427  }
3428  }
3429  if ( fMCEventAOD){
3430  AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(fMCEventAOD->At(firstindexA));
3431  if (aodMCParticle && (fPeriodEnum == kLHC14a1b || fPeriodEnum == kLHC14a1c) ){
3432  if ( aodMCParticle->GetPdgCode() == fAddedSignalPDGCode ){
3433  if (gh->NProduced() > 10){
3434  AliAODMCParticle *aodMCParticle2 = static_cast<AliAODMCParticle*>(fMCEventAOD->At(firstindexA+10));
3435  if ( aodMCParticle2->GetPdgCode() == fAddedSignalPDGCode ){
3436  if (fDebugLevel > 0 ) cout << "cond 1: " << fnHeaders << endl;
3437  fnHeaders++;
3438  continue;
3439  }
3440  }
3441  continue;
3442  }
3443  } else {
3444  if (fDebugLevel > 0 ) cout << "cond 2: " << fnHeaders << endl;
3445  fnHeaders++;
3446  continue;
3447  }
3448  }
3449  continue;
3450  }
3451  if (fDebugLevel > 0 ) cout << "cond 3: "<< fnHeaders << endl;
3452  fnHeaders++;
3453  continue;
3454  }
3455  }
3456  firstindexA = firstindexA + gh->NProduced();
3457  }
3458  }
3459  if (fDebugLevel > 0 ) cout << "number of headers: " <<fnHeaders << endl;
3460 
3464 
3465  if(rejection == 1 || rejection == 3){
3466  fNotRejectedStart[0] = 0;
3467  fNotRejectedEnd[0] = ((AliGenEventHeader*)genHeaders->At(0))->NProduced()-1;
3468  fGeneratorNames[0] = ((AliGenEventHeader*)genHeaders->At(0))->GetName();
3469  if (fDebugLevel > 0 ) cout << 0 << "\t" <<fGeneratorNames[0] << "\t" << fNotRejectedStart[0] << "\t" <<fNotRejectedEnd[0] << endl;
3470  return;
3471  }
3472 
3473  Int_t firstindex = 0;
3474  Int_t lastindex = -1;
3475  Int_t number = 0;
3476 
3477  for(Int_t i = 0; i<genHeaders->GetEntries();i++){
3478  gh = (AliGenEventHeader*)genHeaders->At(i);
3479  TString GeneratorName = gh->GetName();
3480  lastindex = lastindex + gh->NProduced();
3481  for(Int_t j = 0; j<HeaderList->GetEntries();j++){
3482  TString GeneratorInList = ((TObjString*)HeaderList->At(j))->GetString();
3483  if (fDebugLevel > 0 ) cout << i << "\t" << GeneratorName.Data() << endl;
3484  if(GeneratorName.CompareTo(GeneratorInList) == 0){
3485  if (GeneratorInList.CompareTo("PARAM") == 0 || GeneratorInList.CompareTo("BOX") == 0 ){
3486  if(fMCEvent){
3487  if (fPeriodEnum == kLHC14a1b || fPeriodEnum == kLHC14a1c ){
3488  if (fMCEvent->Particle(firstindex)->GetPdgCode() == fAddedSignalPDGCode ) {
3489  if (fDebugLevel > 0 ) cout << "produced " << gh->NProduced() << " with box generator" << endl;
3490  if (gh->NProduced() > 10 && fMCEvent->Particle(firstindex+10)->GetPdgCode() == fAddedSignalPDGCode){
3491  if (fDebugLevel > 0 ) cout << "one of them was a pi0 or eta" << endl;
3492  fNotRejectedStart[number] = firstindex;
3493  fNotRejectedEnd[number] = lastindex;
3494  fGeneratorNames[number] = GeneratorName;
3495  number++;
3496  if (fDebugLevel > 0 ) cout << "Number of particles produced for: " << i << "\t" << GeneratorName.Data() << "\t" << lastindex-firstindex+1 << endl;
3497  continue;
3498  }
3499  }
3500  } else {
3501  fNotRejectedStart[number] = firstindex;
3502  fNotRejectedEnd[number] = lastindex;
3503  fGeneratorNames[number] = GeneratorName;
3504  number++;
3505  continue;
3506  }
3507  }
3508  if ( fMCEventAOD){
3509  AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(fMCEventAOD->At(firstindex));
3510  if (fPeriodEnum == kLHC14a1b || fPeriodEnum == kLHC14a1c ){
3511  if ( aodMCParticle->GetPdgCode() == fAddedSignalPDGCode ){
3512  if (gh->NProduced() > 10) {
3513  AliAODMCParticle *aodMCParticle2 = static_cast<AliAODMCParticle*>(fMCEventAOD->At(firstindex+10));
3514  if ( aodMCParticle2->GetPdgCode() == fAddedSignalPDGCode ){
3515  fNotRejectedEnd[number] = lastindex;
3516  fNotRejectedStart[number] = firstindex;
3517  fGeneratorNames[number] = GeneratorName;
3518  number++;
3519  }
3520  continue;
3521  }
3522  }
3523  } else {
3524  fNotRejectedStart[number] = firstindex;
3525  fNotRejectedEnd[number] = lastindex;
3526  fGeneratorNames[number] = GeneratorName;
3527  number++;
3528  continue;
3529  }
3530  }
3531  continue;
3532  } else {
3533  fNotRejectedStart[number] = firstindex;
3534  fNotRejectedEnd[number] = lastindex;
3535  fGeneratorNames[number] = GeneratorName;
3536  if (fDebugLevel > 0 ) cout << "Number of particles produced for: " << i << "\t" << GeneratorName.Data() << "\t" << lastindex-firstindex+1 << endl;
3537  number++;
3538  continue;
3539  }
3540 
3541  }
3542  }
3543  firstindex = firstindex + gh->NProduced();
3544  }
3545  if (fDebugLevel > 0 ) {
3546  for (Int_t i = 0; i < number; i++){
3547  cout << i << "\t" <<fGeneratorNames[i] << "\t" << fNotRejectedStart[i] << "\t" <<fNotRejectedEnd[i] << endl;
3548  }
3549  }
3550  } else { // No Cocktail Header Found
3551  fNotRejectedStart = new Int_t[1];
3552  fNotRejectedEnd = new Int_t[1];
3553 
3554  fnHeaders = 1;
3555  fNotRejectedStart[0] = 0;
3556  fNotRejectedEnd[0] = static_cast<AliMCEvent*>(event)->GetNumberOfPrimaries()-1;
3557  if (rejection > 1){
3558  fNotRejectedStart[0] = -1;
3559  fNotRejectedEnd[0] = -1;
3560  }
3561 
3562  fGeneratorNames = new TString[1];
3563  fGeneratorNames[0] = "NoCocktailGeneratorFound";
3564 // SetRejectExtraSignalsCut(0);
3565  }
3566 
3567 }
3568 
3569 //_________________________________________________________________________
3570 Int_t AliConvEventCuts::IsParticleFromBGEvent(Int_t index, AliMCEvent *mcEvent, AliVEvent *InputEvent){
3571 
3572  if (fDebugLevel > 2 ) cout << index << endl;
3573  if(index < 0) return 0; // No Particle
3574 
3575  Int_t accepted = 0;
3576  if(!InputEvent || InputEvent->IsA()==AliESDEvent::Class()){
3577  if(!mcEvent) return 0; // no mcEvent available, return 0
3578  if(index >= mcEvent->GetNumberOfPrimaries()){ // initial particle is secondary particle
3579  if( ((TParticle*)mcEvent->Particle(index))->GetMother(0) < 0) return 0; // material particle, return 0
3580  return IsParticleFromBGEvent(((TParticle*)mcEvent->Particle(index))->GetMother(0),mcEvent,InputEvent);
3581  }
3582  for(Int_t i = 0;i<fnHeaders;i++){
3583  if (fDebugLevel > 2 ) cout << "header " << i << ":"<< fNotRejectedStart[i] << "\t" << fNotRejectedEnd[i] << endl;
3584  if(index >= fNotRejectedStart[i] && index <= fNotRejectedEnd[i]){
3585  if (fDebugLevel > 1 ) cout << "accepted:" << index << "\t header " << i << ": "<< fNotRejectedStart[i] << "\t" << fNotRejectedEnd[i] << endl;
3586  accepted = 1;
3587  if(i == 0) accepted = 2; // MB Header
3588  }
3589  }
3590  }
3591  else if(InputEvent->IsA()==AliAODEvent::Class()){
3592  TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(InputEvent->FindListObject(AliAODMCParticle::StdBranchName()));
3593  if (AODMCTrackArray){
3594  AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index));
3595  if(!aodMCParticle) return 0; // no particle
3596  if(!aodMCParticle->IsPrimary()){
3597  if( aodMCParticle->GetMother() < 0) return 0;// material particle, return 0
3598  return IsParticleFromBGEvent(aodMCParticle->GetMother(),mcEvent,InputEvent);
3599  }
3600  index = TMath::Abs(static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index))->GetLabel());
3601  for(Int_t i = 0;i<fnHeaders;i++){
3602  if(index >= fNotRejectedStart[i] && index <= fNotRejectedEnd[i]){
3603  accepted = 1;
3604  if(i == 0) accepted = 2; // MB Header
3605  }
3606  }
3607  }
3608  }
3609 
3610  return accepted;
3611 }
3612 
3613 //_________________________________________________________________________
3614 Int_t AliConvEventCuts::IsEventAcceptedByCut(AliConvEventCuts *ReaderCuts, AliVEvent *event, AliMCEvent *mcEvent, Int_t isHeavyIon, Bool_t isEMCALAnalysis){
3615 
3616  Bool_t isMC = kFALSE;
3617  if (mcEvent){isMC = kTRUE;}
3618 
3619  if ( !IsTriggerSelected(event, isMC) )
3620  return 3;
3621 
3622  if( !(IsCentralitySelected(event,mcEvent)))
3623  return 1; // Check Centrality --> Not Accepted => eventQuality = 1
3624 
3625  Bool_t hasV0And = ReaderCuts->HasV0AND();
3626  Bool_t isSDDFired = ReaderCuts->IsSDDFired();
3627 
3628  if( ( (IsSpecialTrigger() == 0 && IsSpecialSubTrigger() == 1) || (IsSpecialTrigger() == 1 && IsSpecialSubTrigger() == 1) ) && !isSDDFired && !mcEvent)
3629  //if V0OR with SDD requested or V0AND with SDD request but the SDD has not fired
3630  return 7; // V0 with SDD requested but no fired
3631 
3632  if( ( (IsSpecialTrigger() == 1 && IsSpecialSubTrigger() == 0) || (IsSpecialTrigger() == 1 && IsSpecialSubTrigger() == 1) ) && !hasV0And)
3633  //if V0AND (only) or V0AND with SDD requested but V0AND requested but no fired
3634  return 8; // V0AND requested but no fired
3635 
3636 
3637  if( (IsSpecialTrigger() == 2 || IsSpecialTrigger() == 3) && !isSDDFired && !mcEvent)
3638  return 7; // With SDD requested but no fired
3639 
3640  if( (IsSpecialTrigger() == 1 || IsSpecialTrigger() == 3) && !hasV0And)
3641  return 8; // V0AND requested but no fired
3642 
3643  // Special EMCAL checks due to hardware issues in LHC11a
3644  if (isEMCALAnalysis || IsSpecialTrigger() == 5 || IsSpecialTrigger() == 8 || IsSpecialTrigger() == 9 ){
3645  Int_t runnumber = event->GetRunNumber();
3646  if ((runnumber>=144871) && (runnumber<=146860)) {
3647 
3648  AliVCaloCells *cells = event->GetEMCALCells();
3649  const Short_t nCells = cells->GetNumberOfCells();
3650 
3651  if (event->IsA()==AliESDEvent::Class()) AliAnalysisManager::GetAnalysisManager()->LoadBranch("EMCALCells.");
3652 
3653  AliInputEventHandler *fInputHandler=(AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
3654  if (!fInputHandler) return 3;
3655 
3656  // count cells above threshold
3657  Int_t nCellCount[12] = {0,0,0,0,0,0,0,0,0,0,0,0};
3658  for(Int_t iCell=0; iCell<nCells; ++iCell) {
3659  Short_t cellId = cells->GetCellNumber(iCell);
3660  Double_t cellE = cells->GetCellAmplitude(cellId);
3661  Int_t sm = cellId / (24*48);
3662  if (cellE>0.1) ++nCellCount[sm];
3663  }
3664 
3665  Bool_t fIsLedEvent = kFALSE;
3666  if (nCellCount[4] > 100) {
3667  fIsLedEvent = kTRUE;
3668  } else {
3669  if ((runnumber>=146858) && (runnumber<=146860)) {
3670  if ((fInputHandler->IsEventSelected() & AliVEvent::kMB) && (nCellCount[3]>=21))
3671  fIsLedEvent = kTRUE;
3672  else if ((fInputHandler->IsEventSelected() & AliVEvent::kEMC1) && (nCellCount[3]>=35))
3673  fIsLedEvent = kTRUE;
3674  }
3675  }
3676  if (fIsLedEvent) {
3677  return 9;
3678  }
3679  }
3680  }
3681 
3682  // SPD clusters vs tracklets to check for pileup/background
3683  Int_t nClustersLayer0 = event->GetNumberOfITSClusters(0);
3684  Int_t nClustersLayer1 = event->GetNumberOfITSClusters(1);
3685  Int_t nTracklets = event->GetMultiplicity()->GetNumberOfTracklets();
3686  if(hSPDClusterTrackletBackgroundBefore) hSPDClusterTrackletBackgroundBefore->Fill(nTracklets, (nClustersLayer0 + nClustersLayer1));
3687 
3688 
3689  Double_t distZMax = 0;
3690  if(event->IsA()==AliESDEvent::Class()){
3691  Int_t nPileVert = ((AliESDEvent*)event)->GetNumberOfPileupVerticesSPD();
3692  if (hNPileupVertices) hNPileupVertices->Fill(nPileVert);
3693  if (nPileVert > 0){
3694  for(Int_t i=0; i<nPileVert;i++){
3695  const AliESDVertex* pv= ((AliESDEvent*)event)->GetPileupVertexSPD(i);
3696  Int_t nc2 = pv->GetNContributors();
3697  if(nc2>=3){
3698  Double_t z1 = ((AliESDEvent*)event)->GetPrimaryVertexSPD()->GetZ();
3699  Double_t z2 = pv->GetZ();
3700  Double_t distZ = z2-z1;
3701  if (TMath::Abs(distZMax) < TMath::Abs(distZ) ){
3702  distZMax = distZ;
3703  }
3704  }
3705  }
3706  if (hPileupVertexToPrimZ) hPileupVertexToPrimZ->Fill(distZMax);
3707  }
3708  }
3709  if(GetPastFutureLowBC()!=0 && GetPastFutureHighBC()!=0 ){
3710  if(IsOutOfBunchPileupPastFuture(event))
3711  return 12;
3712  }
3713 
3714  if( isHeavyIon != 2 && GetIsFromPileup()){
3715  if(event->IsPileupFromSPD(3,0.8,3.,2.,5.) ){
3717  return 6; // Check Pileup --> Not Accepted => eventQuality = 6
3718  }
3719  if (fUtils->IsSPDClusterVsTrackletBG(event)){
3721  return 11; // Check Pileup --> Not Accepted => eventQuality = 11
3722  }
3723  }
3724  if(isHeavyIon == 2 && GetIsFromPileup()){
3725  if(fUtils->IsPileUpEvent(event) ){
3727  return 6; // Check Pileup --> Not Accepted => eventQuality = 6
3728  }
3729  if (fUtils->IsSPDClusterVsTrackletBG(event)){
3731  return 11; // Check Pileup --> Not Accepted => eventQuality = 11
3732  }
3733  }
3734 
3736  if( IsPileUpV0MTPCout(event) ){
3737  return 13;
3738  }
3739  }
3740 
3741  if(hCentrality)hCentrality->Fill(GetCentrality(event));
3742 
3743  if(hVertexZ)hVertexZ->Fill(event->GetPrimaryVertex()->GetZ());
3744 // if(hCentralityVsNumberOfPrimaryTracks)
3745 // hCentralityVsNumberOfPrimaryTracks->Fill(GetCentrality(event),
3746 // ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()
3747 // ->GetTask(fV0ReaderName.Data()))->GetNumberOfPrimaryTracks());
3748 
3749  if(fIsHeavyIon == 1){
3750  AliEventplane *EventPlane = event->GetEventplane();
3751  fEventPlaneAngle = EventPlane->GetEventplane("V0",event,2);
3752  if(hEventPlaneAngle)hEventPlaneAngle->Fill(TMath::Abs(fEventPlaneAngle));
3753  }
3754  if(hSPDClusterTrackletBackground) hSPDClusterTrackletBackground->Fill(nTracklets, (nClustersLayer0 + nClustersLayer1));
3755 
3756  return 0;
3757 }
3758 
3759 
3760 //_________________________________________________________________________
3762 
3763  AliInfo("Inside the GetWeightForCentralityFlattening function");
3764  Double_t centrality = 0.;
3765  //obtain centrality for ESD or AOD
3766  if(!event || event->IsA()==AliESDEvent::Class()){
3767  AliESDEvent *esdEvent=dynamic_cast<AliESDEvent*>(event);
3768  if(esdEvent){
3769  AliCentrality *fESDCentrality=(AliCentrality*)esdEvent->GetCentrality();
3770  if(fDetectorCentrality==0 && fIsHeavyIon==1){
3771  centrality = fESDCentrality->GetCentralityPercentile("V0M"); // default for PbPb
3772  }
3773  }
3774  } else if(event->IsA()==AliAODEvent::Class()){
3775  AliAODEvent *aodEvent=dynamic_cast<AliAODEvent*>(event);
3776  if(aodEvent){
3777  if(aodEvent->GetHeader()){
3778  centrality = ((AliVAODHeader*)aodEvent->GetHeader())->GetCentrality();
3779  }
3780  }
3781  }
3782 
3783  //Get the maximum vlaue from the reference distribution and interpolated value
3784  Float_t GetValueForWeight = 1.;
3785  Float_t maximum = 1.;
3786  Double_t weightCentrality = 1.;
3787  Bool_t CorrCentrLoop = kFALSE;
3788 
3789  //depending on the value of the flag, flattening in different cent. range
3790  if ( fDoCentralityFlat == 1 && (centrality >= 0. && centrality <= 10.) ){
3791  GetValueForWeight = hCentralityNotFlat->Interpolate(centrality);
3792  maximum = hCentralityNotFlat->GetMaximum();
3793  CorrCentrLoop = kTRUE;
3794  } else if ( fDoCentralityFlat == 2 && (centrality >=10. && centrality <= 20.) ){
3795  GetValueForWeight = hCentralityNotFlat->Interpolate(centrality);
3796  maximum = hCentralityNotFlat->GetMaximum();
3797  CorrCentrLoop = kTRUE;
3798  } else if ( fDoCentralityFlat == 8 ){
3799  GetValueForWeight = hCentralityNotFlat->Interpolate(centrality);
3800  maximum = hCentralityNotFlat->GetMaximum();
3801  CorrCentrLoop = kTRUE;
3802  } else {
3803  CorrCentrLoop = kFALSE;
3804  }
3805 
3806  if (CorrCentrLoop && GetValueForWeight != 0. && maximum !=0. && isfinite(GetValueForWeight) && isfinite(maximum) ){
3807  weightCentrality = maximum/GetValueForWeight;
3808  if (!isfinite(GetValueForWeight)) weightCentrality = 1.;
3809  if (!isfinite(weightCentrality)) weightCentrality = 1.;
3810  }
3811 
3812  return weightCentrality;
3813 }
3814 
3815 //_________________________________________________________________________
3817 
3818  Double_t weightMult = 1.;
3819 
3820  Float_t valueMultData = -1.;
3821  Float_t valueMultMC = -1.;
3822 
3823  if (hReweightMultData == NULL || hReweightMultMC == NULL ) return weightMult;
3824 
3825  valueMultData = hReweightMultData->Interpolate(mult);
3826  valueMultMC = hReweightMultMC->Interpolate(mult);
3827 
3828  Float_t relativeErrorMC = hReweightMultMC->GetBinError(hReweightMultMC->FindBin(mult))/hReweightMultMC->GetBinContent(hReweightMultMC->FindBin(mult));
3829  Float_t relativeErrorData = hReweightMultData->GetBinError(hReweightMultData->FindBin(mult))/hReweightMultData->GetBinContent(hReweightMultData->FindBin(mult));
3830 
3831  if (relativeErrorData < 0.2 && relativeErrorMC < 0.2 ){
3832  if (isfinite(valueMultData) && isfinite(valueMultMC) ){
3833  weightMult = valueMultData/valueMultMC;
3834  }
3835  }
3836 
3837  return weightMult;
3838 }
3839 
3840 
3841 
3842 //_________________________________________________________________________
3843 Float_t AliConvEventCuts::GetWeightForMeson(Int_t index, AliMCEvent *mcEvent, AliVEvent *event){
3844  if (!( fPeriodEnum == kLHC13d2 || fPeriodEnum == kLHC13d2b || // LHC10h MCs
3845  fPeriodEnum == kLHC14a1a || fPeriodEnum == kLHC14a1b || fPeriodEnum == kLHC14a1c || // LHC11h MCs
3846  fPeriodEnum == kLHC13e7 || fPeriodEnum == kLHC13b2_efix || fPeriodEnum == kLHC14b2 || // LHC13bc MCs
3847  fPeriodEnum == kLHC14e2a || fPeriodEnum == kLHC14e2b || fPeriodEnum == kLHC14e2c || // LHC12[a-i] pass 1 MCs
3848  fPeriodEnum == kLHC12f1a || fPeriodEnum == kLHC12f1b || fPeriodEnum == kLHC12i3 // LHC11a MCs
3849  ) ) return 1.;
3850  Int_t kCaseGen = 0;
3851 
3852  if(index < 0) return 0; // No Particle
3853 
3854  if (IsParticleFromBGEvent(index, mcEvent, event)){
3857  kCaseGen = 1;
3858  }
3859  }
3860  if (kCaseGen == 0) return 1;
3861 
3862  Double_t mesonPt = 0;
3863  //Double_t mesonMass = 0;
3864  Int_t PDGCode = 0;
3865  if(!event || event->IsA()==AliESDEvent::Class()){
3866  mesonPt = ((TParticle*)mcEvent->Particle(index))->Pt();
3867  //mesonMass = ((TParticle*)mcEvent->Particle(index))->GetCalcMass();
3868  PDGCode = ((TParticle*)mcEvent->Particle(index))->GetPdgCode();
3869  } else if(event->IsA()==AliAODEvent::Class()){
3870  TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(event->FindListObject(AliAODMCParticle::StdBranchName()));
3871  if (AODMCTrackArray){
3872  AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index));
3873  mesonPt = aodMCParticle->Pt();
3874  //mesonMass = aodMCParticle->GetCalcMass();
3875  PDGCode = aodMCParticle->GetPdgCode();
3876  } else {
3877  return 1;
3878  }
3879  }
3880 
3881  Float_t functionResultMC = 1.;
3882  if ( PDGCode == 111 && fDoReweightHistoMCPi0 && hReweightMCHistPi0!= 0x0){
3883  functionResultMC = hReweightMCHistPi0->Interpolate(mesonPt);
3884  }
3885  if ( PDGCode == 221 && fDoReweightHistoMCEta && hReweightMCHistEta!= 0x0){
3886  functionResultMC = hReweightMCHistEta->Interpolate(mesonPt);
3887  }
3888  if ( PDGCode == 310 && fDoReweightHistoMCK0s && hReweightMCHistK0s!= 0x0){
3889  functionResultMC = hReweightMCHistK0s->Interpolate(mesonPt);
3890  }
3891 
3892  Float_t functionResultData = 1;
3893  if ( PDGCode == 111 && fDoReweightHistoMCPi0 && fFitDataPi0!= 0x0){
3894  functionResultData = fFitDataPi0->Eval(mesonPt);
3895  }
3896  if ( PDGCode == 221 && fDoReweightHistoMCEta && fFitDataEta!= 0x0){
3897  functionResultData = fFitDataEta->Eval(mesonPt);
3898  }
3899  if ( PDGCode == 310 && fDoReweightHistoMCK0s && fFitDataK0s!= 0x0){
3900  functionResultData = fFitDataK0s->Eval(mesonPt);
3901  }
3902 
3903  Double_t weight = 1;
3904  if (PDGCode == 111 || PDGCode == 221){
3905  if (functionResultData != 0. && functionResultMC != 0. && isfinite(functionResultData) && isfinite(functionResultMC)){
3906  weight = functionResultData/functionResultMC;
3907  if ( kCaseGen == 3){
3908  if (PDGCode == 111){
3909  if (!(fDoReweightHistoMCPi0 && hReweightMCHistPi0!= 0x0 && PDGCode == 111)){
3910  weight = 1.;
3911  }
3912  }
3913  if (PDGCode == 221){
3914  if (!(fDoReweightHistoMCEta && hReweightMCHistEta!= 0x0 && PDGCode == 221)){
3915  weight = 1.;
3916  }
3917  }
3918  }
3919  if (!isfinite(functionResultData)) weight = 1.;
3920  if (!isfinite(weight)) weight = 1.;
3921  }
3922  } else if (PDGCode == 310 && functionResultMC != 0 && isfinite(functionResultMC)){
3923  weight = functionResultMC;
3924  }
3925  return weight;
3926 }
3927 
3928 
3931 
3932  if( fPeriodEnum == kLHC13bc || // mainly minimum bias
3933  fPeriodEnum == kLHC13de || // mainly triggered
3934  fPeriodEnum == kLHC13b4_fix || fPeriodEnum == kLHC13b4_plus || // MC Pythia 6 (Jet-Jet), anchor LHC13b-e
3935  fPeriodEnum == kLHC13b2_efix || //MC DPMJET, anchr LHC13b+c
3936  fPeriodEnum == kLHC13e7 || //MC HIJING, anchr LHC13b+c
3937  fPeriodEnum == kLHC14b2 //MC HIJING, anchr LHC13b+c
3938  ){
3939  printf(" Gamma Conversion Cuts %s :: pPb Run doing Eta Shift of %f \n\n",(GetCutNumber()).Data(),-0.465);
3940  SetEtaShift(-0.465);
3941  } else if( fPeriodEnum == kLHC13f ) {
3942  printf(" Gamma Conversion Cuts %s :: Pbp Run doing Eta Shift of %f \n\n",(GetCutNumber()).Data(),0.465);
3943  SetEtaShift(+0.465);
3944  }
3945  else printf(" Gamma Conversion Cuts %s :: Automatic Eta Shift requested but Period is not known -> No Shift \n\n",(GetCutNumber()).Data());
3946 }
3947 
3948 //________________________________________________________________________
3949 AliEMCALTriggerPatchInfo* AliConvEventCuts::GetMainTriggerPatch()
3950 {
3951  //get main trigger match; if not known yet, look for it and cache
3952 
3954  return fMainTriggerPatchEMCAL;
3955 
3956  if (!fTriggerPatchInfo) {
3957  AliError(Form("%s: fTriggerPatchInfo not available",GetName()));
3958  return 0;
3959  }
3960 
3961  //number of patches in event
3962  Int_t nPatch = fTriggerPatchInfo->GetEntries();
3963 
3964  //extract main trigger patch
3965  AliEMCALTriggerPatchInfo *patch;
3966  for (Int_t iPatch = 0; iPatch < nPatch; iPatch++) {
3967  patch = (AliEMCALTriggerPatchInfo*)fTriggerPatchInfo->At( iPatch );
3968  if (patch->IsMainTrigger()) {
3969  fMainTriggerPatchEMCAL = patch;
3970  break;
3971  }
3972  }
3973 
3974  return fMainTriggerPatchEMCAL;
3975 }
3976 
3977 
3978 //________________________________________________________________________
3980 {
3981 // cout << "entered EMCAL trigger initialization" << endl;
3982 
3983  // Init the analysis.
3984  if (fCaloTriggersName.IsNull()){
3985  if (event->IsA()==AliESDEvent::Class()){
3986  fCaloTriggersName = "EMCALTrigger";
3987  } else {
3988  fCaloTriggersName = "emcalTrigger";
3989  }
3990  }
3991 
3992  if (!fCaloTriggersName.IsNull() && !fCaloTriggers) {
3993  fCaloTriggers = dynamic_cast<AliVCaloTrigger*>(event->FindListObject(fCaloTriggersName));
3994  if (!fCaloTriggers) {
3995  AliError(Form("%s: Could not retrieve calo triggers %s!", GetName(), fCaloTriggersName.Data()));
3996  return;
3997  }
3998  }
3999 
4000  if (fCaloTriggerPatchInfoName.IsNull()){
4001  if (event->IsA()==AliESDEvent::Class()){
4002  fCaloTriggerPatchInfoName = "EmcalTriggers";
4003  } else {
4004  fCaloTriggerPatchInfoName = "EmcalTriggers";
4005  }
4006  }
4007 
4008  if (!fCaloTriggerPatchInfoName.IsNull() && !fTriggerPatchInfo) {
4009  fTriggerPatchInfo = GetArrayFromEvent(event, fCaloTriggerPatchInfoName.Data(), "AliEMCALTriggerPatchInfo");
4010  if (!fTriggerPatchInfo) {
4011  AliError(Form("%s: Could not retrieve calo trigger patch info %s!", GetName(), fCaloTriggerPatchInfoName.Data()));
4012  return;
4013  }
4014 
4015  }
4016 
4017  fEMCALTrigInitialized = kTRUE;
4018 }
4019 
4020 //________________________________________________________________________
4022  if (!fTriggerPatchInfo)
4023  return 0;
4024  //number of patches in event
4025  Int_t nPatch = fTriggerPatchInfo->GetEntries();
4026 
4027  //loop over patches to define trigger type of event
4028  Int_t nG1 = 0;
4029  Int_t nG2 = 0;
4030  Int_t nJ1 = 0;
4031  Int_t nJ2 = 0;
4032  Int_t nL0 = 0;
4033  AliEMCALTriggerPatchInfo *patch;
4034 // if (nPatch> 0) {cout << "NEW Triggers in this event*********************************" << endl;}
4035  for (Int_t iPatch = 0; iPatch < nPatch; iPatch++) {
4036  patch = (AliEMCALTriggerPatchInfo*)fTriggerPatchInfo->At( iPatch );
4037 // cout << "Patch energy: "<<patch->GetPatchE() << "\t ADC counts: " << patch->GetADCAmp() << endl;
4038 // cout << "Phi: " << patch->GetPhiMin() << " - " << patch->GetPhiMax() << " delta phi: " <<TMath::Abs(patch->GetPhiMin()-patch->GetPhiMax())<< endl;
4039 // cout << "Eta: " << patch->GetEtaMin() << " - " << patch->GetEtaMax() << " delta eta: " <<TMath::Abs(patch->GetEtaMin()-patch->GetEtaMax())<< endl;
4040  if (patch->IsGammaHigh()){
4041 // cout << "fired L1GA high" << endl;
4042  nG1++;
4043  }
4044  if (patch->IsGammaLow()){
4045 // cout << "fired L1GA low" << endl;
4046  nG2++;
4047  }
4048  if (patch->IsJetHigh()){
4049 // cout << "fired L1JE high" << endl;
4050  nJ1++;
4051  }
4052  if (patch->IsJetLow()){
4053 // cout << "fired L1JE low" << endl;
4054  nJ2++;
4055  }
4056  if (patch->IsLevel0()){
4057 // cout << "fired L0" << endl;
4058  nL0++;
4059  }
4060 // cout << patch->GetPatchE() << "\t" << patch->GetADCAmp() << "\t" << patch->IsGammaHigh() << "\t" << patch->IsGammaLow()
4061 // << "\t" << patch->IsJetHigh() << "\t" << patch->IsJetLow() << "\t" << patch->IsLevel0()
4062 // << "\t" << patch->GetPhiMin() << "\t" << patch->GetPhiMax() << "\t" << TMath::Abs(patch->GetPhiMin()-patch->GetPhiMax())
4063 // << "\t" << patch->GetEtaMin() << "\t" << patch->GetEtaMax() << "\t" << TMath::Abs(patch->GetEtaMin()-patch->GetEtaMax()) << endl;
4064  }
4065 
4066  if (nPatch > 0){
4067  AliDebug(2, "Patch summary: ");
4068  AliDebug(2, Form("Number of patches: %d", nPatch));
4069  AliDebug(2, Form("Level0: [%d]" ,nL0));
4070  AliDebug(2, Form("Jet: low[%d], high[%d]" ,nJ2, nJ1));
4071  AliDebug(2, Form("Gamma: low[%d], high[%d]" ,nG2, nG1));
4072  }
4073 
4074 // if (nPatch > 0){
4075 // cout << Form("Number of patches: %d", nPatch) << endl;
4076 // cout << Form("Level0: [%d]" ,nL0) << endl;
4077 // cout << Form("Jet: low[%d], high[%d]" ,nJ2, nJ1) << endl;
4078 // cout << Form("Gamma: low[%d], high[%d]" ,nG2, nG1) << endl;
4079 // }
4080 
4081  ULong_t triggers(0);
4082  if (nG1>0)
4083  SETBIT(triggers, kG1);
4084  if (nG2>0)
4085  SETBIT(triggers, kG2);
4086  if (nJ1>0)
4087  SETBIT(triggers, kJ1);
4088  if (nJ2>0)
4089  SETBIT(triggers, kJ2);
4090  if (nL0>0)
4091  SETBIT(triggers, kL0);
4092  return triggers;
4093 }
4094 
4095 //________________________________________________________________________
4097  // Check if event has a given trigger type
4098  if(t == kND){
4099  return fTriggersEMCAL == 0;
4100  }
4101  return TESTBIT(fTriggersEMCAL, int(t));
4102 }
4103 
4104 
4105 //________________________________________________________________________
4106 TClonesArray *AliConvEventCuts::GetArrayFromEvent(AliVEvent* event, const char *name, const char *clname)
4107 {
4108  // Get array from event.
4109 
4110  TClonesArray *arr = 0;
4111  TString sname(name);
4112  if (!sname.IsNull()) {
4113  arr = dynamic_cast<TClonesArray*>(event->FindListObject(sname));
4114  if (!arr) {
4115  AliWarning(Form("%s: Could not retrieve array with name %s!", GetName(), name));
4116  return 0;
4117  }
4118  } else {
4119  return 0;
4120  }
4121 
4122  if (!clname)
4123  return arr;
4124 
4125  TString objname(arr->GetClass()->GetName());
4126  TClass cls(objname);
4127  if (!cls.InheritsFrom(clname)) {
4128  AliWarning(Form("%s: Objects of type %s in %s are not inherited from %s!",
4129  GetName(), cls.GetName(), name, clname));
4130  return 0;
4131  }
4132  return arr;
4133 }
4134 
4135 //_________________________________________________________________________
4136 Bool_t AliConvEventCuts::IsConversionPrimaryESD( AliMCEvent *mcEvent, Long_t eventpos, Double_t prodVtxX, Double_t prodVtxY, Double_t prodVtxZ){
4137 
4138  if (eventpos < 0) return kFALSE;
4139  TParticle* particle = (TParticle *)mcEvent->Particle(eventpos);
4140  if (!particle) return kFALSE;
4141  if (TMath::Abs(particle->GetPdgCode()) == 11 ){
4142  if (particle->GetMother(0) != -1){
4143  TParticle* particleMother = (TParticle *)mcEvent->Particle(particle->GetMother(0));
4144  if (particleMother){
4145  if (TMath::Abs(particleMother->GetPdgCode()) == 22)
4146  particle = particleMother;
4147  }
4148  }
4149  }
4150  if (particle->GetMother(0) != -1){
4151  Double_t deltaX = particle->Vx() - prodVtxX;
4152  Double_t deltaY = particle->Vy() - prodVtxY;
4153  Double_t deltaZ = particle->Vz() - prodVtxZ;
4154 
4155  //Double_t realRadius2D = TMath::Sqrt(deltaX*deltaX+deltaY*deltaY);
4156  Double_t realRadius3D = TMath::Sqrt(deltaX*deltaX+deltaY*deltaY+deltaZ*deltaZ);
4157 
4158 
4159  Bool_t dalitzCand = kFALSE;
4160 
4161  TParticle* firstmother = (TParticle *)mcEvent->Particle(particle->GetMother(0));
4162  if (!firstmother) return kFALSE;
4163  Int_t pdgCodeFirstMother = firstmother->GetPdgCode();
4164  Bool_t intDecay = kFALSE;
4165  if ( pdgCodeFirstMother == 111 || pdgCodeFirstMother == 221 ) intDecay = kTRUE;
4166  if ( intDecay && TMath::Abs(particle->GetPdgCode()) == 11 ){
4167  dalitzCand = kTRUE;
4168 // cout << "dalitz candidate found" << endl;
4169  }
4170 
4171  Long_t source = particle->GetMother(0);
4172  Bool_t foundExcludedPart = kFALSE;
4173  Bool_t foundShower = kFALSE;
4174  Int_t pdgCodeMotherPrev = 0;
4175  Int_t pdgCodeMotherPPrevMother = 0;
4176  Int_t depth = 0;
4177  if (dalitzCand || realRadius3D < fSecProdBoundary ){
4178 // if (particle->GetPdgCode() == 22){
4179 // cout << endl << endl << "new particle: " << eventpos <<endl;
4180 // cout << particle->GetPdgCode() << "\t" << particle->R() << "\t" << realRadius2D << "\t" << realRadius3D << endl;
4181 // }
4182  while (depth < 20){
4183  TParticle* mother = (TParticle *)mcEvent->Particle(source);
4184  source = mother->GetMother(0);
4185 // if (particle->GetPdgCode() == 22)cout << "eventposition: "<< source << endl;
4186  Int_t pdgCodeMother = mother->GetPdgCode();
4187 // if (particle->GetPdgCode() == 22)cout << "Previous mothers: " << pdgCodeMother << "\t"<< pdgCodeMotherPrev<< "\t" << pdgCodeMotherPPrevMother << endl;
4188  if (pdgCodeMother == pdgCodeMotherPrev && pdgCodeMother == pdgCodeMotherPPrevMother) depth = 20;
4189  if (TMath::Abs(pdgCodeMother) == 11 && TMath::Abs(pdgCodeMotherPrev) == 22 && TMath::Abs(pdgCodeMotherPPrevMother) == 11 ){
4190  foundShower = kTRUE;
4191  depth =20;
4192  }
4193  if (TMath::Abs(pdgCodeMother) == 22 && TMath::Abs(pdgCodeMotherPrev) == 11 && TMath::Abs(pdgCodeMotherPPrevMother) == 22 ){
4194  foundShower = kTRUE;
4195  depth =20;
4196  }
4197 
4198  // particles to be excluded:
4199  // K0s - 310
4200  // K0l - 130
4201  // K+/- - 321
4202  // Lambda - 3122
4203  // Sigma0 - 3212
4204  // Sigma+/- - 3222, 3112
4205  // Cascades - 3322, 3312
4206  if (TMath::Abs(pdgCodeMother) == 310 || TMath::Abs(pdgCodeMother) == 130 || TMath::Abs(pdgCodeMother) == 321 ||
4207  TMath::Abs(pdgCodeMother) == 3122 || TMath::Abs(pdgCodeMother) == 3212 || TMath::Abs(pdgCodeMother) == 3222 ||
4208  TMath::Abs(pdgCodeMother) == 3112 || TMath::Abs(pdgCodeMother) == 3322 || TMath::Abs(pdgCodeMother) == 3312
4209  ) {
4210  foundExcludedPart = kTRUE;
4211  }
4212 // if (particle->GetPdgCode() == 22)cout << mother->GetPdgCode() << "\t" << source << "\t" << foundExcludedPart<< endl;
4213  pdgCodeMotherPPrevMother = pdgCodeMotherPrev;
4214  pdgCodeMotherPrev = pdgCodeMother;
4215  if (source == -1) depth = 20;
4216 
4217 // if (particle->GetPdgCode() == 22)cout << depth << endl;
4218  depth++;
4219  }
4220  }
4221  if (foundExcludedPart){
4222 // if (particle->GetPdgCode() == 22)cout << "This is definitely a secondary, manually excluded" << endl;
4223  return kFALSE;
4224  } else if (dalitzCand && realRadius3D < fSecProdBoundary ){
4225 // if (particle->GetPdgCode() == 22)cout << "This was a decay via a virtual photon" << endl;
4226  return kTRUE;
4227  } else if (foundShower){
4228 // if (particle->GetPdgCode() == 22)cout << "This is a shower" << endl;
4229  return kFALSE;
4230  } else if (realRadius3D >= fSecProdBoundary){
4231 // cout << "This is a secondary, to large production radius" << endl;
4232  return kFALSE;
4233  }
4234  }
4235 
4236  return kTRUE;
4237 }
4238 
4239 //_________________________________________________________________________
4240 Bool_t AliConvEventCuts::IsConversionPrimaryAOD(AliVEvent *event, AliAODMCParticle* AODMCParticle, Double_t prodVtxX, Double_t prodVtxY, Double_t prodVtxZ){
4241 
4242  TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(event->FindListObject(AliAODMCParticle::StdBranchName()));
4243  if (AODMCTrackArray == NULL) return kFALSE;
4244  AliAODMCParticle* currentParticle = AODMCParticle;
4245  if (TMath::Abs(currentParticle->GetPdgCode()) == 11 ){
4246  if (currentParticle->GetMother() != -1){
4247  AliAODMCParticle* particleMother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(currentParticle->GetMother()));
4248  if (particleMother){
4249  if (TMath::Abs(particleMother->GetPdgCode()) == 22)
4250  currentParticle = particleMother;
4251  }
4252  }
4253  }
4254  if (currentParticle->GetMother() > -1){
4255  Double_t deltaX = currentParticle->Xv() - prodVtxX;
4256  Double_t deltaY = currentParticle->Yv() - prodVtxY;
4257  Double_t deltaZ = currentParticle->Zv() - prodVtxZ;
4258 
4259  //Double_t realRadius2D = TMath::Sqrt(deltaX*deltaX+deltaY*deltaY);
4260  Double_t realRadius3D = TMath::Sqrt(deltaX*deltaX+deltaY*deltaY+deltaZ*deltaZ);
4261 
4262  Bool_t dalitzCand = kFALSE;
4263 
4264  AliAODMCParticle* firstmother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(currentParticle->GetMother()));
4265  if (!firstmother) return kFALSE;
4266  Int_t pdgCodeFirstMother = firstmother->GetPdgCode();
4267  Bool_t intDecay = kFALSE;
4268  if ( pdgCodeFirstMother == 111 || pdgCodeFirstMother == 221 ) intDecay = kTRUE;
4269  if ( intDecay && TMath::Abs(currentParticle->GetPdgCode()) == 11 ){
4270  dalitzCand = kTRUE;
4271 // cout << "dalitz candidate found" << endl;
4272  }
4273 
4274  Long_t source = currentParticle->GetMother();
4275  Bool_t foundExcludedPart = kFALSE;
4276  Bool_t foundShower = kFALSE;
4277  Int_t pdgCodeMotherPrev = 0;
4278  Int_t pdgCodeMotherPPrevMother = 0;
4279  Int_t depth = 0;
4280  if (dalitzCand || realRadius3D < fSecProdBoundary ){
4281 // if (currentParticle->GetPdgCode() == 22){
4282 // cout << endl << endl << "new particle: " << eventpos <<endl;
4283 // cout << currentParticle->GetPdgCode() << "\t" << currentParticle->R() << "\t" << realRadius2D << "\t" << realRadius3D << endl;
4284 // }
4285  while (depth < 20){
4286  AliAODMCParticle* mother = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(source));
4287  source = mother->GetMother();
4288 // if (currentParticle->GetPdgCode() == 22)cout << "eventposition: "<< source << endl;
4289  Int_t pdgCodeMother = mother->GetPdgCode();
4290 // if (currentParticle->GetPdgCode() == 22)cout << "Previous mothers: " << pdgCodeMother << "\t"<< pdgCodeMotherPrev<< "\t" << pdgCodeMotherPPrevMother << endl;
4291  if (pdgCodeMother == pdgCodeMotherPrev && pdgCodeMother == pdgCodeMotherPPrevMother) depth = 20;
4292  if (TMath::Abs(pdgCodeMother) == 11 && TMath::Abs(pdgCodeMotherPrev) == 22 && TMath::Abs(pdgCodeMotherPPrevMother) == 11 ){
4293  foundShower = kTRUE;
4294  depth =20;
4295  }
4296  if (TMath::Abs(pdgCodeMother) == 22 && TMath::Abs(pdgCodeMotherPrev) == 11 && TMath::Abs(pdgCodeMotherPPrevMother) == 22 ){
4297  foundShower = kTRUE;
4298  depth =20;
4299  }
4300 
4301  // particles to be excluded:
4302  // K0s - 310
4303  // K0l - 130
4304  // K+/- - 321
4305  // Lambda - 3122
4306  // Sigma0 - 3212
4307  // Sigma+/- - 3222, 3112
4308  // Cascades - 3322, 3312
4309  if (TMath::Abs(pdgCodeMother) == 310 || TMath::Abs(pdgCodeMother) == 130 || TMath::Abs(pdgCodeMother) == 321 ||
4310  TMath::Abs(pdgCodeMother) == 3122 || TMath::Abs(pdgCodeMother) == 3212 || TMath::Abs(pdgCodeMother) == 3222 ||
4311  TMath::Abs(pdgCodeMother) == 3112 || TMath::Abs(pdgCodeMother) == 3322 || TMath::Abs(pdgCodeMother) == 3312)
4312  {
4313  foundExcludedPart = kTRUE;
4314  }
4315 // if (currentParticle->GetPdgCode() == 22)cout << mother->GetPdgCode() << "\t" << source << "\t" << foundExcludedPart<< endl;
4316  pdgCodeMotherPPrevMother = pdgCodeMotherPrev;
4317  pdgCodeMotherPrev = pdgCodeMother;
4318  if (source == -1) depth = 20;
4319 
4320 // if (currentParticle->GetPdgCode() == 22)cout << depth << endl;
4321  depth++;
4322  }
4323  }
4324  if (foundExcludedPart){
4325 // if (currentParticle->GetPdgCode() == 22)cout << "This is definitely a secondary, manually excluded" << endl;
4326  return kFALSE;
4327  } else if (dalitzCand && realRadius3D < fSecProdBoundary ){
4328 // if (currentParticle->GetPdgCode() == 22)cout << "This was a decay via a virtual photon" << endl;
4329  return kTRUE;
4330  } else if (foundShower){
4331 // if (currentParticle->GetPdgCode() == 22)cout << "This is a shower" << endl;
4332  return kFALSE;
4333  } else if (realRadius3D >= fSecProdBoundary){
4334 // cout << "This is a secondary, too large production radius" << endl;
4335  return kFALSE;
4336  }
4337  }
4338 
4339  return kTRUE;
4340 }
4341 
4342 
4343 //________________________________________________________________________
4344 Int_t AliConvEventCuts::SecondaryClassificationPhoton( TParticle *particle, AliMCEvent* mcEvent, Bool_t isConversion ){
4345  if (particle != NULL && mcEvent != NULL){
4346  Int_t pdgSecondary = 0;
4347  if (!isConversion){
4348  //Bool_t hasMother = kFALSE;
4349  //Bool_t hasGrandMother = kFALSE;
4350  Long_t motherID = particle->GetMother(0);
4351  Long_t grandMotherID = -1;
4352  // is the photon a direct photons, without a mother?
4353  if (motherID > -1){
4354  //hasMother = kTRUE;
4355  grandMotherID = mcEvent->Particle(motherID)->GetMother(0);
4356  // is the meson a primary?
4357  if (grandMotherID > -1){
4358  //hasGrandMother = kTRUE;
4359  pdgSecondary = mcEvent->Particle(grandMotherID)->GetPdgCode();
4360  }
4361  }
4362  } else {
4363  //Bool_t hasMother = kFALSE;
4364  //Bool_t hasGrandMother = kFALSE;
4365  //Bool_t hasGreatGrandMother = kFALSE;
4366  Long_t motherID = particle->GetMother(0);
4367  Long_t grandMotherID = -1;
4368  Long_t greatGrandMotherID = -1;
4369  // is the electron a direct electron, without a mother?
4370  if (motherID > -1){
4371  //hasMother = kTRUE;
4372  grandMotherID = mcEvent->Particle(motherID)->GetMother(0);
4373  // is the photon a direct photons, without a mother?
4374  if (grandMotherID > -1){
4375  //hasGrandMother = kTRUE;
4376  greatGrandMotherID = mcEvent->Particle(grandMotherID)->GetMother(0);
4377  // is the meson a primary?
4378  if (greatGrandMotherID > -1){
4379  //hasGreatGrandMother = kTRUE;
4380  pdgSecondary = mcEvent->Particle(greatGrandMotherID)->GetPdgCode();
4381  }
4382  }
4383  }
4384  }
4385  // is the secondary photon from a lambda
4386  if (TMath::Abs(pdgSecondary) == 3122 )
4387  return 3;
4388  // is the secondary photon from a K0s
4389  else if ( TMath::Abs(pdgSecondary) == 310 )
4390  return 2;
4391  // is the secondary photon from a K0l
4392  else if ( TMath::Abs(pdgSecondary) == 130 )
4393  return 5;
4394  // is the secondary photon from a eta
4395  else if ( TMath::Abs(pdgSecondary) == 221 )
4396  return 4;
4397  // is the secondary photon from something else
4398  else if ( TMath::Abs(pdgSecondary) != 0 )
4399  return 1;
4400 
4401  }
4402 
4403  return 0;
4404 }
4405 
4406 //________________________________________________________________________
4407 Int_t AliConvEventCuts::SecondaryClassificationPhotonAOD( AliAODMCParticle *particle, TClonesArray *aodmcArray, Bool_t isConversion ){
4408  if (particle != NULL && aodmcArray != NULL){
4409  Int_t pdgSecondary = 0;
4410  if (!isConversion){
4411  //Bool_t hasMother = kFALSE;
4412  //Bool_t hasGrandMother = kFALSE;
4413  Long_t motherID = particle->GetMother();
4414  Long_t grandMotherID = -1;
4415  // is the photon a direct photons, without a mother?
4416  if (motherID > -1){
4417  //hasMother = kTRUE;
4418  grandMotherID = ((AliAODMCParticle*)aodmcArray->At(motherID))->GetMother();
4419  // is the meson a primary?
4420  if (grandMotherID > -1){
4421  //hasGrandMother = kTRUE;
4422  pdgSecondary = ((AliAODMCParticle*)aodmcArray->At(grandMotherID))->GetPdgCode();
4423  }
4424  }
4425  } else {
4426  //Bool_t hasMother = kFALSE;
4427  //Bool_t hasGrandMother = kFALSE;
4428  //Bool_t hasGreatGrandMother = kFALSE;
4429  Long_t motherID = particle->GetMother();
4430  Long_t grandMotherID = -1;
4431  Long_t greatGrandMotherID = -1;
4432  // is the electron a direct electron, without a mother?
4433  if (motherID > -1){
4434  //hasMother = kTRUE;
4435  grandMotherID = ((AliAODMCParticle*)aodmcArray->At(motherID))->GetMother();
4436  // is the photon a direct photons, without a mother?
4437  if (grandMotherID > -1){
4438  //hasGrandMother = kTRUE;
4439  greatGrandMotherID = ((AliAODMCParticle*)aodmcArray->At(grandMotherID))->GetMother();
4440  // is the meson a primary?
4441  if (greatGrandMotherID > -1){
4442  //hasGreatGrandMother = kTRUE;
4443  pdgSecondary = ((AliAODMCParticle*)aodmcArray->At(greatGrandMotherID))->GetPdgCode();
4444  }
4445  }
4446  }
4447  }
4448  // is the secondary photon from a lambda
4449  if (TMath::Abs(pdgSecondary) == 3122 )
4450  return 3;
4451  // is the secondary photon from a K0s
4452  else if ( TMath::Abs(pdgSecondary) == 310 )
4453  return 2;
4454  // is the secondary photon from a K0l
4455  else if ( TMath::Abs(pdgSecondary) == 130 )
4456  return 5;
4457  // is the secondary photon from a eta
4458  else if ( TMath::Abs(pdgSecondary) == 221 )
4459  return 4;
4460  // is the secondary photon from something else
4461  else if ( TMath::Abs(pdgSecondary) != 0 )
4462  return 1;
4463 
4464  }
4465 
4466  return 0;
4467 }
4468 
4470 
4471  if (periodName.CompareTo("") == 0){
4472  periodName = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask(fV0ReaderName.Data()))->GetPeriodName();
4473  }
4474 
4475  if (periodName.CompareTo("") == 0) {
4477  fEnergyEnum = kUnset;
4478  AliError("No correct period could be set, periodName string empty");
4479  return;
4480  }
4481 
4482  // Data
4483  if (periodName.CompareTo("LHC10b") == 0 || periodName.CompareTo("LHC10c") == 0 || periodName.CompareTo("LHC10d") == 0 || periodName.CompareTo("LHC10e") == 0 ||
4484  periodName.CompareTo("LHC10f") == 0 || periodName.CompareTo("LHC10g") == 0 || periodName.CompareTo("LHC10bg") == 0
4485  ){
4487  fEnergyEnum = k7TeV;
4488  } else if (periodName.CompareTo("LHC10h") == 0) {
4489  fPeriodEnum = kLHC10h;
4491  } else if (periodName.CompareTo("LHC11a") == 0) {
4492  fPeriodEnum = kLHC11a;
4494  } else if (periodName.CompareTo("LHC11b") == 0) {
4495  fPeriodEnum = kLHC11b;
4496  fEnergyEnum = k7TeV;
4497  } else if (periodName.CompareTo("LHC11c") == 0 || periodName.CompareTo("LHC11d") == 0 || periodName.CompareTo("LHC11e") == 0 || periodName.CompareTo("LHC11f") == 0 ||
4498  periodName.CompareTo("LHC11g") == 0
4499  ) {
4501  fEnergyEnum = k7TeV;
4502  } else if (periodName.CompareTo("LHC11h") == 0) {
4503  fPeriodEnum = kLHC11h;
4505  } else if (periodName.CompareTo("LHC12a") == 0 || periodName.CompareTo("LHC12b") == 0 || periodName.CompareTo("LHC12c") == 0 || periodName.CompareTo("LHC12d") == 0 ||
4506  periodName.CompareTo("LHC12e") == 0 || periodName.CompareTo("LHC12f") == 0 || periodName.CompareTo("LHC12g") == 0 || periodName.CompareTo("LHC12h") == 0 ||
4507  periodName.CompareTo("LHC12i") == 0 || periodName.CompareTo("LHC12ai") == 0
4508  ) {
4509  fPeriodEnum = kLHC12;
4510  fEnergyEnum = k8TeV;
4511  } else if (periodName.CompareTo("LHC13b") == 0 || periodName.CompareTo("LHC13c") == 0 || periodName.CompareTo("LHC13bc") == 0){
4514  } else if (periodName.CompareTo("LHC13d") == 0 || periodName.CompareTo("LHC13e") == 0 || periodName.CompareTo("LHC13de") == 0){
4517  } else if (periodName.CompareTo("LHC13f") == 0 ){
4518  fPeriodEnum = kLHC13f;
4520  } else if (periodName.CompareTo("LHC13g") == 0 ){
4521  fPeriodEnum = kLHC13g;
4523  } else if ( periodName.CompareTo("LHC15f") == 0 || periodName.CompareTo("LHC15g") == 0 || periodName.CompareTo("LHC15h") == 0 || periodName.CompareTo("LHC15i") == 0 ||
4524  periodName.CompareTo("LHC15j") == 0 || periodName.CompareTo("LHC15k") == 0 || periodName.CompareTo("LHC15l") == 0 || periodName.CompareTo("LHC15m") == 0 ||
4525  periodName.CompareTo("LHC15fm") == 0
4526  ) {
4528  fEnergyEnum = k13TeV;
4529  } else if (periodName.CompareTo("LHC15n") == 0 ){
4530  fPeriodEnum = kLHC15n;
4531  fEnergyEnum = k5TeV;
4532  } else if (periodName.CompareTo("LHC15o") == 0 ){
4533  fPeriodEnum = kLHC15o;
4535  } else if ( periodName.CompareTo("LHC16d") == 0 || periodName.CompareTo("LHC16e") == 0 || periodName.CompareTo("LHC16g") == 0 || periodName.CompareTo("LHC16h") == 0 ||
4536  periodName.CompareTo("LHC16i") == 0 || periodName.CompareTo("LHC16j") == 0 || periodName.CompareTo("LHC16k") == 0 || periodName.CompareTo("LHC16l") == 0 ||
4537  periodName.CompareTo("LHC16m") == 0 || periodName.CompareTo("LHC16n") == 0 || periodName.CompareTo("LHC16o") == 0 || periodName.CompareTo("LHC16p") == 0){
4539  fEnergyEnum = k13TeV;
4540  } else if (periodName.CompareTo("LHC16f") == 0 ){
4543  } else if (periodName.CompareTo("LHC16q") == 0 ){
4544  fPeriodEnum = kLHC16q;
4546  } else if (periodName.CompareTo("LHC16r") == 0 ){
4547  fPeriodEnum = kLHC16r;
4549  } else if (periodName.CompareTo("LHC16s") == 0 ){
4550  fPeriodEnum = kLHC16s;
4552  } else if (periodName.CompareTo("LHC16t") == 0 ){
4553  fPeriodEnum = kLHC16t;
4555  } else if ( periodName.CompareTo("LHC17c") == 0 || periodName.CompareTo("LHC17d") == 0 || periodName.CompareTo("LHC17e") == 0 || periodName.CompareTo("LHC17f") == 0 ||
4556  periodName.CompareTo("LHC17h") == 0 || periodName.CompareTo("LHC17i") == 0 || periodName.CompareTo("LHC17j") == 0 || periodName.CompareTo("LHC17k") == 0 ||
4557  periodName.CompareTo("LHC17l") == 0 || periodName.CompareTo("LHC17m") == 0 || periodName.CompareTo("LHC17o") == 0 ){
4559  fEnergyEnum = k13TeV;
4560  } else if ( periodName.CompareTo("LHC17g") == 0 ){
4563  } else if ( periodName.CompareTo("LHC17n") == 0 ){
4564  fPeriodEnum = kLHC17n;
4566 
4567  // LHC10x anchored MCs
4568  } else if (periodName.CompareTo("LHC10d1") == 0){
4570  fEnergyEnum = k7TeV;
4571  } else if (periodName.CompareTo("LHC10d2") == 0){
4573  fEnergyEnum = k7TeV;
4574  } else if (periodName.CompareTo("LHC10d4a") == 0){
4576  fEnergyEnum = k7TeV;
4577  } else if (periodName.CompareTo("LHC10d4") == 0){
4579  fEnergyEnum = k7TeV;
4580  } else if (periodName.CompareTo("LHC10e12") == 0){
4582  fEnergyEnum = k900GeV;
4583  } else if (periodName.CompareTo("LHC10e13") == 0){
4585  fEnergyEnum = k900GeV;
4586  } else if (periodName.CompareTo("LHC10e20") == 0){
4588  fEnergyEnum = k7TeV;
4589  } else if (periodName.CompareTo("LHC10e21") == 0){
4591  fEnergyEnum = k7TeV;
4592  } else if (periodName.CompareTo("LHC10f6a") == 0){
4594  fEnergyEnum = k7TeV;
4595  } else if (periodName.CompareTo("LHC10f6") == 0){
4597  fEnergyEnum = k7TeV;
4598  } else if (periodName.Contains("LHC14j4")){
4600  fEnergyEnum = k7TeV;
4601  } else if (periodName.CompareTo("LHC13d2") == 0){
4604  } else if (periodName.CompareTo("LHC13d2b") == 0){
4607  } else if (periodName.CompareTo("LHC12a11a") == 0){
4610  } else if (periodName.CompareTo("LHC12a11b") == 0){
4613  } else if (periodName.CompareTo("LHC12a11c") == 0){
4616  } else if (periodName.CompareTo("LHC12a11d") == 0){
4619  } else if (periodName.CompareTo("LHC12a11e") == 0){
4622  } else if (periodName.CompareTo("LHC12a11f") == 0){
4625  // LHC11x anchored MCs
4626  } else if (periodName.CompareTo("LHC12a15c") == 0){
4629  } else if (periodName.Contains("LHC12f1a") ){
4632  } else if (periodName.Contains("LHC12f1b") ){
4635  } else if (periodName.Contains("LHC12i3") ){
4638  } else if (periodName.CompareTo("LHC15g1a") == 0){
4641  } else if (periodName.CompareTo("LHC15g1b") == 0){
4644  } else if (periodName.CompareTo("LHC13e4") == 0){
4646  fEnergyEnum = k7TeV;
4647  } else if (periodName.CompareTo("LHC13e5") == 0){
4649  fEnergyEnum = k7TeV;
4650  } else if (periodName.CompareTo("LHC14k1a") == 0){
4652  fEnergyEnum = k7TeV;
4653  } else if (periodName.CompareTo("LHC14k1b") == 0){
4655  fEnergyEnum = k7TeV;
4656  } else if (periodName.CompareTo("LHC12a15f") == 0){
4658  fEnergyEnum = k7TeV;
4659  } else if (periodName.CompareTo("LHC12a15g") == 0){
4661  fEnergyEnum = k7TeV;
4662  } else if (periodName.CompareTo("LHC12f2a") == 0){
4664  fEnergyEnum = k7TeV;
4665  } else if (periodName.CompareTo("LHC14a1a") == 0){
4668  } else if (periodName.CompareTo("LHC14a1b") == 0){
4671  } else if (periodName.CompareTo("LHC14a1c") == 0){
4674  // LHC12x anchored MCs
4675  } else if (periodName.CompareTo("LHC14e2a") == 0){
4677  fEnergyEnum = k8TeV;
4678  } else if (periodName.CompareTo("LHC14e2b") == 0){
4680  fEnergyEnum = k8TeV;
4681  } else if (periodName.CompareTo("LHC14e2c") == 0){
4683  fEnergyEnum = k8TeV;
4684  } else if (periodName.Contains("LHC15h1")){
4686  fEnergyEnum = k8TeV;
4687  } else if (periodName.Contains("LHC15h2")){
4689  fEnergyEnum = k8TeV;
4690  } else if (periodName.CompareTo("LHC16c2") == 0){
4692  fEnergyEnum = k8TeV;
4693  } else if (periodName.CompareTo("LHC16c2_plus") == 0){
4695  fEnergyEnum = k8TeV;
4696  // LHC13x anchored MCs
4697  } else if (periodName.Contains("LHC13b2_efix")){
4700  } else if (periodName.CompareTo("LHC13e7") == 0){
4703  } else if (periodName.CompareTo("LHC14b2") == 0){
4706  } else if (periodName.CompareTo("LHC13b4_fix") == 0){
4709  } else if (periodName.CompareTo("LHC13b4_plus") == 0){
4712  } else if (periodName.CompareTo("LHC16c3a") == 0){
4715  } else if (periodName.CompareTo("LHC16c3b") == 0){
4718  } else if (periodName.CompareTo("LHC16c3c") == 0){
4721  } else if (periodName.CompareTo("LHC15g2") == 0){
4724  } else if (periodName.CompareTo("LHC15a3a") == 0){
4727  } else if (periodName.CompareTo("LHC15a3a_plus") == 0){
4730  } else if (periodName.CompareTo("LHC15a3b") == 0){
4733  } else if (periodName.CompareTo("LHC15d3a") == 0){
4736  } else if (periodName.CompareTo("LHC15d3b") == 0){
4739  // LHC15x anchored MCs
4740  } else if (periodName.CompareTo("LHC15g3a3") == 0){
4742  fEnergyEnum = k13TeV;
4743  } else if (periodName.CompareTo("LHC15g3a") == 0){
4745  fEnergyEnum = k13TeV;
4746  } else if (periodName.CompareTo("LHC15g3c2") == 0){
4748  fEnergyEnum = k13TeV;
4749  } else if (periodName.CompareTo("LHC15g3c3") == 0){
4751  fEnergyEnum = k13TeV;
4752  } else if (periodName.CompareTo("LHC15g3") == 0){
4754  fEnergyEnum = k13TeV;
4755  } else if (periodName.CompareTo("LHC16a2a") == 0){
4757  fEnergyEnum = k13TeV;
4758  } else if (periodName.CompareTo("LHC16a2b") == 0){
4760  fEnergyEnum = k13TeV;
4761  } else if (periodName.CompareTo("LHC16a2c") == 0){
4763  fEnergyEnum = k13TeV;
4764  } else if (periodName.CompareTo("LHC15l1a2") == 0){
4766  fEnergyEnum = k5TeV;
4767  } else if (periodName.CompareTo("LHC15l1b2") == 0){
4769  fEnergyEnum = k5TeV;
4770  } else if (periodName.CompareTo("LHC16h3") == 0){
4772  fEnergyEnum = k5TeV;
4773  } else if (periodName.CompareTo("LHC16h8a") == 0){
4775  fEnergyEnum = k5TeV;
4776  } else if (periodName.CompareTo("LHC16h8b") == 0){
4778  fEnergyEnum = k5TeV;
4779  } else if (periodName.CompareTo("LHC16k5a") == 0){
4781  fEnergyEnum = k5TeV;
4782  } else if (periodName.CompareTo("LHC16k5b") == 0){
4784  fEnergyEnum = k5TeV;
4785  } else if (periodName.CompareTo("LHC16k3a") == 0){
4787  fEnergyEnum = k5TeV;
4788  } else if (periodName.CompareTo("LHC16k3a2") == 0){
4790  fEnergyEnum = k5TeV;
4791  } else if (periodName.CompareTo("LHC17e2") == 0){
4793  fEnergyEnum = k5TeV;
4794  } else if (periodName.Contains("LHC15k1a1")){
4797  } else if (periodName.Contains("LHC15k1a2")){
4800  } else if (periodName.Contains("LHC15k1a3")){
4803  } else if (periodName.Contains("LHC16j7")){
4806  } else if (periodName.Contains("LHC16h4")){
4809  } else if (periodName.Contains("LHC16g1")){
4812  } else if (periodName.Contains("LHC16g1a")){
4815  } else if (periodName.Contains("LHC16g1b")){
4818  } else if (periodName.Contains("LHC16g1c")){
4821  } else if (periodName.Contains("LHC16g2")){
4824  } else if (periodName.Contains("LHC16g3")){
4827  } else if (periodName.Contains("LHC16h2a")){
4830  } else if (periodName.Contains("LHC16h2b")){
4833  } else if (periodName.Contains("LHC16h2c")){
4836  } else if (periodName.CompareTo("LHC16k3b") == 0){
4839  } else if (periodName.CompareTo("LHC16k3b2") == 0){
4842 
4843  // LHC16x anchored MCs
4844  // 13TeV LHC16* anchors full field Pythia 8 MB
4845  } else if ( periodName.CompareTo("LHC16P1Pyt8") == 0 || periodName.CompareTo("LHC17f6") == 0 || periodName.CompareTo("LHC17d17") == 0 || periodName.CompareTo("LHC17f5") == 0 ||
4846  periodName.CompareTo("LHC17d3") == 0 || periodName.CompareTo("LHC17e5") == 0 || periodName.CompareTo("LHC17d20a1") == 0 || periodName.CompareTo("LHC17d20a1_extra") == 0 ||
4847  periodName.CompareTo("LHC17d20a2") == 0 || periodName.CompareTo("LHC17d20a2_extra") == 0 || periodName.CompareTo("LHC17d16") == 0 || periodName.CompareTo("LHC17d18") == 0 ||
4848  periodName.CompareTo("LHC17f9") == 0 || periodName.CompareTo("LHC17f9_test") == 0){
4850  fEnergyEnum = k13TeV;
4851  // 13TeV LHC16* anchors low field Pythia 8 MB
4852  } else if ( periodName.CompareTo("LHC16P1Pyt8LowB") == 0 || periodName.CompareTo("LHC17d1") == 0 ){
4855  // 13TeV LHC16* anchors full field EPOS MB
4856  } else if ( periodName.CompareTo("LHC16P1EPOS") == 0 || periodName.CompareTo("LHC17d20b1") == 0 || periodName.CompareTo("LHC17d20b2") == 0 ){
4858  fEnergyEnum = k13TeV;
4859  // 13TeV LHC16* anchors full field JJ Pythia 8 MB
4860  } else if ( periodName.CompareTo("LHC16P1JJ") == 0 || periodName.CompareTo("LHC17f8a") == 0 || periodName.CompareTo("LHC17f8c") == 0 || periodName.CompareTo("LHC17f8d") == 0 ||
4861  periodName.CompareTo("LHC17f8e") == 0 ){
4863  fEnergyEnum = k13TeV;
4864  // 13TeV LHC16* anchors low field JJ Pythia 8 MB
4865  } else if ( periodName.CompareTo("LHC16P1JJLowB") == 0 || periodName.CompareTo("LHC17f8b") == 0 ){
4868 
4869  // LHC16qt anchored MCs
4870  } else if (periodName.CompareTo("LHC17f2a") == 0){
4873  } else if (periodName.CompareTo("LHC17f2a_fast") == 0){
4876  } else if (periodName.CompareTo("LHC17f2a_cent") == 0){
4879  } else if (periodName.CompareTo("LHC17f2a_cent_woSDD") == 0){
4882  } else if (periodName.CompareTo("LHC17f2a_fast_fix") == 0){
4885  } else if (periodName.CompareTo("LHC17f2a_cent_fix") == 0){
4888  } else if (periodName.CompareTo("LHC17f2a_cent_woSDD_fix") == 0){
4891  } else if (periodName.CompareTo("LHC17f2b") == 0){
4894  } else if (periodName.CompareTo("LHC17f2b_fast") == 0){
4897  } else if (periodName.CompareTo("LHC17f2b_cent") == 0){
4900  } else if (periodName.CompareTo("LHC17f2b_cent_woSDD") == 0){
4903  } else if (periodName.CompareTo("LHC17g8a_fast") == 0){
4906  } else if (periodName.CompareTo("LHC17g8a_cent_woSDD") == 0){
4909  // LHC16r anchored MCs
4910  } else if (periodName.CompareTo("LHC17a3a") == 0){
4912  fEnergyEnum =