|
49 | 49 | #include <TPDGCode.h> |
50 | 50 | #include <TProfile.h> |
51 | 51 |
|
| 52 | +#include <vector> |
| 53 | +#include <unordered_map> |
52 | 54 | #include <array> |
53 | 55 | #include <cmath> |
54 | 56 | #include <cstdlib> |
@@ -130,6 +132,7 @@ struct sigmaanalysis { |
130 | 132 | Configurable<float> mc_rapidityMin{"mc_rapidityMin", -0.5, "Min generated particle rapidity"}; |
131 | 133 | Configurable<float> mc_rapidityMax{"mc_rapidityMax", 0.5, "Max generated particle rapidity"}; |
132 | 134 | } genSelections; |
| 135 | + |
133 | 136 |
|
134 | 137 | // QA |
135 | 138 | Configurable<bool> fdoSigma0QA{"doSigma0QA", false, "if true, perform Sigma0 QA analysis. Only works with MC."}; |
@@ -549,7 +552,8 @@ struct sigmaanalysis { |
549 | 552 | histos.add("Gen/h2dGenAntiSigma0VsMultMC_RecoedEvt", "h2dGenAntiSigma0VsMultMC_RecoedEvt", kTH2D, {axisNch, axisPt}); |
550 | 553 | histos.add("Gen/h2dGenSigma0VsMultMC", "h2dGenSigma0VsMultMC", kTH2D, {axisNch, axisPt}); |
551 | 554 | histos.add("Gen/h2dGenAntiSigma0VsMultMC", "h2dGenAntiSigma0VsMultMC", kTH2D, {axisNch, axisPt}); |
552 | | - } |
| 555 | + |
| 556 | + } |
553 | 557 | if (doprocessPi0GeneratedRun3) { // Pi0 specific |
554 | 558 | histos.add("Gen/h2dGenPi0VsMultMC_RecoedEvt", "h2dGenPi0VsMultMC_RecoedEvt", kTH2D, {axisNch, axisPt}); |
555 | 559 | histos.add("Gen/h2dGenPi0", "h2dGenPi0", kTH2D, {axisCentrality, axisPt}); |
@@ -1221,15 +1225,15 @@ struct sigmaanalysis { |
1221 | 1225 | "TPCTOFPID", "DCADauToPV", "Mass"}; |
1222 | 1226 |
|
1223 | 1227 | if (PDGRequired == 22) { |
1224 | | - if constexpr (selection_index >= 0 && selection_index < (int)std::size(PhotonSelsLocal)) { |
| 1228 | + if constexpr (selection_index >= 0 && selection_index < static_cast<int>(std::size(PhotonSelsLocal))) { |
1225 | 1229 | histos.fill(HIST("Selection/Photon/hCandidateSel"), selection_index); |
1226 | 1230 | histos.fill(HIST("Selection/Photon/h2d") + HIST(PhotonSelsLocal[selection_index]), sigma.photonPt(), sigma.photonMass()); |
1227 | 1231 | histos.fill(HIST("Selection/Sigma0/h2dPhoton") + HIST(PhotonSelsLocal[selection_index]), sigma.pt(), sigma.sigma0Mass()); |
1228 | 1232 | } |
1229 | 1233 | } |
1230 | 1234 |
|
1231 | 1235 | if (PDGRequired == 3122) { |
1232 | | - if constexpr (selection_index >= 0 && selection_index < (int)std::size(LambdaSelsLocal)) { |
| 1236 | + if constexpr (selection_index >= 0 && selection_index < static_cast<int>(std::size(LambdaSelsLocal))) { |
1233 | 1237 | histos.fill(HIST("Selection/Lambda/hCandidateSel"), selection_index); |
1234 | 1238 | histos.fill(HIST("Selection/Lambda/h2d") + HIST(LambdaSelsLocal[selection_index]), sigma.lambdaPt(), sigma.lambdaMass()); |
1235 | 1239 | histos.fill(HIST("Selection/Sigma0/h2dLambda") + HIST(LambdaSelsLocal[selection_index]), sigma.pt(), sigma.sigma0Mass()); |
@@ -1263,9 +1267,8 @@ struct sigmaanalysis { |
1263 | 1267 | auto sigma0mc = fullSigma0s.rawIteratorAt(sigma0Index[mcid]); |
1264 | 1268 | histos.fill(HIST("Sigma0QA/hDuplicates"), NDuplicates); // how many times a mc sigma0 was reconstructed |
1265 | 1269 |
|
1266 | | - if (sigma0mc.isSigma0()) { |
1267 | | - } |
1268 | | - histos.fill(HIST("Sigma0QA/hSigma0Duplicates"), NDuplicates); // how many times a mc sigma0 was reconstructed |
| 1270 | + if (sigma0mc.isSigma0()) |
| 1271 | + histos.fill(HIST("Sigma0QA/hSigma0Duplicates"), NDuplicates); // how many times a mc sigma0 was reconstructed |
1269 | 1272 |
|
1270 | 1273 | if (sigma0mc.isAntiSigma0()) |
1271 | 1274 | histos.fill(HIST("Sigma0QA/hASigma0Duplicates"), NDuplicates); // how many times a mc sigma0 was reconstructed |
@@ -1730,3 +1733,4 @@ WorkflowSpec defineDataProcessing(ConfigContext const& cfgc) |
1730 | 1733 | { |
1731 | 1734 | return WorkflowSpec{adaptAnalysisTask<sigmaanalysis>(cfgc)}; |
1732 | 1735 | } |
| 1736 | + |
0 commit comments