[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / KaonTrain / PbPb2011 / MagFieldMinus / ConfigFemtoAnalysis.C


/*********************************************************************
 *                                                                   *
 * ConfigFemtoAnalysis.C - configuration macro for the femtoscopic   *
 * analysis, meant as a QA process for two-particle effects          *
 *                                                                   *
 * Author: Adam Kisiel (Adam.Kisiel@cern.ch)                         *
 *                                                                   *
 *********************************************************************/

#if !defined(__CINT__) || defined(__MAKECINT_)
#include "AliFemtoManager.h"
#include "AliFemtoEventReaderESDChain.h"
#include "AliFemtoEventReaderESDChainKine.h"
#include "AliFemtoEventReaderAODChain.h"
#include "AliFemtoSimpleAnalysis.h"
#include "AliFemtoBasicEventCut.h"
#include "AliFemtoESDTrackCut.h"
#include "AliFemtoCorrFctn.h"
#include "AliFemtoCutMonitorParticleYPt.h"
#include "AliFemtoCutMonitorParticleVertPos.h"
#include "AliFemtoCutMonitorParticleMomRes.h"
#include "AliFemtoCutMonitorParticlePID.h"
#include "AliFemtoCutMonitorEventMult.h"
#include "AliFemtoCutMonitorEventVertex.h"
#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
#include "AliFemtoPairCutAntiGamma.h"
#include "AliFemtoPairCutRadialDistance.h"
#include "AliFemtoQinvCorrFctn.h"
#include "AliFemtoShareQualityCorrFctn.h"
#include "AliFemtoTPCInnerCorrFctn.h"
#include "AliFemtoVertexMultAnalysis.h"
#include "AliFemtoCorrFctn3DSpherical.h"
#include "AliFemtoChi2CorrFctn.h"
#include "AliFemtoCorrFctnTPCNcls.h"
#include "AliFemtoBPLCMS3DCorrFctn.h"
#include "AliFemtoCorrFctn3DLCMSSym.h"
#include "AliFemtoModelBPLCMSCorrFctn.h"
#include "AliFemtoModelCorrFctn3DSpherical.h"
#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
#include "AliFemtoModelGausRinvFreezeOutGenerator.h"
#include "AliFemtoModelManager.h"
#include "AliFemtoModelWeightGeneratorBasic.h"
#include "AliFemtoModelWeightGeneratorLednicky.h"
#include "AliFemtoCorrFctnDirectYlm.h"
#include "AliFemtoModelCorrFctnDirectYlm.h"
#include "AliFemtoModelCorrFctnSource.h"
#include "AliFemtoKTPairCut.h"
#endif

//________________________________________________________________________
AliFemtoManager* ConfigFemtoAnalysis() {

  double PionMass = 0.13956995;
  double KaonMass = 0.493677;
	
  // Switches for QA analyses
  // int runmults[10] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 1};
  // int multbins[11] = {30, 70, 150, 310, 590, 990, 1570, 2370, 2370, 2370, 6500};

  int runmults[10] = {1, 1, 1, 0, 0, 0, 0, 0, 0, 0};
  int multbins[11] = {0, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};

  int runch[2] = {1, 1};
  const char *chrgs[2] = { "Kp", "Km"};
  
  
  int runktdep = 1;
  double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};

  int run3d = 1;
  int runshlcms = 0;

  int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
  int isrealdata = 1;

  //   AliFemtoEventReaderESDChainKine* Reader=new AliFemtoEventReaderESDChainKine();
  //   Reader->SetConstrained(true);
  //   Reader->SetUseTPCOnly(false);

  double shqmax;
  double shqmaxSH;
  int nbinssh = 100;

//ml  if (runshlcms) shqmax = 0.25;
//  else shqmax = 0.9;
  

  if (runshlcms) shqmaxSH = 0.25;
  shqmax = 0.9;
  

  // AliFemtoEventReaderESDChain* Reader=new AliFemtoEventReaderESDChain();
  // Reader->SetConstrained(true);
  // Reader->SetUseTPCOnly(false);
  // Reader->SetReadTPCInner(false);
  // Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);

  // if (runtype == 0)
  //   Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kGlobal);
  // else if (runtype == 1)
  //   Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kITSOnly);
  // else if (runtype == 2)
  //   Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kTPCOnly);
  // if (isrealdata)
  //   Reader->SetUsePhysicsSelection(kTRUE);
  // else
  //   Reader->SetUsePhysicsSelection(kFALSE);

  // Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);

  AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
    Reader->SetFilterBit(7);
    Reader->SetCentralityPreSelection(0, 200);
    
  AliFemtoManager* Manager=new AliFemtoManager();
  Manager->SetEventReader(Reader);

  AliFemtoVertexMultAnalysis    *anetaphitpc[20];
  AliFemtoBasicEventCut         *mecetaphitpc[20];
  AliFemtoCutMonitorEventMult   *cutPassEvMetaphitpc[20];
  AliFemtoCutMonitorEventMult   *cutFailEvMetaphitpc[20];
  AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[20];
  AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[20];
  AliFemtoESDTrackCut           *dtc1etaphitpc[20];
  AliFemtoESDTrackCut           *dtc2etaphitpc[20];
  AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[20];
  AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[20];
  AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[20];
  AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[20];
  AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[20];
  AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[20];
  AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[20];
  AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[20];
//AliFemtoPairCutAntiGamma      *sqpcetaphitpc[20];
//    AliFemtoShareQualityTPCEntranceSepPairCut      *sqpcetaphitpc[20];
   AliFemtoPairCutRadialDistance      *sqpcetaphitpc[20];
  AliFemtoCorrFctnDirectYlm     *cylmetaphitpc[20];
  AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[20];
  AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[20];
  AliFemtoKTPairCut             *ktpcuts[20*7];
  AliFemtoCorrFctnDirectYlm     *cylmkttpc[20*7];
  AliFemtoQinvCorrFctn          *cqinvkttpc[20*7];
  AliFemtoCorrFctn3DLCMSSym     *cq3dlcmskttpc[20*7];
  AliFemtoCorrFctnTPCNcls       *cqinvnclstpc[20];
  AliFemtoShareQualityCorrFctn  *cqinvsqtpc[20*10];
  AliFemtoChi2CorrFctn          *cqinvchi2tpc[20];
  AliFemtoTPCInnerCorrFctn      *cqinvinnertpc[20*10];

  // *** Third QA task - HBT analysis with all pair cuts off, TPC only ***
  // *** Begin Kaon-Kaon (positive) analysis ***
  int aniter = 0;

  for (int imult=0; imult<10; imult++) {
    if (runmults[imult]) {
      for (int ichg=0; ichg<2; ichg++) {
	if (runch[ichg]) {
	  aniter = ichg*10+imult;

	  anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(4, -8.0, 8.0, 5, multbins[imult], multbins[imult+1]);
	  anetaphitpc[aniter]->SetNumEventsToMix(3);
	  anetaphitpc[aniter]->SetMinSizePartCollection(1);

	  mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
	  mecetaphitpc[aniter]->SetEventMult(0,100000);
	  mecetaphitpc[aniter]->SetVertZPos(-8,8);
	  if (isrealdata)
	     mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
	  
	  //    mecetaphitpc->SetAcceptBadVertex(kTRUE);
	  
	  cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	  cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	  mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
	  
	  cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	  cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	  mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
	  
	  dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
	  //     dtc1etaphitpc[aniter]->SetPidProbPion(0.0,1.001);
	  //     dtc1etaphitpc[aniter]->SetPidProbMuon(0.0,1.0);
	  //     dtc1etaphitpc[aniter]->SetPidProbKaon(0.0,1.0);
	  //     dtc1etaphitpc[aniter]->SetPidProbProton(0.0,1.0);
	  if (ichg == 0)
	    dtc1etaphitpc[aniter]->SetCharge(1.0);
	  else if (ichg == 1)
	    dtc1etaphitpc[aniter]->SetCharge(-1.0);
	    
	  dtc1etaphitpc[aniter]->SetPt(0.14,1.5);
	  //	  dtc1etaphitpc[aniter]->SetEta(-1.2,1.2);
	  dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
	  // 	//    dtc1etaphitpc[aniter]->SetEta(-0.5,0.5);
	  dtc1etaphitpc[aniter]->SetMass(KaonMass);
	  //	  dtc1etaphitpc[aniter]->SetminTPCncls(80);
		  
	  dtc1etaphitpc[aniter]->SetMostProbableKaon();
	  
	  //	  dtc1etaphitpc[aniter]->SetMostProbablePion();
	  // 	// Track quality cuts
	  if (runtype == 0) {
	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
	    //	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	    //    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	    dtc1etaphitpc[aniter]->SetminTPCncls(80);
	    dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	    dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	    //    dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	    dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
	    dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
	    //            dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15);
	    //      dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	  }
	  else if (runtype == 1) {
	    //      dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
	    //    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	    //	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit|AliESDtrack::kITSpureSA);
	    //      dtc1etaphitpc[aniter]->SetminTPCncls(70);
	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	    dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	    dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	    //    dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	    //      dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(6.0);
	    dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(0.25);
	    //      dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	  }
	  else if (runtype == 2) {
	    //	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
	    //	    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	    //    dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	    dtc1etaphitpc[aniter]->SetminTPCncls(80); //was "0"
	    dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	    dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	    //    dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	    dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
	    dtc1etaphitpc[aniter]->SetMaxImpactXY(0.20); //2.4
	    //            dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15);  //3.0
	    //      dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	  }

	  
	  cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%i", chrgs[ichg], imult), 0.493677);
	  cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%i", chrgs[ichg], imult), 0.493677);
	  dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
	  
	  cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass1%stpcM%i", chrgs[ichg], imult),0);
	  cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%i", chrgs[ichg], imult),0);
	  dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);

	  
	 // sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
	//  sqpcetaphitpc[aniter] = new AliFemtoShareQualityTPCEntranceSepPairCut();
	  
          if (ichg < 2) {
    	  sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
	  if (runtype == 0) {
	    sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
	    sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
	    sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	    // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
	    // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
	    	    //ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
//ml	    sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
//ml	    sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);

        
	  sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
          sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
          sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
          sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
	  
	 
	  }
	  else if (runtype == 1) {
	    sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
	    sqpcetaphitpc[aniter]->SetShareFractionMax(1.05);
	    sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	    // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
	    // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
	    	    //ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(5.0);
//	    sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.03);
//	    sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);

//  void SetPhiStarDifferenceMinimum(double dtpc);
//  void SetEtaDifferenceMinimum(double etpc);
//  void SetMinimumRadius(double minrad);
//  void SetMagneticFieldSign(int magsign);

         sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
         sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
         sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
         sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
     //sqpcetaphitpc[aniter]->SetMagneticFieldSign(1.0);


	  }
	  else if (runtype == 2) {
            sqpcetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
	    sqpcetaphitpc[aniter]->SetShareQualityMax(0.0);
	    sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
	    sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	    // sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
	    // sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
	    	//ml    sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.0);
//ml	    sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.045);
//ml	    sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.016);

     	   sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
           sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
           sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
           sqpcetaphitpc[aniter]->SetMagneticFieldSign(-1);
	 

	  }
         }        

	  
	  anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
	  anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
	  anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
	  anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
	  
	 // 	  cylmetaphitpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),3,nbinssh,0.0,shqmaxSH,runshlcms);
	  ///	  anetaphitpc[aniter]->AddCorrFctn(cylmetaphitpc[aniter]);
	  
	  // cqinvnclstpc[aniter] = new AliFemtoCorrFctnTPCNcls(Form("cqinvncls%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	  // anetaphitpc[aniter]->AddCorrFctn(cqinvnclstpc[aniter]);

	  //	  cqinvchi2tpc[aniter] = new AliFemtoChi2CorrFctn(Form("cqinvchi2%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	  //	  anetaphitpc[aniter]->AddCorrFctn(cqinvchi2tpc[aniter]);

	  if (runktdep) {
	    int ktm;
	    for (int ikt=0; ikt<7; ikt++) {
	      ktm = aniter*7 + ikt;
	      ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
	      
////////  	      cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,
///////							     nbinssh, 0.0,
//////							     (imult>6)?shqmaxSH*2.5:shqmaxSH,
/////							     runshlcms);
//////  	      cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
/////   	      anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[ktm]);
	      
              cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
//	      cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,0.5);
	      cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	      anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);

	      cqinvsqtpc[ktm] = new AliFemtoShareQualityCorrFctn(Form("cqinvsq%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
	      cqinvsqtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	      anetaphitpc[aniter]->AddCorrFctn(cqinvsqtpc[ktm]);

	      cqinvinnertpc[ktm] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
	      cqinvinnertpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	      cqinvinnertpc[ktm]->SetRadius(1.2);
	      anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[ktm]);

	      if (run3d) {
		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);
//		cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),50,0.5);
		cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
		anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
	      }
	    }
	  }
	  
	  // cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),24, 24);
	  // anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
	  
	  Manager->AddAnalysis(anetaphitpc[aniter]);	
	}
      }
    }
  }
  // *** End Kaon-Kaon (positive) analysis

  return Manager;
}
Commit	Line	Data
3a32e4dc	1
	2	/*********************************************************************
	3	* *
	4	* ConfigFemtoAnalysis.C - configuration macro for the femtoscopic *
	5	* analysis, meant as a QA process for two-particle effects *
	6	* *
	7	* Author: Adam Kisiel (Adam.Kisiel@cern.ch) *
	8	* *
	9	*********************************************************************/
	10
	11	#if !defined(__CINT__) \|\| defined(__MAKECINT_)
	12	#include "AliFemtoManager.h"
	13	#include "AliFemtoEventReaderESDChain.h"
	14	#include "AliFemtoEventReaderESDChainKine.h"
	15	#include "AliFemtoEventReaderAODChain.h"
	16	#include "AliFemtoSimpleAnalysis.h"
	17	#include "AliFemtoBasicEventCut.h"
	18	#include "AliFemtoESDTrackCut.h"
	19	#include "AliFemtoCorrFctn.h"
	20	#include "AliFemtoCutMonitorParticleYPt.h"
	21	#include "AliFemtoCutMonitorParticleVertPos.h"
	22	#include "AliFemtoCutMonitorParticleMomRes.h"
	23	#include "AliFemtoCutMonitorParticlePID.h"
	24	#include "AliFemtoCutMonitorEventMult.h"
	25	#include "AliFemtoCutMonitorEventVertex.h"
	26	#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
	27	#include "AliFemtoPairCutAntiGamma.h"
	28	#include "AliFemtoPairCutRadialDistance.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 "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 "AliFemtoKTPairCut.h"
	49	#endif
	50
	51	//________________________________________________________________________
	52	AliFemtoManager* ConfigFemtoAnalysis() {
	53
	54	double PionMass = 0.13956995;
	55	double KaonMass = 0.493677;
	56
	57	// Switches for QA analyses
	58	// int runmults[10] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 1};
	59	// int multbins[11] = {30, 70, 150, 310, 590, 990, 1570, 2370, 2370, 2370, 6500};
	60
bb3ceab4	61	int runmults[10] = {1, 1, 1, 0, 0, 0, 0, 0, 0, 0};
3a32e4dc	62	int multbins[11] = {0, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
	63
	64	int runch[2] = {1, 1};
	65	const char *chrgs[2] = { "Kp", "Km"};
	66
	67
	68	int runktdep = 1;
	69	double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};
	70
	71	int run3d = 1;
	72	int runshlcms = 0;
	73
	74	int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
	75	int isrealdata = 1;
	76
	77	// AliFemtoEventReaderESDChainKine* Reader=new AliFemtoEventReaderESDChainKine();
	78	// Reader->SetConstrained(true);
	79	// Reader->SetUseTPCOnly(false);
	80
	81	double shqmax;
	82	double shqmaxSH;
	83	int nbinssh = 100;
	84
	85	//ml if (runshlcms) shqmax = 0.25;
	86	// else shqmax = 0.9;
	87
	88
	89	if (runshlcms) shqmaxSH = 0.25;
bb3ceab4	90	shqmax = 0.9;
3a32e4dc	91
	92
	93	// AliFemtoEventReaderESDChain* Reader=new AliFemtoEventReaderESDChain();
	94	// Reader->SetConstrained(true);
	95	// Reader->SetUseTPCOnly(false);
	96	// Reader->SetReadTPCInner(false);
	97	// Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);
	98
	99	// if (runtype == 0)
	100	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kGlobal);
	101	// else if (runtype == 1)
	102	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kITSOnly);
	103	// else if (runtype == 2)
	104	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kTPCOnly);
	105	// if (isrealdata)
	106	// Reader->SetUsePhysicsSelection(kTRUE);
	107	// else
	108	// Reader->SetUsePhysicsSelection(kFALSE);
	109
	110	// Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);
	111
	112	AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
	113	Reader->SetFilterBit(7);
bb3ceab4	114	Reader->SetCentralityPreSelection(0, 200);
3a32e4dc	115
	116	AliFemtoManager* Manager=new AliFemtoManager();
	117	Manager->SetEventReader(Reader);
	118
	119	AliFemtoVertexMultAnalysis *anetaphitpc[20];
	120	AliFemtoBasicEventCut *mecetaphitpc[20];
	121	AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[20];
	122	AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[20];
	123	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[20];
	124	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[20];
	125	AliFemtoESDTrackCut *dtc1etaphitpc[20];
	126	AliFemtoESDTrackCut *dtc2etaphitpc[20];
	127	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[20];
	128	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[20];
	129	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[20];
	130	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[20];
	131	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[20];
	132	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[20];
	133	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[20];
	134	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[20];
	135	//AliFemtoPairCutAntiGamma *sqpcetaphitpc[20];
	136	// AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpc[20];
	137	AliFemtoPairCutRadialDistance *sqpcetaphitpc[20];
	138	AliFemtoCorrFctnDirectYlm *cylmetaphitpc[20];
	139	AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[20];
	140	AliFemtoChi2CorrFctn *cchiqinvetaphitpc[20];
	141	AliFemtoKTPairCut ktpcuts[207];
	142	AliFemtoCorrFctnDirectYlm cylmkttpc[207];
	143	AliFemtoQinvCorrFctn cqinvkttpc[207];
	144	AliFemtoCorrFctn3DLCMSSym cq3dlcmskttpc[207];
	145	AliFemtoCorrFctnTPCNcls *cqinvnclstpc[20];
	146	AliFemtoShareQualityCorrFctn cqinvsqtpc[2010];
	147	AliFemtoChi2CorrFctn *cqinvchi2tpc[20];
	148	AliFemtoTPCInnerCorrFctn cqinvinnertpc[2010];
	149
	150	// * Third QA task - HBT analysis with all pair cuts off, TPC only *
	151	// * Begin Kaon-Kaon (positive) analysis *
	152	int aniter = 0;
	153
	154	for (int imult=0; imult<10; imult++) {
	155	if (runmults[imult]) {
	156	for (int ichg=0; ichg<2; ichg++) {
	157	if (runch[ichg]) {
	158	aniter = ichg*10+imult;
	159
	160	anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(4, -8.0, 8.0, 5, multbins[imult], multbins[imult+1]);
bb3ceab4	161	anetaphitpc[aniter]->SetNumEventsToMix(3);
3a32e4dc	162	anetaphitpc[aniter]->SetMinSizePartCollection(1);
	163
	164	mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
	165	mecetaphitpc[aniter]->SetEventMult(0,100000);
	166	mecetaphitpc[aniter]->SetVertZPos(-8,8);
	167	if (isrealdata)
	168	mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
	169
	170	// mecetaphitpc->SetAcceptBadVertex(kTRUE);
	171
	172	cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	173	cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	174	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
	175
bb3ceab4	176	cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	177	cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	178	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
3a32e4dc	179
	180	dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
	181	// dtc1etaphitpc[aniter]->SetPidProbPion(0.0,1.001);
	182	// dtc1etaphitpc[aniter]->SetPidProbMuon(0.0,1.0);
	183	// dtc1etaphitpc[aniter]->SetPidProbKaon(0.0,1.0);
	184	// dtc1etaphitpc[aniter]->SetPidProbProton(0.0,1.0);
	185	if (ichg == 0)
	186	dtc1etaphitpc[aniter]->SetCharge(1.0);
	187	else if (ichg == 1)
	188	dtc1etaphitpc[aniter]->SetCharge(-1.0);
	189
bb3ceab4	190	dtc1etaphitpc[aniter]->SetPt(0.14,1.5);
3a32e4dc	191	// dtc1etaphitpc[aniter]->SetEta(-1.2,1.2);
	192	dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
	193	// // dtc1etaphitpc[aniter]->SetEta(-0.5,0.5);
	194	dtc1etaphitpc[aniter]->SetMass(KaonMass);
	195	// dtc1etaphitpc[aniter]->SetminTPCncls(80);
	196
	197	dtc1etaphitpc[aniter]->SetMostProbableKaon();
	198
	199	// dtc1etaphitpc[aniter]->SetMostProbablePion();
	200	// // Track quality cuts
	201	if (runtype == 0) {
	202	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
	203	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	204	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	205	dtc1etaphitpc[aniter]->SetminTPCncls(80);
	206	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	207	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	208	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	209	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
	210	dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
	211	// dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	212	dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15);
	213	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	214	}
	215	else if (runtype == 1) {
	216	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
	217	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	218	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit\|AliESDtrack::kITSpureSA);
	219	// dtc1etaphitpc[aniter]->SetminTPCncls(70);
	220	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	221	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	222	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	223	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	224	// dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(6.0);
	225	dtc1etaphitpc[aniter]->SetMaxImpactXY(0.2);
	226	dtc1etaphitpc[aniter]->SetMaxImpactZ(0.25);
	227	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	228	}
	229	else if (runtype == 2) {
	230	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
	231	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
	232	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
	233	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
	234	dtc1etaphitpc[aniter]->SetminTPCncls(80); //was "0"
	235	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	236	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	237	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
	238	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
	239	dtc1etaphitpc[aniter]->SetMaxImpactXY(0.20); //2.4
	240	// dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	241	dtc1etaphitpc[aniter]->SetMaxImpactZ(0.15); //3.0
	242	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
	243	}
	244
	245
bb3ceab4	246
3a32e4dc	247	cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%i", chrgs[ichg], imult), 0.493677);
	248	cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%i", chrgs[ichg], imult), 0.493677);
	249	dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
	250
	251	cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass1%stpcM%i", chrgs[ichg], imult),0);
	252	cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%i", chrgs[ichg], imult),0);
	253	dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);
bb3ceab4	254
3a32e4dc	255
	256	// sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
	257	// sqpcetaphitpc[aniter] = new AliFemtoShareQualityTPCEntranceSepPairCut();
	258
	259	if (ichg < 2) {
	260	sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
	261	if (runtype == 0) {
	262	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
	263	sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
	264	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	265	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
	266	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
	267	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
	268	//ml sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
	269	//ml sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
	270
	271
	272	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
	273	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
	274	sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
bb3ceab4	275	sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
3a32e4dc	276
	277
	278
	279	}
	280	else if (runtype == 1) {
	281	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
	282	sqpcetaphitpc[aniter]->SetShareFractionMax(1.05);
	283	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	284	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
	285	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
	286	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(5.0);
	287	// sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.03);
	288	// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
	289
	290	// void SetPhiStarDifferenceMinimum(double dtpc);
	291	// void SetEtaDifferenceMinimum(double etpc);
	292	// void SetMinimumRadius(double minrad);
	293	// void SetMagneticFieldSign(int magsign);
	294
	295	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
	296	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
	297	sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
bb3ceab4	298	sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
3a32e4dc	299	//sqpcetaphitpc[aniter]->SetMagneticFieldSign(1.0);
	300
	301
	302	}
	303	else if (runtype == 2) {
	304	sqpcetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
bb3ceab4	305	sqpcetaphitpc[aniter]->SetShareQualityMax(0.0);
3a32e4dc	306	sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
	307	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	308	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
	309	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
	310	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.0);
	311	//ml sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.045);
	312	//ml sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.016);
	313
	314	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.017);
	315	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.015);
	316	sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
	317	sqpcetaphitpc[aniter]->SetMagneticFieldSign(-1);
	318
	319
	320	}
	321	}
	322
	323
	324	anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
	325	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
	326	anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
	327	anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
	328
	329	// cylmetaphitpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),3,nbinssh,0.0,shqmaxSH,runshlcms);
	330	/// anetaphitpc[aniter]->AddCorrFctn(cylmetaphitpc[aniter]);
	331
	332	// cqinvnclstpc[aniter] = new AliFemtoCorrFctnTPCNcls(Form("cqinvncls%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	333	// anetaphitpc[aniter]->AddCorrFctn(cqinvnclstpc[aniter]);
	334
	335	// cqinvchi2tpc[aniter] = new AliFemtoChi2CorrFctn(Form("cqinvchi2%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	336	// anetaphitpc[aniter]->AddCorrFctn(cqinvchi2tpc[aniter]);
	337
	338	if (runktdep) {
	339	int ktm;
	340	for (int ikt=0; ikt<7; ikt++) {
	341	ktm = aniter*7 + ikt;
	342	ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
	343
	344	//////// cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,
	345	/////// nbinssh, 0.0,
	346	////// (imult>6)?shqmaxSH*2.5:shqmaxSH,
	347	///// runshlcms);
	348	////// cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	349	///// anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[ktm]);
	350
bb3ceab4	351	cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
bb3ceab4	352	// cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,0.5);
3a32e4dc	353	cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	354	anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
	355
bb3ceab4	356	cqinvsqtpc[ktm] = new AliFemtoShareQualityCorrFctn(Form("cqinvsq%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
	357	cqinvsqtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	358	anetaphitpc[aniter]->AddCorrFctn(cqinvsqtpc[ktm]);
3a32e4dc	359
bb3ceab4	360	cqinvinnertpc[ktm] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
	361	cqinvinnertpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	362	cqinvinnertpc[ktm]->SetRadius(1.2);
	363	anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[ktm]);
3a32e4dc	364
3a32e4dc	365	if (run3d) {
bb3ceab4	366	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);
bb3ceab4	367	// cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),50,0.5);
3a32e4dc	368	cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	369	anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
	370	}
	371	}
	372	}
	373
	374	// cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),24, 24);
	375	// anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
	376
	377	Manager->AddAnalysis(anetaphitpc[aniter]);
	378	}
	379	}
	380	}
	381	}
	382	// *** End Kaon-Kaon (positive) analysis
	383
	384	return Manager;
	385	}
	386