Skip to content

Commit 7b9b978

Browse files
spucillospucillo
committed
Fix build errors PR 13987
1 parent 5fcbe44 commit 7b9b978

3 files changed

Lines changed: 176 additions & 13 deletions

File tree

‎PWGLF/DataModel/LFInJets.h‎

Lines changed: 13 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -11,6 +11,7 @@
1111
///
1212
/// \brief Derived Data table for LF in jets analysis
1313
/// \author Francesca Ercolessi (francesca.ercolessi@cern.ch)
14+
/// \author Sara Pucillo (sara.pucillo@cern.ch)
1415

1516
#ifndef PWGLF_DATAMODEL_LFINJETS_H_
1617
#define PWGLF_DATAMODEL_LFINJETS_H_
@@ -23,8 +24,8 @@ namespace o2::aod
2324
namespace lfinjets
2425
{
2526

26-
DECLARE_SOA_COLUMN(Sign, sign, int);
2727
DECLARE_SOA_COLUMN(Pt, pt, float);
28+
DECLARE_SOA_COLUMN(Sign, sign, int);
2829
DECLARE_SOA_COLUMN(MassLambda, masslambda, float);
2930
DECLARE_SOA_COLUMN(MassAntiLambda, massantilambda, float);
3031
DECLARE_SOA_COLUMN(MassK0Short, massk0short, float);
@@ -37,6 +38,10 @@ DECLARE_SOA_COLUMN(NTPCSigmaNegPr, ntpcsigmanegpr, float);
3738
DECLARE_SOA_COLUMN(NTPCSigmaPosPr, ntpcsigmapospr, float);
3839
DECLARE_SOA_COLUMN(NTPCSigmaNegPi, ntpcsigmanegpi, float);
3940
DECLARE_SOA_COLUMN(NTPCSigmaPosPi, ntpcsigmapospi, float);
41+
DECLARE_SOA_COLUMN(NTOFSigmaNegPr, ntofsigmanegpr, float);
42+
DECLARE_SOA_COLUMN(NTOFSigmaPosPr, ntofsigmapospr, float);
43+
DECLARE_SOA_COLUMN(NTOFSigmaNegPi, ntofsigmanegpi, float);
44+
DECLARE_SOA_COLUMN(NTOFSigmaPosPi, ntofsigmapospi, float);
4045
DECLARE_SOA_COLUMN(MultFT0M, multft0m, float);
4146
DECLARE_SOA_COLUMN(V0PosTPCCrossedRows, v0postpcCrossedRows, float);
4247
DECLARE_SOA_COLUMN(V0NegTPCCrossedRows, v0negtpcCrossedRows, float);
@@ -46,14 +51,15 @@ DECLARE_SOA_COLUMN(V0PosTPCChi2, v0posTPCChi2, float);
4651
DECLARE_SOA_COLUMN(V0PosITSlayers, v0posITSlayers, int);
4752
DECLARE_SOA_COLUMN(MassXi, massxi, float);
4853
DECLARE_SOA_COLUMN(MassOmega, massomega, float);
49-
DECLARE_SOA_COLUMN(MassLambdaDau, masslambdadau, float);
5054
DECLARE_SOA_COLUMN(CascRadius, cascradius, float);
5155
DECLARE_SOA_COLUMN(CascCosPA, casccospa, float);
5256
DECLARE_SOA_COLUMN(DCABachToPV, dcabachtopv, float);
5357
DECLARE_SOA_COLUMN(DCACascDaughters, dcacascdaughters, float);
5458
DECLARE_SOA_COLUMN(DCAV0ToPV, dcav0topv, float);
5559
DECLARE_SOA_COLUMN(NTPCSigmaBachPi, ntpcsigmabachpi, float);
5660
DECLARE_SOA_COLUMN(NTPCSigmaBachKa, ntpcsigmabachka, float);
61+
DECLARE_SOA_COLUMN(NTOFSigmaBachPi, ntofsigmabachpi, float);
62+
DECLARE_SOA_COLUMN(NTOFSigmaBachKa, ntofsigmabachka, float);
5763
DECLARE_SOA_COLUMN(BachTPCCrossedRows, bachtpcCrossedRows, float);
5864
DECLARE_SOA_COLUMN(BachTPCChi2, bachTPCChi2, float);
5965
DECLARE_SOA_COLUMN(BachITSlayers, bachITSlayers, int);
@@ -67,17 +73,20 @@ DECLARE_SOA_TABLE(V0InJets, "AOD", "V0INJETS",
6773
lfinjets::V0Radius, lfinjets::V0CosPA, lfinjets::V0DCAPosToPV,
6874
lfinjets::V0DCANegToPV, lfinjets::V0DCAV0Daughters,
6975
lfinjets::NTPCSigmaNegPr, lfinjets::NTPCSigmaPosPr, lfinjets::NTPCSigmaNegPi, lfinjets::NTPCSigmaPosPi,
76+
lfinjets::NTOFSigmaNegPr, lfinjets::NTOFSigmaPosPr, lfinjets::NTOFSigmaNegPi, lfinjets::NTOFSigmaPosPi,
7077
lfinjets::MultFT0M, lfinjets::V0PosTPCCrossedRows, lfinjets::V0NegTPCCrossedRows,
7178
lfinjets::V0NegTPCChi2, lfinjets::V0NegITSlayers, lfinjets::V0PosTPCChi2, lfinjets::V0PosITSlayers,
7279
lfinjets::IsUE, lfinjets::IsJC);
7380

7481
DECLARE_SOA_TABLE(CascInJets, "AOD", "CASCINJETS",
75-
lfinjets::Pt, lfinjets::Sign, lfinjets::MassXi, lfinjets::MassOmega, lfinjets::MassLambdaDau,
82+
lfinjets::Pt, lfinjets::Sign, lfinjets::MassXi, lfinjets::MassOmega, lfinjets::MassLambda,
7683
lfinjets::CascRadius, lfinjets::CascCosPA, lfinjets::V0Radius, lfinjets::V0CosPA,
7784
lfinjets::V0DCAPosToPV, lfinjets::V0DCANegToPV, lfinjets::DCABachToPV,
7885
lfinjets::DCACascDaughters, lfinjets::V0DCAV0Daughters, lfinjets::DCAV0ToPV,
7986
lfinjets::NTPCSigmaNegPr, lfinjets::NTPCSigmaPosPr, lfinjets::NTPCSigmaNegPi, lfinjets::NTPCSigmaPosPi,
80-
lfinjets::NTPCSigmaBachPi, lfinjets::NTPCSigmaBachKa, lfinjets::MultFT0M,
87+
lfinjets::NTPCSigmaBachPi, lfinjets::NTPCSigmaBachKa,
88+
lfinjets::NTOFSigmaNegPr, lfinjets::NTOFSigmaPosPr, lfinjets::NTOFSigmaNegPi, lfinjets::NTOFSigmaPosPi,
89+
lfinjets::NTOFSigmaBachPi, lfinjets::NTOFSigmaBachKa, lfinjets::MultFT0M,
8190
lfinjets::V0PosTPCCrossedRows, lfinjets::V0NegTPCCrossedRows, lfinjets::BachTPCCrossedRows,
8291
lfinjets::V0NegTPCChi2, lfinjets::V0NegITSlayers, lfinjets::V0PosTPCChi2, lfinjets::V0PosITSlayers,
8392
lfinjets::BachTPCChi2, lfinjets::BachITSlayers,

‎PWGLF/TableProducer/Strangeness/lfinjets.cxx‎

Lines changed: 10 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -1149,6 +1149,10 @@ struct LFInJets {
11491149
pos.tpcNSigmaPr(),
11501150
neg.tpcNSigmaPi(),
11511151
pos.tpcNSigmaPi(),
1152+
neg.tofNSigmaPr(),
1153+
pos.tofNSigmaPr(),
1154+
neg.tofNSigmaPi(),
1155+
pos.tofNSigmaPi(),
11521156
collision.centFT0M(),
11531157
pos.tpcNClsCrossedRows(),
11541158
neg.tpcNClsCrossedRows(),
@@ -1229,6 +1233,12 @@ struct LFInJets {
12291233
pos.tpcNSigmaPi(),
12301234
bach.tpcNSigmaPi(),
12311235
bach.tpcNSigmaKa(),
1236+
neg.tofNSigmaPr(),
1237+
pos.tofNSigmaPr(),
1238+
neg.tofNSigmaPi(),
1239+
pos.tofNSigmaPi(),
1240+
bach.tofNSigmaPi(),
1241+
bach.tofNSigmaKa(),
12321242
collision.centFT0M(),
12331243
pos.tpcNClsCrossedRows(),
12341244
neg.tpcNClsCrossedRows(),

‎PWGLF/Tasks/Strangeness/strangenessInJets.cxx‎

Lines changed: 153 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -428,17 +428,17 @@ struct StrangenessInJets {
428428
return false;
429429

430430
// Constants for identifying heavy-flavor (charm and bottom) content from PDG codes
431-
static constexpr int kCharmQuark = 4;
432-
static constexpr int kBottomQuark = 5;
433-
static constexpr int hundreds = 100;
434-
static constexpr int thousands = 1000;
431+
static constexpr int CharmQuark = 4;
432+
static constexpr int BottomQuark = 5;
433+
static constexpr int Hundreds = 100;
434+
static constexpr int Thousands = 1000;
435435

436436
// Check if particle is from heavy-flavor decay
437437
bool fromHF = false;
438438
if (particle.has_mothers()) {
439439
auto mother = mcParticles.iteratorAt(particle.mothersIds()[0]);
440440
int motherPdg = std::abs(mother.pdgCode());
441-
fromHF = (motherPdg / hundreds == kCharmQuark || motherPdg / hundreds == kBottomQuark || motherPdg / thousands == kCharmQuark || motherPdg / thousands == kBottomQuark);
441+
fromHF = (motherPdg / Hundreds == CharmQuark || motherPdg / Hundreds == BottomQuark || motherPdg / Thousands == CharmQuark || motherPdg / Thousands == BottomQuark);
442442
}
443443

444444
// Select only physical primary particles or from heavy-flavor
@@ -1255,8 +1255,8 @@ struct StrangenessInJets {
12551255
continue;
12561256

12571257
// Build 4-momentum assuming charged pion mass
1258-
static constexpr float kMassPionChargedSquared = o2::constants::physics::MassPionCharged * o2::constants::physics::MassPionCharged;
1259-
const double energy = std::sqrt(particle.p() * particle.p() + kMassPionChargedSquared);
1258+
static constexpr float MassPionChargedSquared = o2::constants::physics::MassPionCharged * o2::constants::physics::MassPionCharged;
1259+
const double energy = std::sqrt(particle.p() * particle.p() + MassPionChargedSquared);
12601260
fastjet::PseudoJet fourMomentum(particle.px(), particle.py(), particle.pz(), energy);
12611261
fourMomentum.set_user_index(particle.pdgCode());
12621262
fjParticles.emplace_back(fourMomentum);
@@ -1773,7 +1773,7 @@ struct StrangenessInJets {
17731773
// Postprocessing
17741774
void processDerivedAnalysis(aod::V0InJets const& v0s, aod::CascInJets const& cascades)
17751775
{
1776-
for (auto& v0 : v0s) {
1776+
for (const auto& v0 : v0s) {
17771777

17781778
if (v0.v0negITSlayers() < minITSnCls || v0.v0posITSlayers() < minITSnCls)
17791779
continue;
@@ -1791,6 +1791,45 @@ struct StrangenessInJets {
17911791
continue;
17921792
if (std::fabs(v0.v0dcanegtopv()) < dcanegtoPVmin)
17931793
continue;
1794+
// PID selections (TPC) -- K0s
1795+
if (v0.ntpcsigmapospi() < nsigmaTPCmin || v0.ntpcsigmapospi() > nsigmaTPCmax)
1796+
continue;
1797+
if (v0.ntpcsigmanegpi() < nsigmaTPCmin || v0.ntpcsigmanegpi() > nsigmaTPCmax)
1798+
continue;
1799+
1800+
// PID selections (TOF) -- K0s
1801+
if (requireTOF) {
1802+
if (v0.ntofsigmapospi() < nsigmaTOFmin || v0.ntofsigmapospi() > nsigmaTOFmax)
1803+
continue;
1804+
if (v0.ntofsigmanegpi() < nsigmaTOFmin || v0.ntofsigmanegpi() > nsigmaTOFmax)
1805+
continue;
1806+
}
1807+
// PID selections (TPC): positive track = proton, negative track = pion -- Lam
1808+
if (v0.ntpcsigmapospr() < nsigmaTPCmin || v0.ntpcsigmapospr() > nsigmaTPCmax)
1809+
continue;
1810+
if (v0.ntpcsigmanegpi() < nsigmaTPCmin || v0.ntpcsigmanegpi() > nsigmaTPCmax)
1811+
continue;
1812+
1813+
// PID selections (TOF): positive track = proton, negative track = pion -- Lam
1814+
if (requireTOF) {
1815+
if (v0.ntofsigmapospr() < nsigmaTOFmin || v0.ntofsigmapospr() > nsigmaTOFmax)
1816+
continue;
1817+
if (v0.ntofsigmanegpi() < nsigmaTOFmin || v0.ntofsigmanegpi() > nsigmaTOFmax)
1818+
continue;
1819+
}
1820+
// PID selections (TPC): negative track = proton, positive track = pion --- ALam
1821+
if (v0.ntpcsigmapospi() < nsigmaTPCmin || v0.ntpcsigmapospi() > nsigmaTPCmax)
1822+
continue;
1823+
if (v0.ntpcsigmanegpr() < nsigmaTPCmin || v0.ntpcsigmanegpr() > nsigmaTPCmax)
1824+
continue;
1825+
1826+
// PID selections (TOF): negative track = proton, positive track = pion --- ALam
1827+
if (requireTOF) {
1828+
if (v0.ntofsigmapospi() < nsigmaTOFmin || v0.ntofsigmapospi() > nsigmaTOFmax)
1829+
continue;
1830+
if (v0.ntofsigmanegpr() < nsigmaTOFmin || v0.ntofsigmanegpr() > nsigmaTOFmax)
1831+
continue;
1832+
}
17941833

17951834
if (v0.isUE()) {
17961835
registryData.fill(HIST("K0s_in_ue"), v0.multft0m(), v0.pt(), v0.massk0short());
@@ -1803,7 +1842,7 @@ struct StrangenessInJets {
18031842
}
18041843
}
18051844

1806-
for (auto& casc : cascades) {
1845+
for (const auto& casc : cascades) {
18071846

18081847
if (casc.v0negITSlayers() < minITSnCls || casc.v0posITSlayers() < minITSnCls || casc.bachITSlayers() < minITSnCls)
18091848
continue;
@@ -1828,6 +1867,111 @@ struct StrangenessInJets {
18281867
continue;
18291868
if (std::fabs(casc.dcabachtopv()) < dcabachtopvMin)
18301869
continue;
1870+
if (std::fabs(casc.dcav0topv()) < dcaV0topvMin)
1871+
continue;
1872+
if (std::fabs(casc.dcacascdaughters()) > dcaCascDaughtersMax)
1873+
continue;
1874+
//Xi
1875+
// Xi+ selection (Xi+ -> antiL + pi+)
1876+
if (casc.sign() > 0) {
1877+
// PID selections (TPC)
1878+
if (casc.ntpcsigmanegpr() < nsigmaTPCmin || casc.ntpcsigmanegpr() > nsigmaTPCmax)
1879+
continue;
1880+
if (casc.ntpcsigmapospi() < nsigmaTPCmin || casc.ntpcsigmapospi() > nsigmaTPCmax)
1881+
continue;
1882+
1883+
// PID selections (TOF)
1884+
if (requireTOF) {
1885+
if (casc.ntofsigmanegpr() < nsigmaTOFmin || casc.ntofsigmanegpr() > nsigmaTOFmax)
1886+
continue;
1887+
if (casc.ntofsigmapospi() < nsigmaTOFmin || casc.ntofsigmapospi() > nsigmaTOFmax)
1888+
continue;
1889+
}
1890+
}
1891+
// Xi- selection (Xi- -> L + pi-)
1892+
if (casc.sign() < 0) {
1893+
// PID selections (TPC)
1894+
if (casc.ntpcsigmapospr() < nsigmaTPCmin || casc.ntpcsigmapospr() > nsigmaTPCmax)
1895+
continue;
1896+
if (casc.ntpcsigmanegpi() < nsigmaTPCmin || casc.ntpcsigmanegpi() > nsigmaTPCmax)
1897+
continue;
1898+
1899+
// PID selections (TOF)
1900+
if (requireTOF) {
1901+
if (casc.ntofsigmapospr() < nsigmaTOFmin || casc.ntofsigmapospr() > nsigmaTOFmax)
1902+
continue;
1903+
if (casc.ntofsigmanegpi() < nsigmaTOFmin || casc.ntofsigmanegpi() > nsigmaTOFmax)
1904+
continue;
1905+
}
1906+
}
1907+
1908+
// PID selection on bachelor
1909+
if (casc.ntpcsigmabachpi() < nsigmaTPCmin || casc.ntpcsigmabachpi() > nsigmaTPCmax)
1910+
continue;
1911+
1912+
// PID selections (TOF)
1913+
if (requireTOF) {
1914+
if (casc.ntofsigmabachpi() < nsigmaTOFmin || casc.ntofsigmabachpi() > nsigmaTOFmax)
1915+
continue;
1916+
}
1917+
// V0 mass window
1918+
if (std::fabs(casc.masslambda() - o2::constants::physics::MassLambda0) > deltaMassLambda)
1919+
continue;
1920+
// Reject candidates compatible with Omega
1921+
if (std::fabs(casc.massomega() - o2::constants::physics::MassOmegaMinus) < deltaMassOmega)
1922+
continue;
1923+
1924+
//Omega
1925+
// Omega+ selection (Omega+ -> antiL + K+)
1926+
if (casc.sign() > 0) {
1927+
// PID selections (TPC)
1928+
if (casc.ntpcsigmanegpr() < nsigmaTPCmin || casc.ntpcsigmanegpr() > nsigmaTPCmax)
1929+
continue;
1930+
if (casc.ntpcsigmapospi() < nsigmaTPCmin || casc.ntpcsigmapospi() > nsigmaTPCmax)
1931+
continue;
1932+
1933+
// PID selections (TOF)
1934+
if (requireTOF) {
1935+
if (casc.ntofsigmanegpr() < nsigmaTOFmin || casc.ntofsigmanegpr() > nsigmaTOFmax)
1936+
continue;
1937+
if (casc.ntofsigmapospi() < nsigmaTOFmin || casc.ntofsigmapospi() > nsigmaTOFmax)
1938+
continue;
1939+
}
1940+
}
1941+
1942+
// Omega- selection (Omega- -> L + K-)
1943+
if (casc.sign() < 0) {
1944+
// PID selections (TPC)
1945+
if (casc.ntpcsigmapospr() < nsigmaTPCmin || casc.ntpcsigmapospr() > nsigmaTPCmax)
1946+
continue;
1947+
if (casc.ntpcsigmanegpi() < nsigmaTPCmin || casc.ntpcsigmanegpi() > nsigmaTPCmax)
1948+
continue;
1949+
1950+
// PID selections (TOF)
1951+
if (requireTOF) {
1952+
if (casc.ntofsigmapospr() < nsigmaTOFmin || casc.ntofsigmapospr() > nsigmaTOFmax)
1953+
continue;
1954+
if (casc.ntofsigmanegpi() < nsigmaTOFmin || casc.ntofsigmanegpi() > nsigmaTOFmax)
1955+
continue;
1956+
}
1957+
}
1958+
1959+
// PID selection on bachelor
1960+
if (casc.ntpcsigmabachka() < nsigmaTPCmin || casc.ntpcsigmabachka() > nsigmaTPCmax)
1961+
continue;
1962+
1963+
// PID selections (TOF)
1964+
if (requireTOF) {
1965+
if (casc.ntofsigmabachka() < nsigmaTOFmin || casc.ntofsigmabachka() > nsigmaTOFmax)
1966+
continue;
1967+
}
1968+
// V0 mass window
1969+
if (std::fabs(casc.masslambda() - o2::constants::physics::MassLambda0) > deltaMassLambda)
1970+
continue;
1971+
1972+
// Reject candidates compatible with Xi
1973+
if (std::fabs(casc.massxi() - o2::constants::physics::MassXiMinus) < deltaMassXi)
1974+
continue;
18311975

18321976
if (casc.isUE()) {
18331977
if (casc.sign() < 0) {

0 commit comments

Comments
 (0)