PWGCF/FEMTOSCOPY/macros/Train/ProtonLambdaFemto/SemiCentral/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 "AliFemtoSimpleAnalysis.h"
  15 #include "AliFemtoBasicEventCut.h"
  16 #include "AliFemtoESDTrackCut.h"
  17 #include "AliFemtoCorrFctn.h"
  18 #include "AliFemtoCutMonitorParticleYPt.h"
  19 #include "AliFemtoCutMonitorParticleVertPos.h"
  20 #include "AliFemtoCutMonitorParticleMomRes.h"
  21 #include "AliFemtoCutMonitorParticlePID.h"
  22 #include "AliFemtoCutMonitorEventMult.h"
  23 #include "AliFemtoCutMonitorEventVertex.h"
  24 #include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
  25 #include "AliFemtoPairCutAntiGamma.h"
  26 #include "AliFemtoQinvCorrFctn.h"
  27 #include "AliFemtoShareQualityCorrFctn.h"
  28 #include "AliFemtoTPCInnerCorrFctn.h"
  29 #include "AliFemtoVertexMultAnalysis.h"
  30 #include "AliFemtoCorrFctn3DSpherical.h"
  31 #include "AliFemtoChi2CorrFctn.h"
  32 #include "AliFemtoCorrFctnTPCNcls.h"
  33 #include "AliFemtoModelBPLCMSCorrFctn.h"
  34 #include "AliFemtoModelCorrFctn3DSpherical.h"
  35 #include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
  36 #include "AliFemtoModelGausRinvFreezeOutGenerator.h"
  37 #include "AliFemtoModelManager.h"
  38 #include "AliFemtoModelWeightGeneratorBasic.h"
  39 #include "AliFemtoModelWeightGeneratorLednicky.h"
  40 #include "AliFemtoCorrFctnDirectYlm.h"
  41 #include "AliFemtoModelCorrFctnDirectYlm.h"
  42 #include "AliFemtoModelCorrFctnSource.h"
  43 #include "AliFemtoCutMonitorParticlePtPDG.h"
  44 #include "AliFemtoKTPairCut.h"
  45 #include "AliFemtoEnumeration.h"
  46 #include "AliFemtoV0PairCut.h"
  47 #include "AliFemtoV0TrackPairCut.h"
  48 #include "AliFemtoV0TrackCut.h"
  49 #include "AliFemtoEventReaderAODChain.h"
  50 #include "AliFemtoCorrFctnNonIdDR.h"
  51 #include "AliFemtoCutMonitorCollections.h"
  52 #include "AliFemtoPairCutRadialDistance.h"
  53 #include "AliFemtoAvgSepCorrFctn.h"
  54 #endif
  55
  56 //________________________________________________________________________
  57 AliFemtoManager* ConfigFemtoAnalysis() {
  58
  59   double PionMass = 0.13956995;
  60   double KaonMass = 0.493677;
  61   double ProtonMass = 0.938272013;
  62   double LambdaMass = 1.115683;
  63
  64   double psi = TMath::Pi()/2.;
  65   double psid = TMath::Pi()/6.;
  66
  67   // Switches for QA analyses
  68   int runmults[10] = {0, 0, 1, 1, 1, 1, 0, 0, 0, 0};
  69   int multbins[11] = {0.00001, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
  70   int runch[10] = {1, 1, 1, 1, 1, 1, 1, 0, 0, 0};
  71   //int runch[10] = {1, 1, 1, 0, 0, 0, 0, 0, 0, 0};
  72   const char *chrgs[10] = { "V0LL", "V0ALAL", "V0LAL", "V0PL", "V0APL", "V0PAL", "V0APAL","PP","PAP","APAP" };
  73
  74   int runktdep = 0;
  75   double ktrng[3] = {0.01, 0.7, 100};
  76
  77   int numOfMultBins = 10;
  78   int numOfChTypes = 10;
  79   int numOfkTBins = 2;
  80
  81   int runshlcms = 0;
  82
  83   int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
  84
  85   double shqmax;
  86   int nbinssh = 100;
  87
  88   shqmax = 1.0;
  89
  90   AliFemtoEventReaderAODChain* Reader = new AliFemtoEventReaderAODChain();
  91   Reader->SetFilterBit(7);
  92   Reader->SetCentralityPreSelection(0.00001, 710);
  93   Reader->SetReadV0(1); //Read V0
  94   Reader->SetEPVZERO(kTRUE);
  95   Reader->SetCentralityFlattening(kTRUE);
  96
  97   // Reader->SetMagneticFieldSign(-1.0); //-1->field1, 1->field3
  98
  99   AliFemtoManager* Manager=new AliFemtoManager();
 100   Manager->SetEventReader(Reader);
 101
 102
 103   AliFemtoVertexMultAnalysis    *anetaphitpc[320];
 104   AliFemtoBasicEventCut         *mecetaphitpc[320];
 105   AliFemtoCutMonitorEventMult   *cutPassEvMetaphitpc[320];
 106   AliFemtoCutMonitorEventMult   *cutFailEvMetaphitpc[320];
 107   AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[320];
 108   AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[320];
 109   AliFemtoCutMonitorCollections   *cutPassColletaphitpc[320];
 110   AliFemtoCutMonitorCollections   *cutFailColletaphitpc[320];
 111   AliFemtoV0TrackCut           *dtc1etaphitpc[320];
 112   AliFemtoV0TrackCut           *dtc2etaphitpc[320];
 113   AliFemtoESDTrackCut           *dtc3etaphitpc[320];
 114   AliFemtoESDTrackCut           *dtc4etaphitpc[320];
 115   AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[320];
 116   AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[320];
 117   AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[320];
 118   AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[320];
 119   AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[320];
 120   AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[320];
 121   AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[320];
 122   AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[320];
 123   AliFemtoCutMonitorV0          *cutPass1V0[320];
 124   AliFemtoCutMonitorV0          *cutFail1V0[320];
 125   AliFemtoCutMonitorV0          *cutPass2V0[320];
 126   AliFemtoCutMonitorV0          *cutFail2V0[320];
 127   AliFemtoV0PairCut             *sqp1cetaphitpc[320];
 128   AliFemtoV0TrackPairCut        *sqp2cetaphitpc[320];
 129   AliFemtoV0TrackPairCut        *sqp3cetaphitpc[320];
 130   AliFemtoV0TrackPairCut        *sqp4cetaphitpc[320];
 131   //AliFemtoPairCutAntiGamma      *sqp3cetaphitpc[320];
 132   AliFemtoPairCutRadialDistance *sqp5cetaphitpc[320];
 133   AliFemtoCorrFctnDirectYlm     *cylmetaphitpc[320];
 134   AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[320*6];
 135   AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[320];
 136   AliFemtoKTPairCut             *ktpcuts[320*6];
 137   AliFemtoCorrFctnDirectYlm     *cylmkttpc[320*6];
 138   AliFemtoQinvCorrFctn          *cqinvkttpc[320*6];
 139   AliFemtoAvgSepCorrFctn        *avgsepcorr[320*6];
 140   //AliFemtoBPLCMS3DCorrFctn      *cq3dlcmskttpc[320*6];
 141   AliFemtoCorrFctnNonIdDR       *cnonidtpc[320*6];
 142   AliFemtoCorrFctnTPCNcls       *cqinvnclstpc[320];
 143   AliFemtoShareQualityCorrFctn  *cqinvsqtpc[320];
 144   AliFemtoTPCInnerCorrFctn      *cqinvtitpc[320];
 145   //AliFemtoChi2CorrFctn          *cqinvchi2tpc[320];
 146   AliFemtoQinvCorrFctn          *cqinvtpc[100*3];
 147
 148
 149   // *** Third QA task - HBT analysis with all pair cuts off, TPC only ***
 150   // *** Begin pion-pion (positive) analysis ***
 151   int aniter = 0;
 152
 153   for (int imult=0; imult<numOfMultBins; imult++) {
 154     if (runmults[imult]) {
 155       for (int ichg=0; ichg<numOfChTypes; ichg++) {
 156         if (runch[ichg]) {
 157           aniter = imult * numOfChTypes + ichg;
 158
 159           if(ichg==3 || ichg==4 || ichg==5 || ichg==6)
 160             runshlcms = 0;
 161           else
 162             runshlcms = 1;
 163
 164           anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(8, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);
 165           anetaphitpc[aniter]->SetNumEventsToMix(5);
 166           anetaphitpc[aniter]->SetMinSizePartCollection(1);
 167           anetaphitpc[aniter]->SetVerboseMode(kFALSE);
 168
 169           mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
 170           mecetaphitpc[aniter]->SetEventMult(0,100000);
 171           mecetaphitpc[aniter]->SetVertZPos(-8,8);
 172
 173           // cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
 174           // cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
 175           // mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
 176
 177           // cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
 178           // cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
 179           // mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
 180
 181           // cutPassColletaphitpc[aniter] = new AliFemtoCutMonitorCollections(Form("cutPass%stpcM%i", chrgs[ichg], imult));
 182           // cutFailColletaphitpc[aniter] = new AliFemtoCutMonitorCollections(Form("cutFail%stpcM%i", chrgs[ichg], imult));
 183           // mecetaphitpc[aniter]->AddCutMonitor(cutPassColletaphitpc[aniter], cutFailColletaphitpc[aniter]);
 184
 185           //V0 first particle cut -> Lambda ichg 0, 2, 6
 186           dtc1etaphitpc[aniter] = new AliFemtoV0TrackCut();
 187           dtc1etaphitpc[aniter]->SetMass(LambdaMass);
 188           dtc1etaphitpc[aniter]->SetEta(0.8); //0.8
 189           dtc1etaphitpc[aniter]->SetPt(0.5,5.0); //0.4,100
 190           dtc1etaphitpc[aniter]->SetEtaDaughters(0.8); //0.8
 191           dtc1etaphitpc[aniter]->SetPtPosDaughter(0.5,4.0); //0.5
 192           dtc1etaphitpc[aniter]->SetPtNegDaughter(0.16,4.0); //0.16
 193           dtc1etaphitpc[aniter]->SetTPCnclsDaughters(80); //80
 194           dtc1etaphitpc[aniter]->SetNdofDaughters(4.0); //4.0
 195           dtc1etaphitpc[aniter]->SetStatusDaughters(AliESDtrack::kTPCrefit/* | AliESDtrack::kITSrefit*/);
 196           dtc1etaphitpc[aniter]->SetOnFlyStatus(kFALSE);
 197           dtc1etaphitpc[aniter]->SetParticleType(0); //0-Lambda
 198           dtc1etaphitpc[aniter]->SetMinDaughtersToPrimVertex(0.1, 0.3);
 199           dtc1etaphitpc[aniter]->SetMaxDcaV0Daughters(0.4); //1.5 Jai, 0.6 //0.4
 200           dtc1etaphitpc[aniter]->SetMaxDcaV0(0.5); //5.0
 201           dtc1etaphitpc[aniter]->SetMaxV0DecayLength(60.0);
 202           dtc1etaphitpc[aniter]->SetMaxCosPointingAngle(0.9993); //0.99 - Jai //0.998
 203           dtc1etaphitpc[aniter]->SetInvariantMassLambda(LambdaMass-0.0034,LambdaMass+0.0034);
 204
 205           //V0 second particle cut -> AntiLambda ichg 1, 3, 4, 5
 206           dtc2etaphitpc[aniter] = new AliFemtoV0TrackCut();
 207           dtc2etaphitpc[aniter]->SetMass(LambdaMass);
 208           dtc2etaphitpc[aniter]->SetEta(0.8);
 209           dtc2etaphitpc[aniter]->SetPt(0.5,5.0);
 210           dtc2etaphitpc[aniter]->SetEtaDaughters(0.8);
 211           dtc2etaphitpc[aniter]->SetPtPosDaughter(0.16,4.0);
 212           dtc2etaphitpc[aniter]->SetPtNegDaughter(0.3,4.0);
 213           dtc2etaphitpc[aniter]->SetTPCnclsDaughters(80);
 214           dtc2etaphitpc[aniter]->SetNdofDaughters(4.0); //4.0
 215           dtc2etaphitpc[aniter]->SetStatusDaughters(AliESDtrack::kTPCrefit/* | AliESDtrack::kITSrefit*/);
 216           dtc2etaphitpc[aniter]->SetOnFlyStatus(kFALSE); //kTRUE
 217           dtc2etaphitpc[aniter]->SetParticleType(1);
 218           dtc2etaphitpc[aniter]->SetMaxDcaV0Daughters(0.4); //1.5 Jai, 0.6
 219           dtc2etaphitpc[aniter]->SetMaxDcaV0(0.5);
 220           dtc2etaphitpc[aniter]->SetMinDaughtersToPrimVertex(0.3, 0.1);
 221           dtc2etaphitpc[aniter]->SetMaxCosPointingAngle(0.9993); //0.99 - Jai
 222           dtc2etaphitpc[aniter]->SetMaxV0DecayLength(60.0);
 223           dtc2etaphitpc[aniter]->SetInvariantMassLambda(LambdaMass-0.0034,LambdaMass+0.0034);
 224
 225
 226           //ESD first particle cut -> Proton 3, 5; AntiProton 4, 6, 7, 8
 227           dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();
 228           dtc3etaphitpc[aniter]->SetMostProbableProton();
 229           dtc3etaphitpc[aniter]->SetMass(ProtonMass);
 230           dtc3etaphitpc[aniter]->SetCharge(1.0);
 231           dtc3etaphitpc[aniter]->SetPt(0.7,4.0);
 232           dtc3etaphitpc[aniter]->SetEta(-0.8,0.8);
 233
 234           // Track quality cuts
 235           dtc3etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
 236           dtc3etaphitpc[aniter]->SetminTPCncls(80);
 237           dtc3etaphitpc[aniter]->SetRemoveKinks(kTRUE);
 238           dtc3etaphitpc[aniter]->SetLabel(kFALSE);
 239           dtc3etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
 240           dtc3etaphitpc[aniter]->SetMaxImpactXY(0.1);
 241           dtc3etaphitpc[aniter]->SetMaxImpactZ(3.2);
 242           dtc3etaphitpc[aniter]->SetNsigma(3.0);
 243           dtc3etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
 244
 245           //ESD first particle cut -> Proton 3, 5; AntiProton 4, 6, 8, 9
 246           dtc4etaphitpc[aniter] = new AliFemtoESDTrackCut();
 247           dtc4etaphitpc[aniter]->SetMostProbableProton();
 248           dtc4etaphitpc[aniter]->SetMass(ProtonMass);
 249           dtc4etaphitpc[aniter]->SetCharge(-1.0);
 250           dtc4etaphitpc[aniter]->SetPt(0.7,5.0);
 251           dtc4etaphitpc[aniter]->SetEta(-0.8,0.8);
 252
 253           // Track quality cuts
 254           dtc4etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
 255           dtc4etaphitpc[aniter]->SetminTPCncls(80);
 256           dtc4etaphitpc[aniter]->SetRemoveKinks(kTRUE);
 257           dtc4etaphitpc[aniter]->SetLabel(kFALSE);
 258           dtc4etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
 259           dtc4etaphitpc[aniter]->SetMaxImpactXY(0.1);
 260           dtc4etaphitpc[aniter]->SetMaxImpactZ(3.2);
 261           dtc4etaphitpc[aniter]->SetNsigma(3.0);
 262           dtc4etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
 263
 264           // // //V0 monitor
 265           // cutPass1V0[aniter] = new AliFemtoCutMonitorV0(Form("cutPass1%stpcM%i", chrgs[ichg], imult));
 266           // cutFail1V0[aniter] = new AliFemtoCutMonitorV0(Form("cutFail1%stpcM%i", chrgs[ichg], imult));
 267           // dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1V0[aniter], cutFail1V0[aniter]);
 268
 269           // cutPass2V0[aniter] = new AliFemtoCutMonitorV0(Form("cutPass2%stpcM%i", chrgs[ichg], imult));
 270           // cutFail2V0[aniter] = new AliFemtoCutMonitorV0(Form("cutFail2%stpcM%i", chrgs[ichg], imult));
 271           // dtc2etaphitpc[aniter]->AddCutMonitor(cutPass2V0[aniter], cutFail2V0[aniter]);
 272
 273           // cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%i", chrgs[ichg], imult), ProtonMass);
 274           // cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%i", chrgs[ichg], imult), ProtonMass);
 275           // //if(ichg==7 || ichg==8 || ichg==9)
 276           // //dtc3etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
 277
 278           // cutPass2YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass2%stpcM%i", chrgs[ichg], imult), ProtonMass);
 279           // cutFail2YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail2%stpcM%i", chrgs[ichg], imult), ProtonMass);
 280           // //if(ichg==7 || ichg==8 || ichg==9)
 281           // //dtc4etaphitpc[aniter]->AddCutMonitor(cutPass2YPtetaphitpc[aniter], cutFail2YPtetaphitpc[aniter]);
 282
 283           // cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("trackCutPass1%stpcM%i", chrgs[ichg], imult),2);
 284           // cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("trackCutFail1%stpcM%i", chrgs[ichg], imult),2);
 285           // dtc3etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);
 286
 287           // cutPass2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("trackCutPass2%stpcM%i", chrgs[ichg], imult),2);
 288           // cutFail2PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("trackCutFail2%stpcM%i", chrgs[ichg], imult),2);
 289           // // //if(ichg == 4 || ichg == 5 || ichg == 6 || ichg == 7 || ichg == 8 || ichg == 9)
 290           // dtc4etaphitpc[aniter]->AddCutMonitor(cutPass2PIDetaphitpc[aniter], cutFail2PIDetaphitpc[aniter]);
 291
 292
 293           sqp1cetaphitpc[aniter] = new AliFemtoV0PairCut();
 294           sqp1cetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
 295           sqp1cetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.00001);
 296           sqp1cetaphitpc[aniter]->SetTPCExitSepMinimum(-1.);
 297           sqp1cetaphitpc[aniter]->SetMinAvgSeparation(0,3);
 298           sqp1cetaphitpc[aniter]->SetMinAvgSeparation(1,0);
 299           sqp1cetaphitpc[aniter]->SetMinAvgSeparation(2,0);
 300           sqp1cetaphitpc[aniter]->SetMinAvgSeparation(3,3);
 301
 302           sqp2cetaphitpc[aniter] = new AliFemtoV0TrackPairCut(); //lambda-proton
 303           sqp2cetaphitpc[aniter]->SetShareQualityMax(1.0); //between V0 daughter and track
 304           sqp2cetaphitpc[aniter]->SetShareFractionMax(0.05);
 305           sqp2cetaphitpc[aniter]->SetTPCOnly(kTRUE);
 306           sqp2cetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
 307           sqp2cetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.00001);
 308           sqp2cetaphitpc[aniter]->SetTPCExitSepMinimum(-1.);
 309           sqp2cetaphitpc[aniter]->SetKstarCut(0.04,AliFemtoV0TrackPairCut::kLambda,AliFemtoV0TrackPairCut::kProton); //0 - lambda, 2 - proton
 310           sqp2cetaphitpc[aniter]->SetMinAvgSeparation(0,11); //0 - track-pos, 1 - track-neg
 311           sqp2cetaphitpc[aniter]->SetMinAvgSeparation(1,0);
 312
 313           sqp3cetaphitpc[aniter] = new AliFemtoV0TrackPairCut(); //antilambda-antiproton
 314           sqp3cetaphitpc[aniter]->SetShareQualityMax(1.0); //between V0 daughter and track
 315           sqp3cetaphitpc[aniter]->SetShareFractionMax(0.05);
 316           sqp3cetaphitpc[aniter]->SetTPCOnly(kTRUE);
 317           sqp3cetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
 318           sqp3cetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.00001);
 319           sqp3cetaphitpc[aniter]->SetTPCExitSepMinimum(-1.);
 320           sqp3cetaphitpc[aniter]->SetKstarCut(0.04,AliFemtoV0TrackPairCut::kAntiLambda,AliFemtoV0TrackPairCut::kAntiProton); //1 - antilambda, 3 - antiproton
 321           sqp3cetaphitpc[aniter]->SetMinAvgSeparation(0,0); //0 - track-pos, 1 - track-neg
 322           sqp3cetaphitpc[aniter]->SetMinAvgSeparation(1,11);
 323
 324           sqp4cetaphitpc[aniter] = new AliFemtoV0TrackPairCut(); //lambda-antiproton, antilambda-proton
 325           sqp4cetaphitpc[aniter]->SetShareQualityMax(1.0); //between V0 daughter and track
 326           sqp4cetaphitpc[aniter]->SetShareFractionMax(0.05);
 327           sqp4cetaphitpc[aniter]->SetTPCOnly(kTRUE);
 328           sqp4cetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
 329           sqp4cetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.00001);
 330           sqp4cetaphitpc[aniter]->SetTPCExitSepMinimum(-1.);
 331           //SetMinAvgSeparation w if'ach ponizej
 332
 333           //sqp3cetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
 334           sqp5cetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
 335           sqp5cetaphitpc[aniter]->SetMagneticFieldSign(-1.0); //-1->field1, 1->field3
 336           sqp5cetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.012);
 337           sqp5cetaphitpc[aniter]->SetEtaDifferenceMinimum(0.017);
 338           sqp5cetaphitpc[aniter]->SetShareQualityMax(1.0);
 339           sqp5cetaphitpc[aniter]->SetShareFractionMax(0.05);
 340           sqp5cetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
 341           sqp5cetaphitpc[aniter]->SetMaxEEMinv(0.002);
 342           sqp5cetaphitpc[aniter]->SetMaxThetaDiff(0.008);
 343           sqp5cetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
 344           sqp5cetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.0001);
 345
 346           avgsepcorr[aniter] = new AliFemtoAvgSepCorrFctn(Form("Avgsep%stpcM%i", chrgs[ichg], imult),5000,0,500);
 347
 348
 349           if(ichg == 0) //V0LL
 350           {
 351             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 352             anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 353             anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
 354             anetaphitpc[aniter]->SetPairCut(sqp1cetaphitpc[aniter]);
 355             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kV0s);
 356           }
 357           else if(ichg == 1) //V0ALAL
 358           {
 359             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 360             anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
 361             anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
 362             anetaphitpc[aniter]->SetPairCut(sqp1cetaphitpc[aniter]);
 363             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kV0s);
 364           }
 365           else if(ichg == 2) //VOLAL
 366           {
 367             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 368             anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 369             anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
 370             anetaphitpc[aniter]->SetPairCut(sqp1cetaphitpc[aniter]);
 371             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kV0s);
 372           }
 373           else if(ichg == 3) //V0PL
 374           {
 375             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 376             anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 377             anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
 378             anetaphitpc[aniter]->SetPairCut(sqp2cetaphitpc[aniter]);
 379             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTrackV0);
 380           }
 381           else if(ichg == 4) //V0APL
 382           {
 383             sqp4cetaphitpc[aniter]->SetMinAvgSeparation(0,0); //0 - track-pos, 1 - track-neg
 384             sqp4cetaphitpc[aniter]->SetMinAvgSeparation(1,11);
 385             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 386             anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 387             anetaphitpc[aniter]->SetSecondParticleCut(dtc4etaphitpc[aniter]);
 388             anetaphitpc[aniter]->SetPairCut(sqp4cetaphitpc[aniter]);
 389             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTrackV0);
 390           }
 391           else if(ichg == 5) //V0PAL
 392           {
 393             sqp4cetaphitpc[aniter]->SetMinAvgSeparation(0,11); //0 - track-pos, 1 - track-neg
 394             sqp4cetaphitpc[aniter]->SetMinAvgSeparation(1,0);
 395             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 396             anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
 397             anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
 398             anetaphitpc[aniter]->SetPairCut(sqp4cetaphitpc[aniter]);
 399             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTrackV0);
 400           }
 401           else if(ichg == 6) //V0APAL
 402           {
 403             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 404             anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
 405             anetaphitpc[aniter]->SetSecondParticleCut(dtc4etaphitpc[aniter]);
 406             anetaphitpc[aniter]->SetPairCut(sqp3cetaphitpc[aniter]);
 407             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTrackV0);
 408           }
 409           else if(ichg ==7) //PP
 410           {
 411             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 412             anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
 413             anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
 414             anetaphitpc[aniter]->SetPairCut(sqp5cetaphitpc[aniter]);
 415             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTracks);
 416           }
 417           else if(ichg ==8) //PAP
 418           {
 419             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 420             anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
 421             anetaphitpc[aniter]->SetSecondParticleCut(dtc4etaphitpc[aniter]);
 422             anetaphitpc[aniter]->SetPairCut(sqp5cetaphitpc[aniter]);
 423             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTracks);
 424           }
 425           else if(ichg ==9) //APAP
 426           {
 427             anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 428             anetaphitpc[aniter]->SetFirstParticleCut(dtc4etaphitpc[aniter]);
 429             anetaphitpc[aniter]->SetSecondParticleCut(dtc4etaphitpc[aniter]);
 430             anetaphitpc[aniter]->SetPairCut(sqp5cetaphitpc[aniter]);
 431             avgsepcorr[aniter]->SetPairType(AliFemtoAvgSepCorrFctn::kTracks);
 432           }
 433
 434           anetaphitpc[aniter]->AddCorrFctn(avgsepcorr[aniter]);
 435
 436           //cylmetaphitpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),3,nbinssh,0.0,shqmax,runshlcms);
 437           //anetaphitpc[aniter]->AddCorrFctn(cylmetaphitpc[aniter]);
 438
 439           //          cqinvnclstpc[aniter] = new AliFemtoCorrFctnTPCNcls(Form("cqinvncls%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
 440           // anetaphitpc[aniter]->AddCorrFctn(cqinvnclstpc[aniter]);
 441
 442           //cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),45, 45);
 443           //anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
 444
 445
 446           // cqinvtitpc[aniter] = new AliFemtoTPCInnerCorrFctn(Form("cqinvti%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
 447           // anetaphitpc[aniter]->AddCorrFctn(cqinvtitpc[aniter]);
 448
 449           if(ichg==3 || ichg==4 || ichg==5 || ichg==6 || ichg==8) { //PL, APL, PAL, APAL, PAP
 450             cnonidtpc[aniter] = new AliFemtoCorrFctnNonIdDR(Form("cnonid%stpcM%i", chrgs[ichg], imult), nbinssh, 0.0,shqmax);
 451             anetaphitpc[aniter]->AddCorrFctn(cnonidtpc[aniter]);
 452           }
 453           else {
 454             cqinvtpc[aniter] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
 455             anetaphitpc[aniter]->AddCorrFctn(cqinvtpc[aniter]);
 456           }
 457
 458           // cylmkttpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),3, nbinssh, 0.0,shqmax, runshlcms);
 459           // anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[aniter]);
 460
 461           if (runktdep)
 462           {
 463             int ktm;
 464             for (int ikt=0; ikt<numOfkTBins; ikt++) {
 465               ktm = aniter*numOfkTBins + ikt;
 466               ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
 467
 468               cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,nbinssh, 0.0,shqmax,runshlcms);
 469               cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 470               anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[ktm]);
 471
 472               //cdedpetaphi[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ikT%i", chrgs[ichg], imult,ikt),45, 45);
 473               // cdedpetaphi[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 474               //anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[ktm]);
 475
 476               //if(ichg==3 || ichg==4 || ichg==5 || ichg==6 || ichg==8) { //PL, APL, PAL, APAL, PAP
 477               //cnonidtpc[ktm] = new AliFemtoCorrFctnNonIdDR(Form("cnonid%stpcM%ikT%i", chrgs[ichg], imult,ikt),nbinssh, 0.0,shqmax);
 478               //cnonidtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 479               //anetaphitpc[aniter]->AddCorrFctn(cnonidtpc[ktm]);
 480               //}
 481
 482               //else{
 483               //cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
 484               //cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 485               //anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
 486
 487               //}
 488
 489             }
 490           }
 491
 492           Manager->AddAnalysis(anetaphitpc[aniter]);
 493         }
 494       }
 495     }
 496   }
 497   // *** End of analysis
 498
 499   return Manager;
 500 }