[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / DEtaDPhi / 2014 / Train4 / NoCorrections / 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 "AliFemtoCorrFctnNonIdDR.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 "AliFemtoCutMonitorParticlePtPDG.h"
#include "AliFemtoKTPairCut.h"
#include "AliFemtoPairCutPt.h"
#endif

//________________________________________________________________________
AliFemtoManager* ConfigFemtoAnalysis(const char* params="") {

	double PionMass = 0.13956995;
	double KaonMass = 0.493677;
	double ProtonMass = 0.938272013;
	
	const int numOfMultBins = 5;	
	const int numOfChTypes = 13;
	const int numOfkTbins = 6;

	char *parameter[15];
	if(strlen(params)!=0)
	  {
	    parameter[0] = strtok(params, ","); // Splits spaces between words in params
	    cout<<"Parameter [0] (filterbit):"<<parameter[0]<<endl; // Writes first parameter
	    parameter[1] = strtok(NULL, ",");
	    cout<<"Parameter [1] (ktdep):"<<parameter[1]<<" "<<endl;
	    parameter[2] = strtok(NULL, ",");
	    cout<<"Parameter [2] (multdep):"<<parameter[2]<<" "<<endl;	  
	    parameter[3] = strtok(NULL, ",");
	    cout<<"Parameter [3]: (MinPlpContribSPD)"<<parameter[3]<<" "<<endl;
	    parameter[4] = strtok(NULL, ",");
	    cout<<"Parameter [4]: (multbino)"<<parameter[4]<<" "<<endl;
	    parameter[5] = strtok(NULL, ",");
	    cout<<"Parameter [5]: (zvertbino)"<<parameter[5]<<" "<<endl;
	    parameter[6] = strtok(NULL, ",");
	    cout<<"Parameter [6]: (ifGlobalTracks=true/false)"<<parameter[6]<<" "<<endl;
	    parameter[7] = strtok(NULL, ",");
	    cout<<"Parameter [7]: (shareQuality)"<<parameter[7]<<" "<<endl;
	    parameter[8] = strtok(NULL, ",");
	    cout<<"Parameter [8]: (shareFraction)"<<parameter[8]<<" "<<endl;
	    parameter[9] = strtok(NULL, ",");
	    cout<<"Parameter [9]: (ifElectronRejection)"<<parameter[9]<<" "<<endl;
	    parameter[10] = strtok(NULL, ",");
	    cout<<"Parameter [10]: (nSigma)"<<parameter[10]<<" "<<endl;
	    parameter[11] = strtok(NULL, ",");
	    cout<<"Parameter [12]: (etaMin)"<<parameter[11]<<" "<<endl;
	    parameter[12] = strtok(NULL, ",");
	    cout<<"Parameter [12]: (etaMax)"<<parameter[12]<<" "<<endl;
	    parameter[13] = strtok(NULL, ",");
	    cout<<"Parameter [13]: (ispileup)"<<parameter[13]<<" "<<endl;
	  }
	int filterbit = atoi(parameter[0]); //96 / 768 / 128 
	int runktdep = atoi(parameter[1]); //0
	int runmultdep = atoi(parameter[2]); //0
	int minPlpContribSPD = atoi(parameter[3]); //3
	int multbino = atoi(parameter[4]); //30
	int zvertbino = atoi(parameter[5]); //10
	Bool_t ifGlobalTracks=kFALSE; if(atoi(parameter[6]))ifGlobalTracks=kTRUE;//kTRUE 
	double shareQuality = atof(parameter[7]); //0.00
	double shareFraction = atof(parameter[8]); //0.05
	bool ifElectronRejection = atoi(parameter[9]); //true
	double nSigmaVal = atof(parameter[10]); //3.0
	double nEtaMin = atof(parameter[11]); //-0.8
	double nEtaMax = atof(parameter[12]);  //0.8
	bool ifIsPileUp = atoi(parameter[13]); //true

	int runmults[numOfMultBins] = {0, 0, 0, 0, 1};
	if(runmultdep)	  {runmults[0]=1; runmults[1]=1; runmults[2]=1;	  }
	int multbins[numOfMultBins+1] = {2, 20, 50,150,2,150};
	
	int runch[numOfChTypes] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
	const char *chrgs[numOfChTypes] = { "PP", "aPaP", "PaP", "KpKp", "KmKm", "KpKm", "PIpPIp", "PImPIm", "PIpPIm", "all", "plus", "minus", "mixed" };
	
	int runktdep = 1;
	double ktrng[numOfkTbins+1] = {0.0, 0, 0, 0, 0, 0};
	double ktrngAll[numOfkTbins+1] = {0.0, 7.0, 1.2, 2.0, 3.0, 4.0, 100.0};
	double ktrngPion[numOfkTbins+1] = {0.0, 0.7, 0.9, 1.1, 1.4, 2.5, 100.0};
	double ktrngKaon[numOfkTbins+1] = {0.0, 1.5, 2.5, 3.5, 100.0, 0, 0};
	double ktrngProton[numOfkTbins+1] = {0.0, 2.75, 100, 0, 0};

	int runqinv = 1;
	int runshlcms = 1;// 0:PRF(PAP), 1:LCMS(PP,APAP)

	int runtype = 0; // Types 0 - global, 1 - ITS only, 2 - TPC Inner	//global tracks ->mfit ITS+TPC
	int owncuts = 0; 
	int owndca = 0;

	int gammacut = 0;	// cut na ee z gamma 
	
	double shqmax = 0.5; 
	int nbinssh = 100;

	//AliFemtoEventReaderESDChain *Reader = new AliFemtoEventReaderESDChain();
	//Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kGlobalCount);


	AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
	Reader->SetFilterMask(filterbit);
	Reader->SetDCAglobalTrack(ifGlobalTracks); //false for FB7, true for the rest //we do not use DCA at all
	Reader->SetUseMultiplicity(AliFemtoEventReaderAOD::kReference);
	Reader->SetMinPlpContribSPD(minPlpContribSPD);
	Reader->SetIsPileUpEvent(ifIsPileUp);

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

	AliFemtoVertexMultAnalysis		*anetaphitpc[320];
	AliFemtoBasicEventCut				 *mecetaphitpc[320];
	AliFemtoCutMonitorEventMult	 *cutPassEvMetaphitpc[320];
	AliFemtoCutMonitorEventMult	 *cutFailEvMetaphitpc[320];
	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[320];
	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[320];
	AliFemtoESDTrackCut			 *dtc1etaphitpc[320];
	AliFemtoESDTrackCut			 *dtc2etaphitpc[320];
	AliFemtoESDTrackCut			 *dtc3etaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutPass3YPtetaphitpc[320];
	AliFemtoCutMonitorParticleYPt *cutFail3YPtetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutPass3PIDetaphitpc[320];
	AliFemtoCutMonitorParticlePID *cutFail3PIDetaphitpc[320];
	//	 AliFemtoShareQualityTPCEntranceSepPairCut			*sqpcetaphitpcsame[320];
	//AliFemtoPairCutAntiGamma			*sqpcetaphitpc[320];
	//AliFemtoPairCutRadialDistance			*sqpcetaphitpc[320];
	AliFemtoShareQualityPairCut			*sqpcetaphitpc[320];
	//	AliFemtoChi2CorrFctn					*cchiqinvetaphitpc[320];
	AliFemtoPairCutPt						 *ktpcuts[1020];
	AliFemtoQinvCorrFctn					*cqinvkttpc[320];
	AliFemtoQinvCorrFctn					*cqinvtpc[320];
	AliFemtoCorrFctnDEtaDPhi			*cdedpetaphi[320];
	AliFemtoCorrFctnDEtaDPhi			*cdedpetaphiPt[1320];


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

	for (int imult = 0; imult < numOfMultBins; imult++)
	{
		if (runmults[imult])
		{
			for (int ichg = 0; ichg < numOfChTypes; ichg++)
			{
				if (runch[ichg])
				{

					aniter = ichg * numOfMultBins + imult;
					anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(zvertbino, -10.0, 10.0, multbino, multbins[imult], multbins[imult+1]);
					anetaphitpc[aniter]->SetNumEventsToMix(10);
					anetaphitpc[aniter]->SetMinSizePartCollection(1);
					anetaphitpc[aniter]->SetVerboseMode(kFALSE);//~~~~~~~~~~~~~~~~

					//*** Event cut ***
					mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
					mecetaphitpc[aniter]->SetEventMult(0.001,100000);
					mecetaphitpc[aniter]->SetVertZPos(-10,10);//cm

					//****** event monitors **********	
					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]);


					// ***** single particle track cuts *********
					dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
					dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
					dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
					dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
					dtc1etaphitpc[aniter]->SetNsigma(nSigmaVal);
					dtc2etaphitpc[aniter]->SetNsigma(nSigmaVal);
					//dtc3etaphitpc[aniter]->SetNsigma(3.0);

					dtc1etaphitpc[aniter]->SetCharge(1.0);
					dtc2etaphitpc[aniter]->SetCharge(-1.0);

					dtc1etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
					dtc2etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);


					dtc1etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);	
					dtc2etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);


                                        if (ichg == 0 ||ichg == 1 ||ichg == 2)//protons 0-2
                                          {
                                            dtc1etaphitpc[aniter]->SetPt(0.5,4);
                                            dtc2etaphitpc[aniter]->SetPt(0.5,4);
                                           
					    dtc1etaphitpc[aniter]->SetMass(ProtonMass);		
					    dtc1etaphitpc[aniter]->SetMostProbableProton();
					    dtc2etaphitpc[aniter]->SetMass(ProtonMass);		
					    dtc2etaphitpc[aniter]->SetMostProbableProton();
                                          }

					if (ichg == 3 ||ichg == 4 ||ichg == 5)//kaons 3-5
                                          {
                                            dtc1etaphitpc[aniter]->SetPt(0.3,4);
                                            dtc2etaphitpc[aniter]->SetPt(0.3,4);
                                        
					    dtc1etaphitpc[aniter]->SetMass(KaonMass);
					    dtc1etaphitpc[aniter]->SetMostProbableKaon();
					    dtc2etaphitpc[aniter]->SetMass(KaonMass);
					    dtc2etaphitpc[aniter]->SetMostProbableKaon();

                                          }
                                        if (ichg == 6 ||ichg == 7 ||ichg == 8)//pions 6-8
                                          {
                                            dtc1etaphitpc[aniter]->SetPt(0.2,4);
					    dtc2etaphitpc[aniter]->SetPt(0.2,4);

					    dtc1etaphitpc[aniter]->SetMass(PionMass);		
					    dtc1etaphitpc[aniter]->SetMostProbablePion();
					    dtc2etaphitpc[aniter]->SetMass(PionMass);		
					    dtc2etaphitpc[aniter]->SetMostProbablePion();
                                          }
                                        if (ichg == 9)//all
                                          {

					    dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();
					    dtc3etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
					    dtc3etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
					    dtc3etaphitpc[aniter]->SetElectronRejection(ifElectronRejection); 
                                            dtc3etaphitpc[aniter]->SetPt(0.2,4);
                                          }
                                        if (ichg == 10 ||ichg == 11 ||ichg == 12)//plus,minus,mixed
                                          {
                                            dtc1etaphitpc[aniter]->SetPt(0.2,4);
                                            dtc2etaphitpc[aniter]->SetPt(0.2,4);
                                          }

					//****** DCA ******

					if(owndca){
					  dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); 	//	DCA xy
					  //dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
					  dtc1etaphitpc[aniter]->SetMaxImpactZ(2);	//DCA Z
					  dtc2etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); 	//	DCA xy
					  dtc2etaphitpc[aniter]->SetMaxImpactZ(2);	//DCA Z
					  if (ichg == 9){dtc3etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); 	//	DCA xy
					  dtc3etaphitpc[aniter]->SetMaxImpactZ(2);}	//DCA Z
					}
					//****** Track quality cuts ******

					if(owncuts){

					  dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
					  dtc1etaphitpc[aniter]->SetminTPCncls(70);
					  dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
					  dtc1etaphitpc[aniter]->SetLabel(kFALSE);
					  //	dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
					  dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);	// pisac
					  //dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);

					  dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
					  dtc2etaphitpc[aniter]->SetminTPCncls(70);
					  dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
					  dtc2etaphitpc[aniter]->SetLabel(kFALSE);
					  //	dtc2etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
					  dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
					  //	dtc2etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
					  if (ichg == 9){
					    dtc3etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit|AliESDtrack::kITSrefit);
					    dtc3etaphitpc[aniter]->SetminTPCncls(70);
					    dtc3etaphitpc[aniter]->SetRemoveKinks(kTRUE);
					    dtc3etaphitpc[aniter]->SetLabel(kFALSE);
					    //	dtc3etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
					    dtc3etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
					    //	dtc3etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
					  }
	
					}
					//**************** track Monitors ***************

					
					if(1)//ichg>8)
					  {
					    if(ichg<2 || ichg==3||ichg==4 || ichg==6|| ichg==7||ichg==9||ichg==10||ichg==11){ 
					      cutPass3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass%stpcM%i", chrgs[ichg], imult),PionMass);
					      cutFail3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail%stpcM%i", chrgs[ichg], imult),PionMass);
					    }
					    if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
					    if(ichg==0||ichg==3||ichg==6||ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
					    if(ichg==1||ichg==4||ichg==7||ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);

					    if(ichg<2){ //PP, PaP
					      cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
					      cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),2);
					      if(ichg==0) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					      if(ichg==1) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
		    			    }
					    else if(ichg>=3 && ichg<=4){//KpKp, KmKm
					      cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),1);//0-pion,1-kaon,2-proton
					      cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),1);
					      if(ichg==3) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					      if(ichg==4) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					    }
					    else if(ichg==6||ichg==10||ichg==7||ichg==11||ichg==9){ //pions, all
					      cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),0);//0-pion,1-kaon,2-proton
					      cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),0);
					      if(ichg==6||ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					      if(ichg==7||ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					      if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					    }
					  }
					 
					//******** Two - track cuts ************
					//sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
					sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
				
					sqpcetaphitpc[aniter]->SetShareQualityMax(shareQuality);		// two track cuts on splitting and merging  //1- wylaczany 0 -wlaczany   
					sqpcetaphitpc[aniter]->SetShareFractionMax(shareFraction);	//  ile moga miec wspolnych klastrow //1 - wylaczany, 0.05 - wlaczany
					sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
					// sqpcetaphitpc[aniter]->SetMaximumRadius(0.82);
					// sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
					// sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.02);
					// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
					/*
					if (gammacut == 0)
					  {
					    sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
					    sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
					  }
					else if (gammacut == 1)
					  {
					    sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
					    sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
					  }
					*/
					// sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
					// sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
					// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
					  
					
					//***** Setting cuts ***********

		
					// setting event cut
					anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
					//setting single track cuts
					if(ichg==0 || ichg==3  || ichg==6 || ichg==10) //positive like-sign
					{
						anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
						anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
					}
					if(ichg==1 || ichg==4 || ichg==7 || ichg==11)//negative like-sign
					{
						anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
						anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
					}
					if(ichg==2 || ichg==5 || ichg==8 || ichg==12)//unlike-sign
					{
						anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
						anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
					}
					if(ichg==9) //all
                                        {
					  anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
					  anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
                                        }

					//setting two-track cuts
					anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);


					//**** Correlation functions *******

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

		
					if (runktdep)
					{

					  if(ichg<=2){
					  for(int kit=0;kit<=numOfkTbins;kit++)
					    ktrng[kit]=ktrngProton[kit];
					  }
					  else if(ichg>2&&ichg<6){
					    for(int kit=0;kit<=numOfkTbins;kit++)
					      ktrng[kit]=ktrngKaon[kit];
                                          }
					  else if(ichg>=6&&ichg<=8){
					    for(int kit=0;kit<=numOfkTbins;kit++)
					      ktrng[kit]=ktrngPion[kit];
                                          }
					  else if(ichg>=9){
					    for(int kit=0;kit<=numOfkTbins;kit++)
					      ktrng[kit]=ktrngAll[kit];
                                          }


						int ktm;
						for (int ikt=0; ikt<numOfkTbins; ikt++)
						{
						  	if(ktrng[ikt+1]==0) continue;
							ktm = aniter * numOfkTbins + ikt;
							ktpcuts[ktm] = new AliFemtoPairCutPt(ktrng[ikt], ktrng[ikt+1]);
				
							//cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
							//cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
							//anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);

							cdedpetaphiPt[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ipT%i", chrgs[ichg], imult,ikt),35, 35);
							cdedpetaphiPt[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
							anetaphitpc[aniter]->AddCorrFctn(cdedpetaphiPt[ktm]);

						}
					}		
					Manager->AddAnalysis(anetaphitpc[aniter]);	
				}
			}
		}
	}
	return Manager;
}
Commit	Line	Data
9ee51cbc	1	/*********************************************************************
	2	* *
	3	* Configfemtoanalysis.C - configuration macro for the femtoscopic *
	4	* analysis, meant as a QA process for two-particle effects *
	5	* *
	6	* Author: Adam Kisiel (Adam.Kisiel@cern.ch) *
	7	* *
	8	*********************************************************************/
	9
	10	#if !defined(__CINT__) \|\| defined(__MAKECINT_)
	11	#include "AliFemtoManager.h"
	12	#include "AliFemtoEventReaderESDChain.h"
	13	#include "AliFemtoEventReaderESDChainKine.h"
	14	#include "AliFemtoEventReaderAODChain.h"
	15	#include "AliFemtoSimpleAnalysis.h"
	16	#include "AliFemtoBasicEventCut.h"
	17	#include "AliFemtoESDTrackCut.h"
	18	#include "AliFemtoCorrFctn.h"
	19	#include "AliFemtoCutMonitorParticleYPt.h"
	20	#include "AliFemtoCutMonitorParticleVertPos.h"
	21	#include "AliFemtoCutMonitorParticleMomRes.h"
	22	#include "AliFemtoCutMonitorParticlePID.h"
	23	#include "AliFemtoCutMonitorEventMult.h"
	24	#include "AliFemtoCutMonitorEventVertex.h"
	25	#include "AliFemtoShareQualityTPCEntranceSepPairCut.h"
	26	#include "AliFemtoPairCutAntiGamma.h"
	27	#include "AliFemtoPairCutRadialDistance.h"
	28	#include "AliFemtoQinvCorrFctn.h"
	29	#include "AliFemtoCorrFctnNonIdDR.h"
	30	#include "AliFemtoShareQualityCorrFctn.h"
	31	#include "AliFemtoTPCInnerCorrFctn.h"
	32	#include "AliFemtoVertexMultAnalysis.h"
	33	#include "AliFemtoCorrFctn3DSpherical.h"
	34	#include "AliFemtoChi2CorrFctn.h"
	35	#include "AliFemtoCorrFctnTPCNcls.h"
	36	#include "AliFemtoBPLCMS3DCorrFctn.h"
	37	#include "AliFemtoCorrFctn3DLCMSSym.h"
	38	#include "AliFemtoModelBPLCMSCorrFctn.h"
	39	#include "AliFemtoModelCorrFctn3DSpherical.h"
	40	#include "AliFemtoModelGausLCMSFreezeOutGenerator.h"
	41	#include "AliFemtoModelGausRinvFreezeOutGenerator.h"
	42	#include "AliFemtoModelManager.h"
	43	#include "AliFemtoModelWeightGeneratorBasic.h"
	44	#include "AliFemtoModelWeightGeneratorLednicky.h"
	45	#include "AliFemtoCorrFctnDirectYlm.h"
	46	#include "AliFemtoModelCorrFctnDirectYlm.h"
	47	#include "AliFemtoModelCorrFctnSource.h"
	48	#include "AliFemtoCutMonitorParticlePtPDG.h"
	49	#include "AliFemtoKTPairCut.h"
	50	#include "AliFemtoPairCutPt.h"
	51	#endif
	52
	53	//________________________________________________________________________
	54	AliFemtoManager* ConfigFemtoAnalysis(const char* params="") {
	55
	56	double PionMass = 0.13956995;
	57	double KaonMass = 0.493677;
	58	double ProtonMass = 0.938272013;
	59
	60	const int numOfMultBins = 5;
	61	const int numOfChTypes = 13;
	62	const int numOfkTbins = 6;
	63
	64	char *parameter[15];
65	if(strlen(params)!=0)
66	{
67	parameter[0] = strtok(params, ","); // Splits spaces between words in params
68	cout<<"Parameter [0] (filterbit):"<<parameter[0]<<endl; // Writes first parameter
69	parameter[1] = strtok(NULL, ",");
70	cout<<"Parameter [1] (ktdep):"<<parameter[1]<<" "<<endl;
71	parameter[2] = strtok(NULL, ",");
72	cout<<"Parameter [2] (multdep):"<<parameter[2]<<" "<<endl;
73	parameter[3] = strtok(NULL, ",");
74	cout<<"Parameter [3]: (MinPlpContribSPD)"<<parameter[3]<<" "<<endl;
75	parameter[4] = strtok(NULL, ",");
76	cout<<"Parameter [4]: (multbino)"<<parameter[4]<<" "<<endl;
77	parameter[5] = strtok(NULL, ",");
78	cout<<"Parameter [5]: (zvertbino)"<<parameter[5]<<" "<<endl;
79	parameter[6] = strtok(NULL, ",");
80	cout<<"Parameter [6]: (ifGlobalTracks=true/false)"<<parameter[6]<<" "<<endl;
81	parameter[7] = strtok(NULL, ",");
82	cout<<"Parameter [7]: (shareQuality)"<<parameter[7]<<" "<<endl;
83	parameter[8] = strtok(NULL, ",");
84	cout<<"Parameter [8]: (shareFraction)"<<parameter[8]<<" "<<endl;
85	parameter[9] = strtok(NULL, ",");
86	cout<<"Parameter [9]: (ifElectronRejection)"<<parameter[9]<<" "<<endl;
87	parameter[10] = strtok(NULL, ",");
88	cout<<"Parameter [10]: (nSigma)"<<parameter[10]<<" "<<endl;
89	parameter[11] = strtok(NULL, ",");
90	cout<<"Parameter [12]: (etaMin)"<<parameter[11]<<" "<<endl;
91	parameter[12] = strtok(NULL, ",");
92	cout<<"Parameter [12]: (etaMax)"<<parameter[12]<<" "<<endl;
93	parameter[13] = strtok(NULL, ",");
94	cout<<"Parameter [13]: (ispileup)"<<parameter[13]<<" "<<endl;
95	}
96	int filterbit = atoi(parameter[0]); //96 / 768 / 128
97	int runktdep = atoi(parameter[1]); //0
98	int runmultdep = atoi(parameter[2]); //0
99	int minPlpContribSPD = atoi(parameter[3]); //3
100	int multbino = atoi(parameter[4]); //30
101	int zvertbino = atoi(parameter[5]); //10
102	Bool_t ifGlobalTracks=kFALSE; if(atoi(parameter[6]))ifGlobalTracks=kTRUE;//kTRUE
103	double shareQuality = atof(parameter[7]); //0.00
104	double shareFraction = atof(parameter[8]); //0.05
105	bool ifElectronRejection = atoi(parameter[9]); //true
106	double nSigmaVal = atof(parameter[10]); //3.0
107	double nEtaMin = atof(parameter[11]); //-0.8
108	double nEtaMax = atof(parameter[12]); //0.8
109	bool ifIsPileUp = atoi(parameter[13]); //true
110
111	int runmults[numOfMultBins] = {0, 0, 0, 0, 1};
112	if(runmultdep) {runmults[0]=1; runmults[1]=1; runmults[2]=1; }
113	int multbins[numOfMultBins+1] = {2, 20, 50,150,2,150};
114
115	int runch[numOfChTypes] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
116	const char *chrgs[numOfChTypes] = { "PP", "aPaP", "PaP", "KpKp", "KmKm", "KpKm", "PIpPIp", "PImPIm", "PIpPIm", "all", "plus", "minus", "mixed" };
117
118	int runktdep = 1;
119	double ktrng[numOfkTbins+1] = {0.0, 0, 0, 0, 0, 0};
120	double ktrngAll[numOfkTbins+1] = {0.0, 7.0, 1.2, 2.0, 3.0, 4.0, 100.0};
121	double ktrngPion[numOfkTbins+1] = {0.0, 0.7, 0.9, 1.1, 1.4, 2.5, 100.0};
122	double ktrngKaon[numOfkTbins+1] = {0.0, 1.5, 2.5, 3.5, 100.0, 0, 0};
123	double ktrngProton[numOfkTbins+1] = {0.0, 2.75, 100, 0, 0};
124
125	int runqinv = 1;
126	int runshlcms = 1;// 0:PRF(PAP), 1:LCMS(PP,APAP)
127
128	int runtype = 0; // Types 0 - global, 1 - ITS only, 2 - TPC Inner //global tracks ->mfit ITS+TPC
129	int owncuts = 0;
130	int owndca = 0;
131
132	int gammacut = 0; // cut na ee z gamma
133
134	double shqmax = 0.5;
135	int nbinssh = 100;
136
137	//AliFemtoEventReaderESDChain *Reader = new AliFemtoEventReaderESDChain();
138	//Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kGlobalCount);
139
140
141	AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
142	Reader->SetFilterMask(filterbit);
143	Reader->SetDCAglobalTrack(ifGlobalTracks); //false for FB7, true for the rest //we do not use DCA at all
144	Reader->SetUseMultiplicity(AliFemtoEventReaderAOD::kReference);
145	Reader->SetMinPlpContribSPD(minPlpContribSPD);
146	Reader->SetIsPileUpEvent(ifIsPileUp);
147
148	AliFemtoManager* Manager = new AliFemtoManager();
149	Manager->SetEventReader(Reader);
150
151	AliFemtoVertexMultAnalysis *anetaphitpc[320];
152	AliFemtoBasicEventCut *mecetaphitpc[320];
153	AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[320];
154	AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[320];
155	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[320];
156	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[320];
157	AliFemtoESDTrackCut *dtc1etaphitpc[320];
158	AliFemtoESDTrackCut *dtc2etaphitpc[320];
159	AliFemtoESDTrackCut *dtc3etaphitpc[320];
160	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[320];
161	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[320];
162	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[320];
163	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[320];
164	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[320];
165	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[320];
166	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[320];
167	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[320];
168	AliFemtoCutMonitorParticleYPt *cutPass3YPtetaphitpc[320];
169	AliFemtoCutMonitorParticleYPt *cutFail3YPtetaphitpc[320];
170	AliFemtoCutMonitorParticlePID *cutPass3PIDetaphitpc[320];
171	AliFemtoCutMonitorParticlePID *cutFail3PIDetaphitpc[320];
172	// AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpcsame[320];
173	//AliFemtoPairCutAntiGamma *sqpcetaphitpc[320];
174	//AliFemtoPairCutRadialDistance *sqpcetaphitpc[320];
175	AliFemtoShareQualityPairCut *sqpcetaphitpc[320];
176	// AliFemtoChi2CorrFctn *cchiqinvetaphitpc[320];
177	AliFemtoPairCutPt *ktpcuts[1020];
178	AliFemtoQinvCorrFctn *cqinvkttpc[320];
179	AliFemtoQinvCorrFctn *cqinvtpc[320];
180	AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[320];
181	AliFemtoCorrFctnDEtaDPhi *cdedpetaphiPt[1320];
182
183
184
185	// * Third QA task - HBT analysis with all pair cuts off, TPC only *
186	// * Begin pion-pion (positive) analysis *
187	int aniter = 0;
188
189	for (int imult = 0; imult < numOfMultBins; imult++)
190	{
191	if (runmults[imult])
192	{
193	for (int ichg = 0; ichg < numOfChTypes; ichg++)
194	{
195	if (runch[ichg])
196	{
197
198	aniter = ichg * numOfMultBins + imult;
199	anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(zvertbino, -10.0, 10.0, multbino, multbins[imult], multbins[imult+1]);
200	anetaphitpc[aniter]->SetNumEventsToMix(10);
201	anetaphitpc[aniter]->SetMinSizePartCollection(1);
41fa2a11	202	anetaphitpc[aniter]->SetVerboseMode(kFALSE);//~~~~~~~~~~~~~~~~
9ee51cbc	203
	204	//* Event cut *
	205	mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
	206	mecetaphitpc[aniter]->SetEventMult(0.001,100000);
	207	mecetaphitpc[aniter]->SetVertZPos(-10,10);//cm
	208
	209	//**** event monitors ********
	210	cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	211	cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	212	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
	213
	214	cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	215	cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	216	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
	217
	218
	219	// *** single particle track cuts *******
	220	dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
	221	dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
	222	dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	223	dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	224	dtc1etaphitpc[aniter]->SetNsigma(nSigmaVal);
	225	dtc2etaphitpc[aniter]->SetNsigma(nSigmaVal);
	226	//dtc3etaphitpc[aniter]->SetNsigma(3.0);
	227
	228	dtc1etaphitpc[aniter]->SetCharge(1.0);
	229	dtc2etaphitpc[aniter]->SetCharge(-1.0);
	230
	231	dtc1etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
	232	dtc2etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
	233
	234
	235	dtc1etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
	236	dtc2etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
	237
	238
	239
	240	if (ichg == 0 \|\|ichg == 1 \|\|ichg == 2)//protons 0-2
	241	{
	242	dtc1etaphitpc[aniter]->SetPt(0.5,4);
	243	dtc2etaphitpc[aniter]->SetPt(0.5,4);
	244
	245	dtc1etaphitpc[aniter]->SetMass(ProtonMass);
	246	dtc1etaphitpc[aniter]->SetMostProbableProton();
	247	dtc2etaphitpc[aniter]->SetMass(ProtonMass);
	248	dtc2etaphitpc[aniter]->SetMostProbableProton();
	249	}
	250
	251	if (ichg == 3 \|\|ichg == 4 \|\|ichg == 5)//kaons 3-5
	252	{
	253	dtc1etaphitpc[aniter]->SetPt(0.3,4);
	254	dtc2etaphitpc[aniter]->SetPt(0.3,4);
	255
	256	dtc1etaphitpc[aniter]->SetMass(KaonMass);
	257	dtc1etaphitpc[aniter]->SetMostProbableKaon();
	258	dtc2etaphitpc[aniter]->SetMass(KaonMass);
	259	dtc2etaphitpc[aniter]->SetMostProbableKaon();
	260
	261	}
	262	if (ichg == 6 \|\|ichg == 7 \|\|ichg == 8)//pions 6-8
	263	{
	264	dtc1etaphitpc[aniter]->SetPt(0.2,4);
	265	dtc2etaphitpc[aniter]->SetPt(0.2,4);
	266
267	dtc1etaphitpc[aniter]->SetMass(PionMass);
268	dtc1etaphitpc[aniter]->SetMostProbablePion();
269	dtc2etaphitpc[aniter]->SetMass(PionMass);
270	dtc2etaphitpc[aniter]->SetMostProbablePion();
271	}
272	if (ichg == 9)//all
273	{
274
275	dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();
276	dtc3etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
277	dtc3etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
278	dtc3etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
279	dtc3etaphitpc[aniter]->SetPt(0.2,4);
280	}
281	if (ichg == 10 \|\|ichg == 11 \|\|ichg == 12)//plus,minus,mixed
282	{
283	dtc1etaphitpc[aniter]->SetPt(0.2,4);
284	dtc2etaphitpc[aniter]->SetPt(0.2,4);
285	}
286
287	//**** DCA ****
288
289	if(owndca){
290	dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
291	//dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
292	dtc1etaphitpc[aniter]->SetMaxImpactZ(2); //DCA Z
293	dtc2etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
294	dtc2etaphitpc[aniter]->SetMaxImpactZ(2); //DCA Z
295	if (ichg == 9){dtc3etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
296	dtc3etaphitpc[aniter]->SetMaxImpactZ(2);} //DCA Z
297	}
298	//**** Track quality cuts ****
299
300	if(owncuts){
301
302	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
303	dtc1etaphitpc[aniter]->SetminTPCncls(70);
304	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
305	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
306	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
307	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0); // pisac
308	//dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
309
310	dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
311	dtc2etaphitpc[aniter]->SetminTPCncls(70);
312	dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
313	dtc2etaphitpc[aniter]->SetLabel(kFALSE);
314	// dtc2etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
315	dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
316	// dtc2etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
317	if (ichg == 9){
318	dtc3etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
319	dtc3etaphitpc[aniter]->SetminTPCncls(70);
320	dtc3etaphitpc[aniter]->SetRemoveKinks(kTRUE);
321	dtc3etaphitpc[aniter]->SetLabel(kFALSE);
322	// dtc3etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
323	dtc3etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
324	// dtc3etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
325	}
326
327	}
328	//************** track Monitors *************
329
330
331	if(1)//ichg>8)
332	{
333	if(ichg<2 \|\| ichg==3\|\|ichg==4 \|\| ichg==6\|\| ichg==7\|\|ichg==9\|\|ichg==10\|\|ichg==11){
334	cutPass3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass%stpcM%i", chrgs[ichg], imult),PionMass);
335	cutFail3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail%stpcM%i", chrgs[ichg], imult),PionMass);
336	}
337	if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
338	if(ichg==0\|\|ichg==3\|\|ichg==6\|\|ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
339	if(ichg==1\|\|ichg==4\|\|ichg==7\|\|ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
340
341	if(ichg<2){ //PP, PaP
342	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
343	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),2);
344	if(ichg==0) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
345	if(ichg==1) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
346	}
347	else if(ichg>=3 && ichg<=4){//KpKp, KmKm
348	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),1);//0-pion,1-kaon,2-proton
349	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),1);
350	if(ichg==3) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
351	if(ichg==4) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
352	}
353	else if(ichg==6\|\|ichg==10\|\|ichg==7\|\|ichg==11\|\|ichg==9){ //pions, all
354	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),0);//0-pion,1-kaon,2-proton
355	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),0);
356	if(ichg==6\|\|ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
357	if(ichg==7\|\|ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
358	if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
359	}
360	}
361
362	//****** Two - track cuts **********
363	//sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
364	sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
365
366	sqpcetaphitpc[aniter]->SetShareQualityMax(shareQuality); // two track cuts on splitting and merging //1- wylaczany 0 -wlaczany
367	sqpcetaphitpc[aniter]->SetShareFractionMax(shareFraction); // ile moga miec wspolnych klastrow //1 - wylaczany, 0.05 - wlaczany
368	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
369	// sqpcetaphitpc[aniter]->SetMaximumRadius(0.82);
370	// sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
371	// sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.02);
372	// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
373	/*
374	if (gammacut == 0)
375	{
376	sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
377	sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
378	}
379	else if (gammacut == 1)
380	{
381	sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
382	sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
383	}
384	*/
385	// sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
386	// sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
387	// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
388
389
390
391
392	//*** Setting cuts *********
393
394
395	// setting event cut
396	anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
397	//setting single track cuts
398	if(ichg==0 \|\| ichg==3 \|\| ichg==6 \|\| ichg==10) //positive like-sign
399	{
400	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
401	anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
402	}
403	if(ichg==1 \|\| ichg==4 \|\| ichg==7 \|\| ichg==11)//negative like-sign
404	{
405	anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
406	anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
407	}
408	if(ichg==2 \|\| ichg==5 \|\| ichg==8 \|\| ichg==12)//unlike-sign
409	{
410	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
411	anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
412	}
413	if(ichg==9) //all
414	{
415	anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
416	anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
417	}
418
419	//setting two-track cuts
420	anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
421
422
423	//** Correlation functions *****
424
425	cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),35, 35);
426	anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
427
428
429	if (runktdep)
430	{
431
432	if(ichg<=2){
433	for(int kit=0;kit<=numOfkTbins;kit++)
434	ktrng[kit]=ktrngProton[kit];
435	}
436	else if(ichg>2&&ichg<6){
437	for(int kit=0;kit<=numOfkTbins;kit++)
438	ktrng[kit]=ktrngKaon[kit];
439	}
440	else if(ichg>=6&&ichg<=8){
441	for(int kit=0;kit<=numOfkTbins;kit++)
442	ktrng[kit]=ktrngPion[kit];
443	}
444	else if(ichg>=9){
445	for(int kit=0;kit<=numOfkTbins;kit++)
446	ktrng[kit]=ktrngAll[kit];
447	}
448
449
450	int ktm;
451	for (int ikt=0; ikt<numOfkTbins; ikt++)
452	{
453	if(ktrng[ikt+1]==0) continue;
454	ktm = aniter * numOfkTbins + ikt;
455	ktpcuts[ktm] = new AliFemtoPairCutPt(ktrng[ikt], ktrng[ikt+1]);
456
457	//cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
458	//cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
459	//anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
460
461	cdedpetaphiPt[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ipT%i", chrgs[ichg], imult,ikt),35, 35);
462	cdedpetaphiPt[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
463	anetaphitpc[aniter]->AddCorrFctn(cdedpetaphiPt[ktm]);
464
465	}
466	}
467	Manager->AddAnalysis(anetaphitpc[aniter]);
468	}
469	}
470	}
471	}
472	return Manager;
473	}