From 271acb2e2b555ff30d55d107131a9fe8b3b674a9 Mon Sep 17 00:00:00 2001 From: xsanchez Date: Tue, 28 Oct 2014 14:39:58 +0100 Subject: [PATCH] Adding efficiency for the following selections: ctau, ficutial volumen and transverse momentum of daughters --- .../AliAnalysisTaskLambdaOverK0sJets.cxx | 372 +++++++++++++----- .../AliAnalysisTaskLambdaOverK0sJets.h | 27 ++ .../macros/AddTaskLambdaOverK0sJets.C | 8 +- 3 files changed, 312 insertions(+), 95 deletions(-) diff --git a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.cxx b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.cxx index a9675ebd298..6a5f7abf682 100644 --- a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.cxx +++ b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.cxx @@ -125,24 +125,32 @@ fAntiLambdaMass(0), fAntiLambdaMassEmbeded(0), fAntiLambdaMass2(0), fAntiLambdaM // K0s fK0sMCPtRapVtx[i] = 0; fK0sMCPtRapVtxEmbeded[i] = 0; + fK0sMCPtRapPtDaugPt[i] = 0x0; + fK0sMCPtRapPtDaugPtEmbeded[i] = 0x0; fK0sMCPtPhiEta[i] = 0; fK0sAssocPtPhiEta[i] = 0; // -- Natural particles - fK0sAssocPtMassArm[i] = 0; - fK0sAssocMassPtVtx[i] = 0; - fK0sAssocMassPtDCADaug[i] = 0; - fK0sAssocMassPtCPA[i] = 0; - fK0sAssocMassPtDCAPV[i] = 0; - fK0sAssocMassPtDaugNClsTPC[i] = 0; - fK0sAssocMassPtShTPCcls[i] = 0; + fK0sAssocPtMassArm[i] = 0x0; + fK0sAssocMassPtVtx[i] = 0x0; + fK0sAssocMassPtDCADaug[i] = 0x0; + fK0sAssocMassPtCPA[i] = 0x0; + fK0sAssocMassPtDCAPV[i] = 0x0; + fK0sAssocMassPtDaugNClsTPC[i] = 0x0; + fK0sAssocMassPtShTPCcls[i] = 0x0; + fK0sAssocMassPtDaugPt[i] = 0x0; + fK0sAssocMassPtCtau[i] = 0x0; + fK0sAssocMassPtFidVolume[i] = 0x0; // -- Embeded particles - fK0sAssocPtMassArmEmbeded[i] = 0; - fK0sAssocMassPtVtxEmbeded[i] = 0; - fK0sAssocMassPtDCADaug[i] = 0; - fK0sAssocMassPtCPAEmbeded[i] = 0; - fK0sAssocMassPtDCAPVEmbeded[i] = 0; - fK0sAssocMassPtDaugNClsTPCEmbeded[i] = 0; - fK0sAssocMassPtShTPCclsEmbeded[i] = 0; + fK0sAssocPtMassArmEmbeded[i] = 0x0; + fK0sAssocMassPtVtxEmbeded[i] = 0x0; + fK0sAssocMassPtDCADaug[i] = 0x0; + fK0sAssocMassPtCPAEmbeded[i] = 0x0; + fK0sAssocMassPtDCAPVEmbeded[i] = 0x0; + fK0sAssocMassPtDaugNClsTPCEmbeded[i] = 0x0; + fK0sAssocMassPtShTPCclsEmbeded[i] = 0x0; + fK0sAssocMassPtDaugPtEmbeded[i] = 0x0; + fK0sAssocMassPtCtauEmbeded[i] = 0x0; + fK0sAssocMassPtFidVolumeEmbeded[i] = 0x0; // -- Mass vs rapidity vs pt vs centrlaity fK0sMassPtRap[i] = 0; // -- Splitting checks @@ -156,26 +164,34 @@ fAntiLambdaMass(0), fAntiLambdaMassEmbeded(0), fAntiLambdaMass2(0), fAntiLambdaM // Lambda fLambdaMCPtRapVtx[i] = 0; fLambdaMCPtRapVtxEmbeded[i] = 0; + fLambdaMCPtRapPtDaugPt[i] = 0x0; + fLambdaMCPtRapPtDaugPtEmbeded[i] = 0x0; fLambdaMCPtPhiEta[i] = 0; fLambdaAssocPtPhiEta[i] = 0; // -- Natural particles - fLambdaAssocMassPtRap[i] = 0; - fLambdaAssocMassPtRap2[i] = 0; - fLambdaAssocMassPtVtx[i] = 0; - fLambdaAssocMassPtDCADaug[i] = 0; - fLambdaAssocMassPtCPA[i] = 0; - fLambdaAssocMassPtDCAPV[i] = 0; - fLambdaAssocMassPtDaugNClsTPC[i] = 0; - fLambdaAssocMassPtShTPCcls[i] = 0; + fLambdaAssocMassPtRap[i] = 0x0; + fLambdaAssocMassPtRap2[i] = 0x0; + fLambdaAssocMassPtVtx[i] = 0x0; + fLambdaAssocMassPtDCADaug[i] = 0x0; + fLambdaAssocMassPtCPA[i] = 0x0; + fLambdaAssocMassPtDCAPV[i] = 0x0; + fLambdaAssocMassPtDaugNClsTPC[i] = 0x0; + fLambdaAssocMassPtShTPCcls[i] = 0x0; + fLambdaAssocMassPtDaugPt[i] = 0x0; + fLambdaAssocMassPtCtau[i] = 0x0; + fLambdaAssocMassPtFidVolume[i] = 0x0; // -- Embeded particles - fLambdaAssocMassPtRapEmbeded[i] = 0; - fLambdaAssocMassPtRapEmbeded2[i] = 0; - fLambdaAssocMassPtVtxEmbeded[i] = 0; - fLambdaAssocMassPtDCADaug[i] = 0; - fLambdaAssocMassPtCPAEmbeded[i] = 0; - fLambdaAssocMassPtDCAPVEmbeded[i] = 0; - fLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0; - fLambdaAssocMassPtShTPCclsEmbeded[i] = 0; + fLambdaAssocMassPtRapEmbeded[i] = 0x0; + fLambdaAssocMassPtRapEmbeded2[i] = 0x0; + fLambdaAssocMassPtVtxEmbeded[i] = 0x0; + fLambdaAssocMassPtDCADaug[i] = 0x0; + fLambdaAssocMassPtCPAEmbeded[i] = 0x0; + fLambdaAssocMassPtDCAPVEmbeded[i] = 0x0; + fLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0x0; + fLambdaAssocMassPtShTPCclsEmbeded[i] = 0x0; + fLambdaAssocMassPtDaugPtEmbeded[i] = 0x0; + fLambdaAssocMassPtCtauEmbeded[i] = 0x0; + fLambdaAssocMassPtFidVolumeEmbeded[i] = 0x0; // -- Mass vs rapidity vs pt vs centrlaity fLambdaMassPtRap[i] = 0; // -- Splitting checks @@ -189,26 +205,34 @@ fAntiLambdaMass(0), fAntiLambdaMassEmbeded(0), fAntiLambdaMass2(0), fAntiLambdaM // AntiLambda fAntiLambdaMCPtRapVtx[i] = 0; fAntiLambdaMCPtRapVtxEmbeded[i] = 0; + fAntiLambdaMCPtRapPtDaugPt[i] = 0x0; + fAntiLambdaMCPtRapPtDaugPtEmbeded[i] = 0x0; fAntiLambdaMCPtPhiEta[i] = 0; fAntiLambdaAssocPtPhiEta[i] = 0; // -- Natural particles - fAntiLambdaAssocMassPtRap[i] = 0; - fAntiLambdaAssocMassPtRap2[i] = 0; - fAntiLambdaAssocMassPtVtx[i] = 0; - fAntiLambdaAssocMassPtDCADaug[i] = 0; - fAntiLambdaAssocMassPtCPA[i] = 0; - fAntiLambdaAssocMassPtDCAPV[i] = 0; - fAntiLambdaAssocMassPtDaugNClsTPC[i] = 0; - fAntiLambdaAssocMassPtShTPCcls[i] = 0; + fAntiLambdaAssocMassPtRap[i] = 0x0; + fAntiLambdaAssocMassPtRap2[i] = 0x0; + fAntiLambdaAssocMassPtVtx[i] = 0x0; + fAntiLambdaAssocMassPtDCADaug[i] = 0x0; + fAntiLambdaAssocMassPtCPA[i] = 0x0; + fAntiLambdaAssocMassPtDCAPV[i] = 0x0; + fAntiLambdaAssocMassPtDaugNClsTPC[i] = 0x0; + fAntiLambdaAssocMassPtShTPCcls[i] = 0x0; + fAntiLambdaAssocMassPtDaugPt[i] = 0x0; + fAntiLambdaAssocMassPtCtau[i] = 0x0; + fAntiLambdaAssocMassPtFidVolume[i] = 0x0; // -- Embeded particles - fAntiLambdaAssocMassPtRapEmbeded[i] = 0; - fAntiLambdaAssocMassPtRapEmbeded2[i] = 0; - fAntiLambdaAssocMassPtVtxEmbeded[i] = 0; - fAntiLambdaAssocMassPtDCADaug[i] = 0; - fAntiLambdaAssocMassPtCPAEmbeded[i] = 0; - fAntiLambdaAssocMassPtDCAPVEmbeded[i] = 0; - fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0; - fAntiLambdaAssocMassPtShTPCclsEmbeded[i] = 0; + fAntiLambdaAssocMassPtRapEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtRapEmbeded2[i] = 0x0; + fAntiLambdaAssocMassPtVtxEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtDCADaug[i] = 0x0; + fAntiLambdaAssocMassPtCPAEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtDCAPVEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtShTPCclsEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtDaugPtEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtCtauEmbeded[i] = 0x0; + fAntiLambdaAssocMassPtFidVolumeEmbeded[i] = 0x0; // -- Mass vs rapidity vs pt vs centrlaity fAntiLambdaMassPtRap[i] = 0; // -- Splitting checks @@ -545,6 +569,22 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fOutput->Add(fK0sMCPtRapVtxEmbeded[jj]); } + + Int_t binsK0sMC[4] = {nbins,20,100.,100.}; Double_t xminK0sMC[4] = {pMin,-1.0,0.,0.}; Double_t xmaxK0sMC[4] = {pMax,1.0,10.,10.}; // gral efficiency + for(Int_t jj=0;jjAdd(fK0sMCPtRapPtDaugPt[jj]); + + + snprintf(hNameHist,100, "fK0sMCPtRapPtDaugPtEmbeded_Cent_%d",jj); + fK0sMCPtRapPtDaugPtEmbeded[jj] = new THnSparseD(hNameHist, "K^{0}_{S} MC Embeded;#it{p}_{T} (GeV/#it{c});y;#it{p}_{T,Pos Daug} (GeV/#it{c});#it{p}_{T,Neg Daug} (GeV/#it{c});", + 4,binsK0sMC,xminK0sMC,xmaxK0sMC); + fOutput->Add(fK0sMCPtRapPtDaugPtEmbeded[jj]); + } + + for(Int_t jj=0;jjAdd(fK0sAssocPtRapEmbeded); - + for(Int_t jj=0;jjAdd(fK0sAssocMassPtShTPCcls[i]); + snprintf(hNameHist,100, "fK0sAssocMassPtDaugPt_Cent_%d",i); + fK0sAssocMassPtDaugPt[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffKPtDaug,xminEffKPtDaug,xmaxEffKPtDaug); + fOutput->Add(fK0sAssocMassPtDaugPt[i]); + + snprintf(hNameHist,100, "fK0sAssocMassPtCtau_Cent_%d",i); + fK0sAssocMassPtCtau[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau;",4,binsEffKCtau,xminEffKCtau,xmaxEffKCtau); + fOutput->Add(fK0sAssocMassPtCtau[i]); + + snprintf(hNameHist,100, "fK0sAssocMassPtFidVolume_Cent_%d",i); + fK0sAssocMassPtFidVolume[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffKFidVol,xminEffKFidVol,xmaxEffKFidVol); + fOutput->Add(fK0sAssocMassPtFidVolume[i]); + + /// ----- Embeded particles snprintf(hNameHist,100, "fK0sAssocPtMassArmEmbeded_Cent_%d",i); fK0sAssocPtMassArmEmbeded[i] = new THnSparseD(hNameHist,"K^{0}_{S} Assoc Embeded;Mass (GeV/c^{2});#it{p}_{T} (GeV/#it{c});rap;",3,binsEff1,xminEff1,xmaxEff1); @@ -638,9 +694,21 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fK0sAssocMassPtDaugNClsTPCEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos # TPC Cls; Neg # TPC Cls;",5,binsEff6,xminEff6,xmaxEff6); fOutput->Add(fK0sAssocMassPtDaugNClsTPCEmbeded[i]); - snprintf(hNameHist,100, "fK0sAssocMassPtShTPCclsEmbeded_Cent_%d",i); - fK0sAssocMassPtShTPCclsEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos fraction shared TPC Cls; Neg fraction shared TPC Cls;",5,binsEffKsh,xminEffKsh,xmaxEffKsh); - fOutput->Add(fK0sAssocMassPtShTPCclsEmbeded[i]); + snprintf(hNameHist,100, "fK0sAssocMassPtShTPCclsEmbeded_Cent_%d",i); + fK0sAssocMassPtShTPCclsEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos fraction shared TPC Cls; Neg fraction shared TPC Cls;",5,binsEffKsh,xminEffKsh,xmaxEffKsh); + fOutput->Add(fK0sAssocMassPtShTPCclsEmbeded[i]); + + snprintf(hNameHist,100, "fK0sAssocMassPtDaugPtEmbeded_Cent_%d",i); + fK0sAssocMassPtDaugPtEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffKPtDaug,xminEffKPtDaug,xmaxEffKPtDaug); + fOutput->Add(fK0sAssocMassPtDaugPtEmbeded[i]); + + snprintf(hNameHist,100, "fK0sAssocMassPtCtauEmbeded_Cent_%d",i); + fK0sAssocMassPtCtauEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau;",4,binsEffKCtau,xminEffKCtau,xmaxEffKCtau); + fOutput->Add(fK0sAssocMassPtCtauEmbeded[i]); + + snprintf(hNameHist,100, "fK0sAssocMassPtFidVolumeEmbeded_Cent_%d",i); + fK0sAssocMassPtFidVolumeEmbeded[i] = new THnSparseD(hNameHist, "K^{0}_{S}; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffKFidVol,xminEffKFidVol,xmaxEffKFidVol); + fOutput->Add(fK0sAssocMassPtFidVolumeEmbeded[i]); } @@ -696,6 +764,20 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fOutput->Add(fLambdaMCPtRapVtxEmbeded[jj]); } + Int_t binsLambdaMC[4] = {nbins,20,100.,100.}; Double_t xminLambdaMC[4] = {pMin,-1.0,0.,0.}; Double_t xmaxLambdaMC[4] = {pMax,1.0,10.,10.}; // gral efficiency + for(Int_t jj=0;jjAdd(fLambdaMCPtRapPtDaugPt[jj]); + + + snprintf(hNameHist,100, "fLambdaMCPtRapPtDaugPtEmbeded_Cent_%d",jj); + fLambdaMCPtRapPtDaugPtEmbeded[jj] = new THnSparseD(hNameHist, "#Lambda MC Embeded;#it{p}_{T} (GeV/#it{c});y;#it{p}_{T,Pos Daug} (GeV/#it{c});#it{p}_{T,Neg Daug} (GeV/#it{c});", + 4,binsLambdaMC,xminLambdaMC,xmaxLambdaMC); + fOutput->Add(fLambdaMCPtRapPtDaugPtEmbeded[jj]); + } + fLambdaMCFromXi = new TH2F("fLambdaMCFromXi", "#Lambda from Xi MC;#it{p}_{T} (GeV/#it{c});centrality",nbins,pMin,pMax,100,0.,100.); fOutput->Add(fLambdaMCFromXi); @@ -730,6 +812,10 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() Int_t binsEff11[5] = {nbins,nbins,20,99,99}; Double_t xminEff11[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEff11[5] = {1.165,pMax,1.0,3.3,3.3}; // dca to prim. vtx Int_t binsEff12[5] = {nbins,nbins,20,170,170}; Double_t xminEff12[5] = {1.065,pMin,-1.0,0.5,0.5}; Double_t xmaxEff12[5] = {1.165,pMax,1.0,170.5,170.5}; // No. TPC Cls Int_t binsEffLsh[5] = {nbins,nbins,20,50,50}; Double_t xminEffLsh[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEffLsh[5] = {1.165,pMax,1.0,1.0,1.0}; // shared TPC cls + Int_t binsEffLPtDaug[5] = {nbins,nbins,20,100,100}; Double_t xminEffLPtDaug[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEffLPtDaug[5] = {1.165,pMax,1.0,10.,10.}; //PtDaug + Int_t binsEffLCtau[4] = {nbins,nbins,20,60}; Double_t xminEffLCtau[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEffLCtau[4] = {1.165,pMax,1.0,30.0}; //CTau + Int_t binsEffLFidVol[4] = {nbins,nbins,20,100}; Double_t xminEffLFidVol[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEffLFidVol[4] = {1.165,pMax,1.0,100.}; //Fiducial volume + for(Int_t i=0; iAdd(fLambdaAssocMassPtShTPCcls[i]); + snprintf(hNameHist,100, "fLambdaAssocMassPtDaugPt_Cent_%d",i); + fLambdaAssocMassPtDaugPt[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffLPtDaug,xminEffLPtDaug,xmaxEffLPtDaug); + fOutput->Add(fLambdaAssocMassPtDaugPt[i]); + + snprintf(hNameHist,100, "fLambdaAssocMassPtCtau_Cent_%d",i); + fLambdaAssocMassPtCtau[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau;",4,binsEffLCtau,xminEffLCtau,xmaxEffLCtau); + fOutput->Add(fLambdaAssocMassPtCtau[i]); + + snprintf(hNameHist,100, "fLambdaAssocMassPtFidVolume_Cent_%d",i); + fLambdaAssocMassPtFidVolume[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffLFidVol,xminEffLFidVol,xmaxEffLFidVol); + fOutput->Add(fLambdaAssocMassPtFidVolume[i]); // ------------ Embeded particles snprintf(hNameHist,100, "fLambdaAssocMassPtRapEmbeded_Cent_%d",i); @@ -800,6 +897,18 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fLambdaAssocMassPtShTPCclsEmbeded[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos fraction shared TPC Cls; Neg fraction shared TPC Cls;",5,binsEffLsh,xminEffLsh,xmaxEffLsh); fOutput->Add(fLambdaAssocMassPtShTPCclsEmbeded[i]); + snprintf(hNameHist,100, "fLambdaAssocMassPtDaugPtEmbeded_Cent_%d",i); + fLambdaAssocMassPtDaugPtEmbeded[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffLPtDaug,xminEffLPtDaug,xmaxEffLPtDaug); + fOutput->Add(fLambdaAssocMassPtDaugPtEmbeded[i]); + + snprintf(hNameHist,100, "fLambdaAssocMassPtCtauEmbeded_Cent_%d",i); + fLambdaAssocMassPtCtauEmbeded[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau ;",4,binsEffLCtau,xminEffLCtau,xmaxEffLCtau); + fOutput->Add(fLambdaAssocMassPtCtauEmbeded[i]); + + snprintf(hNameHist,100, "fLambdaAssocMassPtFidVolumeEmbeded_Cent_%d",i); + fLambdaAssocMassPtFidVolumeEmbeded[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffLFidVol,xminEffLFidVol,xmaxEffLFidVol); + fOutput->Add(fLambdaAssocMassPtFidVolumeEmbeded[i]); + } fLambdaMCResEta = new TH3F("fLambdaMCResEta","#Lambda Assoc: #eta resolution; #eta_{MC}-#eta_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.); @@ -854,6 +963,20 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fOutput->Add(fAntiLambdaMCPtRapVtxEmbeded[jj]); } + Int_t binsAntiLambdaMC[4] = {nbins,20,100.,100.}; Double_t xminAntiLambdaMC[4] = {pMin,-1.0,0.,0.}; Double_t xmaxAntiLambdaMC[4] = {pMax,1.0,10.,10.}; // gral efficiency + for(Int_t jj=0;jjAdd(fAntiLambdaMCPtRapPtDaugPt[jj]); + + + snprintf(hNameHist,100, "fLambdaMCPtRapPtDaugPtEmbeded_Cent_%d",jj); + fAntiLambdaMCPtRapPtDaugPtEmbeded[jj] = new THnSparseD(hNameHist, "#bar{#Lambda} MC Embeded;#it{p}_{T} (GeV/#it{c});y;#it{p}_{T,Pos Daug} (GeV/#it{c});#it{p}_{T,Neg Daug} (GeV/#it{c});", + 4,binsAntiLambdaMC,xminAntiLambdaMC,xmaxAntiLambdaMC); + fOutput->Add(fAntiLambdaMCPtRapPtDaugPtEmbeded[jj]); + } + fAntiLambdaMCFromXi = new TH2F("fAntiLambdaMCFromXi", "#bar{#Lambda} from Xi MC;#it{p}_{T} (GeV/#it{c});centrality",nbins,pMin,pMax,100,0.,100.); fOutput->Add(fAntiLambdaMCFromXi); @@ -888,7 +1011,10 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() Int_t binsEff17[5] = {nbins,nbins,20,99,99}; Double_t xminEff17[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEff17[5] = {1.165,pMax,1.0,3.3,3.3}; // dca to prim. vtx Int_t binsEff18[5] = {nbins,nbins,20,170,170}; Double_t xminEff18[5] = {1.065,pMin,-1.0,0.5,0.5}; Double_t xmaxEff18[5] = {1.165,pMax,1.0,170.5,170.5}; // No. TPC Cls Int_t binsEffALsh[5] = {nbins,nbins,20,50,50}; Double_t xminEffALsh[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEffALsh[5] = {1.165,pMax,1.0,1.0,1.0}; // shared TPC cls - + Int_t binsEffALPtDaug[5] = {nbins,nbins,20,100,100}; Double_t xminEffALPtDaug[5] = {1.065,pMin,-1.0,0.,0.}; Double_t xmaxEffALPtDaug[5] = {1.065,pMax,1.0,10.,10.}; //PtDaug + Int_t binsEffALCtau[4] = {nbins,nbins,20,60}; Double_t xminEffALCtau[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEffALCtau[4] = {1.065,pMax,1.0,30.0}; //CTau + Int_t binsEffALFidVol[4] = {nbins,nbins,20,100}; Double_t xminEffALFidVol[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEffALFidVol[4] = {1.065,pMax,1.0,100.}; //Fiducial volume + for(Int_t i=0; iAdd(fAntiLambdaAssocMassPtShTPCcls[i]); + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtDaugPt_Cent_%d",i); + fAntiLambdaAssocMassPtDaugPt[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffALPtDaug,xminEffALPtDaug,xmaxEffALPtDaug); + fOutput->Add(fAntiLambdaAssocMassPtDaugPt[i]); + + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtCtau_Cent_%d",i); + fAntiLambdaAssocMassPtCtau[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau;",4,binsEffALCtau,xminEffALCtau,xmaxEffALCtau); + fOutput->Add(fAntiLambdaAssocMassPtCtau[i]); + + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtFidVolume_Cent_%d",i); + fAntiLambdaAssocMassPtFidVolume[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffALFidVol,xminEffALFidVol,xmaxEffALFidVol); + fOutput->Add(fAntiLambdaAssocMassPtFidVolume[i]); + // ------------ Embeded particles snprintf(hNameHist,100, "fAntiLambdaAssocMassPtRapEmbeded_Cent_%d",i); fAntiLambdaAssocMassPtRapEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda} Embeded; Mass (GeV/c^{2}); #it{p}_{T}; rap;",3,binsEff13,xminEff13,xmaxEff13); @@ -952,10 +1090,23 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos # TPC Cls; Neg # TPC Cls;",5,binsEff18,xminEff18,xmaxEff18); fOutput->Add(fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i]); - snprintf(hNameHist,100, "fAntiLambdaAssocMassPtShTPCclsEmbeded_Cent_%d",i); + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtShTPCclsEmbeded_Cent_%d",i); fAntiLambdaAssocMassPtShTPCclsEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos fraction shared TPC Cls; Neg fraction shared TPC Cls;",5,binsEffALsh,xminEffALsh,xmaxEffALsh); fOutput->Add(fAntiLambdaAssocMassPtShTPCclsEmbeded[i]); + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtDaugPtEmbeded_Cent_%d",i); + fAntiLambdaAssocMassPtDaugPtEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; #it{p}_{T,Pos Daug} (GeV/#it{c}); #it{p}_{T,Neg Daug} (GeV/#it{c});",5,binsEffALPtDaug,xminEffALPtDaug,xmaxEffALPtDaug); + fOutput->Add(fAntiLambdaAssocMassPtDaugPtEmbeded[i]); + + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtCtauEmbeded_Cent_%d",i); + fAntiLambdaAssocMassPtCtauEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; c#tau;",4,binsEffALCtau,xminEffALCtau,xmaxEffALCtau); + fOutput->Add(fAntiLambdaAssocMassPtCtauEmbeded[i]); + + snprintf(hNameHist,100, "fAntiLambdaAssocMassPtFidVolumeEmbeded_Cent_%d",i); + fAntiLambdaAssocMassPtFidVolumeEmbeded[i] = new THnSparseD(hNameHist, "#bar{#Lambda}; Mass (GeV/c^{2}); #it{p}_{T}; rap; l_{T} (cm);",4,binsEffALFidVol,xminEffALFidVol,xmaxEffALFidVol); + fOutput->Add(fAntiLambdaAssocMassPtFidVolumeEmbeded[i]); + + } fAntiLambdaMCResEta = new TH3F("fAntiLambdaMCResEta","#bar{#Lambda} Assoc: #eta resolution; #eta_{MC}-#eta_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.); @@ -1014,6 +1165,7 @@ void AliAnalysisTaskLambdaOverK0sJets::UserCreateOutputObjects() // ****** K0s ****** fK0sMass = new TH3F("fK0sMass", "K^{0}_{s}: mass vs #it{p}_{T};Mass (GeV/#it{c}^2);#it{p}_{T} (GeV/#it{c});centrality",nbins,0.398,0.598,nbins,pMin,pMax,100,0.,100.); + fOutput->Add(fK0sMass); fK0sMassEmbeded = @@ -3123,7 +3275,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge } /* - // Daughter momentum cut: ! FIX it in case of AOD ! //MC or REc +3 // Daughter momentum cut: ! FIX it in case of AOD ! //MC or REc if ( (nPart->Pt() < fMinPtDaughter ) || (pPart->Pt() < fMinPtDaughter ) ) goto noas; @@ -3227,7 +3379,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Natural particles if(isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fK0sAssocPt->Fill(ptAs); fK0sAssocPtRap->Fill(ptAs,rapAs,centrality); @@ -3241,6 +3393,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge Double_t effK0sDCA[4] = {v0->MassK0Short(),ptAs,rapAs,dca}; Double_t effK0sCPA[4] = {v0->MassK0Short(),ptAs,rapAs,cpa}; Double_t effK0sShTPCcls[5] = {v0->MassK0Short(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effK0sDaugPt[5] = {v0->MassK0Short(),ptAs,rapAs,lPtPos,lPtNeg}; + Double_t effK0sCtau[4] = {v0->MassK0Short(),ptAs,rapAs,dlK}; + Double_t effK0sFidVol[4] = {v0->MassK0Short(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fK0sAssocPtMassArm[curCentBin]->Fill(effK0sArm); @@ -3248,6 +3403,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fK0sAssocMassPtDCADaug[curCentBin]->Fill(effK0sDCA); fK0sAssocMassPtCPA[curCentBin]->Fill(effK0sCPA); fK0sAssocMassPtShTPCcls[curCentBin]->Fill(effK0sShTPCcls); + fK0sAssocMassPtDaugPt[curCentBin]->Fill(effK0sDaugPt); + fK0sAssocMassPtCtau[curCentBin]->Fill(effK0sCtau); + fK0sAssocMassPtFidVolume[curCentBin]->Fill(effK0sFidVol); } fK0sMCResEta->Fill(resEta,pt,centrality); @@ -3268,7 +3426,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( (lPtArmV0 > TMath::Abs(0.2*lAlphaV0) ) && TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effK0sdcaPV[5] = {v0->MassK0Short(),ptAs,rapAs,dcaPos,dcaNeg}; fK0sAssocMassPtDCAPV[curCentBin]->Fill(effK0sdcaPV); @@ -3288,7 +3446,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Embeded particles if(!isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fK0sAssocPtRapEmbeded->Fill(ptAs,rapAs,centrality); @@ -3299,6 +3457,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge Double_t effK0sDCA[4] = {v0->MassK0Short(),ptAs,rapAs,dca}; Double_t effK0sCPA[4] = {v0->MassK0Short(),ptAs,rapAs,cpa}; Double_t effK0sShTPCcls[5] = {v0->MassK0Short(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effK0sDaugPt[5] = {v0->MassK0Short(),ptAs,rapAs,lPtPos,lPtPos}; + Double_t effK0sCtau[4] = {v0->MassK0Short(),ptAs,rapAs,dlK}; + Double_t effK0sFidVol[4] = {v0->MassK0Short(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fK0sAssocPtMassArmEmbeded[curCentBin]->Fill(effK0sArm); @@ -3306,6 +3467,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fK0sAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effK0sDCA); fK0sAssocMassPtCPAEmbeded[curCentBin]->Fill(effK0sCPA); fK0sAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effK0sShTPCcls); + fK0sAssocMassPtDaugPtEmbeded[curCentBin]->Fill(effK0sDaugPt); + fK0sAssocMassPtCtauEmbeded[curCentBin]->Fill(effK0sCtau); + fK0sAssocMassPtFidVolumeEmbeded[curCentBin]->Fill(effK0sFidVol); } } // End selection in the dca to prim. vtx and the number of clusters @@ -3314,7 +3478,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( (lPtArmV0 > TMath::Abs(0.2*lAlphaV0) ) && TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effK0sdcaPV[5] = {v0->MassK0Short(),ptAs,rapAs,dcaPos,dcaNeg}; fK0sAssocMassPtDCAPVEmbeded[curCentBin]->Fill(effK0sdcaPV); @@ -3339,7 +3503,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Natural particles if(isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fLambdaAssocPt->Fill(ptAs); fLambdaAssocPtRap->Fill(ptAs,rapAs,centrality); @@ -3353,6 +3517,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge Double_t effLambdaDCA[4] = {v0->MassLambda(),ptAs,rapAs,dca}; Double_t effLambdaCPA[4] = {v0->MassLambda(),ptAs,rapAs,cpa}; Double_t effLambdaShTPCcls[5] = {v0->MassLambda(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effLambdaDaugPt[5] = {v0->MassLambda(),ptAs,rapAs,lPtPos,lPtNeg}; + Double_t effLambdaCtau[4] = {v0->MassLambda(),ptAs,rapAs,dlL}; + Double_t effLambdaFidVol[4] = {v0->MassLambda(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fLambdaAssocMassPtRap[curCentBin]->Fill(effLambda); @@ -3360,6 +3527,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fLambdaAssocMassPtDCADaug[curCentBin]->Fill(effLambdaDCA); fLambdaAssocMassPtCPA[curCentBin]->Fill(effLambdaCPA); fLambdaAssocMassPtShTPCcls[curCentBin]->Fill(effLambdaShTPCcls); + fLambdaAssocMassPtDaugPt[curCentBin]->Fill(effLambdaDaugPt); + fLambdaAssocMassPtCtau[curCentBin]->Fill(effLambdaCtau); + fLambdaAssocMassPtFidVolume[curCentBin]->Fill(effLambdaFidVol); if( !isCandidate2K0s && !isCandidate2LambdaBar) fLambdaAssocMassPtRap2[curCentBin]->Fill(effLambda); @@ -3384,7 +3554,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effLambdadcaPV[5] = {v0->MassLambda(),ptAs,rapAs,dcaPos,dcaNeg}; fLambdaAssocMassPtDCAPV[curCentBin]->Fill(effLambdadcaPV); @@ -3403,7 +3573,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Embeded particles if(!isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ if( TMath::Abs(rapAs)MassLambda(),ptAs,rapAs,dca}; Double_t effLambdaCPA[4] = {v0->MassLambda(),ptAs,rapAs,cpa}; Double_t effLambdaShTPCcls[5] = {v0->MassLambda(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effLambdaDaugPt[5] = {v0->MassLambda(),ptAs,rapAs,lPtPos,lPtNeg}; + Double_t effLambdaCtau[4] = {v0->MassLambda(),ptAs,rapAs,dlL}; + Double_t effLambdaFidVol[4] = {v0->MassLambda(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fLambdaAssocMassPtRapEmbeded[curCentBin]->Fill(effLambda); @@ -3419,6 +3592,10 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fLambdaAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effLambdaDCA); fLambdaAssocMassPtCPAEmbeded[curCentBin]->Fill(effLambdaCPA); fLambdaAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effLambdaShTPCcls); + fLambdaAssocMassPtDaugPtEmbeded[curCentBin]->Fill(effLambdaDaugPt); + fLambdaAssocMassPtCtauEmbeded[curCentBin]->Fill(effLambdaCtau); + fLambdaAssocMassPtFidVolumeEmbeded[curCentBin]->Fill(effLambdaFidVol); + if( !isCandidate2K0s && !isCandidate2LambdaBar) fLambdaAssocMassPtRapEmbeded2[curCentBin]->Fill(effLambda); @@ -3430,7 +3607,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effLambdadcaPV[5] = {v0->MassLambda(),ptAs,rapAs,dcaPos,dcaNeg}; fLambdaAssocMassPtDCAPVEmbeded[curCentBin]->Fill(effLambdadcaPV); } @@ -3452,7 +3629,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if(isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fAntiLambdaAssocPt->Fill(ptAs); fAntiLambdaAssocPtRap->Fill(ptAs,rapAs,centrality); @@ -3466,6 +3643,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge Double_t effAntiLambdaDCA[4] = {v0->MassAntiLambda(),ptAs,rapAs,dca}; Double_t effAntiLambdaCPA[4] = {v0->MassAntiLambda(),ptAs,rapAs,cpa}; Double_t effAntiLambdaShTPCcls[5] = {v0->MassAntiLambda(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effAntiLambdaDaugPt[5] = {v0->MassAntiLambda(),ptAs,rapAs,lPtPos,lPtNeg}; + Double_t effAntiLambdaCtau[4] = {v0->MassAntiLambda(),ptAs,rapAs,dlL}; + Double_t effAntiLambdaFidVol[4] = {v0->MassAntiLambda(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fAntiLambdaAssocMassPtRap[curCentBin]->Fill(effAntiLambda); @@ -3473,7 +3653,11 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fAntiLambdaAssocMassPtDCADaug[curCentBin]->Fill(effAntiLambdaDCA); fAntiLambdaAssocMassPtCPA[curCentBin]->Fill(effAntiLambdaCPA); fAntiLambdaAssocMassPtShTPCcls[curCentBin]->Fill(effAntiLambdaShTPCcls); - + fAntiLambdaAssocMassPtDaugPt[curCentBin]->Fill(effAntiLambdaDaugPt); + fAntiLambdaAssocMassPtCtau[curCentBin]->Fill(effAntiLambdaCtau); + fAntiLambdaAssocMassPtFidVolume[curCentBin]->Fill(effAntiLambdaFidVol); + + if( !isCandidate2K0s && !isCandidate2Lambda ) fAntiLambdaAssocMassPtRap2[curCentBin]->Fill(effAntiLambda); } @@ -3496,7 +3680,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effAntiLambdadcaPV[5] = {v0->MassAntiLambda(),ptAs,rapAs,dcaPos,dcaNeg}; fAntiLambdaAssocMassPtDCAPV[curCentBin]->Fill(effAntiLambdadcaPV); @@ -3514,7 +3698,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Embeded particles if(!isNaturalPart){ - if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ if( TMath::Abs(rapAs)MassAntiLambda(),ptAs,rapAs,dca}; Double_t effAntiLambdaCPA[4] = {v0->MassAntiLambda(),ptAs,rapAs,cpa}; Double_t effAntiLambdaShTPCcls[5] = {v0->MassAntiLambda(),ptAs,rapAs,fracPosDaugTPCSharedMap,fracNegDaugTPCSharedMap}; + Double_t effAntiLambdaDaugPt[5] = {v0->MassAntiLambda(),ptAs,rapAs,lPtPos,lPtNeg}; + Double_t effAntiLambdaCtau[4] = {v0->MassAntiLambda(),ptAs,rapAs,dlL}; + Double_t effAntiLambdaFidVol[4] = {v0->MassAntiLambda(),ptAs,rapAs,lt}; // Distributions for the efficiency (systematics chechks) fAntiLambdaAssocMassPtRapEmbeded[curCentBin]->Fill(effAntiLambda); @@ -3530,6 +3717,9 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge fAntiLambdaAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effAntiLambdaDCA); fAntiLambdaAssocMassPtCPAEmbeded[curCentBin]->Fill(effAntiLambdaCPA); fAntiLambdaAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effAntiLambdaShTPCcls); + fAntiLambdaAssocMassPtDaugPtEmbeded[curCentBin]->Fill(effAntiLambdaDaugPt); + fAntiLambdaAssocMassPtCtauEmbeded[curCentBin]->Fill(effAntiLambdaCtau); + fAntiLambdaAssocMassPtFidVolumeEmbeded[curCentBin]->Fill(effAntiLambdaFidVol); if( !isCandidate2K0s && !isCandidate2Lambda ) fAntiLambdaAssocMassPtRapEmbeded2[curCentBin]->Fill(effAntiLambda); @@ -3542,7 +3732,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( TMath::Abs(rapAs)fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ Double_t effAntiLambdadcaPV[5] = {v0->MassAntiLambda(),ptAs,rapAs,dcaPos,dcaNeg}; fAntiLambdaAssocMassPtDCAPVEmbeded[curCentBin]->Fill(effAntiLambdadcaPV); @@ -3600,7 +3790,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge switch(step) { case kTriggerCheck: - if (isCandidate2K0s && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if (isCandidate2K0s && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) ){ if(pt>ptTrig){ fIsV0LP = 1; @@ -3659,7 +3849,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge break; // End K0s selection for TriggerCheck case kReconstruction: - if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) && (pt<10.) ){ + if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) && (pt<10.) ){ if(isNaturalPart) fK0sMass->Fill(massK0s,pt,centrality); else fK0sMassEmbeded->Fill(massK0s,pt,centrality); @@ -3682,14 +3872,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Invariant Mass cut if (TMath::Abs(mK0s-massK0s) < 3*sK0s) { - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fK0sDCAPosDaug->Fill(dcaPos,pt); fK0sDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos > 0.1) && (dcaNeg > 0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC ) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70 ) ){ fK0sPtPosDaug->Fill(pt,lPtPos); fK0sPtNegDaug->Fill(pt,lPtNeg); @@ -3715,14 +3905,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( TMath::Abs(mK0s-massK0s + 6.5*sK0s) < 1.5*sK0s || TMath::Abs(mK0s-massK0s - 6.5*sK0s) < 1.5*sK0s ) { - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fK0sBckgDCAPosDaug->Fill(dcaPos,pt); fK0sBckgDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos > 0.1) && (dcaNeg > 0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fK0sBckgPtPosDaug->Fill(pt,lPtPos); fK0sBckgPtNegDaug->Fill(pt,lPtNeg); @@ -3765,7 +3955,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge switch(step) { case kTriggerCheck: - if (isCandidate2Lambda && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) && !isCandidate2K0s && !isCandidate2LambdaBar ){ + if (isCandidate2Lambda && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) && !isCandidate2K0s && !isCandidate2LambdaBar ){ if(pt>ptTrig) { fIsV0LP = 1; @@ -3821,7 +4011,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge break; // End Lambda selection for TriggerCheck case kReconstruction: - if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) && (pt<10.) ){ + if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) && (pt<10.) ){ if(isNaturalPart) fLambdaMass->Fill(massLambda,pt,centrality); else fLambdaMassEmbeded->Fill(massLambda,pt,centrality); @@ -3851,14 +4041,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Invariant Mass cut if (TMath::Abs(mLambda-massLambda) < 3*sL) { - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fLambdaDCAPosDaug->Fill(dcaPos,pt); fLambdaDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos > 0.1) && (dcaNeg > 0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fLambdaPtPosDaug->Fill(pt,lPtPos); fLambdaPtNegDaug->Fill(pt,lPtNeg); @@ -3885,14 +4075,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( (TMath::Abs(mLambda-massLambda + 6.5*sL) < 1.5*sL) || (TMath::Abs(mLambda-massLambda - 6.5*sL) < 1.5*sL) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fLambdaBckgDCAPosDaug->Fill(dcaPos,pt); fLambdaBckgDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos > 0.1) && (dcaNeg > 0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fLambdaBckgPtPosDaug->Fill(pt,lPtPos); fLambdaBckgPtNegDaug->Fill(pt,lPtNeg); @@ -3934,7 +4124,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge switch(step) { case kTriggerCheck: - if (isCandidate2LambdaBar && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) && !isCandidate2K0s && !isCandidate2Lambda ){ + if (isCandidate2LambdaBar && (dcaPos>0.1) && (dcaNeg>0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) && !isCandidate2K0s && !isCandidate2Lambda ){ if(pt>ptTrig) { fIsV0LP = 1; @@ -3991,7 +4181,7 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge break; // End AntiLambda selection for CheckTrigger case kReconstruction: - if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) && (pt<10.) ) { + if( (dcaPos > 0.1) && (dcaNeg > 0.1) && (nClsTPCPos>70) && (nClsTPCNeg>70) && (pt<10.) ) { if(isNaturalPart) fAntiLambdaMass->Fill(massAntiLambda,pt,centrality); else fAntiLambdaMassEmbeded->Fill(massAntiLambda,pt,centrality); @@ -4019,14 +4209,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge // Invariant Mass cut if (TMath::Abs(mLambda-massAntiLambda) < 3*sAL) { - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fAntiLambdaDCAPosDaug->Fill(dcaPos,pt); fAntiLambdaDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos>0.1) && (dcaNeg>0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fAntiLambdaPtPosDaug->Fill(pt,lPtPos); fAntiLambdaPtNegDaug->Fill(pt,lPtNeg); @@ -4052,14 +4242,14 @@ void AliAnalysisTaskLambdaOverK0sJets::V0Loop(V0LoopStep_t step, Bool_t isTrigge if( (TMath::Abs(mLambda-massAntiLambda + 6.5*sAL) < 1.5*sAL) || (TMath::Abs(mLambda-massAntiLambda - 6.5*sAL) < 1.5*sAL) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fAntiLambdaBckgDCAPosDaug->Fill(dcaPos,pt); fAntiLambdaBckgDCANegDaug->Fill(dcaNeg,pt); } if( (dcaPos>0.1) && (dcaNeg>0.1) ){ - if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){ + if( (nClsTPCPos>70) && (nClsTPCNeg>70) ){ fAntiLambdaBckgPtPosDaug->Fill(pt,lPtPos); fAntiLambdaBckgPtNegDaug->Fill(pt,lPtNeg); diff --git a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.h b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.h index f8aa84b50d3..cd927afe960 100644 --- a/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.h +++ b/PWGLF/STRANGENESS/Correlations/AliAnalysisTaskLambdaOverK0sJets.h @@ -235,6 +235,8 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { TH3F* fK0sMCPtRapVtx[kNCent]; //! K0s MC: pt vs Z vtx position vs centrality TH3F* fK0sMCPtRapEmbeded; //! K0s MC: pt vs rapidity (embeded particles) TH3F* fK0sMCPtRapVtxEmbeded[kNCent]; //! K0s MC: pt vs Z vtx position rapidity vs centrality (embeded particles) + THnSparse* fK0sMCPtRapPtDaugPt[kNCent]; //! K0s MC: pt vs rapidity vs transverse momemtum of daughters + THnSparse* fK0sMCPtRapPtDaugPtEmbeded[kNCent]; //! K0s MC: pt vs rapidity vs transverse momemtum of daughters TH3F* fK0sMCPtPhiEta[kNCent]; //! K0s MC: pt vs pseudo-rapidity TH1F* fK0sAssocPt; //! K0s Assoc: pt @@ -250,6 +252,9 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fK0sAssocMassPtDCAPV[kNCent]; //! K0s Assoc: mass vs pt vs dca to prim. vtx THnSparse* fK0sAssocMassPtDaugNClsTPC[kNCent]; //! K0s Assoc: mass vs pt vs num. of tpc clusters THnSparse* fK0sAssocMassPtShTPCcls[kNCent]; //! K0s Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fK0sAssocMassPtDaugPt[kNCent]; //! K0s Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fK0sAssocMassPtCtau[kNCent]; //! K0s Assoc: mass vs pt vs proper time of life + THnSparse* fK0sAssocMassPtFidVolume[kNCent]; //! K0s Assoc: mass vs pt vs fiducial volume THnSparse* fK0sAssocPtMassArmEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs rapidity (embeded particles) THnSparse* fK0sAssocMassPtVtxEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs Z vertex position (embeded particles) @@ -258,6 +263,10 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fK0sAssocMassPtDCAPVEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs dca to prim. vtx (embeded particles) THnSparse* fK0sAssocMassPtDaugNClsTPCEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs num. o ftpc clusters (embeded particles) THnSparse* fK0sAssocMassPtShTPCclsEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fK0sAssocMassPtDaugPtEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fK0sAssocMassPtCtauEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs proper time of life + THnSparse* fK0sAssocMassPtFidVolumeEmbeded[kNCent]; //! K0s Assoc: mass vs pt vs fiducial volume + TH3F* fK0sMCResEta; //! K0s Assoc: eta resolution TH3F* fK0sMCResPhi; //! K0s Assoc: phi resolution @@ -276,6 +285,8 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { TH3F* fLambdaMCPtRapVtx[kNCent]; //! Lambda MC: pt vs Z vtx position rapidity vs centrality TH3F* fLambdaMCPtRapEmbeded; //! Lambda MC: pt vs rapidity (embeded particles) TH3F* fLambdaMCPtRapVtxEmbeded[kNCent]; //! Lambda MC: pt vs Z vtx position vs centrality (embeded particles) + THnSparse* fLambdaMCPtRapPtDaugPt[kNCent]; //! Lambda MC: pt vs rapidity vs transverse momemtum of daughters + THnSparse* fLambdaMCPtRapPtDaugPtEmbeded[kNCent]; //! Lambda MC: pt vs rapidity vs transverse momemtum of daughters TH2F* fLambdaMCFromXi; //! Lambda MC: coming from Xi TH3F* fLambdaMCPtPhiEta[kNCent]; //! Lambda MC: pt vs pseudo-rapidity @@ -292,6 +303,9 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fLambdaAssocMassPtDCAPV[kNCent]; //! Lambda Assoc: mass vs pt vs dca to prim vtx THnSparse* fLambdaAssocMassPtDaugNClsTPC[kNCent]; //! Lambda Assoc: mass vs pt vs num.of tpc clusters THnSparse* fLambdaAssocMassPtShTPCcls[kNCent]; //! Lambda Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fLambdaAssocMassPtDaugPt[kNCent]; //! Lambda Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fLambdaAssocMassPtCtau[kNCent]; //! Lambda Assoc: mass vs pt vs proper time of life + THnSparse* fLambdaAssocMassPtFidVolume[kNCent]; //! Lambda Assoc: mass vs pt vs fiducial volume THnSparse* fLambdaAssocMassPtRapEmbeded[kNCent]; //! Lambda Assoc: pt vs rapidity vs mass (embeded) THnSparse* fLambdaAssocMassPtRapEmbeded2[kNCent]; //! Lambda Assoc: pt vs rapidity vs mass (wo Cross contamination) (embeded) @@ -301,6 +315,10 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fLambdaAssocMassPtDCAPVEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs dca to prim vtx (embeded particles) THnSparse* fLambdaAssocMassPtDaugNClsTPCEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs num. of tpc clusters (embeded particles) THnSparse* fLambdaAssocMassPtShTPCclsEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fLambdaAssocMassPtDaugPtEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fLambdaAssocMassPtCtauEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs proper time of life + THnSparse* fLambdaAssocMassPtFidVolumeEmbeded[kNCent]; //! Lambda Assoc: mass vs pt vs fiducial volume + TH3F* fLambdaMCResEta; //! Lambda Assoc: eta resolution TH3F* fLambdaMCResPhi; //! Lambda Assoc: phi resolution @@ -319,6 +337,8 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { TH3F* fAntiLambdaMCPtRapVtx[kNCent]; //! AntiLambda MC: pt vs rapidity vs Z vtx position TH3F* fAntiLambdaMCPtRapEmbeded; //! AntiLambda MC: pt vs rapidity (embeded particles) TH3F* fAntiLambdaMCPtRapVtxEmbeded[kNCent]; //! AntiLambda MC: pt vs rapidity vs Z vtx position + THnSparse* fAntiLambdaMCPtRapPtDaugPt[kNCent]; //! AntiLambda MC: pt vs rapidity vs transverse momemtum of daughters + THnSparse* fAntiLambdaMCPtRapPtDaugPtEmbeded[kNCent]; //! AntiLambda MC: pt vs rapidity vs transverse momemtum of daughters TH2F* fAntiLambdaMCFromXi; //! AntiLambda MC: coming from Xi TH3F* fAntiLambdaMCPtPhiEta[kNCent]; //! AntiLambda MC: pt vs pseudo-rapidity @@ -335,6 +355,9 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fAntiLambdaAssocMassPtDCAPV[kNCent]; //! AntiLambda Assoc: mass vs pt vs dca to prim. vtx THnSparse* fAntiLambdaAssocMassPtDaugNClsTPC[kNCent]; //! AntiLambda Assoc: mass vs pt vs num. of tpc clusters THnSparse* fAntiLambdaAssocMassPtShTPCcls[kNCent]; //! AntiLambda Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fAntiLambdaAssocMassPtDaugPt[kNCent]; //! AntiLambda Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fAntiLambdaAssocMassPtCtau[kNCent]; //! AntiLambda Assoc: mass vs pt vs proper time of life + THnSparse* fAntiLambdaAssocMassPtFidVolume[kNCent]; //! AntiLambda Assoc: mass vs pt vs fiducial volume THnSparse* fAntiLambdaAssocMassPtRapEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs rapidity (embeded) THnSparse* fAntiLambdaAssocMassPtRapEmbeded2[kNCent]; //! AntiLambda Assoc: mass vs pt vs rapidity (wo Cross contamination) (embeded) @@ -344,6 +367,10 @@ class AliAnalysisTaskLambdaOverK0sJets : public AliAnalysisTaskSE { THnSparse* fAntiLambdaAssocMassPtDCAPVEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs dca to prim. vtx (embeded particles) THnSparse* fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs num. of tpc clusters (embeded particles) THnSparse* fAntiLambdaAssocMassPtShTPCclsEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs fraction of shared TPC cls + THnSparse* fAntiLambdaAssocMassPtDaugPtEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs transverse momemtum of daughters + THnSparse* fAntiLambdaAssocMassPtCtauEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs proper time of life + THnSparse* fAntiLambdaAssocMassPtFidVolumeEmbeded[kNCent]; //! AntiLambda Assoc: mass vs pt vs fiducial volume + TH3F* fAntiLambdaMCResEta; //! AntiLambda Assoc: eta resolution TH3F* fAntiLambdaMCResPhi; //! AntiLambda Assoc: phi resolution diff --git a/PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C b/PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C index 34062059900..8784f5062b6 100644 --- a/PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C +++ b/PWGLF/STRANGENESS/Correlations/macros/AddTaskLambdaOverK0sJets.C @@ -20,7 +20,7 @@ AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString name = Float_t nSigmaPID = 3.0; Float_t dcaDaug = 1.0; Float_t dca2PrmVtx = 0.095; // tighter cut - Float_t nclsDaug = 70; + Float_t nclsDaug = 50; Float_t minPtDaughter = 0.16; Float_t radiusTPC = 125.; @@ -47,10 +47,10 @@ AliAnalysisTaskLambdaOverK0sJets *AddTaskLambdaOverK0sJets( TString name = task->SetEtaCut(useEtaCut); task->SetMaxY(rapMaxV0); task->SetMaxDCADaughter(dcaDaug); // Added to perform systematics - task->SetDCAToPrimVtx(dca2PrmVtx); // Added to perform systematics - task->SetNSigmaPID(nSigmaPID); + //task->SetDCAToPrimVtx(dca2PrmVtx); // Added to perform systematics + //task->SetNSigmaPID(nSigmaPID); task->SetNClsTPC(nclsDaug); // Added to perform systematics - task->SetMinPtDaughter(minPtDaughter); + //task->SetMinPtDaughter(minPtDaughter); // PID task->SetSeparateInjectedPart(sepInjec); // MC -- 2.39.3