Skip to content

Commit 14859b3

Browse files
romainschotteralibuild
romainschotter
and
alibuild
authored
[PWGLF] Add QA histograms (#17440)
Co-authored-by: ALICE Builder <alibuild@users.noreply.github.com>
1 parent 9da4f3d commit 14859b3

1 file changed

Lines changed: 85 additions & 13 deletions

File tree

‎PWGLF/Tasks/Strangeness/forwardlambdakzeroanalysis.cxx‎

Lines changed: 85 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -707,8 +707,8 @@ struct forwardlambdakzeroanalysis {
707707
if (doPlainTopoQA) {
708708
// All candidates received
709709
histos.add("hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
710-
histos.add("hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
711-
histos.add("hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
710+
histos.add("hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
711+
histos.add("hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
712712
histos.add("hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
713713
histos.add("hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
714714
histos.add("hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -727,6 +727,16 @@ struct forwardlambdakzeroanalysis {
727727
histos.add("hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
728728
histos.add("hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
729729
histos.add("hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
730+
if (doCompleteTopoQA) {
731+
histos.add("hV0InvMassK0sVsLambda", "hV0InvMassK0sVsLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
732+
histos.add("hV0InvMassK0sVsAntiLambda", "hV0InvMassK0sVsAntiLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
733+
histos.add("hV0InvMassK0sVsD0", "hV0InvMassK0sVsD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
734+
histos.add("hV0InvMassK0sVsAntiD0", "hV0InvMassK0sVsAntiD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
735+
histos.add("hV0InvMassLambdaVsD0", "hV0InvMassLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
736+
histos.add("hV0InvMassAntiLambdaVsD0", "hV0InvMassAntiLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
737+
histos.add("hV0InvMassLambdaVsAntiD0", "hV0InvMassLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
738+
histos.add("hV0InvMassAntiLambdaVsAntiD0", "hV0InvMassAntiLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
739+
}
730740

731741
histos.add("hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
732742
histos.add("hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -741,7 +751,7 @@ struct forwardlambdakzeroanalysis {
741751
if (analyseK0Short) {
742752
histos.add("K0Short/hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
743753
histos.add("K0Short/hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
744-
histos.add("K0Short/hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
754+
histos.add("K0Short/hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVxy});
745755
histos.add("K0Short/hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
746756
histos.add("K0Short/hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
747757
histos.add("K0Short/hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -760,6 +770,12 @@ struct forwardlambdakzeroanalysis {
760770
histos.add("K0Short/hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
761771
histos.add("K0Short/hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
762772
histos.add("K0Short/hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
773+
if (doCompleteTopoQA) {
774+
histos.add("K0Short/hV0InvMassK0sVsLambda", "hV0InvMassK0sVsLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
775+
histos.add("K0Short/hV0InvMassK0sVsAntiLambda", "hV0InvMassK0sVsAntiLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
776+
histos.add("K0Short/hV0InvMassK0sVsD0", "hV0InvMassK0sVsD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
777+
histos.add("K0Short/hV0InvMassK0sVsAntiD0", "hV0InvMassK0sVsAntiD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
778+
}
763779

764780
histos.add("K0Short/hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
765781
histos.add("K0Short/hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -775,7 +791,7 @@ struct forwardlambdakzeroanalysis {
775791
if (analyseLambda) {
776792
histos.add("Lambda/hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
777793
histos.add("Lambda/hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
778-
histos.add("Lambda/hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
794+
histos.add("Lambda/hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVxy});
779795
histos.add("Lambda/hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
780796
histos.add("Lambda/hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
781797
histos.add("Lambda/hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -794,6 +810,11 @@ struct forwardlambdakzeroanalysis {
794810
histos.add("Lambda/hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
795811
histos.add("Lambda/hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
796812
histos.add("Lambda/hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
813+
if (doCompleteTopoQA) {
814+
histos.add("Lambda/hV0InvMassK0sVsLambda", "hV0InvMassK0sVsLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
815+
histos.add("Lambda/hV0InvMassLambdaVsD0", "hV0InvMassLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
816+
histos.add("Lambda/hV0InvMassLambdaVsAntiD0", "hV0InvMassLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
817+
}
797818

798819
histos.add("Lambda/hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
799820
histos.add("Lambda/hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -809,7 +830,7 @@ struct forwardlambdakzeroanalysis {
809830
if (analyseAntiLambda) {
810831
histos.add("AntiLambda/hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
811832
histos.add("AntiLambda/hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
812-
histos.add("AntiLambda/hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
833+
histos.add("AntiLambda/hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVxy});
813834
histos.add("AntiLambda/hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
814835
histos.add("AntiLambda/hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
815836
histos.add("AntiLambda/hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -828,6 +849,11 @@ struct forwardlambdakzeroanalysis {
828849
histos.add("AntiLambda/hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
829850
histos.add("AntiLambda/hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
830851
histos.add("AntiLambda/hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
852+
if (doCompleteTopoQA) {
853+
histos.add("AntiLambda/hV0InvMassK0sVsAntiLambda", "hV0InvMassK0sVsAntiLambda", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisLambdaMass});
854+
histos.add("AntiLambda/hV0InvMassAntiLambdaVsD0", "hV0InvMassAntiLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
855+
histos.add("AntiLambda/hV0InvMassAntiLambdaVsAntiD0", "hV0InvMassAntiLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
856+
}
831857

832858
histos.add("AntiLambda/hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
833859
histos.add("AntiLambda/hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -843,7 +869,7 @@ struct forwardlambdakzeroanalysis {
843869
if (analyseD0) {
844870
histos.add("D0/hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
845871
histos.add("D0/hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
846-
histos.add("D0/hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
872+
histos.add("D0/hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVxy});
847873
histos.add("D0/hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
848874
histos.add("D0/hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
849875
histos.add("D0/hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -862,6 +888,11 @@ struct forwardlambdakzeroanalysis {
862888
histos.add("D0/hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
863889
histos.add("D0/hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
864890
histos.add("D0/hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
891+
if (doCompleteTopoQA) {
892+
histos.add("D0/hV0InvMassK0sVsD0", "hV0InvMassK0sVsD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
893+
histos.add("D0/hV0InvMassLambdaVsD0", "hV0InvMassLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
894+
histos.add("D0/hV0InvMassAntiLambdaVsD0", "hV0InvMassAntiLambdaVsD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
895+
}
865896

866897
histos.add("D0/hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
867898
histos.add("D0/hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -877,7 +908,7 @@ struct forwardlambdakzeroanalysis {
877908
if (analyseAntiD0) {
878909
histos.add("AntiD0/hPosDCAToPVxy", "hPosDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
879910
histos.add("AntiD0/hNegDCAToPVxy", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVz});
880-
histos.add("AntiD0/hPosDCAToPVz", "hNegDCAToPVxy", kTH1D, {axisConfigurations.axisDCAtoPVxy});
911+
histos.add("AntiD0/hPosDCAToPVz", "hPosDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVxy});
881912
histos.add("AntiD0/hNegDCAToPVz", "hNegDCAToPVz", kTH1D, {axisConfigurations.axisDCAtoPVz});
882913
histos.add("AntiD0/hDCADaughters", "hDCADaughters", kTH1D, {axisConfigurations.axisDCAdau});
883914
histos.add("AntiD0/hCosPA", "hCosPA", kTH1D, {axisConfigurations.axisCosPA});
@@ -896,6 +927,11 @@ struct forwardlambdakzeroanalysis {
896927
histos.add("AntiD0/hV0InvMassAntiLambda", "hV0InvMassAntiLambda", kTH1D, {axisConfigurations.axisLambdaMass});
897928
histos.add("AntiD0/hV0InvMassD0", "hV0InvMassD0", kTH1D, {axisConfigurations.axisD0Mass});
898929
histos.add("AntiD0/hV0InvMassAntiD0", "hV0InvMassAntiD0", kTH1D, {axisConfigurations.axisD0Mass});
930+
if (doCompleteTopoQA) {
931+
histos.add("AntiD0/hV0InvMassK0sVsAntiD0", "hV0InvMassK0sVsAntiD0", kTH2D, {axisConfigurations.axisK0Mass, axisConfigurations.axisD0Mass});
932+
histos.add("AntiD0/hV0InvMassLambdaVsAntiD0", "hV0InvMassLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
933+
histos.add("AntiD0/hV0InvMassAntiLambdaVsAntiD0", "hV0InvMassAntiLambdaVsAntiD0", kTH2D, {axisConfigurations.axisLambdaMass, axisConfigurations.axisD0Mass});
934+
}
899935

900936
histos.add("AntiD0/hPositiveMFTcls", "hPositiveMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
901937
histos.add("AntiD0/hNegativeMFTcls", "hNegativeMFTcls", kTH1D, {axisConfigurations.axisMFTclus});
@@ -1313,6 +1349,16 @@ struct forwardlambdakzeroanalysis {
13131349
histos.fill(HIST("hV0InvMassAntiLambda"), invMassAntiLambda);
13141350
histos.fill(HIST("hV0InvMassD0"), invMassD0);
13151351
histos.fill(HIST("hV0InvMassAntiD0"), invMassAntiD0);
1352+
if (doCompleteTopoQA) {
1353+
histos.fill(HIST("hV0InvMassK0sVsLambda"), invMassK0Short, invMassLambda);
1354+
histos.fill(HIST("hV0InvMassK0sVsAntiLambda"), invMassK0Short, invMassAntiLambda);
1355+
histos.fill(HIST("hV0InvMassK0sVsD0"), invMassK0Short, invMassD0);
1356+
histos.fill(HIST("hV0InvMassK0sVsAntiD0"), invMassK0Short, invMassAntiD0);
1357+
histos.fill(HIST("hV0InvMassLambdaVsD0"), invMassLambda, invMassD0);
1358+
histos.fill(HIST("hV0InvMassAntiLambdaVsD0"), invMassAntiLambda, invMassD0);
1359+
histos.fill(HIST("hV0InvMassLambdaVsAntiD0"), invMassLambda, invMassAntiD0);
1360+
histos.fill(HIST("hV0InvMassAntiLambdaVsAntiD0"), invMassAntiLambda, invMassAntiD0);
1361+
}
13161362

13171363
histos.fill(HIST("hPositiveMFTcls"), v0.posNclusters);
13181364
histos.fill(HIST("hNegativeMFTcls"), v0.negNclusters);
@@ -1374,6 +1420,12 @@ struct forwardlambdakzeroanalysis {
13741420
histos.fill(HIST("K0Short/hV0InvMassAntiLambda"), invMassAntiLambda);
13751421
histos.fill(HIST("K0Short/hV0InvMassD0"), invMassD0);
13761422
histos.fill(HIST("K0Short/hV0InvMassAntiD0"), invMassAntiD0);
1423+
if (doCompleteTopoQA) {
1424+
histos.fill(HIST("K0Short/hV0InvMassK0sVsLambda"), invMassK0Short, invMassLambda);
1425+
histos.fill(HIST("K0Short/hV0InvMassK0sVsAntiLambda"), invMassK0Short, invMassAntiLambda);
1426+
histos.fill(HIST("K0Short/hV0InvMassK0sVsD0"), invMassK0Short, invMassD0);
1427+
histos.fill(HIST("K0Short/hV0InvMassK0sVsAntiD0"), invMassK0Short, invMassAntiD0);
1428+
}
13771429

13781430
histos.fill(HIST("K0Short/hPositiveMFTcls"), v0.posNclusters);
13791431
histos.fill(HIST("K0Short/hNegativeMFTcls"), v0.negNclusters);
@@ -1425,6 +1477,11 @@ struct forwardlambdakzeroanalysis {
14251477
histos.fill(HIST("Lambda/hV0InvMassAntiLambda"), invMassAntiLambda);
14261478
histos.fill(HIST("Lambda/hV0InvMassD0"), invMassD0);
14271479
histos.fill(HIST("Lambda/hV0InvMassAntiD0"), invMassAntiD0);
1480+
if (doCompleteTopoQA) {
1481+
histos.fill(HIST("Lambda/hV0InvMassK0sVsLambda"), invMassK0Short, invMassLambda);
1482+
histos.fill(HIST("Lambda/hV0InvMassLambdaVsD0"), invMassLambda, invMassD0);
1483+
histos.fill(HIST("Lambda/hV0InvMassLambdaVsAntiD0"), invMassLambda, invMassAntiD0);
1484+
}
14281485

14291486
histos.fill(HIST("Lambda/hPositiveMFTcls"), v0.posNclusters);
14301487
histos.fill(HIST("Lambda/hNegativeMFTcls"), v0.negNclusters);
@@ -1476,6 +1533,11 @@ struct forwardlambdakzeroanalysis {
14761533
histos.fill(HIST("AntiLambda/hV0InvMassAntiLambda"), invMassAntiLambda);
14771534
histos.fill(HIST("AntiLambda/hV0InvMassD0"), invMassD0);
14781535
histos.fill(HIST("AntiLambda/hV0InvMassAntiD0"), invMassAntiD0);
1536+
if (doCompleteTopoQA) {
1537+
histos.fill(HIST("AntiLambda/hV0InvMassK0sVsAntiLambda"), invMassK0Short, invMassAntiLambda);
1538+
histos.fill(HIST("AntiLambda/hV0InvMassAntiLambdaVsD0"), invMassAntiLambda, invMassD0);
1539+
histos.fill(HIST("AntiLambda/hV0InvMassAntiLambdaVsAntiD0"), invMassAntiLambda, invMassAntiD0);
1540+
}
14791541

14801542
histos.fill(HIST("AntiLambda/hPositiveMFTcls"), v0.posNclusters);
14811543
histos.fill(HIST("AntiLambda/hNegativeMFTcls"), v0.negNclusters);
@@ -1527,6 +1589,11 @@ struct forwardlambdakzeroanalysis {
15271589
histos.fill(HIST("D0/hV0InvMassAntiLambda"), invMassAntiLambda);
15281590
histos.fill(HIST("D0/hV0InvMassD0"), invMassD0);
15291591
histos.fill(HIST("D0/hV0InvMassAntiD0"), invMassAntiD0);
1592+
if (doCompleteTopoQA) {
1593+
histos.fill(HIST("D0/hV0InvMassK0sVsD0"), invMassK0Short, invMassD0);
1594+
histos.fill(HIST("D0/hV0InvMassLambdaVsD0"), invMassLambda, invMassD0);
1595+
histos.fill(HIST("D0/hV0InvMassAntiLambdaVsD0"), invMassAntiLambda, invMassD0);
1596+
}
15301597

15311598
histos.fill(HIST("D0/hPositiveMFTcls"), v0.posNclusters);
15321599
histos.fill(HIST("D0/hNegativeMFTcls"), v0.negNclusters);
@@ -1578,6 +1645,11 @@ struct forwardlambdakzeroanalysis {
15781645
histos.fill(HIST("AntiD0/hV0InvMassAntiLambda"), invMassAntiLambda);
15791646
histos.fill(HIST("AntiD0/hV0InvMassD0"), invMassD0);
15801647
histos.fill(HIST("AntiD0/hV0InvMassAntiD0"), invMassAntiD0);
1648+
if (doCompleteTopoQA) {
1649+
histos.fill(HIST("AntiD0/hV0InvMassK0sVsAntiD0"), invMassK0Short, invMassAntiD0);
1650+
histos.fill(HIST("AntiD0/hV0InvMassLambdaVsAntiD0"), invMassLambda, invMassAntiD0);
1651+
histos.fill(HIST("AntiD0/hV0InvMassAntiLambdaVsAntiD0"), invMassAntiLambda, invMassAntiD0);
1652+
}
15811653

15821654
histos.fill(HIST("AntiD0/hPositiveMFTcls"), v0.posNclusters);
15831655
histos.fill(HIST("AntiD0/hNegativeMFTcls"), v0.negNclusters);
@@ -2069,9 +2141,9 @@ struct forwardlambdakzeroanalysis {
20692141

20702142
// get decay vertex coordinates
20712143
Vec3D vtx = fitter.getPCACandidate();
2072-
pairInfo.X = vtx[0];
2073-
pairInfo.Y = vtx[1];
2074-
pairInfo.Z = vtx[2];
2144+
pairInfo.X = vtx[0] - collision.posX();
2145+
pairInfo.Y = vtx[1] - collision.posY();
2146+
pairInfo.Z = vtx[2] - collision.posZ();
20752147

20762148
// get daughter DCA to PV
20772149
pairInfo.dcaPosToPVxy = std::sqrt(dcaPosToPVx * dcaPosToPVx + dcaPosToPVy * dcaPosToPVy);
@@ -2101,16 +2173,16 @@ struct forwardlambdakzeroanalysis {
21012173
pairInfo.OpAngle = track1Momentum.Angle(track2Momentum);
21022174

21032175
// Radius
2104-
pairInfo.Radius = std::sqrt(vtx[0] * vtx[0] + vtx[1] * vtx[1]);
2176+
pairInfo.Radius = std::sqrt(pairInfo.X * pairInfo.X + pairInfo.Y * pairInfo.Y);
21052177

21062178
// Dist over tot mom.
21072179
float px = pairInfo.positiveMomentum[0] + pairInfo.negativeMomentum[0];
21082180
float py = pairInfo.positiveMomentum[1] + pairInfo.negativeMomentum[1];
21092181
float pz = pairInfo.positiveMomentum[2] + pairInfo.negativeMomentum[2];
2110-
pairInfo.DistOverTotMom = std::sqrt(vtx[0] * vtx[0] + vtx[1] * vtx[1] + vtx[2] * vtx[2]) / std::sqrt(px * px + py * py + pz * pz);
2182+
pairInfo.DistOverTotMom = std::sqrt(pairInfo.X * pairInfo.X + pairInfo.Y * pairInfo.Y + pairInfo.Z * pairInfo.Z) / std::sqrt(px * px + py * py + pz * pz);
21112183

21122184
// Z dist over pz
2113-
pairInfo.ZdistOverPz = vtx[2] / pz;
2185+
pairInfo.ZdistOverPz = pairInfo.Z / pz;
21142186

21152187
// V0 Momenta
21162188
pairInfo.pT = std::sqrt(px * px + py * py);

0 commit comments

Comments
 (0)