Commit for Simone - Added phi-eta histos

diff --git a/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx b/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx
index 488f1e7fd69686bc5fc1c601a9ab595cceed6446..6c88ef49489833d1034ba0888990bef8cee8a1c2 100644 (file)
--- a/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx
+++ b/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx
@@ -1111,28 +1111,28 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   // ------------------- add always ---------------------------//
   //THnF
   Double_t piForAx = 2.0*TMath::Pi();
-  Int_t binsTHnV0K0s[4] = {150,100,50,18};
-  Int_t binsTHnV0L[4] = {100,100,50,18};
-
-  Int_t binsTHnV0DauEtaK0s[4] = {150,100,40,40};
-  Int_t binsTHnV0DauEtaL[4] = {100,100,40,40};
-
-  Int_t binsTHnV0DauPhiK0s[5] = {150, 18,18, 7,7};
-  Int_t binsTHnV0DauPhiL[5] = {100, 18,18, 7,7};
-
-  Double_t minK0s[4] = {0.35,0.0,-1.0,0.0};
-  Double_t maxK0s[4] = {0.65,50.0,1.0,piForAx};
-  Double_t minK0sDauEta[4] = {0.35, 0.0,-0.8,-0.8};
-  Double_t maxK0sDauEta[4] = {0.65,50.0, 0.8, 0.8};
-  Double_t minK0sDauPhi[5] = {0.35,0.0,0.0,-0.5,-0.5};
-  Double_t maxK0sDauPhi[5] = {0.65,piForAx,piForAx,6.5,6.5};
-
-  Double_t minL[4] = {1.07, 0.0,-1.0,0.0};
-  Double_t maxL[4] = {1.17,50.0, 1.0,piForAx};
-  Double_t minLDauEta[4] = {1.07, 0.0,-0.8,-0.8};
-  Double_t maxLDauEta[4] = {1.17,50.0, 0.8, 0.8};
-  Double_t minLDauPhi[5] = {1.07,0.0,0.0,-0.5,-0.5};
-  Double_t maxLDauPhi[5] = {1.17,piForAx,piForAx,6.5, 6.5};
+  Int_t binsTHnV0K0s[4] = {150,100,80,18};
+  Int_t binsTHnV0L[4] = {100,100,80,18};
+
+  Int_t binsTHnV0DauEtaK0s[4] = {150,100,40,18};
+  Int_t binsTHnV0DauEtaL[4] = {100,100,40,18};
+
+  Int_t binsTHnV0DauPhiK0s[4] = {150, 100,40,18};
+  Int_t binsTHnV0DauPhiL[4] = {100, 100,40,18};
+
+  Double_t minK0s[4] = {0.35,0.0,0.0,0.0};
+  Double_t maxK0s[4] = {0.65,50.0,40.0,piForAx};
+  Double_t minK0sDauEta[4] = {0.35, 0.0,-0.8,0.0};
+  Double_t maxK0sDauEta[4] = {0.65,50.0, 0.8,piForAx};
+  Double_t minK0sDauPhi[4] = {0.35,0.0,-0.8,0.0};
+  Double_t maxK0sDauPhi[4] = {0.65,50.0,0.8,piForAx};
+
+  Double_t minL[4] = {1.07, 0.0,0.0,0.0};
+  Double_t maxL[4] = {1.17,50.0, 40.0,piForAx};
+  Double_t minLDauEta[4] = {1.07, 0.0,-0.8,0.0};
+  Double_t maxLDauEta[4] = {1.17,50.0, 0.8, piForAx};
+  Double_t minLDauPhi[4] = {1.07,0.0,-0.8,0.0};
+  Double_t maxLDauPhi[4] = {1.17,50.0,0.8,piForAx};
   
   
   char histTitK0s[255];
@@ -1159,27 +1159,27 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   snprintf(histTitALDauPhi,255,"fTHnFALDauPhi");
  
   char axTitK0s[255];
-  snprintf(axTitK0s,255,"K^{0}_{s};m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta(V0);#phi(V0)");
+  snprintf(axTitK0s,255,"K^{0}_{s};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
   char axTitL[255];
-  snprintf(axTitL,255,"#Lambda;m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta(V0);#phi(V0)");
+  snprintf(axTitL,255,"#Lambda;m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
   char axTitAL[255];
-  snprintf(axTitAL,255,"#bar{#Lambda};m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta(V0);#phi(V0)");
+  snprintf(axTitAL,255,"#bar{#Lambda};m_{inv} (GeV/c^{2});p_{T} (GeV/c);c#tau(V0);#phi(V0)");
 
 
   char axTitK0sDauEta[255];
-  snprintf(axTitK0sDauEta,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});p_{T} (Gev/c);#eta_{pos};#eta_{neg}");
+  snprintf(axTitK0sDauEta,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});p_{T} (Gev/c);#eta_{pos};#phi_{pos}");
   char axTitLDauEta[255];
-  snprintf(axTitLDauEta,255,"#Lambda daughter;m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta_{pos};#eta_{neg}");
+  snprintf(axTitLDauEta,255,"#Lambda daughter;m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta_{pos};#phi_{pos}");
   char axTitALDauEta[255];
-  snprintf(axTitALDauEta,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta_{pos};#eta_{neg}");
+  snprintf(axTitALDauEta,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});p_{T} (GeV/c);#eta_{pos};#phi_{pos}");
 
 
   char axTitK0sDauPhi[255];
-  snprintf(axTitK0sDauPhi,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});#phi_{pos};#phi_{neg};ITS hits (pos);ITS hits (neg)");
+  snprintf(axTitK0sDauPhi,255,"K^{0}_{s} daughter;m_{inv} (GeV/c^{2});p_{T}(GeV/c);#eta_{neg};#phi_{neg}");
   char axTitLDauPhi[255];
-  snprintf(axTitLDauPhi,255,"#Lambda daughter;m_{inv} (GeV/c^{2});#phi_{pos};#phi_{neg};ITS hits (pos);ITS hits (neg)");
+  snprintf(axTitLDauPhi,255,"#Lambda daughter;m_{inv} (GeV/c^{2});p_{T}(GeV/c);#eta_{neg};#phi_{neg}");
   char axTitALDauPhi[255];
-  snprintf(axTitALDauPhi,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});#phi_{pos};#phi_{neg};ITS hits (pos);ITS hits (neg)");
+  snprintf(axTitALDauPhi,255,"#bar{#Lambda} daughter;m_{inv} (GeV/c^{2});p_{T}(GeV/c);#eta_{neg};#phi_{neg}");
 
 
   fTHnFK0s = new       THnF(histTitK0s,axTitK0s,4,binsTHnV0K0s,minK0s,maxK0s);
@@ -1197,11 +1197,11 @@ void AliAnalysisTaskV0ForRAA::UserCreateOutputObjects(){
   fTHnFALDauEta  = new         THnF(histTitALDauEta ,axTitALDauEta ,4,binsTHnV0DauEtaL,minLDauEta  ,maxLDauEta);
   // fTHnFALDauEta->Sumw2();
 
-  fTHnFK0sDauPhi = new         THnF(histTitK0sDauPhi,axTitK0sDauPhi,5,binsTHnV0DauPhiK0s,minK0sDauPhi,maxK0sDauPhi);
+  fTHnFK0sDauPhi = new         THnF(histTitK0sDauPhi,axTitK0sDauPhi,4,binsTHnV0DauPhiK0s,minK0sDauPhi,maxK0sDauPhi);
   // fTHnFK0sDauPhi->Sumw2();
-  fTHnFLDauPhi   = new         THnF(histTitLDauPhi  ,axTitLDauPhi  ,5,binsTHnV0DauPhiL,minLDauPhi  ,maxLDauPhi);
+  fTHnFLDauPhi   = new         THnF(histTitLDauPhi  ,axTitLDauPhi  ,4,binsTHnV0DauPhiL,minLDauPhi  ,maxLDauPhi);
   // fTHnFLDauPhi->Sumw2();
-  fTHnFALDauPhi  = new         THnF(histTitALDauPhi ,axTitALDauPhi ,5,binsTHnV0DauPhiL,minLDauPhi  ,maxLDauPhi); 
+  fTHnFALDauPhi  = new         THnF(histTitALDauPhi ,axTitALDauPhi ,4,binsTHnV0DauPhiL,minLDauPhi  ,maxLDauPhi); 
   //fTHnFALDauPhi->Sumw2();
 
  
@@ -2658,11 +2658,11 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
 
     
       Double_t phiMC =   p0->Phi(); 
-      Double_t etaMC =   p0->Eta(); 
+      //Double_t etaMC =   p0->Eta(); 
 
-      Double_t valTHnMC[4] = {massV0MC,ptV0MC,etaMC,phiMC};
-      Double_t valTHnMCDauEta[4] = {massV0MC,ptV0MC,etaMCPos,etaMCNeg};
-      Double_t valTHnMCDauPhi[5] = {massV0MC,phiMCPos,phiMCNeg,0.0,0.0};
+      Double_t valTHnMC[4]  = {massV0MC,ptV0MC,0.0,phiMC}; 
+      Double_t valTHnMCDauEta[4] = {massV0MC,ptV0MC,etaMCPos,phiMCPos};
+      Double_t valTHnMCDauPhi[4] = {massV0MC,ptV0MC,etaMCNeg,phiMCNeg};
 
       //-- Fill Particle histos --//
       if (pdgCode==310){//K0s
@@ -2690,7 +2690,8 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
        fHistV0RadiusZVSPt[isSecd]->Fill(ptV0MC,xyzMC[2]);
        // fHistPiPiPhiPosVsPtPosVsMass->Fill(massV0MC,ctTK0s,ptV0MC);//,ctK0s);//phiPosMC);//xxx
        fHistPiPiK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
-       fHistPiPiK0sVsALambdaMass->Fill(calcALambdamass,calcK0smass); 
+       fHistPiPiK0sVsALambdaMass->Fill(calcALambdamass,calcK0smass);
+       valTHnMC[2] = ctTK0s; 
        fTHnFK0s->Fill(valTHnMC);
        fTHnFK0sDauEta->Fill(valTHnMCDauEta);
        fTHnFK0sDauPhi->Fill(valTHnMCDauPhi);
@@ -2722,7 +2723,9 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
        fHistV0RadiusZVSPt[isSecd]->Fill(ptV0MC,xyzMC[2]);
        
        fHistPiPK0sVsLambdaMass->Fill(calcLambdamass,calcK0smass);
-       fTHnFL->Fill(valTHnMC);
+
+              valTHnMC[2] = ctTL;
+fTHnFL->Fill(valTHnMC);
        fTHnFLDauEta->Fill(valTHnMCDauEta);
        fTHnFLDauPhi->Fill(valTHnMCDauPhi);
       }
@@ -2755,6 +2758,7 @@ void AliAnalysisTaskV0ForRAA::V0MCTruthLoop(){
        // if(isSecd <1) fHistPiPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
        //else fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx       
        //fHistPiAPPhiPosVsPtPosVsMass->Fill(massV0MC,ctTL,ptV0MC);//,ctK0s);//phiPosMC);//xxx
+       valTHnMC[2] = ctTAL;
        fTHnFAL->Fill(valTHnMC);
        fTHnFALDauEta->Fill(valTHnMCDauEta);
        fTHnFALDauPhi->Fill(valTHnMCDauPhi);
@@ -3473,7 +3477,7 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
     Double_t px = v0K0.Px();
     Double_t py = v0K0.Py();
     Double_t phi  = TMath::Pi()+TMath::ATan2(-py, -px);
-    Double_t eta =  v0K0.Eta();
+    //    Double_t eta =  v0K0.Eta();
 
     /*     
     //introduce more histo
@@ -3668,9 +3672,9 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            fHistPiPiPtVSY->Fill(rapK0s,ptK0s);
            fHistPiPiDecayLengthVsMass->Fill(massK0s,dim2V0Radius);//decayLength);
            // fHistPiPiPhiPosVsPtPosVsMass->Fill(massK0s,ctTK0,ptV0MC);//,ctK0);//posDaughterPhi);//xxx
-           Double_t valTHnK0s[4]= {massK0s,ptV0MC,eta,phi};
-           Double_t valTHnK0sDauEta[4]= {massK0s,ptV0MC,posDaughterEta,negDaughterEta};
-           Double_t valTHnK0sDauPhi[5]= {massK0s,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
+           Double_t valTHnK0s[4]= {massK0s,ptV0MC,ctTK0,phi};
+           Double_t valTHnK0sDauEta[4]= {massK0s,ptV0MC,posDaughterEta,posDaughterPhi};
+           Double_t valTHnK0sDauPhi[4]= {massK0s,ptV0MC,negDaughterEta,negDaughterPhi};
            if(massK0s >=0.35 && massK0s <= 0.65){
              fTHnFK0s->Fill(valTHnK0s);            
              fTHnFK0sDauEta->Fill(valTHnK0sDauEta);
@@ -3809,9 +3813,15 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            fHistPiPMassVSY[isSecd]->Fill(massLambda,rapL);
            fHistPiPPtVSY[isSecd]->Fill(rapL,ptLambda);
            //fHistPiPDecayLengthVsPt[isSecd]->Fill(ptLambda,ctL);
-           Double_t valTHnL[4]= {massLambda,ptV0MC,eta,phi};
-           Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,negDaughterEta};
-           Double_t valTHnLDauPhi[5]= {massLambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
+           //Double_t valTHnL[4]= {massLambda,ptV0MC,eta,phi};
+           //Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,negDaughterEta};
+           //      Double_t valTHnLDauPhi[5]= {massLambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
+
+           Double_t valTHnL[4]= {massLambda,ptV0MC,ctTL,phi};
+           Double_t valTHnLDauEta[4]= {massLambda,ptV0MC,posDaughterEta,posDaughterPhi};
+           Double_t valTHnLDauPhi[4]= {massLambda,ptV0MC,negDaughterEta,negDaughterPhi};
+
+
            if(massLambda >=1.07 && massLambda <= 1.17){
              fTHnFL->Fill(valTHnL);
              fTHnFLDauEta->Fill(valTHnLDauEta);
@@ -3963,10 +3973,15 @@ void AliAnalysisTaskV0ForRAA::V0RecoLoop(Int_t id0,Int_t id1,Int_t isSecd,Int_t
            // else {
            //  if(fMCTruthMode) fHistPiAPPhiPosVsPtPosVsMass->Fill(massLambda,ctTL,ptV0MC);
            // }              
-
+           /*
            Double_t valTHnAL[4]= {massALambda,ptV0MC,eta,phi};
            Double_t valTHnALDauEta[4]={massALambda,ptV0MC,posDaughterEta,negDaughterEta};
            Double_t valTHnALDauPhi[5]={massALambda,posDaughterPhi,negDaughterPhi,double(nclsITSPos),double(nclsITSNeg)};
+           */
+           Double_t valTHnAL[4]= {massALambda,ptV0MC,ctTAL,phi};
+           Double_t valTHnALDauEta[4]= {massALambda,ptV0MC,posDaughterEta,posDaughterPhi};
+           Double_t valTHnALDauPhi[4]= {massALambda,ptV0MC,negDaughterEta,negDaughterPhi};
+
            if(massALambda >=1.07 && massALambda <= 1.17){
              fTHnFALDauEta->Fill(valTHnALDauEta);
              fTHnFALDauPhi->Fill(valTHnALDauPhi);