Skip to content

Commit fb2bf5f

Browse files
committed
Please consider the following formatting changes
1 parent 2475d4d commit fb2bf5f

File tree

1 file changed

+143
-142
lines changed

1 file changed

+143
-142
lines changed

PWGLF/Tasks/Strangeness/strange-yield-pbpb.cxx

Lines changed: 143 additions & 142 deletions
Original file line numberDiff line numberDiff line change
@@ -63,7 +63,6 @@ using cascadeCandidatesMC = soa::Join<aod::CascCollRefs, aod::CascCores, aod::Ca
6363

6464
using cascMCCoresFull = soa::Join<aod::CascMCCores, aod::CascMCCollRefs>;
6565

66-
6766
using straCollisonFull = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraStamps>::iterator;
6867
using straCollisonFullMC = soa::Join<aod::StraCollisions, aod::StraCents, aod::StraEvSels, aod::StraStamps, aod::StraCollLabels>::iterator;
6968

@@ -1017,50 +1016,50 @@ struct strangeYieldPbPb {
10171016
};
10181017

10191018
const std::array<selectionCheck, 13> checks = {{
1020-
{true, true, 0.0}, // All collisions
1021-
{requireIsTriggerTVX, collision.selection_bit(aod::evsel::kIsTriggerTVX), 1.0}, // Triggered by FT0M
1022-
{true, fabs(collision.posZ()) <= 10.f, 2.0}, // Vertex-Z selected
1023-
{rejectITSROFBorder, collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder), 3.0}, // Not at ITS ROF border
1024-
{rejectTFBorder, collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder), 4.0}, // Not at TF border
1025-
{requireIsVertexITSTPC, collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC), 5.0}, // At least one ITS-TPC track
1026-
{requireIsGoodZvtxFT0VsPV, collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV), 6.0}, // PV position consistency
1027-
{requireIsVertexTOFmatched, collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched), 7.0}, // PV with TOF match
1028-
{requireIsVertexTRDmatched, collision.selection_bit(o2::aod::evsel::kIsVertexTRDmatched), 8.0}, // PV with TRD match
1029-
{rejectSameBunchPileup, collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup), 9.0}, // No same-bunch pileup
1030-
{requireNoCollInTimeRangeStd, collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard), 10.0}, // No collision within +-10 µs
1031-
{requireNoCollInTimeRangeNarrow, collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeNarrow), 11.0}, // No collision within +-4 µs
1032-
{minOccupancy > 0, collision.trackOccupancyInTimeRange() >= minOccupancy, 12.0}, // Above min occupancy
1019+
{true, true, 0.0}, // All collisions
1020+
{requireIsTriggerTVX, collision.selection_bit(aod::evsel::kIsTriggerTVX), 1.0}, // Triggered by FT0M
1021+
{true, fabs(collision.posZ()) <= 10.f, 2.0}, // Vertex-Z selected
1022+
{rejectITSROFBorder, collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder), 3.0}, // Not at ITS ROF border
1023+
{rejectTFBorder, collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder), 4.0}, // Not at TF border
1024+
{requireIsVertexITSTPC, collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC), 5.0}, // At least one ITS-TPC track
1025+
{requireIsGoodZvtxFT0VsPV, collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV), 6.0}, // PV position consistency
1026+
{requireIsVertexTOFmatched, collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched), 7.0}, // PV with TOF match
1027+
{requireIsVertexTRDmatched, collision.selection_bit(o2::aod::evsel::kIsVertexTRDmatched), 8.0}, // PV with TRD match
1028+
{rejectSameBunchPileup, collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup), 9.0}, // No same-bunch pileup
1029+
{requireNoCollInTimeRangeStd, collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeStandard), 10.0}, // No collision within +-10 µs
1030+
{requireNoCollInTimeRangeNarrow, collision.selection_bit(o2::aod::evsel::kNoCollInTimeRangeNarrow), 11.0}, // No collision within +-4 µs
1031+
{minOccupancy > 0, collision.trackOccupancyInTimeRange() >= minOccupancy, 12.0}, // Above min occupancy
10331032
}};
10341033

10351034
for (const auto& check : checks) {
1036-
if (check.selection && !check.condition) {
1037-
return false;
1038-
}
1039-
if (fillQA && check.selection) {
1040-
histos.fill(HIST("eventQA/hEventSelection"), check.qaBin);
1041-
}
1035+
if (check.selection && !check.condition) {
1036+
return false;
1037+
}
1038+
if (fillQA && check.selection) {
1039+
histos.fill(HIST("eventQA/hEventSelection"), check.qaBin);
1040+
}
10421041
}
10431042

10441043
// Additional checks
10451044
if (maxOccupancy > 0 && collision.trackOccupancyInTimeRange() > maxOccupancy) {
1046-
return false;
1045+
return false;
10471046
}
10481047
if (fillQA) {
1049-
histos.fill(HIST("eventQA/hEventSelection"), 13.0); // Below max occupancy
1048+
histos.fill(HIST("eventQA/hEventSelection"), 13.0); // Below max occupancy
10501049
}
10511050

10521051
if (studyUPConly && !collision.isUPC()) {
1053-
return false;
1052+
return false;
10541053
}
10551054
if (collision.isUPC() && fillQA) {
1056-
histos.fill(HIST("eventQA/hEventSelection"), 14.0); // UPC compatible
1055+
histos.fill(HIST("eventQA/hEventSelection"), 14.0); // UPC compatible
10571056
}
10581057

10591058
if (useUPCflag && collision.flags() < 1) {
1060-
return false;
1059+
return false;
10611060
}
10621061
if (collision.flags() >= 1 && fillQA) {
1063-
histos.fill(HIST("eventQA/hEventSelection"), 15.0); // UPC event
1062+
histos.fill(HIST("eventQA/hEventSelection"), 15.0); // UPC event
10641063
}
10651064

10661065
return true;
@@ -1227,8 +1226,8 @@ struct strangeYieldPbPb {
12271226
double bachDcaXYLimit = 0.0105f + 0.035f / pow(bachpt, 1.1f);
12281227

12291228
// TODO: separate xy and z //
1230-
if ((std::abs(casc.dcapostopv()) > posDcaXYLimit) &&
1231-
(std::abs(casc.dcanegtopv()) > negDcaXYLimit) &&
1229+
if ((std::abs(casc.dcapostopv()) > posDcaXYLimit) &&
1230+
(std::abs(casc.dcanegtopv()) > negDcaXYLimit) &&
12321231
(std::abs(casc.dcabachtopv()) > bachDcaXYLimit)) {
12331232
bitMap.set(selDauDCA);
12341233
}
@@ -1386,7 +1385,7 @@ struct strangeYieldPbPb {
13861385
double negDcaXYLimit = 0.0105f + 0.035f / pow(negpt, 1.1f);
13871386

13881387
// TODO: separate xy and z //
1389-
if ((std::abs(v0.dcapostopv()) > posDcaXYLimit) &&
1388+
if ((std::abs(v0.dcapostopv()) > posDcaXYLimit) &&
13901389
(std::abs(v0.dcanegtopv()) > negDcaXYLimit)) {
13911390
bitMap.set(selDauDCA);
13921391
}
@@ -1518,38 +1517,41 @@ struct strangeYieldPbPb {
15181517
}
15191518

15201519
template <typename TV0>
1521-
void computeV0MCAssociation(const TV0& v0, std::bitset<selNum>& bitMap)
1520+
void computeV0MCAssociation(const TV0& v0, std::bitset<selNum>& bitMap)
15221521
{
1523-
const int pdgPos = v0.pdgCodePositive();
1524-
const int pdgNeg = v0.pdgCodeNegative();
1525-
const int pdgV0 = v0.pdgCode();
1526-
const bool isPhysPrim = v0.isPhysicalPrimary();
1527-
1528-
const bool isPositiveProton = (pdgPos == 2212);
1529-
const bool isPositivePion = (pdgPos == 211) || (doTreatPiToMuon && pdgPos == -13);
1530-
const bool isNegativeProton = (pdgNeg == -2212);
1531-
const bool isNegativePion = (pdgNeg == -211) || (doTreatPiToMuon && pdgNeg == 13);
1532-
1533-
switch (pdgV0) {
1534-
case 310: // K0Short
1535-
if (isPositivePion && isNegativePion) {
1536-
bitMap.set(selConsiderK0Short);
1537-
if (isPhysPrim) bitMap.set(selPhysPrimK0Short);
1538-
}
1539-
break;
1540-
case 3122: // Lambda
1541-
if (isPositiveProton && isNegativePion) {
1542-
bitMap.set(selConsiderLambda);
1543-
if (isPhysPrim) bitMap.set(selPhysPrimLambda);
1544-
}
1545-
break;
1546-
case -3122: // AntiLambda
1547-
if (isPositivePion && isNegativeProton) {
1548-
bitMap.set(selConsiderAntiLambda);
1549-
if (isPhysPrim) bitMap.set(selPhysPrimAntiLambda);
1550-
}
1551-
break;
1552-
}
1522+
const int pdgPos = v0.pdgCodePositive();
1523+
const int pdgNeg = v0.pdgCodeNegative();
1524+
const int pdgV0 = v0.pdgCode();
1525+
const bool isPhysPrim = v0.isPhysicalPrimary();
1526+
1527+
const bool isPositiveProton = (pdgPos == 2212);
1528+
const bool isPositivePion = (pdgPos == 211) || (doTreatPiToMuon && pdgPos == -13);
1529+
const bool isNegativeProton = (pdgNeg == -2212);
1530+
const bool isNegativePion = (pdgNeg == -211) || (doTreatPiToMuon && pdgNeg == 13);
1531+
1532+
switch (pdgV0) {
1533+
case 310: // K0Short
1534+
if (isPositivePion && isNegativePion) {
1535+
bitMap.set(selConsiderK0Short);
1536+
if (isPhysPrim)
1537+
bitMap.set(selPhysPrimK0Short);
1538+
}
1539+
break;
1540+
case 3122: // Lambda
1541+
if (isPositiveProton && isNegativePion) {
1542+
bitMap.set(selConsiderLambda);
1543+
if (isPhysPrim)
1544+
bitMap.set(selPhysPrimLambda);
1545+
}
1546+
break;
1547+
case -3122: // AntiLambda
1548+
if (isPositivePion && isNegativeProton) {
1549+
bitMap.set(selConsiderAntiLambda);
1550+
if (isPhysPrim)
1551+
bitMap.set(selPhysPrimAntiLambda);
1552+
}
1553+
break;
1554+
}
15531555
}
15541556

15551557
template <typename TV0, typename TCollision>
@@ -1590,83 +1592,82 @@ struct strangeYieldPbPb {
15901592

15911593
PresliceUnsorted<straCollisonsFullMC> perMcCollision = aod::v0data::straMCCollisionId;
15921594

1593-
std::vector<int> getListOfRecoCollIndices(straMCCollisionsFull const& mcCollisions,
1595+
std::vector<int> getListOfRecoCollIndices(straMCCollisionsFull const& mcCollisions,
15941596
straCollisonsFullMC const& collisions)
15951597
{
1596-
std::vector<int> listBestCollisionIdx(mcCollisions.size(), -1);
1597-
1598-
for (auto const& mcCollision : mcCollisions) {
1599-
auto groupedCollisions = collisions.sliceBy(perMcCollision, mcCollision.globalIndex());
1600-
// Find the collision with the biggest nbr of PV contributors
1601-
// Follows what was done here: https://github.com/AliceO2Group/O2Physics/blob/master/Common/TableProducer/mcCollsExtra.cxx#L93
1602-
int biggestNContribs = -1;
1603-
int bestCollisionIndex = -1;
1604-
for (auto const& collision : groupedCollisions) {
1605-
if (biggestNContribs < collision.multPVTotalContributors()) {
1606-
biggestNContribs = collision.multPVTotalContributors();
1607-
bestCollisionIndex = collision.globalIndex();
1608-
}
1598+
std::vector<int> listBestCollisionIdx(mcCollisions.size(), -1);
1599+
1600+
for (auto const& mcCollision : mcCollisions) {
1601+
auto groupedCollisions = collisions.sliceBy(perMcCollision, mcCollision.globalIndex());
1602+
// Find the collision with the biggest nbr of PV contributors
1603+
// Follows what was done here: https://github.com/AliceO2Group/O2Physics/blob/master/Common/TableProducer/mcCollsExtra.cxx#L93
1604+
int biggestNContribs = -1;
1605+
int bestCollisionIndex = -1;
1606+
for (auto const& collision : groupedCollisions) {
1607+
if (biggestNContribs < collision.multPVTotalContributors()) {
1608+
biggestNContribs = collision.multPVTotalContributors();
1609+
bestCollisionIndex = collision.globalIndex();
16091610
}
1610-
listBestCollisionIdx[mcCollision.globalIndex()] = bestCollisionIndex;
16111611
}
1612+
listBestCollisionIdx[mcCollision.globalIndex()] = bestCollisionIndex;
1613+
}
16121614

1613-
return listBestCollisionIdx;
1615+
return listBestCollisionIdx;
16141616
}
16151617

16161618
void fillGenMCHistogramsQA(straMCCollisionsFull const& mcCollisions, straCollisonsFullMC const& collisions)
16171619
{
1618-
for (auto const& mcCollision : mcCollisions) {
1619-
histos.fill(HIST("eventQA/mc/hEventSelectionMC"), 0.0);
1620-
histos.fill(HIST("eventQA/mc/hMCNParticlesEta10"), mcCollision.multMCNParticlesEta10(), 0 /* all gen. events*/);
1621-
1622-
if (TMath::Abs(mcCollision.posZ()) > 10.f)
1623-
continue;
1624-
1625-
histos.fill(HIST("eventQA/mc/hEventSelectionMC"), 1.0);
1626-
1627-
// Group collisions by MC collision index
1628-
auto groupedCollisions = collisions.sliceBy(perMcCollision, mcCollision.globalIndex());
1629-
1630-
bool atLeastOne = false;
1631-
float centrality = -1.f;
1632-
int nCollisions = 0;
1633-
int biggestNContribs = -1;
1634-
int nTracksGlobal = -1;
1635-
1636-
// Find the max contributors and count accepted collisions
1637-
for (auto const& collision : groupedCollisions) {
1638-
if (!acceptEvent(collision, false)) {
1639-
continue;
1640-
}
1641-
1642-
int selGapSide = collision.isUPC() ? getGapSide(collision) : -1;
1643-
if (studyUPConly && (selGapSide < -0.5))
1644-
continue;
1645-
1646-
if (biggestNContribs < collision.multPVTotalContributors()) {
1647-
biggestNContribs = collision.multPVTotalContributors();
1648-
centrality = collision.centFT0C();
1649-
nTracksGlobal = collision.multNTracksGlobal();
1650-
}
1651-
1652-
++nCollisions;
1653-
atLeastOne = true;
1654-
}
1655-
1656-
// Fill histograms
1657-
histos.fill(HIST("eventQA/mc/hTracksGlobalVsNcoll_beforeEvSel"), nTracksGlobal, groupedCollisions.size());
1658-
histos.fill(HIST("eventQA/mc/hTracksGlobalVsNcoll_afterEvSel"), nTracksGlobal, nCollisions);
1659-
histos.fill(HIST("eventQA/mc/hTracksGlobalvsMCNParticlesEta10gen"), nTracksGlobal, mcCollision.multMCNParticlesEta10());
1660-
histos.fill(HIST("eventQA/mc/hTracksGlobalVsPVzMC"), nTracksGlobal, mcCollision.posZ());
1661-
histos.fill(HIST("eventQA/mc/hEventPVzMC"), mcCollision.posZ());
1662-
1663-
if (atLeastOne) {
1664-
histos.fill(HIST("eventQA/mc/hMCNParticlesEta10"), mcCollision.multMCNParticlesEta10(), 1 /* at least 1 rec. event*/);
1665-
histos.fill(HIST("eventQA/mc/hGenEventCentrality"), centrality);
1666-
}
1620+
for (auto const& mcCollision : mcCollisions) {
1621+
histos.fill(HIST("eventQA/mc/hEventSelectionMC"), 0.0);
1622+
histos.fill(HIST("eventQA/mc/hMCNParticlesEta10"), mcCollision.multMCNParticlesEta10(), 0 /* all gen. events*/);
1623+
1624+
if (TMath::Abs(mcCollision.posZ()) > 10.f)
1625+
continue;
1626+
1627+
histos.fill(HIST("eventQA/mc/hEventSelectionMC"), 1.0);
1628+
1629+
// Group collisions by MC collision index
1630+
auto groupedCollisions = collisions.sliceBy(perMcCollision, mcCollision.globalIndex());
1631+
1632+
bool atLeastOne = false;
1633+
float centrality = -1.f;
1634+
int nCollisions = 0;
1635+
int biggestNContribs = -1;
1636+
int nTracksGlobal = -1;
1637+
1638+
// Find the max contributors and count accepted collisions
1639+
for (auto const& collision : groupedCollisions) {
1640+
if (!acceptEvent(collision, false)) {
1641+
continue;
1642+
}
1643+
1644+
int selGapSide = collision.isUPC() ? getGapSide(collision) : -1;
1645+
if (studyUPConly && (selGapSide < -0.5))
1646+
continue;
1647+
1648+
if (biggestNContribs < collision.multPVTotalContributors()) {
1649+
biggestNContribs = collision.multPVTotalContributors();
1650+
centrality = collision.centFT0C();
1651+
nTracksGlobal = collision.multNTracksGlobal();
1652+
}
1653+
1654+
++nCollisions;
1655+
atLeastOne = true;
16671656
}
1668-
}
16691657

1658+
// Fill histograms
1659+
histos.fill(HIST("eventQA/mc/hTracksGlobalVsNcoll_beforeEvSel"), nTracksGlobal, groupedCollisions.size());
1660+
histos.fill(HIST("eventQA/mc/hTracksGlobalVsNcoll_afterEvSel"), nTracksGlobal, nCollisions);
1661+
histos.fill(HIST("eventQA/mc/hTracksGlobalvsMCNParticlesEta10gen"), nTracksGlobal, mcCollision.multMCNParticlesEta10());
1662+
histos.fill(HIST("eventQA/mc/hTracksGlobalVsPVzMC"), nTracksGlobal, mcCollision.posZ());
1663+
histos.fill(HIST("eventQA/mc/hEventPVzMC"), mcCollision.posZ());
1664+
1665+
if (atLeastOne) {
1666+
histos.fill(HIST("eventQA/mc/hMCNParticlesEta10"), mcCollision.multMCNParticlesEta10(), 1 /* at least 1 rec. event*/);
1667+
histos.fill(HIST("eventQA/mc/hGenEventCentrality"), centrality);
1668+
}
1669+
}
1670+
}
16701671

16711672
void processV0s(straCollisonFull const& collision, v0Candidates const& fullV0s, dauTracks const&)
16721673
{
@@ -1695,9 +1696,9 @@ struct strangeYieldPbPb {
16951696
}
16961697

16971698
void processV0sMC(straCollisonFullMC const& collision,
1698-
v0CandidatesMC const& fullV0s,
1699+
v0CandidatesMC const& fullV0s,
16991700
dauTracks const&,
1700-
aod::MotherMCParts const&,
1701+
aod::MotherMCParts const&,
17011702
straMCCollisionsFull const&,
17021703
v0MCCoresFull const&)
17031704
{
@@ -1716,11 +1717,11 @@ struct strangeYieldPbPb {
17161717
fillHistogramsQA(collision, selGapSide);
17171718

17181719
if (collision.has_straMCCollision()) {
1719-
const auto& mcCollision = collision.straMCCollision_as<straMCCollisionsFull>();
1720-
histos.fill(HIST("eventQA/mc/hNTracksGlobalvsMCNParticlesEta10rec"), collision.multNTracksGlobal(), mcCollision.multMCNParticlesEta10());
1721-
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vsMCNParticlesEta10rec"), collision.multNTracksPVeta1(), mcCollision.multMCNParticlesEta10());
1722-
histos.fill(HIST("eventQA/mc/hNTracksGlobalvstotalMultMCParticles"), collision.multNTracksGlobal(), mcCollision.totalMultMCParticles());
1723-
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vstotalMultMCParticles"), collision.multNTracksPVeta1(), mcCollision.totalMultMCParticles());
1720+
const auto& mcCollision = collision.straMCCollision_as<straMCCollisionsFull>();
1721+
histos.fill(HIST("eventQA/mc/hNTracksGlobalvsMCNParticlesEta10rec"), collision.multNTracksGlobal(), mcCollision.multMCNParticlesEta10());
1722+
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vsMCNParticlesEta10rec"), collision.multNTracksPVeta1(), mcCollision.multMCNParticlesEta10());
1723+
histos.fill(HIST("eventQA/mc/hNTracksGlobalvstotalMultMCParticles"), collision.multNTracksGlobal(), mcCollision.totalMultMCParticles());
1724+
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vstotalMultMCParticles"), collision.multNTracksPVeta1(), mcCollision.totalMultMCParticles());
17241725
}
17251726

17261727
for (auto& v0 : fullV0s) {
@@ -1744,7 +1745,7 @@ struct strangeYieldPbPb {
17441745
} // end v0 loop
17451746
}
17461747

1747-
void processCascadesMC(straCollisonFullMC const& collision,
1748+
void processCascadesMC(straCollisonFullMC const& collision,
17481749
cascadeCandidatesMC const& fullCascades,
17491750
dauTracks const&,
17501751
straMCCollisionsFull const&,
@@ -1765,17 +1766,17 @@ struct strangeYieldPbPb {
17651766
fillHistogramsQA(collision, selGapSide);
17661767

17671768
if (collision.has_straMCCollision()) {
1768-
const auto& mcCollision = collision.straMCCollision_as<straMCCollisionsFull>();
1769-
histos.fill(HIST("eventQA/mc/hNTracksGlobalvsMCNParticlesEta10rec"), collision.multNTracksGlobal(), mcCollision.multMCNParticlesEta10());
1770-
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vsMCNParticlesEta10rec"), collision.multNTracksPVeta1(), mcCollision.multMCNParticlesEta10());
1771-
histos.fill(HIST("eventQA/mc/hNTracksGlobalvstotalMultMCParticles"), collision.multNTracksGlobal(), mcCollision.totalMultMCParticles());
1772-
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vstotalMultMCParticles"), collision.multNTracksPVeta1(), mcCollision.totalMultMCParticles());
1769+
const auto& mcCollision = collision.straMCCollision_as<straMCCollisionsFull>();
1770+
histos.fill(HIST("eventQA/mc/hNTracksGlobalvsMCNParticlesEta10rec"), collision.multNTracksGlobal(), mcCollision.multMCNParticlesEta10());
1771+
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vsMCNParticlesEta10rec"), collision.multNTracksPVeta1(), mcCollision.multMCNParticlesEta10());
1772+
histos.fill(HIST("eventQA/mc/hNTracksGlobalvstotalMultMCParticles"), collision.multNTracksGlobal(), mcCollision.totalMultMCParticles());
1773+
histos.fill(HIST("eventQA/mc/hNTracksPVeta1vstotalMultMCParticles"), collision.multNTracksPVeta1(), mcCollision.totalMultMCParticles());
17731774
}
17741775
}
17751776

1776-
void processGenerated(straMCCollisionsFull const& mcCollisions,
1777-
v0MCCoresFull const& V0MCCores,
1778-
cascMCCoresFull const& CascMCCores,
1777+
void processGenerated(straMCCollisionsFull const& mcCollisions,
1778+
v0MCCoresFull const& V0MCCores,
1779+
cascMCCoresFull const& CascMCCores,
17791780
straCollisonsFullMC const& collisions)
17801781
{
17811782
fillGenMCHistogramsQA(mcCollisions, collisions);

0 commit comments

Comments
 (0)