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