AliAnalysisTaskLambdaOverK0sJets::AliAnalysisTaskLambdaOverK0sJets(const char *name) :
AliAnalysisTaskSE(name),
-fAOD(0), fCollision("PbPb2010"), fIsMC(kFALSE), fUsePID(kFALSE), fCentMin(0.), fCentMax(90.), fDoQA(kFALSE), fDoMixEvt(kFALSE), fTriggerFB(128), fTrigPtMin(5.), fTrigPtMax(10.), fTrigPtMCMin(5.), fTrigPtMCMax(10000.), fTrigEtaMax(0.8), fCheckIDTrig(kFALSE), fSeparateInjPart(kTRUE), fEndOfHijingEvent(-1), fPIDResponse(0),
+fAOD(0), fCollision("PbPb2010"), fIsMC(kFALSE), fUsePID(kFALSE), fCentMin(0.), fCentMax(90.), fDoQA(kFALSE), fDoMixEvt(kFALSE), fTriggerFB(768), fTrigPtMin(5.), fTrigPtMax(10.), fTrigPtMCMin(5.), fTrigPtMCMax(10000.), fTrigEtaMax(0.8), fCheckIDTrig(kFALSE), fSeparateInjPart(kTRUE), fEndOfHijingEvent(-1), fPIDResponse(0),
fMinPtDaughter(0.160), fMaxEtaDaughter(0.8), fMaxDCADaughter(1.0), fUseEtaCut(kFALSE), fYMax(0.7), fDCAToPrimVtx(0.1), fMinCPA(0.998), fNSigma(3.0), fDaugNClsTPC(70.), fMinCtau(0.), fMaxCtau(3.), fIdTrigger(-1), fIsV0LP(0), fPtV0LP(0.), fIsSndCheck(0),
fOutput(0), fOutputQA(0), fOutputME(0), fMEList(0x0), fTriggerParticles(0x0), fTriggerPartMC(0x0), fAssocParticles(0x0), fAssocPartMC(0x0), fEvents(0), fEvtPerCent(0), fCentrality(0), fCentrality2(0), fCentralityTrig(0), fPrimayVtxGlobalvsSPD(0), fPrimaryVertexX(0), fPrimaryVertexY(0), fPrimaryVertexZ(0),
-fTriggerEventPlane(0), fTriggerMCPtCent(0), fTriggerMCResPt(0), fTriggerMCResEta(0), fTriggerMCResPhi(0), fTriggerPtCent(0), fTriggerPtCentCh(0), fNTrigPerEvt(0), fTriggerWiSPDHit(0), fTriggerEtaPhi(0), fTriggerDCA(0), fCheckTriggerFromV0Daug(0), fTriggerComingFromDaug(0), fTriggerIsV0(0), fCheckIDTrigPtK0s(0), fCheckIDTrigPhiK0s(0), fCheckIDTrigEtaK0s(0), fCheckIDTrigNclsK0s(0), fCheckIDTrigPtLambda(0), fCheckIDTrigPhiLambda(0), fCheckIDTrigEtaLambda(0), fCheckIDTrigNclsLambda(0), fCheckIDTrigPtAntiLambda(0), fCheckIDTrigPhiAntiLambda(0), fCheckIDTrigEtaAntiLambda(0), fCheckIDTrigNclsAntiLambda(0),
+fTriggerEventPlane(0), fTriggerMCPtCent(0), fTriggerMCResPt(0), fTriggerMCResEta(0), fTriggerMCResPhi(0), fTriggerPtCent(0), fTriggerPtCentCh(0), fNTrigPerEvt(0), fTriggerWiSPDHit(0), fTriggerEtaPhi(0), fTrigFracShTPCcls(0), fTriggerDCA(0), fCheckTriggerFromV0Daug(0), fTriggerComingFromDaug(0), fTriggerIsV0(0), fCheckIDTrigPtK0s(0), fCheckIDTrigPhiK0s(0), fCheckIDTrigEtaK0s(0), fCheckIDTrigNclsK0s(0), fCheckIDTrigPtLambda(0), fCheckIDTrigPhiLambda(0), fCheckIDTrigEtaLambda(0), fCheckIDTrigNclsLambda(0), fCheckIDTrigPtAntiLambda(0), fCheckIDTrigPhiAntiLambda(0), fCheckIDTrigEtaAntiLambda(0), fCheckIDTrigNclsAntiLambda(0),
fInjectedParticles(0),
fHistArmenterosPodolanski(0), fHistArmPodBckg(0),
- fK0sMass(0), fK0sMassEmbeded(0), fK0sMassPtEta(0), fK0sMassPtPhi(0), fK0sDaughtersPt(0), fK0sDCADaugToPrimVtx(0), fK0sSpatialRes(0), fK0sBckgDecLength(0), fK0sBckgDCADaugToPrimVtx(0), fK0sBckgEtaPhi(0), fK0sBckgPhiRadio(0), fK0sBckgDCANegDaugToPrimVtx(0), fK0sBckgDCAPosDaugToPrimVtx(0), fV0MassCascade(0),
+fK0sMass(0), fK0sMassEmbeded(0), fK0sMassPtEta(0), fK0sMassPtPhi(0), fK0sPosDaugFracShTPCcls(0), fK0sNegDaugFracShTPCcls(0), fK0sDaughtersPt(0), fK0sPosDaugFracShTPCclsTrig(0), fK0sNegDaugFracShTPCclsTrig(0), fK0sDCADaugToPrimVtx(0), fK0sSpatialRes(0), fK0sBckgDecLength(0), fK0sBckgDCADaugToPrimVtx(0), fK0sBckgEtaPhi(0), fK0sBckgPhiRadio(0), fK0sBckgDCANegDaugToPrimVtx(0), fK0sBckgDCAPosDaugToPrimVtx(0), fV0MassCascade(0),
- fLambdaMass(0), fLambdaMassEmbeded(0), fLambdaMass2(0), fLambdaMass2Embeded(0), fLambdaMassPtEta(0), fLambdaMassPtPhi(0), fLambdaDaughtersPt(0), fLambdaDCADaugToPrimVtx(0), fLambdaSpatialRes(0), fLambdaBckgDecLength(0), fLambdaBckgDCADaugToPrimVtx(0), fLambdaBckgEtaPhi(0), fLambdaBckgPhiRadio(0), fLambdaBckgDCANegDaugToPrimVtx(0), fLambdaBckgDCAPosDaugToPrimVtx(0),
+fLambdaMass(0), fLambdaMassEmbeded(0), fLambdaMass2(0), fLambdaMass2Embeded(0), fLambdaMassPtEta(0), fLambdaMassPtPhi(0), fLambdaPosDaugFracShTPCcls(0), fLambdaNegDaugFracShTPCcls(0), fLambdaDaughtersPt(0), fLambdaPosDaugFracShTPCclsTrig(0), fLambdaNegDaugFracShTPCclsTrig(0), fLambdaDCADaugToPrimVtx(0), fLambdaSpatialRes(0), fLambdaBckgDecLength(0), fLambdaBckgDCADaugToPrimVtx(0), fLambdaBckgEtaPhi(0), fLambdaBckgPhiRadio(0), fLambdaBckgDCANegDaugToPrimVtx(0), fLambdaBckgDCAPosDaugToPrimVtx(0),
- fAntiLambdaMass(0), fAntiLambdaMassEmbeded(0), fAntiLambdaMass2(0), fAntiLambdaMass2Embeded(0), fAntiLambdaMassPtEta(0), fAntiLambdaMassPtPhi(0), fAntiLambdaDaughtersPt(0), fAntiLambdaDCADaugToPrimVtx(0), fAntiLambdaSpatialRes(0), fAntiLambdaBckgDecLength(0), fAntiLambdaBckgDCADaugToPrimVtx(0), fAntiLambdaBckgEtaPhi(0), fAntiLambdaBckgPhiRadio(0), fAntiLambdaBckgDCANegDaugToPrimVtx(0), fAntiLambdaBckgDCAPosDaugToPrimVtx(0),
+fAntiLambdaMass(0), fAntiLambdaMassEmbeded(0), fAntiLambdaMass2(0), fAntiLambdaMass2Embeded(0), fAntiLambdaMassPtEta(0), fAntiLambdaMassPtPhi(0), fAntiLambdaPosDaugFracShTPCcls(0), fAntiLambdaNegDaugFracShTPCcls(0), fAntiLambdaDaughtersPt(0), fAntiLambdaPosDaugFracShTPCclsTrig(0), fAntiLambdaNegDaugFracShTPCclsTrig(0), fAntiLambdaDCADaugToPrimVtx(0), fAntiLambdaSpatialRes(0), fAntiLambdaBckgDecLength(0), fAntiLambdaBckgDCADaugToPrimVtx(0), fAntiLambdaBckgEtaPhi(0), fAntiLambdaBckgPhiRadio(0), fAntiLambdaBckgDCANegDaugToPrimVtx(0), fAntiLambdaBckgDCAPosDaugToPrimVtx(0),
fK0sPtPosDaug(0), fK0sPtNegDaug(0), fK0sBckgPtPosDaug(0), fK0sBckgPtNegDaug(0), fK0sPhiEtaPosDaug(0), fK0sPhiEtaNegDaug(0), fK0sBckgPhiEtaPosDaug(0), fK0sBckgPhiEtaNegDaug(0), fK0sDCAPosDaug(0), fK0sDCANegDaug(0), fK0sBckgDCAPosDaug(0), fK0sBckgDCANegDaug(0), fK0sDecayPos(0), fK0sBckgDecayPos(0), fK0sDecayVertex(0), fK0sBckgDecayVertex(0), fK0sCPA(0), fK0sBckgCPA(0), fK0sDCAV0Daug(0), fK0sBckgDCAV0Daug(0), fK0sNClustersTPC(0), fK0sBckgNClustersTPC(0), fK0sNClustersITSPos(0), fK0sNClustersITSNeg(0), fK0sBckgNClustersITSPos(0), fK0sBckgNClustersITSNeg(0),
fK0sAssocMassPtCPA[i] = 0;
fK0sAssocMassPtDCAPV[i] = 0;
fK0sAssocMassPtDaugNClsTPC[i] = 0;
+ fK0sAssocMassPtShTPCcls[i] = 0;
// -- Embeded particles
fK0sAssocPtMassArmEmbeded[i] = 0;
fK0sAssocMassPtVtxEmbeded[i] = 0;
fK0sAssocMassPtCPAEmbeded[i] = 0;
fK0sAssocMassPtDCAPVEmbeded[i] = 0;
fK0sAssocMassPtDaugNClsTPCEmbeded[i] = 0;
+ fK0sAssocMassPtShTPCclsEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fK0sMassPtRap[i] = 0;
// -- Splitting checks
fLambdaAssocMassPtCPA[i] = 0;
fLambdaAssocMassPtDCAPV[i] = 0;
fLambdaAssocMassPtDaugNClsTPC[i] = 0;
+ fLambdaAssocMassPtShTPCcls[i] = 0;
// -- Embeded particles
fLambdaAssocMassPtRapEmbeded[i] = 0;
fLambdaAssocMassPtRapEmbeded2[i] = 0;
fLambdaAssocMassPtCPAEmbeded[i] = 0;
fLambdaAssocMassPtDCAPVEmbeded[i] = 0;
fLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0;
+ fLambdaAssocMassPtShTPCclsEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fLambdaMassPtRap[i] = 0;
// -- Splitting checks
fAntiLambdaAssocMassPtCPA[i] = 0;
fAntiLambdaAssocMassPtDCAPV[i] = 0;
fAntiLambdaAssocMassPtDaugNClsTPC[i] = 0;
+ fAntiLambdaAssocMassPtShTPCcls[i] = 0;
// -- Embeded particles
fAntiLambdaAssocMassPtRapEmbeded[i] = 0;
fAntiLambdaAssocMassPtRapEmbeded2[i] = 0;
fAntiLambdaAssocMassPtCPAEmbeded[i] = 0;
fAntiLambdaAssocMassPtDCAPVEmbeded[i] = 0;
fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0;
+ fAntiLambdaAssocMassPtShTPCclsEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fAntiLambdaMassPtRap[i] = 0;
// -- Splitting checks
fOutput->Add(fEvtPerCent);
// Centrality:
- fCentrality = new TH1F("fCentrality","Centrality;Centrality (%);Events",100,0.,100.);
+ fCentrality = new TH1F("fCentrality","Centrality;Centrality (%);Events",500,0.,100.);
fOutput->Add(fCentrality);
- fCentrality2 = new TH1F("fCentrality2","Centrality in events with |VtxZ|<10 cm;Centrality (%);Events",100,0.,100.);
+ fCentrality2 = new TH1F("fCentrality2","Centrality in events with |VtxZ|<10 cm;Centrality (%);Events",500,0.,100.);
fOutput->Add(fCentrality2);
fCentralityTrig = new TH2F("fCentralityTrig","Centrality in events per trigger selection;Centrality (%);Triger Selection",100,0.,100.,3,0.5,3.5);
fTriggerDCA->GetYaxis()->SetBinLabel(2,"Z");
fOutput->Add(fTriggerDCA);
+ // Fraction of shared TPC cls
+ fTrigFracShTPCcls =
+ new TH2F("fTrigFracShTPCcls","Trigger particle: vs #it{p}_{T} vs fraction Shared TPC cls;#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbinPtLP,pMin,ptMaxLP,50,0,1.);
+ fOutput->Add(fTrigFracShTPCcls);
+
// Check if Trigger particle comes from a V0 daughter:
fCheckTriggerFromV0Daug =
new TH1F("fCheckTriggerFromV0Daug","Trigger particle from a V0 daughter;;Counts",4,-0.5,3.5);
// Histogramas para estudios sistematicos de la eficiencia
Int_t binsEff1[3] = {nbins,nbins,20}; Double_t xminEff1[3] = {0.398,pMin,-1.0}; Double_t xmaxEff1[3] = {0.598,pMax,1.0}; // gral efficiency
- Int_t binsEff2[4] = {nbins,nbins,20,10}; Double_t xminEff2[4] = {0.398,pMin,-1.0,-10.}; Double_t xmaxEff2[4] = {0.598,pMax,1.0,10.}; // vtx cut
+ Int_t binsEff2[4] = {nbins,nbins,20,20}; Double_t xminEff2[4] = {0.398,pMin,-1.0,-10.}; Double_t xmaxEff2[4] = {0.598,pMax,1.0,10.}; // vtx cut
Int_t binsEff3[4] = {nbins,nbins,20,60}; Double_t xminEff3[4] = {0.398,pMin,-1.0,0.}; Double_t xmaxEff3[4] = {0.598,pMax,1.0,1.2}; // dca between daughters
Int_t binsEff4[4] = {nbins,nbins,20,50}; Double_t xminEff4[4] = {0.398,pMin,-1.0,0.9975}; Double_t xmaxEff4[4] = {0.598,pMax,1.0,1.}; // cpa
Int_t binsEff5[5] = {nbins,nbins,20,99,99}; Double_t xminEff5[5] = {0.398,pMin,-1.0,0.,0.}; Double_t xmaxEff5[5] = {0.598,pMax,1.0,3.3,3.3}; // dca to prim. vtx
Int_t binsEff6[5] = {nbins,nbins,20,170,170}; Double_t xminEff6[5] = {0.398,pMin,-1.0,0.5,0.5}; Double_t xmaxEff6[5] = {0.598,pMax,1.0,170.5,170}; // No. TPC Cls
+ Int_t binsEffKsh[5] = {nbins,nbins,20,50,50}; Double_t xminEffKsh[5] = {0.398,pMin,-1.0,0.,0.}; Double_t xmaxEffKsh[5] = {0.598,pMax,1.0,1.0,1.0}; //shared TPC cls
for(Int_t i=0; i<kNCent; i++){
fK0sAssocMassPtDaugNClsTPC[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(fK0sAssocMassPtDaugNClsTPC[i]);
+ snprintf(hNameHist,100, "fK0sAssocMassPtShTPCcls_Cent_%d",i);
+ fK0sAssocMassPtShTPCcls[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(fK0sAssocMassPtShTPCcls[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);
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]);
+
}
fK0sMCResEta = new TH3F("fK0sMCResEta","K^{0}_{S} 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.);
}
// Histogramas para estudios sistematicos de la eficiencia
- Int_t binsEff7[3] = {nbins,nbins,20}; Double_t xminEff7[3] = {1.065,pMin,-1.0}; Double_t xmaxEff7[3] = {1.165,pMax,1.0}; // gral efficiency
- Int_t binsEff8[4] = {nbins,nbins,20,10}; Double_t xminEff8[4] = {1.065,pMin,-1.0,-10.}; Double_t xmaxEff8[4] = {1.165,pMax,1.0,10.}; // vtx
- Int_t binsEff9[4] = {nbins,nbins,20,60}; Double_t xminEff9[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEff9[4] = {1.165,pMax,1.0,1.2}; // dca between daughters
- Int_t binsEff10[4] = {nbins,nbins,20,50}; Double_t xminEff10[4] = {1.065,pMin,-1.0,0.9975}; Double_t xmaxEff10[4] = {1.165,pMax,1.0,1.}; // cpa
- 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 binsEff7[3] = {nbins,nbins,20}; Double_t xminEff7[3] = {1.065,pMin,-1.0}; Double_t xmaxEff7[3] = {1.165,pMax,1.0}; // gral efficiency
+ Int_t binsEff8[4] = {nbins,nbins,20,20}; Double_t xminEff8[4] = {1.065,pMin,-1.0,-10.}; Double_t xmaxEff8[4] = {1.165,pMax,1.0,10.}; // vtx
+ Int_t binsEff9[4] = {nbins,nbins,20,60}; Double_t xminEff9[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEff9[4] = {1.165,pMax,1.0,1.2}; // dca between daughters
+ Int_t binsEff10[4] = {nbins,nbins,20,50}; Double_t xminEff10[4] = {1.065,pMin,-1.0,0.9975}; Double_t xmaxEff10[4] = {1.165,pMax,1.0,1.}; // cpa
+ 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
for(Int_t i=0; i<kNCent; i++){
fLambdaAssocMassPtDaugNClsTPC[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos # TPC Cls; Neg # TPC Cls;",5,binsEff12,xminEff12,xmaxEff12);
fOutput->Add(fLambdaAssocMassPtDaugNClsTPC[i]);
+ snprintf(hNameHist,100, "fLambdaAssocMassPtShTPCcls_Cent_%d",i);
+ fLambdaAssocMassPtShTPCcls[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(fLambdaAssocMassPtShTPCcls[i]);
+
+
// ------------ Embeded particles
snprintf(hNameHist,100, "fLambdaAssocMassPtRapEmbeded_Cent_%d",i);
fLambdaAssocMassPtRapEmbeded[i] = new THnSparseD(hNameHist, "#Lambda Embeded; Mass (GeV/c^{2}); #it{p}_{T}; rap;",3,binsEff7,xminEff7,xmaxEff7);
snprintf(hNameHist,100, "fLambdaAssocMassPtDaugNClsTPCEmbeded_Cent_%d",i);
fLambdaAssocMassPtDaugNClsTPCEmbeded[i] = new THnSparseD(hNameHist, "#Lambda; Mass (GeV/c^{2}); #it{p}_{T}; rap; Pos # TPC Cls; Neg # TPC Cls;",5,binsEff12,xminEff12,xmaxEff12);
fOutput->Add(fLambdaAssocMassPtDaugNClsTPCEmbeded[i]);
+
+ snprintf(hNameHist,100, "fLambdaAssocMassPtShTPCclsEmbeded_Cent_%d",i);
+ 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]);
+
}
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.);
// Histogramas para estudios sistematicos de la eficiencia
Int_t binsEff13[3] = {nbins,nbins,20}; Double_t xminEff13[3] = {1.065,pMin,-1.0}; Double_t xmaxEff13[3] = {1.165,pMax,1.0}; // gral efficiency
- Int_t binsEff14[4] = {nbins,nbins,20,10}; Double_t xminEff14[4] = {1.065,pMin,-1.0,-10.}; Double_t xmaxEff14[4] = {1.165,pMax,1.0,10.}; // vtx
+ Int_t binsEff14[4] = {nbins,nbins,20,20}; Double_t xminEff14[4] = {1.065,pMin,-1.0,-10.}; Double_t xmaxEff14[4] = {1.165,pMax,1.0,10.}; // vtx
Int_t binsEff15[4] = {nbins,nbins,20,60}; Double_t xminEff15[4] = {1.065,pMin,-1.0,0.}; Double_t xmaxEff15[4] = {1.165,pMax,1.0,1.2}; // dca between daug
Int_t binsEff16[4] = {nbins,nbins,20,50}; Double_t xminEff16[4] = {1.065,pMin,-1.0,0.9975}; Double_t xmaxEff16[4] = {1.165,pMax,1.0,1.}; // cpa
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
+
for(Int_t i=0; i<kNCent; i++){
// --------- Natural particles
snprintf(hNameHist,100, "fAntiLambdaAssocMassPtRap_Cent_%d",i);
fAntiLambdaAssocMassPtDaugNClsTPC[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(fAntiLambdaAssocMassPtDaugNClsTPC[i]);
+ snprintf(hNameHist,100, "fAntiLambdaAssocMassPtShTPCcls_Cent_%d",i);
+ fAntiLambdaAssocMassPtShTPCcls[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(fAntiLambdaAssocMassPtShTPCcls[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);
snprintf(hNameHist,100, "fAntiLambdaAssocMassPtDaugNClsTPCEmbeded_Cent_%d",i);
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);
+ 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]);
+
+ }
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.);
fOutput->Add(fAntiLambdaMCResEta);
fK0sMassEmbeded->GetZaxis()->SetTitle("centrality");
fOutput->Add(fK0sMassEmbeded);
+
fK0sMassPtEta =
new TH3F("fK0sMassPtEta","K^{0}_{s}: Mass vs #it{p}_{T} vs #eta;Mass (GeV/C^{2});#it{p}_{T} (GeV/#it{c});#eta",
nbins,0.398,0.598,nbins,pMin,pMax,20,-1.0,1.0);
nbins,0.398,0.598,nbins,pMin,pMax,nbinsPhi,0.,2.*TMath::Pi());
fOutput->Add(fK0sMassPtPhi);
+ // ----- Fraction of shared TPC cls
+ fK0sPosDaugFracShTPCcls =
+ new TH3F("fK0sPosDaugFracShTPCcls","K^{0}_{s}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,0.398,0.598,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fK0sPosDaugFracShTPCcls);
+
+ fK0sNegDaugFracShTPCcls =
+ new TH3F("fK0sNegDaugFracShTPCcls","K^{0}_{s}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,0.398,0.598,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fK0sNegDaugFracShTPCcls);
+
+
// ================== Correlations =================
// ------------------------ Splitting:
}
+ // ----- Fraction of shared TPC cls
+ fK0sPosDaugFracShTPCclsTrig =
+ new TH3F("fK0sPosDaugFracShTPCclsTrig","K^{0}_{s}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,0.398,0.598,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fK0sPosDaugFracShTPCclsTrig);
+
+ fK0sNegDaugFracShTPCclsTrig =
+ new TH3F("fK0sNegDaugFracShTPCclsTrig","K^{0}_{s}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,0.398,0.598,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fK0sNegDaugFracShTPCclsTrig);
+
+
// DCA to prim vertex
fK0sDCADaugToPrimVtx
= new TH3F("fK0sDCADaugToPrimVtx","K^{0}_{S} Bckg: dca daughter vs. #it{p}_{T,l};DCA Pos daug (cm);DCA Neg daug (cm);#it{p}_{T,l} (GeV/#it{c})",
nbins,1.065,1.165,nbins,pMin,pMax,nbinsPhi,0.,2.*TMath::Pi());
fOutput->Add(fLambdaMassPtPhi);
+
+ // ----- Fraction of shared TPC cls
+ fLambdaPosDaugFracShTPCcls =
+ new TH3F("fLambdaPosDaugFracShTPCcls","#Lambda: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fLambdaPosDaugFracShTPCcls);
+
+ fLambdaNegDaugFracShTPCcls =
+ new TH3F("fLambdaNegDaugFracShTPCcls","#Lambda: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fLambdaNegDaugFracShTPCcls);
+
+
+
// ================== Correlations =================
// ----------------Splitting:
}
+ // ----- Fraction of shared TPC cls
+ fLambdaPosDaugFracShTPCclsTrig =
+ new TH3F("fLambdaPosDaugFracShTPCclsTrig","#Lambda: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fLambdaPosDaugFracShTPCclsTrig);
+
+ fLambdaNegDaugFracShTPCclsTrig =
+ new TH3F("fLambdaNegDaugFracShTPCclsTrig","#Lambda: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fLambdaNegDaugFracShTPCclsTrig);
+
+
// DCA to prim vertex
fLambdaDCADaugToPrimVtx
= new TH3F("fLambdaDCADaugToPrimVtx","#Lambda Bckg: dca daughter vs. #it{p}_{T,l}",
nbins,1.065,1.165,nbins,pMin,pMax,nbinsPhi,0.,2.*TMath::Pi());
fOutput->Add(fAntiLambdaMassPtPhi);
+ fAntiLambdaPosDaugFracShTPCcls =
+ new TH3F("fAntiLambdaPosDaugFracShTPCcls","#bar{#Lambda}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fAntiLambdaPosDaugFracShTPCcls);
+
+ fAntiLambdaNegDaugFracShTPCcls =
+ new TH3F("fAntiLambdaNegDaugFracShTPCcls","#bar{#Lambda}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fAntiLambdaNegDaugFracShTPCcls);
// ================== Correlations =================
}
+ // ----- Fraction of shared TPC cls
+ fAntiLambdaPosDaugFracShTPCclsTrig =
+ new TH3F("fAntiLambdaPosDaugFracShTPCclsTrig","#bar{#Lambda}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fAntiLambdaPosDaugFracShTPCclsTrig);
+
+ fAntiLambdaNegDaugFracShTPCclsTrig =
+ new TH3F("fAntiLambdaNegDaugFracShTPCclsTrig","#bar{#Lambda}: mass vs #it{p}_{T} vs fraction Shared TPC cls;Mass (GeV/c^2);#it{p}_{T} (GeV/#it{c});fraction Shared TPC cls",
+ nbins,1.065,1.165,nbins,pMin,pMax,50,0,1.);
+ fOutput->Add(fAntiLambdaNegDaugFracShTPCclsTrig);
+
// DCA to prim vertex
fAntiLambdaDCADaugToPrimVtx
for(Int_t j=0;j<kNVtxZ;j++){
snprintf(hNameHist,100,"fK0sdPhidEtaME_%.2f_%.2f_%.0f_%.0f_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[ll],kBinCent[ll+1],j);
- fK0sdPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH2F(hNameHist,"K^{0}_{S}: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
+ fK0sdPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH3F(hNameHist,"K^{0}_{S}: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
- nbinsdEta,-1.5,1.5);
+ nbinsdEta,-1.5,1.5,nbins,0.398,0.598);
fOutputME->Add(fK0sdPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j]);
}
}
for(Int_t j=0;j<kNVtxZ;j++){
snprintf(hNameHist,100,"fLambdadPhidEtaME_%.2f_%.2f_%.0lf_%.0lf_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[ll],kBinCent[ll+1],j);
- fLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH2F(hNameHist,"#Lambda: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
+ fLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH3F(hNameHist,"#Lambda: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
- nbinsdEta,-1.5,1.5);
+ nbinsdEta,-1.5,1.5,nbins,1.065,1.165);
fOutputME->Add(fLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j]);
}
}
for(Int_t j=0;j<kNVtxZ;j++){
snprintf(hNameHist,100,"fAntiLambdadPhidEtaME_%.2f_%.2f_%.0lf_%.0lf_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[ll],kBinCent[ll+1],j);
- fAntiLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH2F(hNameHist,"#bar{#Lambda}: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
+ fAntiLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j] = new TH3F(hNameHist,"#bar{#Lambda}: #Delta#varphi vs #Delta#eta in ME;#Delta#varphi (rad);#Delta#eta",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
- nbinsdEta,-1.5,1.5);
+ nbinsdEta,-1.5,1.5,nbins,1.065,1.165);
fOutputME->Add(fAntiLambdadPhidEtaME[ll*kN1*kNVtxZ + k*kNVtxZ + j]);
}
}
Float_t lPhi = v0->Phi();
//lPhi = ( (lPhi < 0) ? lPhi + 2*TMath::Pi() : lPhi );
+ // Fraction of TPC Shared Cluster
+ Float_t fracPosDaugTPCSharedMap = GetFractionTPCSharedCls(ptrack);
+ Float_t fracNegDaugTPCSharedMap = GetFractionTPCSharedCls(ntrack);
+
// rapidity
Float_t rapK0s = v0->RapK0Short();
Float_t rapLambda = v0->RapLambda();
Double_t effK0sVtx[4] = {v0->MassK0Short(),ptAs,rapAs,zv};
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};
+
// Distributions for the efficiency (systematics chechks)
fK0sAssocPtMassArm[curCentBin]->Fill(effK0sArm);
fK0sAssocMassPtVtx[curCentBin]->Fill(effK0sVtx);
fK0sAssocMassPtDCADaug[curCentBin]->Fill(effK0sDCA);
fK0sAssocMassPtCPA[curCentBin]->Fill(effK0sCPA);
+ fK0sAssocMassPtShTPCcls[curCentBin]->Fill(effK0sShTPCcls);
}
fK0sMCResEta->Fill(resEta,pt,centrality);
Double_t effK0sVtx[4] = {v0->MassK0Short(),ptAs,rapAs,zv};
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};
+
// Distributions for the efficiency (systematics chechks)
fK0sAssocPtMassArmEmbeded[curCentBin]->Fill(effK0sArm);
fK0sAssocMassPtVtxEmbeded[curCentBin]->Fill(effK0sVtx);
fK0sAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effK0sDCA);
fK0sAssocMassPtCPAEmbeded[curCentBin]->Fill(effK0sCPA);
+ fK0sAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effK0sShTPCcls);
}
} // End selection in the dca to prim. vtx and the number of clusters
Double_t effLambdaVtx[4] = {v0->MassLambda(),ptAs,rapAs,zv};
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};
// Distributions for the efficiency (systematics chechks)
fLambdaAssocMassPtRap[curCentBin]->Fill(effLambda);
fLambdaAssocMassPtVtx[curCentBin]->Fill(effLambdaVtx);
fLambdaAssocMassPtDCADaug[curCentBin]->Fill(effLambdaDCA);
fLambdaAssocMassPtCPA[curCentBin]->Fill(effLambdaCPA);
+ fLambdaAssocMassPtShTPCcls[curCentBin]->Fill(effLambdaShTPCcls);
if( !isCandidate2K0s && !isCandidate2LambdaBar)
fLambdaAssocMassPtRap2[curCentBin]->Fill(effLambda);
Double_t effLambdaVtx[4] = {v0->MassLambda(),ptAs,rapAs,zv};
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};
// Distributions for the efficiency (systematics chechks)
fLambdaAssocMassPtRapEmbeded[curCentBin]->Fill(effLambda);
fLambdaAssocMassPtVtxEmbeded[curCentBin]->Fill(effLambdaVtx);
fLambdaAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effLambdaDCA);
fLambdaAssocMassPtCPAEmbeded[curCentBin]->Fill(effLambdaCPA);
+ fLambdaAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effLambdaShTPCcls);
if( !isCandidate2K0s && !isCandidate2LambdaBar)
fLambdaAssocMassPtRapEmbeded2[curCentBin]->Fill(effLambda);
Double_t effAntiLambdaVtx[4] = {v0->MassAntiLambda(),ptAs,rapAs,zv};
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};
// Distributions for the efficiency (systematics chechks)
fAntiLambdaAssocMassPtRap[curCentBin]->Fill(effAntiLambda);
fAntiLambdaAssocMassPtVtx[curCentBin]->Fill(effAntiLambdaVtx);
fAntiLambdaAssocMassPtDCADaug[curCentBin]->Fill(effAntiLambdaDCA);
fAntiLambdaAssocMassPtCPA[curCentBin]->Fill(effAntiLambdaCPA);
-
+ fAntiLambdaAssocMassPtShTPCcls[curCentBin]->Fill(effAntiLambdaShTPCcls);
+
if( !isCandidate2K0s && !isCandidate2Lambda )
fAntiLambdaAssocMassPtRap2[curCentBin]->Fill(effAntiLambda);
}
Double_t effAntiLambdaVtx[4] = {v0->MassAntiLambda(),ptAs,rapAs,zv};
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};
// Distributions for the efficiency (systematics chechks)
fAntiLambdaAssocMassPtRapEmbeded[curCentBin]->Fill(effAntiLambda);
fAntiLambdaAssocMassPtVtxEmbeded[curCentBin]->Fill(effAntiLambdaVtx);
fAntiLambdaAssocMassPtDCADaugEmbeded[curCentBin]->Fill(effAntiLambdaDCA);
fAntiLambdaAssocMassPtCPAEmbeded[curCentBin]->Fill(effAntiLambdaCPA);
+ fAntiLambdaAssocMassPtShTPCclsEmbeded[curCentBin]->Fill(effAntiLambdaShTPCcls);
if( !isCandidate2K0s && !isCandidate2Lambda )
fAntiLambdaAssocMassPtRapEmbeded2[curCentBin]->Fill(effAntiLambda);
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 3) );
+ fK0sPosDaugFracShTPCcls->Fill(massK0s,pt,fracPosDaugTPCSharedMap);
+ fK0sNegDaugFracShTPCcls->Fill(massK0s,pt,fracNegDaugTPCSharedMap);
}
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 4) );
+ fLambdaPosDaugFracShTPCcls->Fill(massLambda,pt,fracPosDaugTPCSharedMap);
+ fLambdaNegDaugFracShTPCcls->Fill(massLambda,pt,fracNegDaugTPCSharedMap);
+
}
// Invariant Mass cut
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 5) );
-
+
+ fAntiLambdaPosDaugFracShTPCcls->Fill(massAntiLambda,pt,fracPosDaugTPCSharedMap);
+ fAntiLambdaNegDaugFracShTPCcls->Fill(massAntiLambda,pt,fracNegDaugTPCSharedMap);
+
}
if( fDoQA && lCheckMcAntiLambda && isNaturalPart && (pt<10.) ){ // Quality Assurance
// ---------------- Fraction of TPC Shared Cluster:
fracTrigTPCSharedMap = GetFractionTPCSharedCls(tTrig);
-
+ fTrigFracShTPCcls->Fill(ptTrig,fracTrigTPCSharedMap);
+
for(Int_t j=0; j<fAssocParticles->GetEntriesFast(); j++){
AliMiniParticle* trackAssocME = (AliMiniParticle*) (fAssocParticles->At(j));
AliAODv0 *tAssoc=fAOD->GetV0(trackAssocME->ID());
if( trackAssocME->WhichCandidate() == 4 ){
splitCont2[1] = massL;
// Negative daughter
- fLambdaNegDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ fLambdaNegDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
}
// ---- AntiLambda
if( trackAssocME->WhichCandidate() == 5 ){
// +++++++++++ Pt bin & centrality
fK0sdPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massK0s,weight);
+ fK0sPosDaugFracShTPCclsTrig->Fill(massK0s,pt,fracPosDaugTPCSharedMap);
+ fK0sNegDaugFracShTPCclsTrig->Fill(massK0s,pt,fracNegDaugTPCSharedMap);
+
// ==== Correlations K0s invariant mass peak ==== //
if (TMath::Abs(mK0s-massK0s) < 3*sK0s) {
// ==== Correlations Lambda invariant mass peak ==== //
// +++++++++++ Pt bin & centrality
fLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massL,weight);
+
+ fLambdaPosDaugFracShTPCclsTrig->Fill(massL,pt,fracPosDaugTPCSharedMap);
+ fLambdaNegDaugFracShTPCclsTrig->Fill(massL,pt,fracNegDaugTPCSharedMap);
// ==== Correlations Lambda invariant mass peak ==== //
if (TMath::Abs(mL-massL) < 3*sL) {
// +++++++++++ Pt bin & centrality
fAntiLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massAL,weight);
+ fAntiLambdaPosDaugFracShTPCclsTrig->Fill(massAL,pt,fracPosDaugTPCSharedMap);
+ fAntiLambdaNegDaugFracShTPCclsTrig->Fill(massAL,pt,fracNegDaugTPCSharedMap);
+
// ==== Correlations AntiLambda invariant mass peak ==== //
if (TMath::Abs(mL-massAL) < 3*sL) {
pt = tAssoc->Pt();
+ massK0s = tAssoc->MassK0Short();
+ massL = tAssoc->MassLambda();
+ massAL = tAssoc->MassAntiLambda();
+
+ /*
Bool_t IsSelected = kFALSE;
// K0s
if( trackAssocME->WhichCandidate() == 3 ){
if (TMath::Abs(mL-massAL) < 3*sL) IsSelected = kTRUE;
}
- if(!IsSelected) continue;
+ if(!IsSelected) continue;*/
phiAssocME = trackAssocME->Phi();
etaAssocME = trackAssocME->Eta();
if(binPtv0==-1) continue;
if( trackAssocME->WhichCandidate() == 3 ) {
- fK0sdPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);}
+ fK0sdPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta,massK0s);}
else if( trackAssocME->WhichCandidate() == 4 )
- fLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
+ fLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta,massL);
else if( trackAssocME->WhichCandidate() == 5 )
- fAntiLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
+ fAntiLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta,massAL);
} // End loop over V0's