Skip to content

Commit 0f3e933

Browse files
joachimckhjoachimckh
authored
[PWGJE] add h2 ft0c/ft0a qnred, psi44+psi42 (#16831)
1 parent e883f61 commit 0f3e933

1 file changed

Lines changed: 94 additions & 42 deletions

File tree

‎PWGJE/Tasks/jetSpectraEseTask.cxx‎

Lines changed: 94 additions & 42 deletions
Original file line numberDiff line numberDiff line change
@@ -100,11 +100,8 @@ struct JetSpectraEseTask {
100100
Configurable<float> cfgChi2PrTPCcls{"cfgChi2PrTPCcls", 2.5, "cut for chi2 per TPC cluster for tracks"};
101101
Configurable<float> cfgChi2PrITScls{"cfgChi2PrITScls", 36, "cut for chi2 per ITS cluster for tracks"};
102102

103-
Configurable<float> trackPtMinRhoPhi{"trackPtMinRhoPhi", 0.2, "minimum pT acceptance for tracks used in rho(phi) calculation"};
104-
Configurable<float> trackPtMaxRhoPhi{"trackPtMaxRhoPhi", 5.0, "maximum pT acceptance for tracks used in rho(phi) calculation"};
105-
106-
Configurable<float> jetEtaMin{"jetEtaMin", -0.7, "minimum jet pseudorapidity"};
107-
Configurable<float> jetEtaMax{"jetEtaMax", 0.7, "maximum jet pseudorapidity"};
103+
Configurable<std::vector<float>> trackPtRhoPhi{"trackPtRhoPhi", {0.2, 5.0}, "pT range for tracks used in rho(phi) calculation"};
104+
Configurable<std::vector<float>> cfgJetEta{"cfgJetEta", {-0.7, 0.7}, "jet eta range for analysis"};
108105

109106
Configurable<std::string> eventSelections{"eventSelections", "sel8FullPbPb", "choose event selection"};
110107
Configurable<std::string> trackSelections{"trackSelections", "globalTracks", "set track selections"};
@@ -117,6 +114,7 @@ struct JetSpectraEseTask {
117114

118115
Configurable<int> cfgnTotalSystem{"cfgnTotalSystem", 7, "total qvector number // look in Qvector table for this number"};
119116
Configurable<int> cfgnCorrLevel{"cfgnCorrLevel", 3, "QVector step: 0 = no corr, 1 = rect, 2 = twist, 3 = full"};
117+
Configurable<int> cfgnredCorrLevel{"cfgnredCorrLevel", 2, "For the correlation of the reduced q vector, QVector step: 0 = no corr, 1 = rect, 2 = twist, 3 = full"};
120118

121119
Configurable<std::string> cfgEPRefA{"cfgEPRefA", "FT0A", "EP reference A"};
122120
Configurable<std::string> cfgEPRefB{"cfgEPRefB", "TPCpos", "EP reference B"};
@@ -460,13 +458,23 @@ struct JetSpectraEseTask {
460458
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi3FT0C");
461459
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi3FT0A");
462460

463-
registry.add("eventQA/hCosPsi2AmC", ";Centrality;cos(2(#Psi_{2}^{A}-#Psi_{2}^{B}));#it{q}_{2}", {HistType::kTH3F, {{centAxis}, {cosAxis}, {eseAxis}}});
464-
registry.addClone("eventQA/hCosPsi2AmC", "eventQA/hCosPsi2AmB");
465-
registry.addClone("eventQA/hCosPsi2AmC", "eventQA/hCosPsi2BmC");
461+
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi4FT0C");
462+
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi4FT0A");
463+
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi4FV0A");
464+
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi4TPCpos");
465+
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hPsi4TPCneg");
466+
467+
registry.add("eventQA/hCos2Psi2AmC", ";Centrality;cos(2(#Psi_{2}^{A}-#Psi_{2}^{B}));#it{q}_{2}", {HistType::kTH3F, {{centAxis}, {cosAxis}, {eseAxis}}});
468+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos2Psi2AmB");
469+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos2Psi2BmC");
466470

467-
registry.addClone("eventQA/hCosPsi2AmC", "eventQA/hCos4PsiAmC");
468-
registry.addClone("eventQA/hCosPsi2AmC", "eventQA/hCos4PsiAmB");
469-
registry.addClone("eventQA/hCosPsi2AmC", "eventQA/hCos4PsiBmC");
471+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi2AmC");
472+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi2AmB");
473+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi2BmC");
474+
475+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi4AmC");
476+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi4AmB");
477+
registry.addClone("eventQA/hCos2Psi2AmC", "eventQA/hCos4Psi4BmC");
470478

471479
registry.add("eventQA/hQvecUncorV2", ";Centrality;Q_x;Q_y", {HistType::kTH3F, {{centAxis}, {qvecAxis}, {qvecAxis}}});
472480
registry.addClone("eventQA/hQvecUncorV2", "eventQA/hQvecRectrV2");
@@ -475,6 +483,8 @@ struct JetSpectraEseTask {
475483
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hEPUncorV2");
476484
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hEPRectrV2");
477485
registry.addClone("eventQA/hPsi2FT0C", "eventQA/hEPTwistV2");
486+
487+
registry.add("eventQA/h3Centq2FT0Cq2FT0A", ";Centrality;#it{q}_{2}^{FT0C};#it{q}_{2}^{FT0A}", {HistType::kTH3F, {{centAxis}, {250, 0, 35}, {250, 0, 35}}});
478488
}
479489
if (doprocessESEBackground) {
480490
LOGF(info, "JetSpectraEseTask::init() - Background Process");
@@ -674,7 +684,7 @@ struct JetSpectraEseTask {
674684

675685
auto corrL = [&](const auto& j) { return j.pt() - evalRho(rhoFit.get(), jetR, j.phi(), collision.rho()) * j.area(); };
676686
for (const auto& jet : jets) {
677-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax))
687+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax))
678688
continue;
679689
// if (!isAcceptedJet<aod::JetTracks>(jet)) {
680690
if (!isAcceptedJet<soa::Join<aod::JetTracks, aod::JTrackPIs>>(jet)) {
@@ -809,7 +819,7 @@ struct JetSpectraEseTask {
809819

810820
auto corrL = [&](const auto& j) { return j.pt() - evalRho(rhoFit.get(), jetR, j.phi(), c1.rho()) * j.area(); };
811821
for (const auto& jet : jets1) {
812-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax))
822+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax))
813823
continue;
814824
if (!isAcceptedJet<soa::Join<aod::JetTracks, aod::JTrackPIs>>(jet)) {
815825
continue;
@@ -896,6 +906,7 @@ struct JetSpectraEseTask {
896906
return;
897907

898908
[[maybe_unused]] const auto psi{procEP<PsiFillerEP>(collision)};
909+
detCorrelation(collision);
899910
}
900911
PROCESS_SWITCH(JetSpectraEseTask, processESEEPData, "process ese collisions for filling EP and EPR", false);
901912

@@ -939,7 +950,7 @@ struct JetSpectraEseTask {
939950
registry.fill(HIST("hTrackPhi"), collision.centFT0M(), track.phi(), occupancy);
940951
}
941952
for (const auto& jet : jets) {
942-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax))
953+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax))
943954
continue;
944955
if (!isAcceptedJet<aod::JetTracks>(jet)) {
945956
continue;
@@ -1203,6 +1214,27 @@ struct JetSpectraEseTask {
12031214
}
12041215
PROCESS_SWITCH(JetSpectraEseTask, processMCRecoTrack, "jet MC process: Reconstructed track", false);
12051216

1217+
float redqn(const std::vector<float>& vec)
1218+
{
1219+
float redq = -999.0f;
1220+
if (vec[2] > LowFT0Cut) {
1221+
redq = std::sqrt(vec[0] * vec[0] + vec[1] * vec[1]) * std::sqrt(vec[2]);
1222+
}
1223+
return redq;
1224+
}
1225+
template <typename Col>
1226+
void detCorrelation(Col const& col)
1227+
{
1228+
1229+
auto qC = qVecNoESE<DetID::FT0C, false>(col, 2, cfgnredCorrLevel);
1230+
auto qA = qVecNoESE<DetID::FT0A, false>(col, 2, cfgnredCorrLevel);
1231+
1232+
auto redFT0C = redqn(qC);
1233+
auto redFT0A = redqn(qA);
1234+
1235+
registry.fill(HIST("eventQA/h3Centq2FT0Cq2FT0A"), col.centFT0M(), redFT0C, redFT0A);
1236+
}
1237+
12061238
static constexpr float InvalidValue = 999.;
12071239

12081240
template <EventPlaneFiller P, typename EPCol>
@@ -1213,52 +1245,66 @@ struct JetSpectraEseTask {
12131245
};
12141246

12151247
std::map<std::string, float> epMap{
1216-
{"FT0A", computeEP(qVecNoESE<DetID::FT0A, P.hist>(vec, 2), LowFT0Cut, 2.0f)},
1217-
{"FT0C", computeEP(qVecNoESE<DetID::FT0C, false>(vec, 2), LowFT0Cut, 2.0f)},
1218-
{"FV0A", computeEP(qVecNoESE<DetID::FV0A, false>(vec), LowFT0Cut, 2.0f)},
1219-
{"TPCpos", computeEP(qVecNoESE<DetID::TPCpos, false>(vec), 0.0f, 2.0f)},
1220-
{"TPCneg", computeEP(qVecNoESE<DetID::TPCneg, false>(vec), 0.0f, 2.0f)}};
1248+
{"FT0A", computeEP(qVecNoESE<DetID::FT0A, P.hist>(vec, 2, cfgnCorrLevel), LowFT0Cut, 2.0f)},
1249+
{"FT0C", computeEP(qVecNoESE<DetID::FT0C, false>(vec, 2, cfgnCorrLevel), LowFT0Cut, 2.0f)},
1250+
{"FV0A", computeEP(qVecNoESE<DetID::FV0A, false>(vec, 2, cfgnCorrLevel), LowFT0Cut, 2.0f)},
1251+
{"TPCpos", computeEP(qVecNoESE<DetID::TPCpos, false>(vec, 2, cfgnCorrLevel), 0.0f, 2.0f)},
1252+
{"TPCneg", computeEP(qVecNoESE<DetID::TPCneg, false>(vec, 2, cfgnCorrLevel), 0.0f, 2.0f)}};
12211253
std::map<std::string, float> ep3Map{
1222-
{"FT0A", computeEP(qVecNoESE<DetID::FT0A, false>(vec, 3), LowFT0Cut, 3.0f)},
1223-
{"FT0C", computeEP(qVecNoESE<DetID::FT0C, false>(vec, 3), LowFT0Cut, 3.0f)}};
1254+
{"FT0A", computeEP(qVecNoESE<DetID::FT0A, false>(vec, 3, cfgnCorrLevel), LowFT0Cut, 3.0f)},
1255+
{"FT0C", computeEP(qVecNoESE<DetID::FT0C, false>(vec, 3, cfgnCorrLevel), LowFT0Cut, 3.0f)}};
1256+
std::map<std::string, float> ep4Map{
1257+
{"FT0A", computeEP(qVecNoESE<DetID::FT0A, false>(vec, 4, cfgnCorrLevel), LowFT0Cut, 4.0f)},
1258+
{"FT0C", computeEP(qVecNoESE<DetID::FT0C, false>(vec, 4, cfgnCorrLevel), LowFT0Cut, 4.0f)},
1259+
{"FV0A", computeEP(qVecNoESE<DetID::FV0A, false>(vec, 4, cfgnCorrLevel), LowFT0Cut, 4.0f)},
1260+
{"TPCpos", computeEP(qVecNoESE<DetID::TPCpos, false>(vec, 4, cfgnCorrLevel), 0.0f, 4.0f)},
1261+
{"TPCneg", computeEP(qVecNoESE<DetID::TPCneg, false>(vec, 4, cfgnCorrLevel), 0.0f, 4.0f)}};
12241262

12251263
if constexpr (P.psi) {
12261264
if constexpr (P.hist) {
1227-
fillEP(vec, epMap, ep3Map);
1265+
fillEP(vec, epMap, ep3Map, ep4Map);
12281266
}
12291267

12301268
auto cosPsi = [](float psiX, float psiY, const float harmonic = 2.0f) {
12311269
return psiX == InvalidValue || psiY == InvalidValue ? InvalidValue : std::cos(harmonic * (psiX - psiY));
12321270
};
1233-
const std::array<float, 3> epCorrContainer{
1271+
const std::array<float, 3> epCorrContainer22{
12341272
cosPsi(epMap.at(cfgEPRefA), epMap.at(cfgEPRefC)),
12351273
cosPsi(epMap.at(cfgEPRefA), epMap.at(cfgEPRefB)),
12361274
cosPsi(epMap.at(cfgEPRefB), epMap.at(cfgEPRefC))};
1237-
const std::array<float, 3> epCorrContainer4{
1275+
const std::array<float, 3> epCorrContainer24{
12381276
cosPsi(epMap.at(cfgEPRefA), epMap.at(cfgEPRefC), 4.0f),
12391277
cosPsi(epMap.at(cfgEPRefA), epMap.at(cfgEPRefB), 4.0f),
12401278
cosPsi(epMap.at(cfgEPRefB), epMap.at(cfgEPRefC), 4.0f)};
1279+
const std::array<float, 3> epCorrContainer44{
1280+
cosPsi(ep4Map.at(cfgEPRefA), ep4Map.at(cfgEPRefC), 4.0f),
1281+
cosPsi(ep4Map.at(cfgEPRefA), ep4Map.at(cfgEPRefB), 4.0f),
1282+
cosPsi(ep4Map.at(cfgEPRefB), ep4Map.at(cfgEPRefC), 4.0f)};
12411283

12421284
if constexpr (P.hist) {
1243-
fillEPCos(vec, epCorrContainer, epCorrContainer4);
1285+
fillEPCos(vec, epCorrContainer22, epCorrContainer24, epCorrContainer44);
12441286
}
12451287
}
12461288
return {epMap.at(cfgEPRefA), ep3Map.at(cfgEPRefA)};
12471289
}
12481290
template <typename collision>
1249-
void fillEPCos(const collision& col, const std::array<float, 3>& Corr, const std::array<float, 3>& Corr4)
1291+
void fillEPCos(const collision& col, const std::array<float, 3>& Corr22, const std::array<float, 3>& Corr42, const std::array<float, 3>& Corr44)
12501292
{
1251-
registry.fill(HIST("eventQA/hCosPsi2AmC"), col.centFT0M(), Corr[0], col.qPERCFT0C()[0]);
1252-
registry.fill(HIST("eventQA/hCosPsi2AmB"), col.centFT0M(), Corr[1], col.qPERCFT0C()[0]);
1253-
registry.fill(HIST("eventQA/hCosPsi2BmC"), col.centFT0M(), Corr[2], col.qPERCFT0C()[0]);
1293+
registry.fill(HIST("eventQA/hCos2Psi2AmC"), col.centFT0M(), Corr22[0], col.qPERCFT0C()[0]);
1294+
registry.fill(HIST("eventQA/hCos2Psi2AmB"), col.centFT0M(), Corr22[1], col.qPERCFT0C()[0]);
1295+
registry.fill(HIST("eventQA/hCos2Psi2BmC"), col.centFT0M(), Corr22[2], col.qPERCFT0C()[0]);
1296+
1297+
registry.fill(HIST("eventQA/hCos4Psi2AmC"), col.centFT0M(), Corr42[0], col.qPERCFT0C()[0]);
1298+
registry.fill(HIST("eventQA/hCos4Psi2AmB"), col.centFT0M(), Corr42[1], col.qPERCFT0C()[0]);
1299+
registry.fill(HIST("eventQA/hCos4Psi2BmC"), col.centFT0M(), Corr42[2], col.qPERCFT0C()[0]);
12541300

1255-
registry.fill(HIST("eventQA/hCos4PsiAmC"), col.centFT0M(), Corr4[0], col.qPERCFT0C()[0]);
1256-
registry.fill(HIST("eventQA/hCos4PsiAmB"), col.centFT0M(), Corr4[1], col.qPERCFT0C()[0]);
1257-
registry.fill(HIST("eventQA/hCos4PsiBmC"), col.centFT0M(), Corr4[2], col.qPERCFT0C()[0]);
1301+
registry.fill(HIST("eventQA/hCos4Psi4AmC"), col.centFT0M(), Corr44[0], col.qPERCFT0C()[0]);
1302+
registry.fill(HIST("eventQA/hCos4Psi4AmB"), col.centFT0M(), Corr44[1], col.qPERCFT0C()[0]);
1303+
registry.fill(HIST("eventQA/hCos4Psi4BmC"), col.centFT0M(), Corr44[2], col.qPERCFT0C()[0]);
12581304
}
12591305

12601306
template <typename collision>
1261-
void fillEP(const collision& col, const std::map<std::string, float>& epMap, const std::map<std::string, float>& ep3Map)
1307+
void fillEP(const collision& col, const std::map<std::string, float>& epMap, const std::map<std::string, float>& ep3Map, const std::map<std::string, float>& ep4Map)
12621308
{
12631309
registry.fill(HIST("eventQA/hPsi2FT0A"), col.centFT0M(), epMap.at("FT0A"));
12641310
registry.fill(HIST("eventQA/hPsi2FV0A"), col.centFT0M(), epMap.at("FV0A"));
@@ -1268,6 +1314,12 @@ struct JetSpectraEseTask {
12681314

12691315
registry.fill(HIST("eventQA/hPsi3FT0A"), col.centFT0M(), ep3Map.at("FT0A"));
12701316
registry.fill(HIST("eventQA/hPsi3FT0C"), col.centFT0M(), ep3Map.at("FT0C"));
1317+
1318+
registry.fill(HIST("eventQA/hPsi4FT0A"), col.centFT0M(), ep4Map.at("FT0A"));
1319+
registry.fill(HIST("eventQA/hPsi4FT0C"), col.centFT0M(), ep4Map.at("FT0C"));
1320+
registry.fill(HIST("eventQA/hPsi4FV0A"), col.centFT0M(), ep4Map.at("FV0A"));
1321+
registry.fill(HIST("eventQA/hPsi4TPCpos"), col.centFT0M(), ep4Map.at("TPCpos"));
1322+
registry.fill(HIST("eventQA/hPsi4TPCneg"), col.centFT0M(), ep4Map.at("TPCneg"));
12711323
}
12721324
constexpr int detIDN(const DetID id)
12731325
{
@@ -1292,7 +1344,7 @@ struct JetSpectraEseTask {
12921344

12931345
const int secondHarmonic{2};
12941346
template <DetID id, bool fill, typename Col>
1295-
std::vector<float> qVecNoESE(Col collision, int nmode = 2)
1347+
std::vector<float> qVecNoESE(Col collision, int nmode = 2, int corrLevel = 3)
12961348
{
12971349
int detId{detIDN(id)};
12981350
int detInd{detId * 4 + cfgnTotalSystem * 4 * (nmode - 2)};
@@ -1308,8 +1360,8 @@ struct JetSpectraEseTask {
13081360
}
13091361
}
13101362
std::vector<float> qVec{};
1311-
qVec.push_back(collision.qvecRe()[detInd + cfgnCorrLevel]);
1312-
qVec.push_back(collision.qvecIm()[detInd + cfgnCorrLevel]);
1363+
qVec.push_back(collision.qvecRe()[detInd + corrLevel]);
1364+
qVec.push_back(collision.qvecIm()[detInd + corrLevel]);
13131365
qVec.push_back(collision.qvecAmp()[detId]);
13141366
return qVec;
13151367
}
@@ -1370,7 +1422,7 @@ struct JetSpectraEseTask {
13701422
for (const auto& track : tracks) {
13711423
if constexpr (fillHist)
13721424
registry.fill(HIST("trackQA/hRhoTrackCounter"), 0.5);
1373-
if (jetderiveddatautilities::selectTrack(track, trackSelection) && (std::fabs(track.eta() - leadingJetEta) > jetR) && track.pt() >= trackPtMinRhoPhi && track.pt() <= trackPtMaxRhoPhi) {
1425+
if (jetderiveddatautilities::selectTrack(track, trackSelection) && (std::fabs(track.eta() - leadingJetEta) > jetR) && track.pt() >= trackPtRhoPhi->at(0) && track.pt() <= trackPtRhoPhi->at(1)) {
13741426
nTrk++;
13751427
if constexpr (fillHist)
13761428
registry.fill(HIST("trackQA/hRhoTrackCounter"), 1.5);
@@ -1382,7 +1434,7 @@ struct JetSpectraEseTask {
13821434

13831435
auto hPhiPt = std::unique_ptr<TH1F>(new TH1F("h_ptsum_sumpt_fit", "h_ptsum_sumpt fit use", TMath::CeilNint(std::sqrt(nTrk)), 0., o2::constants::math::TwoPI));
13841436
for (const auto& track : tracks) {
1385-
if (jetderiveddatautilities::selectTrack(track, trackSelection) && (std::fabs(track.eta() - leadingJetEta) > jetR) && track.pt() >= trackPtMinRhoPhi && track.pt() <= trackPtMaxRhoPhi) {
1437+
if (jetderiveddatautilities::selectTrack(track, trackSelection) && (std::fabs(track.eta() - leadingJetEta) > jetR) && track.pt() >= trackPtRhoPhi->at(0) && track.pt() <= trackPtRhoPhi->at(1)) {
13861438
hPhiPt->Fill(track.phi(), track.pt());
13871439
if constexpr (fillHist) {
13881440
registry.fill(HIST("trackQA/hRhoTrackCounter"), 2.5);
@@ -1617,7 +1669,7 @@ struct JetSpectraEseTask {
16171669
{
16181670
float weight = 1.0;
16191671
for (const auto& jet : jets) {
1620-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1672+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax)) {
16211673
continue;
16221674
}
16231675
if (!isAcceptedJet<JTracks>(jet)) {
@@ -1640,7 +1692,7 @@ struct JetSpectraEseTask {
16401692
bool mcLevelIsParticleLevel = true;
16411693
float weight = 1.0;
16421694
for (const auto& jet : jets) {
1643-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1695+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax)) {
16441696
continue;
16451697
}
16461698
if (!isAcceptedJet<JTracks>(jet, mcLevelIsParticleLevel)) {
@@ -1662,7 +1714,7 @@ struct JetSpectraEseTask {
16621714
{
16631715
float weight = 1.0;
16641716
for (const auto& jet : jets) {
1665-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax)) {
1717+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax)) {
16661718
continue;
16671719
}
16681720
if (!isAcceptedJet<JTracks>(jet)) {
@@ -1711,7 +1763,7 @@ struct JetSpectraEseTask {
17111763
float leadJetEta = -999;
17121764
bool hasLeadingJet = false;
17131765
for (const auto& jet : jets) {
1714-
if (!jetfindingutilities::isInEtaAcceptance(jet, jetEtaMin, jetEtaMax, trackEtaMin, trackEtaMax))
1766+
if (!jetfindingutilities::isInEtaAcceptance(jet, cfgJetEta->at(0), cfgJetEta->at(1), trackEtaMin, trackEtaMax))
17151767
continue;
17161768
if (!isAcceptedJet<soa::Join<aod::JetTracks, aod::JTrackPIs>>(jet)) {
17171769
continue;

0 commit comments

Comments
 (0)