Skip to content

Commit 9ddec13

Browse files
imsurajprasad786imsurajprasad786
authored
[PWGLF] Adding configurables for inv. mass endpoints (AliceO2Group#10542)
1 parent c84e29d commit 9ddec13

File tree

1 file changed

+34
-30
lines changed

1 file changed

+34
-30
lines changed

‎PWGLF/Tasks/Strangeness/lambdak0sflattenicity.cxx‎

Lines changed: 34 additions & 30 deletions
Original file line numberDiff line numberDiff line change
@@ -101,10 +101,14 @@ struct Lambdak0sflattenicity {
101101

102102
// Configurable for histograms
103103
Configurable<int> nBinsVz{"nBinsVz", 100, "N bins in Vz"};
104-
Configurable<int> nBinsK0sMass{"nBinsK0sMass", 200, "N bins in K0sMass"};
105-
Configurable<int> nBinsLambdaMass{"nBinsLambdaMass", 200,
104+
Configurable<int> nBinsK0sMass{"nBinsK0sMass", 400, "N bins in K0sMass"};
105+
Configurable<int> nBinsLambdaMass{"nBinsLambdaMass", 400,
106106
"N bins in LambdaMass"};
107-
Configurable<int> nBinsXiMass{"nBinsXiMass", 200, "N bins in XiMass"};
107+
Configurable<int> nBinsXiMass{"nBinsXiMass", 400, "N bins in XiMass"};
108+
109+
Configurable<float> kK0sEPshiftfromMass{"kK0sEPshiftfromMass", 0.1, "distance of K0s Inv mass histogram start and end points from PDG mass"};
110+
Configurable<float> kLambdaEPshiftfromMass{"kLambdaEPshiftfromMass", 0.05, "distance of Lambda Inv mass histogram start and end points from PDG mass"};
111+
Configurable<float> kXiEPshiftfromMass{"kXiEPshiftfromMass", 0.05, "distance of Xi Inv mass histogram start and end points from PDG mass"};
108112

109113
Configurable<int> nBinspT{"nBinspT", 250, "N bins in pT"};
110114
Configurable<int> nBinsFlattenicity{"nBinsFlattenicity", 100, "N bins in Flattenicity"};
@@ -215,13 +219,13 @@ struct Lambdak0sflattenicity {
215219
void init(InitContext const&)
216220
{
217221
// Axes
218-
AxisSpec k0sMassAxis = {nBinsK0sMass, 0.45f, 0.55f,
222+
AxisSpec k0sMassAxis = {nBinsK0sMass, 0.49f - kK0sEPshiftfromMass, 0.49f + kK0sEPshiftfromMass,
219223
"#it{M}_{#pi^{+}#pi^{-}} [GeV/#it{c}^{2}]"};
220-
AxisSpec lambdaMassAxis = {nBinsLambdaMass, 1.015f, 1.215f,
224+
AxisSpec lambdaMassAxis = {nBinsLambdaMass, 1.115f - kLambdaEPshiftfromMass, 1.115f + kLambdaEPshiftfromMass,
221225
"#it{M}_{p#pi^{-}} [GeV/#it{c}^{2}]"};
222-
AxisSpec antilambdaMassAxis = {nBinsLambdaMass, 1.015f, 1.215f,
226+
AxisSpec antilambdaMassAxis = {nBinsLambdaMass, 1.115f - kLambdaEPshiftfromMass, 1.115f + kLambdaEPshiftfromMass,
223227
"#it{M}_{#pi^{+}#bar{p}} [GeV/#it{c}^{2}]"};
224-
AxisSpec xiMassAxis = {nBinsXiMass, 1.3f, 1.34f,
228+
AxisSpec xiMassAxis = {nBinsXiMass, 1.32f - kXiEPshiftfromMass, 1.32f + kXiEPshiftfromMass,
225229
"#it{M}_{#Lambda#pi} [GeV/#it{c}^{2}]"};
226230
AxisSpec vertexZAxis = {nBinsVz, -15., 15., "vrtx_{Z} [cm]"};
227231
AxisSpec ptAxis = {nBinspT, 0.0f, 25.0f, "#it{p}_{T} (GeV/#it{c})"};
@@ -1524,7 +1528,7 @@ struct Lambdak0sflattenicity {
15241528
auto v0mcParticle = v0.mcParticle();
15251529
// Cut on dynamic columns for K0s
15261530

1527-
if (v0mcParticle.pdgCode() == 310 && v0.v0cosPA() >= v0settingCosPAK0s &&
1531+
if (v0mcParticle.pdgCode() == PDG_t::kK0Short && v0.v0cosPA() >= v0settingCosPAK0s &&
15281532
v0.v0radius() >= v0settingRadiusK0s &&
15291533
std::abs(posDaughterTrack.tpcNSigmaPi()) <= nSigmaTPCPion &&
15301534
std::abs(negDaughterTrack.tpcNSigmaPi()) <= nSigmaTPCPion &&
@@ -1557,7 +1561,7 @@ struct Lambdak0sflattenicity {
15571561
}
15581562

15591563
// Cut on dynamic columns for Lambda
1560-
if (v0mcParticle.pdgCode() == 3122 &&
1564+
if (v0mcParticle.pdgCode() == PDG_t::kLambda0 &&
15611565
v0.v0cosPA() >= v0settingCosPALambda &&
15621566
v0.v0radius() >= v0settingRadiusLambda &&
15631567
std::abs(posDaughterTrack.tpcNSigmaPr()) <= nSigmaTPCProton &&
@@ -1587,7 +1591,7 @@ struct Lambdak0sflattenicity {
15871591
}
15881592

15891593
// Cut on dynamic columns for AntiLambda
1590-
if (v0mcParticle.pdgCode() == -3122 &&
1594+
if (v0mcParticle.pdgCode() == PDG_t::kLambda0Bar &&
15911595
v0.v0cosPA() >= v0settingCosPALambda &&
15921596
v0.v0radius() >= v0settingRadiusLambda &&
15931597
std::abs(posDaughterTrack.tpcNSigmaPi()) <= nSigmaTPCPion &&
@@ -1633,13 +1637,13 @@ struct Lambdak0sflattenicity {
16331637
continue;
16341638
}
16351639

1636-
if (mcParticle.pdgCode() == 310) {
1640+
if (mcParticle.pdgCode() == PDG_t::kK0Short) {
16371641
rKzeroShort.fill(HIST("Generated_MCRecoCollCheck_INEL_K0Short"), mcParticle.pt(), flattenicity); // K0s
16381642
}
1639-
if (mcParticle.pdgCode() == 3122) {
1643+
if (mcParticle.pdgCode() == PDG_t::kLambda0) {
16401644
rLambda.fill(HIST("Generated_MCRecoCollCheck_INEL_Lambda"), mcParticle.pt(), flattenicity); // Lambda
16411645
}
1642-
if (mcParticle.pdgCode() == -3122) {
1646+
if (mcParticle.pdgCode() == PDG_t::kLambda0Bar) {
16431647
rAntiLambda.fill(HIST("Generated_MCRecoCollCheck_INEL_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
16441648
}
16451649
}
@@ -1700,25 +1704,25 @@ struct Lambdak0sflattenicity {
17001704
continue;
17011705
}
17021706

1703-
if (mcParticle.pdgCode() == 310) {
1707+
if (mcParticle.pdgCode() == PDG_t::kK0Short) {
17041708
rKzeroShort.fill(HIST("pGen_MCGenColl_INEL_K0Short"), mcParticle.pt(), flattenicity); // K0s
17051709
if (isINELgt0true) {
17061710
rKzeroShort.fill(HIST("pGen_MCGenColl_INELgt0_K0Short"), mcParticle.pt(), flattenicity); // K0s
17071711
}
17081712
}
1709-
if (mcParticle.pdgCode() == 3122) {
1713+
if (mcParticle.pdgCode() == PDG_t::kLambda0) {
17101714
rLambda.fill(HIST("pGen_MCGenColl_INEL_Lambda"), mcParticle.pt(), flattenicity); // Lambda
17111715
if (isINELgt0true) {
17121716
rLambda.fill(HIST("pGen_MCGenColl_INELgt0_Lambda"), mcParticle.pt(), flattenicity); // Lambda
17131717
}
17141718
}
1715-
if (mcParticle.pdgCode() == -3122) {
1719+
if (mcParticle.pdgCode() == PDG_t::kLambda0Bar) {
17161720
rAntiLambda.fill(HIST("pGen_MCGenColl_INEL_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
17171721
if (isINELgt0true) {
17181722
rAntiLambda.fill(HIST("pGen_MCGenColl_INELgt0_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
17191723
}
17201724
}
1721-
if (std::abs(mcParticle.pdgCode()) == 3312) {
1725+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kXiPlusBar) {
17221726
rXi.fill(HIST("pGen_MCGenColl_INEL_Xi"), mcParticle.pt(), flattenicity); // Xi
17231727
if (isINELgt0true) {
17241728
rXi.fill(HIST("pGen_MCGenColl_INELgt0_Xi"), mcParticle.pt(), flattenicity); // Xi
@@ -1765,25 +1769,25 @@ struct Lambdak0sflattenicity {
17651769
continue;
17661770
}
17671771

1768-
if (mcParticle.pdgCode() == 310) {
1772+
if (mcParticle.pdgCode() == PDG_t::kK0Short) {
17691773
rKzeroShort.fill(HIST("Generated_MCRecoColl_INEL_K0Short"), mcParticle.pt(), flattenicity); // K0s
17701774
if (recoCollIndexINELgt0 > 0) {
17711775
rKzeroShort.fill(HIST("Generated_MCRecoColl_INELgt0_K0Short"), mcParticle.pt(), flattenicity); // K0s
17721776
}
17731777
}
1774-
if (mcParticle.pdgCode() == 3122) {
1778+
if (mcParticle.pdgCode() == PDG_t::kLambda0) {
17751779
rLambda.fill(HIST("Generated_MCRecoColl_INEL_Lambda"), mcParticle.pt(), flattenicity); // Lambda
17761780
if (recoCollIndexINELgt0 > 0) {
17771781
rLambda.fill(HIST("Generated_MCRecoColl_INELgt0_Lambda"), mcParticle.pt(), flattenicity); // Lambda
17781782
}
17791783
}
1780-
if (mcParticle.pdgCode() == -3122) {
1784+
if (mcParticle.pdgCode() == PDG_t::kLambda0Bar) {
17811785
rAntiLambda.fill(HIST("Generated_MCRecoColl_INEL_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
17821786
if (recoCollIndexINELgt0 > 0) {
17831787
rAntiLambda.fill(HIST("Generated_MCRecoColl_INELgt0_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
17841788
}
17851789
}
1786-
if (std::abs(mcParticle.pdgCode()) == 3312) {
1790+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kXiPlusBar) {
17871791
rXi.fill(HIST("Generated_MCRecoColl_INEL_Xi"), mcParticle.pt(), flattenicity); // Xi
17881792
if (recoCollIndexINELgt0 > 0) {
17891793
rXi.fill(HIST("Generated_MCRecoColl_INELgt0_Xi"), mcParticle.pt(), flattenicity); // Xi
@@ -1823,25 +1827,25 @@ struct Lambdak0sflattenicity {
18231827
continue;
18241828
}
18251829

1826-
if (mcParticle.pdgCode() == 310) {
1830+
if (mcParticle.pdgCode() == PDG_t::kK0Short) {
18271831
rKzeroShort.fill(HIST("pGen_MCGenRecoColl_INEL_K0Short"), mcParticle.pt(), flattenicity); // K0s
18281832
if (recoCollIndexINELgt0 > 0) {
18291833
rKzeroShort.fill(HIST("pGen_MCGenRecoColl_INELgt0_K0Short"), mcParticle.pt(), flattenicity); // K0s
18301834
}
18311835
}
1832-
if (mcParticle.pdgCode() == 3122) {
1836+
if (mcParticle.pdgCode() == PDG_t::kLambda0) {
18331837
rLambda.fill(HIST("pGen_MCGenRecoColl_INEL_Lambda"), mcParticle.pt(), flattenicity); // Lambda
18341838
if (recoCollIndexINELgt0 > 0) {
18351839
rLambda.fill(HIST("pGen_MCGenRecoColl_INELgt0_Lambda"), mcParticle.pt(), flattenicity); // Lambda
18361840
}
18371841
}
1838-
if (mcParticle.pdgCode() == -3122) {
1842+
if (mcParticle.pdgCode() == PDG_t::kLambda0Bar) {
18391843
rAntiLambda.fill(HIST("pGen_MCGenRecoColl_INEL_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
18401844
if (recoCollIndexINELgt0 > 0) {
18411845
rAntiLambda.fill(HIST("pGen_MCGenRecoColl_INELgt0_AntiLambda"), mcParticle.pt(), flattenicity); // AntiLambda
18421846
}
18431847
}
1844-
if (std::abs(mcParticle.pdgCode()) == 3312) {
1848+
if (std::abs(mcParticle.pdgCode()) == PDG_t::kXiPlusBar) {
18451849
rXi.fill(HIST("pGen_MCGenRecoColl_INEL_Xi"), mcParticle.pt(), flattenicity); // Xi
18461850
if (recoCollIndexINELgt0 > 0) {
18471851
rXi.fill(HIST("pGen_MCGenRecoColl_INELgt0_Xi"), mcParticle.pt(), flattenicity); // Xi
@@ -1853,8 +1857,8 @@ struct Lambdak0sflattenicity {
18531857
// Cascade Analysis Starts here
18541858
using DauTracks = soa::Join<aod::TracksIU, aod::TrackSelection, aod::TracksExtra, aod::TracksDCA, aod::pidTPCFullPi, aod::pidTPCFullPr, aod::pidTPCFullKa, aod::pidTOFPi, aod::pidTOFPr, aod::pidTOFKa>;
18551859
using LabeledCascades = soa::Join<aod::CascDataExt, aod::McCascLabels>;
1856-
float ctauxiPDG = 4.91; // from PDG
1857-
float ctauomegaPDG = 2.461; // from PDG
1860+
// float ctauxiPDG = 4.91; // from PDG
1861+
// float ctauomegaPDG = 2.461; // from PDG
18581862

18591863
void processDataRun3Cascade(soa::Join<aod::Collisions, aod::EvSels,
18601864
aod::PVMults>::iterator const& collision,
@@ -1902,7 +1906,7 @@ struct Lambdak0sflattenicity {
19021906
casc.dcapostopv() > v0settingDCApostopv && casc.dcanegtopv() > v0settingDCAnegtopv && casc.dcabachtopv() > v0settingDCAbactopv && casc.dcaV0daughters() < v0settingDCAv0dau && dcav0pv > cascsettingDCAv0toPV &&
19031907
casc.dcacascdaughters() < cascsettingDCAv0bach && casc.bachBaryonDCAxyToPV() < cascsettingDCAxybaryonbach && cosPAcasc > cascsettingCosPAcascPV && cosPAv0 > cascsettingCosPAv0PV &&
19041908
casc.cascradius() > cascsettingcascradius && casc.v0radius() > cascsettingv0radius &&
1905-
std::abs(casc.yXi()) < cascsettingRapidity && ctauXi < ctauxiPDG * cascsettingproplifetime &&
1909+
std::abs(casc.yXi()) < cascsettingRapidity && ctauXi < 4.91 * cascsettingproplifetime &&
19061910
std::abs(casc.mLambda() - pdgmassLambda) < cascsettingMassRejectionLambdaXi && std::abs(casc.mOmega() - pdgmassOmega) > cascsettingMassRejectioOmegaXi) {
19071911

19081912
rXi.fill(HIST("hMassXiSelected"), massXi);
@@ -1969,7 +1973,7 @@ struct Lambdak0sflattenicity {
19691973
casc.dcapostopv() > v0settingDCApostopv && casc.dcanegtopv() > v0settingDCAnegtopv && casc.dcabachtopv() > v0settingDCAbactopv && casc.dcaV0daughters() < v0settingDCAv0dau && dcav0pv > cascsettingDCAv0toPV &&
19701974
casc.dcacascdaughters() < cascsettingDCAv0bach && casc.bachBaryonDCAxyToPV() < cascsettingDCAxybaryonbach && cosPAcasc > cascsettingCosPAcascPV && cosPAv0 > cascsettingCosPAv0PV &&
19711975
casc.cascradius() > cascsettingcascradius && casc.v0radius() > cascsettingv0radius &&
1972-
std::abs(casc.yXi()) < cascsettingRapidity && ctauXi < ctauxiPDG * cascsettingproplifetime &&
1976+
std::abs(casc.yXi()) < cascsettingRapidity && ctauXi < 4.91 * cascsettingproplifetime &&
19731977
std::abs(casc.mLambda() - pdgmassLambda) < cascsettingMassRejectionLambdaXi && std::abs(casc.mOmega() - pdgmassOmega) > cascsettingMassRejectioOmegaXi) {
19741978

19751979
rXi.fill(HIST("hMassXiSelected"), massXi);
@@ -1988,7 +1992,7 @@ struct Lambdak0sflattenicity {
19881992
rEventSelection.fill(HIST("hFlattenicity_Corr_Gen_vs_Rec"), flattenicityMCGen, flattenicity);
19891993

19901994
for (const auto& mcParticle : particlesInCollision) {
1991-
if (mcParticle.isPhysicalPrimary() && std::abs(mcParticle.y()) < 0.5f && std::abs(mcParticle.pdgCode()) == 3312) {
1995+
if (mcParticle.isPhysicalPrimary() && std::abs(mcParticle.y()) < 0.5f && std::abs(mcParticle.pdgCode()) == PDG_t::kXiPlusBar) {
19921996
rXi.fill(HIST("Generated_MCRecoCollCheck_INEL_Xi"), mcParticle.pt(), flattenicity); // K0s
19931997
}
19941998
}

0 commit comments

Comments
 (0)