[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / KchHBT_Eta0015Phi0017_DCAxy2.4_DCAz3_ITSTPC_Update / 0010 / 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 "AliFemtoKKTrackCut.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 "AliFemtoPairCutRadialDistanceKK.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;
  const int cMu=4;
  const int cKt=8;

  //-------Single track cuts------------------------------------------------->
  double DCAxy=2.4;//cm // our standard is 0.20 cm; super narrow was 0.015cm
  double DCAz =3.0;//cm // our standard is 0.15 cm;
  //-------Single track cuts-------------------------------------------------<
  //=======Double track cuts=================================================>
  //Dhevan's : PhiStarDifferenceMinimum=0.06; EtaDifferenceMinimum=0.02;
  double PhiStarDifferenceMinimum=0.017; //[radian]
  double EtaDifferenceMinimum=0.015; //[radian]
  //=======Double track cuts=================================================<

  // 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[cMu] = {1, 0, 0, 0};
  //int multbins[11] = {0, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
  int multbins[cMu+1] = {0, 100, 300, 500, 900};

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

  int run3d = 0;
  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, 100);
    /*On 02/26/2014 06:09 PM, Maciej Szymański wrote:
> Dear Konstantin,
>
> as I remember correctly, you were using the DCA cut on the global tracks
> corrsponding to chosen TPC-only tracks (SetDCAglobalTrack(kTRUE), right?
> The update I made regards cutting on the DCA on the TPC-only track, so
> you need to change in your config to: SetDCAglobalTrack(kFALSE).*/
    //Reader->SetDCAglobalTrack(kTRUE);//option the DCA information from global tracks (ITS+TPC)
    Reader->SetDCAglobalTrack(kFALSE);
    
  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];
  AliFemtoKKTrackCut           *dtc1etaphitpc[20];
  AliFemtoKKTrackCut           *dtc2etaphitpc[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];//AliFemto dphi* cut
  AliFemtoPairCutRadialDistanceKK      *sqpcetaphitpc[20];//Dhevan's dphi* cut
  AliFemtoCorrFctnDirectYlm     *cylmetaphitpc[20];
  AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[20*10];//20->20*10 due to kT
  AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[20];
  AliFemtoKTPairCut             *ktpcuts[20*8];
  AliFemtoCorrFctnDirectYlm     *cylmkttpc[20*8];
  AliFemtoQinvCorrFctn          *cqinvkttpc[20*8];
  AliFemtoCorrFctn3DLCMSSym     *cq3dlcmskttpc[20*8];
  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<cMu/*4*/; imult++) {
    if (runmults[imult]) {
      for (int ichg=0; ichg<2; ichg++) {
	if (runch[ichg]) {
	  aniter = ichg*5+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.0,8.0);
	  /* //was in aliroot 5.03.76
	  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 AliFemtoKKTrackCut();
//	  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(PionMass);
	  dtc1etaphitpc[aniter]->SetMass(KaonMass);
	  
	  
	  ////	  dtc1etaphitpc[aniter]->SetminTPCncls(80);
		  
	  //--- Choose Kaon as Most Probable (switch on all cuts: TPC, TOF)---
	  dtc1etaphitpc[aniter]->SetMostProbableKaon();  
	  //------------------- November 2013 -----------------------------------< 
	  // new cuts to remove electron (do not take into analysis if 400<p<500) 
	  //remove 0.0 to include electrons
	  dtc1etaphitpc[aniter]->SetNsigmaTPCle250(1.5);
	  dtc1etaphitpc[aniter]->SetNsigmaTPC250_400(2.0);
	  dtc1etaphitpc[aniter]->SetNsigmaTPC400_450(2.0);//0.0);
	  dtc1etaphitpc[aniter]->SetNsigmaTPC450_500(2.0);//0.0);
	  dtc1etaphitpc[aniter]->SetNsigmaTPCge500(3.0);    
	  // new cuts are stronger, better separation of pion in TOF 
	  // when momentum is greater then 800 MeV/c
	  dtc1etaphitpc[aniter]->SetNsigmaTOF500_800(2.0);
	  dtc1etaphitpc[aniter]->SetNsigmaTOF800_1000(1.5);
	  dtc1etaphitpc[aniter]->SetNsigmaTOFge1000(1.0);
	  //------------------- November 2013 ----------------------------------->
	  
////	  	  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(DCAxy);
	    //Poland: dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz);
	    //      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(DCAxy);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz);
	    //      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(DCAxy);
	    //dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
	    dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz);  //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),1);
	  cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%i", chrgs[ichg], imult),1);
	  dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);

	  
	 // sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
	//  sqpcetaphitpc[aniter] = new AliFemtoShareQualityTPCEntranceSepPairCut();
	  
          if (ichg < 2) {
	    //sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();//AliFemto dphi* cut
    	  sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistanceKK();  //Dhevan's dphi* cut
	  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);

        
	   //--------- km:  eta-phi* Dhevan's custs ----------->>>>
     	    sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
            sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
           //sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
	   //--------- km:  eta-phi* Dhevan's custs -----------<<<

          //////////////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);

	   //--------- km:  eta-phi* Dhevan's custs ----------->>>>
     	    sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
            sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
           //sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
	   //--------- km:  eta-phi* Dhevan's custs -----------<<<

         /////////sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
     //sqpcetaphitpc[aniter]->SetMagneticFieldSign(1.0);


	  }
	  else if (runtype == 2) {
            sqpcetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
	    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(0.0);
//ml	    sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.045);
//ml	    sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.016);

	   //--------- km:  eta-phi* Dhevan's custs ----------->>>>
     	    sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
            sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
           //sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
	   //--------- km:  eta-phi* Dhevan's custs -----------<<<

           ////////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<cKt/*8*/; ikt++) {
	      ktm = aniter*cKt/*8*/ + 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]);

//---- Correlation Function vs Delta_Eta and Delta_Phi (not Phi*)---->>>
	      cdedpetaphi[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ikT%i", chrgs[ichg], imult, ikt),100,100);
	      anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[ktm]);
//---- Correlation Function vs Delta_Eta and Delta_Phi (not Phi*)----<<<

	      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
0b2b8550	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 "AliFemtoKKTrackCut.h"
	20	#include "AliFemtoCorrFctn.h"
	21	#include "AliFemtoCutMonitorParticleYPt.h"
	22	#include "AliFemtoCutMonitorParticleVertPos.h"
	23	#include "AliFemtoCutMonitorParticleMomRes.h"
	24	#include "AliFemtoCutMonitorParticlePID.h"
	25	#include "AliFemtoCutMonitorEventMult.h"
	26	#include "AliFemtoCutMonitorEventVertex.h"
	27	#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
	28	#include "AliFemtoPairCutAntiGamma.h"
	29	#include "AliFemtoPairCutRadialDistance.h"
	30	#include "AliFemtoPairCutRadialDistanceKK.h"
	31	#include "AliFemtoQinvCorrFctn.h"
	32	#include "AliFemtoShareQualityCorrFctn.h"
	33	#include "AliFemtoTPCInnerCorrFctn.h"
	34	#include "AliFemtoVertexMultAnalysis.h"
	35	#include "AliFemtoCorrFctn3DSpherical.h"
	36	#include "AliFemtoChi2CorrFctn.h"
	37	#include "AliFemtoCorrFctnTPCNcls.h"
	38	#include "AliFemtoBPLCMS3DCorrFctn.h"
	39	#include "AliFemtoCorrFctn3DLCMSSym.h"
	40	#include "AliFemtoModelBPLCMSCorrFctn.h"
	41	#include "AliFemtoModelCorrFctn3DSpherical.h"
	42	#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
	43	#include "AliFemtoModelGausRinvFreezeOutGenerator.h"
	44	#include "AliFemtoModelManager.h"
	45	#include "AliFemtoModelWeightGeneratorBasic.h"
	46	#include "AliFemtoModelWeightGeneratorLednicky.h"
	47	#include "AliFemtoCorrFctnDirectYlm.h"
	48	#include "AliFemtoModelCorrFctnDirectYlm.h"
	49	#include "AliFemtoModelCorrFctnSource.h"
	50	#include "AliFemtoKTPairCut.h"
	51	#endif
	52
	53	//________________________________________________________________________
	54	AliFemtoManager* ConfigFemtoAnalysis() {
	55
	56	double PionMass = 0.13956995;
	57	double KaonMass = 0.493677;
	58	const int cMu=4;
	59	const int cKt=8;
	60
	61	//-------Single track cuts------------------------------------------------->
	62	double DCAxy=2.4;//cm // our standard is 0.20 cm; super narrow was 0.015cm
	63	double DCAz =3.0;//cm // our standard is 0.15 cm;
	64	//-------Single track cuts-------------------------------------------------<
65	//=======Double track cuts=================================================>
66	//Dhevan's : PhiStarDifferenceMinimum=0.06; EtaDifferenceMinimum=0.02;
67	double PhiStarDifferenceMinimum=0.017; //[radian]
68	double EtaDifferenceMinimum=0.015; //[radian]
69	//=======Double track cuts=================================================<
70
71	// Switches for QA analyses
72	// int runmults[10] = {0, 0, 0, 0, 0, 0, 0, 0, 0, 1};
73	// int multbins[11] = {30, 70, 150, 310, 590, 990, 1570, 2370, 2370, 2370, 6500};
74
75	int runmults[cMu] = {1, 0, 0, 0};
76	//int multbins[11] = {0, 50, 100, 200, 300, 400, 500, 600, 700, 800, 900};
77	int multbins[cMu+1] = {0, 100, 300, 500, 900};
78
79	int runch[2] = {1, 1};
80	const char *chrgs[2] = { "Kp", "Km"};
81
82
83	int runktdep = 1;
84	double ktrng[cKt+1] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0, 1.3};
85	// double ktrng[8] = {0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 2.0};
86
87	int run3d = 0;
88	int runshlcms = 0;
89
90	int runtype = 2; // Types 0 - global, 1 - ITS only, 2 - TPC Inner
91	int isrealdata = 1;
92
93	// AliFemtoEventReaderESDChainKine* Reader=new AliFemtoEventReaderESDChainKine();
94	// Reader->SetConstrained(true);
95	// Reader->SetUseTPCOnly(false);
96
97	double shqmax;
98	double shqmaxSH;
99	int nbinssh = 100;
100
101	//ml if (runshlcms) shqmax = 0.25;
102	// else shqmax = 0.9;
103
104
105	if (runshlcms) shqmaxSH = 0.25;
106	shqmax = 0.9;
107
108
109	// AliFemtoEventReaderESDChain* Reader=new AliFemtoEventReaderESDChain();
110	// Reader->SetConstrained(true);
111	// Reader->SetUseTPCOnly(false);
112	// Reader->SetReadTPCInner(false);
113	// Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);
114
115	// if (runtype == 0)
116	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kGlobal);
117	// else if (runtype == 1)
118	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kITSOnly);
119	// else if (runtype == 2)
120	// Reader->SetReadTrackType(AliFemtoEventReaderESDChain::kTPCOnly);
121	// if (isrealdata)
122	// Reader->SetUsePhysicsSelection(kTRUE);
123	// else
124	// Reader->SetUsePhysicsSelection(kFALSE);
125
126	// Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kV0Centrality);
127
128	AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
129	Reader->SetFilterBit(7);
130	Reader->SetCentralityPreSelection(0, 100);
131	/*On 02/26/2014 06:09 PM, Maciej Szymański wrote:
132	> Dear Konstantin,
133	>
134	> as I remember correctly, you were using the DCA cut on the global tracks
135	> corrsponding to chosen TPC-only tracks (SetDCAglobalTrack(kTRUE), right?
136	> The update I made regards cutting on the DCA on the TPC-only track, so
137	> you need to change in your config to: SetDCAglobalTrack(kFALSE).*/
138	//Reader->SetDCAglobalTrack(kTRUE);//option the DCA information from global tracks (ITS+TPC)
139	Reader->SetDCAglobalTrack(kFALSE);
140
141	AliFemtoManager* Manager=new AliFemtoManager();
142	Manager->SetEventReader(Reader);
143
144	AliFemtoVertexMultAnalysis *anetaphitpc[20];
145	AliFemtoBasicEventCut *mecetaphitpc[20];
146	AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[20];
147	AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[20];
148	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[20];
149	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[20];
150	AliFemtoKKTrackCut *dtc1etaphitpc[20];
151	AliFemtoKKTrackCut *dtc2etaphitpc[20];
152	// AliFemtoESDTrackCut *dtc1etaphitpc[20];
153	// AliFemtoESDTrackCut *dtc2etaphitpc[20];
154	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[20];
155	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[20];
156	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[20];
157	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[20];
158	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[20];
159	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[20];
160	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[20];
161	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[20];
162	//AliFemtoPairCutAntiGamma *sqpcetaphitpc[20];
163	// AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpc[20];
164	//AliFemtoPairCutRadialDistance sqpcetaphitpc[20];//AliFemto dphi cut
165	AliFemtoPairCutRadialDistanceKK sqpcetaphitpc[20];//Dhevan's dphi cut
166	AliFemtoCorrFctnDirectYlm *cylmetaphitpc[20];
167	AliFemtoCorrFctnDEtaDPhi cdedpetaphi[2010];//20->20*10 due to kT
168	AliFemtoChi2CorrFctn *cchiqinvetaphitpc[20];
169	AliFemtoKTPairCut ktpcuts[208];
170	AliFemtoCorrFctnDirectYlm cylmkttpc[208];
171	AliFemtoQinvCorrFctn cqinvkttpc[208];
172	AliFemtoCorrFctn3DLCMSSym cq3dlcmskttpc[208];
173	AliFemtoCorrFctnTPCNcls *cqinvnclstpc[20];
174	AliFemtoShareQualityCorrFctn cqinvsqtpc[2010];
175	AliFemtoChi2CorrFctn *cqinvchi2tpc[20];
176	AliFemtoTPCInnerCorrFctn cqinvinnertpc[2010];
177
178	// * Third QA task - HBT analysis with all pair cuts off, TPC only *
179	// * Begin Kaon-Kaon (positive) analysis *
180	int aniter = 0;
181
182	for (int imult=0; imult<cMu/4/; imult++) {
183	if (runmults[imult]) {
184	for (int ichg=0; ichg<2; ichg++) {
185	if (runch[ichg]) {
186	aniter = ichg*5+imult;
187
188	anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(4, -8.0, 8.0, 5, multbins[imult], multbins[imult+1]);
189	anetaphitpc[aniter]->SetNumEventsToMix(3);
190	anetaphitpc[aniter]->SetMinSizePartCollection(1);
191
192	mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
193	mecetaphitpc[aniter]->SetEventMult(0,100000);
194	mecetaphitpc[aniter]->SetVertZPos(-8.0,8.0);
195	/* //was in aliroot 5.03.76
196	if (isrealdata)
197	mecetaphitpc[aniter]->SetAcceptOnlyPhysics(kTRUE);
198	*/
199	// mecetaphitpc->SetAcceptBadVertex(kTRUE);
200
201	cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
202	cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
203	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
204
205	cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
206	cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
207	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
208
209	dtc1etaphitpc[aniter] = new AliFemtoKKTrackCut();
210	// dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
211	// dtc1etaphitpc[aniter]->SetPidProbPion(0.0,1.001);
212	// dtc1etaphitpc[aniter]->SetPidProbMuon(0.0,1.0);
213	// dtc1etaphitpc[aniter]->SetPidProbKaon(0.0,1.0);
214	// dtc1etaphitpc[aniter]->SetPidProbProton(0.0,1.0);
215	if (ichg == 0)
216	dtc1etaphitpc[aniter]->SetCharge(1.0);
217	else if (ichg == 1)
218	dtc1etaphitpc[aniter]->SetCharge(-1.0);
219
220	dtc1etaphitpc[aniter]->SetPt(0.14,1.5);
221	// dtc1etaphitpc[aniter]->SetEta(-1.2,1.2);
222	dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
223	// // dtc1etaphitpc[aniter]->SetEta(-0.5,0.5);
224	/// dtc1etaphitpc[aniter]->SetMass(PionMass);
225	dtc1etaphitpc[aniter]->SetMass(KaonMass);
226
227
228	//// dtc1etaphitpc[aniter]->SetminTPCncls(80);
229
230	//--- Choose Kaon as Most Probable (switch on all cuts: TPC, TOF)---
231	dtc1etaphitpc[aniter]->SetMostProbableKaon();
232	//------------------- November 2013 -----------------------------------<
233	// new cuts to remove electron (do not take into analysis if 400<p<500)
234	//remove 0.0 to include electrons
235	dtc1etaphitpc[aniter]->SetNsigmaTPCle250(1.5);
236	dtc1etaphitpc[aniter]->SetNsigmaTPC250_400(2.0);
237	dtc1etaphitpc[aniter]->SetNsigmaTPC400_450(2.0);//0.0);
238	dtc1etaphitpc[aniter]->SetNsigmaTPC450_500(2.0);//0.0);
239	dtc1etaphitpc[aniter]->SetNsigmaTPCge500(3.0);
240	// new cuts are stronger, better separation of pion in TOF
241	// when momentum is greater then 800 MeV/c
242	dtc1etaphitpc[aniter]->SetNsigmaTOF500_800(2.0);
243	dtc1etaphitpc[aniter]->SetNsigmaTOF800_1000(1.5);
244	dtc1etaphitpc[aniter]->SetNsigmaTOFge1000(1.0);
245	//------------------- November 2013 ----------------------------------->
246
247	//// dtc1etaphitpc[aniter]->SetMostProbablePion();
248	// // Track quality cuts
249	if (runtype == 0) {
250	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
251	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
252	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
253	dtc1etaphitpc[aniter]->SetminTPCncls(80);
254	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
255	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
256	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
257	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
258	dtc1etaphitpc[aniter]->SetMaxImpactXY(DCAxy);
259	//Poland: dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
260	dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz);
261	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
262	}
263	else if (runtype == 1) {
264	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
265	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
266	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit\|AliESDtrack::kITSpureSA);
267	// dtc1etaphitpc[aniter]->SetminTPCncls(70);
268	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
269	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
270	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
271	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
272	// dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(6.0);
273	dtc1etaphitpc[aniter]->SetMaxImpactXY(DCAxy);
274	dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz);
275	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
276	}
277	else if (runtype == 2) {
278	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
279	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
280	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit);
281	// dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kITSrefit);
282	dtc1etaphitpc[aniter]->SetminTPCncls(80); //was "0"
283	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
284	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
285	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
286	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
287	dtc1etaphitpc[aniter]->SetMaxImpactXY(DCAxy);
288	//dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
289	dtc1etaphitpc[aniter]->SetMaxImpactZ(DCAz); //3.0
290	// dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
291	}
292
293
294
295	cutPass1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass1%stpcM%i", chrgs[ichg], imult), 0.493677);
296	cutFail1YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail1%stpcM%i", chrgs[ichg], imult), 0.493677);
297	dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1YPtetaphitpc[aniter], cutFail1YPtetaphitpc[aniter]);
298
299	cutPass1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass1%stpcM%i", chrgs[ichg], imult),1);
300	cutFail1PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail1%stpcM%i", chrgs[ichg], imult),1);
301	dtc1etaphitpc[aniter]->AddCutMonitor(cutPass1PIDetaphitpc[aniter], cutFail1PIDetaphitpc[aniter]);
302
303
304	// sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
305	// sqpcetaphitpc[aniter] = new AliFemtoShareQualityTPCEntranceSepPairCut();
306
307	if (ichg < 2) {
308	//sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();//AliFemto dphi* cut
309	sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistanceKK(); //Dhevan's dphi* cut
310	if (runtype == 0) {
311	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
312	sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
313	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
314	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
315	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
316	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
317	//ml sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
318	//ml sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
319
320
321	//--------- km: eta-phi* Dhevan's custs ----------->>>>
322	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
323	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
324	//sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
325	//--------- km: eta-phi* Dhevan's custs -----------<<<
326
327	//////////////sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
328
329
330
331	}
332	else if (runtype == 1) {
333	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
334	sqpcetaphitpc[aniter]->SetShareFractionMax(1.05);
335	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
336	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
337	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
338	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(5.0);
339	// sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.03);
340	// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
341
342	// void SetPhiStarDifferenceMinimum(double dtpc);
343	// void SetEtaDifferenceMinimum(double etpc);
344	// void SetMinimumRadius(double minrad);
345	// void SetMagneticFieldSign(int magsign);
346
347	//--------- km: eta-phi* Dhevan's custs ----------->>>>
348	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
349	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
350	//sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
351	//--------- km: eta-phi* Dhevan's custs -----------<<<
352
353	/////////sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
354	//sqpcetaphitpc[aniter]->SetMagneticFieldSign(1.0);
355
356
357	}
358	else if (runtype == 2) {
359	sqpcetaphitpc[aniter]->SetDataType(AliFemtoPairCut::kAOD);
360	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
361	sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
362	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
363	// sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
364	// sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
365	//ml sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(0.0);
366	//ml sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(1.2, 0.045);
367	//ml sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.016);
368
369	//--------- km: eta-phi* Dhevan's custs ----------->>>>
370	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(PhiStarDifferenceMinimum);
371	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(EtaDifferenceMinimum);
372	//sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);//not need for AliFemtoPairCutRadialDistanceKK()
373	//--------- km: eta-phi* Dhevan's custs -----------<<<
374
375	////////sqpcetaphitpc[aniter]->SetMagneticFieldSign(1);
376
377
378	}
379	}
380
381
382	anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
383	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
384	anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
385	anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
386
387	// cylmetaphitpc[aniter] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%i", chrgs[ichg], imult),3,nbinssh,0.0,shqmaxSH,runshlcms);
388	/// anetaphitpc[aniter]->AddCorrFctn(cylmetaphitpc[aniter]);
389
390	// cqinvnclstpc[aniter] = new AliFemtoCorrFctnTPCNcls(Form("cqinvncls%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
391	// anetaphitpc[aniter]->AddCorrFctn(cqinvnclstpc[aniter]);
392
393	// cqinvchi2tpc[aniter] = new AliFemtoChi2CorrFctn(Form("cqinvchi2%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
394	// anetaphitpc[aniter]->AddCorrFctn(cqinvchi2tpc[aniter]);
395
396	if (runktdep) {
397	int ktm;
398	for (int ikt=0; ikt<cKt/8/; ikt++) {
399	ktm = anitercKt/8*/ + ikt;
400	ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
401
402	//////// cylmkttpc[ktm] = new AliFemtoCorrFctnDirectYlm(Form("cylm%stpcM%ikT%i", chrgs[ichg], imult, ikt),3,
403	/////// nbinssh, 0.0,
404	////// (imult>6)?shqmaxSH*2.5:shqmaxSH,
405	///// runshlcms);
406	////// cylmkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
407	///// anetaphitpc[aniter]->AddCorrFctn(cylmkttpc[ktm]);
408
409	cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
410	// cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,0.5);
411	cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
412	anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
413
414	cqinvsqtpc[ktm] = new AliFemtoShareQualityCorrFctn(Form("cqinvsq%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
415	cqinvsqtpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
416	anetaphitpc[aniter]->AddCorrFctn(cqinvsqtpc[ktm]);
417
418	cqinvinnertpc[ktm] = new AliFemtoTPCInnerCorrFctn(Form("cqinvinner%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,shqmax);
419	cqinvinnertpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
420	cqinvinnertpc[ktm]->SetRadius(1.2);
421	anetaphitpc[aniter]->AddCorrFctn(cqinvinnertpc[ktm]);
422
423	//---- Correlation Function vs Delta_Eta and Delta_Phi (not Phi*)---->>>
424	cdedpetaphi[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ikT%i", chrgs[ichg], imult, ikt),100,100);
425	anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[ktm]);
426	//---- Correlation Function vs Delta_Eta and Delta_Phi (not Phi*)----<<<
427
428	if (run3d) {
429	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);
430	// cq3dlcmskttpc[ktm] = new AliFemtoCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),50,0.5);
431	cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
432	anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
433	}
434	}
435	}
436
437	// cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),24, 24);
438	// anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
439
440	Manager->AddAnalysis(anetaphitpc[aniter]);
441	}
442	}
443	}
444	}
445	// *** End Kaon-Kaon (positive) analysis
446
447	return Manager;
448	}
449