]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGGA/GammaConv/AliConversionCuts.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
changed standard cut for pp & PbPb to have the proper length
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / AliConversionCuts.cxx
CommitLineData
3b77b2d1 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
92efd725 3 * *
4 * Authors: Svein Lindal, Daniel Lohner *
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. *
3b77b2d1 14 **************************************************************************/
15
16////////////////////////////////////////////////
e5b6e8a6 17//---------------------------------------------
3b77b2d1 18// Class handling all kinds of selection cuts for
19// Gamma Conversion analysis
20//---------------------------------------------
21////////////////////////////////////////////////
22
23#include "AliConversionCuts.h"
24
25#include "AliKFVertex.h"
26#include "AliAODTrack.h"
27#include "AliESDtrack.h"
28#include "AliAnalysisManager.h"
29#include "AliInputEventHandler.h"
92efd725 30#include "AliMCEventHandler.h"
3b77b2d1 31#include "AliAODHandler.h"
32#include "AliPIDResponse.h"
33#include "TH1.h"
34#include "TH2.h"
35#include "AliStack.h"
36#include "AliAODConversionMother.h"
37#include "TObjString.h"
38#include "AliAODEvent.h"
39#include "AliESDEvent.h"
40#include "AliCentrality.h"
41#include "TList.h"
ae947965 42#include "TFile.h"
92efd725 43#include "AliLog.h"
ca91a3e1 44#include "AliGenCocktailEventHeader.h"
a280ac15 45#include "AliGenDPMjetEventHeader.h"
46#include "AliGenPythiaEventHeader.h"
47#include "AliGenHijingEventHeader.h"
e5b6e8a6 48#include "AliTriggerAnalysis.h"
49#include "AliV0ReaderV1.h"
ae947965 50#include "AliAODMCParticle.h"
51#include "AliAODMCHeader.h"
92efd725 52
3b77b2d1 53class iostream;
54
55using namespace std;
56
57ClassImp(AliConversionCuts)
58
59
60const char* AliConversionCuts::fgkCutNames[AliConversionCuts::kNCuts] = {
4803eb1f 61 "HeavyIon",//0
62 "CentralityMin",//1
63 "CentralityMax",//2
64 "SelectV0AND",//3
65 "MultiplicityMethod",//4
66 "RemovePileUp",//5
67 "RejectExtraSignals",//6
68 "V0FinderType",//7
69 "EtaCut",//8
70 "MinRCut",//9
71 "SinglePtCut",//10
72 "ClsTPCCut", //11
73 "ededxSigmaCut",//12
74 "pidedxSigmaCut",//13
75 "piMomdedxSigmaCut",//14
76 "piMaxMomdedxSigmaCut",//15
77 "LowPRejectionSigmaCut",//16
78 "TOFelectronPID",//17
79 "QtMaxCut",//18
80 "Chi2GammaCut", //19
81 "PsiPair", //20
82 "DoPhotonAsymmetryCut",//21
83 "CosinePointingAngle", //22
84 "SharedElectronCuts", //23
85 "RejectToCloseV0s", //24
86 "DcaRPrimVtx", //25
87 "DcaZPrimVtx" //26
3b77b2d1 88};
89
90
91//________________________________________________________________________
4a0aab28 92AliConversionCuts::AliConversionCuts(const char *name,const char *title) :
93 AliAnalysisCuts(name,title),
94 fHistograms(NULL),
e5b6e8a6 95 fHeaderList(NULL),
4a0aab28 96 fPIDResponse(NULL),
97 fEventQuality(-1),
98 fMaxR(200),
99 fMinR(0),
100 fEtaCut(0.9),
101 fEtaCutMin(-0.1),
102 fPtCut(0.02),
103 fSinglePtCut(0),
104 fMaxZ(1000),
105 fMinClsTPC(0.),
106 fMinClsTPCToF(0.),
107 fLineCutZRSlope(0.),
108 fLineCutZValue(0),
109 fLineCutZRSlopeMin(0.),
110 fLineCutZValueMin(0),
111 fChi2CutConversion(1000),
112 fPIDProbabilityCutNegativeParticle(0),
113 fPIDProbabilityCutPositiveParticle(0),
114 fDodEdxSigmaCut(kTRUE),
e5b6e8a6 115 fDoTOFsigmaCut(kFALSE),
4a0aab28 116 fPIDTRDEfficiency(1),
117 fDoTRDPID(kFALSE),
118 fPIDnSigmaAboveElectronLine(100),
119 fPIDnSigmaBelowElectronLine(-100),
e5b6e8a6 120 fTofPIDnSigmaAboveElectronLine(100),
121 fTofPIDnSigmaBelowElectronLine(-100),
4a0aab28 122 fPIDnSigmaAbovePionLine(0),
123 fPIDnSigmaAbovePionLineHighPt(-100),
124 fPIDMinPnSigmaAbovePionLine(0),
125 fPIDMaxPnSigmaAbovePionLine(0),
126 fDoKaonRejectionLowP(kFALSE),
127 fDoProtonRejectionLowP(kFALSE),
128 fDoPionRejectionLowP(kFALSE),
129 fPIDnSigmaAtLowPAroundKaonLine(0),
130 fPIDnSigmaAtLowPAroundProtonLine(0),
131 fPIDnSigmaAtLowPAroundPionLine(0),
e5b6e8a6 132 fPIDMinPKaonRejectionLowP(1.5),
133 fPIDMinPProtonRejectionLowP(2),
4a0aab28 134 fPIDMinPPionRejectionLowP(0),
135 fDoQtGammaSelection(kTRUE),
e5b6e8a6 136 fDoHighPtQtGammaSelection(kFALSE),
4a0aab28 137 fQtMax(100),
e5b6e8a6 138 fHighPtQtMax(0.),
139 fPtBorderForQt(0),
4a0aab28 140 fXVertexCut(0.),
141 fYVertexCut(0.),
142 fZVertexCut(0.),
143 fNSigmaMass(0.),
144 fUseEtaMinCut(kFALSE),
145 fUseOnFlyV0Finder(kTRUE),
146 fDoPhotonAsymmetryCut(kTRUE),
147 fMinPPhotonAsymmetryCut(100.),
148 fMinPhotonAsymmetry(0.),
0a2b2b4b 149 fIsHeavyIon(0),
4a0aab28 150 fDetectorCentrality(0),
151 fModCentralityClass(0),
152 fMaxVertexZ(10),
153 fCentralityMin(0),
154 fCentralityMax(0),
155 fUseCorrectedTPCClsInfo(kFALSE),
156 fUseTOFpid(kFALSE),
157 fMultiplicityMethod(0),
e5b6e8a6 158 fSpecialTrigger(0),
4a0aab28 159 fRemovePileUp(kFALSE),
160 fOpeningAngle(0.005),
161 fPsiPairCut(10000),
162 fCosPAngleCut(10000),
163 fDoToCloseV0sCut(kFALSE),
164 fRejectExtraSignals(0),
165 fminV0Dist(200.),
166 fDoSharedElecCut(kFALSE),
167 fOfflineTriggerMask(0),
e5b6e8a6 168 fHasV0AND(kTRUE),
169 fIsSDDFired(kTRUE),
4a0aab28 170 fRandom(0),
e5b6e8a6 171 fElectronArraySize(500),
4a0aab28 172 fElectronLabelArray(NULL),
4803eb1f 173 fDCAZPrimVtxCut(1000),
174 fDCARPrimVtxCut(1000),
e5b6e8a6 175 fConversionPointXArray(0.0),
176 fConversionPointYArray(0.0),
177 fConversionPointZArray(0.0),
178 fnHeaders(0),
4a0aab28 179 fNotRejectedStart(NULL),
180 fNotRejectedEnd(NULL),
181 fGeneratorNames(NULL),
182 fCutString(NULL),
0a2b2b4b 183 fUtils(NULL),
11c1e680 184 fEtaShift(0.0),
185 fDoEtaShift(kFALSE),
ae947965 186 fDoReweightHistoMCPi0(kFALSE),
187 fDoReweightHistoMCEta(kFALSE),
188 fDoReweightHistoMCK0s(kFALSE),
189 fPathTrFReweighting(""),
190 fNameHistoReweightingPi0(""),
191 fNameHistoReweightingEta(""),
192 fNameHistoReweightingK0s(""),
4a0aab28 193 hdEdxCuts(NULL),
194 hTPCdEdxbefore(NULL),
195 hTPCdEdxafter(NULL),
e5b6e8a6 196 hTPCdEdxSigbefore(NULL),
197 hTPCdEdxSigafter(NULL),
4a0aab28 198 hTOFbefore(NULL),
e5b6e8a6 199 hTOFSigbefore(NULL),
200 hTOFSigafter(NULL),
4a0aab28 201 hTrackCuts(NULL),
202 hPhotonCuts(NULL),
203 hInvMassbefore(NULL),
204 hArmenterosbefore(NULL),
205 hInvMassafter(NULL),
206 hArmenterosafter(NULL),
207 hAcceptanceCuts(NULL),
208 hCutIndex(NULL),
209 hV0EventCuts(NULL),
210 hCentrality(NULL),
4803eb1f 211 hCentralityVsNumberOfPrimaryTracks(NULL),
4a0aab28 212 hVertexZ(NULL),
a280ac15 213 hTriggerClass(NULL),
ae947965 214 hTriggerClassSelected(NULL),
215 hReweightMCHistPi0(NULL),
216 hReweightMCHistEta(NULL),
217 hReweightMCHistK0s(NULL),
218 fPreSelCut(kFALSE),
219 fTriggerSelectedManually(kFALSE),
220 fSpecialTriggerName("")
221
3b77b2d1 222{
4a0aab28 223 InitPIDResponse();
224 for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=0;}
225 fCutString=new TObjString((GetCutNumber()).Data());
3b77b2d1 226
e5b6e8a6 227 fElectronLabelArray = new Int_t[fElectronArraySize];
0a2b2b4b 228 fUtils = new AliAnalysisUtils();
229 //if you do not want to apply the cut on the distance between the SPD and TRK vertex:
11c1e680 230 //fUtils->SetCutOnZVertexSPD(kFALSE);
231
ca91a3e1 232
3b77b2d1 233}
234
1d9e6011 235//________________________________________________________________________
236AliConversionCuts::AliConversionCuts(const AliConversionCuts &ref) :
237 AliAnalysisCuts(ref),
238 fHistograms(NULL),
239 fHeaderList(ref.fHeaderList),
240 fPIDResponse(NULL),
241 fEventQuality(ref.fEventQuality),
242 fMaxR(ref.fMaxR),
243 fMinR(ref.fMinR),
244 fEtaCut(ref.fEtaCut),
245 fEtaCutMin(ref.fEtaCutMin),
246 fPtCut(ref.fPtCut),
247 fSinglePtCut(ref.fSinglePtCut),
248 fMaxZ(ref.fMaxZ),
249 fMinClsTPC(ref.fMinClsTPC),
250 fMinClsTPCToF(ref.fMinClsTPCToF),
251 fLineCutZRSlope(ref.fLineCutZRSlope),
252 fLineCutZValue(ref.fLineCutZValue),
253 fLineCutZRSlopeMin(ref.fLineCutZRSlopeMin),
254 fLineCutZValueMin(ref.fLineCutZValueMin),
255 fChi2CutConversion(ref.fChi2CutConversion),
256 fPIDProbabilityCutNegativeParticle(ref.fPIDProbabilityCutNegativeParticle),
257 fPIDProbabilityCutPositiveParticle(ref.fPIDProbabilityCutPositiveParticle),
258 fDodEdxSigmaCut(ref. fDodEdxSigmaCut),
259 fDoTOFsigmaCut(ref.fDoTOFsigmaCut),
260 fPIDTRDEfficiency(ref.fPIDTRDEfficiency),
261 fDoTRDPID(ref.fDoTRDPID),
262 fPIDnSigmaAboveElectronLine(ref.fPIDnSigmaAboveElectronLine),
263 fPIDnSigmaBelowElectronLine(ref.fPIDnSigmaBelowElectronLine),
264 fTofPIDnSigmaAboveElectronLine(ref.fTofPIDnSigmaAboveElectronLine),
265 fTofPIDnSigmaBelowElectronLine(ref.fTofPIDnSigmaBelowElectronLine),
266 fPIDnSigmaAbovePionLine(ref.fPIDnSigmaAbovePionLine),
267 fPIDnSigmaAbovePionLineHighPt(ref.fPIDnSigmaAbovePionLineHighPt),
268 fPIDMinPnSigmaAbovePionLine(ref.fPIDMinPnSigmaAbovePionLine),
269 fPIDMaxPnSigmaAbovePionLine(ref.fPIDMaxPnSigmaAbovePionLine),
270 fDoKaonRejectionLowP(ref.fDoKaonRejectionLowP),
271 fDoProtonRejectionLowP(ref.fDoProtonRejectionLowP),
272 fDoPionRejectionLowP(ref.fDoPionRejectionLowP),
273 fPIDnSigmaAtLowPAroundKaonLine(ref.fPIDnSigmaAtLowPAroundKaonLine),
274 fPIDnSigmaAtLowPAroundProtonLine(ref.fPIDnSigmaAtLowPAroundProtonLine),
275 fPIDnSigmaAtLowPAroundPionLine(ref.fPIDnSigmaAtLowPAroundPionLine),
276 fPIDMinPKaonRejectionLowP(ref.fPIDMinPKaonRejectionLowP),
277 fPIDMinPProtonRejectionLowP(ref.fPIDMinPProtonRejectionLowP),
278 fPIDMinPPionRejectionLowP(ref.fPIDMinPPionRejectionLowP),
279 fDoQtGammaSelection(ref.fDoQtGammaSelection),
280 fDoHighPtQtGammaSelection(ref.fDoHighPtQtGammaSelection),
281 fQtMax(ref.fQtMax),
282 fHighPtQtMax(ref.fHighPtQtMax),
283 fPtBorderForQt(ref.fPtBorderForQt),
284 fXVertexCut(ref.fXVertexCut),
285 fYVertexCut(ref.fYVertexCut),
286 fZVertexCut(ref.fZVertexCut),
287 fNSigmaMass(ref.fNSigmaMass),
288 fUseEtaMinCut(ref.fUseEtaMinCut),
289 fUseOnFlyV0Finder(ref.fUseOnFlyV0Finder),
0a2b2b4b 290 fDoPhotonAsymmetryCut(ref.fDoPhotonAsymmetryCut),
1d9e6011 291 fMinPPhotonAsymmetryCut(ref.fMinPPhotonAsymmetryCut),
292 fMinPhotonAsymmetry(ref.fMinPhotonAsymmetry),
293 fIsHeavyIon(ref.fIsHeavyIon),
294 fDetectorCentrality(ref.fDetectorCentrality),
295 fModCentralityClass(ref.fModCentralityClass),
296 fMaxVertexZ(ref.fMaxVertexZ),
297 fCentralityMin(ref.fCentralityMin),
298 fCentralityMax(ref.fCentralityMax),
299 fUseCorrectedTPCClsInfo(ref.fUseCorrectedTPCClsInfo),
300 fUseTOFpid(ref.fUseTOFpid),
301 fMultiplicityMethod(ref.fMultiplicityMethod),
302 fSpecialTrigger(ref.fSpecialTrigger),
303 fRemovePileUp(ref.fRemovePileUp),
304 fOpeningAngle(ref.fOpeningAngle),
305 fPsiPairCut(ref.fPsiPairCut),
306 fCosPAngleCut(ref.fCosPAngleCut),
307 fDoToCloseV0sCut(ref.fDoToCloseV0sCut),
308 fRejectExtraSignals(ref.fRejectExtraSignals),
309 fminV0Dist(ref.fminV0Dist),
310 fDoSharedElecCut(ref.fDoSharedElecCut),
311 fOfflineTriggerMask(ref.fOfflineTriggerMask),
312 fHasV0AND(ref.fHasV0AND),
313 fIsSDDFired(ref.fIsSDDFired),
314 fRandom(ref.fRandom),
315 fElectronArraySize(ref.fElectronArraySize),
316 fElectronLabelArray(NULL),
4803eb1f 317 fDCAZPrimVtxCut(ref.fDCAZPrimVtxCut),
318 fDCARPrimVtxCut(ref.fDCAZPrimVtxCut),
1d9e6011 319 fConversionPointXArray(ref.fConversionPointXArray),
320 fConversionPointYArray(ref.fConversionPointYArray),
321 fConversionPointZArray(ref.fConversionPointZArray),
322 fnHeaders(ref.fnHeaders),
323 fNotRejectedStart(NULL),
324 fNotRejectedEnd(NULL),
325 fGeneratorNames(ref.fGeneratorNames),
326 fCutString(NULL),
0a2b2b4b 327 fUtils(NULL),
11c1e680 328 fEtaShift(ref.fEtaShift),
329 fDoEtaShift(ref.fDoEtaShift),
ae947965 330 fDoReweightHistoMCPi0(ref.fDoReweightHistoMCPi0),
331 fDoReweightHistoMCEta(ref.fDoReweightHistoMCEta),
332 fDoReweightHistoMCK0s(ref.fDoReweightHistoMCK0s),
333 fPathTrFReweighting(ref.fPathTrFReweighting),
334 fNameHistoReweightingPi0(ref.fNameHistoReweightingPi0),
335 fNameHistoReweightingEta(ref.fNameHistoReweightingEta),
336 fNameHistoReweightingK0s(ref.fNameHistoReweightingK0s),
1d9e6011 337 hdEdxCuts(NULL),
338 hTPCdEdxbefore(NULL),
339 hTPCdEdxafter(NULL),
340 hTPCdEdxSigbefore(NULL),
341 hTPCdEdxSigafter(NULL),
342 hTOFbefore(NULL),
343 hTOFSigbefore(NULL),
344 hTOFSigafter(NULL),
345 hTrackCuts(NULL),
346 hPhotonCuts(NULL),
347 hInvMassbefore(NULL),
348 hArmenterosbefore(NULL),
349 hInvMassafter(NULL),
350 hArmenterosafter(NULL),
351 hAcceptanceCuts(NULL),
352 hCutIndex(NULL),
353 hV0EventCuts(NULL),
354 hCentrality(NULL),
4803eb1f 355 hCentralityVsNumberOfPrimaryTracks(NULL),
1d9e6011 356 hVertexZ(NULL),
a280ac15 357 hTriggerClass(NULL),
ae947965 358 hTriggerClassSelected(NULL),
359 hReweightMCHistPi0(NULL),
360 hReweightMCHistEta(NULL),
361 hReweightMCHistK0s(NULL),
362 fPreSelCut(ref.fPreSelCut),
363 fTriggerSelectedManually(ref.fTriggerSelectedManually),
364 fSpecialTriggerName(ref.fSpecialTriggerName)
1d9e6011 365{
366 // Copy Constructor
a280ac15 367 for(Int_t jj=0;jj<kNCuts;jj++){fCuts[jj]=ref.fCuts[jj];}
368 fCutString=new TObjString((GetCutNumber()).Data());
369 fElectronLabelArray = new Int_t[fElectronArraySize];
0a2b2b4b 370 fUtils = new AliAnalysisUtils();
a280ac15 371 // dont copy histograms (if you like histograms, call InitCutHistograms())
11c1e680 372
1d9e6011 373}
374
a280ac15 375
3b77b2d1 376//________________________________________________________________________
377AliConversionCuts::~AliConversionCuts() {
4a0aab28 378 // Destructor
379 //Deleting fHistograms leads to seg fault it it's added to output collection of a task
380 // if(fHistograms)
381 // delete fHistograms;
382 // fHistograms = NULL;
92efd725 383 if(fCutString != NULL){
384 delete fCutString;
385 fCutString = NULL;
386 }
92efd725 387 if(fElectronLabelArray){
388 delete fElectronLabelArray;
389 fElectronLabelArray = NULL;
390 }
ca91a3e1 391 if(fNotRejectedStart){
392 delete[] fNotRejectedStart;
393 fNotRejectedStart = NULL;
394 }
395 if(fNotRejectedEnd){
396 delete[] fNotRejectedEnd;
397 fNotRejectedEnd = NULL;
398 }
399 if(fGeneratorNames){
400 delete[] fGeneratorNames;
401 fGeneratorNames = NULL;
402 }
0a2b2b4b 403 if(fUtils){
404 delete fUtils;
405 fUtils = NULL;
406 }
407
3b77b2d1 408}
409
410//________________________________________________________________________
92efd725 411void AliConversionCuts::InitCutHistograms(TString name, Bool_t preCut){
3b77b2d1 412
4a0aab28 413 // Initialize Cut Histograms for QA (only initialized and filled if function is called)
ccfa8c0d 414 TH1::AddDirectory(kFALSE);
4a0aab28 415
416 if(fHistograms != NULL){
417 delete fHistograms;
418 fHistograms=NULL;
419 }
420 if(fHistograms==NULL){
421 fHistograms=new TList();
ccfa8c0d 422 fHistograms->SetOwner(kTRUE);
4a0aab28 423 if(name=="")fHistograms->SetName(Form("ConvCuts_%s",GetCutNumber().Data()));
424 else fHistograms->SetName(Form("%s_%s",name.Data(),GetCutNumber().Data()));
425 }
426
ae947965 427 if (hReweightMCHistPi0) fHistograms->Add(hReweightMCHistPi0);
428 if (hReweightMCHistEta) fHistograms->Add(hReweightMCHistEta);
429 if (hReweightMCHistK0s) fHistograms->Add(hReweightMCHistK0s);
430
4a0aab28 431 // IsPhotonSelected
432 hCutIndex=new TH1F(Form("IsPhotonSelected %s",GetCutNumber().Data()),"IsPhotonSelected",10,-0.5,9.5);
433 hCutIndex->GetXaxis()->SetBinLabel(kPhotonIn+1,"in");
434 hCutIndex->GetXaxis()->SetBinLabel(kOnFly+1,"onfly");
435 hCutIndex->GetXaxis()->SetBinLabel(kNoTracks+1,"no tracks");
436 hCutIndex->GetXaxis()->SetBinLabel(kdEdxCuts+1,"dEdx");
437 hCutIndex->GetXaxis()->SetBinLabel(kTrackCuts+1,"Track cuts");
438 hCutIndex->GetXaxis()->SetBinLabel(kConvPointFail+1,"ConvPoint fail");
439 hCutIndex->GetXaxis()->SetBinLabel(kPhotonCuts+1,"PhotonCuts");
440 hCutIndex->GetXaxis()->SetBinLabel(kPhotonOut+1,"out");
441 fHistograms->Add(hCutIndex);
442
443 // Track Cuts
e5b6e8a6 444 hTrackCuts=new TH1F(Form("TrackCuts %s",GetCutNumber().Data()),"TrackCuts",9,-0.5,8.5);
4a0aab28 445 hTrackCuts->GetXaxis()->SetBinLabel(1,"in");
446 hTrackCuts->GetXaxis()->SetBinLabel(2,"likesign");
447 hTrackCuts->GetXaxis()->SetBinLabel(3,"ntpccl");
448 hTrackCuts->GetXaxis()->SetBinLabel(4,"acceptance");
449 hTrackCuts->GetXaxis()->SetBinLabel(5,"singlept");
450 hTrackCuts->GetXaxis()->SetBinLabel(6,"TPCrefit");
451 hTrackCuts->GetXaxis()->SetBinLabel(7,"kink");
452 hTrackCuts->GetXaxis()->SetBinLabel(8,"out");
453 fHistograms->Add(hTrackCuts);
454
455 // Photon Cuts
4803eb1f 456 hPhotonCuts=new TH1F(Form("PhotonCuts %s",GetCutNumber().Data()),"PhotonCuts",14,-0.5,13.5);
4a0aab28 457 hPhotonCuts->GetXaxis()->SetBinLabel(1,"in");
458 hPhotonCuts->GetXaxis()->SetBinLabel(2,"qtcut");
459 hPhotonCuts->GetXaxis()->SetBinLabel(3,"chi2");
460 hPhotonCuts->GetXaxis()->SetBinLabel(4,"acceptance");
461 hPhotonCuts->GetXaxis()->SetBinLabel(5,"asymmetry");
462 hPhotonCuts->GetXaxis()->SetBinLabel(6,"pidprob");
463 hPhotonCuts->GetXaxis()->SetBinLabel(7,"cortpcclinfo");
464 hPhotonCuts->GetXaxis()->SetBinLabel(8,"PsiPair");
465 hPhotonCuts->GetXaxis()->SetBinLabel(9,"CosPAngle");
4803eb1f 466 hPhotonCuts->GetXaxis()->SetBinLabel(10,"DCA R");
467 hPhotonCuts->GetXaxis()->SetBinLabel(11,"DCA Z");
468 hPhotonCuts->GetXaxis()->SetBinLabel(12,"out");
4a0aab28 469 fHistograms->Add(hPhotonCuts);
470
471 if(preCut){
e5b6e8a6 472 hInvMassbefore=new TH1F(Form("InvMass_before %s",GetCutNumber().Data()),"InvMass_before",1000,0,0.3);
4a0aab28 473 fHistograms->Add(hInvMassbefore);
e5b6e8a6 474 hArmenterosbefore=new TH2F(Form("Armenteros_before %s",GetCutNumber().Data()),"Armenteros_before",200,-1,1,1000,0,1.);
4a0aab28 475 fHistograms->Add(hArmenterosbefore);
476 }
e5b6e8a6 477 hInvMassafter=new TH1F(Form("InvMass_after %s",GetCutNumber().Data()),"InvMass_after",1000,0,0.3);
4a0aab28 478 fHistograms->Add(hInvMassafter);
479 hArmenterosafter=new TH2F(Form("Armenteros_after %s",GetCutNumber().Data()),"Armenteros_after",200,-1,1,250,0,0.25);
480 fHistograms->Add(hArmenterosafter);
481
482 hAcceptanceCuts=new TH1F(Form("PhotonAcceptanceCuts %s",GetCutNumber().Data()),"PhotonAcceptanceCuts",10,-0.5,9.5);
483 hAcceptanceCuts->GetXaxis()->SetBinLabel(1,"in");
484 hAcceptanceCuts->GetXaxis()->SetBinLabel(2,"maxR");
485 hAcceptanceCuts->GetXaxis()->SetBinLabel(3,"minR");
486 hAcceptanceCuts->GetXaxis()->SetBinLabel(4,"line");
487 hAcceptanceCuts->GetXaxis()->SetBinLabel(5,"maxZ");
488 hAcceptanceCuts->GetXaxis()->SetBinLabel(6,"eta");
489 hAcceptanceCuts->GetXaxis()->SetBinLabel(7,"minpt");
490 hAcceptanceCuts->GetXaxis()->SetBinLabel(8,"out");
491 fHistograms->Add(hAcceptanceCuts);
492
493 // dEdx Cuts
494 hdEdxCuts=new TH1F(Form("dEdxCuts %s",GetCutNumber().Data()),"dEdxCuts",10,-0.5,9.5);
495 hdEdxCuts->GetXaxis()->SetBinLabel(1,"in");
496 hdEdxCuts->GetXaxis()->SetBinLabel(2,"TPCelectron");
497 hdEdxCuts->GetXaxis()->SetBinLabel(3,"TPCpion");
498 hdEdxCuts->GetXaxis()->SetBinLabel(4,"TPCpionhighp");
499 hdEdxCuts->GetXaxis()->SetBinLabel(5,"TPCkaonlowprej");
500 hdEdxCuts->GetXaxis()->SetBinLabel(6,"TPCprotonlowprej");
501 hdEdxCuts->GetXaxis()->SetBinLabel(7,"TPCpionlowprej");
502 hdEdxCuts->GetXaxis()->SetBinLabel(8,"TOFelectron");
503 hdEdxCuts->GetXaxis()->SetBinLabel(9,"TRDelectron");
504 hdEdxCuts->GetXaxis()->SetBinLabel(10,"out");
505 fHistograms->Add(hdEdxCuts);
e5b6e8a6 506
4a0aab28 507 TAxis *AxisBeforedEdx = NULL;
e5b6e8a6 508 TAxis *AxisBeforedEdxSig = NULL;
4a0aab28 509 TAxis *AxisBeforeTOF = NULL;
e5b6e8a6 510 TAxis *AxisBeforeTOFSig = NULL;
4a0aab28 511 if(preCut){
e5b6e8a6 512 hTPCdEdxbefore=new TH2F(Form("Gamma_dEdx_before %s",GetCutNumber().Data()),"dEdx Gamma before" ,150,0.03,20,800,0,200);
4a0aab28 513 fHistograms->Add(hTPCdEdxbefore);
e5b6e8a6 514 AxisBeforedEdx = hTPCdEdxbefore->GetXaxis();
515 hTPCdEdxSigbefore=new TH2F(Form("Gamma_dEdxSig_before %s",GetCutNumber().Data()),"dEdx Sigma Gamma before" ,150,0.03,20,400,-10,10);
516 fHistograms->Add(hTPCdEdxSigbefore);
517 AxisBeforedEdxSig = hTPCdEdxSigbefore->GetXaxis();
4a0aab28 518
e5b6e8a6 519 hTOFbefore=new TH2F(Form("Gamma_TOF_before %s",GetCutNumber().Data()),"TOF Gamma before" ,150,0.03,20,11000,-1000,10000);
4a0aab28 520 fHistograms->Add(hTOFbefore);
e5b6e8a6 521 AxisBeforeTOF = hTOFbefore->GetXaxis();
522 hTOFSigbefore=new TH2F(Form("Gamma_TOFSig_before %s",GetCutNumber().Data()),"TOF Sigma Gamma before" ,150,0.03,20,400,-6,10);
523 fHistograms->Add(hTOFSigbefore);
524 AxisBeforeTOFSig = hTOFSigbefore->GetXaxis();
525
4a0aab28 526 }
e5b6e8a6 527 hTPCdEdxSigafter=new TH2F(Form("Gamma_dEdxSig_after %s",GetCutNumber().Data()),"dEdx Sigma Gamma after" ,150,0.03,20,400, -10,10);
528 fHistograms->Add(hTPCdEdxSigafter);
529
530 hTPCdEdxafter=new TH2F(Form("Gamma_dEdx_after %s",GetCutNumber().Data()),"dEdx Gamma after" ,150,0.03,20,800,0,200);
4a0aab28 531 fHistograms->Add(hTPCdEdxafter);
532
e5b6e8a6 533 hTOFSigafter=new TH2F(Form("Gamma_TOFSig_after %s",GetCutNumber().Data()),"TOF Sigma Gamma after" ,150,0.03,20,400,-6,10);
534 fHistograms->Add(hTOFSigafter);
4a0aab28 535
e5b6e8a6 536 TAxis *AxisAfter = hTPCdEdxSigafter->GetXaxis();
4a0aab28 537 Int_t bins = AxisAfter->GetNbins();
538 Double_t from = AxisAfter->GetXmin();
539 Double_t to = AxisAfter->GetXmax();
540 Double_t *newBins = new Double_t[bins+1];
541 newBins[0] = from;
542 Double_t factor = TMath::Power(to/from, 1./bins);
543 for(Int_t i=1; i<=bins; ++i) newBins[i] = factor * newBins[i-1];
544 AxisAfter->Set(bins, newBins);
e5b6e8a6 545 AxisAfter = hTOFSigafter->GetXaxis();
546 AxisAfter->Set(bins, newBins);
547 AxisAfter = hTPCdEdxafter->GetXaxis();
4a0aab28 548 AxisAfter->Set(bins, newBins);
549 if(preCut){
550 AxisBeforedEdx->Set(bins, newBins);
551 AxisBeforeTOF->Set(bins, newBins);
e5b6e8a6 552 AxisBeforedEdxSig->Set(bins, newBins);
553 AxisBeforeTOFSig->Set(bins, newBins);
4a0aab28 554 }
555 delete [] newBins;
e5b6e8a6 556
4a0aab28 557 // Event Cuts and Info
558 if(preCut){
e5b6e8a6 559 hV0EventCuts=new TH1F(Form("ESD_EventCuts %s",GetCutNumber().Data()),"Event Cuts",7,-0.5,6.5);
4a0aab28 560 hV0EventCuts->GetXaxis()->SetBinLabel(1,"in");
561 hV0EventCuts->GetXaxis()->SetBinLabel(2,"OfflineTrigger");
e5b6e8a6 562 hV0EventCuts->GetXaxis()->SetBinLabel(3,"nvtxcontr");
563 hV0EventCuts->GetXaxis()->SetBinLabel(4,"VertexZ");
4a0aab28 564 hV0EventCuts->GetXaxis()->SetBinLabel(5,"pileup");
565 hV0EventCuts->GetXaxis()->SetBinLabel(6,"centrsel");
566 hV0EventCuts->GetXaxis()->SetBinLabel(7,"out");
567 fHistograms->Add(hV0EventCuts);
e5b6e8a6 568
4803eb1f 569 hCentrality=new TH1F(Form("Centrality %s",GetCutNumber().Data()),"Centrality",1000,0,100);
4a0aab28 570 fHistograms->Add(hCentrality);
4803eb1f 571 hCentralityVsNumberOfPrimaryTracks=new TH2F(Form("Centrality vs Primary Tracks %s",GetCutNumber().Data()),"Centrality vs Primary Tracks ",100,0,100,3000,0,3000);
572 fHistograms->Add(hCentralityVsNumberOfPrimaryTracks);
4a0aab28 573 hVertexZ=new TH1F(Form("VertexZ %s",GetCutNumber().Data()),"VertexZ",1000,-50,50);
574 fHistograms->Add(hVertexZ);
e5b6e8a6 575
11c1e680 576 hTriggerClass= new TH1F(Form("OfflineTrigger %s",GetCutNumber().Data()),"OfflineTrigger",35,-0.5,34.5);
e5b6e8a6 577 hTriggerClass->GetXaxis()->SetBinLabel( 1,"kMB");
578 hTriggerClass->GetXaxis()->SetBinLabel( 2,"kINT7");
579 hTriggerClass->GetXaxis()->SetBinLabel( 3,"kMUON");
580 hTriggerClass->GetXaxis()->SetBinLabel( 4,"kHighMult");
581 hTriggerClass->GetXaxis()->SetBinLabel( 5,"kKEMC1");
582 hTriggerClass->GetXaxis()->SetBinLabel( 6,"kCINT5");
583 hTriggerClass->GetXaxis()->SetBinLabel( 7,"kCMUS5/kMUSPB");
584 hTriggerClass->GetXaxis()->SetBinLabel( 8,"kMUSH7/kMUSHPB");
585 hTriggerClass->GetXaxis()->SetBinLabel( 9,"kMUL7/kMuonLikePB");
586 hTriggerClass->GetXaxis()->SetBinLabel(10,"kMUU7/kMuonUnlikePB");
587 hTriggerClass->GetXaxis()->SetBinLabel(11,"kEMC7/kEMC8");
588 hTriggerClass->GetXaxis()->SetBinLabel(12,"kMUS7");
589 hTriggerClass->GetXaxis()->SetBinLabel(13,"kPHI1");
590 hTriggerClass->GetXaxis()->SetBinLabel(14,"kPHI7/kPHI8/kPHOSPb");
591 hTriggerClass->GetXaxis()->SetBinLabel(15,"kEMCEJE");
592 hTriggerClass->GetXaxis()->SetBinLabel(16,"kEMCEGA");
593 hTriggerClass->GetXaxis()->SetBinLabel(17,"kCentral");
594 hTriggerClass->GetXaxis()->SetBinLabel(18,"kSemiCentral");
595 hTriggerClass->GetXaxis()->SetBinLabel(19,"kDG5");
596 hTriggerClass->GetXaxis()->SetBinLabel(20,"kZED");
597 hTriggerClass->GetXaxis()->SetBinLabel(21,"kSPI7/kSPI");
598 hTriggerClass->GetXaxis()->SetBinLabel(22,"kINT8");
599 hTriggerClass->GetXaxis()->SetBinLabel(23,"kMuonSingleLowPt8");
600 hTriggerClass->GetXaxis()->SetBinLabel(24,"kMuonSingleHighPt8");
601 hTriggerClass->GetXaxis()->SetBinLabel(25,"kMuonLikeLowPt8");
602 hTriggerClass->GetXaxis()->SetBinLabel(26,"kMuonUnlikeLowPt8");
603 hTriggerClass->GetXaxis()->SetBinLabel(27,"kMuonUnlikeLowPt0");
604 hTriggerClass->GetXaxis()->SetBinLabel(28,"kUserDefined");
11c1e680 605 hTriggerClass->GetXaxis()->SetBinLabel(29,"kTRD");
606 hTriggerClass->GetXaxis()->SetBinLabel(30,"kFastOnly");
607 hTriggerClass->GetXaxis()->SetBinLabel(31,"kAnyINT");
608 hTriggerClass->GetXaxis()->SetBinLabel(32,"kAny");
609 hTriggerClass->GetXaxis()->SetBinLabel(33,"V0AND");
610 hTriggerClass->GetXaxis()->SetBinLabel(34,"NOT kFastOnly");
611 hTriggerClass->GetXaxis()->SetBinLabel(35,"failed Physics Selection");
4a0aab28 612 fHistograms->Add(hTriggerClass);
ae947965 613 }
614 if(!preCut){
11c1e680 615 hTriggerClassSelected= new TH1F(Form("OfflineTriggerSelected %s",GetCutNumber().Data()),"OfflineTriggerSelected",34,-0.5,33.5);
a280ac15 616 hTriggerClassSelected->GetXaxis()->SetBinLabel( 1,"kMB");
617 hTriggerClassSelected->GetXaxis()->SetBinLabel( 2,"kINT7");
618 hTriggerClassSelected->GetXaxis()->SetBinLabel( 3,"kMUON");
619 hTriggerClassSelected->GetXaxis()->SetBinLabel( 4,"kHighMult");
620 hTriggerClassSelected->GetXaxis()->SetBinLabel( 5,"kKEMC1");
621 hTriggerClassSelected->GetXaxis()->SetBinLabel( 6,"kCINT5");
622 hTriggerClassSelected->GetXaxis()->SetBinLabel( 7,"kCMUS5/kMUSPB");
623 hTriggerClassSelected->GetXaxis()->SetBinLabel( 8,"kMUSH7/kMUSHPB");
624 hTriggerClassSelected->GetXaxis()->SetBinLabel( 9,"kMUL7/kMuonLikePB");
625 hTriggerClassSelected->GetXaxis()->SetBinLabel(10,"kMUU7/kMuonUnlikePB");
626 hTriggerClassSelected->GetXaxis()->SetBinLabel(11,"kEMC7/kEMC8");
627 hTriggerClassSelected->GetXaxis()->SetBinLabel(12,"kMUS7");
628 hTriggerClassSelected->GetXaxis()->SetBinLabel(13,"kPHI1");
629 hTriggerClassSelected->GetXaxis()->SetBinLabel(14,"kPHI7/kPHI8/kPHOSPb");
630 hTriggerClassSelected->GetXaxis()->SetBinLabel(15,"kEMCEJE");
631 hTriggerClassSelected->GetXaxis()->SetBinLabel(16,"kEMCEGA");
632 hTriggerClassSelected->GetXaxis()->SetBinLabel(17,"kCentral");
633 hTriggerClassSelected->GetXaxis()->SetBinLabel(18,"kSemiCentral");
634 hTriggerClassSelected->GetXaxis()->SetBinLabel(19,"kDG5");
635 hTriggerClassSelected->GetXaxis()->SetBinLabel(20,"kZED");
636 hTriggerClassSelected->GetXaxis()->SetBinLabel(21,"kSPI7/kSPI");
637 hTriggerClassSelected->GetXaxis()->SetBinLabel(22,"kINT8");
638 hTriggerClassSelected->GetXaxis()->SetBinLabel(23,"kMuonSingleLowPt8");
639 hTriggerClassSelected->GetXaxis()->SetBinLabel(24,"kMuonSingleHighPt8");
640 hTriggerClassSelected->GetXaxis()->SetBinLabel(25,"kMuonLikeLowPt8");
641 hTriggerClassSelected->GetXaxis()->SetBinLabel(26,"kMuonUnlikeLowPt8");
642 hTriggerClassSelected->GetXaxis()->SetBinLabel(27,"kMuonUnlikeLowPt0");
643 hTriggerClassSelected->GetXaxis()->SetBinLabel(28,"kUserDefined");
11c1e680 644 hTriggerClassSelected->GetXaxis()->SetBinLabel(29,"kTRD");
645 hTriggerClassSelected->GetXaxis()->SetBinLabel(30,"kFastOnly");
646 hTriggerClassSelected->GetXaxis()->SetBinLabel(31,"kAnyINT");
647 hTriggerClassSelected->GetXaxis()->SetBinLabel(32,"kAny");
648 hTriggerClassSelected->GetXaxis()->SetBinLabel(33,"V0AND");
649 hTriggerClassSelected->GetXaxis()->SetBinLabel(34,"NOT kFastOnly");
a280ac15 650 fHistograms->Add(hTriggerClassSelected);
4a0aab28 651 }
ccfa8c0d 652 TH1::AddDirectory(kTRUE);
3b77b2d1 653}
654
655//________________________________________________________________________
656Bool_t AliConversionCuts::InitPIDResponse(){
4a0aab28 657 // Set Pointer to AliPIDResponse
3b77b2d1 658
4a0aab28 659 AliAnalysisManager *man=AliAnalysisManager::GetAnalysisManager();
e5b6e8a6 660 if(man) {
4a0aab28 661 AliInputEventHandler* inputHandler = (AliInputEventHandler*) (man->GetInputEventHandler());
662 fPIDResponse = (AliPIDResponse*)inputHandler->GetPIDResponse();
663 if(fPIDResponse)return kTRUE;
e5b6e8a6 664
4a0aab28 665 }
e5b6e8a6 666
ca91a3e1 667
4a0aab28 668 return kFALSE;
3b77b2d1 669}
670///________________________________________________________________________
92efd725 671Bool_t AliConversionCuts::EventIsSelected(AliVEvent *fInputEvent, AliVEvent *fMCEvent){
4a0aab28 672 // Process Event Selection
3b77b2d1 673
4a0aab28 674 Int_t cutindex=0;
675 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
676 cutindex++;
e5b6e8a6 677
4a0aab28 678 // Check for MC event
ae947965 679 if(fMCEvent && fInputEvent->IsA()==AliESDEvent::Class()){
4a0aab28 680 // Check if MC event is correctly loaded
681 AliMCEventHandler* mcHandler = dynamic_cast<AliMCEventHandler*> (AliAnalysisManager::GetAnalysisManager()->GetMCtruthEventHandler());
682 if (!mcHandler){
683 fEventQuality = 2;
684 return kFALSE;
685 }
686 if (!mcHandler->InitOk() ){
687 fEventQuality = 2;
688 return kFALSE;
689 }
690 if (!mcHandler->TreeK() ){
691 fEventQuality = 2;
692 return kFALSE;
693 }
694 if (!mcHandler->TreeTR() ) {
695 fEventQuality = 2;
696 return kFALSE;
697 }
698 }
4a0aab28 699
e5b6e8a6 700 // Event Trigger
4803eb1f 701// cout << "before event trigger" << endl;
1186afd2 702 if(!IsTriggerSelected(fInputEvent)){
4a0aab28 703 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
704 fEventQuality = 3;
705 return kFALSE;
706 }
707 cutindex++;
a280ac15 708
e5b6e8a6 709 if(fInputEvent->IsA()==AliESDEvent::Class()){
710 AliTriggerAnalysis fTriggerAnalysis;// = new AliTriggerAnalysis;
711 fHasV0AND = fTriggerAnalysis.IsOfflineTriggerFired((AliESDEvent*)fInputEvent, AliTriggerAnalysis::kV0AND);
11c1e680 712 if(fHasV0AND&&hTriggerClass)hTriggerClass->Fill(32);
e5b6e8a6 713 }
ccfa8c0d 714// cout << "event number " << ((AliESDEvent*)fInputEvent)->GetEventNumberInFile() << " entered"<< endl;
715
4a0aab28 716
e5b6e8a6 717 // Number of Contributors Cut
718 if(GetNumberOfContributorsVtx(fInputEvent)<=0) {
4a0aab28 719 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
e5b6e8a6 720 fEventQuality = 5;
4a0aab28 721 return kFALSE;
722 }
723 cutindex++;
724
e5b6e8a6 725 // Z Vertex Position Cut
726 if(!VertexZCut(fInputEvent)){
4a0aab28 727 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
e5b6e8a6 728 fEventQuality = 4;
4a0aab28 729 return kFALSE;
730 }
731 cutindex++;
92efd725 732
4a0aab28 733 // Pile Up Rejection
3b77b2d1 734
4a0aab28 735 if(fRemovePileUp){
736 if(fInputEvent->IsPileupFromSPD(3,0.8,3.,2.,5.)){
737 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
738 fEventQuality = 6;
739 return kFALSE;
740 }
741 }
742 cutindex++;
743
744 // Centrality Selection
e5b6e8a6 745 if(!IsCentralitySelected(fInputEvent,fMCEvent)){
4a0aab28 746 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
747 fEventQuality = 1;
748 return kFALSE;
749 }
750 cutindex++;
751
752 // Fill Event Histograms
753 if(hV0EventCuts)hV0EventCuts->Fill(cutindex);
754 if(hVertexZ)hVertexZ->Fill(fInputEvent->GetPrimaryVertex()->GetZ());
755 if(hCentrality)hCentrality->Fill(GetCentrality(fInputEvent));
4803eb1f 756 if(hCentralityVsNumberOfPrimaryTracks)
757 hCentralityVsNumberOfPrimaryTracks->Fill(GetCentrality(fInputEvent),
758 ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()
759 ->GetTask("V0ReaderV1"))->GetNumberOfPrimaryTracks());
4a0aab28 760 fEventQuality = 0;
761 return kTRUE;
3b77b2d1 762}
763
764///________________________________________________________________________
92efd725 765Bool_t AliConversionCuts::PhotonIsSelectedMC(TParticle *particle,AliStack *fMCStack,Bool_t checkForConvertedGamma){
4a0aab28 766 // MonteCarlo Photon Selection
3b77b2d1 767
4a0aab28 768 if(!fMCStack)return kFALSE;
3b77b2d1 769
4a0aab28 770 if (particle->GetPdgCode() == 22){
e5b6e8a6 771
4a0aab28 772
11c1e680 773 if( particle->Eta() > (fEtaCut + fEtaShift) || particle->Eta() < (-fEtaCut + fEtaShift) )
774 return kFALSE;
775 if(fEtaCutMin>-0.1){
776 if( particle->Eta() < (fEtaCutMin + fEtaShift) && particle->Eta() > (-fEtaCutMin + fEtaShift) )
777 return kFALSE;
778 }
ccfa8c0d 779
4a0aab28 780 if(particle->GetMother(0) >-1 && fMCStack->Particle(particle->GetMother(0))->GetPdgCode() == 22){
781 return kFALSE; // no photon as mothers!
782 }
783
784 if(particle->GetMother(0) >= fMCStack->GetNprimary()){
785 return kFALSE; // the gamma has a mother, and it is not a primary particle
786 }
787
788 if(!checkForConvertedGamma) return kTRUE; // return in case of accepted gamma
789
790 // looking for conversion gammas (electron + positron from pairbuilding (= 5) )
791 TParticle* ePos = NULL;
792 TParticle* eNeg = NULL;
793
794 if(particle->GetNDaughters() >= 2){
795 for(Int_t daughterIndex=particle->GetFirstDaughter();daughterIndex<=particle->GetLastDaughter();daughterIndex++){
796 TParticle *tmpDaughter = fMCStack->Particle(daughterIndex);
797 if(tmpDaughter->GetUniqueID() == 5){
798 if(tmpDaughter->GetPdgCode() == 11){
799 eNeg = tmpDaughter;
800 } else if(tmpDaughter->GetPdgCode() == -11){
801 ePos = tmpDaughter;
802 }
803 }
804 }
805 }
806
807 if(ePos == NULL || eNeg == NULL){ // means we do not have two daughters from pair production
808 return kFALSE;
809 }
e5b6e8a6 810
4a0aab28 811 if(ePos->Pt()<fSinglePtCut || eNeg->Pt()<fSinglePtCut){
812 return kFALSE; // no reconstruction below the Pt cut
813 }
e5b6e8a6 814
11c1e680 815 if( ePos->Eta() > (fEtaCut + fEtaShift) || ePos->Eta() < (-fEtaCut + fEtaShift) ||
816 eNeg->Eta() > (fEtaCut + fEtaShift) || eNeg->Eta() < (-fEtaCut + fEtaShift) )
4a0aab28 817 return kFALSE;
11c1e680 818
819 if(fEtaCutMin > -0.1){
820 if( (ePos->Eta() < (fEtaCutMin + fEtaShift) && ePos->Eta() > (-fEtaCutMin + fEtaShift)) ||
821 (eNeg->Eta() < (fEtaCutMin + fEtaShift) && eNeg->Eta() > (-fEtaCutMin + fEtaShift)) )
822 return kFALSE;
4a0aab28 823 }
e5b6e8a6 824
4a0aab28 825 if(ePos->R()>fMaxR){
826 return kFALSE; // cuts on distance from collision point
827 }
828
e5b6e8a6 829 if(abs(ePos->Vz()) > fMaxZ){
4a0aab28 830 return kFALSE; // outside material
831 }
e5b6e8a6 832 if(abs(eNeg->Vz()) > fMaxZ){
4a0aab28 833 return kFALSE; // outside material
834 }
835
e5b6e8a6 836 if( ePos->R() <= ((abs(ePos->Vz()) * fLineCutZRSlope) - fLineCutZValue)){
4a0aab28 837 return kFALSE; // line cut to exclude regions where we do not reconstruct
e5b6e8a6 838 } else if ( fEtaCutMin != -0.1 && ePos->R() >= ((abs(ePos->Vz()) * fLineCutZRSlopeMin) - fLineCutZValueMin)){
4a0aab28 839 return kFALSE;
840 }
e5b6e8a6 841
842 if( eNeg->R() <= ((abs(eNeg->Vz()) * fLineCutZRSlope) - fLineCutZValue)){
4a0aab28 843 return kFALSE; // line cut to exclude regions where we do not reconstruct
e5b6e8a6 844 } else if ( fEtaCutMin != -0.1 && eNeg->R() >= ((abs(eNeg->Vz()) * fLineCutZRSlopeMin) - fLineCutZValueMin)){
4a0aab28 845 return kFALSE;
846 }
e5b6e8a6 847
4a0aab28 848 return kTRUE;
849 //if(AcceptanceCut(particle,ePos,eNeg))return kTRUE;
850 }
851 return kFALSE;
3b77b2d1 852}
ae947965 853///________________________________________________________________________
854Bool_t AliConversionCuts::PhotonIsSelectedAODMC(AliAODMCParticle *particle,TClonesArray *aodmcArray,Bool_t checkForConvertedGamma){
855 // MonteCarlo Photon Selection
856
857 if(!aodmcArray)return kFALSE;
858
859 if (particle->GetPdgCode() == 22){
860 if( particle->Eta() > (fEtaCut + fEtaShift) || particle->Eta() < (-fEtaCut + fEtaShift) )
861 return kFALSE;
862 if(fEtaCutMin>-0.1){
863 if( particle->Eta() < (fEtaCutMin + fEtaShift) && particle->Eta() > (-fEtaCutMin + fEtaShift) )
864 return kFALSE;
865 }
866
867 if(particle->GetMother() > -1){
868 if((static_cast<AliAODMCParticle*>(aodmcArray->At(particle->GetMother())))->GetPdgCode() == 22){
869 return kFALSE; // no photon as mothers!
870 }
871 if(!(static_cast<AliAODMCParticle*>(aodmcArray->At(particle->GetMother()))->IsPrimary())){
872 return kFALSE; // the gamma has a mother, and it is not a primary particle
873 }
874 }
875
876 if(!checkForConvertedGamma) return kTRUE; // return in case of accepted gamma
877
878 // looking for conversion gammas (electron + positron from pairbuilding (= 5) )
879 AliAODMCParticle* ePos = NULL;
880 AliAODMCParticle* eNeg = NULL;
881
882 if(particle->GetNDaughters() >= 2){
883 for(Int_t daughterIndex=particle->GetDaughter(0);daughterIndex<=particle->GetDaughter(1);daughterIndex++){
884 AliAODMCParticle *tmpDaughter = static_cast<AliAODMCParticle*>(aodmcArray->At(daughterIndex));
885 if(!tmpDaughter) continue;
886 if(((tmpDaughter->GetMCProcessCode())) == 5){ // STILL A BUG IN ALIROOT >>8 HAS TPO BE REMOVED AFTER FIX
887 if(tmpDaughter->GetPdgCode() == 11){
888 eNeg = tmpDaughter;
889 } else if(tmpDaughter->GetPdgCode() == -11){
890 ePos = tmpDaughter;
891 }
892 }
893 }
894 }
895
896 if(ePos == NULL || eNeg == NULL){ // means we do not have two daughters from pair production
897 return kFALSE;
898 }
899
900 if(ePos->Pt()<fSinglePtCut || eNeg->Pt()<fSinglePtCut){
901 return kFALSE; // no reconstruction below the Pt cut
902 }
903
904 if( ePos->Eta() > (fEtaCut + fEtaShift) || ePos->Eta() < (-fEtaCut + fEtaShift) ||
905 eNeg->Eta() > (fEtaCut + fEtaShift) || eNeg->Eta() < (-fEtaCut + fEtaShift) )
906 return kFALSE;
907
908 if(fEtaCutMin > -0.1){
909 if( (ePos->Eta() < (fEtaCutMin + fEtaShift) && ePos->Eta() > (-fEtaCutMin + fEtaShift)) ||
910 (eNeg->Eta() < (fEtaCutMin + fEtaShift) && eNeg->Eta() > (-fEtaCutMin + fEtaShift)) )
911 return kFALSE;
912 }
913
914 Double_t rPos = sqrt( (ePos->Xv()*ePos->Xv()) + (ePos->Yv()*ePos->Yv()) );
915 Double_t rNeg = sqrt( (eNeg->Xv()*eNeg->Xv()) + (eNeg->Yv()*eNeg->Yv()) );
916
917 if(rPos>fMaxR){
918 return kFALSE; // cuts on distance from collision point
919 }
920 if(abs(ePos->Zv()) > fMaxZ){
921 return kFALSE; // outside material
922 }
923 if(abs(eNeg->Zv()) > fMaxZ){
924 return kFALSE; // outside material
925 }
3b77b2d1 926
ae947965 927 if( rPos <= ((abs(ePos->Zv()) * fLineCutZRSlope) - fLineCutZValue)){
928 return kFALSE; // line cut to exclude regions where we do not reconstruct
929 } else if ( fEtaCutMin != -0.1 && rPos >= ((abs(ePos->Zv()) * fLineCutZRSlopeMin) - fLineCutZValueMin)){
930 return kFALSE;
931 }
932
933 if( rNeg <= ((abs(eNeg->Zv()) * fLineCutZRSlope) - fLineCutZValue)){
934 return kFALSE; // line cut to exclude regions where we do not reconstruct
935 } else if ( fEtaCutMin != -0.1 && rNeg >= ((abs(eNeg->Zv()) * fLineCutZRSlopeMin) - fLineCutZValueMin)){
936 return kFALSE;
937 }
2bb2434e 938
ae947965 939 return kTRUE;
940 //if(AcceptanceCut(particle,ePos,eNeg))return kTRUE;
941 }
942 return kFALSE;
943}
3b77b2d1 944///________________________________________________________________________
945Bool_t AliConversionCuts::PhotonCuts(AliConversionPhotonBase *photon,AliVEvent *event)
946{ // Specific Photon Cuts
947
4a0aab28 948 Int_t cutIndex = 0;
949 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex);
950 cutIndex++;
951
952 // Fill Histos before Cuts
953 if(hInvMassbefore)hInvMassbefore->Fill(photon->GetMass());
e5b6e8a6 954
4a0aab28 955 if(hArmenterosbefore)hArmenterosbefore->Fill(photon->GetArmenterosAlpha(),photon->GetArmenterosQt());
956
957 // Gamma selection based on QT from Armenteros
958 if(fDoQtGammaSelection == kTRUE){
959 if(!ArmenterosQtCut(photon)){
960 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //1
961 return kFALSE;
962 }
963 }
964 cutIndex++; //2
965
966 // Chi Cut
967 if(photon->GetChi2perNDF() > fChi2CutConversion || photon->GetChi2perNDF() <=0){
968 {
969 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //2
970 return kFALSE;
971 }
972 }
973 cutIndex++;//3
974
975 // Reconstruction Acceptance Cuts
976 if(!AcceptanceCuts(photon)){
977 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //3
978 return kFALSE;
979 }
980
981 cutIndex++; //4
982 // Asymmetry Cut
983 if(fDoPhotonAsymmetryCut == kTRUE){
984 if(!AsymmetryCut(photon,event)){
985 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //4
986 return kFALSE;
987 }
988 }
989
990 //Check the pid probability
991 cutIndex++; //5
992 if(!PIDProbabilityCut(photon, event)) {
993 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //5
994 return kFALSE;
995 }
996
997 cutIndex++; //6
998 if(!CorrectedTPCClusterCut(photon, event)) {
999 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //6
1000 return kFALSE;
1001 }
1002
1003
1004 cutIndex++; //7
1005 if(!PsiPairCut(photon)) {
1006 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //7
1007 return kFALSE;
1008 }
1009
1010 cutIndex++; //8
1011 if(!CosinePAngleCut(photon, event)) {
1012 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //8
1013 return kFALSE;
1014 }
1015
4803eb1f 1016 AliAODConversionPhoton* photonAOD = dynamic_cast<AliAODConversionPhoton*>(photon);
1017 if (photonAOD){
1018 photonAOD->CalculateDistanceOfClossetApproachToPrimVtx(event->GetPrimaryVertex());
1019
1020 cutIndex++; //9
1021 if(photonAOD->GetDCArToPrimVtx() > fDCARPrimVtxCut) { //DCA R cut of photon to primary vertex
1022 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //9
1023 return kFALSE;
1024 }
1025
1026 cutIndex++; //10
1027 if(abs(photonAOD->GetDCAzToPrimVtx()) > fDCAZPrimVtxCut) { //DCA Z cut of photon to primary vertex
1028 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //10
1029 return kFALSE;
1030 }
1031 } else {
1032 cutIndex++; //9
1033 cutIndex++; //10
1034 }
1035 cutIndex++; //11
1036 if(hPhotonCuts)hPhotonCuts->Fill(cutIndex); //11
4a0aab28 1037
1038 // Histos after Cuts
1039 if(hInvMassafter)hInvMassafter->Fill(photon->GetMass());
1040 if(hArmenterosafter)hArmenterosafter->Fill(photon->GetArmenterosAlpha(),photon->GetArmenterosQt());
4803eb1f 1041
4a0aab28 1042 return kTRUE;
3b77b2d1 1043
1044}
1045
1046///________________________________________________________________________
1047Bool_t AliConversionCuts::CorrectedTPCClusterCut(AliConversionPhotonBase *photon, AliVEvent * event)
1048{ //Cut on corrected TPC Cluster Info
1049
4a0aab28 1050 AliVTrack * negTrack = GetTrack(event, photon->GetTrackLabelNegative());
1051 AliVTrack * posTrack = GetTrack(event, photon->GetTrackLabelPositive());
3b77b2d1 1052
4a0aab28 1053 if(!negTrack||!posTrack)return kFALSE;
3b77b2d1 1054
4a0aab28 1055 Double_t negclsToF=0;
3b77b2d1 1056
4a0aab28 1057 if (!fUseCorrectedTPCClsInfo ){
1058 if(negTrack->GetTPCNclsF()!=0){
1059 negclsToF = (Double_t)negTrack->GetNcls(1)/(Double_t)negTrack->GetTPCNclsF();}// Ncluster/Nfindablecluster
1060 }
1061 else {
1062 negclsToF = negTrack->GetTPCClusterInfo(2,0,GetFirstTPCRow(photon->GetConversionRadius()));
1063 }
3b77b2d1 1064
4a0aab28 1065 Double_t posclsToF = 0.;
1066 if (!fUseCorrectedTPCClsInfo ){
e5b6e8a6 1067 if(posTrack->GetTPCNclsF()!=0){
4a0aab28 1068 posclsToF = (Double_t)posTrack->GetNcls(1)/(Double_t)posTrack->GetTPCNclsF();
1069 }
1070 }else{
1071 posclsToF = posTrack->GetTPCClusterInfo(2,0,GetFirstTPCRow(photon->GetConversionRadius()));
1072 }
3b77b2d1 1073
e5b6e8a6 1074 if( negclsToF < fMinClsTPCToF || posclsToF < fMinClsTPCToF ){
4a0aab28 1075 return kFALSE;
1076 }
3b77b2d1 1077
4a0aab28 1078 return kTRUE;
3b77b2d1 1079}
1080
3b77b2d1 1081///________________________________________________________________________
d008165e 1082Bool_t AliConversionCuts::PhotonIsSelected(AliConversionPhotonBase *photon, AliVEvent * event)
1083{
4a0aab28 1084 //Selection of Reconstructed Photons
3b77b2d1 1085
4a0aab28 1086 FillPhotonCutIndex(kPhotonIn);
11c1e680 1087
0a2b2b4b 1088 if(event->IsA()==AliESDEvent::Class()) {
ccfa8c0d 1089 if(!SelectV0Finder( ( ((AliESDEvent*)event)->GetV0(photon->GetV0Index()))->GetOnFlyStatus() ) ){
0a2b2b4b 1090 FillPhotonCutIndex(kOnFly);
1091 return kFALSE;
1092 }
1093 }
1094 // else if(event->IsA()==AliAODEvent::Class()) {
1095 // if(!SelectV0Finder( ( ((AliAODEvent*)event)->GetV0(photon->GetV0Index())) ) ){
1096 // FillPhotonCutIndex(kOnFly);
1097 // return kFALSE;
1098 // }
1099 // }
3b77b2d1 1100
4a0aab28 1101 // Get Tracks
1102 AliVTrack * negTrack = GetTrack(event, photon->GetTrackLabelNegative());
1103 AliVTrack * posTrack = GetTrack(event, photon->GetTrackLabelPositive());
3b77b2d1 1104
4a0aab28 1105 if(!negTrack || !posTrack) {
1106 FillPhotonCutIndex(kNoTracks);
1107 return kFALSE;
1108 }
4803eb1f 1109 photon->DeterminePhotonQuality(negTrack,posTrack);
4a0aab28 1110 // Track Cuts
1111 if(!TracksAreSelected(negTrack, posTrack)){
1112 FillPhotonCutIndex(kTrackCuts);
1113 return kFALSE;
1114 }
4803eb1f 1115
1186afd2 1116 // dEdx Cuts
1117 if(!dEdxCuts(negTrack) || !dEdxCuts(posTrack)) {
1118 FillPhotonCutIndex(kdEdxCuts);
1119 return kFALSE;
1120 }
1121
4a0aab28 1122 // Photon Cuts
1123 if(!PhotonCuts(photon,event)){
1124 FillPhotonCutIndex(kPhotonCuts);
1125 return kFALSE;
1126 }
3b77b2d1 1127
4a0aab28 1128 // Photon passed cuts
1129 FillPhotonCutIndex(kPhotonOut);
1130 return kTRUE;
3b77b2d1 1131}
1132
3b77b2d1 1133///________________________________________________________________________
1134Bool_t AliConversionCuts::ArmenterosQtCut(AliConversionPhotonBase *photon)
1135{ // Armenteros Qt Cut
1136
4a0aab28 1137 if(fDoHighPtQtGammaSelection){
1138 if(photon->GetPhotonPt() < fPtBorderForQt){
1139 if(photon->GetArmenterosQt()>fQtMax){
1140 return kFALSE;
1141 }
1142 } else {
1143 if(photon->GetArmenterosQt()>fHighPtQtMax){
1144 return kFALSE;
1145 }
1146 }
1147 } else {
3b77b2d1 1148
4a0aab28 1149 if(photon->GetArmenterosQt()>fQtMax){
1150 return kFALSE;
1151 }
1152 }
1153 return kTRUE;
3b77b2d1 1154}
1155
1156
1157///________________________________________________________________________
1158Bool_t AliConversionCuts::AcceptanceCuts(AliConversionPhotonBase *photon) {
4a0aab28 1159 // Exclude certain areas for photon reconstruction
3b77b2d1 1160
4a0aab28 1161 Int_t cutIndex=0;
1162 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1163 cutIndex++;
3b77b2d1 1164
4a0aab28 1165 if(photon->GetConversionRadius()>fMaxR){ // cuts on distance from collision point
1166 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1167 return kFALSE;
1168 }
1169 cutIndex++;
3b77b2d1 1170
4a0aab28 1171 if(photon->GetConversionRadius()<fMinR){ // cuts on distance from collision point
1172 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1173 return kFALSE;
1174 }
1175 cutIndex++;
3b77b2d1 1176
e5b6e8a6 1177 if(photon->GetConversionRadius() <= ((abs(photon->GetConversionZ())*fLineCutZRSlope)-fLineCutZValue)){
3b77b2d1 1178 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1179 return kFALSE;
4a0aab28 1180 }
e5b6e8a6 1181 else if (fUseEtaMinCut && photon->GetConversionRadius() >= ((abs(photon->GetConversionZ())*fLineCutZRSlopeMin)-fLineCutZValueMin )){
4a0aab28 1182 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1183 return kFALSE;
1184 }
1185 cutIndex++;
3b77b2d1 1186
e5b6e8a6 1187 if(abs(photon->GetConversionZ()) > fMaxZ ){ // cuts out regions where we do not reconstruct
3b77b2d1 1188 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1189 return kFALSE;
4a0aab28 1190 }
1191 cutIndex++;
3b77b2d1 1192
1193
11c1e680 1194 if( photon->GetPhotonEta() > (fEtaCut + fEtaShift) || photon->GetPhotonEta() < (-fEtaCut + fEtaShift) ){
3b77b2d1 1195 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1196 return kFALSE;
4a0aab28 1197 }
11c1e680 1198 if(fEtaCutMin>-0.1){
1199 if( photon->GetPhotonEta() < (fEtaCutMin + fEtaShift) && photon->GetPhotonEta() > (-fEtaCutMin + fEtaShift) ){
1200 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1201 return kFALSE;
1202 }
1203 }
4a0aab28 1204 cutIndex++;
3b77b2d1 1205
4a0aab28 1206 if(photon->GetPhotonPt()<fPtCut){
3b77b2d1 1207 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
1208 return kFALSE;
4a0aab28 1209 }
1210 cutIndex++;
3b77b2d1 1211
4a0aab28 1212 if(hAcceptanceCuts)hAcceptanceCuts->Fill(cutIndex);
e5b6e8a6 1213
4a0aab28 1214 return kTRUE;
3b77b2d1 1215}
1216
1217
1218///________________________________________________________________________
1219Bool_t AliConversionCuts::SpecificTrackCuts(AliAODTrack * negTrack, AliAODTrack * posTrack,Int_t &cutIndex) {
4a0aab28 1220 // Track Cuts which require AOD/ESD specific implementation
3b77b2d1 1221
4a0aab28 1222 if( !negTrack->IsOn(AliESDtrack::kTPCrefit) || !posTrack->IsOn(AliESDtrack::kTPCrefit) ) {
3b77b2d1 1223 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1224 return kFALSE;
4a0aab28 1225 }
1226 cutIndex++;
3b77b2d1 1227
4a0aab28 1228 AliAODVertex * NegVtxType=negTrack->GetProdVertex();
1229 AliAODVertex * PosVtxType=posTrack->GetProdVertex();
4803eb1f 1230 if( (NegVtxType->GetType())==AliAODVertex::kKink || (PosVtxType->GetType())==AliAODVertex::kKink) {
3b77b2d1 1231 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1232 return kFALSE;
4a0aab28 1233 }
1234 return kTRUE;
3b77b2d1 1235
1236}
1237
1238
1239///________________________________________________________________________
1240Bool_t AliConversionCuts::SpecificTrackCuts(AliESDtrack * negTrack, AliESDtrack * posTrack,Int_t &cutIndex) {
4a0aab28 1241 // Track Cuts which require AOD/ESD specific implementation
3b77b2d1 1242
92efd725 1243 if( !negTrack->IsOn(AliESDtrack::kTPCrefit) || !posTrack->IsOn(AliESDtrack::kTPCrefit) ) {
3b77b2d1 1244 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1245 return kFALSE;
4a0aab28 1246 }
1247 cutIndex++;
3b77b2d1 1248
4a0aab28 1249 if(negTrack->GetKinkIndex(0) > 0 || posTrack->GetKinkIndex(0) > 0 ) {
3b77b2d1 1250 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1251 return kFALSE;
4a0aab28 1252 }
1253 return kTRUE;
3b77b2d1 1254}
1255
1256
1257
1258///________________________________________________________________________
1259Bool_t AliConversionCuts::TracksAreSelected(AliVTrack * negTrack, AliVTrack * posTrack) {
4a0aab28 1260 // Track Selection for Photon Reconstruction
3b77b2d1 1261
4a0aab28 1262 Int_t cutIndex=0;
1263 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1264 cutIndex++;
3b77b2d1 1265
4a0aab28 1266 // avoid like sign
1267 if(negTrack->Charge() == posTrack->Charge()) {
1268 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1269 return kFALSE;
1270 }
1271 cutIndex++;
3b77b2d1 1272
4a0aab28 1273 // Number of TPC Clusters
4803eb1f 1274
1275
4a0aab28 1276 if( negTrack->GetNcls(1) < fMinClsTPC || posTrack->GetNcls(1) < fMinClsTPC ) {
3b77b2d1 1277 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
4a0aab28 1278 return kFALSE;
1279 }
1280 cutIndex++;
ccfa8c0d 1281
4a0aab28 1282 // Acceptance
11c1e680 1283 if( posTrack->Eta() > (fEtaCut + fEtaShift) || posTrack->Eta() < (-fEtaCut + fEtaShift) ||
1284 negTrack->Eta() > (fEtaCut + fEtaShift) || negTrack->Eta() < (-fEtaCut + fEtaShift) ){
4a0aab28 1285 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
3b77b2d1 1286 return kFALSE;
4a0aab28 1287 }
11c1e680 1288 if(fEtaCutMin>-0.1){
1289 if( (posTrack->Eta() < (fEtaCutMin + fEtaShift) && posTrack->Eta() > (-fEtaCutMin + fEtaShift)) ||
1290 (negTrack->Eta() < (fEtaCutMin + fEtaShift) && negTrack->Eta() > (-fEtaCutMin + fEtaShift)) ){
1291 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1292 return kFALSE;
1293 }
1294 }
4a0aab28 1295 cutIndex++;
3b77b2d1 1296
4a0aab28 1297 // Single Pt Cut
1298 if( negTrack->Pt()< fSinglePtCut || posTrack->Pt()< fSinglePtCut){
3b77b2d1 1299 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1300 return kFALSE;
4a0aab28 1301 }
1302 cutIndex++;
3b77b2d1 1303
4a0aab28 1304 // AOD ESD specific cuts
1305 Bool_t passCuts = kTRUE;
3b77b2d1 1306
4a0aab28 1307 if(negTrack->IsA()==AliAODTrack::Class()) {
1308 passCuts = passCuts * SpecificTrackCuts(static_cast<AliAODTrack*>(negTrack), static_cast<AliAODTrack*>(posTrack),cutIndex);
e5b6e8a6 1309 } else {
4a0aab28 1310 passCuts = passCuts * SpecificTrackCuts(static_cast<AliESDtrack*>(negTrack), static_cast<AliESDtrack*>(posTrack),cutIndex);
e5b6e8a6 1311 }
3b77b2d1 1312
4a0aab28 1313 if(!passCuts){
3b77b2d1 1314 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
1315 return kFALSE;
4a0aab28 1316 }
1317 cutIndex++;
3b77b2d1 1318
4a0aab28 1319 if(hTrackCuts)hTrackCuts->Fill(cutIndex);
3b77b2d1 1320
4a0aab28 1321 return kTRUE;
e5b6e8a6 1322
3b77b2d1 1323}
1324
1325///________________________________________________________________________
1326Bool_t AliConversionCuts::dEdxCuts(AliVTrack *fCurrentTrack){
4a0aab28 1327 // Electron Identification Cuts for Photon reconstruction
3b77b2d1 1328
4a0aab28 1329 if(!fPIDResponse){InitPIDResponse();}// Try to reinitialize PID Response
1330 if(!fPIDResponse){AliError("No PID Response"); return kTRUE;}// if still missing fatal error
3b77b2d1 1331
4a0aab28 1332 Int_t cutIndex=0;
1333 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
e5b6e8a6 1334 if(hTPCdEdxSigbefore)hTPCdEdxSigbefore->Fill(fCurrentTrack->P(),fPIDResponse->NumberOfSigmasTPC(fCurrentTrack, AliPID::kElectron));
1335 if(hTPCdEdxbefore)hTPCdEdxbefore->Fill(fCurrentTrack->P(),fCurrentTrack->GetTPCsignal());
4a0aab28 1336 cutIndex++;
e5b6e8a6 1337
3b77b2d1 1338
4a0aab28 1339 if(fDodEdxSigmaCut == kTRUE){
3b77b2d1 1340 // TPC Electron Line
1341 if( fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)<fPIDnSigmaBelowElectronLine ||
4a0aab28 1342 fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)>fPIDnSigmaAboveElectronLine){
3b77b2d1 1343
4a0aab28 1344 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1345 return kFALSE;
3b77b2d1 1346 }
1347 cutIndex++;
1348
1349 // TPC Pion Line
4a0aab28 1350 if( fCurrentTrack->P()>fPIDMinPnSigmaAbovePionLine && fCurrentTrack->P()<fPIDMaxPnSigmaAbovePionLine ){
1351 if(fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)>fPIDnSigmaBelowElectronLine &&
1352 fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)<fPIDnSigmaAboveElectronLine&&
1353 fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kPion)<fPIDnSigmaAbovePionLine){
1354
1355 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1356 return kFALSE;
1357 }
1358 }
1359 cutIndex++;
e5b6e8a6 1360
4a0aab28 1361 // High Pt Pion rej
1362 if( fCurrentTrack->P()>fPIDMaxPnSigmaAbovePionLine ){
1363 if(fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)>fPIDnSigmaBelowElectronLine &&
e5b6e8a6 1364 fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kElectron)<fPIDnSigmaAboveElectronLine &&
4a0aab28 1365 fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kPion)<fPIDnSigmaAbovePionLineHighPt){
1366
1367 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1368 return kFALSE;
1369 }
1370 }
1371 cutIndex++;
1372 }
1373 else{cutIndex+=3;}
1374
1375 if(fDoKaonRejectionLowP == kTRUE){
1376 if(fCurrentTrack->P()<fPIDMinPKaonRejectionLowP ){
e5b6e8a6 1377 if( abs(fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kKaon))<fPIDnSigmaAtLowPAroundKaonLine){
4a0aab28 1378
1379 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1380 return kFALSE;
1381 }
1382 }
1383 }
1384 cutIndex++;
e5b6e8a6 1385
4a0aab28 1386 if(fDoProtonRejectionLowP == kTRUE){
1387 if( fCurrentTrack->P()<fPIDMinPProtonRejectionLowP ){
e5b6e8a6 1388 if( abs(fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kProton))<fPIDnSigmaAtLowPAroundProtonLine){
4a0aab28 1389
1390 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1391 return kFALSE;
1392 }
1393 }
1394 }
3b77b2d1 1395 cutIndex++;
e5b6e8a6 1396
4a0aab28 1397 if(fDoPionRejectionLowP == kTRUE){
1398 if( fCurrentTrack->P()<fPIDMinPPionRejectionLowP ){
e5b6e8a6 1399 if( abs(fPIDResponse->NumberOfSigmasTPC(fCurrentTrack,AliPID::kPion))<fPIDnSigmaAtLowPAroundPionLine){
4a0aab28 1400
1401 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1402 return kFALSE;
1403 }
1404 }
1405 }
1406 cutIndex++;
e5b6e8a6 1407
1408
1409 // cout<<"Start"<<endl;
1410 // AliPIDResponse::EDetPidStatus status=fPIDResponse->CheckPIDStatus(AliPIDResponse::kTOF,fCurrentTrack);
a280ac15 1411
e5b6e8a6 1412 // if( ( (status & AliVTrack::kTOFout) == AliVTrack::kTOFout ) && ( (status & AliVTrack::kTIME) == AliVTrack::kTIME ))
1413 // {cout<<"TOF DA"<<endl;}
a280ac15 1414 // if(status == AliPIDResponse::kDetPidOk){
e5b6e8a6 1415 // Float_t probMis = fPIDResponse->GetTOFMismatchProbability(fCurrentTrack);
1416 // cout<<"--> "<<probMis<<endl;
1417 // if(probMis > 0.01){
a280ac15 1418
e5b6e8a6 1419 // }
1420 // }
a280ac15 1421
4a0aab28 1422 if((fCurrentTrack->GetStatus() & AliESDtrack::kTOFpid) && !(fCurrentTrack->GetStatus() & AliESDtrack::kTOFmismatch)){
e5b6e8a6 1423 if(hTOFbefore){
1424 Double_t t0 = fPIDResponse->GetTOFResponse().GetStartTime(fCurrentTrack->P());
1425 Double_t times[5];
1426 fCurrentTrack->GetIntegratedTimes(times);
1427 Double_t TOFsignal = fCurrentTrack->GetTOFsignal();
1428 Double_t dT = TOFsignal - t0 - times[0];
1429 hTOFbefore->Fill(fCurrentTrack->P(),dT);
1430 }
1431 if(hTOFSigbefore) hTOFSigbefore->Fill(fCurrentTrack->P(),fPIDResponse->NumberOfSigmasTOF(fCurrentTrack, AliPID::kElectron));
4a0aab28 1432 if(fUseTOFpid){
1433 if(fPIDResponse->NumberOfSigmasTOF(fCurrentTrack, AliPID::kElectron)>fTofPIDnSigmaAboveElectronLine ||
1434 fPIDResponse->NumberOfSigmasTOF(fCurrentTrack, AliPID::kElectron)<fTofPIDnSigmaBelowElectronLine ){
1435 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1436 return kFALSE;
1437 }
1438 }
e5b6e8a6 1439 if(hTOFSigafter)hTOFSigafter->Fill(fCurrentTrack->P(),fPIDResponse->NumberOfSigmasTOF(fCurrentTrack, AliPID::kElectron));
4a0aab28 1440 }
1441 cutIndex++;
e5b6e8a6 1442
4a0aab28 1443 // Apply TRD PID
1444 if(fDoTRDPID){
1445 if(!fPIDResponse->IdentifiedAsElectronTRD(fCurrentTrack,fPIDTRDEfficiency)){
4a0aab28 1446 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
1447 return kFALSE;
1448 }
1449 }
1450 cutIndex++;
3b77b2d1 1451
4a0aab28 1452 if(hdEdxCuts)hdEdxCuts->Fill(cutIndex);
e5b6e8a6 1453 if(hTPCdEdxSigafter)hTPCdEdxSigafter->Fill(fCurrentTrack->P(),fPIDResponse->NumberOfSigmasTPC(fCurrentTrack, AliPID::kElectron));
1454 if(hTPCdEdxafter)hTPCdEdxafter->Fill(fCurrentTrack->P(),fCurrentTrack->GetTPCsignal());
4a0aab28 1455 return kTRUE;
3b77b2d1 1456}
1457
1458///________________________________________________________________________
1459Bool_t AliConversionCuts::AsymmetryCut(AliConversionPhotonBase * photon,AliVEvent *event) {
4a0aab28 1460 // Cut on Energy Assymetry
3b77b2d1 1461
4a0aab28 1462 for(Int_t ii=0;ii<2;ii++){
3b77b2d1 1463
4a0aab28 1464 AliVTrack *track=GetTrack(event,photon->GetTrackLabel(ii));
3b77b2d1 1465
4a0aab28 1466 if( track->P() > fMinPPhotonAsymmetryCut ){
1467 Double_t trackNegAsy=0;
1468 if (photon->GetPhotonP()!=0.){
1469 trackNegAsy= track->P()/photon->GetPhotonP();
1470 }
4803eb1f 1471
4a0aab28 1472 if( trackNegAsy<fMinPhotonAsymmetry ||trackNegAsy>(1.- fMinPhotonAsymmetry)){
1473 return kFALSE;
1474 }
1475 }
3b77b2d1 1476
4a0aab28 1477 }
1478 return kTRUE;
3b77b2d1 1479}
1480
1481///________________________________________________________________________
92efd725 1482AliVTrack *AliConversionCuts::GetTrack(AliVEvent * event, Int_t label){
4a0aab28 1483 //Returns pointer to the track with given ESD label
1484 //(Important for AOD implementation, since Track array in AOD data is different
1485 //from ESD array, but ESD tracklabels are stored in AOD Tracks)
1486
1487 AliESDEvent * esdEvent = dynamic_cast<AliESDEvent*>(event);
1488 if(esdEvent) {
1489 if(label > event->GetNumberOfTracks() ) return NULL;
1490 AliESDtrack * track = esdEvent->GetTrack(label);
1491 return track;
e5b6e8a6 1492
1493 } else {
ae947965 1494 AliVTrack * track = 0x0;
1186afd2 1495 if(((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask("V0ReaderV1"))->AreAODsRelabeled()){
1496 track = dynamic_cast<AliVTrack*>(event->GetTrack(label));
1497 return track;
1498 }
1499 else{
1500 for(Int_t ii=0; ii<event->GetNumberOfTracks(); ii++) {
1501 track = dynamic_cast<AliVTrack*>(event->GetTrack(ii));
1502 if(track){
1503 if(track->GetID() == label) {
ae947965 1504 return track;
1505 }
1506 }
1507 }
1508 }
4a0aab28 1509 }
e5b6e8a6 1510 //AliDebug(5,(Form("track not found %d %d",label,event->GetNumberOfTracks()));
1511 return NULL;
1512}
1513
1514///________________________________________________________________________
1515AliESDtrack *AliConversionCuts::GetESDTrack(AliESDEvent * event, Int_t label){
1516 //Returns pointer to the track with given ESD label
1517 //(Important for AOD implementation, since Track array in AOD data is different
1518 //from ESD array, but ESD tracklabels are stored in AOD Tracks)
1519
1520 if(event) {
1521 if(label > event->GetNumberOfTracks() ) return NULL;
1522 AliESDtrack * track = event->GetTrack(label);
1523 return track;
1524 }
4a0aab28 1525 //AliDebug(5,(Form("track not found %d %d",label,event->GetNumberOfTracks()));
1526 return NULL;
3b77b2d1 1527}
1528
1529
e5b6e8a6 1530
3b77b2d1 1531///________________________________________________________________________
1532Bool_t AliConversionCuts::PIDProbabilityCut(AliConversionPhotonBase *photon, AliVEvent * event){
4a0aab28 1533 // Cut on Electron Probability for Photon Reconstruction
3b77b2d1 1534
4a0aab28 1535 AliESDEvent * esdEvent = dynamic_cast<AliESDEvent*>(event);
3b77b2d1 1536
4a0aab28 1537 if(esdEvent){
e5b6e8a6 1538
4a0aab28 1539 Bool_t iResult=kFALSE;
e5b6e8a6 1540
4a0aab28 1541 Double_t *posProbArray = new Double_t[AliPID::kSPECIES];
1542 Double_t *negProbArray = new Double_t[AliPID::kSPECIES];
e5b6e8a6 1543
4a0aab28 1544 AliESDtrack* negTrack = esdEvent->GetTrack(photon->GetTrackLabelNegative());
1545 AliESDtrack* posTrack = esdEvent->GetTrack(photon->GetTrackLabelPositive());
e5b6e8a6 1546
4a0aab28 1547 if(negProbArray && posProbArray){
e5b6e8a6 1548
4a0aab28 1549 negTrack->GetTPCpid(negProbArray);
1550 posTrack->GetTPCpid(posProbArray);
3b77b2d1 1551
4a0aab28 1552 if(negProbArray[AliPID::kElectron]>=fPIDProbabilityCutNegativeParticle && posProbArray[AliPID::kElectron]>=fPIDProbabilityCutPositiveParticle){
1553 iResult=kTRUE;
1554 }
1555 }
e5b6e8a6 1556
4a0aab28 1557 delete [] posProbArray;
1558 delete [] negProbArray;
1559 return iResult;
3b77b2d1 1560
4a0aab28 1561 } else {
3b77b2d1 1562 ///Not possible for AODs
1563 return kTRUE;
4a0aab28 1564 }
3b77b2d1 1565
1566
1567
e5b6e8a6 1568
3b77b2d1 1569}
1570
1571
1572///________________________________________________________________________
1573Bool_t AliConversionCuts::AcceptanceCut(TParticle *particle, TParticle * ePos,TParticle* eNeg){
4a0aab28 1574 // MC Acceptance Cuts
1575 //(Certain areas were excluded for photon reconstruction)
3b77b2d1 1576
4a0aab28 1577 if(particle->R()>fMaxR){
1578 return kFALSE;}
3b77b2d1 1579
4a0aab28 1580 if(ePos->R()>fMaxR){
1581 return kFALSE;
1582 }
3b77b2d1 1583
4a0aab28 1584 if(ePos->R()<fMinR){
1585 return kFALSE;
1586 }
3b77b2d1 1587
e5b6e8a6 1588 if( ePos->R() <= ((abs(ePos->Vz())*fLineCutZRSlope)-fLineCutZValue)){
4a0aab28 1589 return kFALSE;
1590 }
e5b6e8a6 1591 else if (fUseEtaMinCut && ePos->R() >= ((abs(ePos->Vz())*fLineCutZRSlopeMin)-fLineCutZValueMin )){
3b77b2d1 1592 return kFALSE;
4a0aab28 1593 }
3b77b2d1 1594
e5b6e8a6 1595 if(abs(eNeg->Vz()) > fMaxZ ){ // cuts out regions where we do not reconstruct
4a0aab28 1596 return kFALSE;
1597 }
3b77b2d1 1598
4a0aab28 1599 if(eNeg->Vz()!=ePos->Vz()||eNeg->R()!=ePos->R()){
1600 return kFALSE;
1601 }
3b77b2d1 1602
e5b6e8a6 1603 if(abs(ePos->Vz()) > fMaxZ ){ // cuts out regions where we do not reconstruct
4a0aab28 1604 return kFALSE;
1605 }
3b77b2d1 1606
ccfa8c0d 1607
11c1e680 1608 if( particle->Eta() > (fEtaCut + fEtaShift) || particle->Eta() < (-fEtaCut + fEtaShift) ){
4a0aab28 1609 return kFALSE;
1610 }
11c1e680 1611 if( ePos->Eta() > (fEtaCut + fEtaShift) || ePos->Eta() < (-fEtaCut + fEtaShift) ){
4a0aab28 1612 return kFALSE;
1613 }
11c1e680 1614 if( eNeg->Eta() > (fEtaCut + fEtaShift) || eNeg->Eta() < (-fEtaCut + fEtaShift) ){
4a0aab28 1615 return kFALSE;
1616 }
11c1e680 1617 if(fEtaCutMin>-0.1){
1618 if( particle->Eta() < (fEtaCutMin + fEtaShift) && particle->Eta() > (-fEtaCutMin + fEtaShift) ){
1619 return kFALSE;
1620 }
1621 if( ePos->Eta() < (fEtaCutMin + fEtaShift) && ePos->Eta() > (-fEtaCutMin + fEtaShift) ){
1622 return kFALSE;
1623 }
1624 if( eNeg->Eta() < (fEtaCutMin + fEtaShift) && eNeg->Eta() > (-fEtaCutMin + fEtaShift) ){
1625 return kFALSE;
1626 }
1627 }
ccfa8c0d 1628
4a0aab28 1629 if( ePos->Pt()< fSinglePtCut || eNeg->Pt()< fSinglePtCut){
1630 return kFALSE;
1631 }
3b77b2d1 1632
4a0aab28 1633 if(particle->Pt()<fPtCut){
1634 return kFALSE;
1635 }
3b77b2d1 1636
4a0aab28 1637 return kTRUE;
3b77b2d1 1638}
1639///________________________________________________________________________
e5b6e8a6 1640Bool_t AliConversionCuts::UpdateCutString() {
4a0aab28 1641 ///Update the cut string (if it has been created yet)
3b77b2d1 1642
4a0aab28 1643 if(fCutString && fCutString->GetString().Length() == kNCuts) {
1644 fCutString->SetString(GetCutNumber());
1645 } else {
1646 return kFALSE;
1647 }
1648 return kTRUE;
3b77b2d1 1649}
1650
ae947965 1651void AliConversionCuts::LoadReweightingHistosMCFromFile() {
1652
1653 AliInfo("Entering loading of histograms for weighting");
1654 TFile *f = TFile::Open(fPathTrFReweighting.Data());
1655 if(!f){
1656 AliError(Form("file for weighting %s not found",fPathTrFReweighting.Data()));
1657 return;
1658 }
1659 if (fNameHistoReweightingPi0.CompareTo("") != 0 && fDoReweightHistoMCPi0 ){
1660 hReweightMCHistPi0 = (TH1D*)f->Get(fNameHistoReweightingPi0.Data());
1661 if (hReweightMCHistPi0) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingPi0.Data(),fPathTrFReweighting.Data() ));
1662 else AliWarning(Form("%s not found in %s",fPathTrFReweighting.Data(), fNameHistoReweightingPi0.Data() ));
1663 }
1664 if (fNameHistoReweightingEta.CompareTo("") != 0 && fDoReweightHistoMCEta){
1665 hReweightMCHistEta = (TH1D*)f->Get(fNameHistoReweightingEta.Data());
1666 if (hReweightMCHistEta) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingEta.Data(),fPathTrFReweighting.Data() ));
1667 else AliWarning(Form("%s not found in %s",fPathTrFReweighting.Data(), fNameHistoReweightingEta.Data() ));
1668
1669 }
1670 if (fNameHistoReweightingK0s.CompareTo("") != 0 && fDoReweightHistoMCK0s){
1671 hReweightMCHistK0s = (TH1D*)f->Get(fNameHistoReweightingK0s.Data());
1672 if (hReweightMCHistK0s) AliInfo(Form("%s has been loaded from %s", fNameHistoReweightingK0s.Data(),fPathTrFReweighting.Data() ));
1673 else AliWarning(Form("%s not found in %s",fPathTrFReweighting.Data(), fNameHistoReweightingK0s.Data() ));
1674
1675 }
1676
1677}
1678
1679
3b77b2d1 1680///________________________________________________________________________
1681Bool_t AliConversionCuts::InitializeCutsFromCutString(const TString analysisCutSelection ) {
1682 // Initialize Cuts from a given Cut string
ae947965 1683 if(fDoReweightHistoMCPi0 || fDoReweightHistoMCEta || fDoReweightHistoMCK0s) LoadReweightingHistosMCFromFile();
1684
4a0aab28 1685 AliInfo(Form("Set Photoncut Number: %s",analysisCutSelection.Data()));
1686 if(analysisCutSelection.Length()!=kNCuts) {
1687 AliError(Form("Cut selection has the wrong length! size is %d, number of cuts is %d", analysisCutSelection.Length(), kNCuts));
1688 return kFALSE;
1689 }
1690 if(!analysisCutSelection.IsDigit()){
1691 AliError("Cut selection contains characters");
1692 return kFALSE;
1693 }
e5b6e8a6 1694
4a0aab28 1695 const char *cutSelection = analysisCutSelection.Data();
1696#define ASSIGNARRAY(i) fCuts[i] = cutSelection[i] - '0'
1697 for(Int_t ii=0;ii<kNCuts;ii++){
3b77b2d1 1698 ASSIGNARRAY(ii);
4a0aab28 1699 }
3b77b2d1 1700
4a0aab28 1701 // Set Individual Cuts
1702 for(Int_t ii=0;ii<kNCuts;ii++){
3b77b2d1 1703 if(!SetCut(cutIds(ii),fCuts[ii]))return kFALSE;
4a0aab28 1704 }
3b77b2d1 1705
4a0aab28 1706 //PrintCuts();
92efd725 1707
4a0aab28 1708 return kTRUE;
3b77b2d1 1709}
1710///________________________________________________________________________
1711Bool_t AliConversionCuts::SetCut(cutIds cutID, const Int_t value) {
4a0aab28 1712 ///Set individual cut ID
3b77b2d1 1713
4a0aab28 1714 switch (cutID) {
e5b6e8a6 1715
4a0aab28 1716 case kv0FinderType:
1717 if( SetV0Finder(value)) {
1718 fCuts[kv0FinderType] = value;
e5b6e8a6 1719 UpdateCutString();
4a0aab28 1720 return kTRUE;
1721 } else return kFALSE;
1722
1723 case kededxSigmaCut:
1724 if( SetTPCdEdxCutElectronLine(value)) {
1725 fCuts[kededxSigmaCut] = value;
e5b6e8a6 1726 UpdateCutString();
4a0aab28 1727 return kTRUE;
1728 } else return kFALSE;
1729
1730 case kpidedxSigmaCut:
1731 if( SetTPCdEdxCutPionLine(value)) {
1732 fCuts[kpidedxSigmaCut] = value;
e5b6e8a6 1733 UpdateCutString();
4a0aab28 1734 return kTRUE;
1735 } else return kFALSE;
1736
1737 case kpiMomdedxSigmaCut:
1738 if( SetMinMomPiondEdxCut(value)) {
1739 fCuts[kpiMomdedxSigmaCut] = value;
e5b6e8a6 1740 UpdateCutString();
4a0aab28 1741 return kTRUE;
1742 } else return kFALSE;
1743
1744 case kchi2GammaCut:
1745 if( SetChi2GammaCut(value)) {
1746 fCuts[kchi2GammaCut] = value;
e5b6e8a6 1747 UpdateCutString();
4a0aab28 1748 return kTRUE;
1749 } else return kFALSE;
1750
1751 case ksinglePtCut:
1752 if( SetSinglePtCut(value)) {
1753 fCuts[ksinglePtCut] = value;
e5b6e8a6 1754 UpdateCutString();
4a0aab28 1755 return kTRUE;
1756 } else return kFALSE;
1757
1758 case kclsTPCCut:
1759 if( SetTPCClusterCut(value)) {
1760 fCuts[kclsTPCCut] = value;
e5b6e8a6 1761 UpdateCutString();
4a0aab28 1762 return kTRUE;
1763 } else return kFALSE;
1764
1765 case ketaCut:
1766 if( SetEtaCut(value)) {
1767 fCuts[ketaCut] = value;
e5b6e8a6 1768 UpdateCutString();
4a0aab28 1769 return kTRUE;
1770 } else return kFALSE;
1771
1772 case kLowPRejectionSigmaCut:
1773 if( SetLowPRejectionCuts(value)) {
1774 fCuts[kLowPRejectionSigmaCut] = value;
e5b6e8a6 1775 UpdateCutString();
4a0aab28 1776 return kTRUE;
1777 } else return kFALSE;
1778
1779 case kQtMaxCut:
1780 if( SetQtMaxCut(value)) {
1781 fCuts[kQtMaxCut] = value;
e5b6e8a6 1782 UpdateCutString();
4a0aab28 1783 return kTRUE;
1784 } else return kFALSE;
1785
1786 case kpiMaxMomdedxSigmaCut:
1787 if( SetMaxMomPiondEdxCut(value)) {
1788 fCuts[kpiMaxMomdedxSigmaCut] = value;
e5b6e8a6 1789 UpdateCutString();
4a0aab28 1790 return kTRUE;
1791 } else return kFALSE;
1792
1793 case kRCut:
1794 if( SetRCut(value)) {
1795 fCuts[kRCut] = value;
e5b6e8a6 1796 UpdateCutString();
4a0aab28 1797 return kTRUE;
1798 } else return kFALSE;
1799
1800 case kremovePileUp:
1801 if( SetRemovePileUp(value)) {
1802 fCuts[kremovePileUp] = value;
e5b6e8a6 1803 UpdateCutString();
4a0aab28 1804 return kTRUE;
1805 } else return kFALSE;
1806
1807 case kselectV0AND:
e5b6e8a6 1808 if( SetSelectSpecialTrigger(value)) {
4a0aab28 1809 fCuts[kselectV0AND] = value;
e5b6e8a6 1810 UpdateCutString();
4a0aab28 1811 return kTRUE;
1812 } else return kFALSE;
1813
1814 case kmultiplicityMethod:
1815 if( SetMultiplicityMethod(value)) {
1816 fCuts[kmultiplicityMethod] = value;
e5b6e8a6 1817 UpdateCutString();
4a0aab28 1818 return kTRUE;
1819 } else return kFALSE;
1820
1821 case kisHeavyIon:
1822 if( SetIsHeavyIon(value)) {
1823 fCuts[kisHeavyIon] = value;
e5b6e8a6 1824 UpdateCutString();
4a0aab28 1825 return kTRUE;
1826 } else return kFALSE;
1827
1828 case kCentralityMin:
1829 if( SetCentralityMin(value)) {
1830 fCuts[kCentralityMin] = value;
e5b6e8a6 1831 UpdateCutString();
4a0aab28 1832 return kTRUE;
1833 } else return kFALSE;
1834
1835 case kCentralityMax:
1836 if( SetCentralityMax(value)) {
1837 fCuts[kCentralityMax] = value;
e5b6e8a6 1838 UpdateCutString();
4a0aab28 1839 return kTRUE;
1840 } else return kFALSE;
1841
1842 case kTOFelectronPID:
1843 if( SetTOFElectronPIDCut(value)) {
1844 fCuts[kTOFelectronPID] = value;
e5b6e8a6 1845 UpdateCutString();
4a0aab28 1846 return kTRUE;
1847 } else return kFALSE;
1848
1849 case kdoPhotonAsymmetryCut:
1850 if( SetPhotonAsymmetryCut(value)) {
1851 fCuts[kdoPhotonAsymmetryCut] = value;
e5b6e8a6 1852 UpdateCutString();
4a0aab28 1853 return kTRUE;
1854 } else return kFALSE;
1855
1856 case kPsiPair:
1857 if( SetPsiPairCut(value)) {
1858 fCuts[kPsiPair] = value;
e5b6e8a6 1859 UpdateCutString();
4a0aab28 1860 return kTRUE;
1861 } else return kFALSE;
1862
1863 case kCosPAngle:
1864 if( SetCosPAngleCut(value)) {
1865 fCuts[kCosPAngle] = value;
e5b6e8a6 1866 UpdateCutString();
4a0aab28 1867 return kTRUE;
1868 } else return kFALSE;
1869
1870
1871 case kElecShare:
1872 if( SetSharedElectronCut(value)) {
1873 fCuts[kElecShare] = value;
e5b6e8a6 1874 UpdateCutString();
4a0aab28 1875 return kTRUE;
1876 } else return kFALSE;
1877
4a0aab28 1878 case kToCloseV0s:
1879 if( SetToCloseV0sCut(value)) {
1880 fCuts[kToCloseV0s] = value;
e5b6e8a6 1881 UpdateCutString();
4a0aab28 1882 return kTRUE;
1883 } else return kFALSE;
1884
1885 case kExtraSignals:
1886 if( SetRejectExtraSignalsCut(value)) {
1887 fCuts[kExtraSignals] = value;
e5b6e8a6 1888 UpdateCutString();
4a0aab28 1889 return kTRUE;
1890 } else return kFALSE;
1891
4803eb1f 1892 case kDcaRPrimVtx:
1893 if( SetDCARPhotonPrimVtxCut(value)) {
1894 fCuts[kDcaRPrimVtx] = value;
1895 UpdateCutString();
1896 return kTRUE;
1897 } else return kFALSE;
1898
1899 case kDcaZPrimVtx:
1900 if( SetDCAZPhotonPrimVtxCut(value)) {
1901 fCuts[kDcaZPrimVtx] = value;
1902 UpdateCutString();
1903 return kTRUE;
1904 } else return kFALSE;
1905
1906
4a0aab28 1907 case kNCuts:
92efd725 1908 AliError("Cut id out of range");
4a0aab28 1909 return kFALSE;
1910 }
d008165e 1911
4a0aab28 1912 AliError("Cut id %d not recognized");
1913 return kFALSE;
3b77b2d1 1914
e5b6e8a6 1915
3b77b2d1 1916}
1917///________________________________________________________________________
e5b6e8a6 1918void AliConversionCuts::PrintCuts() {
1919 // Print out current Cut Selection
1920 for(Int_t ic = 0; ic < kNCuts; ic++) {
1921 printf("%-30s : %d \n", fgkCutNames[ic], fCuts[ic]);
1922 }
1923}
1924///________________________________________________________________________
1925Bool_t AliConversionCuts::SetIsHeavyIon(Int_t isHeavyIon)
1926{ // Set Cut
1927 switch(isHeavyIon){
4a0aab28 1928 case 0:
e5b6e8a6 1929 fIsHeavyIon=0;
4a0aab28 1930 break;
1931 case 1:
e5b6e8a6 1932 fIsHeavyIon=1;
1933 fDetectorCentrality=0;
1934 break;
1935 case 2:
1936 fIsHeavyIon=1;
1937 fDetectorCentrality=1;
1938 break;
1939 case 3: //allows to select centrality 0-45% in steps of 5% for V0 Multiplicity
1940 fIsHeavyIon=1;
1941 fDetectorCentrality=0;
1942 fModCentralityClass=1;
1943 break;
1944 case 4: //allows to select centrality 45-90% in steps of 5% for V0 Multiplicity
1945 fIsHeavyIon=1;
1946 fDetectorCentrality=0;
1947 fModCentralityClass=2;
1948 break;
1949 case 5: //strict cut on v0 tracks for MC
1950 fIsHeavyIon=1;
1951 fDetectorCentrality=0;
1952 fModCentralityClass=3;
1953 break;
4803eb1f 1954 case 6: //allows to select centrality 0-45% in steps of 5% for track mult
e5b6e8a6 1955 //strict cut on v0 tracks for MC
1956 fIsHeavyIon=1;
1957 fDetectorCentrality=0;
1958 fModCentralityClass=4;
1959 break;
1960 case 7: //allows to select centrality 45-90% in steps of 5% for V0 Multiplicity
1961 //strict cut on v0 tracks for MC
1962 fIsHeavyIon=1;
1963 fDetectorCentrality=0;
1964 fModCentralityClass=5;
4a0aab28 1965 break;
0a2b2b4b 1966 case 8:
1967 fIsHeavyIon=2;
1968 fDetectorCentrality=0;
1969 break;
1970 case 9:
1971 fIsHeavyIon=2;
1972 fDetectorCentrality=1;
11c1e680 1973 break;
4a0aab28 1974 default:
e5b6e8a6 1975 AliError(Form("SetHeavyIon not defined %d",isHeavyIon));
4a0aab28 1976 return kFALSE;
1977 }
1978 return kTRUE;
3b77b2d1 1979}
0a2b2b4b 1980
e5b6e8a6 1981//___________________________________________________________________
1982Bool_t AliConversionCuts::SetCentralityMin(Int_t minCentrality)
1983{
1984 // Set Cut
1985 if(minCentrality<0||minCentrality>9){
1986 AliError(Form("minCentrality not defined %d",minCentrality));
1987 return kFALSE;
1988 }
3b77b2d1 1989
e5b6e8a6 1990 fCentralityMin=minCentrality;
1991 return kTRUE;
1992}
1993//___________________________________________________________________
1994Bool_t AliConversionCuts::SetCentralityMax(Int_t maxCentrality)
1995{
1996 // Set Cut
1997 if(maxCentrality<0||maxCentrality>9){
1998 AliError(Form("maxCentrality not defined %d",maxCentrality));
1999 return kFALSE;
2000 }
2001 fCentralityMax=maxCentrality;
2002 return kTRUE;
2003}
3b77b2d1 2004///________________________________________________________________________
e5b6e8a6 2005Int_t AliConversionCuts::SetSelectSpecialTrigger(Int_t selectSpecialTrigger)
3b77b2d1 2006{// Set Cut
e5b6e8a6 2007
2008 switch(selectSpecialTrigger){
4a0aab28 2009 case 0:
e5b6e8a6 2010 fSpecialTrigger=0; // dont care
3b77b2d1 2011 break;
4a0aab28 2012 case 1:
e5b6e8a6 2013 fSpecialTrigger=1; // V0AND
2014 break;
2015 case 2:
2016 fSpecialTrigger=2; // with SDD requested
2017 break;
2018 case 3:
2019 fSpecialTrigger=3; // V0AND plus with SDD requested
3b77b2d1 2020 break;
ae947965 2021 // allows to run MB & 6 other different trigger classes in parallel with the same photon cut
2022 case 4:
2023 fSpecialTrigger=4; // different trigger class as MB
2024 fTriggerSelectedManually = kTRUE;
2025 break;
2026 case 5:
2027 fSpecialTrigger=4; // different trigger class as MB
2028 fTriggerSelectedManually = kTRUE;
2029 break;
2030 case 6:
2031 fSpecialTrigger=4; // different trigger class as MB
2032 fTriggerSelectedManually = kTRUE;
2033 break;
2034 case 7:
2035 fSpecialTrigger=4; // different trigger class as MB
2036 fTriggerSelectedManually = kTRUE;
2037 break;
2038 case 8:
2039 fSpecialTrigger=4; // different trigger class as MB
2040 fTriggerSelectedManually = kTRUE;
2041 break;
2042 case 9:
2043 fSpecialTrigger=4; // different trigger class as MB
2044 fTriggerSelectedManually = kTRUE;
2045 break;
4a0aab28 2046 default:
e5b6e8a6 2047 AliError("Warning: Special Trigger Not known");
3b77b2d1 2048 return kFALSE;
4a0aab28 2049 }
2050 return kTRUE;
3b77b2d1 2051}
3b77b2d1 2052///________________________________________________________________________
92efd725 2053Bool_t AliConversionCuts::SetMultiplicityMethod(Int_t multiplicityMethod)
2054{
4a0aab28 2055 // Set Cut
2056 fMultiplicityMethod=multiplicityMethod;
92efd725 2057
4a0aab28 2058 // 0 Photon Multiplicity
2059 // 1 TPC Track multiplicity
2060 // 2 V0 Mult
2061 // 3 SPD Mult
92efd725 2062
4a0aab28 2063 return kTRUE;
3b77b2d1 2064}
3b77b2d1 2065///________________________________________________________________________
e5b6e8a6 2066Bool_t AliConversionCuts::SetRemovePileUp(Int_t removePileUp)
2067{// Set Cut
2068 switch(removePileUp){
2069 case 0:
2070 fRemovePileUp=kFALSE;
2071 break;
2072 case 1:
2073 fRemovePileUp=kTRUE;
2074 break;
2075 default:
2076 AliError("RemovePileUpCut not defined");
2077 return kFALSE;
4a0aab28 2078 }
e5b6e8a6 2079 return kTRUE;
3b77b2d1 2080}
3b77b2d1 2081///________________________________________________________________________
e5b6e8a6 2082Bool_t AliConversionCuts::SetRejectExtraSignalsCut(Int_t extraSignal) {
2083
2084 switch(extraSignal){
4a0aab28 2085 case 0:
e5b6e8a6 2086 fRejectExtraSignals = 0;
2087 break; // No Rejection
4a0aab28 2088 case 1:
e5b6e8a6 2089 fRejectExtraSignals = 1;
2090 break; // MinBias Header
4a0aab28 2091 case 2:
e5b6e8a6 2092 fRejectExtraSignals = 2;
2093 break; // User String Array
4a0aab28 2094 case 3:
e5b6e8a6 2095 fRejectExtraSignals = 3;
2096 break; // Rejection for Gamma Correction only
2097 default:
2098 AliError(Form("Extra Signal Rejection not defined %d",extraSignal));
2099 return kFALSE;
2100 }
2101 return kTRUE;
2102}
2103///________________________________________________________________________
2104Bool_t AliConversionCuts::SetV0Finder(Int_t v0FinderType)
2105{ // Set Cut
2106 switch (v0FinderType){
2107 case 0: // on fly V0 finder
ccfa8c0d 2108 cout << "have chosen onfly V0" << endl;
e5b6e8a6 2109 fUseOnFlyV0Finder=kTRUE;
4a0aab28 2110 break;
e5b6e8a6 2111 case 1: // offline V0 finder
ccfa8c0d 2112 cout << "have chosen offline V0" << endl;
e5b6e8a6 2113 fUseOnFlyV0Finder=kFALSE;
4a0aab28 2114 break;
4a0aab28 2115 default:
e5b6e8a6 2116 AliError(Form(" v0FinderType not defined %d",v0FinderType));
4a0aab28 2117 return kFALSE;
2118 }
2119 return kTRUE;
3b77b2d1 2120}
3b77b2d1 2121///________________________________________________________________________
e5b6e8a6 2122Bool_t AliConversionCuts::SetEtaCut(Int_t etaCut)
3b77b2d1 2123{ // Set Cut
e5b6e8a6 2124
2125 //Set Standard LineCutZValues
2126 fLineCutZValueMin = -2;
2127 fLineCutZValue = 7.;
2128
2129 switch(etaCut){
2130 case 0: // 0.9
2131 fEtaCut = 0.9;
2132 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2133 fEtaCutMin = -0.1;
2134 fLineCutZRSlopeMin = 0.;
4a0aab28 2135 break;
ae947965 2136 case 1: // 1.2 // changed from 1.2 to 0.6 on 2013.06.10
2137 fEtaCut = 0.6;
e5b6e8a6 2138 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2139 fEtaCutMin = -0.1;
2140 fLineCutZRSlopeMin = 0.;
4a0aab28 2141 break;
e5b6e8a6 2142 case 2: // 1.4
2143 fEtaCut = 1.4;
2144 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2145 fEtaCutMin = -0.1;
2146 fLineCutZRSlopeMin = 0.;
4a0aab28 2147 break;
e5b6e8a6 2148 case 3: // 0.8
2149 fEtaCut = 0.8;
2150 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2151 fEtaCutMin = -0.1;
2152 fLineCutZRSlopeMin = 0.;
4a0aab28 2153 break;
a280ac15 2154 case 4: // 0.75
2155 fEtaCut = 0.75;
e5b6e8a6 2156 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2157 fEtaCutMin = -0.1;
2158 fLineCutZRSlopeMin = 0.;
4a0aab28 2159 break;
11c1e680 2160 case 5: // 0.5
2161 fEtaCut = 0.5;
e5b6e8a6 2162 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
11c1e680 2163 fEtaCutMin = -0.1;
2164 fLineCutZRSlopeMin = 0.;
4a0aab28 2165 break;
e5b6e8a6 2166 case 6: // 5.
2167 fEtaCut = 5.;
2168 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2169 fEtaCutMin = -0.1;
2170 fLineCutZRSlopeMin = 0.;
4a0aab28 2171 break;
ccfa8c0d 2172 case 7:
11c1e680 2173 fEtaCut = 0.3;
e5b6e8a6 2174 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
11c1e680 2175 fEtaCutMin = -0.1;
2176 fLineCutZRSlopeMin = 0.;
e5b6e8a6 2177 break;
11c1e680 2178 // case 8: // 0.1 - 0.8
2179 // fEtaCut = 0.9;
2180 // fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2181 // fEtaCutMin = 0.1;
2182 // fLineCutZRSlopeMin = tan(2*atan(exp(-fEtaCutMin)));
2183 // break;
2184 case 8: // 0.4
2185 fEtaCut = 0.4;
e5b6e8a6 2186 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
11c1e680 2187 fEtaCutMin = -0.1;
2188 fLineCutZRSlopeMin = 0.;
e5b6e8a6 2189 break;
2190 case 9: // 10
2191 fEtaCut = 10;
2192 fLineCutZRSlope = tan(2*atan(exp(-fEtaCut)));
2193 fEtaCutMin = -0.1;
2194 fLineCutZRSlopeMin = 0.;
4a0aab28 2195 break;
2196 default:
e5b6e8a6 2197 AliError(Form(" EtaCut not defined %d",etaCut));
4a0aab28 2198 return kFALSE;
4a0aab28 2199 }
2200 return kTRUE;
3b77b2d1 2201}
3b77b2d1 2202///________________________________________________________________________
e5b6e8a6 2203Bool_t AliConversionCuts::SetRCut(Int_t RCut){
2204 // Set Cut
2205 switch(RCut){
2206 case 0:
2207 fMinR=0;
2208 fMaxR = 180.;
4a0aab28 2209 break;
e5b6e8a6 2210 case 1:
2211 fMinR=2.8;
2212 fMaxR = 180.;
4a0aab28 2213 break;
e5b6e8a6 2214 case 2:
2215 fMinR=5.;
2216 fMaxR = 180.;
4a0aab28 2217 break;
e5b6e8a6 2218 case 3:
2219 fMaxR = 70.;
2220 fMinR = 10.;
4a0aab28 2221 break;
e5b6e8a6 2222 case 4:
2223 fMaxR = 70.;
2224 fMinR = 5.;
4a0aab28 2225 break;
e5b6e8a6 2226 // High purity cuts for PbPb (remove first layers of material)
2227 case 5:
2228 fMaxR = 180.;
2229 fMinR = 10.;
4a0aab28 2230 break;
e5b6e8a6 2231 case 6:
2232 fMaxR = 180.;
2233 fMinR = 20.;
4a0aab28 2234 break;
e5b6e8a6 2235 case 7:
2236 fMaxR = 180.;
2237 fMinR = 26.;
4a0aab28 2238 break;
e5b6e8a6 2239 case 8:
2240 fMaxR = 180.;
2241 fMinR = 35.;
4a0aab28 2242 break;
2243 case 9:
e5b6e8a6 2244 fMaxR = 35.;
2245 fMinR = 5.;
4a0aab28 2246 break;
e5b6e8a6 2247
4a0aab28 2248 default:
e5b6e8a6 2249 AliError("RCut not defined");
4a0aab28 2250 return kFALSE;
2251 }
2252 return kTRUE;
92efd725 2253}
92efd725 2254///________________________________________________________________________
e5b6e8a6 2255Bool_t AliConversionCuts::SetSinglePtCut(Int_t singlePtCut)
92efd725 2256{ // Set Cut
e5b6e8a6 2257 switch(singlePtCut){
2258 case 0: // 0.050 GeV
2259 fSinglePtCut = 0.050;
4a0aab28 2260 break;
e5b6e8a6 2261 case 1: // 0.100 GeV
2262 fSinglePtCut = 0.100;
4a0aab28 2263 break;
e5b6e8a6 2264 case 2: // 0.150 GeV
2265 fSinglePtCut = 0.150;
4a0aab28 2266 break;
e5b6e8a6 2267 case 3: // 0.200 GeV
2268 fSinglePtCut = 0.200;
4a0aab28 2269 break;
e5b6e8a6 2270 case 4: // 0.075 GeV
2271 fSinglePtCut = 0.075;
4a0aab28 2272 break;
e5b6e8a6 2273 case 5: // 0.125 GeV
2274 fSinglePtCut = 0.125;
4a0aab28 2275 break;
e5b6e8a6 2276 case 6: // 0.04 GeV
2277 fSinglePtCut = 0.040;
4a0aab28 2278 break;
e5b6e8a6 2279 case 7: // 0.0 GeV
2280 fSinglePtCut = 0.0;
4a0aab28 2281 break;
2282 default:
e5b6e8a6 2283 AliError(Form("singlePtCut not defined %d",singlePtCut));
4a0aab28 2284 return kFALSE;
2285 }
2286 return kTRUE;
3b77b2d1 2287}
3b77b2d1 2288///________________________________________________________________________
e5b6e8a6 2289Bool_t AliConversionCuts::SetTPCClusterCut(Int_t clsTPCCut)
3b77b2d1 2290{ // Set Cut
e5b6e8a6 2291 switch(clsTPCCut){
2292 case 0: // 0
2293 fMinClsTPC= 0.;
4a0aab28 2294 break;
ccfa8c0d 2295 case 1: // 60
2296 fMinClsTPC= 60.;
4a0aab28 2297 break;
e5b6e8a6 2298 case 2: // 80
2299 fMinClsTPC= 80.;
4a0aab28 2300 break;
e5b6e8a6 2301 case 3: // 100
2302 fMinClsTPC= 100.;
4a0aab28 2303 break;
ccfa8c0d 2304 case 4: // 95% of findable clusters
2305 fMinClsTPCToF= 0.95;
2306 fUseCorrectedTPCClsInfo=1;
4a0aab28 2307 break;
e5b6e8a6 2308 case 5: // 0% of findable clusters
2309 fMinClsTPCToF= 0.0;
2310 fUseCorrectedTPCClsInfo=1;
4a0aab28 2311 break;
e5b6e8a6 2312 case 6: // 70% of findable clusters
2313 fMinClsTPCToF= 0.7;
2314 fUseCorrectedTPCClsInfo=1;
4a0aab28 2315 break;
e5b6e8a6 2316 case 7: // 0% of findable clusters
2317 fMinClsTPCToF= 0.35;
2318 fUseCorrectedTPCClsInfo=0;
4a0aab28 2319 break;
2320 case 8:
e5b6e8a6 2321 fMinClsTPCToF= 0.35;
2322 fUseCorrectedTPCClsInfo=1;
4a0aab28 2323 break;
2324 case 9:
e5b6e8a6 2325 fMinClsTPCToF= 0.6;
2326 fUseCorrectedTPCClsInfo=1;
4a0aab28 2327 break;
2328 default:
e5b6e8a6 2329 AliError(Form("Warning: clsTPCCut not defined %d",clsTPCCut));
4a0aab28 2330 return kFALSE;
2331 }
2332 return kTRUE;
3b77b2d1 2333}
3b77b2d1 2334///________________________________________________________________________
e5b6e8a6 2335Bool_t AliConversionCuts::SetTPCdEdxCutElectronLine(Int_t ededxSigmaCut)
3b77b2d1 2336{ // Set Cut
e5b6e8a6 2337 switch(ededxSigmaCut){
2338 case 0: // -10,10
2339 fPIDnSigmaBelowElectronLine=-10;
2340 fPIDnSigmaAboveElectronLine=10;
4a0aab28 2341 break;
e5b6e8a6 2342 case 1: // -5,5
2343 fPIDnSigmaBelowElectronLine=-5;
2344 fPIDnSigmaAboveElectronLine=5;
4a0aab28 2345 break;
e5b6e8a6 2346 case 2: // -3,5
2347 fPIDnSigmaBelowElectronLine=-3;
2348 fPIDnSigmaAboveElectronLine=5;
4a0aab28 2349 break;
e5b6e8a6 2350 case 3: // -4,5
2351 fPIDnSigmaBelowElectronLine=-4;
2352 fPIDnSigmaAboveElectronLine=5;
4a0aab28 2353 break;
e5b6e8a6 2354 case 4: // -6,7
2355 fPIDnSigmaBelowElectronLine=-6;
2356 fPIDnSigmaAboveElectronLine=7;
4a0aab28 2357 break;
e5b6e8a6 2358 case 5: // -4,4
2359 fPIDnSigmaBelowElectronLine=-4;
2360 fPIDnSigmaAboveElectronLine=4;
4a0aab28 2361 break;
e5b6e8a6 2362 case 6: // -2.5,4
2363 fPIDnSigmaBelowElectronLine=-2.5;
2364 fPIDnSigmaAboveElectronLine=4;
4a0aab28 2365 break;
e5b6e8a6 2366 case 7: // -2,3.5
2367 fPIDnSigmaBelowElectronLine=-2;
2368 fPIDnSigmaAboveElectronLine=3.5;
4a0aab28 2369 break;
2370 default:
e5b6e8a6 2371 AliError("TPCdEdxCutElectronLine not defined");
4a0aab28 2372 return kFALSE;
e5b6e8a6 2373
4a0aab28 2374 }
2375 return kTRUE;
3b77b2d1 2376}
2377///________________________________________________________________________
e5b6e8a6 2378Bool_t AliConversionCuts::SetTPCdEdxCutPionLine(Int_t pidedxSigmaCut)
3b77b2d1 2379{ // Set Cut
e5b6e8a6 2380
2381 switch(pidedxSigmaCut){
2382 case 0: // -10
2383 fPIDnSigmaAbovePionLine=-10;
2384 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2385 break;
e5b6e8a6 2386 case 1: // 0
2387 fPIDnSigmaAbovePionLine=0;
2388 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2389 break;
e5b6e8a6 2390 case 2: // 1
2391 fPIDnSigmaAbovePionLine=1;
2392 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2393 break;
e5b6e8a6 2394 case 3: // 1
2395 fPIDnSigmaAbovePionLine=2.5;
2396 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2397 break;
e5b6e8a6 2398 case 4: // 1
2399 fPIDnSigmaAbovePionLine=0.5;
2400 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2401 break;
e5b6e8a6 2402 case 5: // 1
2403 fPIDnSigmaAbovePionLine=2.;
2404 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2405 break;
e5b6e8a6 2406 case 6: // 1
2407 fPIDnSigmaAbovePionLine=2.;
2408 fPIDnSigmaAbovePionLineHighPt=0.5;
4a0aab28 2409 break;
e5b6e8a6 2410 case 7: // 1
2411 fPIDnSigmaAbovePionLine=3.5;
2412 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2413 break;
e5b6e8a6 2414 case 8: // 1
2415 fPIDnSigmaAbovePionLine=2.;
2416 fPIDnSigmaAbovePionLineHighPt=1.;
4a0aab28 2417 break;
2418 case 9:
e5b6e8a6 2419 fPIDnSigmaAbovePionLine=3.0; // We need a bit less tight cut on dE/dx
2420 fPIDnSigmaAbovePionLineHighPt=-10;
4a0aab28 2421 break;
2422 default:
e5b6e8a6 2423 AliError(Form("Warning: pidedxSigmaCut not defined %d",pidedxSigmaCut));
4a0aab28 2424 return kFALSE;
2425 }
2426 return kTRUE;
3b77b2d1 2427}
3b77b2d1 2428///________________________________________________________________________
e5b6e8a6 2429Bool_t AliConversionCuts::SetMinMomPiondEdxCut(Int_t piMomdedxSigmaCut)
3b77b2d1 2430{ // Set Cut
e5b6e8a6 2431 switch(piMomdedxSigmaCut){
2432 case 0: // 0.5 GeV
2433 fPIDMinPnSigmaAbovePionLine=0.5;
4a0aab28 2434 break;
e5b6e8a6 2435 case 1: // 1. GeV
2436 fPIDMinPnSigmaAbovePionLine=1.;
4a0aab28 2437 break;
e5b6e8a6 2438 case 2: // 1.5 GeV
2439 fPIDMinPnSigmaAbovePionLine=1.5;
4a0aab28 2440 break;
e5b6e8a6 2441 case 3: // 20.0 GeV
2442 fPIDMinPnSigmaAbovePionLine=20.;
4a0aab28 2443 break;
e5b6e8a6 2444 case 4: // 50.0 GeV
2445 fPIDMinPnSigmaAbovePionLine=50.;
4a0aab28 2446 break;
e5b6e8a6 2447 case 5: // 0.3 GeV
2448 fPIDMinPnSigmaAbovePionLine=0.3;
4a0aab28 2449 break;
e5b6e8a6 2450 case 6: // 0.25 GeV
2451 fPIDMinPnSigmaAbovePionLine=0.25;
4a0aab28 2452 break;
e5b6e8a6 2453 case 7: // 0.4 GeV
2454 fPIDMinPnSigmaAbovePionLine=0.4;
4a0aab28 2455 break;
e5b6e8a6 2456 case 8: // 0.2 GeV
2457 fPIDMinPnSigmaAbovePionLine=0.2;
4a0aab28 2458 break;
2459 default:
e5b6e8a6 2460 AliError(Form("piMomdedxSigmaCut not defined %d",piMomdedxSigmaCut));
4a0aab28 2461 return kFALSE;
2462 }
2463 return kTRUE;
3b77b2d1 2464}
3b77b2d1 2465///________________________________________________________________________
2466Bool_t AliConversionCuts::SetMaxMomPiondEdxCut(Int_t piMaxMomdedxSigmaCut)
2467{ // Set Cut
4a0aab28 2468 switch(piMaxMomdedxSigmaCut){
2469 case 0: // 100. GeV
2470 fPIDMaxPnSigmaAbovePionLine=100.;
2471 break;
2472 case 1: // 5. GeV
2473 fPIDMaxPnSigmaAbovePionLine=5.;
2474 break;
2475 case 2: // 4. GeV
2476 fPIDMaxPnSigmaAbovePionLine=4.;
2477 break;
2478 case 3: // 3.5 GeV
2479 fPIDMaxPnSigmaAbovePionLine=3.5;
2480 break;
2481 case 4: // 3. GeV
2482 fPIDMaxPnSigmaAbovePionLine=3.;
2483 break;
e5b6e8a6 2484 case 5: // 7. GeV
2485 fPIDMaxPnSigmaAbovePionLine=7.;
4a0aab28 2486 break;
2487 default:
e5b6e8a6 2488 AliError(Form("piMaxMomdedxSigmaCut not defined %d",piMaxMomdedxSigmaCut));
4a0aab28 2489 return kFALSE;
92efd725 2490 }
2491 return kTRUE;
3b77b2d1 2492}
3b77b2d1 2493///________________________________________________________________________
2494Bool_t AliConversionCuts::SetLowPRejectionCuts(Int_t LowPRejectionSigmaCut)
2495{ // Set Cut
4a0aab28 2496 switch(LowPRejectionSigmaCut){
2497 case 0: //
2498 fPIDnSigmaAtLowPAroundKaonLine=0;
2499 fPIDnSigmaAtLowPAroundProtonLine=0;
2500 fPIDnSigmaAtLowPAroundPionLine=0;
e5b6e8a6 2501 fDoKaonRejectionLowP = kFALSE;
2502 fDoProtonRejectionLowP = kFALSE;
2503 fDoPionRejectionLowP = kFALSE;
2504 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2505 break;
2506 case 1: //
2507 fPIDnSigmaAtLowPAroundKaonLine=0.5;
2508 fPIDnSigmaAtLowPAroundProtonLine=0.5;
2509 fPIDnSigmaAtLowPAroundPionLine=0.5;
e5b6e8a6 2510 fDoKaonRejectionLowP = kTRUE;
2511 fDoProtonRejectionLowP = kTRUE;
2512 fDoPionRejectionLowP = kTRUE;
2513 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2514 break;
2515 case 2: //
2516 fPIDnSigmaAtLowPAroundKaonLine=1;
2517 fPIDnSigmaAtLowPAroundProtonLine=1;
2518 fPIDnSigmaAtLowPAroundPionLine=1;
e5b6e8a6 2519 fDoKaonRejectionLowP = kTRUE;
2520 fDoProtonRejectionLowP = kTRUE;
2521 fDoPionRejectionLowP = kTRUE;
2522 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2523 break;
2524 case 3: //
2525 fPIDnSigmaAtLowPAroundKaonLine=2.;
2526 fPIDnSigmaAtLowPAroundProtonLine=2.;
2527 fPIDnSigmaAtLowPAroundPionLine=2.;
e5b6e8a6 2528 fDoKaonRejectionLowP = kTRUE;
2529 fDoProtonRejectionLowP = kTRUE;
2530 fDoPionRejectionLowP = kTRUE;
2531 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2532 break;
2533 case 4: //
2534 fPIDnSigmaAtLowPAroundKaonLine=0.;
2535 fPIDnSigmaAtLowPAroundProtonLine=0.;
2536 fPIDnSigmaAtLowPAroundPionLine=1;
e5b6e8a6 2537 fDoKaonRejectionLowP = kFALSE;
2538 fDoProtonRejectionLowP = kFALSE;
2539 fDoPionRejectionLowP = kTRUE;
2540 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2541 break;
2542 case 5: //
2543 fPIDnSigmaAtLowPAroundKaonLine=0.;
2544 fPIDnSigmaAtLowPAroundProtonLine=0.;
2545 fPIDnSigmaAtLowPAroundPionLine=1.5;
e5b6e8a6 2546 fDoKaonRejectionLowP = kFALSE;
2547 fDoProtonRejectionLowP = kFALSE;
2548 fDoPionRejectionLowP = kTRUE;
2549 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2550 break;
2551 case 6: //
2552 fPIDnSigmaAtLowPAroundKaonLine=0.;
2553 fPIDnSigmaAtLowPAroundProtonLine=0.;
2554 fPIDnSigmaAtLowPAroundPionLine=2.;
e5b6e8a6 2555 fDoKaonRejectionLowP = kFALSE;
2556 fDoProtonRejectionLowP = kFALSE;
2557 fDoPionRejectionLowP = kTRUE;
2558 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
2559 break;
2560 case 7: //
2561 fPIDnSigmaAtLowPAroundKaonLine=0.;
2562 fPIDnSigmaAtLowPAroundProtonLine=0.;
2563 fPIDnSigmaAtLowPAroundPionLine=0.5;
2564 fDoKaonRejectionLowP = kFALSE;
2565 fDoProtonRejectionLowP = kFALSE;
2566 fDoPionRejectionLowP = kTRUE;
2567 fPIDMinPPionRejectionLowP = fPIDMinPnSigmaAbovePionLine;
4a0aab28 2568 break;
2569 default:
2570 AliError(Form("LowPRejectionSigmaCut not defined %d",LowPRejectionSigmaCut));
2571 return kFALSE;
2572 }
2573 return kTRUE;
3b77b2d1 2574}
3b77b2d1 2575///________________________________________________________________________
2576Bool_t AliConversionCuts::SetTOFElectronPIDCut(Int_t TOFelectronPID){
4a0aab28 2577 // Set Cut
e5b6e8a6 2578 switch(TOFelectronPID){
4a0aab28 2579 case 0: // no cut
2580 fUseTOFpid = kFALSE;
2581 fTofPIDnSigmaBelowElectronLine=-100;
2582 fTofPIDnSigmaAboveElectronLine=100;
2583 break;
2584 case 1: // -7,7
2585 fUseTOFpid = kTRUE;
2586 fTofPIDnSigmaBelowElectronLine=-7;
2587 fTofPIDnSigmaAboveElectronLine=7;
2588 break;
2589 case 2: // -5,5
2590 fUseTOFpid = kTRUE;
2591 fTofPIDnSigmaBelowElectronLine=-5;
2592 fTofPIDnSigmaAboveElectronLine=5;
2593 break;
2594 case 3: // -3,5
2595 fUseTOFpid = kTRUE;
2596 fTofPIDnSigmaBelowElectronLine=-3;
2597 fTofPIDnSigmaAboveElectronLine=5;
2598 break;
2599 case 4: // -2,3
2600 fUseTOFpid = kTRUE;
2601 fTofPIDnSigmaBelowElectronLine=-2;
2602 fTofPIDnSigmaAboveElectronLine=3;
2603 break;
e5b6e8a6 2604 case 5: // -3,3
2605 fUseTOFpid = kTRUE;
2606 fTofPIDnSigmaBelowElectronLine=-3;
2607 fTofPIDnSigmaAboveElectronLine=3;
4a0aab28 2608 break;
2609 default:
e5b6e8a6 2610 AliError(Form("TOFElectronCut not defined %d",TOFelectronPID));
4a0aab28 2611 return kFALSE;
2612 }
4a0aab28 2613 return kTRUE;
3b77b2d1 2614}
2615///________________________________________________________________________
2616Bool_t AliConversionCuts::SetQtMaxCut(Int_t QtMaxCut)
2617{ // Set Cut
4a0aab28 2618 switch(QtMaxCut){
2619 case 0: //
2620 fQtMax=1.;
e5b6e8a6 2621 fDoQtGammaSelection=kFALSE;
2622 fDoHighPtQtGammaSelection=kFALSE;
2623 fHighPtQtMax=100.;
2624 fPtBorderForQt=100.;
4a0aab28 2625 break;
2626 case 1:
2627 fQtMax=0.1;
e5b6e8a6 2628 fDoHighPtQtGammaSelection=kFALSE;
2629 fHighPtQtMax=100.;
2630 fPtBorderForQt=100.;
4a0aab28 2631 break;
2632 case 2:
2633 fQtMax=0.07;
e5b6e8a6 2634 fDoHighPtQtGammaSelection=kFALSE;
2635 fHighPtQtMax=100.;
2636 fPtBorderForQt=100.;
4a0aab28 2637 break;
2638 case 3:
2639 fQtMax=0.05;
e5b6e8a6 2640 fDoHighPtQtGammaSelection=kFALSE;
2641 fHighPtQtMax=100.;
2642 fPtBorderForQt=100.;
4a0aab28 2643 break;
2644 case 4:
2645 fQtMax=0.03;
e5b6e8a6 2646 fDoHighPtQtGammaSelection=kFALSE;
2647 fHighPtQtMax=100.;
2648 fPtBorderForQt=100.;
4a0aab28 2649 break;
e5b6e8a6 2650 case 5:
4a0aab28 2651 fQtMax=0.02;
e5b6e8a6 2652 fDoHighPtQtGammaSelection=kFALSE;
2653 fHighPtQtMax=100.;
2654 fPtBorderForQt=100.;
2655 break;
2656 case 6:
4a0aab28 2657 fQtMax=0.02;
2658 fDoHighPtQtGammaSelection=kTRUE;
2659 fHighPtQtMax=0.06;
2660 fPtBorderForQt=2.5;
e5b6e8a6 2661 break;
4a0aab28 2662 case 7:
2663 fQtMax=0.15;
e5b6e8a6 2664 fDoHighPtQtGammaSelection=kFALSE;
2665 fHighPtQtMax=100.;
2666 fPtBorderForQt=100.;
4a0aab28 2667 break;
2668 default:
2669 AliError(Form("Warning: QtMaxCut not defined %d",QtMaxCut));
2670 return kFALSE;
2671 }
2672 return kTRUE;
3b77b2d1 2673}
3b77b2d1 2674///________________________________________________________________________
e5b6e8a6 2675Bool_t AliConversionCuts::SetChi2GammaCut(Int_t chi2GammaCut)
2676{ // Set Cut
2677
2678 switch(chi2GammaCut){
2679 case 0: // 100
2680 fChi2CutConversion = 100.;
4a0aab28 2681 break;
e5b6e8a6 2682 case 1: // 50
2683 fChi2CutConversion = 50.;
4a0aab28 2684 break;
e5b6e8a6 2685 case 2: // 30
2686 fChi2CutConversion = 30.;
2687 break;
2688 case 3:
2689 fChi2CutConversion = 200.;
2690 break;
2691 case 4:
2692 fChi2CutConversion = 500.;
2693 break;
2694 case 5:
2695 fChi2CutConversion = 100000.;
2696 break;
2697 case 6:
2698 fChi2CutConversion = 5.;
2699 break;
2700 case 7:
2701 fChi2CutConversion = 10.;
2702 break;
2703 case 8:
2704 fChi2CutConversion = 20.;
2705 break;
2706 case 9:
2707 fChi2CutConversion = 15.;
4a0aab28 2708 break;
2709 default:
e5b6e8a6 2710 AliError(Form("Warning: Chi2GammaCut not defined %d",chi2GammaCut));
4a0aab28 2711 return kFALSE;
2712 }
4a0aab28 2713 return kTRUE;
3b77b2d1 2714}
3b77b2d1 2715///________________________________________________________________________
2716Bool_t AliConversionCuts::SetPsiPairCut(Int_t psiCut) {
3b77b2d1 2717
4a0aab28 2718 switch(psiCut) {
2719 case 0:
e5b6e8a6 2720 fPsiPairCut = 10000; //
4a0aab28 2721 break;
2722 case 1:
e5b6e8a6 2723 fPsiPairCut = 0.1; //
4a0aab28 2724 break;
2725 case 2:
2726 fPsiPairCut = 0.05; // Standard
2727 break;
2728 case 3:
e5b6e8a6 2729 fPsiPairCut = 0.035; //
4a0aab28 2730 break;
2731 case 4:
e5b6e8a6 2732 fPsiPairCut = 0.15; //
4a0aab28 2733 break;
2734 case 5:
e5b6e8a6 2735 fPsiPairCut = 0.2; //
4a0aab28 2736 break;
2737 case 6:
e5b6e8a6 2738 fPsiPairCut = 0.03; //
4a0aab28 2739 break;
2740 case 7:
e5b6e8a6 2741 fPsiPairCut = 0.025; //
4a0aab28 2742 break;
2743 case 8:
e5b6e8a6 2744 fPsiPairCut = 0.01; //
2745 break;
2746 case 9:
2747 fPsiPairCut = 0.5; //
4a0aab28 2748 break;
2749 default:
92efd725 2750 AliError(Form("PsiPairCut not defined %d",psiCut));
3b77b2d1 2751 return kFALSE;
4a0aab28 2752 }
3b77b2d1 2753
4a0aab28 2754 return kTRUE;
3b77b2d1 2755}
e5b6e8a6 2756///________________________________________________________________________
2757Bool_t AliConversionCuts::SetPhotonAsymmetryCut(Int_t doPhotonAsymmetryCut){
2758 // Set Cut
2759 switch(doPhotonAsymmetryCut){
2760 case 0:
2761 fDoPhotonAsymmetryCut=0;
2762 fMinPPhotonAsymmetryCut=100.;
2763 fMinPhotonAsymmetry=0.;
2764 break;
2765 case 1:
2766 fDoPhotonAsymmetryCut=1;
2767 fMinPPhotonAsymmetryCut=3.5;
2768 fMinPhotonAsymmetry=0.04;
2769 break;
2770 case 2:
2771 fDoPhotonAsymmetryCut=1;
2772 fMinPPhotonAsymmetryCut=3.5;
2773 fMinPhotonAsymmetry=0.06;
2774 break;
2775 default:
2776 AliError(Form("PhotonAsymmetryCut not defined %d",doPhotonAsymmetryCut));
2777 return kFALSE;
2778 }
2779 fCuts[kdoPhotonAsymmetryCut]=doPhotonAsymmetryCut;
2780 return kTRUE;
2781}
3b77b2d1 2782///________________________________________________________________________
2783Bool_t AliConversionCuts::SetCosPAngleCut(Int_t cosCut) {
2784
4a0aab28 2785 switch(cosCut){
2786 case 0:
e5b6e8a6 2787 fCosPAngleCut = TMath::Pi(); // -1
4a0aab28 2788 break;
2789 case 1:
e5b6e8a6 2790 fCosPAngleCut = 0.1; // 0.99500
4a0aab28 2791 break;
2792 case 2:
e5b6e8a6 2793 fCosPAngleCut = 0.05; // 0.99875
4a0aab28 2794 break;
2795 case 3:
e5b6e8a6 2796 fCosPAngleCut = 0.025; // 0.99969
4a0aab28 2797 break;
2798 case 4:
e5b6e8a6 2799 fCosPAngleCut = 0.01; // 0.99995
2800 break;
2801 case 5:
2802 fCosPAngleCut = 0.2; // 0.98007
2803 break;
2804 case 6:
2805 fCosPAngleCut = 0.5; // 0.87758
2806 break;
2807 case 7:
2808 fCosPAngleCut = 0.075; // 0.73169
4a0aab28 2809 break;
2810 default:
2811 AliError(Form("Cosine Pointing Angle cut not defined %d",cosCut));
2812 return kFALSE;
2813 }
3b77b2d1 2814
4a0aab28 2815 return kTRUE;
3b77b2d1 2816}
92efd725 2817///________________________________________________________________________
2818Bool_t AliConversionCuts::SetSharedElectronCut(Int_t sharedElec) {
2819
2820 switch(sharedElec){
4a0aab28 2821 case 0:
2822 fDoSharedElecCut = kFALSE;
2823 break;
2824 case 1:
2825 fDoSharedElecCut = kTRUE;
2826 break;
2827 default:
2828 AliError(Form("Shared Electron Cut not defined %d",sharedElec));
2829 return kFALSE;
2830 }
e5b6e8a6 2831
4a0aab28 2832 return kTRUE;
92efd725 2833}
92efd725 2834///________________________________________________________________________
2835Bool_t AliConversionCuts::SetToCloseV0sCut(Int_t toClose) {
2836
2837 switch(toClose){
2838 case 0:
2839 fDoToCloseV0sCut = kFALSE;
2840 fminV0Dist = 250;
2841 break;
2842 case 1:
2843 fDoToCloseV0sCut = kTRUE;
2844 fminV0Dist = 1;
2845 break;
2846 case 2:
2847 fDoToCloseV0sCut = kTRUE;
2848 fminV0Dist = 2;
2849 break;
2850 case 3:
2851 fDoToCloseV0sCut = kTRUE;
2852 fminV0Dist = 3;
2853 break;
2854 default:
4a0aab28 2855 AliError(Form("Shared Electron Cut not defined %d",toClose));
2856 return kFALSE;
92efd725 2857 }
2858 return kTRUE;
2859}
ca91a3e1 2860///________________________________________________________________________
e5b6e8a6 2861Bool_t AliConversionCuts::SetTRDElectronCut(Int_t TRDElectronCut)
2862{ // Set Cut
2863 switch(TRDElectronCut){
ca91a3e1 2864 case 0:
e5b6e8a6 2865 fDoTRDPID=kFALSE;
2866 break;
ca91a3e1 2867 case 1:
e5b6e8a6 2868 fDoTRDPID=kTRUE;
2869 fPIDTRDEfficiency=0.1;
2870 break;
2871 case 8:
2872 fDoTRDPID=kTRUE;
2873 fPIDTRDEfficiency=0.8;
2874 break;
2875 case 9:
2876 fDoTRDPID=kTRUE;
2877 fPIDTRDEfficiency=0.9;
2878 break;
ca91a3e1 2879 default:
e5b6e8a6 2880 AliError(Form("TRDElectronCut not defined %d",TRDElectronCut));
ca91a3e1 2881 return kFALSE;
2882 }
e5b6e8a6 2883
ca91a3e1 2884 return kTRUE;
2885}
4803eb1f 2886
2887///________________________________________________________________________
2888Bool_t AliConversionCuts::SetDCAZPhotonPrimVtxCut(Int_t DCAZPhotonPrimVtx){
2889 // Set Cut
2890 switch(DCAZPhotonPrimVtx){
2891 case 0: //
2892 fDCAZPrimVtxCut = 1000;
2893 break;
2894 case 1: //
2895 fDCAZPrimVtxCut = 10;
2896 break;
2897 case 2: //
2898 fDCAZPrimVtxCut = 5;
2899 break;
2900 case 3: //
2901 fDCAZPrimVtxCut = 4;
2902 break;
2903 case 4: //
2904 fDCAZPrimVtxCut = 3;
2905 break;
2906 case 5: //
2907 fDCAZPrimVtxCut = 2.5;
2908 break;
2909 case 6: //
2910 fDCAZPrimVtxCut = 2;
2911 break;
2912 case 7: //
2913 fDCAZPrimVtxCut = 1.5;
2914 break;
2915 case 8: //
2916 fDCAZPrimVtxCut = 1;
2917 break;
2918 case 9: //
2919 fDCAZPrimVtxCut = 0.5;
2920 break;
2921 default:
2922 cout<<"Warning: DCAZPhotonPrimVtx not defined "<<DCAZPhotonPrimVtx<<endl;
2923 return kFALSE;
2924 }
2925 return kTRUE;
2926}
2927
2928///________________________________________________________________________
2929Bool_t AliConversionCuts::SetDCARPhotonPrimVtxCut(Int_t DCARPhotonPrimVtx){
2930 // Set Cut
2931 switch(DCARPhotonPrimVtx){
2932 case 0: //
2933 fDCARPrimVtxCut = 1000;
2934 break;
2935 case 1: //
2936 fDCARPrimVtxCut = 10;
2937 break;
2938 case 2: //
2939 fDCARPrimVtxCut = 5;
2940 break;
2941 case 3: //
2942 fDCARPrimVtxCut = 4;
2943 break;
2944 case 4: //
2945 fDCARPrimVtxCut = 3;
2946 break;
2947 case 5: //
2948 fDCARPrimVtxCut = 2.5;
2949 break;
2950 case 6: //
2951 fDCARPrimVtxCut = 2;
2952 break;
2953 case 7: //
2954 fDCARPrimVtxCut = 1.5;
2955 break;
2956 case 8: //
2957 fDCARPrimVtxCut = 1;
2958 break;
2959 case 9: //
2960 fDCARPrimVtxCut = 0.5;
2961 break;
2962 default:
2963 cout<<"Warning: DCARPhotonPrimVtx not defined "<<DCARPhotonPrimVtx<<endl;
2964 return kFALSE;
2965 }
2966 return kTRUE;
2967}
2968
2969
e5b6e8a6 2970//-------------------------------------------------------------
2971Double_t AliConversionCuts::GetCentrality(AliVEvent *event)
2972{ // Get Event Centrality
2973
2974 AliESDEvent *esdEvent=dynamic_cast<AliESDEvent*>(event);
2975 if(esdEvent){
2976 AliCentrality *fESDCentrality=(AliCentrality*)esdEvent->GetCentrality();
2977
2978 if(fDetectorCentrality==0){
ae947965 2979 if (fIsHeavyIon==2){
2980 return fESDCentrality->GetCentralityPercentile("V0A"); // default for pPb
2981 } else{
2982 return fESDCentrality->GetCentralityPercentile("V0M"); // default
2983 }
e5b6e8a6 2984 }
2985 if(fDetectorCentrality==1){
2986 return fESDCentrality->GetCentralityPercentile("CL1");
2987 }
2988 }
2989
2990 AliAODEvent *aodEvent=dynamic_cast<AliAODEvent*>(event);
2991 if(aodEvent){
2992 if(aodEvent->GetHeader()){return aodEvent->GetHeader()->GetCentrality();}
2993 }
2994
2995 return -1;
2996}
2997//-------------------------------------------------------------
2998Bool_t AliConversionCuts::IsCentralitySelected(AliVEvent *event, AliVEvent *fMCEvent)
2999{ // Centrality Selection
3000 if(!fIsHeavyIon)return kTRUE;
3001
4803eb1f 3002 if(fCentralityMin == fCentralityMax ) return kTRUE;//0-100%
3003 else if(fCentralityMax==0) fCentralityMax=10; //CentralityRange = fCentralityMin-100%
3004
e5b6e8a6 3005 Double_t centrality=GetCentrality(event);
3006 if(centrality<0)return kFALSE;
3007
3008 Int_t centralityC=0;
3009 if (fModCentralityClass == 0){
3010 centralityC= Int_t(centrality/10);
3011 if(centralityC >= fCentralityMin && centralityC < fCentralityMax)
3012 return kTRUE;
3013 else return kFALSE;
3014 }
3015 else if (fModCentralityClass ==1){
3016 centralityC= Int_t(centrality);
3017 if(centralityC >= fCentralityMin*5 && centralityC < fCentralityMax*5){
3018 return kTRUE;
3019 } else return kFALSE;
3020 }
3021 else if (fModCentralityClass ==2){
3022 centralityC= Int_t(centrality);
3023 if(centralityC >= ((fCentralityMin*5)+45) && centralityC < ((fCentralityMax*5)+45))
3024 return kTRUE;
3025 else return kFALSE;
3026 }
ccfa8c0d 3027
4803eb1f 3028 Int_t nprimaryTracks = ((AliV0ReaderV1*)AliAnalysisManager::GetAnalysisManager()->GetTask("V0ReaderV1"))->GetNumberOfPrimaryTracks();
3029 Int_t PrimaryTracks10[10][2] =
3030 {
3031 {9999,9999}, // 0
3032 {1210,1210}, // 10
3033 { 817, 817}, // 20
3034 { 536, 536}, // 30
3035 { 337, 337}, // 40
3036 { 197, 197}, // 50
3037 { 106, 106}, // 60
3038 { 51, 51}, // 70
3039 { 21, 21}, // 80
3040 { 0, 0} // 90
3041 };
3042 Int_t PrimaryTracks5a[10][2] =
3043 {
3044 {9999,9999}, // 0
3045 {1485,2560}, // 5
3046 {1210,2070}, // 10
3047 { 995,1730}, // 15
3048 { 817,1415}, // 20
3049 { 666,1145}, // 25
3050 { 536, 925}, // 30
3051 { 428, 735}, // 35
3052 { 337, 557}, // 40
3053 { 260, 445} // 45
3054 };
3055 Int_t PrimaryTracks5b[10][2] =
3056 {
3057 { 260, 445}, // 45
3058 { 197, 333}, // 50 // Below still ESDS for AODs
3059 { 147, 147}, // 55
3060 { 106, 106}, // 60
3061 { 75, 75}, // 65
3062 { 51, 51}, // 70
3063 { 34, 34}, // 75
3064 { 21, 21}, // 80
3065 { 13, 13}, // 85
3066 { 0, 0} // 90
3067 };
3068
3069 Int_t column = -1;
3070 if(event->IsA()==AliESDEvent::Class()) column = 0;
3071 if(event->IsA()==AliAODEvent::Class()) column = 1;
3072
e5b6e8a6 3073 if (fModCentralityClass == 3){
3074 if(fMCEvent){
4803eb1f 3075 if(nprimaryTracks > PrimaryTracks10[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks10[fCentralityMin][column])
e5b6e8a6 3076 return kTRUE;
3077 else return kFALSE;
3078 }
3079 else{
3080 centralityC= Int_t(centrality/10);
3081 if(centralityC >= fCentralityMin && centralityC < fCentralityMax)
3082 return kTRUE;
3083 else return kFALSE;
3084 }
3085 }
3086 else if (fModCentralityClass ==4){
3087 if(fMCEvent){
4803eb1f 3088 if(nprimaryTracks > PrimaryTracks5a[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks5a[fCentralityMin][column])
e5b6e8a6 3089 return kTRUE;
3090 else return kFALSE;
3091 }
3092 else{
3093 centralityC= Int_t(centrality);
3094 if(centralityC >= fCentralityMin*5 && centralityC < fCentralityMax*5){
3095 return kTRUE;
3096 } else return kFALSE;
3097 }
3098 }
3099 else if (fModCentralityClass ==5){
3100 if(fMCEvent){
4803eb1f 3101 if(nprimaryTracks > PrimaryTracks5b[fCentralityMax][column] && nprimaryTracks <= PrimaryTracks5b[fCentralityMin][column])
e5b6e8a6 3102 return kTRUE;
3103 else return kFALSE;
3104 }
3105 else{
3106 centralityC= Int_t(centrality);
3107 if(centralityC >= ((fCentralityMin*5)+45) && centralityC < ((fCentralityMax*5)+45))
3108 return kTRUE;
3109 else return kFALSE;
3110 }
3111 }
3112
3113 return kFALSE;
3114}
3115///________________________________________________________________________
3116Bool_t AliConversionCuts::VertexZCut(AliVEvent *event){
3117 // Cut on z position of primary vertex
3118 Double_t fVertexZ=event->GetPrimaryVertex()->GetZ();
3b77b2d1 3119
e5b6e8a6 3120 if(abs(fVertexZ)>fMaxVertexZ)return kFALSE;
11c1e680 3121
0a2b2b4b 3122 if (fIsHeavyIon == 2){
3123 if(fUtils->IsFirstEventInChunk(event)) return kFALSE;
3124 if(!fUtils->IsVertexSelected2013pA(event)) return kFALSE;
3125
3126 }
11c1e680 3127
e5b6e8a6 3128 return kTRUE;
3129}
3b77b2d1 3130///________________________________________________________________________
3131
3132Int_t AliConversionCuts::GetNumberOfContributorsVtx(AliVEvent *event){
4a0aab28 3133 // returns number of contributors to the vertex
e5b6e8a6 3134
4a0aab28 3135 AliESDEvent *fESDEvent=dynamic_cast<AliESDEvent*>(event);
3136 if(fESDEvent){
e5b6e8a6 3137 if (fESDEvent->GetPrimaryVertex() != NULL){
3138 if(fESDEvent->GetPrimaryVertex()->GetNContributors()>0) {
ccfa8c0d 3139// cout << "accepted global" << fESDEvent->GetEventNumberInFile() << " with NCont: " << fESDEvent->GetPrimaryVertex()->GetNContributors() << endl;
e5b6e8a6 3140 return fESDEvent->GetPrimaryVertex()->GetNContributors();
3141 }
4a0aab28 3142 }
a280ac15 3143
e5b6e8a6 3144 if(fESDEvent->GetPrimaryVertexSPD() !=NULL){
4a0aab28 3145 if(fESDEvent->GetPrimaryVertexSPD()->GetNContributors()>0) {
ccfa8c0d 3146// cout << "accepted SPD" << fESDEvent->GetEventNumberInFile() << " with NCont: " << fESDEvent->GetPrimaryVertexSPD()->GetNContributors() << endl;
4a0aab28 3147 return fESDEvent->GetPrimaryVertexSPD()->GetNContributors();
e5b6e8a6 3148 } else {
3149 AliWarning(Form("Number of contributors from bad vertex type:: %s",fESDEvent->GetPrimaryVertex()->GetName()));
ccfa8c0d 3150// cout << "rejected " << fESDEvent->GetEventNumberInFile() << endl;
4a0aab28 3151 return 0;
3152 }
3153 }
3154 }
3155
3156 AliAODEvent *fAODEvent=dynamic_cast<AliAODEvent*>(event);
3157 if(fAODEvent){
e5b6e8a6 3158 if (fAODEvent->GetPrimaryVertex() != NULL){
3159 if(fAODEvent->GetPrimaryVertex()->GetNContributors()>0) {
3160 return fAODEvent->GetPrimaryVertex()->GetNContributors();
3161 }
a280ac15 3162 }
e5b6e8a6 3163 if(fAODEvent->GetPrimaryVertexSPD() !=NULL){
4a0aab28 3164 if(fAODEvent->GetPrimaryVertexSPD()->GetNContributors()>0) {
3165 return fAODEvent->GetPrimaryVertexSPD()->GetNContributors();
e5b6e8a6 3166 } else {
4a0aab28 3167 AliWarning(Form("Number of contributors from bad vertex type:: %s",fAODEvent->GetPrimaryVertex()->GetName()));
3168 return 0;
3169 }
3170 }
3171 }
ccfa8c0d 3172 // cout << "rejected " << fESDEvent->GetEventNumberInFile() << endl;
4a0aab28 3173 return 0;
3b77b2d1 3174}
3175
92efd725 3176///________________________________________________________________________
3b77b2d1 3177
1186afd2 3178Bool_t AliConversionCuts::IsTriggerSelected(AliVEvent *fInputEvent)
92efd725 3179{
e5b6e8a6 3180
4a0aab28 3181 AliInputEventHandler *fInputHandler=(AliInputEventHandler*)(AliAnalysisManager::GetAnalysisManager()->GetInputEventHandler());
3182
3183 UInt_t isSelected = AliVEvent::kAny;
0a2b2b4b 3184 if (fInputHandler==NULL) return kFALSE;
1186afd2 3185 if( fInputHandler->GetEventSelection() || fInputEvent->IsA()==AliAODEvent::Class()) {
ae947965 3186 if (!fTriggerSelectedManually){
3187 if (fPreSelCut) fOfflineTriggerMask = AliVEvent::kAny;
3188 else {
3189 if (fIsHeavyIon == 1) fOfflineTriggerMask = AliVEvent::kMB | AliVEvent::kCentral | AliVEvent::kSemiCentral;
3190 else if (fIsHeavyIon == 2) fOfflineTriggerMask = AliVEvent::kINT7;
3191 else fOfflineTriggerMask = AliVEvent::kMB;
3192 }
3193 }
4a0aab28 3194 // Get the actual offline trigger mask for the event and AND it with the
3195 // requested mask. If no mask requested select by default the event.
ae947965 3196// if (fPreSelCut) cout << "Trigger selected from outside: "<< fTriggerSelectedManually <<"\t Offline Trigger mask for Precut: " << fOfflineTriggerMask << endl;
3197// else cout << "Trigger selected from outside: "<< fTriggerSelectedManually <<"\t Offline Trigger mask: " << fOfflineTriggerMask << endl;
3198
4a0aab28 3199 if (fOfflineTriggerMask)
3200 isSelected = fOfflineTriggerMask & fInputHandler->IsEventSelected();
3201 }
e5b6e8a6 3202 fIsSDDFired = !(fInputHandler->IsEventSelected() & AliVEvent::kFastOnly);
4a0aab28 3203
3204 // Fill Histogram
3205 if(hTriggerClass){
11c1e680 3206 if (fIsSDDFired) hTriggerClass->Fill(33);
e5b6e8a6 3207 if (fInputHandler->IsEventSelected() & AliVEvent::kMB)hTriggerClass->Fill(0);
3208 if (fInputHandler->IsEventSelected() & AliVEvent::kINT7)hTriggerClass->Fill(1);
3209 if (fInputHandler->IsEventSelected() & AliVEvent::kMUON)hTriggerClass->Fill(2);
3210 if (fInputHandler->IsEventSelected() & AliVEvent::kHighMult)hTriggerClass->Fill(3);
3211 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC1)hTriggerClass->Fill(4);
4a0aab28 3212 if (fInputHandler->IsEventSelected() & AliVEvent::kCINT5)hTriggerClass->Fill(5);
e5b6e8a6 3213 if (fInputHandler->IsEventSelected() & AliVEvent::kCMUS5)hTriggerClass->Fill(6);
3214 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSPB)hTriggerClass->Fill(6);
e5b6e8a6 3215 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSH7)hTriggerClass->Fill(7);
3216 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSHPB)hTriggerClass->Fill(7);
3217 if (fInputHandler->IsEventSelected() & AliVEvent::kMUL7)hTriggerClass->Fill(8);
3218 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikePB)hTriggerClass->Fill(8);
3219 if (fInputHandler->IsEventSelected() & AliVEvent::kMUU7)hTriggerClass->Fill(9);
3220 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikePB)hTriggerClass->Fill(9);
3221 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7)hTriggerClass->Fill(10);
3222 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC8)hTriggerClass->Fill(10);
3223 if (fInputHandler->IsEventSelected() & AliVEvent::kMUS7)hTriggerClass->Fill(11);
3224 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI1)hTriggerClass->Fill(12);
3225 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI7)hTriggerClass->Fill(13);
3226 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI8)hTriggerClass->Fill(13);
3227 if (fInputHandler->IsEventSelected() & AliVEvent::kPHOSPb)hTriggerClass->Fill(13);
3228 if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEJE)hTriggerClass->Fill(14);
3229 if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEGA)hTriggerClass->Fill(15);
3230 if (fInputHandler->IsEventSelected() & AliVEvent::kCentral)hTriggerClass->Fill(16);
3231 if (fInputHandler->IsEventSelected() & AliVEvent::kSemiCentral)hTriggerClass->Fill(17);
3232 if (fInputHandler->IsEventSelected() & AliVEvent::kDG5)hTriggerClass->Fill(18);
3233 if (fInputHandler->IsEventSelected() & AliVEvent::kZED)hTriggerClass->Fill(19);
3234 if (fInputHandler->IsEventSelected() & AliVEvent::kSPI7)hTriggerClass->Fill(20);
3235 if (fInputHandler->IsEventSelected() & AliVEvent::kSPI)hTriggerClass->Fill(20);
3236 if (fInputHandler->IsEventSelected() & AliVEvent::kINT8)hTriggerClass->Fill(21);
3237 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleLowPt8)hTriggerClass->Fill(22);
3238 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleHighPt8)hTriggerClass->Fill(23);
3239 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikeLowPt8)hTriggerClass->Fill(24);
3240 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt8)hTriggerClass->Fill(25);
3241 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt0)hTriggerClass->Fill(26);
3242 if (fInputHandler->IsEventSelected() & AliVEvent::kUserDefined)hTriggerClass->Fill(27);
11c1e680 3243 if (fInputHandler->IsEventSelected() & AliVEvent::kTRD)hTriggerClass->Fill(28);
3244 if (fInputHandler->IsEventSelected() & AliVEvent::kFastOnly)hTriggerClass->Fill(29);
3245 if (fInputHandler->IsEventSelected() & AliVEvent::kAnyINT)hTriggerClass->Fill(30);
3246 if (fInputHandler->IsEventSelected() & AliVEvent::kAny)hTriggerClass->Fill(31);
3247 if (!fInputHandler->IsEventSelected()) hTriggerClass->Fill(34);
a280ac15 3248 }
3249
3250 if(hTriggerClassSelected && isSelected){
11c1e680 3251 if (!fIsSDDFired) hTriggerClassSelected->Fill(33);
a280ac15 3252 if (fInputHandler->IsEventSelected() & AliVEvent::kMB)hTriggerClassSelected->Fill(0);
3253 if (fInputHandler->IsEventSelected() & AliVEvent::kINT7)hTriggerClassSelected->Fill(1);
3254 if (fInputHandler->IsEventSelected() & AliVEvent::kMUON)hTriggerClassSelected->Fill(2);
3255 if (fInputHandler->IsEventSelected() & AliVEvent::kHighMult)hTriggerClassSelected->Fill(3);
3256 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC1)hTriggerClassSelected->Fill(4);
3257 if (fInputHandler->IsEventSelected() & AliVEvent::kCINT5)hTriggerClassSelected->Fill(5);
3258 if (fInputHandler->IsEventSelected() & AliVEvent::kCMUS5)hTriggerClassSelected->Fill(6);
3259 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSPB)hTriggerClassSelected->Fill(6);
3260 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSH7)hTriggerClassSelected->Fill(7);
3261 if (fInputHandler->IsEventSelected() & AliVEvent::kMUSHPB)hTriggerClassSelected->Fill(7);
3262 if (fInputHandler->IsEventSelected() & AliVEvent::kMUL7)hTriggerClassSelected->Fill(8);
3263 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikePB)hTriggerClassSelected->Fill(8);
3264 if (fInputHandler->IsEventSelected() & AliVEvent::kMUU7)hTriggerClassSelected->Fill(9);
3265 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikePB)hTriggerClassSelected->Fill(9);
3266 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC7)hTriggerClassSelected->Fill(10);
3267 if (fInputHandler->IsEventSelected() & AliVEvent::kEMC8)hTriggerClassSelected->Fill(10);
3268 if (fInputHandler->IsEventSelected() & AliVEvent::kMUS7)hTriggerClassSelected->Fill(11);
3269 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI1)hTriggerClassSelected->Fill(12);
3270 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI7)hTriggerClassSelected->Fill(13);
3271 if (fInputHandler->IsEventSelected() & AliVEvent::kPHI8)hTriggerClassSelected->Fill(13);
3272 if (fInputHandler->IsEventSelected() & AliVEvent::kPHOSPb)hTriggerClassSelected->Fill(13);
3273 if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEJE)hTriggerClassSelected->Fill(14);
3274 if (fInputHandler->IsEventSelected() & AliVEvent::kEMCEGA)hTriggerClassSelected->Fill(15);
3275 if (fInputHandler->IsEventSelected() & AliVEvent::kCentral)hTriggerClassSelected->Fill(16);
3276 if (fInputHandler->IsEventSelected() & AliVEvent::kSemiCentral)hTriggerClassSelected->Fill(17);
3277 if (fInputHandler->IsEventSelected() & AliVEvent::kDG5)hTriggerClassSelected->Fill(18);
3278 if (fInputHandler->IsEventSelected() & AliVEvent::kZED)hTriggerClassSelected->Fill(19);
3279 if (fInputHandler->IsEventSelected() & AliVEvent::kSPI7)hTriggerClassSelected->Fill(20);
3280 if (fInputHandler->IsEventSelected() & AliVEvent::kSPI)hTriggerClassSelected->Fill(20);
3281 if (fInputHandler->IsEventSelected() & AliVEvent::kINT8)hTriggerClassSelected->Fill(21);
3282 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleLowPt8)hTriggerClassSelected->Fill(22);
3283 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonSingleHighPt8)hTriggerClassSelected->Fill(23);
3284 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonLikeLowPt8)hTriggerClassSelected->Fill(24);
3285 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt8)hTriggerClassSelected->Fill(25);
3286 if (fInputHandler->IsEventSelected() & AliVEvent::kMuonUnlikeLowPt0)hTriggerClassSelected->Fill(26);
3287 if (fInputHandler->IsEventSelected() & AliVEvent::kUserDefined)hTriggerClassSelected->Fill(27);
11c1e680 3288 if (fInputHandler->IsEventSelected() & AliVEvent::kTRD)hTriggerClassSelected->Fill(28);
3289 if (fInputHandler->IsEventSelected() & AliVEvent::kFastOnly)hTriggerClassSelected->Fill(29);
3290 if (fInputHandler->IsEventSelected() & AliVEvent::kAnyINT)hTriggerClassSelected->Fill(30);
3291 if (fInputHandler->IsEventSelected() & AliVEvent::kAny)hTriggerClassSelected->Fill(31);
4a0aab28 3292 }
3b77b2d1 3293
e5b6e8a6 3294 if(!isSelected)return kFALSE;
a280ac15 3295
4a0aab28 3296 return kTRUE;
3b77b2d1 3297
92efd725 3298}
3b77b2d1 3299
3300///________________________________________________________________________
3301Int_t AliConversionCuts::GetFirstTPCRow(Double_t radius){
4a0aab28 3302 // Get first TPC row
e5b6e8a6 3303 Int_t firstTPCRow = 0;
3304 Double_t radiusI = 84.8;
3305 Double_t radiusO = 134.6;
4a0aab28 3306 Double_t radiusOB = 198.;
e5b6e8a6 3307 Double_t rSizeI = 0.75;
3308 Double_t rSizeO = 1.;
3309 Double_t rSizeOB = 1.5;
3310 Int_t nClsI = 63;
3311 Int_t nClsIO = 127;
4a0aab28 3312
3313 if(radius <= radiusI){
3314 return firstTPCRow;
3315 }
3316 if(radius>radiusI && radius<=radiusO){
3317 firstTPCRow = (Int_t)((radius-radiusI)/rSizeI);
3318 }
3319 if(radius>radiusO && radius<=radiusOB){
3320 firstTPCRow = (Int_t)(nClsI+(radius-radiusO)/rSizeO);
3321 }
3322
3323 if(radius>radiusOB){
3324 firstTPCRow =(Int_t)(nClsIO+(radius-radiusOB)/rSizeOB);
3325 }
3326
3327 return firstTPCRow;
3b77b2d1 3328}
3329
92efd725 3330Bool_t AliConversionCuts::CosinePAngleCut(const AliConversionPhotonBase * photon, AliVEvent * event) const {
4a0aab28 3331 ///Check if passes cosine of pointing angle cut
92efd725 3332 if(GetCosineOfPointingAngle(photon, event) < (TMath::Cos(fCosPAngleCut))){
3333 return kFALSE;
3334 }
3335 return kTRUE;
3b77b2d1 3336}
3337
92efd725 3338Double_t AliConversionCuts::GetCosineOfPointingAngle( const AliConversionPhotonBase * photon, AliVEvent * event) const{
e5b6e8a6 3339 // calculates the pointing angle of the recalculated V0
3b77b2d1 3340
92efd725 3341 Double_t momV0[3] = {0,0,0};
3342 if(event->IsA()==AliESDEvent::Class()){
3343 AliESDEvent *esdEvent = dynamic_cast<AliESDEvent*>(event);
e5b6e8a6 3344 if(!esdEvent) return -999;
92efd725 3345 AliESDv0 *v0 = esdEvent->GetV0(photon->GetV0Index());
e5b6e8a6 3346 if(!v0) return -999;
92efd725 3347 v0->GetPxPyPz(momV0[0],momV0[1],momV0[2]);
3348 }
3349 if(event->IsA()==AliAODEvent::Class()){
3350 momV0[0] = photon->GetPx();
3351 momV0[1] = photon->GetPy();
3352 momV0[2] = photon->GetPz();
3353 }
e5b6e8a6 3354
92efd725 3355 //Double_t momV0[3] = { photon->GetPx(), photon->GetPy(), photon->GetPz() }; //momentum of the V0
e5b6e8a6 3356 Double_t PosV0[3] = { photon->GetConversionX() - event->GetPrimaryVertex()->GetX(),
3357 photon->GetConversionY() - event->GetPrimaryVertex()->GetY(),
92efd725 3358 photon->GetConversionZ() - event->GetPrimaryVertex()->GetZ() }; //Recalculated V0 Position vector
e5b6e8a6 3359
3b77b2d1 3360 Double_t momV02 = momV0[0]*momV0[0] + momV0[1]*momV0[1] + momV0[2]*momV0[2];
3361 Double_t PosV02 = PosV0[0]*PosV0[0] + PosV0[1]*PosV0[1] + PosV0[2]*PosV0[2];
3362
e5b6e8a6 3363
3364 Double_t cosinePointingAngle = -999;
3365 if(momV02*PosV02 > 0.0)
3366 cosinePointingAngle = (PosV0[0]*momV0[0] + PosV0[1]*momV0[1] + PosV0[2]*momV0[2] ) / TMath::Sqrt(momV02 * PosV02);
3367
3b77b2d1 3368 return cosinePointingAngle;
3369}
3370
e5b6e8a6 3371///________________________________________________________________________
92efd725 3372Bool_t AliConversionCuts::PsiPairCut(const AliConversionPhotonBase * photon) const {
3b77b2d1 3373
4a0aab28 3374 if(photon->GetPsiPair() > fPsiPairCut){
3375 return kFALSE;}
3376 else{return kTRUE;}
3b77b2d1 3377
3b77b2d1 3378}
3379
3380///________________________________________________________________________
3381TString AliConversionCuts::GetCutNumber(){
4a0aab28 3382 // returns TString with current cut number
3383 TString a(kNCuts);
3384 for(Int_t ii=0;ii<kNCuts;ii++){
92efd725 3385 a.Append(Form("%d",fCuts[ii]));
4a0aab28 3386 }
3387 return a;
3b77b2d1 3388}
3389
92efd725 3390///________________________________________________________________________
3391void AliConversionCuts::FillElectonLabelArray(AliAODConversionPhoton* photon, Int_t nV0){
e5b6e8a6 3392
92efd725 3393 Int_t posLabel = photon->GetTrackLabelPositive();
3394 Int_t negLabel = photon->GetTrackLabelNegative();
e5b6e8a6 3395
92efd725 3396 fElectronLabelArray[nV0*2] = posLabel;
3397 fElectronLabelArray[(nV0*2)+1] = negLabel;
3398}
3399///________________________________________________________________________
3400Bool_t AliConversionCuts::RejectSharedElectronV0s(AliAODConversionPhoton* photon, Int_t nV0, Int_t nV0s){
3b77b2d1 3401
92efd725 3402 Int_t posLabel = photon->GetTrackLabelPositive();
3403 Int_t negLabel = photon->GetTrackLabelNegative();
e5b6e8a6 3404
92efd725 3405 for(Int_t i = 0; i<nV0s*2;i++){
3406 if(i==nV0*2) continue;
3407 if(i==(nV0*2)+1) continue;
3408 if(fElectronLabelArray[i] == posLabel){
3409 return kFALSE;}
3410 if(fElectronLabelArray[i] == negLabel){
3411 return kFALSE;}
3412 }
3b77b2d1 3413
92efd725 3414 return kTRUE;
3415}
3b77b2d1 3416///________________________________________________________________________
92efd725 3417Bool_t AliConversionCuts::RejectToCloseV0s(AliAODConversionPhoton* photon, TList *photons, Int_t nV0){
3418
3419
3420 Double_t posX = photon->GetConversionX();
3421 Double_t posY = photon->GetConversionY();
3422 Double_t posZ = photon->GetConversionZ();
3423
3424 for(Int_t i = 0;i<photons->GetEntries();i++){
3425 if(nV0 == i) continue;
3426 AliAODConversionPhoton *photonComp = (AliAODConversionPhoton*) photons->At(i);
3427 Double_t posCompX = photonComp->GetConversionX();
3428 Double_t posCompY = photonComp->GetConversionY();
3429 Double_t posCompZ = photonComp->GetConversionZ();
e5b6e8a6 3430
92efd725 3431 Double_t dist = pow((posX - posCompX),2)+pow((posY - posCompY),2)+pow((posZ - posCompZ),2);
3432
3433 if(dist < fminV0Dist*fminV0Dist){
3434 if(photon->GetChi2perNDF() < photonComp->GetChi2perNDF()) return kTRUE;
3435 else {
3436 return kFALSE;}
3437 }
e5b6e8a6 3438
92efd725 3439 }
3440 return kTRUE;
3b77b2d1 3441}
ca91a3e1 3442///________________________________________________________________________
ae947965 3443void AliConversionCuts::GetNotRejectedParticles(Int_t rejection, TList *HeaderList, AliVEvent *MCEvent){
e5b6e8a6 3444
4803eb1f 3445
3446
ca91a3e1 3447 if(fNotRejectedStart){
3448 delete[] fNotRejectedStart;
3449 fNotRejectedStart = NULL;
3450 }
3451 if(fNotRejectedEnd){
3452 delete[] fNotRejectedEnd;
3453 fNotRejectedEnd = NULL;
3454 }
3455 if(fGeneratorNames){
3456 delete[] fGeneratorNames;
3457 fGeneratorNames = NULL;
3458 }
e5b6e8a6 3459
ca91a3e1 3460 if(rejection == 0) return; // No Rejection
ae947965 3461
3462 AliGenCocktailEventHeader *cHeader = 0x0;
3463 AliAODMCHeader *cHeaderAOD = 0x0;
3464 Bool_t headerFound = kFALSE;
3465
3466 if(MCEvent->IsA()==AliMCEvent::Class()){
3467 cHeader = dynamic_cast<AliGenCocktailEventHeader*>(dynamic_cast<AliMCEvent*>(MCEvent)->GenEventHeader());
3468 if(cHeader) headerFound = kTRUE;
3469 }
4803eb1f 3470 if(MCEvent->IsA()==AliAODEvent::Class()){ // MCEvent is a AODEvent in case of AOD
ae947965 3471 cHeaderAOD = dynamic_cast<AliAODMCHeader*>(MCEvent->FindListObject(AliAODMCHeader::StdBranchName()));
3472 if(cHeaderAOD) headerFound = kTRUE;
3473 }
3474
3475 if(headerFound){
3476 TList *genHeaders = 0x0;
3477 if(cHeader) genHeaders = cHeader->GetHeaders();
3478 if(cHeaderAOD){
3479 genHeaders = cHeaderAOD->GetCocktailHeaders();
3480 if(genHeaders->GetEntries()==1){
3481 SetRejectExtraSignalsCut(0);
3482 return;
3483 }
3484 }
ca91a3e1 3485 AliGenEventHeader* gh = 0;
3486 fnHeaders = 0;
e5b6e8a6 3487 if(rejection == 1 || rejection == 3) fnHeaders = 1; // MinBiasHeader
ca91a3e1 3488 if(rejection == 2){ // TList of Headers Names
3489 for(Int_t i = 0; i<genHeaders->GetEntries();i++){
3490 gh = (AliGenEventHeader*)genHeaders->At(i);
3491 TString GeneratorName = gh->GetName();
3492 for(Int_t j = 0; j<HeaderList->GetEntries();j++){
3493 TString GeneratorInList = ((TObjString*)HeaderList->At(j))->GetString();
3494 if(GeneratorName.CompareTo(GeneratorInList) == 0){
3495 fnHeaders++;
3496 continue;
3497 }
3498 }
3499 }
3500 }
3501
3502 fNotRejectedStart = new Int_t[fnHeaders];
3503 fNotRejectedEnd = new Int_t[fnHeaders];
3504 fGeneratorNames = new TString[fnHeaders];
3505
e5b6e8a6 3506 if(rejection == 1 || rejection == 3){
ca91a3e1 3507 fNotRejectedStart[0] = 0;
e5b6e8a6 3508 fNotRejectedEnd[0] = ((AliGenEventHeader*)genHeaders->At(0))->NProduced()-1;
ca91a3e1 3509 fGeneratorNames[0] = ((AliGenEventHeader*)genHeaders->At(0))->GetName();
ca91a3e1 3510 return;
3511 }
e5b6e8a6 3512
ca91a3e1 3513 Int_t firstindex = 0;
e5b6e8a6 3514 Int_t lastindex = -1;
ca91a3e1 3515 Int_t nummer = 0;
3516 for(Int_t i = 0; i<genHeaders->GetEntries();i++){
3517 gh = (AliGenEventHeader*)genHeaders->At(i);
3518 TString GeneratorName = gh->GetName();
3519 lastindex = lastindex + gh->NProduced();
3520 for(Int_t j = 0; j<HeaderList->GetEntries();j++){
3521 TString GeneratorInList = ((TObjString*)HeaderList->At(j))->GetString();
3522 if(GeneratorName.CompareTo(GeneratorInList) == 0){
3523 fNotRejectedStart[nummer] = firstindex;
3524 fNotRejectedEnd[nummer] = lastindex;
3525 fGeneratorNames[nummer] = GeneratorName;
4803eb1f 3526 //cout << "Number of particles produced for: " << i << "\t" << GeneratorName.Data() << "\t" << lastindex-firstindex+1 << endl;
ca91a3e1 3527 nummer++;
3528 continue;
3529 }
3530 }
3531 firstindex = firstindex + gh->NProduced();
3532 }
a280ac15 3533 } else { // No Cocktail Header Found
ca91a3e1 3534 fNotRejectedStart = new Int_t[1];
3535 fNotRejectedEnd = new Int_t[1];
e5b6e8a6 3536
3537 fnHeaders = 1;
ca91a3e1 3538 fNotRejectedStart[0] = 0;
ae947965 3539 fNotRejectedEnd[0] = static_cast<AliMCEvent*>(MCEvent)->Stack()->GetNprimary()-1;
a280ac15 3540 fGeneratorNames = new TString[1];
3541 fGeneratorNames[0] = "NoCocktailGeneratorFound";
11c1e680 3542
ae947965 3543 AliGenPythiaEventHeader *mcHeaderPythia = dynamic_cast<AliGenPythiaEventHeader*>(static_cast<AliMCEvent*>(MCEvent)->GenEventHeader());
a280ac15 3544 if (mcHeaderPythia) fGeneratorNames[0] = "NoCocktailGeneratorFound_Pythia";
ae947965 3545 AliGenDPMjetEventHeader *mcHeaderPhojet = dynamic_cast<AliGenDPMjetEventHeader*>(static_cast<AliMCEvent*>(MCEvent)->GenEventHeader());
a280ac15 3546 if (mcHeaderPhojet) fGeneratorNames[0] = "NoCocktailGeneratorFound_Phojet";
ae947965 3547 AliGenHijingEventHeader *mcHeaderHijing = dynamic_cast<AliGenHijingEventHeader*>(static_cast<AliMCEvent*>(MCEvent)->GenEventHeader());
a280ac15 3548 if (mcHeaderHijing) fGeneratorNames[0] = "NoCocktailGeneratorFound_Hijing";
3549
e5b6e8a6 3550 SetRejectExtraSignalsCut(0);
ca91a3e1 3551 }
e5b6e8a6 3552
ca91a3e1 3553}
3554//_________________________________________________________________________
ae947965 3555Int_t AliConversionCuts::IsParticleFromBGEvent(Int_t index, AliStack *MCStack, AliVEvent *InputEvent){
a280ac15 3556
3557 // Not Accepted == kFALSE == 0
3558 // Accepted == kTRUE == 1
3559 // FirstHeader == kTRUE == 3
a280ac15 3560 if(index < 0) return 0; // No Particle
e5b6e8a6 3561
a280ac15 3562 Int_t accepted = 0;
ae947965 3563 if(!InputEvent || InputEvent->IsA()==AliESDEvent::Class()){
3564 if( index >= MCStack->GetNprimary()){ // Secondary Particle
3565 if( ((TParticle*)MCStack->Particle(index))->GetMother(0) < 0) return 1; // Secondary Particle without Mother??
3566 return IsParticleFromBGEvent(((TParticle*)MCStack->Particle(index))->GetMother(0),MCStack,InputEvent);
3567 }
3568 for(Int_t i = 0;i<fnHeaders;i++){
3569 if(index >= fNotRejectedStart[i] && index <= fNotRejectedEnd[i]){
3570 accepted = 1;
3571 if(i == 0) accepted = 2; // MB Header
3572 }
3573 }
ca91a3e1 3574 }
ae947965 3575 else if(InputEvent->IsA()==AliAODEvent::Class()){
3576 TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(InputEvent->FindListObject(AliAODMCParticle::StdBranchName()));
3577 AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index));
4803eb1f 3578 if(!aodMCParticle) return 1; // Photon Without a Mother ? --> Accepted
ae947965 3579 if(!aodMCParticle->IsPrimary()){
3580 if( aodMCParticle->GetMother() < 0) return 1;// Secondary Particle without Mother??
3581 return IsParticleFromBGEvent(aodMCParticle->GetMother(),MCStack,InputEvent);
3582 }
3583 index = abs(static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index))->GetLabel());
3584 for(Int_t i = 0;i<fnHeaders;i++){
3585 if(index >= fNotRejectedStart[i] && index <= fNotRejectedEnd[i]){
3586 accepted = 1;
3587 if(i == 0) accepted = 2; // MB Header
3588 }
e5b6e8a6 3589 }
ca91a3e1 3590 }
e5b6e8a6 3591
ca91a3e1 3592 return accepted;
3593}
e5b6e8a6 3594//_________________________________________________________________________
3595Int_t AliConversionCuts::IsEventAcceptedByConversionCut(AliConversionCuts *ReaderCuts, AliVEvent *InputEvent, AliMCEvent *MCEvent, Bool_t isHeavyIon){
3596
1186afd2 3597 if ( !IsTriggerSelected(InputEvent) )
ae947965 3598 return 3;
3599
e5b6e8a6 3600 if(isHeavyIon && !(IsCentralitySelected(InputEvent,MCEvent)))
3601 return 1; // Check Centrality --> Not Accepted => eventQuality = 1
a280ac15 3602
e5b6e8a6 3603 if(!isHeavyIon && GetIsFromPileup()){
3604 if(InputEvent->IsPileupFromSPD(3,0.8,3.,2.,5.)){
ccfa8c0d 3605
e5b6e8a6 3606 return 6; // Check Pileup --> Not Accepted => eventQuality = 6
3607 }
3608 }
ae947965 3609
e5b6e8a6 3610 Bool_t hasV0And = ReaderCuts->HasV0AND();
3611 Bool_t isSDDFired = ReaderCuts->IsSDDFired();
3612 if( (IsSpecialTrigger() == 2 || IsSpecialTrigger() == 3) && !isSDDFired && !MCEvent)
3613 return 7; // With SDD requested but no fired
a280ac15 3614
e5b6e8a6 3615 if( (IsSpecialTrigger() == 1 || IsSpecialTrigger() == 3) && !hasV0And)
3616 return 8; // V0AND requested but no fired
ae947965 3617
e5b6e8a6 3618
3619 return 0;
3620}
ccfa8c0d 3621//_________________________________________________________________________
ae947965 3622Float_t AliConversionCuts::GetWeightForMeson(TString period, Int_t index, AliStack *MCStack, AliVEvent *InputEvent){
3623 if (!(period.CompareTo("LHC12f1a") == 0 || period.CompareTo("LHC12f1b") == 0 || period.CompareTo("LHC12i3") == 0 || period.CompareTo("LHC11a10a") == 0 || period.CompareTo("LHC11a10b") == 0 || period.CompareTo("LHC11a10b_bis") == 0 || period.CompareTo("LHC11a10a_bis") == 0 || period.CompareTo("LHC11a10b_plus") == 0 || period.CompareTo("LHC13d2") == 0)) return 1.;
3624
a280ac15 3625 Int_t kCaseGen = 0;
3626 for (Int_t i = 0; i < fnHeaders; i++){
3627 if (index >= fNotRejectedStart[i] && index < fNotRejectedEnd[i]+1){
a280ac15 3628 if (fGeneratorNames[i].CompareTo("Pythia") == 0){
3629 kCaseGen = 1;
3630 } else if (fGeneratorNames[i].CompareTo("DPMJET") == 0){
3631 kCaseGen = 2;
ae947965 3632 } else if (fGeneratorNames[i].CompareTo("HIJING") == 0 ||
3633 fGeneratorNames[i].CompareTo("Hijing") == 0 ||
3634 fGeneratorNames[i].Contains("hijing")){
a280ac15 3635 kCaseGen = 3;
3636 } else if (fGeneratorNames[i].CompareTo("BOX") == 0){
3637 kCaseGen = 4;
3638 } else if (fGeneratorNames[i].CompareTo("PARAM") == 0){
3639 kCaseGen = 5;
3640 } else if (fGeneratorNames[i].CompareTo("NoCocktailGeneratorFound") == 0){
3641 kCaseGen = 6;
3642 } else if (fGeneratorNames[i].CompareTo("NoCocktailGeneratorFound_Pythia") == 0){
3643 kCaseGen = 1;
3644 } else if (fGeneratorNames[i].CompareTo("NoCocktailGeneratorFound_Phojet") == 0){
3645 kCaseGen = 2;
3646 } else if (fGeneratorNames[i].CompareTo("NoCocktailGeneratorFound_Hijing") == 0){
3647 kCaseGen = 3;
11c1e680 3648 }
a280ac15 3649 }
3650 }
ae947965 3651 if (kCaseGen == 0) return 1;
3652
11c1e680 3653
ae947965 3654 Double_t mesonPt = 0;
3655 Double_t mesonMass = 0;
3656 Int_t PDGCode = 0;
3657 if(!InputEvent || InputEvent->IsA()==AliESDEvent::Class()){
3658 mesonPt = ((TParticle*)MCStack->Particle(index))->Pt();
3659 mesonMass = ((TParticle*)MCStack->Particle(index))->GetCalcMass();
3660 PDGCode = ((TParticle*)MCStack->Particle(index))->GetPdgCode();
3661 }
3662 else if(InputEvent->IsA()==AliAODEvent::Class()){
3663 TClonesArray *AODMCTrackArray = dynamic_cast<TClonesArray*>(InputEvent->FindListObject(AliAODMCParticle::StdBranchName()));
3664 AliAODMCParticle *aodMCParticle = static_cast<AliAODMCParticle*>(AODMCTrackArray->At(index));
3665 mesonPt = aodMCParticle->Pt();
3666 mesonMass = aodMCParticle->GetCalcMass();
3667 PDGCode = aodMCParticle->GetPdgCode();
3668 }
3669
3670 Float_t functionResultMC = 1.;
3671 if (kCaseGen == 1){ // Pythia 6
a280ac15 3672 Float_t dNdyMC = 2.1462;
3673 Float_t nMC = 7.06055;
3674 Float_t tMC = 0.12533;
ae947965 3675 if ( PDGCode == 111){
a280ac15 3676 dNdyMC = 2.1462;
3677 nMC = 7.06055;
3678 tMC = 0.12533;
ae947965 3679 } else if ( PDGCode == 221){
a280ac15 3680 dNdyMC = 0.2357;
3681 nMC = 5.9105;
3682 tMC = 0.1525;
3683 }
ae947965 3684 functionResultMC = dNdyMC / ( 2 * TMath::Pi())*(nMC-1.)*(nMC-2.) / (nMC*tMC*(nMC*tMC+mesonMass*(nMC-2.))) * TMath::Power(1.+(TMath::Sqrt(mesonPt*mesonPt+mesonMass*mesonMass)-mesonMass)/(nMC*tMC), -nMC);
3685 } else if (kCaseGen == 2){ // Phojet
a280ac15 3686 Float_t dNdyMC = 2.35978;
3687 Float_t nMC = 6.81795;
3688 Float_t tMC = 0.11492;
ae947965 3689 if ( PDGCode == 111){
a280ac15 3690 dNdyMC = 2.35978;
3691 nMC = 6.81795;
3692 tMC = 0.11492;
ae947965 3693 } else if ( PDGCode == 221){
a280ac15 3694 dNdyMC = 0.3690;
3695 nMC = 5.55809;
3696 tMC = 0.13387;
3697 }
ae947965 3698 functionResultMC = dNdyMC / ( 2 * TMath::Pi())*(nMC-1.)*(nMC-2.) / (nMC*tMC*(nMC*tMC+mesonMass*(nMC-2.))) * TMath::Power(1.+(TMath::Sqrt(mesonPt*mesonPt+mesonMass*mesonMass)-mesonMass)/(nMC*tMC), -nMC);
3699 } else if (kCaseGen == 4){ // BOX generators pp
3700// functionResultMC = 1./sqrt(1.-mesonMass*mesonMass/((mesonMass*mesonMass+mesonPt*mesonPt)*cosh(mesonY)*cosh(mesonY)));
a280ac15 3701 Float_t a = 0.23437;
3702 Float_t b = 5.6661;
3703 Float_t c = -1430.5863;
3704 Float_t d = -0.6966624;
3705 Float_t e = 252.3742;
ae947965 3706 if ( PDGCode == 111){
a280ac15 3707 a = 0.23437;
3708 b = 5.6661;
3709 c = -1430.5863;
3710 d = -0.6966624;
3711 e = 252.3742;
ae947965 3712 } else if ( PDGCode == 221){
a280ac15 3713 a = 0.10399;
3714 b = 4.35311;
3715 c = -12.17723;
3716 d = -0.01172;
3717 e =1.85140;
3718 }
ae947965 3719 functionResultMC = a*TMath::Power(mesonPt,-1.*(b+c/(TMath::Power(mesonPt,d)+e)))*1./mesonPt *1./1.6 *1./(2.* TMath::Pi());
4803eb1f 3720// cout << functionResultMC << endl;
ae947965 3721 } else if (kCaseGen == 3 ){ // HIJING
3722 if ( PDGCode == 111 && fDoReweightHistoMCPi0 && hReweightMCHistPi0!= 0x0){
3723 functionResultMC = hReweightMCHistPi0->Interpolate(mesonPt);
3724 }
3725 if ( PDGCode == 310 && fDoReweightHistoMCK0s && hReweightMCHistK0s!= 0x0){
3726 functionResultMC = hReweightMCHistK0s->Interpolate(mesonPt);
3727 }
a280ac15 3728 }
11c1e680 3729
ae947965 3730 Float_t functionResultData = 1;
3731 if (kCaseGen == 1 || kCaseGen == 2 || kCaseGen == 4 ){
3732 Float_t dNdyData = 2.2328;
3733 Float_t nData = 7.1473;
3734 Float_t tData = 0.1346;
3735 if ( PDGCode == 111){
a280ac15 3736 dNdyData = 2.2328;
3737 nData = 7.1473;
3738 tData = 0.1346;
ae947965 3739 } else if ( PDGCode == 221){
a280ac15 3740 dNdyData = 0.38992; //be careful this fit is not optimal, eta in data still has problems
3741 nData = 5.72778;
3742 tData = 0.13835;
3743 }
ae947965 3744 functionResultData = dNdyData / ( 2 * TMath::Pi())*(nData-1.)*(nData-2.) / (nData*tData*(nData*tData+mesonMass*(nData-2.))) * TMath::Power(1.+(TMath::Sqrt(mesonPt*mesonPt+mesonMass*mesonMass)-mesonMass)/(nData*tData), -nData);
4803eb1f 3745// cout << functionResultData << endl;
ae947965 3746 } else {
3747 Float_t a = 0.;
3748 Float_t b = 0.;
3749 Float_t c = 0.;
3750 Float_t d = 0.;
3751 Float_t e = 0.;
3752 if ( PDGCode == 111 ){
3753 if (fModCentralityClass == 1 && fCentralityMin == 0 && fCentralityMax == 1 ){ // 0-5 % PbPb
3754 a = 25.8747458223;
3755 b = 5.8761820045;
3756 c = -33.9928191673;
3757 d = 3.0731850142;
3758 e = 13.2500447620;
3759 } else if (fModCentralityClass == 1 && fCentralityMin == 1 && fCentralityMax == 2){ // 5-10% PbPb
3760 a = 21.7518148922;
3761 b = 5.8441200081;
3762 c = -17.1497051691;
3763 d = 2.3799090842;
3764 e = 5.4346404718;
3765 } else if (fModCentralityClass == 0 && fCentralityMin == 0 && fCentralityMax == 1){ // 0-10% PbPb
3766 a = 22.9852133622;
3767 b = 5.8602063916;
3768 c = -17.0992478654;
3769 d = 2.4426218039;
3770 e = 5.1194526345;
3771 } else if (fModCentralityClass == 0 && fCentralityMin == 1 && fCentralityMax == 2){ // 10-20% PbPb
3772 a = 19.3237333776;
3773 b = 5.8145906958;
3774 c = -13.8316665424;
3775 d = 2.3737630637;
3776 e = 4.7690300693;
3777 } else if (fModCentralityClass == 0 && fCentralityMin == 2 && fCentralityMax == 4){ // 20-40% PbPb
3778 a = 11.2656032751;
3779 b = 5.8003194354;
3780 c = -13.3936105929;
3781 d = 2.3371452334;
3782 e = 4.4726244958;
3783 } else if (fModCentralityClass == 0 && fCentralityMin == 4 && fCentralityMax == 6){ // 40-60% PbPb
3784 a = 4.1578154081;
3785 b = 5.6450005163;
3786 c = -8.4309375240;
3787 d = 1.8918308704;
3788 e = 2.9429194709;
3789 } else if (fModCentralityClass == 0 && fCentralityMin == 6 && fCentralityMax == 8){ // 60-80% PbPb
3790 a = 1.0635443810;
3791 b = 5.1337469970;
3792 c = -8.5906997238;
3793 d = 2.9794995997;
3794 e = 3.9294980048;
3795 } else if (fModCentralityClass == 0 && fCentralityMin == 0 && fCentralityMax == 2){ // 0-20% PbPb
3796 a = 21.7018745556;
3797 b = 5.9019352094;
3798 c = -14.2295510326;
3799 d = 2.2104490688;
3800 e = 4.2969671500;
3801 } else if (fModCentralityClass == 0 && fCentralityMin == 0 && fCentralityMax == 4){ // 0-40% PbPb
3802 a = 16.8227412106;
3803 b = 5.8660502207;
3804 c = -12.0978551215;
3805 d = 2.1695068981;
3806 e = 3.5349621182;
3807 } else if (fModCentralityClass == 0 && fCentralityMin == 0 && fCentralityMax == 8){ // 0-80% PbPb
3808 a = 9.4675681080;
3809 b = 5.8114944205;
3810 c = -10.4901523616;
3811 d = 2.0607982712;
3812 e = 2.9262259130;
3813 } else if (fModCentralityClass == 0 && fCentralityMin == 4 && fCentralityMax == 8){ // 60-80% PbPb
3814 a = 2.5985551785;
3815 b = 5.4118895738;
3816 c = -8.2510958428;
3817 d = 2.2551249190;
3818 e = 3.0700919491;
3819 }
3820
3821 functionResultData = a*TMath::Power(mesonPt,-1*(b+c/(TMath::Power(mesonPt,d)+e)));
3822 }
3823
3824 }
3825
3826 Double_t weight = 1;
3827 if (PDGCode == 111 || PDGCode == 221){
3828 if (functionResultData != 0. && functionResultMC != 0. && isfinite(functionResultData) && isfinite(functionResultMC)){
3829 weight = functionResultData/functionResultMC;
4803eb1f 3830 if ( kCaseGen == 3){
3831 if (!(fDoReweightHistoMCPi0 && hReweightMCHistPi0!= 0x0 && PDGCode == 111)){
3832 weight = 1.;
3833 }
ae947965 3834 }
3835 if (!isfinite(functionResultData)) weight = 1.;
3836 if (!isfinite(weight)) weight = 1.;
3837 }
3838 } else if (PDGCode == 310 && functionResultMC != 0 && isfinite(functionResultMC)){
3839 weight = functionResultMC;
3840 }
3841
3842// if (fModCentralityClass == 0 && fCentralityMin == 4 && fCentralityMax == 6 && PDGCode == 111){
4803eb1f 3843// cout << period.Data() << "\t" << kCaseGen << "\t" <<fModCentralityClass<< "\t" <<fCentralityMin<< "\t" <<fCentralityMax << "\t" << mesonPt << "\t" <<mesonMass<< "\t"<<functionResultData << "\t"<< functionResultMC << "\t" << weight <<endl;
ae947965 3844// }
3845 return weight;
a280ac15 3846}
a280ac15 3847///________________________________________________________________________
3848AliConversionCuts* AliConversionCuts::GetStandardCuts2010PbPb(){
3849 //Create and return standard 2010 PbPb cuts
3850 AliConversionCuts *cuts=new AliConversionCuts("StandardCuts2010PbPb","StandardCuts2010PbPb");
430b9c94 3851 if(!cuts->InitializeCutsFromCutString("100000204209297002322000000")){
a280ac15 3852 cout<<"Warning: Initialization of Standardcuts2010PbPb failed"<<endl;}
3853 return cuts;
3854}
3855
3856///________________________________________________________________________
3857AliConversionCuts* AliConversionCuts::GetStandardCuts2010pp(){
3858 //Create and return standard 2010 PbPb cuts
3859 AliConversionCuts *cuts=new AliConversionCuts("StandardCuts2010pp","StandardCuts2010pp");
430b9c94 3860 if(!cuts->InitializeCutsFromCutString("000001100209366300380000000")){
a280ac15 3861 cout<<"Warning: Initialization of Standardcuts2010pp failed"<<endl;}
3862 return cuts;
3863}
ccfa8c0d 3864///________________________________________________________________________
3865void AliConversionCuts::GetCorrectEtaShiftFromPeriod(TString periodName){
3866
3867 if(periodName.CompareTo("LHC12g") == 0 || //pilot run 2012
3868 periodName.CompareTo("LHC13b") == 0 || //mainly minimum bias
3869 periodName.CompareTo("LHC13c") == 0 || //mainly minimum bias
3870 periodName.CompareTo("LHC13d") == 0 || //mainly triggered
3871 periodName.CompareTo("LHC13e") == 0 || //mainly triggered
3872 periodName.CompareTo("LHC13c3") == 0 || //MC Starlight, anchor LHC13d+e
3873 periodName.CompareTo("LHC13c2") == 0 || //MC Starlight, coherent J/Psi, UPC muon anchor LHC13d+e
3874 periodName.CompareTo("LHC13b4") == 0 || //MC Pythia 6 (Jet-Jet), anchor LHC13b
3875 periodName.CompareTo("LHC13b2_fix_1") == 0 || //MC DPMJET, anchr LHC13b+c
3876 periodName.CompareTo("LHC13b3") == 0 || //MC HIJING, weighted to number of events per run, anchor LHC13b
3877 periodName.CompareTo("LHC13b2") == 0 || // MC DPMJET, wrong energy, anchor LHC13b
3878 periodName.CompareTo("LHC13b2_plus") == 0 || // MC DPMJET, weighted to number event per run, anchor LHC13b
3879 periodName.CompareTo("LHC13c1_bis") == 0 || // MC AMPT fast generation, pT hardbin, anchor ?
3880 periodName.CompareTo("LHC13c1") == 0 || // MC AMPT fast generation, anchor ?
3881 periodName.CompareTo("LHC13b1") == 0 || // MC DPMJET, fragments, with fixed label 0, anchor LHC12g
3882 periodName.CompareTo("LHC12g4b_fix") == 0 || // MC DPMJET, with fixed label 0, anchor LHC12g
3883 periodName.CompareTo("LHC12g1_fix") == 0 || // MC ?, with fixed label 0, anchor LHC12g
3884 periodName.CompareTo("LHC12g4c") == 0 || // MC DPMJET, shifted vertex runs, anchor LHC12g
3885 periodName.CompareTo("LHC12h6") == 0 || // MC muon cocktail, anchor LHC12g
3886 periodName.CompareTo("LHC12g4b") == 0 || // MC DPMJET 3rd iteration, anchor LHC12g
3887 periodName.CompareTo("LHC12g4a") == 0 || // MC DPMJET improved, anchor LHC12g
3888 periodName.CompareTo("LHC12g4") == 0 || // MC DPMJET, anchor LHC12g
3889 periodName.CompareTo("LHC12g5") == 0 || // MC PHOJET, anchor LHC12g
3890 periodName.CompareTo("LHC12g2") == 0 || // MC Starlight background, anchor LHC12g
3891 periodName.CompareTo("LHC12g1") == 0 ) // MC ?, anchor LHC12g
3892 {
3893 printf(" Gamma Conversion Cuts %s :: pPb Run doing Eta Shift of %f \n\n",(GetCutNumber()).Data(),-0.465);
3894 SetEtaShift(-0.465);
3895 }
3896 else if(periodName.CompareTo("LHC13f") == 0 ||
3897 periodName.CompareTo("LHC13c6b") == 0 ||// MC Jpsi -> mumu, anchor LHC13f
3898 periodName.CompareTo("LHC13c5") == 0 || //MC Starlight, gamma gamma UPC muon, anchor LHC13f
3899 periodName.CompareTo("LHC13c4") == 0 )//MC Starlight, coherent JPsi, UPC muon, anchor LHC13f
3900 {
3901 printf(" Gamma Conversion Cuts %s :: Pbp Run doing Eta Shift of %f \n\n",(GetCutNumber()).Data(),0.465);
3902 SetEtaShift(+0.465);
3903 }
3904 else printf(" Gamma Conversion Cuts %s :: Automatic Eta Shift requested but Period is not known -> No Shift \n\n",(GetCutNumber()).Data());
3905}