PWGCF/FEMTOSCOPY/macros/Train/DataTrain1/ConfigFemtoAnalysis.C

   1 /*********************************************************************
   2  *                                                                   *
   3  * ConfigFemtoAnalysis.C - configuration macro for the femtoscopic   *
   4  * analysis, to be run in the analysis train.                        *
   5  * Assumed input data: large (>10M) sample of MC pp events           *
   6  * Inluded analysis:                                                 *
   7  *    - positive pion HBT, 3 kt bins, 1D+3D functions                *
   8  *    - negative pion HBT, 3 kt bins, 1D+3D functions                *
   9  *    - positive kaon HBT, 1 kt bin,  1D+3D functions                *
  10  *                                                                   *
  11  * Author: Adam Kisiel (Adam.Kisiel@cern.ch)                         *
  12  *                                                                   *
  13  *********************************************************************/
  14
  15 #if !defined(__CINT__) || defined(__MAKECINT_)
  16 #include "AliFemtoManager.h"
  17 #include "AliFemtoEventReaderESDChain.h"
  18 #include "AliFemtoEventReaderESDChainKine.h"
  19 #include "AliFemtoSimpleAnalysis.h"
  20 #include "AliFemtoBasicEventCut.h"
  21 #include "AliFemtoESDTrackCut.h"
  22 #include "AliFemtoCorrFctn.h"
  23 #include "AliFemtoCutMonitorParticleYPt.h"
  24 #include "AliFemtoCutMonitorParticleVertPos.h"
  25 #include "AliFemtoCutMonitorParticleMomRes.h"
  26 #include "AliFemtoCutMonitorEventMult.h"
  27 #include "AliFemtoCutMonitorEventVertex.h"
  28 #include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
  29 #include "AliFemtoQinvCorrFctn.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 "AliFemtoModelBPLCMSCorrFctn.h"
  37 #include "AliFemtoModelCorrFctn3DSpherical.h"
  38 #include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
  39 #include "AliFemtoModelGausRinvFreezeOutGenerator.h"
  40 #include "AliFemtoModelManager.h"
  41 #include "AliFemtoModelWeightGeneratorBasic.h"
  42 #include "AliFemtoModelWeightGeneratorLednicky.h"
  43 #include "AliFemtoCorrFctnDirectYlm.h"
  44 #include "AliFemtoModelCorrFctnDirectYlm.h"
  45 #include "AliFemtoModelCorrFctnSource.h"
  46 #include "AliFemtoCutMonitorParticlePtPDG.h"
  47 #include "AliFemtoKTPairCut.h"
  48 #endif
  49
  50 //________________________________________________________________________
  51 AliFemtoManager* ConfigFemtoAnalysis() {
  52
  53   double PionMass = 0.13956995;
  54   double KaonMass = 0.493677;
  55
  56   AliFemtoEventReaderESDChain* Reader=new AliFemtoEventReaderESDChain();
  57   Reader->SetConstrained(true);
  58   Reader->SetUseTPCOnly(false);
  59
  60   AliFemtoManager* Manager=new AliFemtoManager();
  61   Manager->SetEventReader(Reader);
  62
  63   int runPositivePions = 1;
  64   int runNegativePions = 1;
  65   int runPositiveKaons = 1;
  66   int runNegativeKaons = 1;
  67   int runPositiveNegativeKaons = 1;
  68
  69   // Z vertex mixing settings
  70   int nZVertexBins = 4;
  71   Double_t minZVertex = -10.0;
  72   Double_t maxZVertex = 10.0;
  73
  74   if (runPositivePions) {
  75     // *** Begin pion-pion (positive) analysis ***
  76     AliFemtoVertexMultAnalysis *anpip = new AliFemtoVertexMultAnalysis(nZVertexBins, minZVertex, maxZVertex, 1, 2, 200000);
  77     anpip->SetNumEventsToMix(10);
  78     anpip->SetMinSizePartCollection(2);
  79
  80     AliFemtoBasicEventCut* mecpip = new AliFemtoBasicEventCut();
  81     mecpip->SetEventMult(2,100000);
  82     mecpip->SetVertZPos(-1000,1000);
  83
  84     AliFemtoESDTrackCut* dtcpip = new AliFemtoESDTrackCut();
  85     dtcpip->SetPidProbPion(0.2,1.001);
  86     dtcpip->SetPidProbMuon(0.0,1.0);
  87     dtcpip->SetPidProbKaon(0.0,1.0);
  88     dtcpip->SetPidProbProton(0.0,1.0);
  89     dtcpip->SetCharge(1.0);
  90     dtcpip->SetPt(0.05,1.0);
  91     dtcpip->SetRapidity(-0.9,0.9);
  92     dtcpip->SetMass(PionMass);
  93     // Track quality cuts
  94     dtcpip->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
  95     //  dtcpip->SetStatus(AliESDtrack::kTPCrefit);
  96     dtcpip->SetminTPCncls(50);
  97     dtcpip->SetRemoveKinks(kTRUE);
  98     dtcpip->SetLabel(kFALSE);
  99     dtcpip->SetMaxITSChiNdof(2.5);
 100     dtcpip->SetMaxTPCChiNdof(3.0);
 101     dtcpip->SetMaxImpactXY(3.0);
 102     dtcpip->SetMaxImpactZ(3.0);
 103
 104     // Track monitors
 105     AliFemtoCutMonitorParticleYPt *cutPassYPtpip = new AliFemtoCutMonitorParticleYPt("cutPasspip", 0.13957);
 106     AliFemtoCutMonitorParticleYPt *cutFailYPtpip = new AliFemtoCutMonitorParticleYPt("cutFailpip", 0.13957);
 107     dtcpip->AddCutMonitor(cutPassYPtpip, cutFailYPtpip);
 108
 109     AliFemtoCutMonitorEventMult *cutPassEvMpip = new AliFemtoCutMonitorEventMult("cutPasspip");
 110     AliFemtoCutMonitorEventMult *cutFailEvMpip = new AliFemtoCutMonitorEventMult("cutFailpip");
 111     mecpip->AddCutMonitor(cutPassEvMpip, cutFailEvMpip);
 112
 113     AliFemtoCutMonitorEventVertex *cutPassEvVpip = new AliFemtoCutMonitorEventVertex("cutPasspip");
 114     AliFemtoCutMonitorEventVertex *cutFailEvVpip = new AliFemtoCutMonitorEventVertex("cutFailpip");
 115     mecpip->AddCutMonitor(cutPassEvVpip, cutFailEvVpip);
 116
 117     // Pair cut
 118     AliFemtoShareQualityTPCEntranceSepPairCut *sqpcpip = new AliFemtoShareQualityTPCEntranceSepPairCut();
 119     sqpcpip->SetShareQualityMax(0.0);
 120     sqpcpip->SetShareFractionMax(0.02);
 121     sqpcpip->SetRemoveSameLabel(kFALSE);
 122     sqpcpip->SetTPCEntranceSepMinimum(2.0);
 123
 124     anpip->SetEventCut(mecpip);
 125     anpip->SetFirstParticleCut(dtcpip);
 126     anpip->SetSecondParticleCut(dtcpip);
 127     anpip->SetPairCut(sqpcpip);
 128
 129     // Two-track quality monitoring
 130     AliFemtoShareQualityCorrFctn *csqqinvpip= new AliFemtoShareQualityCorrFctn("sqqinvcfpip",40,0.0,0.4);
 131     AliFemtoChi2CorrFctn *cchiqinvpip= new AliFemtoChi2CorrFctn("chicfpip",40,0.0,0.4);
 132     AliFemtoCorrFctnTPCNcls *cqtpcnclspip = new AliFemtoCorrFctnTPCNcls("cqtpcnclspip",40,0.0,0.4);
 133
 134     // Intrdouce kT binning
 135     AliFemtoKTPairCut *ktpairkT1pip = new AliFemtoKTPairCut(0.05,0.27);
 136     AliFemtoKTPairCut *ktpairkT2pip = new AliFemtoKTPairCut(0.27,0.37);
 137     AliFemtoKTPairCut *ktpairkT3pip = new AliFemtoKTPairCut(0.37,1.0);
 138
 139     // Purely experimental correlation function
 140     AliFemtoCorrFctnDirectYlm *cylmkT1pip = new AliFemtoCorrFctnDirectYlm("cylmkT1pip",3,80,0.0,0.8,1);
 141     cylmkT1pip->SetPairSelectionCut(ktpairkT1pip);
 142     anpip->AddCorrFctn(cylmkT1pip);
 143
 144     AliFemtoCorrFctnDirectYlm *cylmkT2pip = new AliFemtoCorrFctnDirectYlm("cylmkT2pip",3,80,0.0,0.8,1);
 145     cylmkT2pip->SetPairSelectionCut(ktpairkT2pip);
 146     anpip->AddCorrFctn(cylmkT2pip);
 147
 148     AliFemtoCorrFctnDirectYlm *cylmkT3pip = new AliFemtoCorrFctnDirectYlm("cylmkT3pip",3,80,0.0,0.8,1);
 149     cylmkT3pip->SetPairSelectionCut(ktpairkT3pip);
 150     anpip->AddCorrFctn(cylmkT3pip);
 151
 152     AliFemtoQinvCorrFctn *cqinvkt1pip = new AliFemtoQinvCorrFctn("qinvcfkt1pip", 100,0.0,1.0);
 153     cqinvkt1pip->SetPairSelectionCut(ktpairkT1pip);
 154     anpip->AddCorrFctn(cqinvkt1pip);
 155
 156     AliFemtoQinvCorrFctn *cqinvkt2pip = new AliFemtoQinvCorrFctn("qinvcfkt2pip", 100,0.0,1.0);
 157     cqinvkt2pip->SetPairSelectionCut(ktpairkT2pip);
 158     anpip->AddCorrFctn(cqinvkt2pip);
 159
 160     AliFemtoQinvCorrFctn *cqinvkt3pip = new AliFemtoQinvCorrFctn("qinvcfkt3pip", 100,0.0,1.0);
 161     cqinvkt3pip->SetPairSelectionCut(ktpairkT3pip);
 162     anpip->AddCorrFctn(cqinvkt3pip);
 163
 164     // Add correlation functions to the analysis
 165     anpip->AddCorrFctn(csqqinvpip);
 166     anpip->AddCorrFctn(cchiqinvpip);
 167     anpip->AddCorrFctn(cqtpcnclspip);
 168
 169     Manager->AddAnalysis(anpip);
 170
 171     // *** End pion-pion (positive) analysis
 172   }
 173
 174   if (runNegativePions) {
 175     // *** Begin pion-pion (negative) analysis ***
 176     AliFemtoVertexMultAnalysis *anpim = new AliFemtoVertexMultAnalysis(nZVertexBins, minZVertex, maxZVertex, 1, 2, 200000);
 177     anpim->SetNumEventsToMix(10);
 178     anpim->SetMinSizePartCollection(2);
 179
 180     AliFemtoBasicEventCut* mecpim = new AliFemtoBasicEventCut();
 181     mecpim->SetEventMult(2,100000);
 182     mecpim->SetVertZPos(-1000,1000);
 183
 184     AliFemtoESDTrackCut* dtcpim = new AliFemtoESDTrackCut();
 185     dtcpim->SetPidProbPion(0.2,1.001);
 186     dtcpim->SetPidProbMuon(0.0,1.0);
 187     dtcpim->SetPidProbKaon(0.0,1.0);
 188     dtcpim->SetPidProbProton(0.0,1.0);
 189     dtcpim->SetCharge(1.0);
 190     dtcpim->SetPt(0.05,1.0);
 191     dtcpim->SetRapidity(-0.9,0.9);
 192     dtcpim->SetMass(PionMass);
 193     // Track quality cuts
 194     dtcpim->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 195     //  dtcpim->SetStatus(AliESDtrack::kTPCrefit);
 196     dtcpim->SetminTPCncls(50);
 197     dtcpim->SetRemoveKinks(kTRUE);
 198     dtcpim->SetLabel(kFALSE);
 199     dtcpim->SetMaxITSChiNdof(2.5);
 200     dtcpim->SetMaxTPCChiNdof(3.0);
 201     dtcpim->SetMaxImpactXY(3.0);
 202     dtcpim->SetMaxImpactZ(3.0);
 203
 204     // Track monitors
 205     AliFemtoCutMonitorParticleYPt *cutPassYPtpim = new AliFemtoCutMonitorParticleYPt("cutPasspim", 0.13957);
 206     AliFemtoCutMonitorParticleYPt *cutFailYPtpim = new AliFemtoCutMonitorParticleYPt("cutFailpim", 0.13957);
 207     dtcpim->AddCutMonitor(cutPassYPtpim, cutFailYPtpim);
 208
 209     AliFemtoCutMonitorEventMult *cutPassEvMpim = new AliFemtoCutMonitorEventMult("cutPasspim");
 210     AliFemtoCutMonitorEventMult *cutFailEvMpim = new AliFemtoCutMonitorEventMult("cutFailpim");
 211     mecpim->AddCutMonitor(cutPassEvMpim, cutFailEvMpim);
 212
 213     AliFemtoCutMonitorEventVertex *cutPassEvVpim = new AliFemtoCutMonitorEventVertex("cutPasspim");
 214     AliFemtoCutMonitorEventVertex *cutFailEvVpim = new AliFemtoCutMonitorEventVertex("cutFailpim");
 215     mecpim->AddCutMonitor(cutPassEvVpim, cutFailEvVpim);
 216
 217     // Pair cut
 218     AliFemtoShareQualityTPCEntranceSepPairCut *sqpcpim = new AliFemtoShareQualityTPCEntranceSepPairCut();
 219     sqpcpim->SetShareQualityMax(0.0);
 220     sqpcpim->SetShareFractionMax(0.02);
 221     sqpcpim->SetRemoveSameLabel(kFALSE);
 222     sqpcpim->SetTPCEntranceSepMinimum(2.0);
 223
 224     anpim->SetEventCut(mecpim);
 225     anpim->SetFirstParticleCut(dtcpim);
 226     anpim->SetSecondParticleCut(dtcpim);
 227     anpim->SetPairCut(sqpcpim);
 228
 229     // Two-track quality monitoring
 230     AliFemtoShareQualityCorrFctn *csqqinvpim= new AliFemtoShareQualityCorrFctn("sqqinvcfpim",40,0.0,0.4);
 231     AliFemtoChi2CorrFctn *cchiqinvpim= new AliFemtoChi2CorrFctn("chicfpim",40,0.0,0.4);
 232     AliFemtoCorrFctnTPCNcls *cqtpcnclspim = new AliFemtoCorrFctnTPCNcls("cqtpcnclspim",40,0.0,0.4);
 233
 234     // Intrdouce kT binning
 235     AliFemtoKTPairCut *ktpairkT1pim = new AliFemtoKTPairCut(0.05,0.27);
 236     AliFemtoKTPairCut *ktpairkT2pim = new AliFemtoKTPairCut(0.27,0.37);
 237     AliFemtoKTPairCut *ktpairkT3pim = new AliFemtoKTPairCut(0.37,1.0);
 238
 239     // Purely experimental correlation function
 240     AliFemtoCorrFctnDirectYlm *cylmkT1pim = new AliFemtoCorrFctnDirectYlm("cylmkT1pim",3,80,0.0,0.8,1);
 241     cylmkT1pim->SetPairSelectionCut(ktpairkT1pim);
 242     anpim->AddCorrFctn(cylmkT1pim);
 243
 244     AliFemtoCorrFctnDirectYlm *cylmkT2pim = new AliFemtoCorrFctnDirectYlm("cylmkT2pim",3,80,0.0,0.8,1);
 245     cylmkT2pim->SetPairSelectionCut(ktpairkT2pim);
 246     anpim->AddCorrFctn(cylmkT2pim);
 247
 248     AliFemtoCorrFctnDirectYlm *cylmkT3pim = new AliFemtoCorrFctnDirectYlm("cylmkT3pim",3,80,0.0,0.8,1);
 249     cylmkT3pim->SetPairSelectionCut(ktpairkT3pim);
 250     anpim->AddCorrFctn(cylmkT3pim);
 251
 252     AliFemtoQinvCorrFctn *cqinvkt1pim = new AliFemtoQinvCorrFctn("qinvcfkt1pim", 100,0.0,1.0);
 253     cqinvkt1pim->SetPairSelectionCut(ktpairkT1pim);
 254     anpim->AddCorrFctn(cqinvkt1pim);
 255
 256     AliFemtoQinvCorrFctn *cqinvkt2pim = new AliFemtoQinvCorrFctn("qinvcfkt2pim", 100,0.0,1.0);
 257     cqinvkt2pim->SetPairSelectionCut(ktpairkT2pim);
 258     anpim->AddCorrFctn(cqinvkt2pim);
 259
 260     AliFemtoQinvCorrFctn *cqinvkt3pim = new AliFemtoQinvCorrFctn("qinvcfkt3pim", 100,0.0,1.0);
 261     cqinvkt3pim->SetPairSelectionCut(ktpairkT3pim);
 262     anpim->AddCorrFctn(cqinvkt3pim);
 263
 264     // Add correlation functions to the analysis
 265     anpim->AddCorrFctn(csqqinvpim);
 266     anpim->AddCorrFctn(cchiqinvpim);
 267     anpim->AddCorrFctn(cqtpcnclspim);
 268
 269     Manager->AddAnalysis(anpim);
 270
 271     // *** End pion-pion (negative) analysis
 272   }
 273
 274   if (runPositiveKaons) {
 275     // *** Begin Kaon-Kaon (positive) analysis
 276     AliFemtoVertexMultAnalysis *ankp = new AliFemtoVertexMultAnalysis(nZVertexBins, minZVertex, maxZVertex, 1, 2, 20000);
 277     ankp->SetNumEventsToMix(5);
 278     ankp->SetMinSizePartCollection(2);
 279
 280     AliFemtoBasicEventCut* meckp = new AliFemtoBasicEventCut();
 281     meckp->SetEventMult(1,100000);
 282     meckp->SetVertZPos(-1000,1000);
 283
 284     AliFemtoESDTrackCut* dtckp = new AliFemtoESDTrackCut();
 285     dtckp->SetPidProbKaon(0.7,1.001);
 286     dtckp->SetPidProbMuon(0.0,0.5);
 287     dtckp->SetPidProbPion(0.0,0.5);
 288     dtckp->SetPidProbProton(0.0,0.5);
 289     dtckp->SetCharge(1.0);
 290     dtckp->SetMostProbableKaon();
 291     dtckp->SetMomRangeTOFpidIs(0.6,10000.);
 292     dtckp->SetPt(0.15,2.0);
 293     dtckp->SetMass(KaonMass);
 294     // Track quality cuts
 295     dtckp->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 296     dtckp->SetminTPCncls(50);
 297     dtckp->SetRemoveKinks(kTRUE);
 298     dtckp->SetLabel(kFALSE);
 299     dtckp->SetMaxITSChiNdof(2.5);
 300     dtckp->SetMaxTPCChiNdof(3.0);
 301     dtckp->SetMaxImpactXY(3.0);
 302     dtckp->SetMaxImpactZ(3.0);
 303
 304     AliFemtoCutMonitorParticleYPt *cutPassYPtkp = new AliFemtoCutMonitorParticleYPt("cutPasskp", 0.493677);
 305     AliFemtoCutMonitorParticleYPt *cutFailYPtkp = new AliFemtoCutMonitorParticleYPt("cutFailkp", 0.493677);
 306     dtckp->AddCutMonitor(cutPassYPtkp, cutFailYPtkp);
 307
 308     AliFemtoCutMonitorParticlePtPDG *cutPassPidkp = new AliFemtoCutMonitorParticlePtPDG("cutPasskp", 0.493677);
 309     AliFemtoCutMonitorParticlePtPDG *cutFailPidkp = new AliFemtoCutMonitorParticlePtPDG("cutFailkp", 0.493677);
 310     dtckp->AddCutMonitor(cutPassPidkp, cutFailPidkp);
 311
 312     AliFemtoCutMonitorParticleMomRes *cutPassMRkp = new AliFemtoCutMonitorParticleMomRes("cutPasskp");
 313     AliFemtoCutMonitorParticleMomRes *cutFailMRkp = new AliFemtoCutMonitorParticleMomRes("cutFailkp");
 314     dtckp->AddCutMonitor(cutPassMRkp, cutFailMRkp);
 315
 316     AliFemtoCutMonitorParticleVertPos *cutPassVPkp = new AliFemtoCutMonitorParticleVertPos("cutPasskp");
 317     AliFemtoCutMonitorParticleVertPos *cutFailVPkp = new AliFemtoCutMonitorParticleVertPos("cutFailkp");
 318     dtckp->AddCutMonitor(cutPassVPkp, cutFailVPkp);
 319
 320     AliFemtoCutMonitorEventMult *cutPassEvMkp = new AliFemtoCutMonitorEventMult("cutPasskp");
 321     AliFemtoCutMonitorEventMult *cutFailEvMkp = new AliFemtoCutMonitorEventMult("cutFailkp");
 322     meckp->AddCutMonitor(cutPassEvMkp, cutFailEvMkp);
 323
 324     AliFemtoCutMonitorEventVertex *cutPassEvVkp = new AliFemtoCutMonitorEventVertex("cutPasskp");
 325     AliFemtoCutMonitorEventVertex *cutFailEvVkp = new AliFemtoCutMonitorEventVertex("cutFailkp");
 326     meckp->AddCutMonitor(cutPassEvVkp, cutFailEvVkp);
 327
 328     AliFemtoShareQualityTPCEntranceSepPairCut *sqpckp = new AliFemtoShareQualityTPCEntranceSepPairCut();
 329     sqpckp->SetShareQualityMax(0.0);
 330     sqpckp->SetShareFractionMax(0.02);
 331     sqpckp->SetRemoveSameLabel(kFALSE);
 332     sqpckp->SetTPCEntranceSepMinimum(3.0);
 333
 334     ankp->SetEventCut(meckp);
 335     ankp->SetFirstParticleCut(dtckp);
 336     ankp->SetSecondParticleCut(dtckp);
 337     ankp->SetPairCut(sqpckp);
 338
 339     AliFemtoQinvCorrFctn *cqinvkp= new AliFemtoQinvCorrFctn("qinvcf",100,0.0,1.0);
 340
 341     //###
 342     ankp->AddCorrFctn(cqinvkp);
 343
 344     Manager->AddAnalysis(ankp);
 345
 346     // *** End Kaon-Kaon (positive) analysis
 347   }
 348
 349   if (runNegativeKaons) {
 350     // *** Begin Kaon-Kaon (negative) analysis
 351     AliFemtoVertexMultAnalysis *ankm = new AliFemtoVertexMultAnalysis(nZVertexBins, minZVertex, maxZVertex, 1, 2, 20000);
 352     ankm->SetNumEventsToMix(5);
 353     ankm->SetMinSizePartCollection(2);
 354
 355     AliFemtoBasicEventCut* meckm = new AliFemtoBasicEventCut();
 356     meckm->SetEventMult(1,100000);
 357     meckm->SetVertZPos(-1000,1000);
 358
 359     AliFemtoESDTrackCut* dtckm = new AliFemtoESDTrackCut();
 360     dtckm->SetPidProbKaon(0.7,1.001);
 361     dtckm->SetPidProbMuon(0.0,0.5);
 362     dtckm->SetPidProbPion(0.0,0.5);
 363     dtckm->SetPidProbProton(0.0,0.5);
 364     dtckp->SetMomRangeTOFpidIs(0.6,10000.);
 365     dtckm->SetCharge(-1.0);
 366     dtckm->SetMostProbableKaon();
 367     dtckm->SetPt(0.15,2.0);
 368     dtckm->SetMass(KaonMass);
 369     // Track quality cuts
 370     dtckm->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 371     dtckm->SetminTPCncls(50);
 372     dtckm->SetRemoveKinks(kTRUE);
 373     dtckm->SetLabel(kFALSE);
 374     dtckm->SetMaxITSChiNdof(2.5);
 375     dtckm->SetMaxTPCChiNdof(3.0);
 376     dtckm->SetMaxImpactXY(3.0);
 377     dtckm->SetMaxImpactZ(3.0);
 378
 379     AliFemtoCutMonitorParticleYPt *cutPassYPtkm = new AliFemtoCutMonitorParticleYPt("cutPasskm", 0.493677);
 380     AliFemtoCutMonitorParticleYPt *cutFailYPtkm = new AliFemtoCutMonitorParticleYPt("cutFailkm", 0.493677);
 381     dtckm->AddCutMonitor(cutPassYPtkm, cutFailYPtkm);
 382
 383     AliFemtoCutMonitorParticlePtPDG *cutPassPidkm = new AliFemtoCutMonitorParticlePtPDG("cutPasskm", 0.493677);
 384     AliFemtoCutMonitorParticlePtPDG *cutFailPidkm = new AliFemtoCutMonitorParticlePtPDG("cutFailkm", 0.493677);
 385     dtckm->AddCutMonitor(cutPassPidkm, cutFailPidkm);
 386
 387     AliFemtoCutMonitorParticleMomRes *cutPassMRkm = new AliFemtoCutMonitorParticleMomRes("cutPasskm");
 388     AliFemtoCutMonitorParticleMomRes *cutFailMRkm = new AliFemtoCutMonitorParticleMomRes("cutFailkm");
 389     dtckm->AddCutMonitor(cutPassMRkm, cutFailMRkm);
 390
 391     AliFemtoCutMonitorParticleVertPos *cutPassVPkm = new AliFemtoCutMonitorParticleVertPos("cutPasskm");
 392     AliFemtoCutMonitorParticleVertPos *cutFailVPkm = new AliFemtoCutMonitorParticleVertPos("cutFailkm");
 393     dtckm->AddCutMonitor(cutPassVPkm, cutFailVPkm);
 394
 395     AliFemtoCutMonitorEventMult *cutPassEvMkm = new AliFemtoCutMonitorEventMult("cutPasskm");
 396     AliFemtoCutMonitorEventMult *cutFailEvMkm = new AliFemtoCutMonitorEventMult("cutFailkm");
 397     meckm->AddCutMonitor(cutPassEvMkm, cutFailEvMkm);
 398
 399     AliFemtoCutMonitorEventVertex *cutPassEvVkm = new AliFemtoCutMonitorEventVertex("cutPasskm");
 400     AliFemtoCutMonitorEventVertex *cutFailEvVkm = new AliFemtoCutMonitorEventVertex("cutFailkm");
 401     meckm->AddCutMonitor(cutPassEvVkm, cutFailEvVkm);
 402
 403     AliFemtoShareQualityTPCEntranceSepPairCut *sqpckm = new AliFemtoShareQualityTPCEntranceSepPairCut();
 404     sqpckm->SetShareQualityMax(0.0);
 405     sqpckm->SetShareFractionMax(0.02);
 406     sqpckm->SetRemoveSameLabel(kFALSE);
 407     sqpckm->SetTPCEntranceSepMinimum(3.0);
 408
 409     ankm->SetEventCut(meckm);
 410     ankm->SetFirstParticleCut(dtckm);
 411     ankm->SetSecondParticleCut(dtckm);
 412     ankm->SetPairCut(sqpckm);
 413
 414     AliFemtoQinvCorrFctn *cqinvkm= new AliFemtoQinvCorrFctn("qinvcf",100,0.0,1.0);
 415
 416
 417     //###
 418     ankm->AddCorrFctn(cqinvkm);
 419
 420     Manager->AddAnalysis(ankm);
 421
 422     // *** End Kaon-Kaon (negative) analysis
 423   }
 424
 425   if (runPositiveNegativeKaons) {
 426     // *** Begin Kaon+Kaon- analysis
 427     AliFemtoVertexMultAnalysis *ankpkm = new AliFemtoVertexMultAnalysis(nZVertexBins, minZVertex, maxZVertex, 1, 2, 20000);
 428     ankpkm->SetNumEventsToMix(5);
 429     ankpkm->SetMinSizePartCollection(2);
 430
 431     AliFemtoBasicEventCut* meckpkm = new AliFemtoBasicEventCut();
 432     meckpkm->SetEventMult(1,100000);
 433     meckpkm->SetVertZPos(-1000,1000);
 434
 435     AliFemtoESDTrackCut* dtckp = new AliFemtoESDTrackCut();
 436     dtckp->SetPidProbKaon(0.7,1.001);
 437     dtckp->SetPidProbMuon(0.0,0.5);
 438     dtckp->SetPidProbPion(0.0,0.5);
 439     dtckp->SetPidProbProton(0.0,0.5);
 440     dtckp->SetMomRangeTOFpidIs(0.6,10000.);
 441     dtckp->SetCharge(1.0);
 442     dtckp->SetMostProbableKaon();
 443     dtckp->SetPt(0.15,2.0);
 444     dtckp->SetMass(KaonMass);
 445     // Track quality cuts
 446     dtckp->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 447     dtckp->SetminTPCncls(50);
 448     dtckp->SetRemoveKinks(kTRUE);
 449     dtckp->SetLabel(kFALSE);
 450     dtckp->SetMaxITSChiNdof(2.5);
 451     dtckp->SetMaxTPCChiNdof(3.0);
 452     dtckp->SetMaxImpactXY(3.0);
 453     dtckp->SetMaxImpactZ(3.0);
 454
 455     AliFemtoESDTrackCut* dtckm = new AliFemtoESDTrackCut();
 456     dtckm->SetPidProbKaon(0.7,1.001);
 457     dtckm->SetPidProbMuon(0.0,0.5);
 458     dtckm->SetPidProbPion(0.0,0.5);
 459     dtckm->SetPidProbProton(0.0,0.5);
 460     dtckm->SetMomRangeTOFpidIs(0.6,10000.);
 461     dtckm->SetCharge(-1.0);
 462     dtckm->SetMostProbableKaon();
 463     dtckm->SetPt(0.15,2.0);
 464     dtckm->SetMass(KaonMass);
 465     // Track quality cuts
 466     dtckm->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
 467     dtckm->SetminTPCncls(50);
 468     dtckm->SetRemoveKinks(kTRUE);
 469     dtckm->SetLabel(kFALSE);
 470     dtckm->SetMaxITSChiNdof(2.5);
 471     dtckm->SetMaxTPCChiNdof(3.0);
 472     dtckm->SetMaxImpactXY(3.0);
 473     dtckm->SetMaxImpactZ(3.0);
 474
 475
 476     AliFemtoCutMonitorParticleYPt *cutPassYPtkp = new AliFemtoCutMonitorParticleYPt("cutPasskp", 0.493677);
 477     AliFemtoCutMonitorParticleYPt *cutFailYPtkp = new AliFemtoCutMonitorParticleYPt("cutFailkp", 0.493677);
 478     dtckp->AddCutMonitor(cutPassYPtkp, cutFailYPtkp);
 479
 480     AliFemtoCutMonitorParticleYPt *cutPassYPtkm = new AliFemtoCutMonitorParticleYPt("cutPasskm", 0.493677);
 481     AliFemtoCutMonitorParticleYPt *cutFailYPtkm = new AliFemtoCutMonitorParticleYPt("cutFailkm", 0.493677);
 482     dtckm->AddCutMonitor(cutPassYPtkm, cutFailYPtkm);
 483
 484     AliFemtoCutMonitorParticlePtPDG *cutPassPidkp = new AliFemtoCutMonitorParticlePtPDG("cutPasskp", 0.493677);
 485     AliFemtoCutMonitorParticlePtPDG *cutFailPidkp = new AliFemtoCutMonitorParticlePtPDG("cutFailkp", 0.493677);
 486     dtckp->AddCutMonitor(cutPassPidkp, cutFailPidkp);
 487
 488     AliFemtoCutMonitorParticlePtPDG *cutPassPidkm = new AliFemtoCutMonitorParticlePtPDG("cutPasskm", 0.493677);
 489     AliFemtoCutMonitorParticlePtPDG *cutFailPidkm = new AliFemtoCutMonitorParticlePtPDG("cutFailkm", 0.493677);
 490     dtckm->AddCutMonitor(cutPassPidkm, cutFailPidkm);
 491
 492     AliFemtoCutMonitorParticleMomRes *cutPassMRkp = new AliFemtoCutMonitorParticleMomRes("cutPasskp");
 493     AliFemtoCutMonitorParticleMomRes *cutFailMRkp = new AliFemtoCutMonitorParticleMomRes("cutFailkp");
 494     dtckp->AddCutMonitor(cutPassMRkp, cutFailMRkp);
 495
 496     AliFemtoCutMonitorParticleMomRes *cutPassMRkm = new AliFemtoCutMonitorParticleMomRes("cutPasskm");
 497     AliFemtoCutMonitorParticleMomRes *cutFailMRkm = new AliFemtoCutMonitorParticleMomRes("cutFailkm");
 498     dtckm->AddCutMonitor(cutPassMRkm, cutFailMRkm);
 499
 500     AliFemtoCutMonitorParticleVertPos *cutPassVPkp = new AliFemtoCutMonitorParticleVertPos("cutPasskp");
 501     AliFemtoCutMonitorParticleVertPos *cutFailVPkp = new AliFemtoCutMonitorParticleVertPos("cutFailkp");
 502     dtckp->AddCutMonitor(cutPassVPkp, cutFailVPkp);
 503
 504     AliFemtoCutMonitorParticleVertPos *cutPassVPkm = new AliFemtoCutMonitorParticleVertPos("cutPasskm");
 505     AliFemtoCutMonitorParticleVertPos *cutFailVPkm = new AliFemtoCutMonitorParticleVertPos("cutFailkm");
 506     dtckm->AddCutMonitor(cutPassVPkm, cutFailVPkm);
 507
 508     AliFemtoCutMonitorEventMult *cutPassEvMkpkm = new AliFemtoCutMonitorEventMult("cutPasskpkm");
 509     AliFemtoCutMonitorEventMult *cutFailEvMkpkm = new AliFemtoCutMonitorEventMult("cutFailkpkm");
 510     meckpkm->AddCutMonitor(cutPassEvMkpkm, cutFailEvMkpkm);
 511
 512     AliFemtoCutMonitorEventVertex *cutPassEvVkpkm = new AliFemtoCutMonitorEventVertex("cutPasskpkm");
 513     AliFemtoCutMonitorEventVertex *cutFailEvVkpkm = new AliFemtoCutMonitorEventVertex("cutFailkpkm");
 514     meckpkm->AddCutMonitor(cutPassEvVkpkm, cutFailEvVkpkm);
 515
 516     AliFemtoShareQualityTPCEntranceSepPairCut *sqpckpkm = new AliFemtoShareQualityTPCEntranceSepPairCut();
 517     sqpckpkm->SetShareQualityMax(0.0);
 518     sqpckpkm->SetShareFractionMax(0.02);
 519     sqpckpkm->SetRemoveSameLabel(kFALSE);
 520     sqpckpkm->SetTPCEntranceSepMinimum(3.0);
 521
 522     ankpkm->SetEventCut(meckpkm);
 523     ankpkm->SetFirstParticleCut(dtckp);
 524     ankpkm->SetSecondParticleCut(dtckm);
 525     ankpkm->SetPairCut(sqpckpkm);
 526
 527     AliFemtoQinvCorrFctn *cqinvkpkm= new AliFemtoQinvCorrFctn("qinvcf",100,0.0,1.0);
 528
 529     //###
 530     ankpkm->AddCorrFctn(cqinvkpkm);
 531
 532     Manager->AddAnalysis(ankpkm);
 533
 534     // *** End Kaon+Kaon-  analysis
 535   }
 536
 537   return Manager;
 538 }
 539