1 /*********************************************************************
3 * Configfemtoanalysis.C - configuration macro for the femtoscopic *
4 * analysis, meant as a QA process for two-particle effects *
6 * Author: Adam Kisiel (Adam.Kisiel@cern.ch) *
8 *********************************************************************/
10 #if !defined(__CINT__) || defined(__MAKECINT_)
11 #include "AliFemtoManager.h"
12 #include "AliFemtoEventReaderESDChain.h"
13 #include "AliFemtoEventReaderESDChainKine.h"
14 #include "AliFemtoEventReaderAODChain.h"
15 #include "AliFemtoSimpleAnalysis.h"
16 #include "AliFemtoBasicEventCut.h"
17 #include "AliFemtoESDTrackCut.h"
18 #include "AliFemtoCorrFctn.h"
19 #include "AliFemtoCutMonitorParticleYPt.h"
20 #include "AliFemtoCutMonitorParticleVertPos.h"
21 #include "AliFemtoCutMonitorParticleMomRes.h"
22 #include "AliFemtoCutMonitorParticlePID.h"
23 #include "AliFemtoCutMonitorEventMult.h"
24 #include "AliFemtoCutMonitorEventVertex.h"
25 #include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
26 #include "AliFemtoPairCutAntiGamma.h"
27 #include "AliFemtoPairCutRadialDistance.h"
28 #include "AliFemtoQinvCorrFctn.h"
29 #include "AliFemtoCorrFctnNonIdDR.h"
30 #include "AliFemtoShareQualityCorrFctn.h"
31 #include "AliFemtoTPCInnerCorrFctn.h"
32 #include "AliFemtoVertexMultAnalysis.h"
33 #include "AliFemtoCorrFctn3DSpherical.h"
34 #include "AliFemtoChi2CorrFctn.h"
35 #include "AliFemtoCorrFctnTPCNcls.h"
36 #include "AliFemtoBPLCMS3DCorrFctn.h"
37 #include "AliFemtoCorrFctn3DLCMSSym.h"
38 #include "AliFemtoModelBPLCMSCorrFctn.h"
39 #include "AliFemtoModelCorrFctn3DSpherical.h"
40 #include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
41 #include "AliFemtoModelGausRinvFreezeOutGenerator.h"
42 #include "AliFemtoModelManager.h"
43 #include "AliFemtoModelWeightGeneratorBasic.h"
44 #include "AliFemtoModelWeightGeneratorLednicky.h"
45 #include "AliFemtoCorrFctnDirectYlm.h"
46 #include "AliFemtoModelCorrFctnDirectYlm.h"
47 #include "AliFemtoModelCorrFctnSource.h"
48 #include "AliFemtoCutMonitorParticlePtPDG.h"
49 #include "AliFemtoKTPairCut.h"
52 //________________________________________________________________________
53 AliFemtoManager* ConfigFemtoAnalysis() {
55 double PionMass = 0.13956995;
56 double KaonMass = 0.493677;
57 double ProtonMass = 0.938272013;
60 int runmults[10] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
61 int multbins[11] = {0.001, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
63 int runch[3] = {1, 1, 1};
64 const char *chrgs[3] = { "PP", "APAP", "PAP" };
68 double ktrng[3] = {0.01, 0.7, 100.0};
69 // double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};
71 int numOfMultBins = 10;
76 int runshlcms = 1;// 0:PRF(PAP), 1:LCMS(PP,APAP)
78 int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
86 AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
87 Reader->SetFilterBit(7);
88 Reader->SetCentralityPreSelection(0.001, 910);
90 AliFemtoManager* Manager = new AliFemtoManager();
91 Manager->SetEventReader(Reader);
93 AliFemtoVertexMultAnalysis *anetaphitpc[10*3];
94 AliFemtoBasicEventCut *mecetaphitpc[10*3];
95 AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[10*3];
96 AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[10*3];
97 AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[10*3];
98 AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[10*3];
99 AliFemtoESDTrackCut *dtc1etaphitpc[10*3];
100 AliFemtoESDTrackCut *dtc2etaphitpc[10*3];
101 AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[10*3];
102 AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[10*3];
103 AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[10*3];
104 AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[10*3];
105 AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[10*3];
106 AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[10*3];
107 AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[10*3];
108 AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[10*3];
109 // AliFemtoPairCutAntiGamma *sqpcetaphitpcdiff[10*3];
110 // AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpcsame[10*3];
111 AliFemtoPairCutAntiGamma *sqpcetaphitpc[10*3];
112 // AliFemtoPairCutRadialDistance *sqpcetaphitpc[10*3];
113 // AliFemtoChi2CorrFctn *cchiqinvetaphitpc[20*2];
114 AliFemtoKTPairCut *ktpcuts[10*3*2];
115 AliFemtoCorrFctnDirectYlm *cylmkttpc[10*3];
116 //AliFemtoCorrFctnDirectYlm *cylmetaphitpc[10*3];
117 AliFemtoQinvCorrFctn *cqinvkttpc[10*3*2];
118 AliFemtoQinvCorrFctn *cqinvtpc[10*3];
119 AliFemtoCorrFctnNonIdDR *ckstartpc[10*3];
120 AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[10*3];
122 // AliFemtoCorrFctn3DLCMSSym *cq3dlcmskttpc[20*2];
123 // AliFemtoCorrFctnTPCNcls *cqinvnclstpc[20];
124 // AliFemtoShareQualityCorrFctn *cqinvsqtpc[20*10];
125 // AliFemtoChi2CorrFctn *cqinvchi2tpc[20];
126 // AliFemtoTPCInnerCorrFctn *cqinvinnertpc[20*10];
128 // *** Third QA task - HBT analysis with all pair cuts off, TPC only ***
129 // *** Begin pion-pion (positive) analysis ***
132 for (int imult = 0; imult < numOfMultBins; imult++) {
133 if (runmults[imult]) {
135 for (int ichg = 0; ichg < numOfChTypes; ichg++) {
138 aniter = ichg * numOfMultBins + imult;
139 anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(8, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);
140 anetaphitpc[aniter]->SetNumEventsToMix(10);
141 anetaphitpc[aniter]->SetMinSizePartCollection(1);
143 mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
144 mecetaphitpc[aniter]->SetEventMult(0.001,100000);
145 mecetaphitpc[aniter]->SetVertZPos(-8,8);
148 mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
150 cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
151 cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
152 mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
154 cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
155 cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
156 mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
158 dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
159 dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
162 dtc1etaphitpc[aniter]->SetCharge(1.0);
164 dtc1etaphitpc[aniter]->SetCharge(-1.0);
165 else if (ichg == 2) {
166 dtc1etaphitpc[aniter]->SetCharge(-1.0);
167 dtc2etaphitpc[aniter]->SetCharge(1.0);
170 dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
171 dtc1etaphitpc[aniter]->SetMass(ProtonMass);
172 dtc1etaphitpc[aniter]->SetMostProbableProton();
175 dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
176 dtc2etaphitpc[aniter]->SetMass(ProtonMass);
177 dtc2etaphitpc[aniter]->SetMostProbableProton();
180 // Track quality cuts
183 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
184 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
185 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
186 dtc1etaphitpc[aniter]->SetminTPCncls(80);
187 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
188 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
189 // dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
190 dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
191 dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
192 // dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
193 dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15);
194 // dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
196 else if (runtype == 1) {
197 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
198 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
199 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit|AliESDtrack::kITSpureSA);
200 // dtc1etaphitpc[aniter]->SetminTPCncls(70);
201 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
202 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
203 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
204 // dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
205 // dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(6.0);
206 dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
207 dtc1etaphitpc[aniter]->SetMaxImpactZ(0.25);
208 // dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
210 else if (runtype == 2) {
211 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
212 dtc1etaphitpc[aniter]->SetminTPCncls(80);
213 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
214 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
215 dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
216 dtc1etaphitpc[aniter]->SetMaxImpactXY(2.4);
217 dtc1etaphitpc[aniter]->SetMaxImpactZ(3.0);
220 dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
221 dtc2etaphitpc[aniter]->SetminTPCncls(80);
222 dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
223 dtc2etaphitpc[aniter]->SetLabel(kFALSE);
224 dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
225 dtc2etaphitpc[aniter]->SetMaxImpactXY(2.4);
226 dtc2etaphitpc[aniter]->SetMaxImpactZ(3.0);
231 cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%i", chrgs[ichg], imult),ProtonMass);
232 cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%i", chrgs[ichg], imult),ProtonMass);
233 dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
235 cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass1%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
236 cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%i", chrgs[ichg], imult),2);
237 dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);
240 cutPass2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass2%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
241 cutFail2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail2%stpcM%i", chrgs[ichg], imult),2);
242 dtc2etaphitpc[aniter]->AddCutMonitor(cutPass2PIDetaphitpc[aniter], cutFail2PIDetaphitpc[aniter]);
245 sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
248 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
249 sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
250 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
251 // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
252 // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
253 // sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
254 //sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
255 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
257 else if (runtype == 1) {
258 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
259 sqpcetaphitpc[aniter]->SetShareFractionMax(1.05);
260 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
261 // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
262 // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
263 // sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(5.0);
264 //sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.03);
265 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
267 else if (runtype == 2) {
268 sqpcetaphitpc[aniter]->SetUseAOD(kTRUE);
269 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
270 sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
271 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
274 sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
275 sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
277 else if (gammacut == 1) {
278 sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
279 sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
282 //phi-star cut - values from Johana
283 // sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.012);
284 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.017);
288 anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
291 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
292 anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
295 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
296 anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
299 anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
303 ckstartpc[aniter] = new AliFemtoCorrFctnNonIdDR(Form("ckstar%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
304 anetaphitpc[aniter]->AddCorrFctn(ckstartpc[aniter]);
307 cqinvtpc[aniter] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
308 anetaphitpc[aniter]->AddCorrFctn(cqinvtpc[aniter]);
311 cylmkttpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),2,nbinssh, 0.0,shqmax,runshlcms);
312 anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[aniter]);
317 for (int ikt=0; ikt<numOfkTbins; ikt++) {
319 ktm = aniter * numOfkTbins + ikt;
320 ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
323 cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,nbinssh,0.0,shqmax,0);
324 cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
325 anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[ktm]);
327 cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
328 cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
329 anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
331 // cqinvsqtpc[ktm] = new AliFemtoShareQualityCorrFctn(Form("cqinvsq%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
332 // cqinvsqtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
333 // anetaphitpc[aniter]->AddCorrFctn(cqinvsqtpc[ktm]);
335 // cqinvinnertpc[ktm] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
336 // cqinvinnertpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
337 // cqinvinnertpc[ktm]->SetRadius(1.2);
338 // anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[ktm]);
341 // cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,(imult>3)?((imult>6)?((imult>7)?0.6:0.4):0.25):0.15);
342 // cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
343 // anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
348 // cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),24, 24);
349 // anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
351 Manager->AddAnalysis(anetaphitpc[aniter]);
356 // *** End pion-pion (positive) analysis