]> git.uio.no Git - u/mrichter/AliRoot.git/blame - PWGCF/FEMTOSCOPY/macros/Train/ProtonFemtoscopy/TTC/ConfigFemtoAnalysis.C
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Merge branch 'feature-movesplit'
[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / ProtonFemtoscopy / TTC / ConfigFemtoAnalysis.C
CommitLineData
59c3fe5c 1/*********************************************************************
2 * *
3 * Configfemtoanalysis.C - configuration macro for the femtoscopic *
4 * analysis, meant as a QA process for two-particle effects *
5 * *
6 * Author: Adam Kisiel (Adam.Kisiel@cern.ch) *
7 * *
8 *********************************************************************/
9
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"
50#endif
51
52//________________________________________________________________________
53AliFemtoManager* ConfigFemtoAnalysis() {
54
55 double PionMass = 0.13956995;
56 double KaonMass = 0.493677;
57 double ProtonMass = 0.938272013;
58
59 double psi = TMath::Pi()/2.;
60 double psid = TMath::Pi()/6.;
61
62 int runepvzero[7] = {0, 0, 0, 0, 0, 0, 1};
63 double epvzerobins[7] = {-psi, -psi+psid, -psi+2*psid, -psi+3*psid, -psi+4*psid, -psi+5*psid, -psi+6*psid};
64
65 int runmults[10] = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
66 int multbins[11] = {0.001, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
67
68 int runch[3] = {1, 0, 1};
69 const char *chrgs[3] = { "PP", "APAP", "PAP" };
70
71 int runktdep = 1;
72 double ktrng[3] = {0.01, 1.0, 100.0};
73
74 int numOfMultBins = 10;
75 int numOfChTypes = 3;
76 int numOfkTbins = 2;
77 int numOfEPvzero = 7;
78
79 int runqinv = 1;
80 int runshlcms = 0;// 0:PRF(PAP), 1:LCMS(PP,APAP)
81
82 int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
83 int isrealdata = 1;
84
85 // int gammacut = 1;
86
87 double shqmax = 2.0;
88 int nbinssh = 200;
89
90 AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
91 Reader->SetFilterBit(7);
92 Reader->SetCentralityPreSelection(0.001, 510);
93
94 AliFemtoManager* Manager = new AliFemtoManager();
95 Manager->SetEventReader(Reader);
96
97 AliFemtoVertexMultAnalysis *anetaphitpc[10*3*2];
98 AliFemtoBasicEventCut *mecetaphitpc[10*3*2];
99 AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[100];
100 AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[100];
101
102 // AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[50];
103 // AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[50];
104 AliFemtoESDTrackCut *dtc1etaphitpc[50];
105 AliFemtoESDTrackCut *dtc2etaphitpc[50];
106 AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[50];
107 AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[50];
108 AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[50];
109 AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[50];
110 AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[50];
111 AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[50];
112 AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[50];
113 AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[50];
114 // AliFemtoPairCutAntiGamma *sqpcetaphitpcdiff[10*3];
115 // AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpcsame[10*3];
116 //AliFemtoPairCutAntiGamma *sqpcetaphitpc[10*3];
117 AliFemtoPairCutRadialDistance *sqpcetaphitpc[50];
118 // AliFemtoChi2CorrFctn *cchiqinvetaphitpc[20*2];
119 AliFemtoKTPairCut *ktpcuts[100*2];
120 // AliFemtoCorrFctnDirectYlm *cylmkttpc[50];
121 //AliFemtoCorrFctnDirectYlm *cylmetaphitpc[10*3];
122 AliFemtoQinvCorrFctn *cqinvkttpc[50*2];
123 AliFemtoQinvCorrFctn *cqinvtpc[50];
124 AliFemtoCorrFctnNonIdDR *ckstartpc[50];
125 AliFemtoCorrFctnNonIdDR *ckstarkttpc[100*2];
126 // AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[50];
127
128 // AliFemtoCorrFctn3DLCMSSym *cq3dlcmskttpc[20*2];
129 // AliFemtoCorrFctnTPCNcls *cqinvnclstpc[20];
130 // AliFemtoShareQualityCorrFctn *cqinvsqtpc[20*10];
131 // AliFemtoChi2CorrFctn *cqinvchi2tpc[20];
132 // AliFemtoTPCInnerCorrFctn *cqinvinnertpc[20*10];
133
134 // *** Third QA task - HBT analysis with all pair cuts off, TPC only ***
135 // *** Begin pion-pion (positive) analysis ***
136 int aniter = 0;
137
138 for (int imult = 0; imult < numOfMultBins; imult++) {
139 if (runmults[imult]) {
140
141 for (int ichg = 0; ichg < numOfChTypes; ichg++) {
142 if (runch[ichg]) {
143
144 for (int iepvzero = 0; iepvzero < numOfEPvzero; iepvzero++) {
145 if (runepvzero[iepvzero]) {
146
147 aniter = ichg * numOfMultBins + imult * numOfEPvzero + iepvzero;
148
149 cout << "aniter = " << aniter << endl;
150 // aniter = ichg * numOfMultBins + imult;
151
152 // if (ichg == 2)
153 // runshlcms = 0;
154 // else
155 // runshlcms = 1;
156
157 anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(8, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);
158 anetaphitpc[aniter]->SetNumEventsToMix(10);
159 anetaphitpc[aniter]->SetMinSizePartCollection(1);
160
161 mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
162 mecetaphitpc[aniter]->SetEventMult(0.001,100000);
163 mecetaphitpc[aniter]->SetVertZPos(-8,8);
164
165 if (iepvzero == 6)
166 mecetaphitpc[aniter]->SetEPVZERO(epvzerobins[0],epvzerobins[6]);
167 else
168 mecetaphitpc[aniter]->SetEPVZERO(epvzerobins[iepvzero],epvzerobins[iepvzero+1]);
169
170
171 // cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%iPsi%i", chrgs[ichg], imult, iepvzero));
172 // cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%iPsi%i", chrgs[ichg], imult, iepvzero));
173 // mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
174
175 // cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
176 // cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
177 // mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
178
179 dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
180 dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
181
182 if (ichg == 0) {
183 dtc1etaphitpc[aniter]->SetCharge(1.0);
184 dtc1etaphitpc[aniter]->SetPt(0.7,5.0);
185 }
186 else if (ichg == 1) {
187 dtc1etaphitpc[aniter]->SetCharge(-1.0);
188 dtc1etaphitpc[aniter]->SetPt(0.7,5.0);
189 }
190 else if (ichg == 2) {
191 dtc1etaphitpc[aniter]->SetCharge(-1.0);
192 dtc2etaphitpc[aniter]->SetCharge(1.0);
193 dtc1etaphitpc[aniter]->SetPt(0.7,5.0);
194 dtc2etaphitpc[aniter]->SetPt(0.7,5.0);
195 }
196
197 dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
198 dtc1etaphitpc[aniter]->SetMass(ProtonMass);
199 dtc1etaphitpc[aniter]->SetMostProbableProton();
200
201 if (ichg == 2) {
202 dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
203 dtc2etaphitpc[aniter]->SetMass(ProtonMass);
204 dtc2etaphitpc[aniter]->SetMostProbableProton();
205 }
206
207
208 // Track quality cuts
209
210 if (runtype == 0) {
211 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
212 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
213 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
214 dtc1etaphitpc[aniter]->SetminTPCncls(80);
215 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
216 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
217 // dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
218 dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
219 dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
220 // dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
221 dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15);
222 // dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
223 }
224 else if (runtype == 1) {
225 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
226 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
227 // dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit|AliESDtrack::kITSpureSA);
228 // dtc1etaphitpc[aniter]->SetminTPCncls(70);
229 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
230 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
231 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
232 // dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
233 // dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(6.0);
234 dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
235 dtc1etaphitpc[aniter]->SetMaxImpactZ(0.25);
236 // dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
237 }
238 else if (runtype == 2) {
239 //dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
240 dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
241 dtc1etaphitpc[aniter]->SetminTPCncls(80);
242 dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
243 dtc1etaphitpc[aniter]->SetLabel(kFALSE);
244 dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
245 dtc1etaphitpc[aniter]->SetMaxImpactXY(0.1); // 2.4 0.1
246 dtc1etaphitpc[aniter]->SetMaxImpactZ(2.0); // 2.0 0.1
247
248 if (ichg == 2) {
249 //dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
250 dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
251 dtc2etaphitpc[aniter]->SetminTPCncls(80);
252 dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
253 dtc2etaphitpc[aniter]->SetLabel(kFALSE);
254 dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
255 dtc2etaphitpc[aniter]->SetMaxImpactXY(0.1); // 2.4 0.1
256 dtc2etaphitpc[aniter]->SetMaxImpactZ(2.0); // 2.0 0.1
257 }
258
259 }
260
261 // cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%iPsi%i", chrgs[ichg], imult, iepvzero),ProtonMass);
262 // cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%iPsi%i", chrgs[ichg], imult, iepvzero),ProtonMass);
263 // dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
264
265 // cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass1%stpcM%iPsi%i", chrgs[ichg], imult, iepvzero),2);//0-pion,1-kaon,2-proton
266 // cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%iPsi%i", chrgs[ichg], imult , iepvzero),2);
267 // dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);
268
269 // if (ichg == 2){
270 // cutPass2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass2%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
271 // cutFail2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail2%stpcM%i", chrgs[ichg], imult),2);
272 // dtc2etaphitpc[aniter]->AddCutMonitor(cutPass2PIDetaphitpc[aniter], cutFail2PIDetaphitpc[aniter]);
273 // }
274
275 // sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
276 sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
277
278 if (runtype == 0) {
279 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
280 sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
281 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
282 // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
283 // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
284 // sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
285 //sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
286 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
287 }
288 else if (runtype == 1) {
289 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
290 sqpcetaphitpc[aniter]->SetShareFractionMax(1.05);
291 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
292 // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
293 // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
294 // sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(5.0);
295 //sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.03);
296 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
297 }
298 else if (runtype == 2) {
299 //sqpcetaphitpc[aniter]->SetUseAOD(kTRUE);
300
301 sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
302 sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
303 sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
304
305
306 // if (gammacut == 0) {
307 //sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
308 //sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
309 //}
310 //else if (gammacut == 1) {
311 //sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
312 //sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
313 //}
314
315 sqpcetaphitpc[aniter]->SetMagneticFieldSign(-1); // field1 -1, field3 +1
316 sqpcetaphitpc[aniter]->SetMinimumRadius(0.8); // biggest inefficiency for R=1.1 m (checked on small sample)
317 sqpcetaphitpc[aniter]->SetPhiStarMin(1.2); // biggest inefficiency for R=1.1 m (checked on small sample)
318 // sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.045); // 0.012 - pions, 0.017 - kaons, 0.018
319 // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.01); // 0.017 - pions, 0.015 - kaons
320
321
322
323 }
324
325 anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
326
327 if (ichg == 2) {
328 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
329 anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
330 }
331 else {
332 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
333 anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
334 }
335
336 anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
337
338
339 // if (ichg == 2) {
340 // ckstartpc[aniter] = new AliFemtoCorrFctnNonIdDR(Form("ckstar%stpcM%iPsi%iTTC08", chrgs[ichg], imult, iepvzero),nbinssh,0.0,shqmax);
341 // anetaphitpc[aniter]->AddCorrFctn(ckstartpc[aniter]);
342
343 // }
344 // else {
345
346 cqinvtpc[aniter] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%iPsi%iTTC", chrgs[ichg], imult, iepvzero),2*nbinssh,0.0,2*shqmax);
347 cqinvtpc[aniter]->CalculateDetaDphis(kTRUE,1.2);
348 anetaphitpc[aniter]->AddCorrFctn(cqinvtpc[aniter]);
349
350 // }
351
352 // cylmkttpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),2,nbinssh, 0.0,shqmax,runshlcms);
353 // anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[aniter]);
354
355
356 // cqinvinnertpc[aniter] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM", chrgs[ichg], imult),nbinssh,0.0,shqmax);
357 // anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[aniter]);
358
359
360 if (runktdep) {
361 int ktm;
362 for (int ikt=0; ikt<numOfkTbins; ikt++) {
363
364 ktm = aniter * numOfkTbins + ikt;
365 ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
366
367
368 //cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,nbinssh,0.0,shqmax,0);
369 //cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
370 //anetaphitpc[ktm]->AddCorrFctn(cylmkttpc[ktm]);
371
372
373
374 // if (ichg == 2) {
375 // ckstarkttpc[ktm] = new AliFemtoCorrFctnNonIdDR(Form("ckstar%stpcM%iPsi%ikT%iTTC08", chrgs[ichg], imult, iepvzero, ikt),nbinssh,0.0,shqmax);
376 // ckstarkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
377 // anetaphitpc[aniter]->AddCorrFctn(ckstarkttpc[ktm]);
378
379 // }
380 // else {
381 cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%iPsi%ikT%iTTC", chrgs[ichg], imult, iepvzero, ikt),nbinssh,0.0,shqmax);
382 cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
383 cqinvkttpc[ktm]->CalculateDetaDphis(kTRUE,1.2);
384 anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
385
386 // }
387
388 // cqinvsqtpc[ktm] = new AliFemtoShareQualityCorrFctn(Form("cqinvsq%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
389 // cqinvsqtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
390 // anetaphitpc[aniter]->AddCorrFctn(cqinvsqtpc[ktm]);
391
392 // cqinvinnertpc[ktm] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
393 // cqinvinnertpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
394 // cqinvinnertpc[ktm]->SetRadius(1.2);
395 // anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[ktm]);
396
397 // if (run3d) {
398 // 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);
399 // cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
400 // anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
401 // }
402 }
403 }
404
405 // cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),24, 24);
406 // anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
407
408 Manager->AddAnalysis(anetaphitpc[aniter]);
409 }
410 }
411 }
412 }
413 }
414 }
415// *** End pion-pion (positive) analysis
416
417 return Manager;
418}