Skip to content

[PWGLF] Apply missing cuts before filling dE/dx vs momentum histograms #16848

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

Open
PaolaVT wants to merge 23 commits into
base: master
Choose a base branch
from
+32 −32
Open

[PWGLF] Apply missing cuts before filling dE/dx vs momentum histograms #16848

Changes from all commits
Commits
Show all changes
23 commits
Select commit Hold shift + click to select a range
1399f53
New hist was added
Apr 24, 2026
83e925b
Merge branch 'AliceO2Group:master' into master
PaolaVT May 5, 2026
7e25b0d
DCA cuts were changed, mass plots were added
May 5, 2026
d33387b
Merge branch 'AliceO2Group:master' into master
PaolaVT May 6, 2026
51f7078
Merge branch 'AliceO2Group:master' into master
PaolaVT May 8, 2026
460b55d
Histograms for the p to pT conversion were added
May 8, 2026
45d71ca
Histograms for the p to pT conversion
May 8, 2026
3de3411
Merge branch 'AliceO2Group:master' into master
PaolaVT May 13, 2026
6f85ddd
Merge branch 'AliceO2Group:master' into master
PaolaVT May 13, 2026
e2acbfb
Histograms for the p to pT conversion
May 13, 2026
53eaa79
Merge branch 'AliceO2Group:master' into master
PaolaVT Jun 2, 2026
ebacc3e
Histograms with a fine binning was added
Jun 2, 2026
9d64323
Added histograms with finer binning
Jun 2, 2026
5a06c3b
Added histograms with finer binning
Jun 2, 2026
fc3b28b
Merge branch 'AliceO2Group:master' into master
PaolaVT Jun 9, 2026
a08db2b
Added histograms with finer binning for MC and make momentum configur…
Jun 9, 2026
e9a6946
Merge branch 'AliceO2Group:master' into master
PaolaVT Jun 12, 2026
b4bccef
Debug analysis
Jun 12, 2026
1d3bced
Applied clang format
Jun 12, 2026
a8d02c8
Merge branch 'AliceO2Group:master' into master
PaolaVT Jun 12, 2026
e4e3f61
correct array bounds
Jun 12, 2026
0039d75
Merge branch 'AliceO2Group:master' into master
PaolaVT Jun 28, 2026
4783643
Apply missing cuts before filling dE/dx vs momentum histograms
Jun 28, 2026
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
64 changes: 32 additions & 32 deletions PWGLF/Tasks/Nuspex/multiplicityPt.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1172,38 +1172,6 @@ struct MultiplicityPt {
if (centIndex == -1)
continue;

// ====================================================================
// DEDX VS MOMENTUM HISTOGRAMS FILLING - ALL TRACKS
// ====================================================================
hDedxVspTMomentumVsCent[10]->Fill(track.pt(), tpcSignal, eta);
if (charge > 0) {
registry.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Pos"), momentum, tpcSignal, eta);
hDedxVsMomentumVsCentPos[centIndex]->Fill(momentum, tpcSignal, eta);
hDedxVspTMomentumVsCent[centIndex]->Fill(track.pt(), tpcSignal, eta);
hMomentumVsEtaPos[centIndex]->Fill(eta, momentum);
hMomentumVsEtaPos[10]->Fill(eta, momentum);
hpTVsEtaPos[centIndex]->Fill(eta, track.pt());
hpTVsEtaPos[10]->Fill(eta, track.pt());
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos"), eta, track.pt(), momentum);
} else {
registry.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Neg"), momentum, tpcSignal, eta);
hDedxVsMomentumVsCentNeg[centIndex]->Fill(momentum, tpcSignal, eta);
hDedxVspTMomentumVsCent[centIndex]->Fill(track.pt(), tpcSignal, eta);
hMomentumVsEtaNeg[centIndex]->Fill(eta, momentum);
hMomentumVsEtaNeg[10]->Fill(eta, momentum);
hpTVsEtaNeg[centIndex]->Fill(eta, track.pt());
hpTVsEtaNeg[10]->Fill(eta, track.pt());
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg"), eta, track.pt(), momentum);
}

if (isPrimary) {
if (charge > 0) {
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos_Pri"), eta, track.pt(), momentum);
} else {
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg_Pri"), eta, track.pt(), momentum);
}
}

registry.fill(HIST("hEta"), track.eta());
registry.fill(HIST("hPhi"), track.phi());
registry.fill(HIST("EtaVsPhi"), track.eta(), track.phi());
Expand Down Expand Up @@ -1318,6 +1286,38 @@ struct MultiplicityPt {
}
}

// ====================================================================
// DEDX VS MOMENTUM HISTOGRAMS FILLING - ALL TRACKS
// ====================================================================
hDedxVspTMomentumVsCent[10]->Fill(track.pt(), tpcSignal, eta);
if (charge > 0) {
registry.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Pos"), momentum, tpcSignal, eta);
hDedxVsMomentumVsCentPos[centIndex]->Fill(momentum, tpcSignal, eta);
hDedxVspTMomentumVsCent[centIndex]->Fill(track.pt(), tpcSignal, eta);
hMomentumVsEtaPos[centIndex]->Fill(eta, momentum);
hMomentumVsEtaPos[10]->Fill(eta, momentum);
hpTVsEtaPos[centIndex]->Fill(eta, track.pt());
hpTVsEtaPos[10]->Fill(eta, track.pt());
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos"), eta, track.pt(), momentum);
} else {
registry.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Neg"), momentum, tpcSignal, eta);
hDedxVsMomentumVsCentNeg[centIndex]->Fill(momentum, tpcSignal, eta);
hDedxVspTMomentumVsCent[centIndex]->Fill(track.pt(), tpcSignal, eta);
hMomentumVsEtaNeg[centIndex]->Fill(eta, momentum);
hMomentumVsEtaNeg[10]->Fill(eta, momentum);
hpTVsEtaNeg[centIndex]->Fill(eta, track.pt());
hpTVsEtaNeg[10]->Fill(eta, track.pt());
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg"), eta, track.pt(), momentum);
}

if (isPrimary) {
if (charge > 0) {
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos_Pri"), eta, track.pt(), momentum);
} else {
registry.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg_Pri"), eta, track.pt(), momentum);
}
}

// ====================================================================
// DEDX VS MOMENTUM HISTOGRAMS FILLING - PARTICLE SPECIFIC
// ====================================================================
Expand Down
Loading