From 21376f5972c261d6a424352525c910a83a77fdf9 Mon Sep 17 00:00:00 2001 From: majanik Date: Sat, 19 Jul 2014 15:47:52 +0200 Subject: [PATCH] Untriggered DEtaDPhi: New Configs for the train --- .../NoCorrections/ConfigFemtoAnalysis.C | 9 +- .../ConfigFemtoAnalysis.C | 553 ++++++++++++++++++ 2 files changed, 560 insertions(+), 2 deletions(-) create mode 100644 PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train7Light/CorrectionsParameters/ConfigFemtoAnalysis.C diff --git a/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train6Light/NoCorrections/ConfigFemtoAnalysis.C b/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train6Light/NoCorrections/ConfigFemtoAnalysis.C index 163d772eca0..4daf4d5552b 100644 --- a/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train6Light/NoCorrections/ConfigFemtoAnalysis.C +++ b/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train6Light/NoCorrections/ConfigFemtoAnalysis.C @@ -92,6 +92,8 @@ AliFemtoManager* ConfigFemtoAnalysis(const char* params) { cout<<"Parameter [12]: (etaMax)"<SetPt(0.3,4); - dtc2etaphitpc[aniter]->SetPt(0.3,4); + dtc1etaphitpc[aniter]->SetPt(0.3,maxPtKaons); + dtc2etaphitpc[aniter]->SetPt(0.3,maxPtKaons); dtc1etaphitpc[aniter]->SetMass(KaonMass); dtc1etaphitpc[aniter]->SetMostProbable(11);//cut on Nsigma in pT not p dtc2etaphitpc[aniter]->SetMass(KaonMass); diff --git a/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train7Light/CorrectionsParameters/ConfigFemtoAnalysis.C b/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train7Light/CorrectionsParameters/ConfigFemtoAnalysis.C new file mode 100644 index 00000000000..737c5ab49f0 --- /dev/null +++ b/PWGCF/FEMTOSCOPY/macros/Train/DEtaDPhi/2014/Train7Light/CorrectionsParameters/ConfigFemtoAnalysis.C @@ -0,0 +1,553 @@ +/********************************************************************* * * + * 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 "AliFemtoMJTrackCut.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 "AliFemtoCorrFctnDEtaDPhiCorrections.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 = 5; + + char *parameter[20]; + if(strlen(params)!=0) + { + parameter[0] = strtok(params, ","); // Splits spaces between words in params + cout<<"Parameter [0] (filterbit):"<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[numOfMultBins*numOfChTypes]; + AliFemtoBasicEventCut *mecetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorEventMult *cutPassEvMetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorEventMult *cutFailEvMetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorEventVertex *cutPassEvVetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorEventVertex *cutFailEvVetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoMJTrackCut *dtc1etaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoMJTrackCut *dtc2etaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoMJTrackCut *dtc3etaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutPass1YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutFail1YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutPass1PIDetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutFail1PIDetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutPass2YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutFail2YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutPass2PIDetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutFail2PIDetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutPass3YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticleYPt *cutFail3YPtetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutPass3PIDetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoCutMonitorParticlePID *cutFail3PIDetaphitpc[numOfMultBins*numOfChTypes]; + // AliFemtoShareQualityTPCEntranceSepPairCut *sqpcetaphitpcsame[numOfMultBins*numOfChTypes]; + //AliFemtoPairCutAntiGamma *sqpcetaphitpc[numOfMultBins*numOfChTypes]; + //AliFemtoPairCutRadialDistance *sqpcetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoShareQualityPairCut *sqpcetaphitpc[numOfMultBins*numOfChTypes]; + // AliFemtoChi2CorrFctn *cchiqinvetaphitpc[numOfMultBins*numOfChTypes]; + AliFemtoPairCutPt *ktpcuts[numOfMultBins*numOfChTypes*numOfkTbins]; + AliFemtoQinvCorrFctn *cqinvkttpc[numOfMultBins*numOfChTypes*numOfkTbins]; + AliFemtoQinvCorrFctn *cqinvtpc[numOfMultBins*numOfChTypes]; + AliFemtoCorrFctnDEtaDPhiCorrections *cdedpetaphi[numOfMultBins*numOfChTypes]; + AliFemtoCorrFctnDEtaDPhiCorrections *cdedpetaphiPt[numOfMultBins*numOfChTypes*numOfkTbins]; + + + + // *** 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 AliFemtoMJTrackCut(); + dtc2etaphitpc[aniter] = new AliFemtoMJTrackCut(); + 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]->SetMostProbable(setMostProb1);//cut on Nsigma in pT not p + dtc2etaphitpc[aniter]->SetMass(ProtonMass); + dtc2etaphitpc[aniter]->SetMostProbable(setMostProb1);//cut on Nsigma in pT not p + } + + if (ichg == 3 ||ichg == 4 ||ichg == 5)//kaons 3-5 + { + dtc1etaphitpc[aniter]->SetPt(0.3,maxPtKaons); + dtc2etaphitpc[aniter]->SetPt(0.3,maxPtKaons); + dtc1etaphitpc[aniter]->SetMass(KaonMass); + dtc1etaphitpc[aniter]->SetMostProbable(setMostProb2);//cut on Nsigma in pT not p + dtc2etaphitpc[aniter]->SetMass(KaonMass); + dtc2etaphitpc[aniter]->SetMostProbable(setMostProb2);//cut on Nsigma in pT not p + + } + 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]->SetMostProbable(setMostProb3);//cut on Nsigma in pT not p + dtc2etaphitpc[aniter]->SetMass(PionMass); + dtc2etaphitpc[aniter]->SetMostProbable(setMostProb3);//cut on Nsigma in pT not p + } + if (ichg == 9)//all + { + + dtc3etaphitpc[aniter] = new AliFemtoMJTrackCut(); + 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) + { + //FULL + if(ichg<2 || ichg==3||ichg==4 || ichg==6|| ichg==7||ichg==9||ichg==10||ichg==11){ + //if(ichg==0 || ichg==3 || ichg==6 || ichg==10){ + 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 AliFemtoCorrFctnDEtaDPhiCorrections(Form("cdedp%stpcM%i", chrgs[ichg], imult),35, 35); + cdedpetaphi[aniter]->SetDoFullAnalysis(kFALSE); + if(ichg==0 || ichg==1 || ichg==2) + cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kProton, AliFemtoCorrFctnDEtaDPhiCorrections::kProton/*,1,1,1,0*/); + //cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kProton, AliFemtoCorrFctnDEtaDPhiCorrections::kProton,1,1,1,1); + else if(ichg==3 || ichg==4 || ichg==5) + cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kKaon, AliFemtoCorrFctnDEtaDPhiCorrections::kKaon/*,1,1,1,0*/); + //cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kKaon, AliFemtoCorrFctnDEtaDPhiCorrections::kKaon,1,1,1,1); + + else if(ichg==6 || ichg==7 || ichg==8) + cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kPion, AliFemtoCorrFctnDEtaDPhiCorrections::kPion/*,1,1,1,0*/); + //cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kPion, AliFemtoCorrFctnDEtaDPhiCorrections::kPion,1,1,1,1); + + else + cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kAll, AliFemtoCorrFctnDEtaDPhiCorrections::kAll/*,1,1,1,0*/); + //cdedpetaphi[aniter]->LoadCorrectionTabFromROOTFile(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kAll, AliFemtoCorrFctnDEtaDPhiCorrections::kAll,1,1,1,1); + + + + + 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; ikt6)?shqmax*2.5:shqmax); + //cqinvkttpc[ktm]->SetPairSelectionCut(ktpcuts[ktm]); + //anetaphitpc[aniter]->AddCorrFctn(cqinvkttpc[ktm]); + + cdedpetaphiPt[ktm] = new AliFemtoCorrFctnDEtaDPhiCorrections(Form("cdedp%stpcM%ipT%i", chrgs[ichg], imult,ikt),35, 35); + cdedpetaphiPt[ktm]->SetPairSelectionCut(ktpcuts[ktm]); + + if(ichg==0 || ichg==1 || ichg==2) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kProton, AliFemtoCorrFctnDEtaDPhiCorrections::kProton/*,1,1,1,1*/); + else if(ichg==3 || ichg==4 || ichg==5) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kKaon, AliFemtoCorrFctnDEtaDPhiCorrections::kKaon); + else if(ichg==6 || ichg==7 || ichg==8) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kPion, AliFemtoCorrFctnDEtaDPhiCorrections::kPion); + else + cdedpetaphiPt[ktm]->LoadCorrectionTabFromROOTFile1D(Form("%s",fileName), AliFemtoCorrFctnDEtaDPhiCorrections::kAll, AliFemtoCorrFctnDEtaDPhiCorrections::kAll); + + + /*if(ichg==0 || ichg==1 || ichg==2) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromFile("alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/pTab.txt","alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/protonCorrTab.txt"); + else if(ichg==3 || ichg==4 || ichg==5) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromFile("alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/pTab.txt","alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/kaonCorrTab.txt"); + else if(ichg==6 || ichg==7 || ichg==8) + cdedpetaphiPt[ktm]->LoadCorrectionTabFromFile("alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/pTab.txt","alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/pionCorrTab.txt"); + else + cdedpetaphiPt[ktm]->LoadCorrectionTabFromFile("alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/pTab.txt","alien:///alice/cern.ch/user/l/lgraczyk/2014/DEtaDPhi/CorrectionTables/allCorrTab.txt"); */ + anetaphitpc[aniter]->AddCorrFctn(cdedpetaphiPt[ktm]); + + } + } + Manager->AddAnalysis(anetaphitpc[aniter]); + } + } + } + } + return Manager; +} -- 2.43.0