PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train1/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 "AliFemtoPairCutPt.h"
  51 #endif
  52
  53 //________________________________________________________________________
  54 AliFemtoManager* ConfigFemtoAnalysis(const char* params="") {
  55
  56         double PionMass = 0.13956995;
  57         double KaonMass = 0.493677;
  58         double ProtonMass = 0.938272013;
  59
  60         const int numOfMultBins = 5;
  61         const int numOfChTypes = 13;
  62         const int numOfkTbins = 3;
  63
  64
  65
  66         char *parameter[10];
  67         if(strlen(params)!=0)
  68           {
  69             parameter[0] = strtok(params, ","); // Splits spaces between words in params
  70             cout<<"Parameter [0] (filterbit):"<<parameter[0]<<endl; // Writes first parameter
  71             parameter[1] = strtok(NULL, ",");
  72             cout<<"Parameter [1] (ktdep):"<<parameter[1]<<" "<<endl;
  73             parameter[2] = strtok(NULL, ",");
  74             cout<<"Parameter [2] (multdep):"<<parameter[2]<<" "<<endl;
  75             parameter[3] = strtok(NULL, ",");
  76             cout<<"Parameter [3]: (MinPlpContribSPD)"<<parameter[3]<<" "<<endl;
  77           }
  78
  79         int runmultdep = atoi(parameter[2]);
  80
  81         int runmults[numOfMultBins] = {0, 0, 0, 0, 1};
  82         if(runmultdep)    {         runmults[0]=1; runmults[1]=1; runmults[2]=1;          }
  83         int multbins[numOfMultBins+1] = {2, 20, 50,150,2,150};
  84
  85
  86         int runch[numOfChTypes] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
  87         const char *chrgs[numOfChTypes] = { "PP", "aPaP", "PaP", "KpKp", "KmKm", "KpKm", "PIpPIp", "PImPIm", "PIpPIm", "all", "plus", "minus", "mixed" };
  88
  89
  90         int runktdep = atoi(parameter[1]);
  91         double ktrng[numOfkTbins+1] = {0.0, 0.75, 2.0, 100.0};
  92
  93         int runqinv = 1;
  94         int runshlcms = 1;// 0:PRF(PAP), 1:LCMS(PP,APAP)
  95
  96         int runtype = 0; // Types 0 - global, 1 - ITS only, 2 - TPC Inner       //global tracks ->mfit ITS+TPC
  97         int owncuts = 0;
  98         int owndca = 0;
  99
 100         int gammacut = 0;       // cut na ee z gamma
 101
 102         double shqmax = 0.5;
 103         int nbinssh = 100;
 104
 105         //AliFemtoEventReaderESDChain *Reader = new AliFemtoEventReaderESDChain();
 106         //Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kGlobalCount);
 107
 108
 109         AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
 110         Reader->SetFilterMask(atoi(parameter[0]));
 111         Reader->SetDCAglobalTrack(kTRUE); //false for FB7, true for the rest //we do not use DCA at all
 112         Reader->SetUseMultiplicity(AliFemtoEventReaderAOD::kReference);
 113         Reader->SetMinPlpContribSPD(atoi(parameter[3]));
 114         Reader->SetIsPileUpEvent(kTRUE);
 115
 116         AliFemtoManager* Manager = new AliFemtoManager();
 117         Manager->SetEventReader(Reader);
 118
 119         AliFemtoVertexMultAnalysis              *anetaphitpc[320];
 120         AliFemtoBasicEventCut                            *mecetaphitpc[320];
 121         AliFemtoCutMonitorEventMult      *cutPassEvMetaphitpc[320];
 122         AliFemtoCutMonitorEventMult      *cutFailEvMetaphitpc[320];
 123         AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[320];
 124         AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[320];
 125         AliFemtoESDTrackCut                      *dtc1etaphitpc[320];
 126         AliFemtoESDTrackCut                      *dtc2etaphitpc[320];
 127         AliFemtoESDTrackCut                      *dtc3etaphitpc[320];
 128         AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[320];
 129         AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[320];
 130         AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[320];
 131         AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[320];
 132         AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[320];
 133         AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[320];
 134         AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[320];
 135         AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[320];
 136         AliFemtoCutMonitorParticleYPt *cutPass3YPtetaphitpc[320];
 137         AliFemtoCutMonitorParticleYPt *cutFail3YPtetaphitpc[320];
 138         AliFemtoCutMonitorParticlePID *cutPass3PIDetaphitpc[320];
 139         AliFemtoCutMonitorParticlePID *cutFail3PIDetaphitpc[320];
 140         //       AliFemtoShareQualityTPCEntranceSepPairCut                      *sqpcetaphitpcsame[320];
 141         //AliFemtoPairCutAntiGamma                      *sqpcetaphitpc[320];
 142         //AliFemtoPairCutRadialDistance                 *sqpcetaphitpc[320];
 143         AliFemtoShareQualityPairCut                     *sqpcetaphitpc[320];
 144         //      AliFemtoChi2CorrFctn                                    *cchiqinvetaphitpc[320];
 145         AliFemtoPairCutPt                                                *ktpcuts[320];
 146         AliFemtoQinvCorrFctn                                    *cqinvkttpc[320];
 147         AliFemtoQinvCorrFctn                                    *cqinvtpc[320];
 148         AliFemtoCorrFctnDEtaDPhi                        *cdedpetaphi[320];
 149
 150
 151
 152         // *** Third QA task - HBT analysis with all pair cuts off, TPC only ***
 153         // *** Begin pion-pion (positive) analysis ***
 154         int aniter = 0;
 155
 156         for (int imult = 0; imult < numOfMultBins; imult++)
 157         {
 158                 if (runmults[imult])
 159                 {
 160                         for (int ichg = 0; ichg < numOfChTypes; ichg++)
 161                         {
 162                                 if (runch[ichg])
 163                                 {
 164
 165                                         aniter = ichg * numOfMultBins + imult;
 166                                         anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(10, -10.0, 10.0, 30, multbins[imult], multbins[imult+1]);
 167                                         anetaphitpc[aniter]->SetNumEventsToMix(10);
 168                                         anetaphitpc[aniter]->SetMinSizePartCollection(1);
 169                                         anetaphitpc[aniter]->SetVerboseMode(kFALSE);
 170
 171                                         //*** Event cut ***
 172                                         mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
 173                                         mecetaphitpc[aniter]->SetEventMult(0.001,100000);
 174                                         mecetaphitpc[aniter]->SetVertZPos(-10,10);//cm
 175
 176                                         //****** event monitors **********
 177                                         cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
 178                                         cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
 179                                         mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
 180
 181                                         cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
 182                                         cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
 183                                         mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
 184
 185
 186                                         // ***** single particle track cuts *********
 187                                         dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
 188                                         dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
 189                                         dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();
 190
 191                                         dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
 192                                         dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
 193                                         dtc3etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
 194                                         dtc1etaphitpc[aniter]->SetNsigma(3.0);
 195                                         dtc2etaphitpc[aniter]->SetNsigma(3.0);
 196                                         dtc3etaphitpc[aniter]->SetNsigma(3.0);
 197
 198                                         dtc1etaphitpc[aniter]->SetCharge(1.0);
 199                                         dtc2etaphitpc[aniter]->SetCharge(-1.0);
 200
 201                                         dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
 202                                         dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
 203                                         dtc3etaphitpc[aniter]->SetEta(-0.8,0.8);
 204
 205                                         dtc1etaphitpc[aniter]->SetElectronRejection(true);
 206                                         dtc2etaphitpc[aniter]->SetElectronRejection(true);
 207                                         dtc3etaphitpc[aniter]->SetElectronRejection(true);
 208
 209                                         if (ichg == 0 ||ichg == 1 ||ichg == 2)//protons 0-2
 210                                           {
 211                                             dtc1etaphitpc[aniter]->SetPt(0.5,4);
 212                                             dtc2etaphitpc[aniter]->SetPt(0.5,4);
 213
 214                                             dtc1etaphitpc[aniter]->SetMass(ProtonMass);
 215                                             dtc1etaphitpc[aniter]->SetMostProbableProton();
 216                                             dtc2etaphitpc[aniter]->SetMass(ProtonMass);
 217                                             dtc2etaphitpc[aniter]->SetMostProbableProton();
 218                                           }
 219
 220                                         if (ichg == 3 ||ichg == 4 ||ichg == 5)//kaons 3-5
 221                                           {
 222                                             dtc1etaphitpc[aniter]->SetPt(0.3,4);
 223                                             dtc2etaphitpc[aniter]->SetPt(0.3,4);
 224
 225                                             dtc1etaphitpc[aniter]->SetMass(KaonMass);
 226                                             dtc1etaphitpc[aniter]->SetMostProbableKaon();
 227                                             dtc2etaphitpc[aniter]->SetMass(KaonMass);
 228                                             dtc2etaphitpc[aniter]->SetMostProbableKaon();
 229
 230                                           }
 231                                         if (ichg == 6 ||ichg == 7 ||ichg == 8)//pions 6-8
 232                                           {
 233                                             dtc1etaphitpc[aniter]->SetPt(0.2,4);
 234                                             dtc2etaphitpc[aniter]->SetPt(0.2,4);
 235
 236                                             dtc1etaphitpc[aniter]->SetMass(PionMass);
 237                                             dtc1etaphitpc[aniter]->SetMostProbablePion();
 238                                             dtc2etaphitpc[aniter]->SetMass(PionMass);
 239                                             dtc2etaphitpc[aniter]->SetMostProbablePion();
 240                                           }
 241                                         if (ichg == 9)//all
 242                                           {
 243                                             dtc3etaphitpc[aniter]->SetPt(0.2,4);
 244                                           }
 245                                         if (ichg == 10 ||ichg == 11 ||ichg == 12)//plus,minus,mixed
 246                                           {
 247                                             dtc1etaphitpc[aniter]->SetPt(0.2,4);
 248                                             dtc2etaphitpc[aniter]->SetPt(0.2,4);
 249                                           }
 250
 251                                         //****** DCA ******
 252
 253                                         if(owndca){
 254                                           dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01);      //      DCA xy
 255                                           //dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
 256                                           dtc1etaphitpc[aniter]->SetMaxImpactZ(2);      //DCA Z
 257                                           dtc2etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01);      //      DCA xy
 258                                           dtc2etaphitpc[aniter]->SetMaxImpactZ(2);      //DCA Z
 259                                           dtc3etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01);      //      DCA xy
 260                                           dtc3etaphitpc[aniter]->SetMaxImpactZ(2);      //DCA Z
 261                                         }
 262                                         //****** Track quality cuts ******
 263
 264                                         if(owncuts){
 265
 266                                           dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 267                                           dtc1etaphitpc[aniter]->SetminTPCncls(70);
 268                                           dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
 269                                           dtc1etaphitpc[aniter]->SetLabel(kFALSE);
 270                                           //    dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
 271                                           dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0); // pisac
 272                                           //dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
 273
 274                                           dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 275                                           dtc2etaphitpc[aniter]->SetminTPCncls(70);
 276                                           dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
 277                                           dtc2etaphitpc[aniter]->SetLabel(kFALSE);
 278                                           //    dtc2etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
 279                                           dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
 280                                           //    dtc2etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
 281
 282                                           dtc3etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 283                                           dtc3etaphitpc[aniter]->SetminTPCncls(70);
 284                                           dtc3etaphitpc[aniter]->SetRemoveKinks(kTRUE);
 285                                           dtc3etaphitpc[aniter]->SetLabel(kFALSE);
 286                                           //    dtc3etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
 287                                           dtc3etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
 288                                           //    dtc3etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
 289
 290                                         }
 291                                         //**************** track Monitors ***************
 292
 293
 294                                         if(1)//ichg>8)
 295                                           {
 296                                             cutPass3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass%stpcM%i", chrgs[ichg], imult),PionMass);
 297                                             cutFail3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail%stpcM%i", chrgs[ichg], imult),PionMass);
 298                                             if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
 299                                             if(ichg==0||ichg==3||ichg==6||ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
 300                                             if(ichg==1||ichg==4||ichg==7||ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
 301
 302                                             if(ichg<3){
 303                                               cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
 304                                               cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),2);
 305                                             }
 306                                             else if(ichg>=3 && ichg<=5){
 307                                               cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),1);//0-pion,1-kaon,2-proton
 308                                               cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),1);
 309                                             }
 310                                             else{
 311                                               cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),0);//0-pion,1-kaon,2-proton
 312                                               cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),0);
 313
 314                                             }
 315                                             if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
 316                                             if(ichg==0||ichg==3||ichg==6||ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
 317                                             if(ichg==1||ichg==4||ichg==7||ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
 318                                           }
 319
 320                                         //******** Two - track cuts ************
 321                                         //sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
 322                                         //sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
 323                                         sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
 324
 325                                         sqpcetaphitpc[aniter]->SetShareQualityMax(0.0);         // two track cuts on splitting and merging  //1- wylaczany 0 -wlaczany
 326                                         sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);       //  ile moga miec wspolnych klastrow //1 - wylaczany, 0.05 - wlaczany
 327                                         sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
 328                                         /*sqpcetaphitpc[aniter]->SetMaximumRadius(2.51);
 329                                         sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
 330                                         sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.02);
 331                                         sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.2);
 332
 333                                         if (gammacut == 0)
 334                                           {
 335                                             sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
 336                                             sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
 337                                           }
 338                                         else if (gammacut == 1)
 339                                           {
 340                                             sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
 341                                             sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
 342                                             }*/
 343
 344                                         // sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
 345                                         // sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
 346                                         // sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
 347
 348
 349
 350
 351                                         //***** Setting cuts ***********
 352
 353
 354                                         // setting event cut
 355                                         anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
 356                                         //setting single track cuts
 357                                         if(ichg==0 || ichg==3  || ichg==6 || ichg==10) //positive like-sign
 358                                         {
 359                                                 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 360                                                 anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
 361                                         }
 362                                         if(ichg==1 || ichg==4 || ichg==7 || ichg==11)//negative like-sign
 363                                         {
 364                                                 anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
 365                                                 anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
 366                                         }
 367                                         if(ichg==2 || ichg==5 || ichg==8 || ichg==12)//unlike-sign
 368                                         {
 369                                                 anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
 370                                                 anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
 371                                         }
 372                                         if(ichg==9) //all
 373                                         {
 374                                           anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
 375                                           anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
 376                                         }
 377
 378                                         //setting two-track cuts
 379                                         anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
 380
 381
 382                                         //**** Correlation functions *******
 383
 384                                         cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),35, 35);
 385                                         anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
 386
 387
 388                                         if (runktdep)
 389                                         {
 390                                                 int ktm;
 391                                                 for (int ikt=0; ikt<numOfkTbins; ikt++)
 392                                                 {
 393                                                         ktm = aniter * numOfkTbins + ikt;
 394                                                         ktpcuts[ktm] = new AliFemtoPairCutPt(ktrng[ikt], ktrng[ikt+1]);
 395
 396                                                         //cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
 397                                                         //cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 398                                                         //anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
 399
 400                                                         cdedpetaphi[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ipT%i", chrgs[ichg], imult,ikt),35, 35);
 401                                                         cdedpetaphi[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
 402                                                         anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[ktm]);
 403
 404                                                 }
 405                                         }
 406                                         Manager->AddAnalysis(anetaphitpc[aniter]);
 407                                 }
 408                         }
 409                 }
 410         }
 411         return Manager;
 412 }