From 71b568564c60f84e9ec7091b0bf1077677a25b53 Mon Sep 17 00:00:00 2001 From: bdoenigu Date: Mon, 30 Jun 2014 12:55:26 +0200 Subject: [PATCH] Commit for Simone - Added phi-eta histos --- .../LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx | 111 ++++++++++-------- 1 file changed, 63 insertions(+), 48 deletions(-) diff --git a/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx b/PWGLF/STRANGENESS/LambdaK0PbPb/AliAnalysisTaskV0ForRAA.cxx index 488f1e7fd69..6c88ef49489 100644 --- 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); -- 2.43.0