Skip to content

Commit 28f6a1a

Browse files
1507183233715071832337
authored
Add files via upload
1 parent f800edc commit 28f6a1a

1 file changed

Lines changed: 30 additions & 31 deletions

File tree

‎PWGLF/Tasks/Strangeness/lambdaspincorrderived.cxx‎

Lines changed: 30 additions & 31 deletions
Original file line numberDiff line numberDiff line change
@@ -43,6 +43,7 @@
4343
#include <TVector3.h>
4444

4545
#include <algorithm>
46+
#include <array>
4647
#include <chrono>
4748
#include <cmath> // for std::abs
4849
#include <cstddef>
@@ -436,8 +437,6 @@ struct lambdaspincorrderived {
436437
LOGF(info, "sysCuts.size()=%zu nSysTotal=%d", sysCuts.size(), nSysTotal);
437438
const AxisSpec thnAxisSys{nSysTotal, -0.5f, float(nSysTotal) - 0.5f, "sysId"};
438439

439-
histos.add("hSysIdTest", "SysId Test", HistType::kTH1D, {thnAxisSys});
440-
441440
if (fillBasicQAHistos) {
442441
histos.add("hPtRadiusV0", "V0 QA;#it{p}_{T}^{V0} (GeV/#it{c});V0 decay radius (cm)", kTH2F, {{100, 0.0, 10.0}, {120, 0.0, 45.0}});
443442
histos.add("hPtYSame", "hPtYSame", kTH2F, {{100, 0.0, 10.0}, {200, -1.0, 1.0}});
@@ -1587,7 +1586,7 @@ struct lambdaspincorrderived {
15871586
std::vector<int> activePair;
15881587
activePair.reserve(activeSys.size());
15891588

1590-
for (int sysId : activeSys) {
1589+
for (const auto& sysId : activeSys) {
15911590
const auto& sc = sysCuts[sysId];
15921591
if (selectionV0Sys(v02, sc)) {
15931592
activePair.push_back(sysId);
@@ -1603,7 +1602,7 @@ struct lambdaspincorrderived {
16031602
lambda2 = ROOT::Math::PtEtaPhiMVector(v02.lambdaPt(), v02.lambdaEta(), v02.lambdaPhi(),
16041603
v02.lambdaMass());
16051604

1606-
for (int sysId : activePair) {
1605+
for (const auto& sysId : activePair) {
16071606

16081607
if (v0.v0Status() == 0 && v02.v0Status() == 0) {
16091608
fillHistogramsSys(0, 0, lambda, lambda2, proton, proton2, 0, 1.0, sysId);
@@ -1766,7 +1765,7 @@ struct lambdaspincorrderived {
17661765
auto nBins = colBinning.getAllBinsCount();
17671766
std::vector<std::deque<std::pair<int, AllTrackCandidates>>> eventPools(nBins);
17681767

1769-
for (auto& collision1 : collisions) {
1768+
for (const auto& collision1 : collisions) {
17701769
const int bin = colBinning.getBin(std::make_tuple(collision1.posz(), collision1.cent()));
17711770
if (bin < 0) {
17721771
continue;
@@ -1783,7 +1782,7 @@ struct lambdaspincorrderived {
17831782
continue;
17841783
}
17851784

1786-
for (auto& [t1, t2] : soa::combinations(o2::soa::CombinationsFullIndexPolicy(poolA, poolA))) {
1785+
for (const auto& [t1, t2] : soa::combinations(o2::soa::CombinationsFullIndexPolicy(poolA, poolA))) {
17871786
if (!selectionV0(t1) || !selectionV0(t2)) {
17881787
continue;
17891788
}
@@ -1818,7 +1817,7 @@ struct lambdaspincorrderived {
18181817
int nRepl1 = 0;
18191818
int nRepl2 = 0;
18201819

1821-
for (auto& tX : poolB) {
1820+
for (const auto& tX : poolB) {
18221821
if (!selectionV0(tX)) {
18231822
continue;
18241823
}
@@ -1874,10 +1873,10 @@ struct lambdaspincorrderived {
18741873
static_cast<float>(totalRepl2) * wBase);
18751874
}
18761875

1877-
for (auto& pv : usable) {
1876+
for (const auto& pv : usable) {
18781877
auto& poolB = *pv.pool;
18791878

1880-
for (auto& tX : poolB) {
1879+
for (const auto& tX : poolB) {
18811880
if (!selectionV0(tX)) {
18821881
continue;
18831882
}
@@ -2697,9 +2696,9 @@ struct lambdaspincorrderived {
26972696
auto collectFromBins = [&](const std::vector<int>& ptUseBins,
26982697
const std::vector<int>& etaUseBins,
26992698
const std::vector<int>& phiUseBins) {
2700-
for (int ptUse : ptUseBins) {
2701-
for (int etaUse : etaUseBins) {
2702-
for (int phiUse : phiUseBins) {
2699+
for (const auto& ptUse : ptUseBins) {
2700+
for (const auto& etaUse : etaUseBins) {
2701+
for (const auto& phiUse : phiUseBins) {
27032702
const auto& vec = buffer[linearKeyR(colBin, status, ptUse, etaUse, phiUse, mB, rB,
27042703
nStat, nPt, nEta, nPhi, nM, nR)];
27052704

@@ -2753,9 +2752,9 @@ struct lambdaspincorrderived {
27532752
collectNeighborBinsClamp(etaB, nEta, nN_eta, etaBinsN);
27542753
collectNeighborBinsPhi(phiB, nPhi, nN_phi, phiBinsN);
27552754

2756-
for (int ptUse : ptBinsN) {
2757-
for (int etaUse : etaBinsN) {
2758-
for (int phiUse : phiBinsN) {
2755+
for (const auto& ptUse : ptBinsN) {
2756+
for (const auto& etaUse : etaBinsN) {
2757+
for (const auto& phiUse : phiBinsN) {
27592758
if (ptUse == ptB && etaUse == etaB && phiUse == phiB)
27602759
continue;
27612760

@@ -3244,9 +3243,9 @@ struct lambdaspincorrderived {
32443243
auto collectFromBins = [&](const std::vector<int>& ptUseBins,
32453244
const std::vector<int>& etaUseBins,
32463245
const std::vector<int>& phiUseBins) {
3247-
for (int ptUse : ptUseBins) {
3248-
for (int etaUse : etaUseBins) {
3249-
for (int phiUse : phiUseBins) {
3246+
for (const auto& ptUse : ptUseBins) {
3247+
for (const auto& etaUse : etaUseBins) {
3248+
for (const auto& phiUse : phiUseBins) {
32503249
const auto& vec = buffer[linearKeyR(colBin, status, ptUse, etaUse, phiUse, mB, rB,
32513250
nStat, nPt, nEta, nPhi, nM, nR)];
32523251

@@ -3308,9 +3307,9 @@ struct lambdaspincorrderived {
33083307
collectNeighborBinsClamp(etaB, nEta, nN_eta, etaBinsN);
33093308
collectNeighborBinsPhi(phiB, nPhi, nN_phi, phiBinsN);
33103309

3311-
for (int ptUse : ptBinsN) {
3312-
for (int etaUse : etaBinsN) {
3313-
for (int phiUse : phiBinsN) {
3310+
for (const auto& ptUse : ptBinsN) {
3311+
for (const auto& etaUse : etaBinsN) {
3312+
for (const auto& phiUse : phiBinsN) {
33143313
if (ptUse == ptB && etaUse == etaB && phiUse == phiB)
33153314
continue;
33163315

@@ -3511,7 +3510,7 @@ struct lambdaspincorrderived {
35113510
auto tX = V0s.iteratorAt(static_cast<uint64_t>(m.rowIndex));
35123511

35133512
bool hasCommon = false;
3514-
for (int sysId : commonSys) {
3513+
for (const auto& sysId : commonSys) {
35153514
if (selectionV0Sys(tX, sysCuts[sysId])) {
35163515
hasCommon = true;
35173516
break;
@@ -3539,7 +3538,7 @@ struct lambdaspincorrderived {
35393538
for (auto const& m : matches2) {
35403539
auto tY = V0s.iteratorAt(static_cast<uint64_t>(m.rowIndex));
35413540
bool hasCommon = false;
3542-
for (int sysId : commonSys) {
3541+
for (const auto& sysId : commonSys) {
35433542
if (selectionV0Sys(tY, sysCuts[sysId])) {
35443543
hasCommon = true;
35453544
break;
@@ -3634,7 +3633,7 @@ struct lambdaspincorrderived {
36343633
const int s1 = tX.v0Status();
36353634
const int s2 = t2.v0Status();
36363635

3637-
for (int sysId : commonSys) {
3636+
for (const auto& sysId : commonSys) {
36383637
if (selectionV0Sys(tX, sysCuts[sysId])) {
36393638
if (s1 == 0 && s2 == 1) {
36403639
fillHistogramsSys(0, 1, lambda, lambda2, proton, proton2, 1, meWeight, sysId, 1, 1);
@@ -3686,7 +3685,7 @@ struct lambdaspincorrderived {
36863685
const int s1 = t1.v0Status();
36873686
const int s2 = tY.v0Status();
36883687

3689-
for (int sysId : commonSys) {
3688+
for (const auto& sysId : commonSys) {
36903689
if (selectionV0Sys(tY, sysCuts[sysId])) {
36913690
if (s1 == 0 && s2 == 1) {
36923691
fillHistogramsSys(0, 1, lambda, lambda2, proton, proton2, 1, meWeight, sysId, 2, 2);
@@ -3818,9 +3817,9 @@ struct lambdaspincorrderived {
38183817
auto collectFromBins = [&](const std::vector<int>& ptUseBins,
38193818
const std::vector<int>& etaUseBins,
38203819
const std::vector<int>& phiUseBins) {
3821-
for (int ptUse : ptUseBins) {
3822-
for (int etaUse : etaUseBins) {
3823-
for (int phiUse : phiUseBins) {
3820+
for (const auto& ptUse : ptUseBins) {
3821+
for (const auto& etaUse : etaUseBins) {
3822+
for (const auto& phiUse : phiUseBins) {
38243823
const auto& vec = buffer[linearKeyR(colBin, status, ptUse, etaUse, phiUse, mB, rB,
38253824
nStat, nPt, nEta, nPhi, nM, nR)];
38263825

@@ -3873,9 +3872,9 @@ struct lambdaspincorrderived {
38733872
collectNeighborBinsClamp(etaB, nEta, nN_eta, etaBinsN);
38743873
collectNeighborBinsPhi(phiB, nPhi, nN_phi, phiBinsN);
38753874

3876-
for (int ptUse : ptBinsN) {
3877-
for (int etaUse : etaBinsN) {
3878-
for (int phiUse : phiBinsN) {
3875+
for (const auto& ptUse : ptBinsN) {
3876+
for (const auto& etaUse : etaBinsN) {
3877+
for (const auto& phiUse : phiBinsN) {
38793878
if (ptUse == ptB && etaUse == etaB && phiUse == phiB) {
38803879
continue;
38813880
}

0 commit comments

Comments
 (0)