#include "AliInputEventHandler.h"
#include "AliAnalysisManager.h"
+#include "AliExternalTrackParam.h"
+
#include "AliAnalysisTaskLambdaOverK0sJets.h"
//extern TROOT *gROOT;
static Float_t lMin = 0.0; // Limits in the histo for fidutial volume
static Float_t lMax = 100.; // Limits in the fidutial volume
-static Int_t nMaxEvMix = 250;
+//static Int_t nMaxEvMix = 250;
//
//
//
AliAnalysisTaskLambdaOverK0sJets::AliAnalysisTaskLambdaOverK0sJets(const char *name) :
- AliAnalysisTaskSE(name),
+AliAnalysisTaskSE(name),
fAOD(0), fCollision("PbPb2010"), fIsMC(kFALSE), fUsePID(kFALSE), fCentMin(0.), fCentMax(90.), fDoQA(kFALSE), fDoMixEvt(kFALSE), 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),
+ fTPCRadius(125.), fDiffTrigDaugFracTPCSharedCls(0.06),
+
fOutput(0), fOutputQA(0), fOutputME(0), fMEList(0x0), fTriggerParticles(0x0), fTriggerPartMC(0x0), fAssocParticles(0x0), fAssocPartMC(0x0), fEvents(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), 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), 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),
+ fInjectedParticles(0),
-fK0sMCPt(0), fK0sMCPtRap(0), fK0sMCPtRap2(0), fK0sMCPtRapEmbeded(0), fK0sAssocPt(0), fK0sAssocPtArm(0), fK0sAssocPtRap(0), fK0sAssocPtRapEmbeded(0), fK0sMCResEta(0), fK0sMCResPhi(0), fLambdaMCPt(0), fLambdaMCPtRap(0), fLambdaMCPtRap2(0), fLambdaMCPtRapEmbeded(0), fLambdaMCFromXi(0), fLambdaAssocPt(0), fLambdaAssocPtRap(0), fLambdaAssocFromXi(0), fLambdaMCResEta(0), fLambdaMCResPhi(0), fAntiLambdaMCPt(0), fAntiLambdaMCPtRap(0), fAntiLambdaMCPtRap2(0), fAntiLambdaMCPtRapEmbeded(0), fAntiLambdaMCFromXi(0), fAntiLambdaAssocPt(0), fAntiLambdaAssocPtRap(0), fAntiLambdaAssocFromXi(0), fAntiLambdaMCResEta(0), fAntiLambdaMCResPhi(0),
+ fK0sMCPt(0), fK0sMCPtRap(0), fK0sMCPtRap2(0), fK0sMCPtRapEmbeded(0), fK0sAssocPt(0), fK0sAssocPtArm(0), fK0sAssocPtRap(0), fK0sAssocPtRapEmbeded(0), fK0sMCResEta(0), fK0sMCResPhi(0), fK0sMCResPt(0), fK0sPosMCResEta(0), fK0sPosMCResPhi(0), fK0sPosMCResPt(0), fK0sNegMCResEta(0), fK0sNegMCResPhi(0), fK0sNegMCResPt(0), fLambdaMCPt(0), fLambdaMCPtRap(0), fLambdaMCPtRap2(0), fLambdaMCPtRapEmbeded(0), fLambdaMCFromXi(0), fLambdaAssocPt(0), fLambdaAssocPtRap(0), fLambdaAssocFromXi(0), fLambdaMCResEta(0), fLambdaMCResPhi(0), fLambdaMCResPt(0), fLambdaPosMCResEta(0), fLambdaPosMCResPhi(0), fLambdaPosMCResPt(0), fLambdaNegMCResEta(0), fLambdaNegMCResPhi(0), fLambdaNegMCResPt(0), fAntiLambdaMCPt(0), fAntiLambdaMCPtRap(0), fAntiLambdaMCPtRap2(0), fAntiLambdaMCPtRapEmbeded(0), fAntiLambdaMCFromXi(0), fAntiLambdaAssocPt(0), fAntiLambdaAssocPtRap(0), fAntiLambdaAssocFromXi(0), fAntiLambdaMCResEta(0), fAntiLambdaMCResPhi(0), fAntiLambdaMCResPt(0), fAntiLambdaPosMCResEta(0), fAntiLambdaPosMCResPhi(0), fAntiLambdaPosMCResPt(0), fAntiLambdaNegMCResEta(0), fAntiLambdaNegMCResPhi(0), fAntiLambdaNegMCResPt(0),
fHistArmenterosPodolanski(0), fHistArmPodBckg(0),
- fK0sMass(0), fK0sMassEmbeded(0), fK0sMassPtEta(0), fK0sMassPtPhi(0), fK0sDaughtersPt(0), fSharedClsTrigDaug(0), fK0sPosDaugSplCheckCovMat(0x0), fK0sNegDaugSplCheckCovMat(0x0), 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), fK0sDaughtersPt(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), fLambdaPosDaugSplCheckCovMat(0x0), fLambdaNegDaugSplCheckCovMat(0x0), 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), fLambdaDaughtersPt(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), fAntiLambdaPosDaugSplCheckCovMat(0x0), fAntiLambdaNegDaugSplCheckCovMat(0x0), 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), fAntiLambdaDaughtersPt(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),
- fLambdaPtPosDaug(0), fLambdaPtNegDaug(0), fLambdaBckgPtPosDaug(0), fLambdaBckgPtNegDaug(0), fLambdaPhiEtaPosDaug(0),fLambdaPhiEtaNegDaug(0), fLambdaBckgPhiEtaPosDaug(0),fLambdaBckgPhiEtaNegDaug(0), fLambdaDCAPosDaug(0),fLambdaDCANegDaug(0), fLambdaBckgDCAPosDaug(0), fLambdaBckgDCANegDaug(0), fLambdaDecayPos(0), fLambdaBckgDecayPos(0), fLambdaDecayVertex(0), fLambdaBckgDecayVertex(0), fLambdaCPA(0), fLambdaBckgCPA(0), fLambdaDCAV0Daug(0), fLambdaBckgDCAV0Daug(0), fLambdaNClustersTPC(0), fLambdaBckgNClustersTPC(0), fLambdaNClustersITSPos(0), fLambdaNClustersITSNeg(0), fLambdaBckgNClustersITSPos(0), fLambdaBckgNClustersITSNeg(0),
+ fLambdaPtPosDaug(0), fLambdaPtNegDaug(0), fLambdaBckgPtPosDaug(0), fLambdaBckgPtNegDaug(0), fLambdaPhiEtaPosDaug(0),fLambdaPhiEtaNegDaug(0), fLambdaBckgPhiEtaPosDaug(0), fLambdaBckgPhiEtaNegDaug(0), fLambdaDCAPosDaug(0),fLambdaDCANegDaug(0), fLambdaBckgDCAPosDaug(0), fLambdaBckgDCANegDaug(0), fLambdaDecayPos(0), fLambdaBckgDecayPos(0), fLambdaDecayVertex(0), fLambdaBckgDecayVertex(0), fLambdaCPA(0), fLambdaBckgCPA(0), fLambdaDCAV0Daug(0), fLambdaBckgDCAV0Daug(0), fLambdaNClustersTPC(0), fLambdaBckgNClustersTPC(0), fLambdaNClustersITSPos(0), fLambdaNClustersITSNeg(0), fLambdaBckgNClustersITSPos(0), fLambdaBckgNClustersITSNeg(0),
fAntiLambdaPtPosDaug(0), fAntiLambdaPtNegDaug(0), fAntiLambdaBckgPtPosDaug(0), fAntiLambdaBckgPtNegDaug(0), fAntiLambdaPhiEtaPosDaug(0),fAntiLambdaPhiEtaNegDaug(0), fAntiLambdaBckgPhiEtaPosDaug(0),fAntiLambdaBckgPhiEtaNegDaug(0), fAntiLambdaDCAPosDaug(0),fAntiLambdaDCANegDaug(0), fAntiLambdaBckgDCAPosDaug(0), fAntiLambdaBckgDCANegDaug(0), fAntiLambdaDecayPos(0), fAntiLambdaBckgDecayPos(0), fAntiLambdaDecayVertex(0), fAntiLambdaBckgDecayVertex(0), fAntiLambdaCPA(0), fAntiLambdaBckgCPA(0), fAntiLambdaDCAV0Daug(0), fAntiLambdaBckgDCAV0Daug(0), fAntiLambdaNClustersTPC(0), fAntiLambdaBckgNClustersTPC(0), fAntiLambdaNClustersITSPos(0), fAntiLambdaNClustersITSNeg(0), fAntiLambdaBckgNClustersITSPos(0), fAntiLambdaBckgNClustersITSNeg(0)
{
// Dummy Constructor
+ // variables for track splitting:
+ // shifted positionf for thw tracks
+ for(Int_t i=0; i<3; i++){
+ fTrigSftR125[i] = -9999.;
+ fDaugSftR125[i] = -9999.;
+ }
+
// Particles properties in MC
for (Int_t i=0; i<kNCent; i++){
fK0sAssocMassPtDaugNClsTPCEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fK0sMassPtRap[i] = 0;
+ // -- Splitting checks
+ fK0sPosDaugSplCheckCovMat[i] = 0x0;
+ fK0sNegDaugSplCheckCovMat[i] = 0x0;
+ fK0sPosDaugdPhiSdEtaS[i] = 0x0;
+ fK0sNegDaugdPhiSdEtaS[i] = 0x0;
+ fK0sPosMCResdEtaSdPhiS[i] = 0x0;
+ fK0sNegMCResdEtaSdPhiS[i] = 0x0;
// Lambda
fLambdaMCPtRapVtx[i] = 0;
fLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fLambdaMassPtRap[i] = 0;
+ // -- Splitting checks
+ fLambdaPosDaugSplCheckCovMat[i] = 0x0;
+ fLambdaNegDaugSplCheckCovMat[i] =0x0;
+ fLambdaPosDaugdPhiSdEtaS[i] = 0x0;
+ fLambdaNegDaugdPhiSdEtaS[i] = 0x0;
+ fLambdaPosMCResdEtaSdPhiS[i] = 0x0;
+ fLambdaNegMCResdEtaSdPhiS[i] = 0x0;
// AntiLambda
fAntiLambdaMCPtRapVtx[i] = 0;
fAntiLambdaAssocMassPtDaugNClsTPCEmbeded[i] = 0;
// -- Mass vs rapidity vs pt vs centrlaity
fAntiLambdaMassPtRap[i] = 0;
+ // -- Splitting checks
+ fAntiLambdaPosDaugSplCheckCovMat[i] = 0x0;
+ fAntiLambdaNegDaugSplCheckCovMat[i] = 0x0;
+ fAntiLambdaPosDaugdPhiSdEtaS[i] = 0x0;
+ fAntiLambdaNegDaugdPhiSdEtaS[i] = 0x0;
+ fAntiLambdaPosMCResdEtaSdPhiS[i] = 0x0;
+ fAntiLambdaNegMCResdEtaSdPhiS[i] = 0x0;
+
}
// Correlations in MC
fTriggerMCPtCent = new TH2F("fTriggerMCPtCent","Trigger particle MC;#it{p}_{T} (GeV/#it{c});centrality (%)",2*nbinPtLP,pMin,2*ptMaxLP,100,0.,100.);
fOutput->Add(fTriggerMCPtCent);
- fTriggerMCResPt = new TH3F("fTriggerMCResPt","Trigger particle MC: p_{t} resolution;(p_{t,MC}-p_{t,Rec})/p_{t,Rec};p_{t} (GeV/#it{c});centrality",60,-0.3,0.3,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
+ fTriggerMCResPt = new TH3F("fTriggerMCResPt","Trigger particle MC: #it{p}_{T} resolution;(p_{t,MC}-p_{t,Rec})/p_{t,Rec};#it{p}_{T} (GeV/#it{c});centrality",60,-0.3,0.3,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
fOutput->Add(fTriggerMCResPt);
- fTriggerMCResEta = new TH3F("fTriggerMCResEta","Trigger particle MC: #eta resolution; #eta_{MC}-#eta_{Rec};p_{t} (GeV/#it{c}); centrality",40,-0.1,0.1,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
+ fTriggerMCResEta = new TH3F("fTriggerMCResEta","Trigger particle MC: #eta resolution; #eta_{MC}-#eta_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
fOutput->Add(fTriggerMCResEta);
- fTriggerMCResPhi = new TH3F("fTriggerMCResPhi","Trigger particle MC: #varphi resolution; #varphi_{MC}-#varphi_{Rec};p_{t} (GeV/#it{c}); centrality",40,-0.1,0.1,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
+ fTriggerMCResPhi = new TH3F("fTriggerMCResPhi","Trigger particle MC: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,2*nbinPtLP,pMin,ptMaxLP,100,0.,100.);
fOutput->Add(fTriggerMCResPhi);
}
fTriggerIsV0 = new TH1F("fTriggerIsV0","V0 candidate is a LP;#it{p}_{T} (GeV/#it{c});Counts",nbinPtLP,pMin,ptMaxLP);
fOutput->Add(fTriggerIsV0);
- // ------------------- > Comaring properties of this trigger with the daughters
+ // ------------------- > Comparing properties of this trigger with the daughters
// K0s
fCheckIDTrigPtK0s = new TH3F("fCheckIDTrigPtK0s","K^{0}_{S};#deltap/p_{tri};;p_{V0}",40,-0.2,0.2,7,-0.5,6.5,100,1.,11.);
fCheckIDTrigPtK0s->GetYaxis()->SetBinLabel(1,"Pos Daug X");
fK0sMCResPhi = new TH3F("fK0sMCResPhi","K^{0}_{S} Assoc: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
fOutput->Add(fK0sMCResPhi);
+ fK0sMCResPt = new TH3F("fK0sMCResPt","K^{0}_{S} Assoc: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fK0sMCResPt);
+
+ fK0sPosMCResEta = new TH3F("fK0sPosMCResEta","K^{0}_{S} Pos. Daug.: #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(fK0sPosMCResEta);
+
+ fK0sPosMCResPhi = new TH3F("fK0sPosMCResPhi","K^{0}_{S} Pos. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fK0sPosMCResPhi);
+
+ fK0sPosMCResPt = new TH3F("fK0sPosMCResPt","K^{0}_{S} Pos. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fK0sPosMCResPt);
+
+ fK0sNegMCResEta = new TH3F("fK0sNegMCResEta","K^{0}_{S} Neg. Daug.: #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(fK0sNegMCResEta);
+
+ fK0sNegMCResPhi = new TH3F("fK0sNegMCResPhi","K^{0}_{S} Neg. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fK0sNegMCResPhi);
+
+ fK0sNegMCResPt = new TH3F("fK0sNegMCResPt","K^{0}_{S} Neg. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fK0sNegMCResPt);
+
// Lambda MC-true:
fLambdaMCPt = new TH1F("fLambdaMCPt","#Lambda MC;#it{p}_{T} (GeV/#it{c});Counts",nbins,pMin,pMax);
fOutput->Add(fLambdaMCPt);
fLambdaMCResPhi = new TH3F("fLambdaMCResPhi","#Lambda Assoc: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
fOutput->Add(fLambdaMCResPhi);
+ fLambdaMCResPt = new TH3F("fLambdaMCResPt","#Lambda Assoc: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fLambdaMCResPt);
+
+ fLambdaPosMCResEta = new TH3F("fLambdaPosMCResEta","#Lambda Pos. Daug.: #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(fLambdaPosMCResEta);
+
+ fLambdaPosMCResPhi = new TH3F("fLambdaPosMCResPhi","#Lambda Pos. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fLambdaPosMCResPhi);
+
+ fLambdaPosMCResPt = new TH3F("fLambdaPosMCResPt","#Lambda Pos. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fLambdaPosMCResPt);
+
+ fLambdaNegMCResEta = new TH3F("fLambdaNegMCResEta","#Lambda Neg. Daug.: #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(fLambdaNegMCResEta);
+
+ fLambdaNegMCResPhi = new TH3F("fLambdaNegMCResPhi","#Lambda Neg. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fLambdaNegMCResPhi);
+
+ fLambdaNegMCResPt = new TH3F("fLambdaNegMCResPt","#Lambda Neg. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fLambdaNegMCResPt);
+
// AntiLambda MC-true:
fAntiLambdaMCPt = new TH1F("fAntiLambdaMCPt","#bar{#Lambda} MC;#it{p}_{T} (GeV/#it{c});Counts",nbins,pMin,pMax);
fOutput->Add(fAntiLambdaMCPt);
fAntiLambdaMCResPhi = new TH3F("fAntiLambdaMCResPhi","#bar{#Lambda} Assoc: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
fOutput->Add(fAntiLambdaMCResPhi);
+ fAntiLambdaMCResPt = new TH3F("fAntiLambdaMCResPt","#bar{#Lambda} Assoc: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fAntiLambdaMCResPt);
+
+ fAntiLambdaPosMCResEta = new TH3F("fAntiLambdaPosMCResEta","#bar{#Lambda} Pos. Daug.: #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(fAntiLambdaPosMCResEta);
+
+ fAntiLambdaPosMCResPhi = new TH3F("fAntiLambdaPosMCResPhi","#bar{#Lambda} Pos. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fAntiLambdaPosMCResPhi);
+
+ fAntiLambdaPosMCResPt = new TH3F("fAntiLambdaPosMCResPt","#bar{#Lambda} Pos. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fAntiLambdaPosMCResPt);
+
+ fAntiLambdaNegMCResEta = new TH3F("fAntiLambdaNegMCResEta","#bar{#Lambda} Neg. Daug.: #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(fAntiLambdaNegMCResEta);
+
+ fAntiLambdaNegMCResPhi = new TH3F("fAntiLambdaNegMCResPhi","#bar{#Lambda} Neg. Daug.: #varphi resolution; #varphi_{MC}-#varphi_{Rec};#it{p}_{T} (GeV/#it{c}); centrality",40,-0.1,0.1,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fAntiLambdaNegMCResPhi);
+
+ fAntiLambdaNegMCResPt = new TH3F("fAntiLambdaNegMCResPt","#bar{#Lambda} Neg. Daug.: pt resolution; #it{p}_{T,MC}-#it{p]_{T,Rec};#it{p}_{T} (GeV/#it{c}); centrality",60,-0.3,0.3,nbins,pMin,pMax,100,0.,100.);
+ fOutput->Add(fAntiLambdaNegMCResPt);
+
} //End MC
// ======================================================== //
// ================== Correlations =================
- // ----------------Splitting:
- // Shared TPC clusters status between trigger and daughter tracks
- fSharedClsTrigDaug = new TH3F("fSharedClsTrigDaug","Splitting chech with TPC map;",100,0.,10.,100,0.,10.,6,0.,6.);
- fOutput->Add(fSharedClsTrigDaug);
-
- //-----
- Double_t binsDev[120];
+ // ------------------------ Splitting:
+ Double_t binsDev[121];
binsDev[0] = 0;
- for (Int_t k=0;k<=11;k++)
+ for (Int_t k=-7;k<=4;k++)
for (Int_t j=1;j<=10;j++)
- binsDev[k*10+j] = j*TMath::Power(10,k);
+ binsDev[(k+7)*10+j] = j*TMath::Power(10,k);
- Int_t binsSplit[4] = {100,100,120,7}; Double_t xminSplit[4] = {pMin,pMin,0,0.}; Double_t xmaxSplit[4] = {pMax,pMax,1e+12,7.0};
+ Int_t binsSplit[9] = {100,nbins,100,2,301,101,101,120,9}; Double_t xminSplit[9] = {pMin,0.398,pMin,-0.5,-0.001,-0.005,-0.005,0,-0.5}; Double_t xmaxSplit[9] = {pMax,0.598,pMax,1.5,0.3,1.005,1.005,10e+4,8.5};
+
+ Int_t binsSplit2[12] = {100,nbins,100,2,10,20,101,101,100,120,99,2};
+ Double_t xminSplit2[12] = {pMin,0.398,pMin,-0.5,0.,-0.1,-0.005,-0.005,-1.,0,0.,-0.5};
+ Double_t xmaxSplit2[12] = {pMax,0.598,pMax,1.5,0.1,0.1,1.005,1.005,1.,10e+4,3.3,1.5};
+
+ Int_t binsSplit3[6] = {100,nbins,100,46,46,2};
+ Double_t xminSplit3[6] = {pMin,0.398,pMin,-0.16,-0.16,-0.5};
+ Double_t xmaxSplit3[6] = {pMax,0.598,pMax,0.16,0.16,1.5};
+
+ for(Int_t j=0; j<kNCent; j++){
+
+ // positive daughter
+ fK0sPosDaugSplCheckCovMat[j] = new THnSparseD( Form("fK0sPosDaugSplCheckCovMat_%d",j), "K^{0}_{S} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fK0sPosDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fK0sPosDaugSplCheckCovMat[j]);
+
+ // negative daughter
+ fK0sNegDaugSplCheckCovMat[j] = new THnSparseD( Form("fK0sNegDaugSplCheckCovMat_%d",j), "K^{0}_{S} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fK0sNegDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fK0sNegDaugSplCheckCovMat[j]);
- // positive daughter
- fK0sPosDaugSplCheckCovMat = new THnSparseD("fK0sPosDaugSplCheckCovMat", "K^{0}_{S} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fK0sPosDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fK0sPosDaugSplCheckCovMat);
+ // Positive daughter:
+ fK0sPosDaugdPhiSdEtaS[j] = new THnSparseD(Form("fK0sPosDaugdPhiSdEtaS_%d",j), "K^{0}_{S} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; #sum_{x,y,z}(#it{p}_{i,Trig}-#it{p}_{i,Daug})^{2}/( #sigma_{i,Trig}^{2} + #sigma_{i,Daug}^{2} ); DCA to prim. vtx; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fK0sPosDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fK0sPosDaugdPhiSdEtaS[j]);
+
+ // Negative daughter:
+ fK0sNegDaugdPhiSdEtaS[j] = new THnSparseD(Form("fK0sNegDaugdPhiSdEtaS_%d",j), "K^{0}_{S} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; #sum_{x,y,z}(#it{p}_{i,Trig}-#it{p}_{i,Daug})^{2}/( #sigma_{i,Trig}^{2} + #sigma_{i,Daug}^{2} ); DCA to prim. vtx; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fK0sNegDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fK0sNegDaugdPhiSdEtaS[j]);
+
+ if(fIsMC){
+ // Positive daughter:
+ fK0sPosMCResdEtaSdPhiS[j] = new THnSparseD(Form("fK0sPosMCResdEtaSdPhiS_%d",j), "K^{0}_{S} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fK0sPosMCResdEtaSdPhiS[j]);
+
+ // Negative daughter:
+ fK0sNegMCResdEtaSdPhiS[j] = new THnSparseD(Form("fK0sNegMCResdEtaSdPhiS_%d",j), "K^{0}_{S} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fK0sNegMCResdEtaSdPhiS[j]);
+ }
+
+ }
- // negative daughter
- fK0sNegDaugSplCheckCovMat = new THnSparseD("fK0sNegDaugSplCheckCovMat", "K^{0}_{S} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fK0sNegDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fK0sNegDaugSplCheckCovMat);
// DCA to prim vertex
fK0sDCADaugToPrimVtx
- = new TH3F("fK0sDCADaugToPrimVtx","K^{0}_{S} Bckg: dca daughter vs. p_{T,l};DCA Pos daug (cm);DCA Neg daug (cm);p_{T,l} (GeV/#it{c})",
+ = 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})",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fOutput->Add(fK0sDCADaugToPrimVtx);
// Monte-Carlo level:
if(fIsMC){
snprintf(hNameHist,100, "fK0sdPhidEtaMC_%.2f_%.2f_Cent_%.0f_%.0f",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1]);
- fK0sdPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"K^{0}_{S} MC: #Delta#varphi vs #Delta#eta vs p_{T,l}",
+ fK0sdPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"K^{0}_{S} MC: #Delta#varphi vs #Delta#eta vs #it{p}_{T,l}",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
nbinsVtx,-10.,10.);
for(Int_t ll=0;ll<kNVtxZ;ll++){
snprintf(hNameHist,100, "fK0sdPhidEtaPtL_%.2f_%.2f_Cent_%.0f_%.0f_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1],ll);
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"K^{0}_{S}: #Delta#varphi vs #Delta#eta vs Inv. Mass",
- nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
- nbinsdEta,-1.5,1.5,
- nbins,0.398,0.598);
+ nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
+ nbinsdEta,-1.5,1.5,
+ nbins,0.398,0.598);
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetXaxis()->SetTitle("#Delta#varphi (rad)");
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetYaxis()->SetTitle("#Delta#eta");
fK0sdPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll]->GetZaxis()->SetTitle("Inv. Mass");
// Correlations (side-band):
fK0sBckgDecLength
- = new TH2F("fK0sBckgDecLength","K^{0}_{S} Bckg: c#tau vs. p_{T,l}",
+ = new TH2F("fK0sBckgDecLength","K^{0}_{S} Bckg: c#tau vs. #it{p}_{T,l}",
100,0.,15.,nbinPtLP,pMin,ptMaxLP);
fK0sBckgDecLength->GetXaxis()->SetTitle("c#tau (cm)");
- fK0sBckgDecLength->GetYaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fK0sBckgDecLength->GetYaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fK0sBckgDecLength);
fK0sBckgDCADaugToPrimVtx
- = new TH3F("fK0sBckgDCADaugToPrimVtx","K^{0}_{S} Bckg: dca daughter vs. p_{T,l}",
+ = new TH3F("fK0sBckgDCADaugToPrimVtx","K^{0}_{S} Bckg: dca daughter vs. #it{p}_{T,l}",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fK0sBckgDCADaugToPrimVtx->GetXaxis()->SetTitle("DCA Pos daug (cm)");
fK0sBckgDCADaugToPrimVtx->GetYaxis()->SetTitle("DCA Neg daug (cm)");
- fK0sBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fK0sBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fK0sBckgDCADaugToPrimVtx);
fK0sBckgEtaPhi =
// ================== Correlations =================
// ----------------Splitting:
+ xminSplit[1] = 1.065; xmaxSplit[1] = 1.165;
+ xminSplit2[1] = 1.065; xmaxSplit2[1] = 1.165;
+ xminSplit3[1] = 1.065; xmaxSplit3[1] = 1.165;
+
+ for(Int_t j=0; j<kNCent; j++){
+
+ // positive daughter:
+ fLambdaPosDaugSplCheckCovMat[j] = new THnSparseD( Form("fLambdaPosDaugSplCheckCovMat_%d",j), "#Lambda Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fLambdaPosDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fLambdaPosDaugSplCheckCovMat[j]);
+
+ // negative daughter:
+ fLambdaNegDaugSplCheckCovMat[j] = new THnSparseD( Form("fLambdaNegDaugSplCheckCovMat_%d",j), "#Lambda Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fLambdaNegDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fLambdaNegDaugSplCheckCovMat[j]);
+
+ // Positive daughter:
+ fLambdaPosDaugdPhiSdEtaS[j] = new THnSparseD(Form("fLambdaPosDaugdPhiSdEtaS_%d",j), "#Lambda Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; #sum_{x,y,z}(#it{p}_{i,Trig}-#it{p}_{i,Daug})^{2}/( #sigma_{i,Trig}^{2} + #sigma_{i,Daug}^{2} ); DCA to prim. vtx; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fLambdaPosDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fLambdaPosDaugdPhiSdEtaS[j]);
+
+ // Negative daughter:
+ fLambdaNegDaugdPhiSdEtaS[j] = new THnSparseD(Form("fLambdaNegDaugdPhiSdEtaS_%d",j), "#Lambda Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; #sum_{x,y,z}(#it{p}_{i,Trig}-#it{p}_{i,Daug})^{2}/( #sigma_{i,Trig}^{2} + #sigma_{i,Daug}^{2} ); DCA to prim. vtx; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fLambdaNegDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fLambdaNegDaugdPhiSdEtaS[j]);
+
+ if(fIsMC){
+ // Positive daughter:
+ fLambdaPosMCResdEtaSdPhiS[j] = new THnSparseD(Form("fLambdaPosMCResdEtaSdPhiS_%d",j), "#Lambda Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fLambdaPosMCResdEtaSdPhiS[j]);
+
+ // Negative daughter:
+ fLambdaNegMCResdEtaSdPhiS[j] = new THnSparseD(Form("fLambdaNegMCResdEtaSdPhiS_%d",j), "#Lambda Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fLambdaNegMCResdEtaSdPhiS[j]);
+ }
- // positive daughter
- fLambdaPosDaugSplCheckCovMat = new THnSparseD("fLambdaPosDaugSplCheckCovMat", "#Lambda Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fLambdaPosDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fLambdaPosDaugSplCheckCovMat);
+ }
- // negative daughter
- fLambdaNegDaugSplCheckCovMat = new THnSparseD("fLambdaNegDaugSplCheckCovMat", "#Lambda Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fLambdaNegDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fLambdaNegDaugSplCheckCovMat);
// DCA to prim vertex
fLambdaDCADaugToPrimVtx
- = new TH3F("fLambdaDCADaugToPrimVtx","#Lambda Bckg: dca daughter vs. p_{T,l}",
+ = new TH3F("fLambdaDCADaugToPrimVtx","#Lambda Bckg: dca daughter vs. #it{p}_{T,l}",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fLambdaDCADaugToPrimVtx->GetXaxis()->SetTitle("DCA Pos daug (cm)");
fLambdaDCADaugToPrimVtx->GetYaxis()->SetTitle("DCA Neg daug (cm)");
- fLambdaDCADaugToPrimVtx->GetZaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fLambdaDCADaugToPrimVtx->GetZaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fLambdaDCADaugToPrimVtx);
// Spatial Resoltuion between trigger- and asosciated- particles
// Monte-Carlo level:
if(fIsMC){
snprintf(hNameHist,100, "fLambdadPhidEtaMC_%.2f_%.2f_Cent_%.0f_%.0f",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1]);
- fLambdadPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"#Lambda MC: #Delta#varphi vs #Delta#eta vs p_{T,l}",
+ fLambdadPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"#Lambda MC: #Delta#varphi vs #Delta#eta vs #it{p}_{T,l}",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
nbinsVtx,-10.,10.);
// Reconstruction level:
for(Int_t ll=0;ll<kNVtxZ;ll++){
snprintf(hNameHist,100, "fLambdadPhidEtaPtL_%.2f_%.2f_Cent_%.0f_%.0f_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1],ll);
- fLambdadPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"#Lambda: #Delta#varphi vs #Delta#eta vs p_{T,l}",
+ fLambdadPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"#Lambda: #Delta#varphi vs #Delta#eta vs #it{p}_{T,l}",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
nbins,1.065,1.165);
// Correlations (side-band):
fLambdaBckgDecLength
- = new TH2F("fLambdaBckgDecLength","#Lambda Bckg: c#tau vs. p_{T,l}",
+ = new TH2F("fLambdaBckgDecLength","#Lambda Bckg: c#tau vs. #it{p}_{T,l}",
100,0.,25.,nbinPtLP,pMin,ptMaxLP);
fLambdaBckgDecLength->GetXaxis()->SetTitle("c#tau (cm)");
- fLambdaBckgDecLength->GetYaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fLambdaBckgDecLength->GetYaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fLambdaBckgDecLength);
fLambdaBckgDCADaugToPrimVtx
- = new TH3F("fLambdaBckgDCADaugToPrimVtx","#Lambda Bckg: dca daughter vs. p_{T,l}",
+ = new TH3F("fLambdaBckgDCADaugToPrimVtx","#Lambda Bckg: dca daughter vs. #it{p}_{T,l}",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fLambdaBckgDCADaugToPrimVtx->GetXaxis()->SetTitle("DCA Pos daug (cm)");
fLambdaBckgDCADaugToPrimVtx->GetYaxis()->SetTitle("DCA Neg daug (cm)");
- fLambdaBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fLambdaBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fLambdaBckgDCADaugToPrimVtx);
fLambdaBckgEtaPhi =
// ================== Correlations =================
// ----------------Splitting:
+ for(Int_t j=0; j<kNCent; j++){
+
+ // positive daughter:
+ fAntiLambdaPosDaugSplCheckCovMat[j] = new THnSparseD(Form("fAntiLambdaPosDaugSplCheckCovMat_%d",j), "#bar{#Lambda} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fAntiLambdaPosDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fAntiLambdaPosDaugSplCheckCovMat[j]);
+
+ // negative daughter:
+ fAntiLambdaNegDaugSplCheckCovMat[j] = new THnSparseD(Form("fAntiLambdaNegDaugSplCheckCovMat_%d",j), "#bar{#Lambda} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; R#Delta#varphi*_{max}; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; (X-X')^{2}/( #sigma^{2} + #sigma'^{2} ); Variables;",9,binsSplit,xminSplit,xmaxSplit);
+ fAntiLambdaNegDaugSplCheckCovMat[j]->SetBinEdges(7,binsDev);
+ fOutput->Add(fAntiLambdaNegDaugSplCheckCovMat[j]);
+
+ // Positive daughter:
+ fAntiLambdaPosDaugdPhiSdEtaS[j] = new THnSparseD(Form("fAntiLambdaPosDaugdPhiSdEtaS_%d",j), "#bar{#Lambda} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fAntiLambdaPosDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fAntiLambdaPosDaugdPhiSdEtaS[j]);
+
+ // Negative daughter:
+ fAntiLambdaNegDaugdPhiSdEtaS[j] = new THnSparseD(Form("fAntiLambdaNegDaugdPhiSdEtaS_%d",j), "#bar{#Lambda} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); Same Sign as Trigger Particle; #Delta#varphi*; #Delta#eta*; Trigger: fraction of TPC shared cls; Daughter: fraction of TPC shared cls; Correlation fraction of shared cls: Trigger - Daughter; same MC label;",12,binsSplit2,xminSplit2,xmaxSplit2);
+ fAntiLambdaNegDaugdPhiSdEtaS[j]->SetBinEdges(9,binsDev);
+ fOutput->Add(fAntiLambdaNegDaugdPhiSdEtaS[j]);
+
+ if(fIsMC){
+ // Positive daughter:
+ fAntiLambdaPosMCResdEtaSdPhiS[j] = new THnSparseD(Form("fAntiLambdaPosMCResdEtaSdPhiS_%d",j), "#bar{#Lambda} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fAntiLambdaPosMCResdEtaSdPhiS[j]);
+
+ // Negative daughter:
+ fAntiLambdaNegMCResdEtaSdPhiS[j] = new THnSparseD(Form("fAntiLambdaNegMCResdEtaSdPhiS_%d",j), "#bar{#Lambda} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); Mass (GeV/c^2); #it{p}_{Daug} (GeV/#it{c}); #Delta#varphi*; #Delta#eta*; Same Sign as Trigger Particle;",6,binsSplit3,xminSplit3,xmaxSplit3);
+ fOutput->Add(fAntiLambdaNegMCResdEtaSdPhiS[j]);
+ }
- // positive daughter
- fAntiLambdaPosDaugSplCheckCovMat = new THnSparseD("fAntiLambdaPosDaugSplCheckCovMat", "#bar{#Lambda} Pos. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fAntiLambdaPosDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fAntiLambdaPosDaugSplCheckCovMat);
+ }
- // negative daughter
- fAntiLambdaNegDaugSplCheckCovMat = new THnSparseD("fAntiLambdaNegDaugSplCheckCovMat", "#bar{#Lambda} Neg. daughter; #it{p}_{T,V0} (GeV/#it{c}); #it{p}_{Daug} (GeV/#it{c}); X-X'/( #sigma^{2} - #sigma'^{2} ); Variables;",4,binsSplit,xminSplit,xmaxSplit);
- fAntiLambdaNegDaugSplCheckCovMat->SetBinEdges(2,binsDev);
- fOutput->Add(fAntiLambdaNegDaugSplCheckCovMat);
+
// DCA to prim vertex
fAntiLambdaDCADaugToPrimVtx
- = new TH3F("fAntiLambdaDCADaugToPrimVtx","#bar{#Lambda} Bckg: dca daughter vs. p_{T,l}",
+ = new TH3F("fAntiLambdaDCADaugToPrimVtx","#bar{#Lambda} Bckg: dca daughter vs. #it{p}_{T,l}",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fAntiLambdaDCADaugToPrimVtx->GetXaxis()->SetTitle("DCA Pos daug (cm)");
fAntiLambdaDCADaugToPrimVtx->GetYaxis()->SetTitle("DCA Neg daug (cm)");
- fAntiLambdaDCADaugToPrimVtx->GetZaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fAntiLambdaDCADaugToPrimVtx->GetZaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fAntiLambdaDCADaugToPrimVtx);
// Spatial Resoltuion between trigger- and asosciated- particles
// Monte-Carlo level:
if(fIsMC){
snprintf(hNameHist,100, "fAntiLambdadPhidEtaMC_%.2f_%.2f_Cent_%.0f_%.0f",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1]);
- fAntiLambdadPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"#bar{#Lambda} MC: #Delta#varphi vs #Delta#eta vs p_{T,l}",
+ fAntiLambdadPhidEtaMC[jj*kN1+k] = new TH3F(hNameHist,"#bar{#Lambda} MC: #Delta#varphi vs #Delta#eta vs #it{p}_{T,l}",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
nbinsVtx,-10.,10.);
// Reconstruction level:
for(Int_t ll=0;ll<kNVtxZ;ll++){
snprintf(hNameHist,100, "fAntiLambdadPhidEtaPtL_%.2f_%.2f_Cent_%.0f_%.0f_%d",kPtBinV0[k],kPtBinV0[k+1],kBinCent[jj],kBinCent[jj+1],ll);
- fAntiLambdadPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"#bar{#Lambda}: #Delta#varphi vs #Delta#eta vs p_{T,l}",
+ fAntiLambdadPhidEtaPtL[jj*kN1*kNVtxZ + k*kNVtxZ + ll] = new TH3F(hNameHist,"#bar{#Lambda}: #Delta#varphi vs #Delta#eta vs #it{p}_{T,l}",
nbinsdPhi,-TMath::PiOver2(),3*TMath::PiOver2(),
nbinsdEta,-1.5,1.5,
nbins,1.065,1.165);
// Correlations (side-band):
fAntiLambdaBckgDecLength
- = new TH2F("fAntiLambdaBckgDecLength","#bar{#Lambda} Bckg: c#tau vs. p_{T,l}",
+ = new TH2F("fAntiLambdaBckgDecLength","#bar{#Lambda} Bckg: c#tau vs. #it{p}_{T,l}",
100,0.,25.,nbinPtLP,pMin,ptMaxLP);
fAntiLambdaBckgDecLength->GetXaxis()->SetTitle("c#tau (cm)");
- fAntiLambdaBckgDecLength->GetYaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fAntiLambdaBckgDecLength->GetYaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fAntiLambdaBckgDecLength);
fAntiLambdaBckgDCADaugToPrimVtx
- = new TH3F("fAntiLambdaBckgDCADaugToPrimVtx","#bar{#Lambda} Bckg: dca daughter vs. p_{T,l}",
+ = new TH3F("fAntiLambdaBckgDCADaugToPrimVtx","#bar{#Lambda} Bckg: dca daughter vs. #it{p}_{T,l}",
90,0.,3.3,90,0.,3.3,nbinPtLP,pMin,ptMaxLP);
fAntiLambdaBckgDCADaugToPrimVtx->GetXaxis()->SetTitle("DCA Pos daug (cm)");
fAntiLambdaBckgDCADaugToPrimVtx->GetYaxis()->SetTitle("DCA Neg daug (cm)");
- fAntiLambdaBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("p_{T,l} (GeV/#it{c})");
+ fAntiLambdaBckgDCADaugToPrimVtx->GetZaxis()->SetTitle("#it{p}_{T,l} (GeV/#it{c})");
fOutput->Add(fAntiLambdaBckgDCADaugToPrimVtx);
fAntiLambdaBckgEtaPhi =
//___________________________________________________________________________________________
-Bool_t AliAnalysisTaskLambdaOverK0sJets::AcceptTrack(AliAODTrack *t)
+Bool_t AliAnalysisTaskLambdaOverK0sJets::AcceptTrack(const AliAODTrack *t)
{
// Track criteria for primaries particles
-
- // TPC only tracks
+
if (TMath::Abs(t->Eta())>0.8 ) return kFALSE;
- if (!(t->TestFilterMask(1<<7))) return kFALSE;
+ // TPC only tracks
+ //if (!(t->TestFilterMask(1<<7))) return kFALSE;
+ // Hybrid tracks
+ if( !(t->TestFilterBit(272)) ) return kFALSE;
Float_t nCrossedRowsTPC = t->GetTPCClusterInfo(2,1);
if (nCrossedRowsTPC < 70) return kFALSE;
-
- // Hybrid tracks
- //if( !(t->TestFilterBit(272)) ) return kFALSE;
- //if(!t->IsHybridGlobalConstrainedGlobal()) return kFALSE;
-
- // Point in the SPD
+ // Point in the SPD
Int_t SPDHits = t->HasPointOnITSLayer(0) + t->HasPointOnITSLayer(1);
if( SPDHits )
fTriggerWiSPDHit->Fill(1.5);
//___________________________________________________________________________________________
-static Int_t SameTrack(AliAODTrack *trk, const AliAODTrack *daugTrk)
+static Int_t SameTrack(const AliAODTrack *trig, const AliAODTrack *daug)
{
// Local method to compaire the momentum between two tracks
isSamePt = 1;
*/
- if( (TMath::Abs(daugTrk->GetID())+1)==(TMath::Abs(trk->GetID())) )
+ if( (TMath::Abs(daug->GetID())+1)==(TMath::Abs(trig->GetID())) )
isSamePt = 1;
//___________________________________________________________________________________________
+static Int_t SameLabel(const AliAODTrack *trig, const AliAODTrack *daug)
+{
+ // Compaire the the label value that points back to the Monte Carlo production
+ //cout << trig->GetLabel() << " " << daug->GetLabel() << endl;
+
+ if( TMath::Abs(trig->GetLabel() ) ==
+ TMath::Abs(daug->GetLabel() ) )
+ return 1.0;
+
+ return 0.;
+
+}
+
+//___________________________________________________________________________________________
+
static Float_t SpatialResolution(Float_t p1x,Float_t p1y,Float_t p2x,Float_t p2y,Float_t dist)
{
// Obtains the spacial resolution between trigger and V0
//___________________________________________________________________________________________
-Bool_t AliAnalysisTaskLambdaOverK0sJets::GoodTPCSharedMap(const AliAODTrack *track){
- // Rejects tracks with shared clusters after filling a control histogram
- // This overload is used for primaries
-
- // Get the shared maps
- const TBits sharedMap = track->GetTPCSharedMap();
- // Fill a control histogram
- //fPriHistShare->Fill(sharedMap.CountBits());
- // Reject shared clusters
- if((sharedMap.CountBits()) >= 1){
- // Bad track, has too many shared clusters!
- return kFALSE;
- }
- return kTRUE;
+static Float_t GetDPhiStar(Float_t phi1, Float_t pt1, Float_t charge1, Float_t phi2, Float_t pt2, Float_t charge2, Float_t radius, Float_t bSign)
+{
+ //
+ // calculates dphistar
+ //
+
+ Float_t dphistar = phi1 - phi2 - charge1 * bSign * TMath::ASin(0.075 * radius / pt1) + charge2 * bSign * TMath::ASin(0.075 * radius / pt2);
+ static const Double_t kPi = TMath::Pi();
+
+ // circularity
+ if (dphistar > kPi)
+ dphistar = kPi * 2 - dphistar;
+ if (dphistar < -kPi)
+ dphistar = -kPi * 2 - dphistar;
+ if (dphistar > kPi) // might look funny but is needed
+ dphistar = kPi * 2 - dphistar;
+
+ return dphistar;
+
}
+
//___________________________________________________________________________________________
-Bool_t AliAnalysisTaskLambdaOverK0sJets::GoodTPCSharedMap(AliAODTrack *track){
+static Float_t TwoTrackEfficiencyCut(Float_t phi1, Float_t eta1, Float_t pt1, Float_t charge1, Float_t phi2, Float_t eta2, Float_t pt2, Float_t charge2, Float_t bSign)
+{
+ // Code taken from the HBT analysis to reject the track splitting
+ // It was modified to provide only the value of kDphiStarMax
+ // and a possible rejection in the kDphiStarMean
+
+ Float_t kRadiousDphiStarMax = -0.0005;
+ Float_t deta = eta1 - eta2;
+ Float_t twoTrackEfficiencyCutValue = 0.02;
+
+ // optimization
+ if (TMath::Abs(deta) < twoTrackEfficiencyCutValue * 2.5 * 3) {
+
+ // check first boundaries to see if is worth to loop and find the minimum
+ Float_t dphistar1 = GetDPhiStar(phi1, pt1, charge1, phi2, pt2, charge2, 0.8, bSign);
+ Float_t dphistar2 = GetDPhiStar(phi1, pt1, charge1, phi2, pt2, charge2, 2.5, bSign);
+
+ const Float_t kLimit = twoTrackEfficiencyCutValue * 3;
+
+ //Float_t dphistarminabs = 1e5;
+ //Float_t dphistarmin = 1e5;
+
+ if (TMath::Abs(dphistar1) < kLimit || TMath::Abs(dphistar2) < kLimit || dphistar1 * dphistar2 < 0){
+
+ kRadiousDphiStarMax = 0;
+ //kDphiStarMean = 0;
+ //Int_t i=0;
+
+ for (Double_t rad=0.8; rad<2.51; rad+=0.01){
+
+ if ( TMath::Abs(0.075 * rad / pt2)>1 ) break;
+
+ Float_t dphistar = GetDPhiStar(phi1, pt1, charge1, phi2, pt2, charge2, rad, bSign);
+ Float_t dphistarabs = TMath::Abs(dphistar);
+
+ if( ( (dphistarabs*rad) > kRadiousDphiStarMax) && ( TMath::Abs(deta) < twoTrackEfficiencyCutValue ) ){
+ kRadiousDphiStarMax = dphistarabs*rad;
+ }
+
+ //kDphiStarMean += dphistarabs;
+ //i++;
+
+ }
+
+ //kDphiStarMean = kDphiStarMean/i;
+ /*if (TMath::Abs(deta) < twoTrackEfficiencyCutValue && kDphiStarMean < twoTrackEfficiencyCutValue ){
+ return kFALSE;
+ }*/
+
+ } // End selection in dphistar
+
+ } // End dEta value
+
+
+ return kRadiousDphiStarMax;
+
+}
+
+//___________________________________________________________________________________________
+/*
+static Bool_t GoodTPCSharedMap(const AliAODTrack *track)
+{
// Rejects tracks with shared clusters after filling a control histogram
// This overload is used for primaries
}
return kTRUE;
}
-
+*/
//___________________________________________________________________________________________
-/* under developing
-Float_t AliAnalysisTaskLambdaOverK0sJets::GetFractionTPCSharedCls(AliAODTrack *track1,const AliAODTrack *track2){
+
+static Float_t GetFractionTPCSharedCls(const AliAODTrack *track)
+{
// Rejects tracks with shared clusters after filling a control histogram
// This overload is used for primaries
// Get the shared maps
- const TBits fTPCClusterMap1 = track1->GetTPCClusterMap();
- const TBits fTPCClusterMap2 = track2->GetTPCClusterMap();
- // Fill a control histogram
- //fPriHistShare->Fill(sharedMap.CountBits());
- // Reject shared clusters
+ const TBits sharedMap = track->GetTPCSharedMap();
+
+ return 1.*sharedMap.CountBits()/track->GetTPCNclsF();
- for(Int_t i=0; i<=159; i++)
- cout << fTPCClusterMap1.CountBits(i) << " " << fTPCClusterMap2.CountBits(i) << endl;
-
- return 1.;
+}
+
+//___________________________________________________________________________________________
+
+Double_t AliAnalysisTaskLambdaOverK0sJets::ThetaS(TString part)
+{
+ // LINES OBTAINED FROM THE FEMTOSCOPY ANALYSIS:
+ // private communication with Hans Beck
+
+ // Returns the longitudinal angle of the particles propagated
+ // position at R=1.25m. See
+ // https://edms.cern.ch/file/406391/2/ALICE-INT-2003-038.pdf
+ // for the ALICE coordinate system. Theta is zero at positive z,
+ // pi/2 at z = 0 aka the xy plane and pi at negative z
+
+ // R^ ^
+ // | /
+ // |θ'/
+ // | / θ
+ // |/____>z
+ //
+ // Let's compute θ' and θ = π/2 - θ'
+ // where θ' can even be and should
+ // sometimes be negative
+ // tan(θ') = z/R
+ // θ' = arctan(z/R)
+ // θ = π/2 - θ'
+ // = π/2 - arctan(z/R)
+ // Note that in the doc above theta
+ // is calculated as arccos(z/sqrt(x^2+y^2+z^2))
+
+ // Array of positions is 85,105,125,..cm,
+ // we take the z position at R=1.25m
+ // return TMath::Pi()/2. - TMath::ATan(fXshifted[2][2]/125.);
+ /*
+ if( part.EqualTo("Trigger") )
+ return TMath::Pi()/2. - TMath::ATan(fTrigSftR125[2]/125.);
+ else if( part.EqualTo("Daughter") )
+ return TMath::Pi()/2. - TMath::ATan(fDaugSftR125[2]/125.);
+ */
+ Double_t thetaS = -100.;
+
+ if( part.EqualTo("Trigger") )
+ thetaS = TMath::Pi()/2. - TMath::ATan(fTrigSftR125[2]/fTPCRadius);
+ if( part.EqualTo("Daughter") )
+ thetaS = TMath::Pi()/2. - TMath::ATan(fDaugSftR125[2]/fTPCRadius);
+
+ return thetaS;
+
}
-*/
+
+//___________________________________________________________________________________________
+
+Double_t AliAnalysisTaskLambdaOverK0sJets::EtaS(TString part)
+{
+ // LINES OBTAINED FROM THE FEMTOSCOPY ANALYSIS:
+ // private communication with Hans Beck
+
+ // Returns the corresponding eta of a pri. part.
+ // with this particles pos at R=1.25m
+
+ // http://en.wikipedia.org/wiki/Pseudorapidity
+ // η = -ln[ tan(θ/2)]
+ // printf("z: %+04.0f, thetaS %+03.2f etaS %+1.2f\n"
+ // ,fXshifted[2][2],ThetaS(),-TMath::Log( TMath::Tan(ThetaS()/2.) ));
+
+ return -TMath::Log( TMath::Tan(ThetaS(part)/2.) );
+}
+
+//___________________________________________________________________________________________
+
+Float_t AliAnalysisTaskLambdaOverK0sJets::dEtaS()
+{
+ // LINES OBTAINED FROM THE FEMTOSCOPY ANALYSIS:
+ // private communication with Hans Beck
+
+ // Returns the pseudorapidity star difference
+
+ // It is important to keep the calculations easy and separated.
+ // The calculation of EtaS is straight forward, one just has to
+ // do it step by step to not get confused.
+ return EtaS("Trigger") - EtaS("Daughter");
+}
+
+//___________________________________________________________________________________________
+
+Float_t AliAnalysisTaskLambdaOverK0sJets::dPhiSAtR125()
+{
+ // LINES OBTAINED FROM THE FEMTOSCOPY ANALYSIS:
+ // private communication with Hans Beck
+
+ // returns delta phi star at R=1.2m
+ // position at R=1.2m is stored as second radius
+ // const Float_t distSft= TMath::Sqrt(TMath::Power(track1.fXshifted[2][0] - track2.fXshifted[2][0],2)
+ // +TMath::Power(track1.fXshifted[2][1] - track2.fXshifted[2][1],2));
+ const Float_t distSft= TMath::Sqrt( TMath::Power(fTrigSftR125[0] - fDaugSftR125[0],2) +
+ TMath::Power(fTrigSftR125[1] - fDaugSftR125[1],2) );
+ //return 2.0 * TMath::ATan(distSft/2./(125.));
+ return 2.0 * TMath::ATan(distSft/2./(fTPCRadius));
+}
+
+
//___________________________________________________________________________________________
-void AliAnalysisTaskLambdaOverK0sJets::RecCascade(AliAODTrack *trk1,const AliAODTrack *trk2,const AliAODTrack *trkBch,TString histo)
+void AliAnalysisTaskLambdaOverK0sJets::SetSftPosR125(const AliAODTrack *track,const Float_t bfield,const Float_t priVtx[3], TString part)
+{
+ // LINES OBTAINED FROM THE FEMTOSCOPY ANALYSIS:
+ // private communication with Hans Beck
+
+ // Sets the spatial position of the track at the radius R=1.25m in the shifted coordinate system
+
+ // Initialize the array to something indicating there was no propagation
+ if(part.EqualTo("Trigger")){
+ fTrigSftR125[0] = -9999.;
+ fTrigSftR125[1] = -9999.;
+ fTrigSftR125[2] = -9999.;
+ }
+ if(part.EqualTo("Daughter")){
+ fDaugSftR125[0] = -9999.;
+ fDaugSftR125[1] = -9999.;
+ fDaugSftR125[2] = -9999.;
+ }
+
+ // Make a copy of the track to not change parameters of the track
+ AliExternalTrackParam etp;
+ etp.CopyFromVTrack(track);
+
+ // The global position of the the track
+ Double_t xyz[3]={-9999.,-9999.,-9999.};
+
+ // The radius we want to propagate to, squared
+ //const Float_t RSquaredWanted(125.*125.);
+ const Float_t RSquaredWanted(fTPCRadius*fTPCRadius);
+
+ // Propagation is done in local x of the track
+ for (Float_t x = 58.; x < 247.; x+=1.){
+ // Starts at 83 / Sqrt(2) and goes outwards. 85/Sqrt(2) is the smallest local x
+ // for global radius 85 cm. x = 245 is the outer radial limit of the TPC when
+ // the track is straight, i.e. has inifinite pt and doesn't get bent.
+ // If the track's momentum is smaller than infinite, it will develop a y-component,
+ // which adds to the global radius
+ // We don't change the propagation steps to not mess up things!
+
+ // Stop if the propagation was not succesful. This can happen for low pt tracks
+ // that don't reach outer radii
+ if(!etp.PropagateTo(x,bfield)) break;
+ etp.GetXYZ(xyz); // GetXYZ returns global coordinates
+
+ // Calculate the shifted radius we are at, squared.
+ // Compare squared radii for faster code
+ Float_t shiftedRadiusSquared = (xyz[0]-priVtx[0])*(xyz[0]-priVtx[0])
+ + (xyz[1]-priVtx[1])*(xyz[1]-priVtx[1]);
+
+ // Roughly reached the radius we want
+ if(shiftedRadiusSquared > RSquaredWanted){
+
+ // Bigger loop has bad precision, we're nearly one centimeter too far,
+ // go back in small steps.
+ while (shiftedRadiusSquared>RSquaredWanted){
+ // Propagate a mm inwards
+ x-=.1;
+ if(!etp.PropagateTo(x,bfield)){
+ // Propagation failed but we're already with a
+ // cm precision at R=1.25m so we only break the
+ // inner loop
+ break;
+ }
+ // Get the global position
+ etp.GetXYZ(xyz);
+ // Calculate shifted radius, squared
+ shiftedRadiusSquared = (xyz[0]-priVtx[0])*(xyz[0]-priVtx[0])
+ + (xyz[1]-priVtx[1])*(xyz[1]-priVtx[1]);
+ }
+
+ // We reached R=1.25m with a precission of a cm to a mm,
+ // set the spatial position
+ if(part.EqualTo("Trigger")){
+ fTrigSftR125[0] = xyz[0] - priVtx[0];
+ fTrigSftR125[1] = xyz[1] - priVtx[1];
+ fTrigSftR125[2] = xyz[2] - priVtx[2];
+
+ /*cout << endl
+ << xyz[0] << " " << xyz[1] << " " << xyz[2] << endl;
+ cout << fTrigSftR125[0] << " " << fTrigSftR125[1] << " " <<fTrigSftR125[2] << endl;*/
+ }
+ if(part.EqualTo("Daughter")){
+ fDaugSftR125[0] = xyz[0] - priVtx[0];
+ fDaugSftR125[1] = xyz[1] - priVtx[1];
+ fDaugSftR125[2] = xyz[2] - priVtx[2];
+
+ /*cout << endl
+ << xyz[0] << " " << xyz[1] << " " << xyz[2] << endl
+ << fDaugSftR125[0] << " " << fDaugSftR125[1] << " " <<fDaugSftR125[2] << endl;*/
+ }
+
+ // Done
+ return;
+
+ } // End of if roughly reached radius
+
+ } // End of coarse propagation loop
+
+}
+
+//___________________________________________________________________________________________
+
+void AliAnalysisTaskLambdaOverK0sJets::RecCascade(const AliAODTrack *trk1,const AliAODTrack *trk2,const AliAODTrack *trkBch,TString histo)
{
// Local method to reconstruct cascades candidates from the combinations of three tracks
// The input tracks correspond to the trigger particle and the daughter tracks of the V0 candidate (correlation step)
// 1) TriggerCheck
// 2) Reconstruction
- AliAODTrack *trkTrig = 0x0;
+ const AliAODTrack *trkTrig = 0x0;
Float_t ptTrig = -100.;
Float_t phiTrig = -100.;
Float_t etaTrig = -100.;
if( pPart->IsPhysicalPrimary() ) lMCAssocPosDaug = 1;
/*
- if ( TMath::Abs(nPart->Eta()) > fMaxEtaDaughter ||
- TMath::Abs(pPart->Eta()) > fMaxEtaDaughter )
+ if ( TMath::Abs(nPart->Eta()) > fMaxEtaDaughter ||
+ TMath::Abs(pPart->Eta()) > fMaxEtaDaughter )
goto noas;
*/
/*
(inMother==ipMother) && (lPDGCodeV0==3122) ){
if ( ( TMath::Abs(lPdgcodeMotherOfMother) == 3212) /*||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3224) ||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3214) ||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3114)*/
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3224) ||
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3214) ||
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3114)*/
) lComeFromSigma = kTRUE;
else lComeFromSigma = kFALSE;
if ( ( TMath::Abs(lPdgcodeMotherOfMother) == 3212) /*||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3224) ||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3214) ||
- ( TMath::Abs(lPdgcodeMotherOfMother) == 3114)*/
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3224) ||
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3214) ||
+ ( TMath::Abs(lPdgcodeMotherOfMother) == 3114)*/
) lComeFromSigma = kTRUE;
else lComeFromSigma = kFALSE;
rapAs = etaAs;
}
- // phi resolution for V0-reconstruction
+ // phi resolution for V0-reconstruction and daughter tracks
Float_t resEta = p0->Eta() - v0->Eta();
Float_t resPhi = p0->Phi() - v0->Phi();
+ Float_t resPt = p0->Pt() - v0->Pt();
+ Float_t resEtaPosDaug = pPart->Eta() - ptrack->Eta();
+ Float_t resPhiPosDaug = pPart->Phi() - ptrack->Phi();
+ Float_t resPtPosDaug = pPart->Pt() - ptrack->Pt();
+
+ Float_t resEtaNegDaug = nPart->Eta() - ntrack->Eta();
+ Float_t resPhiNegDaug = nPart->Phi() - ntrack->Phi();
+ Float_t resPtNegDaug = nPart->Pt() - ntrack->Pt();
+
if ( (l < 0.01) && (ptAs<10.) ) { // Primary V0
// K0s:
fK0sMCResEta->Fill(resEta,pt,centrality);
fK0sMCResPhi->Fill(resPhi,pt,centrality);
-
+ fK0sMCResPt->Fill(resPt,pt,centrality);
+
+ fK0sPosMCResEta->Fill(resEtaPosDaug,pt,centrality);
+ fK0sPosMCResPhi->Fill(resPhiPosDaug,pt,centrality);
+ fK0sPosMCResPt->Fill(resPtPosDaug,pt,centrality);
+
+ fK0sNegMCResEta->Fill(resEtaNegDaug,pt,centrality);
+ fK0sNegMCResPhi->Fill(resPhiNegDaug,pt,centrality);
+ fK0sNegMCResPt->Fill(resPtNegDaug,pt,centrality);
+
} // End selection in the dca to prim. vtx and the number of clusters
// Distributions for the efficiency (Systematic checks)
fLambdaMCResEta->Fill(resEta,pt,centrality);
fLambdaMCResPhi->Fill(resPhi,pt,centrality);
+ fLambdaMCResPt->Fill(resPt,pt,centrality);
+
+ fLambdaPosMCResEta->Fill(resEtaPosDaug,pt,centrality);
+ fLambdaPosMCResPhi->Fill(resPhiPosDaug,pt,centrality);
+ fLambdaPosMCResPt->Fill(resPtPosDaug,pt,centrality);
+
+ fLambdaNegMCResEta->Fill(resEtaNegDaug,pt,centrality);
+ fLambdaNegMCResPhi->Fill(resPhiNegDaug,pt,centrality);
+ fLambdaNegMCResPt->Fill(resPtNegDaug,pt,centrality);
} // End selection in the dca to prim. vtx and the number of clusters
fAntiLambdaMCResEta->Fill(resEta,pt,centrality);
fAntiLambdaMCResPhi->Fill(resPhi,pt,centrality);
+ fAntiLambdaMCResPt->Fill(resPt,pt,centrality);
+
+ fAntiLambdaPosMCResEta->Fill(resEtaPosDaug,pt,centrality);
+ fAntiLambdaPosMCResPhi->Fill(resPhiPosDaug,pt,centrality);
+ fAntiLambdaPosMCResPt->Fill(resPtPosDaug,pt,centrality);
+
+ fAntiLambdaNegMCResEta->Fill(resEtaNegDaug,pt,centrality);
+ fAntiLambdaNegMCResPhi->Fill(resPhiNegDaug,pt,centrality);
+ fAntiLambdaNegMCResPt->Fill(resPtNegDaug,pt,centrality);
} // End selection in the dca to prim. vtx and the number of clusters
} // End embeded particle selection
- } // End AntiLambda
+ } // End AntiLambda
// Xi decay:
if( lComeFromXi && isNaturalPart ){
if(lPDGCodeV0==3122) { fLambdaAssocFromXi->Fill(ptAs,centrality); }
if( (pt>kPtBinV0[0]) && (pt<kPtBinV0[kN1]) && isNaturalPart )
fAssocParticles->Add( new AliMiniParticle(centrality, zv, iV0, pt, lPhi, lEta, lMCAssocNegDaug, lMCAssocPosDaug, 5) );
-
}
if( fDoQA && lCheckMcAntiLambda && isNaturalPart && (pt<10.) ){ // Quality Assurance
if( (dcaPos>0.1) && (dcaNeg>0.1) ){
if( (nClsTPCPos>fDaugNClsTPC) && (nClsTPCNeg>fDaugNClsTPC) ){
- fAntiLambdaPtPosDaug->Fill(pt,lPtPos);
- fAntiLambdaPtNegDaug->Fill(pt,lPtNeg);
+ fAntiLambdaPtPosDaug->Fill(pt,lPtPos);
+ fAntiLambdaPtNegDaug->Fill(pt,lPtNeg);
- fAntiLambdaPhiEtaPosDaug->Fill(phiPos,etaPos,pt);
- fAntiLambdaPhiEtaNegDaug->Fill(phiNeg,etaNeg,pt);
+ fAntiLambdaPhiEtaPosDaug->Fill(phiPos,etaPos,pt);
+ fAntiLambdaPhiEtaNegDaug->Fill(phiNeg,etaNeg,pt);
- fAntiLambdaDecayPos->Fill(dx,dy,pt);
- fAntiLambdaDecayVertex->Fill(lt,pt);
+ fAntiLambdaDecayPos->Fill(dx,dy,pt);
+ fAntiLambdaDecayVertex->Fill(lt,pt);
- fAntiLambdaCPA->Fill(cpa,pt);
- fAntiLambdaDCAV0Daug->Fill(dca,pt);
+ fAntiLambdaCPA->Fill(cpa,pt);
+ fAntiLambdaDCAV0Daug->Fill(dca,pt);
- fAntiLambdaNClustersITSPos->Fill(phiPos,posITSNcls,pt);
- fAntiLambdaNClustersITSNeg->Fill(phiNeg,negITSNcls,pt);
- }
+ fAntiLambdaNClustersITSPos->Fill(phiPos,posITSNcls,pt);
+ fAntiLambdaNClustersITSNeg->Fill(phiNeg,negITSNcls,pt);
+ }
fAntiLambdaNClustersTPC->Fill(phiPos,nClsTPCPos,pt);
fAntiLambdaNClustersTPC->Fill(phiNeg,nClsTPCNeg,-pt);
Float_t resPhi = -1000.;
for (Int_t i=0; i<nTrk; i++) {
- AliAODTrack *t = fAOD->GetTrack(i);
+ const AliAODTrack *t = fAOD->GetTrack(i);
if(!AcceptTrack(t)) continue;
pt=t->Pt();
eta=t->Eta();
fEvents->Fill(6);
Float_t xv=vtx->GetX(), yv=vtx->GetY();
+ const Float_t priVtx[3] = {xv,yv,zv};
if (TMath::Abs(zv) > 10.) return;
fEvents->Fill(7);
eventPlane = ( (eventPlane < 0) ? eventPlane + TMath::Pi() : eventPlane );
eventPlane = ( ( eventPlane > TMath::Pi() ) ? eventPlane - TMath::Pi() : eventPlane );
+ // Magnetic field sign
+ const Float_t bSign = (fAOD->GetMagneticField() > 0) ? 1 : -1;
+
// Getting PID Response
fPIDResponse = hdr->GetPIDResponse();
if(trig->WhichCandidate()==0){
fTriggerComingFromDaug->Fill(trig->Pt());
fCheckTriggerFromV0Daug->Fill(1);
- fTriggerPtCentCh->Fill(trig->Pt(),centrality,zv);
if(fIsV0LP) fCheckTriggerFromV0Daug->Fill(2);
}
else if( trig->WhichCandidate()==1){
- fTriggerEtaPhi->Fill(trig->Phi(),trig->Eta());
- fTriggerPtCent->Fill(trig->Pt(),centrality,zv);
- fTriggerPtCentCh->Fill(trig->Pt(),centrality,zv);
fCheckTriggerFromV0Daug->Fill(0);
+ fTriggerPtCentCh->Fill(trig->Pt(),centrality,zv);
phi2 = ( (trig->Phi() > TMath::Pi()) ? trig->Phi() - TMath::Pi() : trig->Phi() ) ;
fTriggerEventPlane->Fill(phi2);
if (!pDaughter0 || !pDaughter1) continue;
/*
- if ( TMath::Abs(pDaughter0->Eta()) > fMaxEtaDaughter ||
- TMath::Abs(pDaughter1->Eta()) > fMaxEtaDaughter )
+ if ( TMath::Abs(pDaughter0->Eta()) > fMaxEtaDaughter ||
+ TMath::Abs(pDaughter1->Eta()) > fMaxEtaDaughter )
continue;
*/
// Daughter momentum cut: ! FIX it in case of AOD !
/*
- if ( ( pDaughter0->Pt() < fMinPtDaughter ) ||
- ( pDaughter1->Pt() < fMinPtDaughter ) )
- continue;
+ if ( ( pDaughter0->Pt() < fMinPtDaughter ) ||
+ ( pDaughter1->Pt() < fMinPtDaughter ) )
+ continue;
*/
if ((p0->Pt())<pMin || (p0->Pt())>100. ) continue;
// K0s
if (lPdgcodeCurrentPart == kK0Short) {
-
fK0sMCPt->Fill(lPtCurrentPart);
fK0sMCPtRap->Fill(lPtCurrentPart,lRapCurrentPart,centrality);
} // End pt selection
// Xi-
/*
- if(lPdgcodeCurrentPart == kXiMinus || lPdgcodeCurrentPart == 3322){
+ if(lPdgcodeCurrentPart == kXiMinus || lPdgcodeCurrentPart == 3322){
if( isNaturalPart )
- fAssocPartMC->Add( new AliMiniParticle(centrality, zv, iTrkMC, lPtCurrentPart, lPhiCurrentPart, lEtaCurrentPart, 0, 0, 6) );
+ fAssocPartMC->Add( new AliMiniParticle(centrality, zv, iTrkMC, lPtCurrentPart, lPhiCurrentPart, lEtaCurrentPart, 0, 0, 6) );
- } //End Xi
- */
+ } //End Xi
+ */
} // End loop over MC
} // End loop over trigger particles
} // End MC condition
-
+
// *************************************************
// V0 loop - AOD
fAssocParticles = new TObjArray();
V0Loop(kReconstruction,kTRUE,-1,-1);
else
V0Loop(kReconstruction,kFALSE,-1,-1);
-
+
//-------------------------------------------------------------
// Correlations
//-------------------------------------------------------------
+
+ // variables for correlations:
Float_t ptTrig=0., pxTrig=0., pyTrig=0.;
Float_t massK0s=0., mK0s=0., sK0s=0.;
Float_t massL=0., mL=0., sL=0.;
Float_t dx=-100., dy=-100., lt=-100., res=-100.;
Float_t dlK=-100., dlL=-100.;
Float_t dPhi=-100., dEta=-100., radio=-100.;
- Bool_t trigTPCMapOk=kTRUE, posDaugTPCMapOk=kTRUE, negDaugTPCMapOk=kTRUE;
- Float_t fracPosDaugTPCMap =0., fracNegDaugTPCMap =0., den=1.;
+
+ // variables for track splititing checks:
+ Float_t posdPhiS = -9999., posdEtaS = -9999., negdPhiS = -9999., negdEtaS = -9999.;
+ Float_t fracTrigTPCSharedMap=-1., fracPosDaugTPCSharedMap =-1., fracNegDaugTPCSharedMap =-1.;
+ //Bool_t trigTPCMapOk=kTRUE, posDaugTPCMapOk=kTRUE, negDaugTPCMapOk=kTRUE;
+ Float_t RdPhiStarMaxPosDaug=-1., RdPhiStarMaxNegDaug=-1., den=1.;
Double_t trigCov[21], posDaugCov[21], negDaugCov[21];
Double_t trigPos[6], posDaugPos[6], negDaugPos[6];
Double_t trigXYZ[3], posDaugXYZ[3], negDaugXYZ[3];
- Double_t devPosDaugTrig[7], devNegDaugTrig[7], splitCont[4];
+ Double_t devPosDaugTrig[9], devNegDaugTrig[9], splitCont[9], splitCont2[12];
+ Int_t sameSignPosDaug = -1, sameSignNegDaug = -1;
+ Float_t sameLabelPosDaug = 0., sameLabelNegDaug = 0.;
+ Int_t tlab, nlab, plab;
+ Double_t resdEtsSdPhiS[6];
// --------------------------------
// h-V0 correlations
AliMiniParticle* trig = (AliMiniParticle*) fTriggerParticles->At(i);
if( trig->WhichCandidate() == 0 ) continue;
- AliAODTrack *tTrig = (AliAODTrack*)fAOD->GetTrack(trig->ID());
+ const AliAODTrack *tTrig = (AliAODTrack*)fAOD->GetTrack(trig->ID());
ptTrig = tTrig->Pt(); pxTrig = tTrig->Px(); pyTrig = tTrig->Py();
+ // ---------------- Fraction of TPC Shared Cluster:
+ fracTrigTPCSharedMap = GetFractionTPCSharedCls(tTrig);
+
for(Int_t j=0; j<fAssocParticles->GetEntriesFast(); j++){
AliMiniParticle* trackAssocME = (AliMiniParticle*) (fAssocParticles->At(j));
AliAODv0 *tAssoc=fAOD->GetV0(trackAssocME->ID());
Int_t lMCAssocNegDaug = trackAssocME->NegDaugMCLabel();
Int_t lMCAssocPosDaug = trackAssocME->PosDaugMCLabel();
- // ------------------------------------------
- // Track splitting check
- // ------------------------------------------
+ // ----------------------------------------------------------------------------
+
+ // -----------------------------------------------------------------
+ // ****************** Track splitting check ******************
+ // -----------------------------------------------------------------
- trigTPCMapOk=kTRUE; posDaugTPCMapOk=kTRUE; negDaugTPCMapOk=kTRUE;
- if(TMath::Abs(dPhi)<0.06 && TMath::Abs(dEta)<0.06){
+ sameLabelPosDaug = 0.; sameLabelNegDaug = 0.;
+ sameSignPosDaug = -1; sameSignNegDaug = -1;
+ RdPhiStarMaxPosDaug=-1.; RdPhiStarMaxNegDaug=-1.;
+ //trigTPCMapOk=kTRUE; posDaugTPCMapOk=kTRUE; negDaugTPCMapOk=kTRUE;
+ fracPosDaugTPCSharedMap=0; fracNegDaugTPCSharedMap=0;
- // TPC Shared Map
- trigTPCMapOk = GoodTPCSharedMap(tTrig);
- posDaugTPCMapOk = GoodTPCSharedMap(ptrack);
- negDaugTPCMapOk = GoodTPCSharedMap(ntrack);
+ // ---------------- Fraction of TPC Shared Cluster
+ fracPosDaugTPCSharedMap = GetFractionTPCSharedCls(ptrack);
+ fracNegDaugTPCSharedMap = GetFractionTPCSharedCls(ntrack);
+
+ // =========== Classical methods for track-splitting ============= //
+ if( TMath::Abs(dPhi)<0.1 && TMath::Abs(dEta)<0.1 ){
+
+ // --------- Check sign of the trigger and daughter track:
+ if(tTrig->Charge()==1) { sameSignPosDaug = 1; sameSignNegDaug = 0; }
+ else { sameSignPosDaug = 0; sameSignNegDaug = 1; }
- // fraction of shared clusters
- //if(!trigTPCMapOk && !posDaugTPCMapOk) fracPosDaugTPCMap = GetFractionTPCSharedCls(tTrig,ptrack);
- //if(!trigTPCMapOk && !negDaugTPCMapOk) fracNegDaugTPCMap = GetFractionTPCSharedCls(tTrig,ntrack);
+ // -------- Shifting charged tracks to the primary vertex.
+ // -------- See HBT anlayses:
- // get position
+ // Trigger particle:
+ SetSftPosR125(tTrig,bSign,priVtx,"Trigger");
+
+ // Positive daughter: calculating delta(phi)* and delta(eta)*
+ SetSftPosR125(ptrack,bSign,priVtx,"Daughter");
+ posdPhiS = dPhiSAtR125();
+ posdEtaS = dEtaS();
+
+ // Negative daughter: calculating delta(phi)* and delta(eta)*
+ SetSftPosR125(ntrack,bSign,priVtx,"Daughter");
+ negdPhiS = dPhiSAtR125();
+ negdEtaS = dEtaS();
+
+ // ------ Get position:
tTrig->GetXYZ(trigXYZ);
ptrack->GetXYZ(posDaugXYZ);
ntrack->GetXYZ(negDaugXYZ);
- // Covaraince matrix for the tracks
+ // ------ Covaraince matrix for the tracks:
tTrig->GetCovarianceXYZPxPyPz(trigCov);
ptrack->GetCovarianceXYZPxPyPz(posDaugCov);
ntrack->GetCovarianceXYZPxPyPz(negDaugCov);
- // position and momentum
+ // ------- position and momentum:
// trigger particle
trigPos[0] = trigXYZ[0]; trigPos[1] = trigXYZ[1]; trigPos[2] = trigXYZ[2];
trigPos[3] = tTrig->Px(); trigPos[4] = tTrig->Py(); trigPos[5] = tTrig->Pz();
negDaugPos[0] = negDaugXYZ[0]; negDaugPos[1] = negDaugXYZ[1]; negDaugPos[2] = negDaugXYZ[2];
negDaugPos[3] = ntrack->Px(); negDaugPos[4] = ntrack->Py(); negDaugPos[5] = ntrack->Pz();
-
- // deviation between the two tracks
+ // ------- deviation between the two tracks:
+ // positive daughter
for(Int_t ll=0;ll<6;ll++){
- den = trigCov[ll*(ll+1)/2+ll]*trigCov[ll*(ll+1)/2+ll] + posDaugCov[ll*(ll+1)/2+ll]* posDaugCov[ll*(ll+1)/2+ll] ;
+ den = trigCov[ll*(ll+1)/2+ll] + posDaugCov[ll*(ll+1)/2+ll] ;
devPosDaugTrig[ll] = 0.;
if(den>0) devPosDaugTrig[ll] = TMath::Power( trigPos[ll] - posDaugPos[ll] ,2) / den;
- devPosDaugTrig[6] += devPosDaugTrig[ll];
+ if(ll<3) devPosDaugTrig[6] += devPosDaugTrig[ll]; // sum in X,Y,Z
+ if(ll>2) devPosDaugTrig[7] += devPosDaugTrig[ll]; // sum in momemtum
+ devPosDaugTrig[8] += devPosDaugTrig[ll]; // sum in all variables
}
-
+ // negative daughter
for(Int_t ll=0;ll<6;ll++){
- den = trigCov[ll*(ll+1)/2+ll]*trigCov[ll*(ll+1)/2+ll] + negDaugCov[ll*(ll+1)/2+ll]* negDaugCov[ll*(ll+1)/2+ll] ;
+ den = trigCov[ll*(ll+1)/2+ll] + negDaugCov[ll*(ll+1)/2+ll] ;
devNegDaugTrig[ll] = 0;
if(den>0) devNegDaugTrig[ll] = TMath::Power( trigPos[ll] - negDaugPos[ll] ,2) / den;
- devNegDaugTrig[6] += devNegDaugTrig[ll];
+ if(ll<3) devNegDaugTrig[6] += devNegDaugTrig[ll]; // sum in X,Y,Z
+ if(ll>2) devNegDaugTrig[7] += devNegDaugTrig[ll]; // sum in momemtum
+ devNegDaugTrig[8] += devNegDaugTrig[ll]; // sum in all variables
+
}
- }
+
+ // ---------------- Monte Carlo check for track-splitting
+ if(fIsMC){
+
+ TList *lst = fAOD->GetList();
+ stack = (TClonesArray*)lst->FindObject(AliAODMCParticle::StdBranchName());
+ if (!stack) {
+ Printf("ERROR: stack not available");
+ return;
+ }
+
+ sameLabelPosDaug = 1.*SameLabel(tTrig,ptrack);
+ sameLabelNegDaug = 1.*SameLabel(tTrig,ntrack);
+
+ // Resolution of delta(phi)* and delta(eta)*
+ tlab = TMath::Abs(tTrig->GetLabel());
+ plab = TMath::Abs(ptrack->GetLabel());
+ nlab = TMath::Abs(ntrack->GetLabel());
+
+ AliAODMCParticle *tPart=(AliAODMCParticle*)stack->UncheckedAt(tlab);
+ AliAODMCParticle *pPart=(AliAODMCParticle*)stack->UncheckedAt(plab);
+ AliAODMCParticle *nPart=(AliAODMCParticle*)stack->UncheckedAt(nlab);
+
+ resdEtsSdPhiS[0] = pt;
+
+ //positive daughter
+ resdEtsSdPhiS[2] = ptrack->Pt();
+ resdEtsSdPhiS[3] = (tPart->Phi() - pPart->Phi()) - posdPhiS;
+ resdEtsSdPhiS[4] = (tPart->Eta() - pPart->Eta()) - posdEtaS;
+ resdEtsSdPhiS[5] = sameSignPosDaug;
+
+ if( trackAssocME->WhichCandidate() == 3 ){
+ resdEtsSdPhiS[1] = massK0s;
+ fK0sPosMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+ if( trackAssocME->WhichCandidate() == 4 ){
+ resdEtsSdPhiS[1] = massL;
+ fLambdaPosMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+ if( trackAssocME->WhichCandidate() == 5 ){
+ resdEtsSdPhiS[1] = massAL;
+ fAntiLambdaPosMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+
+ // negative daughter
+ resdEtsSdPhiS[2] = ntrack->Pt();
+ resdEtsSdPhiS[3] = (tPart->Phi() - nPart->Phi()) - negdPhiS;
+ resdEtsSdPhiS[4] = (tPart->Eta() - nPart->Eta()) - negdEtaS;
+ resdEtsSdPhiS[5] = sameSignNegDaug;
+
+ if( trackAssocME->WhichCandidate() == 3 ){
+ resdEtsSdPhiS[1] = massK0s;
+ fK0sNegMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+ if( trackAssocME->WhichCandidate() == 4 ){
+ resdEtsSdPhiS[1] = massL;
+ fLambdaNegMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+ if( trackAssocME->WhichCandidate() == 5 ){
+ resdEtsSdPhiS[1] = massAL;
+ fAntiLambdaNegMCResdEtaSdPhiS[curCentBin]->Fill(resdEtsSdPhiS);
+ }
+
+ }
+
+ // ================ Alternative methods for track-splitting ==================
+ if(TMath::Abs(dPhi)<0.02 && TMath::Abs(dEta)<0.02){
+
+ // --------- Calculate TPCRadius*Delta(phi)Star_Max distance:
+ RdPhiStarMaxPosDaug = TwoTrackEfficiencyCut( tTrig->Phi(), tTrig->Eta(), tTrig->Pt(), tTrig->Charge(), ptrack->Phi(), ptrack->Eta(), ptrack->Pt(), 1, bSign);
+ RdPhiStarMaxNegDaug = TwoTrackEfficiencyCut( tTrig->Phi(), tTrig->Eta(), tTrig->Pt(), tTrig->Charge(), ntrack->Phi(), ntrack->Eta(), ntrack->Pt(), -1, bSign);
+
+ // -------- Comparison between trigger and daughter tracks:
+ // -------- Filling deviation of matrix elements
+ splitCont[0] = pt; splitCont[5] = fracTrigTPCSharedMap;
+
+ // ---------------------------
+ // -------- Positive daughter:
+ splitCont[2] = ptrack->Pt(); splitCont[3] = sameSignPosDaug;
+ splitCont[4] = RdPhiStarMaxPosDaug; splitCont[6] = fracPosDaugTPCSharedMap;
+
+ // ----K0s
+ if( trackAssocME->WhichCandidate() == 3 ){
+ splitCont[1] = massK0s;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devPosDaugTrig[ll]; splitCont[8] = ll;
+ fK0sPosDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // ----Lambda
+ if( trackAssocME->WhichCandidate() == 4 ){
+ splitCont[1] = massL;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devPosDaugTrig[ll]; splitCont[8] = ll;
+ fLambdaPosDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // ----AntiLambda
+ if( trackAssocME->WhichCandidate() == 5 ){
+ splitCont[1] = massAL;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devPosDaugTrig[ll]; splitCont[8] = ll;
+ fAntiLambdaPosDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // End: Positive daughter
+
+ // ---------------------------
+ // -------- Negative daughter:
+ splitCont[2] = ntrack->Pt(); splitCont[3] = sameSignNegDaug;
+ splitCont[4] = RdPhiStarMaxNegDaug; splitCont[6] = fracNegDaugTPCSharedMap;
+
+ // ----K0s
+ if( trackAssocME->WhichCandidate() == 3 ){
+ splitCont[1] = massK0s;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devNegDaugTrig[ll]; splitCont[8] = ll;
+ fK0sNegDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // ----Lambda
+ if( trackAssocME->WhichCandidate() == 4 ){
+ splitCont[1] = massL;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devNegDaugTrig[ll]; splitCont[8] = ll;
+ fLambdaNegDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // ----AntiLambda
+ if( trackAssocME->WhichCandidate() == 5 ){
+ splitCont[1] = massAL;
+ for(Int_t ll=0; ll<=8; ll++){
+ splitCont[7] = devNegDaugTrig[ll]; splitCont[8] = ll;
+ fAntiLambdaNegDaugSplCheckCovMat[curCentBin]->Fill(splitCont);
+ }
+
+ }
+ // End: Negative daughter
+
+ } // end selection in |delta(eta)| < 0.02, |delta(phi)| < 0.02
+
+
+ // ================ FILLING THnSparse: Classical track-splitting method: d(phi)* and d(eta)*
+ splitCont2[0] = pt; splitCont2[6] = fracTrigTPCSharedMap;
+
+ // --------------------------
+ // -------- Positive daughter:
+ splitCont2[2] = ptrack->Pt(); splitCont2[3] = sameSignPosDaug; splitCont2[4] = posdPhiS; splitCont2[5] = posdEtaS;
+ splitCont2[7] = fracPosDaugTPCSharedMap; splitCont2[8] = fracTrigTPCSharedMap - fracPosDaugTPCSharedMap;
+ splitCont2[9] = devPosDaugTrig[7]; splitCont2[10] = tAssoc->DcaPosToPrimVertex(); splitCont2[11] = sameLabelPosDaug;
+
+ // ---- K0s
+ if( trackAssocME->WhichCandidate() == 3 ){
+ splitCont2[1] = massK0s;
+ // Positive daughter
+ fK0sPosDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+ // ---- Lambda
+ if( trackAssocME->WhichCandidate() == 4 ){
+ splitCont2[1] = massL;
+ // Positive daughter
+ fLambdaPosDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+ // ---- AntiLambda
+ if( trackAssocME->WhichCandidate() == 5 ){
+ splitCont2[1] = massAL;
+ // Positive daughter
+ fAntiLambdaPosDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+
+ // --------------------------
+ // ------- Negative daughter:
+ splitCont2[2] = ntrack->Pt(); splitCont2[3] = sameSignNegDaug; splitCont2[4] = negdPhiS; splitCont2[5] = negdEtaS;
+ splitCont2[7] = fracNegDaugTPCSharedMap; splitCont2[8] = fracTrigTPCSharedMap - fracNegDaugTPCSharedMap;
+ splitCont2[9] = devNegDaugTrig[7]; splitCont2[10] = tAssoc->DcaNegToPrimVertex(); splitCont2[11] = sameLabelNegDaug;
+
+ // ---- K0s
+ if( trackAssocME->WhichCandidate() == 3 ){
+ splitCont2[1] = massK0s;
+ // Negative daughter
+ fK0sNegDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+ // ---- Lambda
+ if( trackAssocME->WhichCandidate() == 4 ){
+ splitCont2[1] = massL;
+ // Negative daughter
+ fLambdaNegDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+ // ---- AntiLambda
+ if( trackAssocME->WhichCandidate() == 5 ){
+ splitCont2[1] = massAL;
+ // Negative daughter
+ fAntiLambdaNegDaugdPhiSdEtaS[curCentBin]->Fill(splitCont2);
+ }
+ } // end selection in |delta(eta)| < 0.1, |delta(phi)| < 0.1
+
+ // ----------------------------------------------------------------
+ // Reject the 'fake' correlation due to the TPC shared clusters
+ // between trigger particle and one of the daughter tracks
+ // The rejection will affect more the correlations:
+ // - Trigger track - Positive track (from Lambda with pt above 3 GeV/c)
+ // - Trigger track - Negative track (from AntiLambda with pt above 3 GeV/c)
+ /* if( fracTrigTPCSharedMap>0.5 &&
+ ( ( sameSignPosDaug==1 && TMath::Abs(fracTrigTPCSharedMap - fracPosDaugTPCSharedMap) < fDiffTrigDaugFracTPCSharedCls ) ||
+ ( sameSignNegDaug==1 && TMath::Abs(fracTrigTPCSharedMap - fracNegDaugTPCSharedMap) < fDiffTrigDaugFracTPCSharedCls ) ) )*/
+
+ if( (fracTrigTPCSharedMap > fFracTPCcls) || (fracPosDaugTPCSharedMap > fFracTPCcls) || (fracNegDaugTPCSharedMap > fFracTPCcls) )
+ continue;
+
// ----------------------------------------------------------------------------
// *******************
// ==== Correlations K0s invariant mass peak ==== //
// +++++++++++ Pt bin & centrality
- if(trigTPCMapOk && posDaugTPCMapOk && negDaugTPCMapOk) // splitting rejection according to the TPC shared map
- fK0sdPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massK0s);
+ fK0sdPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massK0s);
// ==== Correlations K0s invariant mass peak ==== //
if (TMath::Abs(mK0s-massK0s) < 3*sK0s) {
- // -------------- splitting checks ---------------
- if(TMath::Abs(dPhi)<0.06 && TMath::Abs(dEta)<0.06){
-
- //// -------- TPC Map ------- ////
- if(!trigTPCMapOk && !posDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ptrack->P(),0.5);
- if(!trigTPCMapOk && !negDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ntrack->P(),1.5);
-
- //// -------- Covariance matrix ------- ////
- // positive daughter
- splitCont[0] = pt; splitCont[1] = ptrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devPosDaugTrig[ll]; splitCont[3] = ll;
- fK0sPosDaugSplCheckCovMat->Fill(splitCont);
- }
-
- // negative daughter
- splitCont[0] = pt; splitCont[1] = ntrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devNegDaugTrig[ll]; splitCont[3] = ll;
- fK0sNegDaugSplCheckCovMat->Fill(splitCont);
- }
-
- }
-
- // ----------------------------------------------
-
- if(radio<0.1){
+ if(radio<0.02){
fK0sSpatialRes->Fill(dPhi,res,lt);
}
if(radio < 0.4){
RecCascade(tTrig,ntrack,ptrack,"K0s");
RecCascade(tTrig,ptrack,ntrack,"K0s");
}
-
-
+
}
// ==== Correlations K0s background ==== //
if( TMath::Abs(mK0s-massK0s + 6.5*sK0s) < 1.5*sK0s ||
// ==== Correlations Lambda invariant mass peak ==== //
// +++++++++++ Pt bin & centrality
- if(trigTPCMapOk && posDaugTPCMapOk && negDaugTPCMapOk) // splitting rejection according to the TPC shared map
- fLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massL);
+ fLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massL);
// ==== Correlations Lambda invariant mass peak ==== //
- if (TMath::Abs(mL-massL) < 3*sL) {
-
- // -------------- splitting checks ---------------
- if(TMath::Abs(dPhi)<0.06 && TMath::Abs(dEta)<0.06){
-
- //// -------- TPC Map ------- ////
- if(!trigTPCMapOk && !posDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ptrack->P(),2.5);
- if(!trigTPCMapOk && !negDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ntrack->P(),3.5);
-
- //// -------- Covariance matrix ------- ////
- // positive daughter
- splitCont[0] = pt; splitCont[1] = ptrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devPosDaugTrig[ll]; splitCont[3] = ll;
- fLambdaPosDaugSplCheckCovMat->Fill(splitCont);
- }
+ if (TMath::Abs(mL-massL) < 3*sL) {
- // negative daughter
- splitCont[0] = pt; splitCont[1] = ntrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devNegDaugTrig[ll]; splitCont[3] = ll;
- fLambdaNegDaugSplCheckCovMat->Fill(splitCont);
- }
-
- }
-
- // ----------------------------------------------
-
- if(radio<0.1)
+ if(radio<0.02)
fLambdaSpatialRes->Fill(dPhi,res,lt);
if(radio < 0.4){
fHistArmPodBckg->Fill(lAlphaV0,lPtArmV0,2);
// ==== Correlations Lambda background ==== //
if( TMath::Abs(mL-massL + 6.5*sL) < 1.5*sL ||
TMath::Abs(mL-massL - 6.5*sL) < 1.5*sL ) {
-
+
// Only fills the histograms when it is a triggered event
if(j==0){
fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0,3);
}// End selection in the correlation peak
- } // End bacground selection
+ } // End background selection
}// End Lambda selection
// *******************
else if( fCollision.Contains("PbPb2011") )
sL = kCteAntiLambda2011[curCentBin] + kLinearAntiLambda2011[curCentBin]*pt;
-
// ==== Correlations Lambda invariant mass peak ==== //
// +++++++++++ Pt bin & centrality
- if(trigTPCMapOk && posDaugTPCMapOk && negDaugTPCMapOk) // splitting rejection according to the TPC shared map
- fAntiLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massAL);
+ fAntiLambdadPhidEtaPtL[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(dPhi,dEta,massAL);
// ==== Correlations AntiLambda invariant mass peak ==== //
if (TMath::Abs(mL-massAL) < 3*sL) {
- // -------------- splitting checks ---------------
- if(TMath::Abs(dPhi)<0.06 && TMath::Abs(dEta)<0.06){
-
- //// -------- TPC Map ------- ////
- if(!trigTPCMapOk && !posDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ptrack->P(),4.5);
- if(!trigTPCMapOk && !negDaugTPCMapOk)
- fSharedClsTrigDaug->Fill(pt,ntrack->P(),5.5);
-
- //// -------- Covariance matrix ------- ////
- // positive daughter
- splitCont[0] = pt; splitCont[1] = ptrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devPosDaugTrig[ll]; splitCont[3] = ll;
- fAntiLambdaPosDaugSplCheckCovMat->Fill(splitCont);
- }
-
- // negative daughter
- splitCont[0] = pt; splitCont[1] = ntrack->P();
- for(Int_t ll=0; ll<=6; ll++){
- splitCont[2] = devNegDaugTrig[ll]; splitCont[3] = ll;
- fAntiLambdaNegDaugSplCheckCovMat->Fill(splitCont);
- }
-
- }
-
- // ----------------------------------------------
-
if(radio<0.1)
fAntiLambdaSpatialRes->Fill(dPhi,res,lt);
if(radio < 0.4){
// ==== Correlations AntiLambda background ==== //
if( (TMath::Abs(mL-massAL + 6.5*sL) < 1.5*sL) ||
(TMath::Abs(mL-massAL - 6.5*sL) < 1.5*sL) ){
-
+
+ // ----------------------------------------------
+
// Only fills the histograms when it is a triggered event
if(j==0){
fHistArmenterosPodolanski->Fill(lAlphaV0,lPtArmV0,5);
} // End loop over associated particles
+
+ // Filling information of the trigger particle
+ // after the rejection in the cut of shared TPC cls
+ fTriggerEtaPhi->Fill(trig->Phi(),trig->Eta());
+ fTriggerPtCent->Fill(trig->Pt(),centrality,zv);
+
} // End loop over trigger particles
//-------------------------------------------------------------
// Mixing
//-------------------------------------------------------------
- Double_t phiTrigME=0, etaTrigME=0, phiAssocME=0, etaAssocME=0;
- Double_t deltaPhi=0, deltaEta=0;
+ /*
+ Double_t phiTrigME=0, etaTrigME=0, phiAssocME=0, etaAssocME=0;
+ Double_t deltaPhi=0, deltaEta=0;
- TList *evMixList = fMEList[curCentBin*kNVtxZ+curVtxBin];
- Int_t nMixed = evMixList->GetSize();
+ TList *evMixList = fMEList[curCentBin*kNVtxZ+curVtxBin];
+ Int_t nMixed = evMixList->GetSize();
- if( nMixed>0 && fAssocParticles->GetEntriesFast() >= 0 ){
+ if( nMixed>0 && fAssocParticles->GetEntriesFast() >= 0 ){
for(Int_t ii=0; ii<nMixed; ii++){
- AliMiniParticle* trackTriggerME = (AliMiniParticle*) (evMixList->At(ii));
- phiTrigME = trackTriggerME->Phi();
- etaTrigME = trackTriggerME->Eta();
+ AliMiniParticle* trackTriggerME = (AliMiniParticle*) (evMixList->At(ii));
+ phiTrigME = trackTriggerME->Phi();
+ etaTrigME = trackTriggerME->Eta();
- // --- V0 associated particles
- for(Int_t j=0; j<fAssocParticles->GetEntriesFast(); j++){
+ // --- V0 associated particles
+ for(Int_t j=0; j<fAssocParticles->GetEntriesFast(); j++){
- AliMiniParticle* trackAssocME = (AliMiniParticle*) (fAssocParticles->At(j));
- if( CentBin(trackTriggerME->Centrality()) != CentBin(trackAssocME->Centrality()) ) continue;
- if( VtxBin(trackTriggerME->VtxZ()) != VtxBin(trackAssocME->VtxZ()) ) continue;
- if( trackAssocME->WhichCandidate() == 2 ) continue;
+ AliMiniParticle* trackAssocME = (AliMiniParticle*) (fAssocParticles->At(j));
+ if( CentBin(trackTriggerME->Centrality()) != CentBin(trackAssocME->Centrality()) ) continue;
+ if( VtxBin(trackTriggerME->VtxZ()) != VtxBin(trackAssocME->VtxZ()) ) continue;
+ if( trackAssocME->WhichCandidate() == 2 ) continue;
- AliAODv0 *tAssoc=fAOD->GetV0(trackAssocME->ID());
- pt = tAssoc->Pt();
+ AliAODv0 *tAssoc=fAOD->GetV0(trackAssocME->ID());
+ pt = tAssoc->Pt();
- Bool_t IsSelected = kFALSE;
- // K0s
- if( trackAssocME->WhichCandidate() == 3 ){
- massK0s = tAssoc->MassK0Short();
- mK0s = TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass();
- if( fCollision.Contains("PbPb2010") )
- sK0s = kCteK0s2010[curCentBin] + kLinearK0s2010[curCentBin]*pt;
- else if( fCollision.Contains("PbPb2011") )
- sK0s = kCteK0s2011[curCentBin] + kLinearK0s2011[curCentBin]*pt;
+ Bool_t IsSelected = kFALSE;
+ // K0s
+ if( trackAssocME->WhichCandidate() == 3 ){
+ massK0s = tAssoc->MassK0Short();
+ mK0s = TDatabasePDG::Instance()->GetParticle(kK0Short)->Mass();
+ if( fCollision.Contains("PbPb2010") )
+ sK0s = kCteK0s2010[curCentBin] + kLinearK0s2010[curCentBin]*pt;
+ else if( fCollision.Contains("PbPb2011") )
+ sK0s = kCteK0s2011[curCentBin] + kLinearK0s2011[curCentBin]*pt;
- if (TMath::Abs(mK0s-massK0s) < 3*sK0s) IsSelected = kTRUE;
- }
- // Lambda
- if( trackAssocME->WhichCandidate() == 4 ){
- massL = tAssoc->MassLambda();
- mL = TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass();
- if( fCollision.Contains("PbPb2010") )
- sL = kCteLambda2010[curCentBin] + kLinearLambda2010[curCentBin]*pt;
- else if( fCollision.Contains("PbPb2011") )
- sL = kCteLambda2011[curCentBin] + kLinearLambda2011[curCentBin]*pt;
-
- if (TMath::Abs(mL-massL) < 3*sL) IsSelected = kTRUE;
- }
- // AntiLambda
- if( trackAssocME->WhichCandidate() == 5 ){
- massAL = tAssoc->MassAntiLambda();
- mL = TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass();
- if( fCollision.Contains("PbPb2010") )
- sL = kCteAntiLambda2010[curCentBin] + kLinearAntiLambda2010[curCentBin]*pt;
- else if( fCollision.Contains("PbPb2011") )
- sL = kCteAntiLambda2011[curCentBin] + kLinearAntiLambda2011[curCentBin]*pt;
+ if (TMath::Abs(mK0s-massK0s) < 3*sK0s) IsSelected = kTRUE;
+ }
+ // Lambda
+ if( trackAssocME->WhichCandidate() == 4 ){
+ massL = tAssoc->MassLambda();
+ mL = TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass();
+ if( fCollision.Contains("PbPb2010") )
+ sL = kCteLambda2010[curCentBin] + kLinearLambda2010[curCentBin]*pt;
+ else if( fCollision.Contains("PbPb2011") )
+ sL = kCteLambda2011[curCentBin] + kLinearLambda2011[curCentBin]*pt;
+
+ if (TMath::Abs(mL-massL) < 3*sL) IsSelected = kTRUE;
+ }
+ // AntiLambda
+ if( trackAssocME->WhichCandidate() == 5 ){
+ massAL = tAssoc->MassAntiLambda();
+ mL = TDatabasePDG::Instance()->GetParticle(kLambda0)->Mass();
+ if( fCollision.Contains("PbPb2010") )
+ sL = kCteAntiLambda2010[curCentBin] + kLinearAntiLambda2010[curCentBin]*pt;
+ else if( fCollision.Contains("PbPb2011") )
+ sL = kCteAntiLambda2011[curCentBin] + kLinearAntiLambda2011[curCentBin]*pt;
- if (TMath::Abs(mL-massAL) < 3*sL) IsSelected = kTRUE;
- }
+ if (TMath::Abs(mL-massAL) < 3*sL) IsSelected = kTRUE;
+ }
- if(!IsSelected) continue;
+ if(!IsSelected) continue;
- phiAssocME = trackAssocME->Phi();
- etaAssocME = trackAssocME->Eta();
+ phiAssocME = trackAssocME->Phi();
+ etaAssocME = trackAssocME->Eta();
- deltaPhi = dPHI(phiTrigME,phiAssocME);
- deltaEta = etaTrigME - etaAssocME;
+ deltaPhi = dPHI(phiTrigME,phiAssocME);
+ deltaEta = etaTrigME - etaAssocME;
- Int_t binPtv0 = PtBin( trackAssocME->Pt() );
- if(binPtv0==-1) continue;
+ Int_t binPtv0 = PtBin( trackAssocME->Pt() );
+ if(binPtv0==-1) continue;
- if( trackAssocME->WhichCandidate() == 3 ) {
- fK0sdPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);}
- else if( trackAssocME->WhichCandidate() == 4 )
- fLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
- else if( trackAssocME->WhichCandidate() == 5 )
- fAntiLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
+ if( trackAssocME->WhichCandidate() == 3 ) {
+ fK0sdPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);}
+ else if( trackAssocME->WhichCandidate() == 4 )
+ fLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
+ else if( trackAssocME->WhichCandidate() == 5 )
+ fAntiLambdadPhidEtaME[curCentBin*kN1*kNVtxZ + binPtv0*kNVtxZ + curVtxBin]->Fill(deltaPhi,deltaEta);
- } // End loop over V0's
+ } // End loop over V0's
}
- }
-
+ }
+ */
//--------------------------------------------------------
//Add the current event to the list of events for mixing
//--------------------------------------------------------
-
+ /*
//Add current event to buffer and Remove redundant events
if(fTriggerParticles->GetEntriesFast()>=0){
- for(Int_t ii=0; ii<(fTriggerParticles->GetEntriesFast()); ii++){
- AliMiniParticle* trkTrig = (AliMiniParticle*) fTriggerParticles->At(ii);
- //cout << trkTrig->Pt() << " " << ii << endl;
+ for(Int_t ii=0; ii<(fTriggerParticles->GetEntriesFast()); ii++){
+ AliMiniParticle* trkTrig = (AliMiniParticle*) fTriggerParticles->At(ii);
+ //cout << trkTrig->Pt() << " " << ii << endl;
- if(evMixList->GetSize() < nMaxEvMix)
- evMixList->AddFirst(trkTrig);
- /*
- if(evMixList->GetSize() >= nMaxEvMix) {
- AliMiniParticle *tmp = (AliMiniParticle*) (evMixList->Last()) ;
- evMixList->RemoveLast();
- delete tmp;
- }
- */
+ if(evMixList->GetSize() < nMaxEvMix)
+ evMixList->AddFirst(trkTrig);
+ / *
+ if(evMixList->GetSize() >= nMaxEvMix) {
+ AliMiniParticle *tmp = (AliMiniParticle*) (evMixList->Last()) ;
+ evMixList->RemoveLast();
+ delete tmp;
+ }
+ * /
- }// End loop over fTriggerParticles
+ }// End loop over fTriggerParticles
}// End adding trigger particles to buffers
-
+ */
}
//___________________________________________________________________________________________