Skip to content

[PWGLF] Please consider the following formatting changes to #12443 #79

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
alcaliva merged 1 commit into from
Aug 5, 2025
Merged

[PWGLF] Please consider the following formatting changes to #12443 #79

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
110 changes: 45 additions & 65 deletions PWGLF/Tasks/Nuspex/antinucleiInJets.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -115,7 +115,7 @@
Configurable<bool> applyRandomEventRejection{"applyRandomEventRejection", false, "reject some events for syst"};
Configurable<double> rejectionPercentage{"rejectionPercentage", 3.0, "percentage of events to reject"};
Configurable<int> nSyst{"nSyst", 50, "number of systematic variations"};

// Track quality, kinematic, and PID selection parameters
Configurable<bool> requirePvContributor{"requirePvContributor", false, "require that the track is a PV contributor"};
Configurable<bool> applyItsPid{"applyItsPid", true, "apply ITS PID"};
Expand Down Expand Up @@ -259,7 +259,7 @@
registryMC.add("antiproton_rec_tof_jet", "antiproton_rec_tof_jet", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tpc_ue", "antiproton_rec_tpc_ue", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_rec_tof_ue", "antiproton_rec_tof_ue", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});

// Fraction of primary antiprotons
registryMC.add("antiproton_prim_jet", "antiproton_prim_jet", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
registryMC.add("antiproton_incl_jet", "antiproton_incl_jet", HistType::kTH1F, {{nbins, min, max, "#it{p}_{T} (GeV/#it{c})"}});
Expand Down Expand Up @@ -430,14 +430,14 @@
template <typename JetTrack>
bool passedTrackSelectionForJetReconstruction(const JetTrack& track)
{
static constexpr int minTpcCr = 70;

Check failure on line 433 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double maxChi2Tpc = 4.0;

Check failure on line 434 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double maxChi2Its = 36.0;

Check failure on line 435 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double minPtTrack = 0.1;

Check failure on line 436 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double dcaxyMaxTrackPar0 = 0.0105;

Check failure on line 437 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double dcaxyMaxTrackPar1 = 0.035;

Check failure on line 438 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double dcaxyMaxTrackPar2 = 1.1;

Check failure on line 439 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double dcazMaxTrack = 2.0;

Check failure on line 440 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".

if (!track.hasITS())
return false;
Expand Down Expand Up @@ -498,43 +498,37 @@
4, 3, 6, 3, 7, 5, 4, 6, 5, 7,
6, 5, 3, 5, 4, 3, 6, 6, 4, 7,
3, 4, 3, 5, 7, 6, 6, 4, 3, 5,
4, 7, 3, 6, 4, 5, 6, 3, 7, 5
};
4, 7, 3, 6, 4, 5, 6, 3, 7, 5};
static std::vector<int> minTpcNcrossedRowsSyst = {
72, 99, 87, 66, 83, 71, 78, 74, 84, 94,
89, 85, 66, 80, 95, 74, 76, 77, 67, 88,
61, 100, 97, 78, 81, 91, 64, 63, 93, 69,
87, 99, 65, 62, 96, 60, 75, 92, 82, 89,
98, 70, 90, 85, 73, 79, 68, 86, 97, 61
};
98, 70, 90, 85, 73, 79, 68, 86, 97, 61};
static std::vector<double> maxChiSquareTpcSyst = {
4.28, 4.81, 4.43, 4.02, 3.38, 3.58, 3.11, 4.17, 3.51, 4.53,
4.90, 3.07, 3.20, 4.86, 4.62, 3.91, 3.98, 4.38, 3.66, 3.84,
3.03, 3.14, 4.96, 4.07, 4.75, 4.32, 3.31, 3.78, 4.11, 3.23,
3.87, 3.70, 4.99, 4.48, 4.69, 4.25, 3.93, 3.45, 4.58, 3.35,
3.18, 3.60, 4.21, 3.75, 4.64, 4.35, 3.26, 3.42, 4.15, 3.09
};
3.18, 3.60, 4.21, 3.75, 4.64, 4.35, 3.26, 3.42, 4.15, 3.09};
static std::vector<double> maxChiSquareItsSyst = {
42.84, 48.66, 39.27, 34.09, 43.73, 36.98, 30.23, 49.11, 37.67, 35.10,
44.55, 46.79, 38.92, 40.66, 47.14, 33.46, 30.88, 41.32, 45.90, 39.68,
31.42, 32.71, 43.17, 36.04, 49.80, 33.95, 31.89, 38.37, 48.08, 35.87,
47.61, 44.02, 32.15, 46.21, 34.75, 40.17, 37.14, 30.55, 45.42, 42.30,
41.79, 33.21, 39.12, 47.98, 36.52, 31.58, 49.44, 38.02, 35.56, 43.49
};
41.79, 33.21, 39.12, 47.98, 36.52, 31.58, 49.44, 38.02, 35.56, 43.49};
static std::vector<double> minEtaSyst = {
-0.804, -0.844, -0.751, -0.784, -0.819, -0.823, -0.768, -0.781, -0.845, -0.787,
-0.758, -0.828, -0.776, -0.842, -0.808, -0.763, -0.849, -0.770, -0.799, -0.754,
-0.825, -0.847, -0.806, -0.783, -0.796, -0.835, -0.777, -0.752, -0.838, -0.813,
-0.785, -0.802, -0.795, -0.846, -0.780, -0.829, -0.817, -0.773, -0.765, -0.789,
-0.800, -0.839, -0.758, -0.820, -0.833, -0.762, -0.792, -0.809, -0.827, -0.751
};
-0.800, -0.839, -0.758, -0.820, -0.833, -0.762, -0.792, -0.809, -0.827, -0.751};
static std::vector<double> maxEtaSyst = {
0.804, 0.844, 0.751, 0.784, 0.819, 0.823, 0.768, 0.781, 0.845, 0.787,
0.758, 0.828, 0.776, 0.842, 0.808, 0.763, 0.849, 0.770, 0.799, 0.754,
0.825, 0.847, 0.806, 0.783, 0.796, 0.835, 0.777, 0.752, 0.838, 0.813,
0.785, 0.802, 0.795, 0.846, 0.780, 0.829, 0.817, 0.773, 0.765, 0.789,
0.800, 0.839, 0.758, 0.820, 0.833, 0.762, 0.792, 0.809, 0.827, 0.751
};
0.800, 0.839, 0.758, 0.820, 0.833, 0.762, 0.792, 0.809, 0.827, 0.751};

// Track Selection
if (requirePvContributor && !(track.isPVContributor()))
Expand All @@ -555,7 +549,7 @@
return false;
if (track.pt() < minPt)
return false;

return true;
}

Expand All @@ -567,8 +561,8 @@
double nsigmaTPCPr = track.tpcNSigmaPr();
double nsigmaTOFPr = track.tofNSigmaPr();
double pt = track.pt();
static constexpr double ptThreshold = 0.5;

Check failure on line 564 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".
static constexpr double nsigmaMaxPr = 2.0;

Check failure on line 565 in PWGLF/Tasks/Nuspex/antinucleiInJets.cxx

View workflow job for this annotation

GitHub Actions / O2 linter

[name/constexpr-constant] 

Use UpperCamelCase for names of constexpr constants. Names of special constants may be prefixed with "k".

if (pt < ptThreshold && std::fabs(nsigmaTPCPr) < nsigmaMaxPr)
return true;
Expand Down Expand Up @@ -596,7 +590,7 @@
// Retrieve the bunch crossing information with timestamps from the collision
auto bc = collision.template bc_as<aod::BCsWithTimestamps>();
initCCDB(bc);

// If skimmed processing is enabled, apply Zorro trigger selection
if (cfgSkimmedProcessing && !zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC())) {
return;
Expand Down Expand Up @@ -942,7 +936,7 @@
fastjet::PseudoJet jetMinusBkg = backgroundSub.doRhoAreaSub(jetForSub, rhoPerp, rhoMPerp);
if (jetMinusBkg.pt() < minJetPt)
continue;

// Apply area cut if required
double normalizedJetArea = jet.area() / (PI * rJet * rJet);
if (applyAreaCut && normalizedJetArea > maxNormalizedJetArea)
Expand Down Expand Up @@ -1334,7 +1328,7 @@
for (const auto& particle : jetConstituents) {
if (particle.user_index() != PDG_t::kProtonBar)
continue;

if (particle.eta() < minEta || particle.eta() > maxEta)
continue;

Expand Down Expand Up @@ -1374,7 +1368,7 @@
// Reject tracks that lie outside the maxConeRadius from both UE axes
if (deltaRUe1 > maxConeRadius && deltaRUe2 > maxConeRadius)
continue;

// Select antiprotons based on PDG
if (particle.pdgCode() != PDG_t::kProtonBar)
continue;
Expand Down Expand Up @@ -1405,7 +1399,7 @@

// Event counter: after event selection
registryMC.fill(HIST("recEvents"), 1.5);

// Reject events near the ITS Read-Out Frame border
if (rejectITSROFBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))
continue;
Expand Down Expand Up @@ -1435,7 +1429,7 @@
if (requireIsVertexTOFmatched && !collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched))
continue;
registryMC.fill(HIST("recEvents"), 7.5);

// Loop over reconstructed tracks
int id(-1);
std::vector<fastjet::PseudoJet> fjParticles;
Expand Down Expand Up @@ -1569,7 +1563,7 @@
if (!track.has_mcParticle())
continue;
const auto mcparticle = track.mcParticle();

// Antiproton selection
if (track.sign() > 0 || mcparticle.pdgCode() != PDG_t::kProtonBar)
continue;
Expand All @@ -1580,7 +1574,7 @@
double pt = track.pt();
double dcaxy = track.dcaXY();
double dcaz = track.dcaZ();

// Fill DCA templates
if (mcparticle.pdgCode() == PDG_t::kProtonBar && std::fabs(dcaz) < maxDcaz) {
if (mcparticle.isPhysicalPrimary()) {
Expand Down Expand Up @@ -1621,11 +1615,11 @@

// Fill inclusive antiproton spectrum
registryMC.fill(HIST("antiproton_incl_ue"), pt);

// Select physical primary antiprotons
if (!mcparticle.isPhysicalPrimary())
continue;

// Fill antiproton spectrum for physical primaries
registryMC.fill(HIST("antiproton_prim_ue"), pt);

Expand Down Expand Up @@ -1653,39 +1647,39 @@

// Apply standard event selection
if (!collision.sel8() || std::fabs(collision.posZ()) > zVtx)
return;
return;

// Event counter: after event selection
registryData.fill(HIST("number_of_events_data_syst"), 1.5);

// Reject events near the ITS Read-Out Frame border
if (rejectITSROFBorder && !collision.selection_bit(o2::aod::evsel::kNoITSROFrameBorder))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 2.5);

// Reject events at the Time Frame border
if (rejectTFBorder && !collision.selection_bit(o2::aod::evsel::kNoTimeFrameBorder))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 3.5);

// Require at least one ITS-TPC matched track
if (requireVtxITSTPC && !collision.selection_bit(o2::aod::evsel::kIsVertexITSTPC))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 4.5);

// Reject events with same-bunch pileup
if (rejectSameBunchPileup && !collision.selection_bit(o2::aod::evsel::kNoSameBunchPileup))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 5.5);

// Require consistent FT0 vs PV z-vertex
if (requireIsGoodZvtxFT0VsPV && !collision.selection_bit(o2::aod::evsel::kIsGoodZvtxFT0vsPV))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 6.5);

// Require TOF match for at least one vertex track
if (requireIsVertexTOFmatched && !collision.selection_bit(o2::aod::evsel::kIsVertexTOFmatched))
return;
return;
registryData.fill(HIST("number_of_events_data_syst"), 7.5);

// Cut settings
Expand All @@ -1694,43 +1688,37 @@
0.079, 0.043, 0.067, 0.059, 0.032, 0.070, 0.048, 0.077, 0.062, 0.034,
0.057, 0.055, 0.073, 0.038, 0.050, 0.075, 0.041, 0.061, 0.033, 0.069,
0.035, 0.044, 0.076, 0.049, 0.037, 0.054, 0.072, 0.046, 0.058, 0.040,
0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030
};
0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030};
static std::vector<double> maxDcazSyst = {
0.064, 0.047, 0.032, 0.076, 0.039, 0.058, 0.043, 0.069, 0.050, 0.035,
0.074, 0.061, 0.045, 0.033, 0.068, 0.055, 0.037, 0.071, 0.042, 0.053,
0.077, 0.038, 0.065, 0.049, 0.036, 0.059, 0.044, 0.067, 0.041, 0.034,
0.073, 0.052, 0.040, 0.063, 0.046, 0.031, 0.070, 0.054, 0.037, 0.062,
0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056
};
0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056};
static std::vector<double> nSigmaItsMinSyst = {
-2.91, -2.77, -3.03, -3.40, -2.69, -3.28, -2.96, -3.11, -3.36, -3.14,
-2.99, -2.75, -3.17, -2.64, -2.72, -3.44, -2.87, -2.95, -3.00, -2.66,
-2.93, -3.31, -3.05, -3.12, -3.21, -3.01, -2.89, -2.73, -3.26, -2.97,
-2.81, -3.33, -2.68, -3.30, -2.78, -3.39, -2.84, -3.45, -2.92, -3.15,
-3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82
};
-3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82};
static std::vector<double> nSigmaItsMaxSyst = {
2.91, 2.77, 3.03, 3.40, 2.69, 3.28, 2.96, 3.11, 3.36, 3.14,
2.99, 2.75, 3.17, 2.64, 2.72, 3.44, 2.87, 2.95, 3.00, 2.66,
2.93, 3.31, 3.05, 3.12, 3.21, 3.01, 2.89, 2.73, 3.26, 2.97,
2.81, 3.33, 2.68, 3.30, 2.78, 3.39, 2.84, 3.45, 2.92, 3.15,
3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82
};
3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82};
static std::vector<double> minNsigmaTpcSyst = {
-3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
-3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
-3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
-3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
};
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
static std::vector<double> maxNsigmaTpcSyst = {
3.18, 2.86, 3.12, 2.91, 3.49, 2.57, 3.33, 2.98, 3.46, 2.70,
3.01, 2.65, 3.27, 3.40, 2.81, 3.10, 2.55, 3.22, 3.07, 2.77,
3.35, 2.68, 3.43, 2.88, 3.04, 2.53, 3.30, 2.79, 3.15, 2.66,
3.41, 2.75, 3.26, 2.61, 3.09, 2.54, 3.36, 2.95, 3.20, 2.58,
3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90
};
3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90};

// Initialize ITS PID Response object
o2::aod::ITSResponse itsResponse;
Expand All @@ -1743,12 +1731,12 @@
continue;

// Loop over different cut settings
for (int isyst = 0 ; isyst < nSyst ; isyst++) {
for (int isyst = 0; isyst < nSyst; isyst++) {

// Apply track selection
if (!passedTrackSelectionSyst(track,isyst))
if (!passedTrackSelectionSyst(track, isyst))
continue;

// Define variables
double nsigmaTPCPr = track.tpcNSigmaPr();
double nsigmaTOFPr = track.tofNSigmaPr();
Expand Down Expand Up @@ -1807,57 +1795,49 @@
0.079, 0.043, 0.067, 0.059, 0.032, 0.070, 0.048, 0.077, 0.062, 0.034,
0.057, 0.055, 0.073, 0.038, 0.050, 0.075, 0.041, 0.061, 0.033, 0.069,
0.035, 0.044, 0.076, 0.049, 0.037, 0.054, 0.072, 0.046, 0.058, 0.040,
0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030
};
0.068, 0.042, 0.056, 0.039, 0.047, 0.065, 0.051, 0.053, 0.063, 0.030};
static std::vector<double> maxDcazSyst = {
0.064, 0.047, 0.032, 0.076, 0.039, 0.058, 0.043, 0.069, 0.050, 0.035,
0.074, 0.061, 0.045, 0.033, 0.068, 0.055, 0.037, 0.071, 0.042, 0.053,
0.077, 0.038, 0.065, 0.049, 0.036, 0.059, 0.044, 0.067, 0.041, 0.034,
0.073, 0.052, 0.040, 0.063, 0.046, 0.031, 0.070, 0.054, 0.037, 0.062,
0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056
};
0.048, 0.035, 0.075, 0.051, 0.039, 0.066, 0.043, 0.060, 0.032, 0.056};
static std::vector<double> nSigmaItsMinSyst = {
-2.91, -2.77, -3.03, -3.40, -2.69, -3.28, -2.96, -3.11, -3.36, -3.14,
-2.99, -2.75, -3.17, -2.64, -2.72, -3.44, -2.87, -2.95, -3.00, -2.66,
-2.93, -3.31, -3.05, -3.12, -3.21, -3.01, -2.89, -2.73, -3.26, -2.97,
-2.81, -3.33, -2.68, -3.30, -2.78, -3.39, -2.84, -3.45, -2.92, -3.15,
-3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82
};
-3.22, -2.58, -3.07, -2.86, -3.10, -2.76, -2.94, -3.25, -3.04, -2.82};
static std::vector<double> nSigmaItsMaxSyst = {
2.91, 2.77, 3.03, 3.40, 2.69, 3.28, 2.96, 3.11, 3.36, 3.14,
2.99, 2.75, 3.17, 2.64, 2.72, 3.44, 2.87, 2.95, 3.00, 2.66,
2.93, 3.31, 3.05, 3.12, 3.21, 3.01, 2.89, 2.73, 3.26, 2.97,
2.81, 3.33, 2.68, 3.30, 2.78, 3.39, 2.84, 3.45, 2.92, 3.15,
3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82
};
3.22, 2.58, 3.07, 2.86, 3.10, 2.76, 2.94, 3.25, 3.04, 2.82};
static std::vector<double> minNsigmaTpcSyst = {
-3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
-3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
-3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
-3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
};
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
static std::vector<double> maxNsigmaTpcSyst = {
3.18, 2.86, 3.12, 2.91, 3.49, 2.57, 3.33, 2.98, 3.46, 2.70,
3.01, 2.65, 3.27, 3.40, 2.81, 3.10, 2.55, 3.22, 3.07, 2.77,
3.35, 2.68, 3.43, 2.88, 3.04, 2.53, 3.30, 2.79, 3.15, 2.66,
3.41, 2.75, 3.26, 2.61, 3.09, 2.54, 3.36, 2.95, 3.20, 2.58,
3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90
};
3.44, 2.83, 3.11, 2.62, 3.28, 2.69, 3.23, 2.73, 3.39, 2.90};
static std::vector<double> minNsigmaTofSyst = {
-3.18, -2.86, -3.12, -2.91, -3.49, -2.57, -3.33, -2.98, -3.46, -2.70,
-3.01, -2.65, -3.27, -3.40, -2.81, -3.10, -2.55, -3.22, -3.07, -2.77,
-3.35, -2.68, -3.43, -2.88, -3.04, -2.53, -3.30, -2.79, -3.15, -2.66,
-3.41, -2.75, -3.26, -2.61, -3.09, -2.54, -3.36, -2.95, -3.20, -2.58,
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90
};
-3.44, -2.83, -3.11, -2.62, -3.28, -2.69, -3.23, -2.73, -3.39, -2.90};
static std::vector<double> maxNsigmaTofSyst = {
3.94, 3.62, 3.83, 3.15, 3.23, 3.49, 3.10, 3.78, 3.54, 3.36,
3.91, 3.84, 3.72, 3.00, 3.63, 3.13, 3.68, 3.40, 3.97, 3.01,
3.74, 3.25, 3.89, 3.08, 3.30, 3.48, 3.59, 3.16, 3.47, 3.31,
3.92, 3.03, 3.43, 3.24, 3.11, 3.86, 3.60, 3.07, 3.21, 3.98,
3.14, 3.69, 3.56, 3.12, 3.28, 3.46, 3.34, 3.39, 3.05, 3.76
};
3.14, 3.69, 3.56, 3.12, 3.28, 3.46, 3.34, 3.39, 3.05, 3.76};

// Loop over generated collisions
for (const auto& collision : genCollisions) {
Expand Down Expand Up @@ -1939,10 +1919,10 @@
const auto mcparticle = track.mcParticle();

// Loop over different cut settings
for (int isyst = 0 ; isyst < nSyst ; isyst++) {
for (int isyst = 0; isyst < nSyst; isyst++) {

// Apply track selection
if (!passedTrackSelectionSyst(track,isyst))
if (!passedTrackSelectionSyst(track, isyst))
continue;

// Define variables
Expand Down
Loading