[u/mrichter/AliRoot.git] / PWGCF / FEMTOSCOPY / macros / Train / DEtaDPhi / 2014 / Train2 / 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 = 3;

	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" };
	

	double ktrng[numOfkTbins+1] = {0.0, 0.75, 2.0, 100.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[320];
	AliFemtoQinvCorrFctn					*cqinvkttpc[320];
	AliFemtoQinvCorrFctn					*cqinvtpc[320];
	AliFemtoCorrFctnDEtaDPhi			*cdedpetaphi[320];


	// *** 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();
					dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();

					dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
					dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
					dtc3etaphitpc[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);
					dtc3etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);

					dtc1etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);	
					dtc2etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
					dtc3etaphitpc[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]->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
					  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);
					
					  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)
					  {
					    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<3){
					      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);
		    			    }
					    else if(ichg>=3 && ichg<=5){
					      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);
					    }
					    else{
					      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==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					    if(ichg==0||ichg==3||ichg==6||ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					    if(ichg==1||ichg==4||ichg==7||ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
					  }
					 
					//******** Two - track cuts ************
					//sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
					//sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
					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(2.51);
					sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
					sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.02);
					sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.2);

					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)
					{
						int ktm;
						for (int ikt=0; ikt<numOfkTbins; ikt++)
						{
							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]);

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

						}
					}		
					Manager->AddAnalysis(anetaphitpc[aniter]);	
				}
			}
		}
	}
	return Manager;
}
Commit	Line	Data
a7f41240	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 = 3;
	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
116	int runch[numOfChTypes] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1};
117	const char *chrgs[numOfChTypes] = { "PP", "aPaP", "PaP", "KpKp", "KmKm", "KpKm", "PIpPIp", "PImPIm", "PIpPIm", "all", "plus", "minus", "mixed" };
118
119
120
121	double ktrng[numOfkTbins+1] = {0.0, 0.75, 2.0, 100.0};
122
123	int runqinv = 1;
124	int runshlcms = 1;// 0:PRF(PAP), 1:LCMS(PP,APAP)
125
126	int runtype = 0; // Types 0 - global, 1 - ITS only, 2 - TPC Inner //global tracks ->mfit ITS+TPC
127	int owncuts = 0;
128	int owndca = 0;
129
130	int gammacut = 0; // cut na ee z gamma
131
132	double shqmax = 0.5;
133	int nbinssh = 100;
134
135	//AliFemtoEventReaderESDChain *Reader = new AliFemtoEventReaderESDChain();
136	//Reader->SetUseMultiplicity(AliFemtoEventReaderESDChain::kGlobalCount);
137
138
139	AliFemtoEventReaderAODChain *Reader = new AliFemtoEventReaderAODChain();
140	Reader->SetFilterMask(filterbit);
141	Reader->SetDCAglobalTrack(ifGlobalTracks); //false for FB7, true for the rest //we do not use DCA at all
142	Reader->SetUseMultiplicity(AliFemtoEventReaderAOD::kReference);
143	Reader->SetMinPlpContribSPD(minPlpContribSPD);
144	Reader->SetIsPileUpEvent(ifIsPileUp);
145
146	AliFemtoManager* Manager = new AliFemtoManager();
147	Manager->SetEventReader(Reader);
148
149	AliFemtoVertexMultAnalysis *anetaphitpc[320];
150	AliFemtoBasicEventCut *mecetaphitpc[320];
151	AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[320];
152	AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[320];
153	AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[320];
154	AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[320];
155	AliFemtoESDTrackCut *dtc1etaphitpc[320];
156	AliFemtoESDTrackCut *dtc2etaphitpc[320];
157	AliFemtoESDTrackCut *dtc3etaphitpc[320];
158	AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[320];
159	AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[320];
160	AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[320];
161	AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[320];
162	AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[320];
163	AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[320];
164	AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[320];
165	AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[320];
166	AliFemtoCutMonitorParticleYPt *cutPass3YPtetaphitpc[320];
167	AliFemtoCutMonitorParticleYPt *cutFail3YPtetaphitpc[320];
168	AliFemtoCutMonitorParticlePID *cutPass3PIDetaphitpc[320];
169	AliFemtoCutMonitorParticlePID *cutFail3PIDetaphitpc[320];
170	// AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpcsame[320];
171	//AliFemtoPairCutAntiGamma *sqpcetaphitpc[320];
172	//AliFemtoPairCutRadialDistance *sqpcetaphitpc[320];
173	AliFemtoShareQualityPairCut *sqpcetaphitpc[320];
174	// AliFemtoChi2CorrFctn *cchiqinvetaphitpc[320];
175	AliFemtoPairCutPt *ktpcuts[320];
176	AliFemtoQinvCorrFctn *cqinvkttpc[320];
177	AliFemtoQinvCorrFctn *cqinvtpc[320];
178	AliFemtoCorrFctnDEtaDPhi *cdedpetaphi[320];
179
180
181
182	// * Third QA task - HBT analysis with all pair cuts off, TPC only *
183	// * Begin pion-pion (positive) analysis *
184	int aniter = 0;
185
186	for (int imult = 0; imult < numOfMultBins; imult++)
187	{
188	if (runmults[imult])
189	{
190	for (int ichg = 0; ichg < numOfChTypes; ichg++)
191	{
192	if (runch[ichg])
193	{
194	aniter = ichg * numOfMultBins + imult;
195	anetaphitpc[aniter] = new AliFemtoVertexMultAnalysis(zvertbino, -10.0, 10.0, multbino, multbins[imult], multbins[imult+1]);
196	anetaphitpc[aniter]->SetNumEventsToMix(10);
197	anetaphitpc[aniter]->SetMinSizePartCollection(1);
7887fda2	198	anetaphitpc[aniter]->SetVerboseMode(kFALSE);//~~~~~~~~~~~~~~~~
a7f41240	199
	200	//* Event cut *
	201	mecetaphitpc[aniter] = new AliFemtoBasicEventCut();
	202	mecetaphitpc[aniter]->SetEventMult(0.001,100000);
	203	mecetaphitpc[aniter]->SetVertZPos(-10,10);//cm
	204
	205	//**** event monitors ********
	206	cutPassEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	207	cutFailEvMetaphitpc[aniter] = new AliFemtoCutMonitorEventMult(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	208	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvMetaphitpc[aniter], cutFailEvMetaphitpc[aniter]);
	209
	210	cutPassEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutPass%stpcM%i", chrgs[ichg], imult));
	211	cutFailEvVetaphitpc[aniter] = new AliFemtoCutMonitorEventVertex(Form("cutFail%stpcM%i", chrgs[ichg], imult));
	212	mecetaphitpc[aniter]->AddCutMonitor(cutPassEvVetaphitpc[aniter], cutFailEvVetaphitpc[aniter]);
	213
	214
	215	// *** single particle track cuts *******
	216	dtc1etaphitpc[aniter] = new AliFemtoESDTrackCut();
	217	dtc2etaphitpc[aniter] = new AliFemtoESDTrackCut();
	218	dtc3etaphitpc[aniter] = new AliFemtoESDTrackCut();
	219
	220	dtc1etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	221	dtc2etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	222	dtc3etaphitpc[aniter]->SetNsigmaTPCTOF(kTRUE);
	223	dtc1etaphitpc[aniter]->SetNsigma(nSigmaVal);
	224	dtc2etaphitpc[aniter]->SetNsigma(nSigmaVal);
	225	//dtc3etaphitpc[aniter]->SetNsigma(3.0);
	226
	227	dtc1etaphitpc[aniter]->SetCharge(1.0);
	228	dtc2etaphitpc[aniter]->SetCharge(-1.0);
	229
	230	dtc1etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
	231	dtc2etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
	232	dtc3etaphitpc[aniter]->SetEta(nEtaMin,nEtaMax);
	233
	234	dtc1etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
	235	dtc2etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
	236	dtc3etaphitpc[aniter]->SetElectronRejection(ifElectronRejection);
	237
	238	if (ichg == 0 \|\|ichg == 1 \|\|ichg == 2)//protons 0-2
	239	{
	240	dtc1etaphitpc[aniter]->SetPt(0.5,4);
	241	dtc2etaphitpc[aniter]->SetPt(0.5,4);
	242
	243	dtc1etaphitpc[aniter]->SetMass(ProtonMass);
	244	dtc1etaphitpc[aniter]->SetMostProbableProton();
	245	dtc2etaphitpc[aniter]->SetMass(ProtonMass);
	246	dtc2etaphitpc[aniter]->SetMostProbableProton();
	247	}
	248
	249	if (ichg == 3 \|\|ichg == 4 \|\|ichg == 5)//kaons 3-5
	250	{
	251	dtc1etaphitpc[aniter]->SetPt(0.3,4);
	252	dtc2etaphitpc[aniter]->SetPt(0.3,4);
	253
	254	dtc1etaphitpc[aniter]->SetMass(KaonMass);
	255	dtc1etaphitpc[aniter]->SetMostProbableKaon();
	256	dtc2etaphitpc[aniter]->SetMass(KaonMass);
	257	dtc2etaphitpc[aniter]->SetMostProbableKaon();
	258
	259	}
	260	if (ichg == 6 \|\|ichg == 7 \|\|ichg == 8)//pions 6-8
	261	{
	262	dtc1etaphitpc[aniter]->SetPt(0.2,4);
263	dtc2etaphitpc[aniter]->SetPt(0.2,4);
264
265	dtc1etaphitpc[aniter]->SetMass(PionMass);
266	dtc1etaphitpc[aniter]->SetMostProbablePion();
267	dtc2etaphitpc[aniter]->SetMass(PionMass);
268	dtc2etaphitpc[aniter]->SetMostProbablePion();
269	}
270	if (ichg == 9)//all
271	{
272	dtc3etaphitpc[aniter]->SetPt(0.2,4);
273	}
274	if (ichg == 10 \|\|ichg == 11 \|\|ichg == 12)//plus,minus,mixed
275	{
276	dtc1etaphitpc[aniter]->SetPt(0.2,4);
277	dtc2etaphitpc[aniter]->SetPt(0.2,4);
278	}
279
280	//**** DCA ****
281
282	if(owndca){
283	dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
284	//dtc1etaphitpc[aniter]->SetMaxImpactXYPtDep(0.0182, 0.0350, -1.01);
285	dtc1etaphitpc[aniter]->SetMaxImpactZ(2); //DCA Z
286	dtc2etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
287	dtc2etaphitpc[aniter]->SetMaxImpactZ(2); //DCA Z
288	dtc3etaphitpc[aniter]->SetMaxImpactXYPtDep(0.018, 0.035, -1.01); // DCA xy
289	dtc3etaphitpc[aniter]->SetMaxImpactZ(2); //DCA Z
290	}
291	//**** Track quality cuts ****
292
293	if(owncuts){
294
295	dtc1etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
296	dtc1etaphitpc[aniter]->SetminTPCncls(70);
297	dtc1etaphitpc[aniter]->SetRemoveKinks(kTRUE);
298	dtc1etaphitpc[aniter]->SetLabel(kFALSE);
299	// dtc1etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
300	dtc1etaphitpc[aniter]->SetMaxTPCChiNdof(4.0); // pisac
301	//dtc1etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
302
303	dtc2etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
304	dtc2etaphitpc[aniter]->SetminTPCncls(70);
305	dtc2etaphitpc[aniter]->SetRemoveKinks(kTRUE);
306	dtc2etaphitpc[aniter]->SetLabel(kFALSE);
307	// dtc2etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
308	dtc2etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
309	// dtc2etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
310
311	dtc3etaphitpc[aniter]->SetStatus(AliESDtrack::kTPCrefit\|AliESDtrack::kITSrefit);
312	dtc3etaphitpc[aniter]->SetminTPCncls(70);
313	dtc3etaphitpc[aniter]->SetRemoveKinks(kTRUE);
314	dtc3etaphitpc[aniter]->SetLabel(kFALSE);
315	// dtc3etaphitpc[aniter]->SetMaxITSChiNdof(6.0);
316	dtc3etaphitpc[aniter]->SetMaxTPCChiNdof(4.0);
317	// dtc3etaphitpc[aniter]->SetMaxSigmaToVertex(6.0);
318
319	}
320	//************** track Monitors *************
321
322
323	if(1)//ichg>8)
324	{
325	cutPass3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutPass%stpcM%i", chrgs[ichg], imult),PionMass);
326	cutFail3YPtetaphitpc[aniter] = new AliFemtoCutMonitorParticleYPt(Form("cutFail%stpcM%i", chrgs[ichg], imult),PionMass);
327	if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
328	if(ichg==0\|\|ichg==3\|\|ichg==6\|\|ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
329	if(ichg==1\|\|ichg==4\|\|ichg==7\|\|ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3YPtetaphitpc[aniter], cutFail3YPtetaphitpc[aniter]);
330
331	if(ichg<3){
332	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),2);//0-pion,1-kaon,2-proton
333	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),2);
334	}
335	else if(ichg>=3 && ichg<=5){
336	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),1);//0-pion,1-kaon,2-proton
337	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),1);
338	}
339	else{
340	cutPass3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutPass%stpcM%i", chrgs[ichg], imult),0);//0-pion,1-kaon,2-proton
341	cutFail3PIDetaphitpc[aniter] = new AliFemtoCutMonitorParticlePID(Form("cutFail%stpcM%i", chrgs[ichg], imult),0);
342
343	}
344	if(ichg==9) dtc3etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
345	if(ichg==0\|\|ichg==3\|\|ichg==6\|\|ichg==10) dtc1etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
346	if(ichg==1\|\|ichg==4\|\|ichg==7\|\|ichg==11) dtc2etaphitpc[aniter]->AddCutMonitor(cutPass3PIDetaphitpc[aniter], cutFail3PIDetaphitpc[aniter]);
347	}
348
349	//****** Two - track cuts **********
350	//sqpcetaphitpc[aniter] = new AliFemtoPairCutAntiGamma();
351	//sqpcetaphitpc[aniter] = new AliFemtoPairCutRadialDistance();
352	sqpcetaphitpc[aniter] = new AliFemtoShareQualityPairCut();
353
354	sqpcetaphitpc[aniter]->SetShareQualityMax(shareQuality); // two track cuts on splitting and merging //1- wylaczany 0 -wlaczany
355	sqpcetaphitpc[aniter]->SetShareFractionMax(shareFraction); // ile moga miec wspolnych klastrow //1 - wylaczany, 0.05 - wlaczany
356	sqpcetaphitpc[aniter]->SetRemoveSameLabel(kFALSE);
357	/*sqpcetaphitpc[aniter]->SetMaximumRadius(2.51);
358	sqpcetaphitpc[aniter]->SetMinimumRadius(0.8);
359	sqpcetaphitpc[aniter]->SetPhiStarDifferenceMinimum(0.02);
360	sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.2);
361
362	if (gammacut == 0)
363	{
364	sqpcetaphitpc[aniter]->SetMaxEEMinv(0.0);
365	sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.0);
366	}
367	else if (gammacut == 1)
368	{
369	sqpcetaphitpc[aniter]->SetMaxEEMinv(0.002);
370	sqpcetaphitpc[aniter]->SetMaxThetaDiff(0.008);
371	}*/
372
373	// sqpcetaphitpc[aniter]->SetTPCEntranceSepMinimum(1.5);
374	// sqpcetaphitpc[aniter]->SetRadialDistanceMinimum(0.12, 0.03);
375	// sqpcetaphitpc[aniter]->SetEtaDifferenceMinimum(0.02);
376
377
378
379
380	//*** Setting cuts *********
381
382
383	// setting event cut
384	anetaphitpc[aniter]->SetEventCut(mecetaphitpc[aniter]);
385	//setting single track cuts
386	if(ichg==0 \|\| ichg==3 \|\| ichg==6 \|\| ichg==10) //positive like-sign
387	{
388	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
389	anetaphitpc[aniter]->SetSecondParticleCut(dtc1etaphitpc[aniter]);
390	}
391	if(ichg==1 \|\| ichg==4 \|\| ichg==7 \|\| ichg==11)//negative like-sign
392	{
393	anetaphitpc[aniter]->SetFirstParticleCut(dtc2etaphitpc[aniter]);
394	anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
395	}
396	if(ichg==2 \|\| ichg==5 \|\| ichg==8 \|\| ichg==12)//unlike-sign
397	{
398	anetaphitpc[aniter]->SetFirstParticleCut(dtc1etaphitpc[aniter]);
399	anetaphitpc[aniter]->SetSecondParticleCut(dtc2etaphitpc[aniter]);
400	}
401	if(ichg==9) //all
402	{
403	anetaphitpc[aniter]->SetFirstParticleCut(dtc3etaphitpc[aniter]);
404	anetaphitpc[aniter]->SetSecondParticleCut(dtc3etaphitpc[aniter]);
405	}
406
407	//setting two-track cuts
408	anetaphitpc[aniter]->SetPairCut(sqpcetaphitpc[aniter]);
409
410
411	//** Correlation functions *****
412
413	cdedpetaphi[aniter] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%i", chrgs[ichg], imult),35, 35);
414	anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[aniter]);
415
416
417	if (runktdep)
418	{
419	int ktm;
420	for (int ikt=0; ikt<numOfkTbins; ikt++)
421	{
422	ktm = aniter * numOfkTbins + ikt;
423	ktpcuts[ktm] = new AliFemtoPairCutPt(ktrng[ikt], ktrng[ikt+1]);
424
425	//cqinvkttpc[ktm] = new AliFemtoQinvCorrFctn(Form("cqinv%stpcM%ikT%i", chrgs[ichg], imult, ikt),nbinssh,0.0,(imult>6)?shqmax*2.5:shqmax);
426	//cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
427	//anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]);
428
429	cdedpetaphi[ktm] = new AliFemtoCorrFctnDEtaDPhi(Form("cdedp%stpcM%ipT%i", chrgs[ichg], imult,ikt),35, 35);
430	cdedpetaphi[ktm]->SetPairSelectionCut(ktpcuts[ktm]);
431	anetaphitpc[aniter]->AddCorrFctn(cdedpetaphi[ktm]);
432
433	}
434	}
435	Manager->AddAnalysis(anetaphitpc[aniter]);
436	}
437	}
438	}
439	}
440	return Manager;
441	}