[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / KchKch_MC_ESD / 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 "AliFemtoModelCorrFctn.h"
#include "AliFemtoCutMonitorParticlePtPDG.h"
#include "AliFemtoKTPairCut.h"
#include "AliFemtoCutMonitorCollections.h"
#endif
 
//_  6 fm  pure kaons _(0.2-0.4)______________________________________________________________________
AliFemtoManager* ConfigFemtoAnalysis() {

  double PionMass = 0.13956995;
  double KaonMass = 0.493677;
  double ProtonMass = 0.938272013;

  //multiplicity bins
  int runmults[10] = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
  int multbins[11] = {0.01, 200000, 400, 600, 900, 950, 500, 600, 700, 800, 900};

  int runch[3] = {1, 1, 0};
  const char *chrgs[3] = { "kk", "akak", "kak" };

  int runktdep = 1;
 // double ktrng[8] = {0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 1.0};
  double ktrng[4] = {0.2, 0.4, 0.6, 1.0};

  int run3d = 0; // Do 3D cartesian analysis?
  int runshlcms = 1;

  int nbinssh = 100;
  double shqmax = 0.5;


  AliFemtoEventReaderESDChainKine *Reader = new AliFemtoEventReaderESDChainKine();
  Reader->SetKaonAnalysis(kTRUE);
  //Reader->SetUseMultiplicity(AliFemtoEventReaderESDChainKine::kCentrality);
  //Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);

  // AliFemtoEventReaderKinematicsChain *Reader = new AliFemtoEventReaderKinematicsChain();
  // Reader->SetUseMultiplicity(AliFemtoEventReaderKinematicsChain::kGlobalCount);
  // // // Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);

  AliFemtoModelGausRinvFreezeOutGenerator *tFreeze = new AliFemtoModelGausRinvFreezeOutGenerator();
  tFreeze->SetSizeInv(5.0);

 // AliFemtoModelWeightGeneratorLednicky *tWeight = new AliFemtoModelWeightGeneratorLednicky();
  // tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
   
    AliFemtoModelWeightGeneratorBasic *tWeight = new AliFemtoModelWeightGeneratorBasic();
        tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
   
 //////// tWeight->SetPairType(5);
  //tWeight->SetPairType(1);

  AliFemtoModelManager *tModelManager = new AliFemtoModelManager();
  tModelManager->AcceptFreezeOutGenerator(tFreeze);
  tModelManager->AcceptWeightGenerator(tWeight);
  tModelManager->CreateCopyHiddenInfo(kFALSE);
  //tModelManager->CreateCopyHiddenInfo(kTRUE);

  AliFemtoManager* Manager=new AliFemtoManager();
  Manager->SetEventReader(Reader);

  AliFemtoVertexMultAnalysis    *anetaphitpc[30];
  AliFemtoBasicEventCut         *mecetaphitpc[30];
  AliFemtoCutMonitorEventMult   *cutPassEvMetaphitpc[30];
  AliFemtoCutMonitorEventMult   *cutFailEvMetaphitpc[30];
  AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[30];
  AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[30];
  AliFemtoCutMonitorCollections   *cutPassColletaphitpc[30];
  AliFemtoCutMonitorCollections   *cutFailColletaphitpc[30];
  AliFemtoESDTrackCut           *dtc1etaphitpc[30];
  AliFemtoESDTrackCut           *dtc2etaphitpc[30];
  AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[30];
  AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[30];
  AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[30];
  AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[30];
  AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[30];
  AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[30];
  AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[30];
  AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[30];
  //  AliFemtoPairCutAntiGamma      *sqpcetaphitpc[30];
  AliFemtoShareQualityPairCut      *sqpcetaphitpc[30];
  //AliFemtoPairCutRadialDistance      *sqpcetaphitpc[30];
  AliFemtoModelCorrFctnDirectYlm     *cylmetaphitpc[30];
  AliFemtoCorrFctnDEtaDPhi      *cdedpetaphi[30];
  AliFemtoChi2CorrFctn          *cchiqinvetaphitpc[30];
  AliFemtoKTPairCut             *ktpcuts[30*7];
  AliFemtoModelCorrFctnDirectYlm     *cylmkttpc[30*7];
  //AliFemtoModelQinvCorrFctn          *cqinvkttpc[30*7];
  AliFemtoModelCorrFctnSource          *cqinvkttpc[30*7];
//  AliFemtoModelCorrFctn          *cqinvkttpc[30*7];
  //AliFemtoModelCorrFctn3DLCMSSym     *cq3dlcmskttpc[30*7];
  AliFemtoCorrFctnTPCNcls       *cqinvnclstpc[30];
  AliFemtoShareQualityCorrFctn  *cqinvsqtpc[30*10];
  AliFemtoChi2CorrFctn          *cqinvchi2tpc[30];
  AliFemtoTPCInnerCorrFctn      *cqinvinnertpc[30*10];


  // *** Begin pion-pion analysis ***
  int aniter = 0;

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

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

          mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
          mecetaphitpc[aniter]->SetEventMult(0,200000);
          mecetaphitpc[aniter]->SetVertZPos(-8,8);


          dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
          dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();


          if (ichg == 0) {
            dtc1etaphitpc[aniter]->SetCharge(1.0);
            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
          }
          else if (ichg == 1) {
            dtc1etaphitpc[aniter]->SetCharge(-1.0);
            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
          }
          else if (ichg == 2) {
            dtc1etaphitpc[aniter]->SetCharge(-1.0);
            dtc2etaphitpc[aniter]->SetCharge(1.0);
            dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
            dtc2etaphitpc[aniter]->SetPt(0.2,1.5);
          }

          dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
          dtc1etaphitpc[aniter]->SetMass(KaonMass);
//test          dtc1etaphitpc[aniter]->SetMostProbableKaon();
//test          dtc1etaphitpc[aniter]->SetNsigma(3.0);
          //dtc1etaphitpc[aniter]->SetNsigma(2.0);
//test          dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
          //dtc1etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);

          if (ichg == 2) {
            dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
            dtc2etaphitpc[aniter]->SetMass(KaonMass);
//test            dtc2etaphitpc[aniter]->SetMostProbableKaon();
//test            dtc2etaphitpc[aniter]->SetNsigma(3.0);
            //dtc2etaphitpc[aniter]->SetNsigma(2.0);
//test            dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
            //dtc2etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);

          }


           dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
           dtc1etaphitpc[aniter]->SetminTPCncls(80);
           dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
           dtc1etaphitpc[aniter]->SetLabel(kFALSE);
           dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
           dtc1etaphitpc[aniter]->SetMaxImpactXY(2.4); // 2.4 0.1
           dtc1etaphitpc[aniter]->SetMaxImpactZ(3.2); // 2.0 0.1

           AliFemtoCutMonitorParticleYPt *cutPassYPtkp = new AliFemtoCutMonitorParticleYPt("cutPasskp", 0.493677);
           AliFemtoCutMonitorParticleYPt *cutFailYPtkp = new AliFemtoCutMonitorParticleYPt("cutFailkp", 0.493677);
           dtc1etaphitpc[aniter]->AddCutMonitor(cutPassYPtkp, cutFailYPtkp);

           AliFemtoCutMonitorParticlePtPDG *cutPassPidkp = new AliFemtoCutMonitorParticlePtPDG("cutPasskp", 0.493677);
           AliFemtoCutMonitorParticlePtPDG *cutFailPidkp = new AliFemtoCutMonitorParticlePtPDG("cutFailkp", 0.493677);
           dtc1etaphitpc[aniter]->AddCutMonitor(cutPassPidkp, cutFailPidkp);

           AliFemtoCutMonitorParticleMomRes *cutPassMRkp = new AliFemtoCutMonitorParticleMomRes("cutPasskp");
           AliFemtoCutMonitorParticleMomRes *cutFailMRkp = new AliFemtoCutMonitorParticleMomRes("cutFailkp");
           dtc1etaphitpc[aniter]->AddCutMonitor(cutPassMRkp, cutFailMRkp);


          sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
          //sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
          sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
          sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
          sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);


          anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
          if (ichg == 2) {
            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
            anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
          }
          else {
            anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
            anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
          }

          anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);


          //Qinv (without kT bins)
          cqinvkttpc[aniter] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
          //cqinvkttpc[aniter] = new AliFemtoModelCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
        //  ktpcuts[0] = new AliFemtoKTPairCut(0.2,0.4);
         // cqinvkttpc[aniter]->SetPairSelectionCut(ktpcuts[0]);
          cqinvkttpc[aniter]->ConnectToManager(tModelManager);
          anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[aniter]);


          if (runktdep) {
            int ktm;
            for (int ikt=0; ikt<3; ikt++) {
              ktm = aniter*3 + ikt;
              ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);

              cqinvkttpc[ktm] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0, shqmax);
              cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
              cqinvkttpc[ktm]->ConnectToManager(tModelManager);
              anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[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 AliFemtoModelCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,0.5);
              //cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
              //anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
              //}
            }
          }

          //cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),39, 39);
          //anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);

          Manager->AddAnalysis(anetaphitpc[aniter]);
        }
      }
    }
  }
  // *** End pion-pion analysis

  return Manager;
}
Commit	Line	Data
1a3b11ce	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 "AliFemtoModelCorrFctn.h"
	49	#include "AliFemtoCutMonitorParticlePtPDG.h"
	50	#include "AliFemtoKTPairCut.h"
	51	#include "AliFemtoCutMonitorCollections.h"
	52	#endif
	53
	54	//_ 6 fm pure kaons _(0.2-0.4)______________________________________________________________________
	55	AliFemtoManager* ConfigFemtoAnalysis() {
	56
	57	double PionMass = 0.13956995;
	58	double KaonMass = 0.493677;
	59	double ProtonMass = 0.938272013;
	60
	61	//multiplicity bins
	62	int runmults[10] = {1, 0, 0, 0, 0, 0, 0, 0, 0, 0};
	63	int multbins[11] = {0.01, 200000, 400, 600, 900, 950, 500, 600, 700, 800, 900};
	64
65	int runch[3] = {1, 1, 0};
66	const char *chrgs[3] = { "kk", "akak", "kak" };
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};
1e1dd217	70	double ktrng[4] = {0.2, 0.4, 0.6, 1.0};
1a3b11ce	71
	72	int run3d = 0; // Do 3D cartesian analysis?
	73	int runshlcms = 1;
	74
	75	int nbinssh = 100;
	76	double shqmax = 0.5;
	77
	78
	79	AliFemtoEventReaderESDChainKine *Reader = new AliFemtoEventReaderESDChainKine();
1e1dd217	80	Reader->SetKaonAnalysis(kTRUE);
1a3b11ce	81	//Reader->SetUseMultiplicity(AliFemtoEventReaderESDChainKine::kCentrality);
	82	//Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);
	83
	84	// AliFemtoEventReaderKinematicsChain *Reader = new AliFemtoEventReaderKinematicsChain();
	85	// Reader->SetUseMultiplicity(AliFemtoEventReaderKinematicsChain::kGlobalCount);
	86	// // // Reader->SetReadTrackType(AliFemtoEventReaderESDChainKine::kGlobal);
	87
	88	AliFemtoModelGausRinvFreezeOutGenerator *tFreeze = new AliFemtoModelGausRinvFreezeOutGenerator();
	89	tFreeze->SetSizeInv(5.0);
	90
	91	// AliFemtoModelWeightGeneratorLednicky *tWeight = new AliFemtoModelWeightGeneratorLednicky();
1e1dd217	92	// tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
	93
	94	AliFemtoModelWeightGeneratorBasic *tWeight = new AliFemtoModelWeightGeneratorBasic();
	95	tWeight->SetPairType(AliFemtoModelWeightGenerator::KaonPlusKaonPlus());
	96
1a3b11ce	97	//////// tWeight->SetPairType(5);
	98	//tWeight->SetPairType(1);
	99
	100	AliFemtoModelManager *tModelManager = new AliFemtoModelManager();
	101	tModelManager->AcceptFreezeOutGenerator(tFreeze);
	102	tModelManager->AcceptWeightGenerator(tWeight);
	103	tModelManager->CreateCopyHiddenInfo(kFALSE);
	104	//tModelManager->CreateCopyHiddenInfo(kTRUE);
	105
	106	AliFemtoManager* Manager=new AliFemtoManager();
	107	Manager->SetEventReader(Reader);
	108
	109	AliFemtoVertexMultAnalysis *anetaphitpc[30];
	110	AliFemtoBasicEventCut *mecetaphitpc[30];
	111	AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[30];
	112	AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[30];
	113	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[30];
	114	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[30];
	115	AliFemtoCutMonitorCollections *cutPassColletaphitpc[30];
	116	AliFemtoCutMonitorCollections *cutFailColletaphitpc[30];
	117	AliFemtoESDTrackCut *dtc1etaphitpc[30];
	118	AliFemtoESDTrackCut *dtc2etaphitpc[30];
	119	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[30];
	120	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[30];
	121	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[30];
	122	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[30];
	123	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[30];
	124	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[30];
	125	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[30];
	126	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[30];
	127	// AliFemtoPairCutAntiGamma *sqpcetaphitpc[30];
	128	AliFemtoShareQualityPairCut *sqpcetaphitpc[30];
	129	//AliFemtoPairCutRadialDistance *sqpcetaphitpc[30];
	130	AliFemtoModelCorrFctnDirectYlm *cylmetaphitpc[30];
	131	AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[30];
	132	AliFemtoChi2CorrFctn *cchiqinvetaphitpc[30];
	133	AliFemtoKTPairCut ktpcuts[307];
	134	AliFemtoModelCorrFctnDirectYlm cylmkttpc[307];
	135	//AliFemtoModelQinvCorrFctn cqinvkttpc[307];
	136	AliFemtoModelCorrFctnSource cqinvkttpc[307];
	137	// AliFemtoModelCorrFctn cqinvkttpc[307];
	138	//AliFemtoModelCorrFctn3DLCMSSym cq3dlcmskttpc[307];
	139	AliFemtoCorrFctnTPCNcls *cqinvnclstpc[30];
	140	AliFemtoShareQualityCorrFctn cqinvsqtpc[3010];
	141	AliFemtoChi2CorrFctn *cqinvchi2tpc[30];
	142	AliFemtoTPCInnerCorrFctn cqinvinnertpc[3010];
	143
	144
	145	// * Begin pion-pion analysis *
	146	int aniter = 0;
	147
	148	for (int imult=0; imult<10; imult++) {
	149	if (runmults[imult]) {
	150	for (int ichg=0; ichg<3; ichg++) {
	151	if (runch[ichg]) {
	152	aniter = ichg*10+imult;
	153
	154	anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(10, -8.0, 8.0, 4, multbins[imult], multbins[imult+1]);
	155	anetaphitpc[aniter]->SetNumEventsToMix(5);
	156	anetaphitpc[aniter]->SetMinSizePartCollection(1);
	157
	158	mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
	159	mecetaphitpc[aniter]->SetEventMult(0,200000);
	160	mecetaphitpc[aniter]->SetVertZPos(-8,8);
161
162
163	dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
164	dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
1e1dd217	165
	166
	167
1a3b11ce	168	if (ichg == 0) {
	169	dtc1etaphitpc[aniter]->SetCharge(1.0);
	170	dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
	171	}
	172	else if (ichg == 1) {
	173	dtc1etaphitpc[aniter]->SetCharge(-1.0);
	174	dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
	175	}
	176	else if (ichg == 2) {
	177	dtc1etaphitpc[aniter]->SetCharge(-1.0);
	178	dtc2etaphitpc[aniter]->SetCharge(1.0);
	179	dtc1etaphitpc[aniter]->SetPt(0.2,1.5);
	180	dtc2etaphitpc[aniter]->SetPt(0.2,1.5);
	181	}
	182
	183	dtc1etaphitpc[aniter]->SetEta(-0.8,0.8);
	184	dtc1etaphitpc[aniter]->SetMass(KaonMass);
	185	//test dtc1etaphitpc[aniter]->SetMostProbableKaon();
	186	//test dtc1etaphitpc[aniter]->SetNsigma(3.0);
	187	//dtc1etaphitpc[aniter]->SetNsigma(2.0);
	188	//test dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	189	//dtc1etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);
	190
	191	if (ichg == 2) {
	192	dtc2etaphitpc[aniter]->SetEta(-0.8,0.8);
	193	dtc2etaphitpc[aniter]->SetMass(KaonMass);
	194	//test dtc2etaphitpc[aniter]->SetMostProbableKaon();
	195	//test dtc2etaphitpc[aniter]->SetNsigma(3.0);
	196	//dtc2etaphitpc[aniter]->SetNsigma(2.0);
	197	//test dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	198	//dtc2etaphitpc[aniter]->SetNsigmaTPConly(kTRUE);
	199
	200	}
	201
	202
	203	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCin);
	204	dtc1etaphitpc[aniter]->SetminTPCncls(80);
	205	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
	206	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
	207	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
	208	dtc1etaphitpc[aniter]->SetMaxImpactXY(2.4); // 2.4 0.1
	209	dtc1etaphitpc[aniter]->SetMaxImpactZ(3.2); // 2.0 0.1
1e1dd217	210
	211	AliFemtoCutMonitorParticleYPt *cutPassYPtkp = new AliFemtoCutMonitorParticleYPt("cutPasskp", 0.493677);
	212	AliFemtoCutMonitorParticleYPt *cutFailYPtkp = new AliFemtoCutMonitorParticleYPt("cutFailkp", 0.493677);
	213	dtc1etaphitpc[aniter]->AddCutMonitor(cutPassYPtkp, cutFailYPtkp);
	214
	215	AliFemtoCutMonitorParticlePtPDG *cutPassPidkp = new AliFemtoCutMonitorParticlePtPDG("cutPasskp", 0.493677);
	216	AliFemtoCutMonitorParticlePtPDG *cutFailPidkp = new AliFemtoCutMonitorParticlePtPDG("cutFailkp", 0.493677);
	217	dtc1etaphitpc[aniter]->AddCutMonitor(cutPassPidkp, cutFailPidkp);
	218
	219	AliFemtoCutMonitorParticleMomRes *cutPassMRkp = new AliFemtoCutMonitorParticleMomRes("cutPasskp");
	220	AliFemtoCutMonitorParticleMomRes *cutFailMRkp = new AliFemtoCutMonitorParticleMomRes("cutFailkp");
	221	dtc1etaphitpc[aniter]->AddCutMonitor(cutPassMRkp, cutFailMRkp);
	222
	223
1a3b11ce	224	sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
	225	//sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
	226	sqpcetaphitpc[aniter]->SetShareQualityMax(1.0);
	227	sqpcetaphitpc[aniter]->SetShareFractionMax(0.05);
	228	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
	229
	230
	231	anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
	232	if (ichg == 2) {
	233	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
	234	anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
	235	}
	236	else {
	237	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
	238	anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
	239	}
	240
	241	anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
	242
	243
	244	//Qinv (without kT bins)
	245	cqinvkttpc[aniter] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	246	//cqinvkttpc[aniter] = new AliFemtoModelCorrFctn(Form("cqinv%stpcM%i", chrgs[ichg], imult),nbinssh,0.0,shqmax);
	247	// ktpcuts[0] = new AliFemtoKTPairCut(0.2,0.4);
	248	// cqinvkttpc[aniter]->SetPairSelectionCut(ktpcuts[0]);
	249	cqinvkttpc[aniter]->ConnectToManager(tModelManager);
	250	anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[aniter]);
	251
	252
	253	if (runktdep) {
	254	int ktm;
	255	for (int ikt=0; ikt<3; ikt++) {
	256	ktm = aniter*3 + ikt;
	257	ktpcuts[ktm] = new AliFemtoKTPairCut(ktrng[ikt], ktrng[ikt+1]);
	258
	259	cqinvkttpc[ktm] = new AliFemtoModelCorrFctnSource(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0, shqmax);
	260	cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	261	cqinvkttpc[ktm]->ConnectToManager(tModelManager);
	262	anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
	263
	264	//if (run3d) {
	265	// 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);
	266	//cq3dlcmskttpc[ktm] = new AliFemtoModelCorrFctn3DLCMSSym(Form("cq3d%stpcM%ikT%i", chrgs[ichg], imult, ikt),60,0.5);
	267	//cq3dlcmskttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
	268	//anetaphitpc[aniter]->AddCorrFctn(cq3dlcmskttpc[ktm]);
	269	//}
	270	}
	271	}
	272
	273	//cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),39, 39);
	274	//anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
	275
	276	Manager->AddAnalysis(anetaphitpc[aniter]);
	277	}
	278	}
	279	}
	280	}
	281	// *** End pion-pion analysis
	282
	283	return Manager;
	284	}