Skip to content

Commit 1af0d57

Browse files
Cristian-MoscatelliCristian Moscatelli
Cristian-Moscatelli
and
Cristian Moscatelli
authored
[PWGLF] Modified track selection criteria in deuteron analysis (#17742)
Co-authored-by: Cristian Moscatelli <cmoscate@gr3srv.ts.infn.it>
1 parent 4754e51 commit 1af0d57

3 files changed

Lines changed: 56 additions & 63 deletions

File tree

‎PWGLF/TableProducer/Nuspex/deuteronInTriggeredEvents.cxx‎

Lines changed: 1 addition & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -283,7 +283,6 @@ struct DeuteronInTriggeredEvents {
283283
Configurable<bool> rapidityToggle{"rapidityToggle", false, "If true, use rapidity cuts"};
284284
Configurable<float> tpcChi2ClusMax{"tpcChi2ClusMax", 4.f, "Max TPC Chi2 per cluster"};
285285
Configurable<int> tpcNCrossedRowsMin{"tpcNCrossedRowsMin", 70, "Minimum number of TPC crossed rows"};
286-
Configurable<float> tpcNCrossedRowsOverFindableMin{"tpcNCrossedRowsOverFindableMin", 0.8f, "Minimum ratio of crossed rows over findable clusters"};
287286
Configurable<int> tpcNClsMin{"tpcNClsMin", 80, "Minimum number of TPC clusters"};
288287
Configurable<float> tpcRigidityMin{"tpcRigidityMin", 0.5f, "Minimum TPC rigidity for tracks"};
289288
Configurable<LabeledArray<double>> tpcNSigmaMax{"tpcNSigmaMax", {nuclei::NSigmaTPCdefault[0], 5, 2, nuclei::names, nuclei::nSigmaConfigName}, "TPC nsigma selection for light nuclei"};
@@ -777,7 +776,6 @@ struct DeuteronInTriggeredEvents {
777776
track.itsNCls() < cfgTrackCut.itsNClusMin ||
778777
track.tpcNClsFound() < cfgTrackCut.tpcNClsMin ||
779778
track.tpcNClsCrossedRows() < cfgTrackCut.tpcNCrossedRowsMin ||
780-
track.tpcNClsCrossedRows() < cfgTrackCut.tpcNCrossedRowsOverFindableMin * track.tpcNClsFindable() ||
781779
track.tpcChi2NCl() > cfgTrackCut.tpcChi2ClusMax ||
782780
track.itsChi2NCl() > cfgTrackCut.itsChi2ClusMax) {
783781
continue;
@@ -993,7 +991,7 @@ struct DeuteronInTriggeredEvents {
993991
}
994992

995993
std::vector<bool> isReconstructed(particlesMC.size(), false);
996-
for (auto& c : nuclei::candidates) {
994+
for (auto& c : nuclei::candidates) { // o2-linter: disable=const-ref-in-for-loop (candidate is modified in loop)
997995
auto label = tracks.iteratorAt(c.globalIndex);
998996
if (label.mcParticleId() < -1 || label.mcParticleId() >= particlesMC.size()) {
999997
continue;

‎PWGLF/Tasks/Nuspex/CMakeLists.txt‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -202,7 +202,7 @@ o2physics_add_dpl_workflow(multiplicity-pt
202202
COMPONENT_NAME Analysis)
203203

204204
o2physics_add_dpl_workflow(deuteron-in-jets-trg-pt
205-
SOURCES DeuteronInJetsTrgPt.cxx
205+
SOURCES deuteronInJetsTrgPt.cxx
206206
PUBLIC_LINK_LIBRARIES O2Physics::AnalysisCore O2Physics::PWGJECore FastJet::FastJet FastJet::Contrib O2Physics::EventFilteringUtils
207207
COMPONENT_NAME Analysis)
208208

PWGLF/Tasks/Nuspex/DeuteronInJetsTrgPt.cxx renamed to PWGLF/Tasks/Nuspex/deuteronInJetsTrgPt.cxx

Lines changed: 54 additions & 59 deletions
Original file line numberDiff line numberDiff line change
@@ -9,7 +9,12 @@
99
// granted to it by virtue of its status as an Intergovernmental Organization
1010
// or submit itself to any jurisdiction.
1111
//
12-
// Task for analysing (anti)deuteron production in jets using pT-triggered data - update: 05-02-2026
12+
/// \file deuteronInJetsTrgPt.cxx
13+
/// \brief d and antid production in and out of jets using pTjet triggered data
14+
///
15+
/// \author Cristian Moscatelli <cristian.moscatelli@cern.ch>
16+
/// \since 02/2026
17+
// ================
1318
//
1419
// Executable : o2-analysis-lf-deuteron-in-jets-trg-pt
1520

@@ -73,58 +78,50 @@ struct DeuteronInJetsTrgPt {
7378

7479
// Setting default selection criteria to select tracks. May be changed when configuring the analysis.
7580
struct : o2::framework::ConfigurableGroup {
76-
std::string prefix{"cfgTrackCut"};
7781
// General specific
7882
Configurable<bool> requirePvContributor{"requirePvContributor", false, "Require that the track is a PV contributor"};
79-
Configurable<float> EtaMax{"EtaMax", 0.9f, "Max Eta for track acceptance"};
83+
Configurable<float> etaMax{"etaMax", 0.9f, "Max Eta for track acceptance"};
8084
Configurable<double> minPt{"minPt", 0.3, "Minimum pt of the tracks"};
8185
Configurable<double> maxDcaxy{"maxDcaxy", 0.05, "Maximum DCAxy"};
8286
Configurable<double> maxDcaz{"maxDcaz", 0.05, "Maximum DCAz"};
8387
// Part relative to ITS
84-
Configurable<int> ITSnClusMin{"ITSnClsMin", 6, "Minimum number of ITS clusters"};
85-
Configurable<float> ITSchi2ClusMax{"ITSchi2ClusMax", 36.f, "Max ITS Chi2 per cluster"};
88+
Configurable<int> itsNClusMin{"itsNClusMin", 6, "Minimum number of ITS clusters"};
89+
Configurable<float> itsChi2ClusMax{"itsChi2ClusMax", 36.f, "Max ITS Chi2 per cluster"};
8690
Configurable<bool> applyItsPid{"applyItsPid", false, "apply ITS PID"};
8791
Configurable<bool> setMCDefaultItsParams{"setMCDefaultItsParams", true, "Set MC default parameters for ITS PID"};
8892
Configurable<double> nSigmaItsMin{"nSigmaItsMin", -3.0, "nSigmaITS min"};
8993
Configurable<double> nSigmaItsMax{"nSigmaItsMax", +3.0, "nSigmaITS max"};
9094
Configurable<double> ptMaxItsPidProt{"ptMaxItsPidProt", 1.0, "maximum pt for ITS PID for protons"};
9195
Configurable<double> ptMaxItsPidDeut{"ptMaxItsPidDeut", 1.0, "maximum pt for ITS PID for deuterons"};
9296
// Part relative to TPC
93-
Configurable<int> TPCnClsMin{"TPCnClsMin", 100, "Minimum number of TPC clusters"};
94-
Configurable<float> TPCchi2ClusMin{"TPCchi2ClusMin", 0.f, "Min TPC Chi2 per cluster"};
95-
Configurable<float> TPCchi2ClusMax{"TPCchi2ClusMax", 4.f, "Max TPC Chi2 per cluster"};
96-
Configurable<int> TPCnCrossedRowsMin{"TPCnCrossedRowsMin", 100, "Minimum number of TPC crossed rows"};
97-
Configurable<double> Rtpc{"minRtpc", 0.8, "Minimum value of TPC crossed rows/TPC n cluster findable"};
98-
Configurable<float> TPCrigidityMin{"TPCrigidityMin", 0.3f, "Minimum TPC rigidity (p/Z) for track"};
99-
Configurable<double> minNsigmaTpc{"minNsigmaTpc", -3.0, "Minimum nsigma TPC"};
100-
Configurable<double> maxNsigmaTpc{"maxNsigmaTpc", +3.0, "Maximum nsigma TPC"};
101-
// Part relatuive to TOF
102-
Configurable<double> minNsigmaTof{"minNsigmaTof", -3.0, "Minimum nsigma TOF"};
103-
Configurable<double> maxNsigmaTof{"maxNsigmaTof", +3.5, "Maximum nsigma TOF"};
97+
Configurable<int> tpcNClsMin{"tpcNClsMin", 100, "Minimum number of TPC clusters"};
98+
Configurable<float> tpcChi2ClusMin{"tpcChi2ClusMin", 0.f, "Min TPC Chi2 per cluster"};
99+
Configurable<float> tpcChi2ClusMax{"tpcChi2ClusMax", 4.f, "Max TPC Chi2 per cluster"};
100+
Configurable<int> tpcNCrossedRowsMin{"tpcNCrossedRowsMin", 100, "Minimum number of TPC crossed rows"};
101+
Configurable<double> minNsigmaTpc{"minNsigmaTpc", -3.0, "Minimum nsigma TPC for TOF analysis"};
102+
Configurable<double> maxNsigmaTpc{"maxNsigmaTpc", +3.0, "Maximum nsigma TPC for TOF analysis"};
104103
} cfgTrackCut;
105104

106105
// Setting default selection criteria for events. May be changes when configuring the analysis.
107106
struct : o2::framework::ConfigurableGroup {
108-
std::string prefix{"cgfEventCut"};
109107
Configurable<double> zVtx{"zVtx", 10.0, "Maximum z vertex"};
110-
} cfgEvCut;
108+
} cfgEventCut;
111109

112110
// Skimmed data flag and list of active triggers for processing
113111
Configurable<bool> cfgSkimmedProcessing{"cfgSkimmedProcessing", false, "Skimmed dataset processing"};
114112
Configurable<std::string> triggerList{"triggerList", "fJetFullLowPt", "Trigger list"};
115113

116114
// Setting default selection criteria fr jet identification. May be changes when configuring the analysis.
117115
struct : o2::framework::ConfigurableGroup {
118-
std::string prefix{"cgfJetCut"};
119116
Configurable<double> minJetPt{"minJetPt", 10.0, "Minimum pt of the jet after bkg subtraction"};
120117
Configurable<double> rJet{"rJet", 0.3, "Jet parameter R"};
121118
Configurable<double> deltaEtaEdge{"deltaEtaEdge", 0.05, "eta gap from the edge"};
122119
} cfgJetCut;
123120

124121
// Setting the number of bins and min and max value for the nsigma distribution
125-
Configurable<int> cfgNbins{"Nbins", 120, "Number of pT-bins"};
126-
Configurable<double> cfgpt_min{"pt_min", 0.0, "Min pT value of pT-axis"};
127-
Configurable<double> cfgpt_max{"pt_max", 6.0, "Max pT value of pT-axis"};
122+
Configurable<int> cfgNbins{"cfgNbins", 120, "Number of pT-bins"};
123+
Configurable<double> cfgPtMin{"cfgPtMin", 0.0, "Min pT value of pT-axis"};
124+
Configurable<double> cfgPtMax{"cfgPtMax", 6.0, "Max pT value of pT-axis"};
128125

129126
// CCDB manager service for accessing condition data
130127
Service<o2::ccdb::BasicCCDBManager> ccdb;
@@ -159,45 +156,45 @@ struct DeuteronInJetsTrgPt {
159156
itsResponse.setMCDefaultParameters();
160157

161158
// Initialize random seed using high-resolution clock to ensure unique sequences across parallel Grid jobs
162-
auto time_seed = std::chrono::high_resolution_clock::now().time_since_epoch().count();
163-
mRand.SetSeed(time_seed);
159+
auto timeSeed = std::chrono::high_resolution_clock::now().time_since_epoch().count();
160+
mRand.SetSeed(timeSeed);
164161

165162
// Histrograms for real data
166163
if (doprocessData) {
167164
registryData.add("number_of_events_data", "number of events in data", HistType::kTH1F, {{4, 0, 4, "counter"}}); // Event counters
168-
registryData.add("settingData", "settingData", HistType::kTH2F, {{100, 0.0, 50.0, "min #it{p}^{jet}_{T} [GeV/#it{c}]"}, {20, 0.0, 1.0, "#it{R}_{jet}"}}); // Configuration
165+
registryData.add("settingData", "settingData", HistType::kTH2F, {{100, 0.0, 50.0, "min #it{p}^{jet}_{T} (GeV/#it{c})"}, {20, 0.0, 1.0, "#it{R}_{jet}"}}); // Configuration
169166
registryData.add("jetEffectiveAreaOverPiR2", "jet effective area / piR^2", HistType::kTH1F, {{2000, 0, 2, "Area/#piR^{2}"}}); // Jet effective area over piR^2
170167

171168
// Antiprotons
172-
registryData.add("antiproton_jet_tpc", "antiproton_jet_tpc", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
173-
registryData.add("antiproton_jet_tof", "antiproton_jet_tof", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
174-
registryData.add("antiproton_ue_tpc", "antiproton_ue_tpc", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
175-
registryData.add("antiproton_ue_tof", "antiproton_ue_tof", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
176-
registryData.add("antiproton_dca_jet", "antiproton_dca_jet", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {200, -1.0, 1.0, "DCA_{xy} [cm]"}});
177-
registryData.add("antiproton_dca_ue", "antiproton_dca_ue", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {200, -1.0, 1.0, "DCA_{xy} [cm]"}});
169+
registryData.add("antiproton_jet_tpc", "antiproton_jet_tpc", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
170+
registryData.add("antiproton_jet_tof", "antiproton_jet_tof", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
171+
registryData.add("antiproton_ue_tpc", "antiproton_ue_tpc", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
172+
registryData.add("antiproton_ue_tof", "antiproton_ue_tof", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
173+
registryData.add("antiproton_dca_jet", "antiproton_dca_jet", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {2000, -1.0, 1.0, "DCA_{xy} (cm)"}});
174+
registryData.add("antiproton_dca_ue", "antiproton_dca_ue", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {2000, -1.0, 1.0, "DCA_{xy} (cm)"}});
178175

179176
// protons
180-
registryData.add("proton_jet_tpc", "proton_jet_tpc", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
181-
registryData.add("proton_jet_tof", "proton_jet_tof", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
182-
registryData.add("proton_ue_tpc", "proton_ue_tpc", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
183-
registryData.add("proton_ue_tof", "proton_ue_tof", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
184-
registryData.add("proton_dca_jet", "proton_dca_jet", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {200, -1.0, 1.0, "DCA_{xy} [cm]"}});
185-
registryData.add("proton_dca_ue", "proton_dca_ue", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {200, -1.0, 1.0, "DCA_{xy} [cm]"}});
177+
registryData.add("proton_jet_tpc", "proton_jet_tpc", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
178+
registryData.add("proton_jet_tof", "proton_jet_tof", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
179+
registryData.add("proton_ue_tpc", "proton_ue_tpc", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
180+
registryData.add("proton_ue_tof", "proton_ue_tof", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
181+
registryData.add("proton_dca_jet", "proton_dca_jet", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {2000, -1.0, 1.0, "DCA_{xy} (cm)"}});
182+
registryData.add("proton_dca_ue", "proton_dca_ue", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {2000, -1.0, 1.0, "DCA_{xy} (cm)"}});
186183

187184
// Antideuterons
188-
registryData.add("antideuteron_jet_tpc", "antideuteron_jet_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
189-
registryData.add("antideuteron_jet_tof", "antideuteron_jet_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
190-
registryData.add("antideuteron_ue_tpc", "antideuteron_ue_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
191-
registryData.add("antideuteron_ue_tof", "antideuteron_ue_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
185+
registryData.add("antideuteron_jet_tpc", "antideuteron_jet_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
186+
registryData.add("antideuteron_jet_tof", "antideuteron_jet_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
187+
registryData.add("antideuteron_ue_tpc", "antideuteron_ue_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
188+
registryData.add("antideuteron_ue_tof", "antideuteron_ue_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
192189

193190
// Deuterons
194-
registryData.add("deuteron_jet_tpc", "deuteron_jet_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
195-
registryData.add("deuteron_jet_tof", "deuteron_jet_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
196-
registryData.add("deuteron_ue_tpc", "deuteron_ue_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TPC}"}});
197-
registryData.add("deuteron_ue_tof", "deuteron_ue_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgpt_min, 2 * cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{TOF}"}});
191+
registryData.add("deuteron_jet_tpc", "deuteron_jet_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
192+
registryData.add("deuteron_jet_tof", "deuteron_jet_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
193+
registryData.add("deuteron_ue_tpc", "deuteron_ue_tpc", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TPC}"}});
194+
registryData.add("deuteron_ue_tof", "deuteron_ue_tof", HistType::kTH2F, {{cfgNbins, 2 * cfgPtMin, 2 * cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{TOF}"}});
198195

199196
// nsigmaITS for antiproton candidates
200-
registryData.add("antiproton_nsigma_its_data", "antiproton_nsigma_its_data", HistType::kTH2F, {{cfgNbins, cfgpt_min, cfgpt_max, "#it{p}_{T} [GeV/#it{c}]"}, {400, -20.0, 20.0, "n#sigma_{ITS}"}});
197+
registryData.add("antiproton_nsigma_its_data", "antiproton_nsigma_its_data", HistType::kTH2F, {{cfgNbins, cfgPtMin, cfgPtMax, "#it{p}_{T} (GeV/#it{c})"}, {1000, -20.0, 20.0, "n#sigma_{ITS}"}});
201198
}
202199
}
203200

@@ -286,7 +283,7 @@ struct DeuteronInJetsTrgPt {
286283
static constexpr double DcazMaxTrack = 2.0;
287284

288285
// General part
289-
if (std::fabs(track.eta()) > cfgTrackCut.EtaMax)
286+
if (std::fabs(track.eta()) > cfgTrackCut.etaMax)
290287
return false;
291288
if (track.pt() < MinPtTrack)
292289
return false;
@@ -319,7 +316,7 @@ struct DeuteronInJetsTrgPt {
319316
// General part
320317
if (cfgTrackCut.requirePvContributor && !(track.isPVContributor()))
321318
return false; // Flag to check if the track contributed to the collision vertex fit
322-
if (std::fabs(track.eta()) > cfgTrackCut.EtaMax)
319+
if (std::fabs(track.eta()) > cfgTrackCut.etaMax)
323320
return false; // Eta
324321
if (track.pt() < cfgTrackCut.minPt)
325322
return false;
@@ -328,23 +325,21 @@ struct DeuteronInJetsTrgPt {
328325
return false; // Flag to check if track has a ITS match
329326
if ((!hasHitITS(track, 1)) && (!hasHitITS(track, 2)) && (!hasHitITS(track, 3)))
330327
return false; // Require IB hit
331-
if (track.itsNCls() < cfgTrackCut.ITSnClusMin)
328+
if (track.itsNCls() < cfgTrackCut.itsNClusMin)
332329
return false; // Minimum number of ITS cluster
333-
if (track.itsChi2NCl() > cfgTrackCut.ITSchi2ClusMax)
330+
if (track.itsChi2NCl() > cfgTrackCut.itsChi2ClusMax)
334331
return false; // Minimum chi2 per cluster in ITS
335332
// Part relative to TPC
336333
if (!track.hasTPC())
337334
return false; // Flag to check if track has a ITS match
338-
if (track.tpcNClsFound() < cfgTrackCut.TPCnClsMin)
335+
if (track.tpcNClsFound() < cfgTrackCut.tpcNClsMin)
339336
return false; // Minimum number of TPC cluster
340-
if (track.tpcNClsCrossedRows() < cfgTrackCut.TPCnCrossedRowsMin)
337+
if (track.tpcNClsCrossedRows() < cfgTrackCut.tpcNCrossedRowsMin)
341338
return false; // Minimum number of crossed rows in TPC
342-
if (track.tpcChi2NCl() < cfgTrackCut.TPCchi2ClusMin)
339+
if (track.tpcChi2NCl() < cfgTrackCut.tpcChi2ClusMin)
343340
return false; // Minimum chi2 per cluster in TPC
344-
if (track.tpcChi2NCl() > cfgTrackCut.TPCchi2ClusMax)
341+
if (track.tpcChi2NCl() > cfgTrackCut.tpcChi2ClusMax)
345342
return false; // Maximum chi2 per cluster in TPC
346-
if (track.tpcCrossedRowsOverFindableCls() < cfgTrackCut.Rtpc)
347-
return false; // R_{TPC} > 0.8
348343

349344
return true;
350345
}
@@ -359,13 +354,13 @@ struct DeuteronInJetsTrgPt {
359354
auto bc = collision.template bc_as<aod::BCsWithTimestamps>();
360355
initCCDB(bc);
361356

362-
// If skimmed processing is enabled, aplly Zorro trigger selection
357+
// If skimmed processing is enabled, apply Zorro trigger selection
363358
if (cfgSkimmedProcessing && !zorro.isSelected(collision.template bc_as<aod::BCsWithTimestamps>().globalBC()))
364359
return;
365360
registryData.fill(HIST("number_of_events_data"), 1.5);
366361

367362
// Apply standard event selection
368-
if (!collision.sel8() || std::fabs(collision.posZ()) >= cfgEvCut.zVtx)
363+
if (!collision.sel8() || std::fabs(collision.posZ()) >= cfgEventCut.zVtx)
369364
return;
370365
registryData.fill(HIST("number_of_events_data"), 2.5); // Save number of collisions that passed standard selections
371366

@@ -397,7 +392,7 @@ struct DeuteronInJetsTrgPt {
397392
bool isAtLeastOneJetSelected = false;
398393
for (auto const& jet : jets) {
399394

400-
if ((std::fabs(jet.eta()) + cfgJetCut.rJet) > (cfgTrackCut.EtaMax - cfgJetCut.deltaEtaEdge))
395+
if ((std::fabs(jet.eta()) + cfgJetCut.rJet) > (cfgTrackCut.etaMax - cfgJetCut.deltaEtaEdge))
401396
continue; // Jet must be fully contained in the acceptance
402397

403398
// Jet pt must be larger than threshold

0 commit comments

Comments
 (0)