Skip to content

Commit 1556a17

Browse files
mjkim525mjkim525
committed
Use DeltaR checks in Xi1820
1 parent 6681d7a commit 1556a17

1 file changed

Lines changed: 38 additions & 36 deletions

File tree

‎PWGLF/Tasks/Resonances/xi1820Analysis.cxx‎

Lines changed: 38 additions & 36 deletions
Original file line numberDiff line numberDiff line change
@@ -155,7 +155,7 @@ struct Xi1820Analysis {
155155
Configurable<bool> cfgRotKaon{"cfgRotKaon", true, "Rotate Kaon"};
156156
Configurable<int> cfgNrotBkg{"cfgNrotBkg", 10, "Number of rotated copies (background) per each original candidate"};
157157

158-
Configurable<bool> cfgDelCheck{"cfgDelCheck", true, "Check Delta phi distribution of candidates"};
158+
Configurable<bool> cfgDelCheck{"cfgDelCheck", true, "Check Delta R distribution of candidates"};
159159
Configurable<bool> cfgPVContributor{"cfgPVContributor", true, "Require PV contributor tracks"};
160160
Configurable<bool> cfgPrimaryTrack{"cfgPrimaryTrack", true, "Require primary tracks flags (just for check)"};
161161

@@ -180,7 +180,7 @@ struct Xi1820Analysis {
180180
AxisSpec armenterosAlphaAxis = {200, -1.0, 1.0, "Armenteros alpha"};
181181
AxisSpec armenterosQtAxis = {500, 0.0, 0.5, "Armenteros qt (GeV/c)"};
182182
AxisSpec axisRot = {360, -o2::constants::math::PI, o2::constants::math::PI};
183-
AxisSpec axisDel = {360, -o2::constants::math::PI, o2::constants::math::PI, "#Delta#varphi"};
183+
AxisSpec axisDel = {400, -0.5, 7.0, "#DeltaR"};
184184

185185
// Event QA histograms
186186
histos.add("Event/posZ", "Event vertex Z position", kTH1F, {{200, -20., 20., "V_{z} (cm)"}});
@@ -288,22 +288,22 @@ struct Xi1820Analysis {
288288
histos.add("xi1820/kminus_antilambda/hMassPtCentKminusAntiLambda", "K^{-} + #bar{#Lambda} mass vs pT vs cent", kTH3D, {invMassAxis, ptAxis, centAxis});
289289
histos.add("xi1820/kminus_antilambda/hMassPtCentKminusAntiLambda_Mix", "Mixed event K^{-} + #bar{#Lambda} mass vs pT vs cent", kTH3D, {invMassAxis, ptAxis, centAxis});
290290
} else {
291-
histos.add("xi1820/kplus_lambda/hMassPtCentDelKplusLambda", "K^{+} + #Lambda mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
292-
histos.add("xi1820/kplus_lambda/hMassPtCentDelKplusLambda_Mix", "Mixed event K^{+} + #Lambda mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
291+
histos.add("xi1820/kplus_lambda/hMassPtCentDelKplusLambda", "K^{+} + #Lambda mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
292+
histos.add("xi1820/kplus_lambda/hMassPtCentDelKplusLambda_Mix", "Mixed event K^{+} + #Lambda mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
293293

294294
// K+ Anti-Lambda
295-
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda", "K^{+} + #bar{#Lambda} mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
296-
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Mix", "Mixed event K^{+} + #bar{#Lambda} mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
297-
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Rot", "Rotated background K^{+} + #bar{#Lambda} mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
295+
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda", "K^{+} + #bar{#Lambda} mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
296+
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Mix", "Mixed event K^{+} + #bar{#Lambda} mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
297+
histos.add("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Rot", "Rotated background K^{+} + #bar{#Lambda} mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
298298

299299
// K- Lambda
300-
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda", "K^{-} + #Lambda mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
301-
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Mix", "Mixed event K^{-} + #Lambda mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
302-
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Rot", "Rotated background K^{-} + #Lambda mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
300+
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda", "K^{-} + #Lambda mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
301+
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Mix", "Mixed event K^{-} + #Lambda mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
302+
histos.add("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Rot", "Rotated background K^{-} + #Lambda mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
303303

304304
// K- Anti-Lambda
305-
histos.add("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda", "K^{-} + #bar{#Lambda} mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
306-
histos.add("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda_Mix", "Mixed event K^{-} + #bar{#Lambda} mass vs pT vs cent vs #Delta#varphi", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
305+
histos.add("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda", "K^{-} + #bar{#Lambda} mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
306+
histos.add("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda_Mix", "Mixed event K^{-} + #bar{#Lambda} mass vs pT vs cent vs #DeltaR", kTHnSparseF, {invMassAxis, ptAxis, centAxis, axisDel});
307307
}
308308
}
309309

@@ -365,14 +365,14 @@ struct Xi1820Analysis {
365365
histos.add("xi1820/k0s_antilambda/hMassPtCentK0sAntiLambda_Mix", "Mixed event Xi(1820) mass vs pT vs cent (K^{0}_{S}#bar{#Lambda})", kTH3D, {invMassAxis, ptAxis, centAxis});
366366
histos.add("xi1820/k0s_antilambda/hMassPtCentK0sAntiLambda_Rot", "Rotated background Xi(1820) mass vs pT vs cent (K^{0}_{S}#bar{#Lambda})", kTH3D, {invMassAxis, ptAxis, centAxis});
367367
} else {
368-
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda", "Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
369-
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Mix", "Mixed event Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
370-
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Rot", "Rotated background Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
368+
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda", "Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
369+
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Mix", "Mixed event Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
370+
histos.add("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Rot", "Rotated background Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#Lambda)", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
371371

372372
// K0s + Anti-Lambda
373-
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda", "Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
374-
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Mix", "Mixed event Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
375-
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Rot", "Rotated background Xi(1820) mass vs pT vs cent vs #Delta#varphi (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
373+
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda", "Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
374+
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Mix", "Mixed event Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
375+
histos.add("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Rot", "Rotated background Xi(1820) mass vs pT vs cent vs #DeltaR (K^{0}_{S}#bar{#Lambda})", kTHnSparseF, {{invMassAxis, ptAxis, centAxis, axisDel}});
376376
}
377377
}
378378

@@ -436,9 +436,11 @@ struct Xi1820Analysis {
436436
}
437437

438438
template <typename FirstVectorType, typename SecondVectorType>
439-
static float deltaPhi(const FirstVectorType& first, const SecondVectorType& second)
439+
static float deltaR(const FirstVectorType& first, const SecondVectorType& second)
440440
{
441-
return RecoDecay::constrainAngle(first.Phi() - second.Phi(), -constants::math::PI);
441+
const auto dPhi = RecoDecay::constrainAngle(first.Phi() - second.Phi(), -constants::math::PI);
442+
const auto dEta = first.Eta() - second.Eta();
443+
return std::sqrt(dPhi * dPhi + dEta * dEta);
442444
}
443445

444446
// Lambda/Anti-Lambda selection
@@ -941,15 +943,15 @@ struct Xi1820Analysis {
941943
if (!(additionalConfig.cfgDelCheck))
942944
histos.fill(HIST("xi1820/kplus_lambda/hMassPtCentKplusLambda"), pRes.M(), pRes.Pt(), cent);
943945
else {
944-
auto pCandDel = deltaPhi(pLambda, pKaon);
946+
auto pCandDel = deltaR(pLambda, pKaon);
945947
histos.fill(HIST("xi1820/kplus_lambda/hMassPtCentDelKplusLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
946948
}
947949
} else {
948950
histos.fill(HIST("xi1820/kplus_lambda/hInvMassKplusLambda_Mix"), pRes.M());
949951
if (!(additionalConfig.cfgDelCheck))
950952
histos.fill(HIST("xi1820/kplus_lambda/hMassPtCentKplusLambda_Mix"), pRes.M(), pRes.Pt(), cent);
951953
else {
952-
auto pCandDel = deltaPhi(pLambda, pKaon);
954+
auto pCandDel = deltaR(pLambda, pKaon);
953955
histos.fill(HIST("xi1820/kplus_lambda/hMassPtCentDelKplusLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
954956
}
955957
}
@@ -967,7 +969,7 @@ struct Xi1820Analysis {
967969
if (!(additionalConfig.cfgDelCheck))
968970
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentKplusAntiLambda"), pRes.M(), pRes.Pt(), cent);
969971
else {
970-
auto pCandDel = deltaPhi(pLambda, pKaon);
972+
auto pCandDel = deltaR(pLambda, pKaon);
971973
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
972974
}
973975
if (additionalConfig.cfgFillRotBkg) {
@@ -990,7 +992,7 @@ struct Xi1820Analysis {
990992
if (!(additionalConfig.cfgDelCheck))
991993
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentKplusAntiLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent);
992994
else {
993-
auto pRotDel = additionalConfig.cfgRotKaon ? deltaPhi(lDaughterRot, pLambda) : deltaPhi(lDaughterRot, pKaon);
995+
auto pRotDel = additionalConfig.cfgRotKaon ? deltaR(lDaughterRot, pLambda) : deltaR(lDaughterRot, pKaon);
994996
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent, pRotDel);
995997
}
996998
}
@@ -1016,7 +1018,7 @@ struct Xi1820Analysis {
10161018
if (!(additionalConfig.cfgDelCheck))
10171019
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentKplusAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent);
10181020
else {
1019-
auto pCandDel = deltaPhi(pLambda, pKaon);
1021+
auto pCandDel = deltaR(pLambda, pKaon);
10201022
histos.fill(HIST("xi1820/kplus_antilambda/hMassPtCentDelKplusAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
10211023
}
10221024
}
@@ -1034,7 +1036,7 @@ struct Xi1820Analysis {
10341036
if (!(additionalConfig.cfgDelCheck))
10351037
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentKminusLambda"), pRes.M(), pRes.Pt(), cent);
10361038
else {
1037-
auto pCandDel = deltaPhi(pLambda, pKaon);
1039+
auto pCandDel = deltaR(pLambda, pKaon);
10381040
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentDelKminusLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
10391041
}
10401042

@@ -1058,7 +1060,7 @@ struct Xi1820Analysis {
10581060
if (!(additionalConfig.cfgDelCheck))
10591061
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentKminusLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent);
10601062
else {
1061-
auto pCandDel = additionalConfig.cfgRotKaon ? deltaPhi(lDaughterRot, pLambda) : deltaPhi(lDaughterRot, pKaon);
1063+
auto pCandDel = additionalConfig.cfgRotKaon ? deltaR(lDaughterRot, pLambda) : deltaR(lDaughterRot, pKaon);
10621064
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent, pCandDel);
10631065
}
10641066
}
@@ -1083,7 +1085,7 @@ struct Xi1820Analysis {
10831085
if (!(additionalConfig.cfgDelCheck))
10841086
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentKminusLambda_Mix"), pRes.M(), pRes.Pt(), cent);
10851087
else {
1086-
auto pCandDel = deltaPhi(pLambda, pKaon);
1088+
auto pCandDel = deltaR(pLambda, pKaon);
10871089
histos.fill(HIST("xi1820/kminus_lambda/hMassPtCentDelKminusLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
10881090
}
10891091
}
@@ -1101,15 +1103,15 @@ struct Xi1820Analysis {
11011103
if (!(additionalConfig.cfgDelCheck))
11021104
histos.fill(HIST("xi1820/kminus_antilambda/hMassPtCentKminusAntiLambda"), pRes.M(), pRes.Pt(), cent);
11031105
else {
1104-
auto pCandDel = deltaPhi(pLambda, pKaon);
1106+
auto pCandDel = deltaR(pLambda, pKaon);
11051107
histos.fill(HIST("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
11061108
}
11071109
} else {
11081110
histos.fill(HIST("xi1820/kminus_antilambda/hInvMassKminusAntiLambda_Mix"), pRes.M());
11091111
if (!(additionalConfig.cfgDelCheck))
11101112
histos.fill(HIST("xi1820/kminus_antilambda/hMassPtCentKminusAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent);
11111113
else {
1112-
auto pCandDel = deltaPhi(pLambda, pKaon);
1114+
auto pCandDel = deltaR(pLambda, pKaon);
11131115
histos.fill(HIST("xi1820/kminus_antilambda/hMassPtCentDelKminusAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
11141116
}
11151117
}
@@ -1296,7 +1298,7 @@ struct Xi1820Analysis {
12961298
if (!(additionalConfig.cfgDelCheck))
12971299
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentK0sLambda"), pRes.M(), pRes.Pt(), cent);
12981300
else {
1299-
auto pCandDel = deltaPhi(pLambda, pK0s);
1301+
auto pCandDel = deltaR(pLambda, pK0s);
13001302
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentDelK0sLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
13011303
}
13021304

@@ -1312,7 +1314,7 @@ struct Xi1820Analysis {
13121314
if (!(additionalConfig.cfgDelCheck))
13131315
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentK0sLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent);
13141316
else {
1315-
auto lRotDel = deltaPhi(lDaughterRot, pLambda);
1317+
auto lRotDel = deltaR(lDaughterRot, pLambda);
13161318
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent, lRotDel);
13171319
}
13181320
}
@@ -1337,7 +1339,7 @@ struct Xi1820Analysis {
13371339
if (!(additionalConfig.cfgDelCheck))
13381340
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentK0sLambda_Mix"), pRes.M(), pRes.Pt(), cent);
13391341
else {
1340-
auto pCandDel = deltaPhi(pLambda, pK0s);
1342+
auto pCandDel = deltaR(pLambda, pK0s);
13411343
histos.fill(HIST("xi1820/k0s_lambda/hMassPtCentDelK0sLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
13421344
}
13431345
}
@@ -1354,7 +1356,7 @@ struct Xi1820Analysis {
13541356
if (!(additionalConfig.cfgDelCheck))
13551357
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentK0sAntiLambda"), pRes.M(), pRes.Pt(), cent);
13561358
else {
1357-
auto pCandDel = deltaPhi(pLambda, pK0s);
1359+
auto pCandDel = deltaR(pLambda, pK0s);
13581360
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda"), pRes.M(), pRes.Pt(), cent, pCandDel);
13591361
}
13601362

@@ -1370,7 +1372,7 @@ struct Xi1820Analysis {
13701372
if (!(additionalConfig.cfgDelCheck))
13711373
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentK0sAntiLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent);
13721374
else {
1373-
auto lRotDel = deltaPhi(lDaughterRot, pLambda);
1375+
auto lRotDel = deltaR(lDaughterRot, pLambda);
13741376
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Rot"), lResonanceRot.M(), lResonanceRot.Pt(), cent, lRotDel);
13751377
}
13761378
}
@@ -1395,7 +1397,7 @@ struct Xi1820Analysis {
13951397
if (!(additionalConfig.cfgDelCheck))
13961398
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentK0sAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent);
13971399
else {
1398-
auto pCandDel = deltaPhi(pLambda, pK0s);
1400+
auto pCandDel = deltaR(pLambda, pK0s);
13991401
histos.fill(HIST("xi1820/k0s_antilambda/hMassPtCentDelK0sAntiLambda_Mix"), pRes.M(), pRes.Pt(), cent, pCandDel);
14001402
}
14011403
}

0 commit comments

Comments
 (0)