Skip to content

Commit fe46cc4

Browse files
authored
[PWGJE] fixed rho local bug (#10263)
1 parent c83005d commit fe46cc4

File tree

1 file changed

+10
-13
lines changed

1 file changed

+10
-13
lines changed

PWGJE/Tasks/jetSpectraEseTask.cxx

Lines changed: 10 additions & 13 deletions
Original file line numberDiff line numberDiff line change
@@ -287,16 +287,18 @@ struct JetSpectraEseTask {
287287

288288
auto rho = collision.rho();
289289
std::unique_ptr<TF1> rhoFit{nullptr};
290-
if (cfgrhoPhi)
290+
if (cfgrhoPhi) {
291291
rhoFit = fitRho<true>(collision, psi, tracks);
292+
}
292293

293294
registry.fill(HIST("hEventCounter"), counter++);
294295
registry.fill(HIST("hRho"), rho);
295296
registry.fill(HIST("hCentralityAnalyzed"), collision.centrality());
296297
for (auto const& jet : jets) {
297-
float jetpTbkgsub = jet.pt() - (rho * jet.area());
298-
if (cfgrhoPhi)
298+
if (cfgrhoPhi) {
299299
rho = evalRho(rhoFit.get(), jet.phi(), rho);
300+
}
301+
float jetpTbkgsub = jet.pt() - (rho * jet.area());
300302
registry.fill(HIST("hJetPt"), jet.pt());
301303
registry.fill(HIST("hJetPt_bkgsub"), jetpTbkgsub);
302304
registry.fill(HIST("hJetEta"), jet.eta());
@@ -631,17 +633,12 @@ struct JetSpectraEseTask {
631633

632634
hPhiPt->Fit(modulationFit.get(), "Q", "", 0, o2::constants::math::TwoPI);
633635

634-
std::array<double, 5> fitParams{};
635-
for (int i{0}; i < 5; ++i) {
636-
fitParams[i] = modulationFit->GetParameter(i);
637-
}
638-
639636
if constexpr (fillHist) {
640-
registry.fill(HIST("hfitPar0"), col.centrality(), fitParams[0]);
641-
registry.fill(HIST("hfitPar1"), col.centrality(), fitParams[1]);
642-
registry.fill(HIST("hfitPar2"), col.centrality(), fitParams[2]);
643-
registry.fill(HIST("hfitPar3"), col.centrality(), fitParams[3]);
644-
registry.fill(HIST("hfitPar4"), col.centrality(), fitParams[4]);
637+
registry.fill(HIST("hfitPar0"), col.centrality(), modulationFit->GetParameter(0));
638+
registry.fill(HIST("hfitPar1"), col.centrality(), modulationFit->GetParameter(1));
639+
registry.fill(HIST("hfitPar2"), col.centrality(), modulationFit->GetParameter(2));
640+
registry.fill(HIST("hfitPar3"), col.centrality(), modulationFit->GetParameter(3));
641+
registry.fill(HIST("hfitPar4"), col.centrality(), modulationFit->GetParameter(4));
645642
}
646643

647644
double chi2{0.};

0 commit comments

Comments
 (0)