update V0-h correlations (Jayarathna, Sandun Y <sypahula@Central.UH.EDU>)

diff --git a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx
index 70f37b29f9be071523df3141b62a95c8d2737790..35a0396f2773b100bba7d2afb16d4c251a36e59a 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx
+++ b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.cxx
@@ -37,10 +37,8 @@
 
 #include "AliLeadingV0Correlation.h"
 
-#define CorrBinsX 240
-#define CorrBinsY 260
-using std::endl;
-using std::cout;
+#define CorrBinsX 24
+#define CorrBinsY 26
 
 
 Double_t PI =TMath::Pi();
@@ -61,6 +59,7 @@ AliLeadingV0Correlation::AliLeadingV0Correlation()
        fNzVtxBins(0), 
        fNCentBins(0),
        fcollidingSys(""),
+       ftriggertype(""),
        fpvzcut(0),
        fTrackEtaCut(0),
        fFilterBit(128),
@@ -126,6 +125,7 @@ AliLeadingV0Correlation::AliLeadingV0Correlation(const char *name)
        fNzVtxBins(0), 
        fNCentBins(0),
        fcollidingSys(""),
+       ftriggertype(""),
        fpvzcut(0),
     fTrackEtaCut(0),
     fFilterBit(128),
@@ -199,10 +199,10 @@ void AliLeadingV0Correlation::UserCreateOutputObjects()
        
        //---------------------------------------------- Events histograms -----------------------------------------------------//
        //0-PVx,1-PVy,2-PVz,3-MULT,4-CENT
-       const Int_t ndimsEV = 5;
-       Int_t    binsEV[ndimsEV] = { 100,  100, 200, 1000,100};
-       Double_t xminEV[ndimsEV] = {-0.5, -0.5,-20 ,    0,  0};
-       Double_t xmaxEV[ndimsEV] = { 0.5,  0.5, 20 ,20000,100};
+       const Int_t ndimsEV = 3;
+       Int_t    binsEV[ndimsEV] = { 200, 1000,100};
+       Double_t xminEV[ndimsEV] = {-20 ,    0,  0};
+       Double_t xmaxEV[ndimsEV] = { 20 ,20000,100};
        
        fHistEventViceGen= new THnSparseD("fHistEventViceGen", "fHistEventViceGen", ndimsEV, binsEV, xminEV, xmaxEV);
        fOutputList->Add(fHistEventViceGen);
@@ -229,33 +229,37 @@ void AliLeadingV0Correlation::UserCreateOutputObjects()
        fHistMCGenK0   = new THnSparseD("fHistMCGenK0"  , "fHistMCGenK0"  , ndimsGenMC, binsGenMCK0, xminGenMCK0, xmaxGenMCK0);
        fOutputList->Add(fHistMCGenK0);
        
-       const Int_t ndims=10;    //MK0  mLA  MALA DCAP DCAN  RV0 DVD   CPA  PT   cent
-       Int_t    bins[ndims] = {  200,  140,  140, 500, 500, 110,110,   200,240  ,100};
-       Double_t xmin[ndims] = {  0.4, 1.06, 1.06,   0,   0,   0,  0, 0.997,  0  ,  0};
-       Double_t xmax[ndims] = {  0.6,  1.2,  1.2,  10,  10, 110,1.1, 1.007, 12  ,100}; 
+       const Int_t ndims=3;    //MK0  mLA  MALA PT   cent
+       Int_t    binsK0[ndims] = {  200, 240  ,100};
+       Double_t xminK0[ndims] = {  0.4,  0   ,  0};
+       Double_t xmaxK0[ndims] = {  0.6, 12   ,100};
+        
+       Int_t    binsLA[ndims] = {  140, 240  ,100};
+       Double_t xminLA[ndims] = { 1.06,   0  ,  0};
+       Double_t xmaxLA[ndims] = {  1.2,  12  ,100};
        
        
-       fHistReconstK0= new THnSparseD("fHistReconstK0"  , "fHistReconstK0",  ndims, bins, xmin, xmax);
+       fHistReconstK0= new THnSparseD("fHistReconstK0"  , "fHistReconstK0",  ndims, binsK0, xminK0, xmaxK0);
        fHistReconstK0->Sumw2();
        fOutputList->Add(fHistReconstK0);
        
-       fHistReconstLA= new THnSparseD("fHistReconstLA"  , "fHistReconstLA",  ndims, bins, xmin, xmax);
+       fHistReconstLA= new THnSparseD("fHistReconstLA"  , "fHistReconstLA",  ndims, binsLA, xminLA, xmaxLA);
        fHistReconstLA->Sumw2();
        fOutputList->Add(fHistReconstLA);
        
-       fHistReconstALA= new THnSparseD("fHistReconstALA", "fHistReconstALA", ndims, bins, xmin, xmax);
+       fHistReconstALA= new THnSparseD("fHistReconstALA", "fHistReconstALA", ndims, binsLA, xminLA, xmaxLA);
        fHistReconstALA->Sumw2();
        fOutputList->Add(fHistReconstALA);
        
-       fHistMCAssoK0= new THnSparseD("fHistMCAssoK0"   , "fHistMCAssoK0"   , ndims, bins, xmin, xmax);
+       fHistMCAssoK0= new THnSparseD("fHistMCAssoK0"   , "fHistMCAssoK0"   , ndims, binsK0, xminK0, xmaxK0);
        fHistMCAssoK0->Sumw2();
        fOutputList->Add(fHistMCAssoK0);
        
-       fHistMCAssoLA= new THnSparseD("fHistMCAssoLA"   , "fHistMCAssoLA"   , ndims, bins, xmin, xmax);
+       fHistMCAssoLA= new THnSparseD("fHistMCAssoLA"   , "fHistMCAssoLA"   , ndims, binsLA, xminLA, xmaxLA);
        fHistMCAssoLA->Sumw2();
        fOutputList->Add(fHistMCAssoLA);
        
-       fHistMCAssoALA= new THnSparseD("fHistMCAssoALA" , "fHistMCAssoALA" ,  ndims, bins, xmin, xmax);
+       fHistMCAssoALA= new THnSparseD("fHistMCAssoALA" , "fHistMCAssoALA" ,  ndims, binsLA, xminLA, xmaxLA);
        fHistMCAssoALA->Sumw2();
        fOutputList->Add(fHistMCAssoALA);
        
@@ -342,7 +346,7 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
        Double_t CentPecentMC = -1;
        Double_t CentPecentAfterPhySel   = -1;
        
-       if(fcollidingSys=="PbPb2010" || fcollidingSys=="PbPb2011"){   
+       if(fcollidingSys=="PbPb"){   
        AliCentrality *centralityObjMC = fAODEvent->GetHeader()->GetCentralityP();
                CentPecentMC  = centralityObjMC->GetCentralityPercentileUnchecked("V0M");
                if ((CentPecentMC < 0.)||(CentPecentMC > 90)) return;
@@ -353,7 +357,7 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
        Double_t poolmax    = CentBins[fNCentBins];
        
        
-       Double_t  dimEventviceMC[5];
+       Double_t  dimEventviceMC[3];
        if(fAnalysisMC)    //Efficency denomenator comes before the physics selection
        {
                AliAODMCHeader *aodMCheader = (AliAODMCHeader*)fAODEvent->FindListObject(AliAODMCHeader::StdBranchName());
@@ -361,27 +365,28 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
                
                if (TMath::Abs(mcZv) >= fpvzcut) return;
                
-               dimEventviceMC[0]=aodMCheader->GetVtxX();
-               dimEventviceMC[1]=aodMCheader->GetVtxY();
-               dimEventviceMC[2]=aodMCheader->GetVtxZ();
+               dimEventviceMC[0]=aodMCheader->GetVtxZ();
                
                TClonesArray *mcArray = (TClonesArray*)fAODEvent->FindListObject(AliAODMCParticle::StdBranchName());
                if(!mcArray)return;
                
-               TObjArray *selectedTracksLeadingMC=fAnalyseUE->FindLeadingObjects(mcArray);
+               TObjArray* selectedTracksLeadingMC = 0;
+               if(ftriggertype=="Leading")selectedTracksLeadingMC=fAnalyseUE->FindLeadingObjects(mcArray);
+               if(ftriggertype=="Alltrigs")selectedTracksLeadingMC = fAnalyseUE->GetAcceptedParticles(mcArray, 0, kTRUE, -1, kTRUE);
                if(!selectedTracksLeadingMC) return;
                selectedTracksLeadingMC->SetOwner(kTRUE);
                
+               
                TObjArray * selectedV0sMC =new TObjArray;
                selectedV0sMC->SetOwner(kTRUE);
                
                Int_t nMCTracks = mcArray->GetEntriesFast();
                
-               if(fcollidingSys=="PbPb2010" || fcollidingSys=="PbPb2011") multiplicityMC=CentPecentMC;
+               if(fcollidingSys=="PbPb") multiplicityMC=CentPecentMC;
                if(fcollidingSys=="PP")   multiplicityMC=nMCTracks;
                
-               dimEventviceMC[3]=nMCTracks;
-               dimEventviceMC[4]=CentPecentMC;
+               dimEventviceMC[1]=nMCTracks;
+               dimEventviceMC[2]=CentPecentMC;
                fHistEventViceGen->Fill(dimEventviceMC);
                
                for (Int_t iMC = 0; iMC<nMCTracks; iMC++)
@@ -448,12 +453,12 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
        if(fRemovePileUP)
        if(fAODEvent->IsPileupFromSPD()) return;
        
-       Double_t  dimEventviceReal[5];
+       Double_t  dimEventviceReal[3];
        Double_t  lPrimaryVtxPosition[3];
        Double_t  lV0Position[3];
        
        
-       if(fcollidingSys=="PbPb2010" || fcollidingSys=="PbPb2011"){  //
+       if(fcollidingSys=="PbPb"){  //
        AliCentrality *centralityObj = fAODEvent->GetHeader()->GetCentralityP();
                CentPecentAfterPhySel  = centralityObj->GetCentralityPercentileUnchecked("V0M");
                if ((CentPecentAfterPhySel < 0.)||(CentPecentAfterPhySel > 90)) return;
@@ -469,20 +474,18 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
        
        if ((TMath::Abs(lPVz)) >= fpvzcut) return ;
        if (TMath::Abs(lPVx)<10e-5 && TMath::Abs(lPVy)<10e-5 && TMath::Abs(lPVz)<10e-5) return;
-       
-       dimEventviceReal[0]=lPVx;
-    dimEventviceReal[1]=lPVy;
-    dimEventviceReal[2]=lPVz;
+
+    dimEventviceReal[0]=lPVz;
        
        multiplicity = fAODEvent->GetNTracks();
        
-       dimEventviceReal[3]=multiplicity;
-       dimEventviceReal[4]=CentPecentAfterPhySel;
+       dimEventviceReal[1]=multiplicity;
+       dimEventviceReal[2]=CentPecentAfterPhySel;
        
        fHistEventViceReconst->Fill(dimEventviceReal);
        
        if(fcollidingSys=="PP")MultipOrCent=multiplicity;
-       if(fcollidingSys=="PbPb2010" || fcollidingSys=="PbPb2011")MultipOrCent=CentPecentAfterPhySel;
+       if(fcollidingSys=="PbPb")MultipOrCent=CentPecentAfterPhySel;
 
        //---------------------------------------------------------------------------------------------
        
@@ -497,7 +500,9 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
        Double_t lPzV0s       = 0; Double_t lAlphaV0 = 0, lPtArmV0           = 0;
        Double_t lPV0s        = 0;
        
-       TObjArray *selectedTracksLeading=fAnalyseUE->FindLeadingObjects(fAODEvent);
+       TObjArray *selectedTracksLeading=0;
+       if(ftriggertype=="Leading")selectedTracksLeading=fAnalyseUE->FindLeadingObjects(fAODEvent);
+       if(ftriggertype=="Alltrigs")selectedTracksLeading = fAnalyseUE->GetAcceptedParticles(fAODEvent, 0, kTRUE, -1, kTRUE);
        if(!selectedTracksLeading) return;
        selectedTracksLeading->SetOwner(kTRUE);
        
@@ -574,7 +579,7 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
         if (pPid)
         {
             Double_t pdMom = pPid->GetTPCmomentum();
-            if (pdMom<1.0 && (fcollidingSys=="PbPb2010"||fcollidingSys=="PbPb2011"))
+            if (pdMom<1.0 && (fcollidingSys=="PbPb"))
             {
                 nSigmaPosPion   = fPIDResponse->NumberOfSigmasTPC(myTrackPos, AliPID::kPion);
                 nSigmaPosProton = fPIDResponse->NumberOfSigmasTPC(myTrackPos, AliPID::kProton);
@@ -590,7 +595,7 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
         if (nPid)
         {
             Double_t ndMom = nPid->GetTPCmomentum();
-            if (ndMom<1.0 && (fcollidingSys=="PbPb2010"||fcollidingSys=="PbPb2011"))
+            if (ndMom<1.0 && (fcollidingSys=="PbPb"))
             {
                 nSigmaNegPion   = fPIDResponse->NumberOfSigmasTPC(myTrackNeg, AliPID::kPion);
                 nSigmaNegProton = fPIDResponse->NumberOfSigmasTPC(myTrackNeg, AliPID::kProton);
@@ -628,14 +633,15 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
                Bool_t laRapcut = (TMath::Abs(lRapLambda)     < fRapidityCut);
                Bool_t alaRapcut= (TMath::Abs(lRapAntiLambda) < fRapidityCut);
                
+               if(fcollidingSys=="PbPb")if(lV0Radius>=100) continue;
                
                Bool_t k0cutset = IsAcseptedK0(lV0Radius,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lDcaV0Daughters,lV0cosPointAngle,lInvMassLambda,lInvMassAntiLambda);
                Bool_t lacutset = IsAcseptedLA(lV0Radius,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lDcaV0Daughters,lV0cosPointAngle,lInvMassK0);
                Bool_t alacutset= IsAcseptedLA(lV0Radius,lDcaNegToPrimVertex,lDcaPosToPrimVertex,lDcaV0Daughters,lV0cosPointAngle,lInvMassK0);
                
-               Double_t spK0[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
-               Double_t spLa[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
-               Double_t spAl[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
+               Double_t spK0[3] = {lInvMassK0, lPtV0s,MultipOrCent};
+               Double_t spLa[3] = {lInvMassLambda,lPtV0s,MultipOrCent};
+               Double_t spAl[3] = {lInvMassAntiLambda,lPtV0s,MultipOrCent};
        
                switch (fCase) {
                        case 1:
@@ -720,9 +726,9 @@ void AliLeadingV0Correlation::UserExec(Option_t *)
                        Int_t MotherPdg  = mcPosMother->GetPdgCode();
                        Bool_t IsPrime   = mcPosMother->IsPhysicalPrimary();
                        
-                       Double_t rcK0[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
-                       Double_t rcLa[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
-                       Double_t rcAl[10] = {lInvMassK0, lInvMassLambda,lInvMassAntiLambda,lDcaPosToPrimVertex,lDcaNegToPrimVertex,lV0Radius,lDcaV0Daughters,lV0cosPointAngle,lPtV0s,MultipOrCent};
+                       Double_t rcK0[3] = {lInvMassK0, lPtV0s,MultipOrCent};
+                       Double_t rcLa[3] = {lInvMassLambda,lPtV0s,MultipOrCent};
+                       Double_t rcAl[3] = {lInvMassAntiLambda,lPtV0s,MultipOrCent};
                        
                        switch (fCase) {
                                case 1:
@@ -962,6 +968,8 @@ void AliLeadingV0Correlation::FillCorrelationSibling(Double_t MultipOrCent,
                
                if(triggerPt<6.0||triggerPt>12.0)continue;
                counterSibMCA++;
+               if(ftriggertype=="Leading")
+               {
                if(counterSibMCA==triggerArray->GetEntriesFast()){
                        
                        binsTrigSib[0]=triggerPt;
@@ -971,6 +979,36 @@ void AliLeadingV0Correlation::FillCorrelationSibling(Double_t MultipOrCent,
                        
                        if(triggerHist)triggerHist->Fill(binsTrigSib);
                        
+                       for (Int_t j=0; j<selectedV0Array->GetEntriesFast(); j++){
+                               
+                               V0Correlationparticle* associate = (V0Correlationparticle*) selectedV0Array->At(j);
+                               if(!associate)continue;
+                               
+                               binsv0CORR[0]= associate->Pt();
+                               binsv0CORR[1]= associate->Phi();
+                               binsv0CORR[2]= associate->Eta();
+                               
+                               if(binsv0CORR[0]>triggerPt) continue;
+                               
+                               binsv0CORR[3]=RangePhi(triggerPhi-binsv0CORR[1]);
+                               binsv0CORR[4]=triggerEta-binsv0CORR[2];
+                               binsv0CORR[5]= associate->WhichCandidate();
+                               binsv0CORR[6]= MultipOrCent;
+                               
+                               associateHist->Fill(binsv0CORR);
+                       }
+                 }
+               }
+               
+               if(ftriggertype=="Alltrigs")
+               {
+                       binsTrigSib[0]=triggerPt;
+                       binsTrigSib[1]=triggerPhi;
+                       binsTrigSib[2]=triggerEta;
+                       binsTrigSib[3]=MultipOrCent;
+                       
+                       if(triggerHist)triggerHist->Fill(binsTrigSib);
+                       
                        for (Int_t j=0; j<selectedV0Array->GetEntriesFast(); j++){
                                
                                V0Correlationparticle* associate = (V0Correlationparticle*) selectedV0Array->At(j);
@@ -1005,7 +1043,7 @@ void AliLeadingV0Correlation::FillCorrelationMixing(Double_t MultipOrCentMix,
        if(TMath::Abs(pvxMix)>=fpvzcut || MultipOrCentMix>poolmax || MultipOrCentMix < poolmin)
        {
                if(fcollidingSys=="PP")AliInfo(Form("pp Event with Zvertex = %.2f cm and multiplicity = %.0f out of pool bounds, SKIPPING",pvxMix,MultipOrCentMix));
-               if(fcollidingSys=="PbPb2010" || fcollidingSys=="PbPb2011") AliInfo(Form("PbPb Event with Zvertex = %.2f cm and centrality = %.1f  out of pool bounds, SKIPPING",pvxMix,MultipOrCentMix));
+               if(fcollidingSys=="PbPb") AliInfo(Form("PbPb Event with Zvertex = %.2f cm and centrality = %.1f  out of pool bounds, SKIPPING",pvxMix,MultipOrCentMix));
                return;
        }
        
@@ -1036,7 +1074,8 @@ void AliLeadingV0Correlation::FillCorrelationMixing(Double_t MultipOrCentMix,
                                
                                if(trigPt<6.0||trigPt>12.0)continue;
                                counterMix++;
-                               
+                               if(ftriggertype=="Leading")
+                               {
                                if(counterMix==triggerArray->GetEntriesFast()){
                                        
                                        binsTrigMix[0]=trigPt;
@@ -1046,6 +1085,36 @@ void AliLeadingV0Correlation::FillCorrelationMixing(Double_t MultipOrCentMix,
                                        
                                        if(triggerHist)triggerHist->Fill(binsTrigMix);
                                        
+                                       for (Int_t j=0; j<mixEvents->GetEntriesFast(); j++){
+                                               
+                                               V0Correlationparticle* associate = (V0Correlationparticle*) mixEvents->At(j);
+                                               if(!associate)continue;
+                                               
+                                               binsv0CORRMix[0]= associate->Pt();
+                                               binsv0CORRMix[1]= associate->Phi();
+                                               binsv0CORRMix[2]= associate->Eta();
+                                               
+                                               if(binsv0CORRMix[0]>trigPt) continue;
+                                               
+                                               binsv0CORRMix[3]=RangePhi(trigPhi-binsv0CORRMix[1]);
+                                               binsv0CORRMix[4]=trigEta-binsv0CORRMix[2];
+                                               binsv0CORRMix[5]=associate->WhichCandidate();
+                                               binsv0CORRMix[6]=MultipOrCentMix;
+                                               
+                                               associateHist->Fill(binsv0CORRMix);
+                                       }
+                                 }
+                               }
+                               
+                               if(ftriggertype=="Alltrigs")
+                               {
+                                       binsTrigMix[0]=trigPt;
+                                       binsTrigMix[1]=trigPhi;
+                                       binsTrigMix[2]=trigEta;
+                                       binsTrigMix[3]=MultipOrCentMix;
+                                       
+                                       if(triggerHist)triggerHist->Fill(binsTrigMix);
+                                       
                                        for (Int_t j=0; j<mixEvents->GetEntriesFast(); j++){
                                                
                                                V0Correlationparticle* associate = (V0Correlationparticle*) mixEvents->At(j);
@@ -1086,4 +1155,4 @@ void AliLeadingV0Correlation::FillCorrelationMixing(Double_t MultipOrCentMix,
                                                                
                                        
                                                                                                
-                                                                                               
+                                                                                               
\ No newline at end of file

diff --git a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.h b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.h
index 7c324f45e37ff1b4e53331907b190d0f9b6cb029..a2c79ae995ac37fe7bc79bcc822495cfe9120c71 100644 (file)
--- a/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.h
+++ b/PWGCF/Correlations/DPhi/AliLeadingV0Correlation.h
@@ -54,6 +54,7 @@ public:
        }
        
        void SetCollidingSystem(TString system){fcollidingSys = system;}
+       void SetTriggerType(TString triggertype){ftriggertype = triggertype;}
        void SetPrimeryVertexCut(Double_t pvzcut){fpvzcut = pvzcut;}
        void SetEatCut(Double_t  TrackEtaCut){fTrackEtaCut = TrackEtaCut;}
        void SetFilterBit(UInt_t  filterBit){fFilterBit = filterBit;}
@@ -125,7 +126,8 @@ private:
        Int_t fNCentBins;                            // number of centrality bins
        Double_t fCentBins[100];                     // [fNCentBinsDim]
        
-       TString         fcollidingSys;               // "PP" or "PbPb2010,2011"
+       TString         fcollidingSys;               // "PP" or "PbPb"
+       TString         ftriggertype;                // "Leading" or "Alltrigs"
        Double_t        fpvzcut;                     // PVz cut of event
     Double_t           fTrackEtaCut;                // Eta cut on particles
     UInt_t             fFilterBit;                  // Select tracks from an specific track cut (default 0xFF all track selected)

diff --git a/PWGCF/Correlations/macros/AddTaskLV0Correlation.C b/PWGCF/Correlations/macros/AddTaskLV0Correlation.C
index adb03b66cd41113fc46ff35b79655edfea6ee0ac..5500f459bc10d62a143f20da2b05a01d77dfeffa 100644 (file)
--- a/PWGCF/Correlations/macros/AddTaskLV0Correlation.C
+++ b/PWGCF/Correlations/macros/AddTaskLV0Correlation.C
@@ -1,5 +1,6 @@
 AliLeadingV0Correlation* AddTaskLV0Correlation(TString  fListName                = "LV0Correlation",
                                                                                           TString      fCollisiontype           = "PP",
+                                                                                          TString  ftriggertype             = "Alltrigs",
                                                                                           Bool_t       fAnalysisMC              = 0,
                                                                                           Int_t    fCase                    = 2,
                                                                                           Bool_t   fRemovePileUP            = 0,
@@ -30,7 +31,7 @@ AliLeadingV0Correlation* AddTaskLV0Correlation(TString  fListName
     if (!mgr) {Error("AddTaskLV0Correlation.C", "No Analysis Manager ");return 0;}
        
        //PVz Binning for pool PP or PbPb       
-       Double_t pvzbinlimits[] = {-10,-9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10};
+       Double_t pvzbinlimits[] = {-12,-11,-10,-9,-8,-7,-6,-5,-4,-3,-2,-1,0,1,2,3,4,5,6,7,8,9,10,11,12};
        Int_t pvzbinnumb = sizeof(pvzbinlimits)/sizeof(Double_t) - 1;
        
        //Mult Binning for pool pp 
@@ -56,6 +57,7 @@ AliLeadingV0Correlation* AddTaskLV0Correlation(TString  fListName
        AliLeadingV0Correlation *myTask = new AliLeadingV0Correlation(fListName.Data());
        
        myTask->SetCollidingSystem(fCollisiontype);
+       myTask->SetTriggerType(ftriggertype);
        myTask->SetMCAnalysis(fAnalysisMC);
        myTask->SetCase(fCase);
        myTask->SetRemovePileUp(fRemovePileUP);
@@ -65,7 +67,7 @@ AliLeadingV0Correlation* AddTaskLV0Correlation(TString  fListName
        myTask->SetMinEventsToMix(fMinEventsToMix);
        myTask->SetPoolPVzBinLimits(pvzbinnumb,pvzbinlimits);
        if(fCollisiontype=="PP")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPP,cent_mult_binlimitsPP);
-       if(fCollisiontype=="PbPb2010"||fCollisiontype=="PbPb2011")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPbPb,cent_mult_binlimitsPbPb);
+       if(fCollisiontype=="PbPb"||fCollisiontype=="PbPb")myTask->SetPoolCentBinLimits(cent_mult_bin_numbPbPb,cent_mult_binlimitsPbPb);
        myTask->SetPrimeryVertexCut(fPVzCut);
        myTask->SetEatCut(fEtaCut);
        myTask->SetFilterBit(fFilterBit);