Skip to content

Commit 555261c

Browse files
PaolaVTPaola Vargas Torres
PaolaVT
and
Paola Vargas Torres
authored
[PWGLF] Added histograms with finer binning (#16513)
Co-authored-by: Paola Vargas Torres <paovargas@MacBook-Air-de-Paola.local>
1 parent 2dd7aca commit 555261c

2 files changed

Lines changed: 63 additions & 10 deletions

File tree

‎PWGLF/Tasks/Nuspex/dedxPidAnalysis.cxx‎

Lines changed: 25 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -69,6 +69,10 @@ static constexpr int NCentHists{10};
6969
std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentPos{};
7070
std::array<std::shared_ptr<TH3>, NCentHists> hDedxVsMomentumVsCentNeg{};
7171
std::array<std::shared_ptr<TH3>, NCentHists + 1> hDedxVspTMomentumVsCent{};
72+
std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaPos{};
73+
std::array<std::shared_ptr<TH2>, NCentHists + 1> hMomentumVsEtaNeg{};
74+
std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaPos{};
75+
std::array<std::shared_ptr<TH2>, NCentHists + 1> hpTVsEtaNeg{};
7276

7377
struct DedxPidAnalysis {
7478

@@ -247,6 +251,11 @@ struct DedxPidAnalysis {
247251
static constexpr std::string_view DedxvsMomentumvsCentPos[CentralityClasses] = {"dEdx_vs_Momentum_Cent0_1_Pos", "dEdx_vs_Momentum_Cent1_5_Pos", "dEdx_vs_Momentum_Cent5_10_Pos", "dEdx_vs_Momentum_Cent10_15_Pos", "dEdx_vs_Momentum_Cent15_20_Pos", "dEdx_vs_Momentum_Cent20_30_Pos", "dEdx_vs_Momentum_Cent30_40_Pos", "dEdx_vs_Momentum_Cent40_50_Pos", "dEdx_vs_Momentum_Cent50_70_Pos", "dEdx_vs_Momentum_Cent70_100_Pos"};
248252
static constexpr std::string_view DedxvsMomentumvsCentNeg[CentralityClasses] = {"dEdx_vs_Momentum_Cent0_1_Neg", "dEdx_vs_Momentum_Cent1_5_Neg", "dEdx_vs_Momentum_Cent5_10_Neg", "dEdx_vs_Momentum_Cent10_15_Neg", "dEdx_vs_Momentum_Cent15_20_Neg", "dEdx_vs_Momentum_Cent20_30_Neg", "dEdx_vs_Momentum_Cent30_40_Neg", "dEdx_vs_Momentum_Cent40_50_Neg", "dEdx_vs_Momentum_Cent50_70_Neg", "dEdx_vs_Momentum_Cent70_100_Neg"};
249253
static constexpr std::string_view DedxvspTMomentumvsCent[CentralityClasses + 1] = {"dEdx_vs_pTMomentum_Cent0_1", "dEdx_vs_pTMomentum_Cent1_5", "dEdx_vs_pTMomentum_Cent5_10", "dEdx_vs_pTMomentum_Cent10_15", "dEdx_vs_pTMomentum_Cent15_20", "dEdx_vs_pTMomentum_Cent20_30", "dEdx_vs_pTMomentum_Cent30_40", "dEdx_vs_pTMomentum_Cent40_50", "dEdx_vs_pTMomentum_Cent50_70", "dEdx_vs_pTMomentum_Cent70_100", "dEdx_vs_pTMomentum_all_Pos"};
254+
// Fine binnind
255+
static constexpr std::string_view CentpPos[CentralityClasses + 1] = {"p_vs_eta_Cent0_1_Pos", "p_vs_eta_Cent1_5_Pos", "p_vs_eta_Cent5_10_Pos", "p_vs_eta_Cent10_15_Pos", "p_vs_eta_Cent15_20_Pos", "p_vs_eta_Cent20_30_Pos", "p_vs_eta_Cent30_40_Pos", "p_vs_eta_Cent40_50_Pos", "p_vs_eta_Cent50_70_Pos", "p_vs_eta_Cent70_100_Pos", "p_vs_eta_MB_Pos"};
256+
static constexpr std::string_view CentpNeg[CentralityClasses + 1] = {"p_vs_eta_Cent0_1_Neg", "p_vs_eta_Cent1_5_Neg", "p_vs_eta_Cent5_10_Neg", "p_vs_eta_Cent10_15_Neg", "p_vs_eta_Cent15_20_Neg", "p_vs_eta_Cent20_30_Neg", "p_vs_eta_Cent30_40_Neg", "p_vs_eta_Cent40_50_Neg", "p_vs_eta_Cent50_70_Neg", "p_vs_eta_Cent70_100_Neg", "p_vs_eta_MB_Neg"};
257+
static constexpr std::string_view CentpTPos[CentralityClasses + 1] = {"pT_vs_eta_Cent0_1_Pos", "pT_vs_eta_Cent1_5_Pos", "pT_vs_eta_Cent5_10_Pos", "pT_vs_eta_Cent10_15_Pos", "pT_vs_eta_Cent15_20_Pos", "pT_vs_eta_Cent20_30_Pos", "pT_vs_eta_Cent30_40_Pos", "pT_vs_eta_Cent40_50_Pos", "pT_vs_eta_Cent50_70_Pos", "pT_vs_eta_Cent70_100_Pos", "pT_vs_eta_MB_Pos"};
258+
static constexpr std::string_view CentpTNeg[CentralityClasses + 1] = {"pT_vs_eta_Cent0_1_Neg", "pT_vs_eta_Cent1_5_Neg", "pT_vs_eta_Cent5_10_Neg", "pT_vs_eta_Cent10_15_Neg", "pT_vs_eta_Cent15_20_Neg", "pT_vs_eta_Cent20_30_Neg", "pT_vs_eta_Cent30_40_Neg", "pT_vs_eta_Cent40_50_Neg", "pT_vs_eta_Cent50_70_Neg", "pT_vs_eta_Cent70_100_Neg", "pT_vs_eta_MB_Neg"};
250259
// Ncl TPC
251260
static constexpr std::string_view NclTPCDedxMomentumNegBefore[EtaIntervals] = {"Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_1_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_2_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_3_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_4_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_5_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_6_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_7_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Neg_8_Before"};
252261
static constexpr std::string_view NclTPCDedxMomentumPosBefore[EtaIntervals] = {"Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_1_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_2_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_3_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_4_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_5_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_6_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_7_Before", "Ncl_FoundTPC_vs_dEdx_vs_Momentum_Pos_8_Before"};
@@ -263,6 +272,7 @@ struct DedxPidAnalysis {
263272
Configurable<std::vector<float>> calibrationFactorPos{"calibrationFactorPos", {50.5157, 50.6359, 50.3198, 49.3345, 48.9197, 49.4931, 50.0188, 50.1406}, "positive calibration factors"};
264273
ConfigurableAxis binP{"binP", {VARIABLE_WIDTH, 0.1, 0.12, 0.14, 0.16, 0.18, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.5, 2.6, 2.7, 2.8, 2.9, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 18.0, 20.0}, ""};
265274
ConfigurableAxis centBins{"centBins", {100, 0, 100}, "Binning for centralidad"};
275+
ConfigurableAxis pFineBins{"pFineBins", {1995, 0.1, 40}, "Binning for momentum"};
266276

267277
// phi cut fits
268278
TF1* fphiCutHigh = nullptr;
@@ -351,11 +361,13 @@ struct DedxPidAnalysis {
351361
LOGF(info, "Centrality clases between %.1f - %.1f", CentClasses[0], CentClasses[10]);
352362

353363
AxisSpec dedxAxis{100, 0.0, 100.0, "dE/dx (a. u.)"};
354-
AxisSpec ptAxis = {binP, "pT (GeV/c)"};
364+
AxisSpec ptAxis = {binP, "#it{p}_{T} (GeV/c)"};
355365
AxisSpec etaAxis{8, -0.8, 0.8, "#eta"};
356366
AxisSpec pAxis = {binP, "#it{p}/Z (GeV/c)"};
357367
AxisSpec pAxisTrack = {binP, "#it{p} (GeV/c)"};
358368
AxisSpec centAxis{centBins, "Undefined multiplicity estimator"};
369+
AxisSpec pFineAxis{pFineBins, "#it{p} (GeV/c)"};
370+
AxisSpec pTFineAxis{pFineBins, "#it{p}_{T} (GeV/c)"};
359371
switch (multiplicityEstimator) {
360372
case MultSelectionMode::NoMultiplicity: // No multiplicity
361373
LOGF(info, "No multiplicity estimator applied");
@@ -488,6 +500,10 @@ struct DedxPidAnalysis {
488500

489501
for (int i = 0; i < CentralityClasses + 1; ++i) {
490502
hDedxVspTMomentumVsCent[i] = registryDeDx.add<TH3>(DedxvspTMomentumvsCent[i].data(), "dE/dx", HistType::kTH3F, {{ptAxis}, {dedxAxis}, {etaAxis}});
503+
hMomentumVsEtaPos[i] = registryDeDx.add<TH2>(CentpPos[i].data(), "p vs eta", HistType::kTH2F, {{etaAxis}, {pFineAxis}});
504+
hMomentumVsEtaNeg[i] = registryDeDx.add<TH2>(CentpNeg[i].data(), "p vs eta", HistType::kTH2F, {{etaAxis}, {pFineAxis}});
505+
hpTVsEtaPos[i] = registryDeDx.add<TH2>(CentpTPos[i].data(), "pT vs eta", HistType::kTH2F, {{etaAxis}, {pTFineAxis}});
506+
hpTVsEtaNeg[i] = registryDeDx.add<TH2>(CentpTNeg[i].data(), "pT vs eta", HistType::kTH2F, {{etaAxis}, {pTFineAxis}});
491507
}
492508

493509
// Invariant Mass
@@ -1704,12 +1720,20 @@ struct DedxPidAnalysis {
17041720
hDedxVsMomentumVsCentPos[centIndex]->Fill(signedP, trk.tpcSignal() * 50 / calibrationFactorPos->at(i), trk.eta());
17051721
hDedxVspTMomentumVsCent[centIndex]->Fill(signedpT, trk.tpcSignal() * 50 / calibrationFactorPos->at(i), trk.eta());
17061722
hDedxVspTMomentumVsCent[10]->Fill(signedpT, trk.tpcSignal() * 50 / calibrationFactorPos->at(i), trk.eta());
1723+
hMomentumVsEtaPos[centIndex]->Fill(trk.eta(), signedP);
1724+
hMomentumVsEtaPos[10]->Fill(trk.eta(), signedP);
1725+
hpTVsEtaPos[centIndex]->Fill(trk.eta(), signedpT);
1726+
hpTVsEtaPos[10]->Fill(trk.eta(), signedpT);
17071727
} else {
17081728
registryDeDx.fill(HIST(DedxvsMomentumNeg[0]), std::abs(signedP), trk.tpcSignal() * 50 / calibrationFactorNeg->at(i), trk.eta());
17091729
registryDeDx.fill(HIST("heta_vs_pt_vs_p_all_Neg"), trk.eta(), trk.pt(), trk.p());
17101730
hDedxVsMomentumVsCentNeg[centIndex]->Fill(std::abs(signedP), trk.tpcSignal() * 50 / calibrationFactorNeg->at(i), trk.eta());
17111731
hDedxVspTMomentumVsCent[centIndex]->Fill(std::abs(signedpT), trk.tpcSignal() * 50 / calibrationFactorNeg->at(i), trk.eta());
17121732
hDedxVspTMomentumVsCent[10]->Fill(std::abs(signedpT), trk.tpcSignal() * 50 / calibrationFactorNeg->at(i), trk.eta());
1733+
hMomentumVsEtaNeg[centIndex]->Fill(trk.eta(), std::abs(signedP));
1734+
hMomentumVsEtaNeg[10]->Fill(trk.eta(), std::abs(signedP));
1735+
hpTVsEtaNeg[centIndex]->Fill(trk.eta(), std::abs(signedpT));
1736+
hpTVsEtaNeg[10]->Fill(trk.eta(), std::abs(signedpT));
17131737
}
17141738
}
17151739
}

‎PWGLF/Tasks/Nuspex/multiplicityPt.cxx‎

Lines changed: 38 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -69,6 +69,7 @@ struct MultiplicityPt {
6969
static constexpr int MultBinMax = 200;
7070
static constexpr int RecMultBinMax = 100;
7171
static constexpr int ParticlesType = 4;
72+
static constexpr int ResponseMatrixTypes = 7;
7273

7374
enum INELCutSelection : int {
7475
INEL = 0,
@@ -116,6 +117,9 @@ struct MultiplicityPt {
116117
// Histogram names for V0s dE/dx analysis
117118
static constexpr std::string_view DedxvsMomentumPos[ParticlesType] = {"dEdx_vs_Momentum_all_Pos", "dEdx_vs_Momentum_Pi_v0_Pos", "dEdx_vs_Momentum_Pr_v0_Pos", "dEdx_vs_Momentum_El_v0_Pos"};
118119
static constexpr std::string_view DedxvsMomentumNeg[ParticlesType] = {"dEdx_vs_Momentum_all_Neg", "dEdx_vs_Momentum_Pi_v0_Neg", "dEdx_vs_Momentum_Pr_v0_Neg", "dEdx_vs_Momentum_El_v0_Neg"};
120+
121+
static constexpr std::string_view EtavspvspTPosPart[ResponseMatrixTypes] = {"heta_vs_pt_vs_p_all_Pos", "heta_vs_pt_vs_p_all_Pos_Pri", "heta_vs_pt_vs_p_all_Pos_Pri_MC", "heta_vs_pt_vs_p_all_Pos_Pri_MC_Part", "heta_vs_pt_vs_p_Pi_Pos", "heta_vs_pt_vs_p_K_Pos", "heta_vs_pt_vs_p_Pr_Pos"};
122+
static constexpr std::string_view EtavspvspTNegPart[ResponseMatrixTypes] = {"heta_vs_pt_vs_p_all_Neg", "heta_vs_pt_vs_p_all_Neg_Pri", "heta_vs_pt_vs_p_all_Neg_Pri_MC", "heta_vs_pt_vs_p_all_Neg_Pri_MC_Part", "heta_vs_pt_vs_p_Pi_Neg", "heta_vs_pt_vs_p_K_Neg", "heta_vs_pt_vs_p_Pr_Neg"};
119123
// Particle fractions histograms
120124
static constexpr std::string_view ParticleFractionsVsMomentumPos[ParticlesType + 1] = {"hFractionVsMomentum_Pion_Pos", "hFractionVsMomentum_Kaon_Pos", "hFractionVsMomentum_Proton_Pos", "hFractionVsMomentum_Electron_Pos", "hFractionVsMomentum_Muon_Pos"};
121125

@@ -482,6 +486,15 @@ void MultiplicityPt::init(InitContext const&)
482486
"dE/dx vs Momentum Negative", HistType::kTH3F,
483487
{{pAxis}, {dedxAxis}, {etaAxis}});
484488
}
489+
// pt vs p
490+
for (int i = 0; i < ResponseMatrixTypes; ++i) {
491+
ue.add(("ResponseMatrix/" + std::string(EtavspvspTPosPart[i])).c_str(),
492+
"eta vs pT vs p Positive", HistType::kTH3F,
493+
{{etaAxis}, {ptAxis}, {pAxis}});
494+
ue.add(("ResponseMatrix/" + std::string(EtavspvspTNegPart[i])).c_str(),
495+
"eta vs pT vs p Negative", HistType::kTH3F,
496+
{{etaAxis}, {ptAxis}, {pAxis}});
497+
}
485498

486499
// ===== Particle Fractions as function of p and pT =====
487500
ue.add("ParticleFractions/hTotalCountsVsMomentumPos", "Total counts vs momentum;#it{p} (GeV/#it{c});Counts", HistType::kTH2D, {{etaAxis}, {pAxis}});
@@ -499,13 +512,6 @@ void MultiplicityPt::init(InitContext const&)
499512
ue.add(("ParticleFractions/" + std::string(ParticleFractionsVsPtNeg[i])).c_str(),
500513
"Particle fraction vs pT", HistType::kTH2D, {{etaAxis}, {ptAxis}});
501514
}
502-
// pt vs p
503-
ue.add(
504-
"heta_vs_pt_vs_p_all_Neg", "eta_vs_pT_vs_p", HistType::kTH3F,
505-
{{etaAxis}, {ptAxis}, {pAxis}});
506-
ue.add(
507-
"heta_vs_pt_vs_p_all_Pos", "eta_vs_pT_vs_p", HistType::kTH3F,
508-
{{etaAxis}, {ptAxis}, {pAxis}});
509515

510516
LOG(info) << "=== Initialization complete ===";
511517
}
@@ -754,24 +760,41 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
754760
// dedx for all particles
755761
if (charge > 0) {
756762
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Pos"), momentum, tpcSignal, eta);
757-
ue.fill(HIST("heta_vs_pt_vs_p_all_Pos"), eta, track.pt(), momentum);
763+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos"), eta, track.pt(), momentum);
758764
} else {
759765
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_all_Neg"), momentum, tpcSignal, eta);
760-
ue.fill(HIST("heta_vs_pt_vs_p_all_Neg"), eta, track.pt(), momentum);
766+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg"), eta, track.pt(), momentum);
767+
}
768+
769+
if (track.mcParticle().isPhysicalPrimary()) {
770+
if (charge > 0) {
771+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos_Pri"), eta, track.pt(), momentum);
772+
} else {
773+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg_Pri"), eta, track.pt(), momentum);
774+
}
761775
}
762776

763777
if (track.has_mcParticle()) {
764778
const auto& particle = track.mcParticle();
765779
int pdgCode = std::abs(particle.pdgCode());
780+
if (particle.isPhysicalPrimary()) {
781+
if (charge > 0) {
782+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos_Pri_MC"), eta, track.pt(), momentum);
783+
} else {
784+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg_Pri_MC"), eta, track.pt(), momentum);
785+
}
786+
}
766787
if (pdgCode == PDG_t::kPiPlus || pdgCode == PDG_t::kKPlus || pdgCode == PDG_t::kProton || pdgCode == PDG_t::kElectron || pdgCode == PDG_t::kMuonPlus) {
767788
if (particle.isPhysicalPrimary()) {
768789
// Fill total counts for fractions
769790
if (charge > 0) {
770791
ue.fill(HIST("ParticleFractions/hTotalCountsVsMomentumPos"), eta, momentum);
771792
ue.fill(HIST("ParticleFractions/hTotalCountsVsPtPos"), eta, track.pt());
793+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Pos_Pri_MC_Part"), eta, track.pt(), momentum);
772794
} else {
773795
ue.fill(HIST("ParticleFractions/hTotalCountsVsMomentumNeg"), eta, momentum);
774796
ue.fill(HIST("ParticleFractions/hTotalCountsVsPtNeg"), eta, track.pt());
797+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_all_Neg_Pri_MC_Part"), eta, track.pt(), momentum);
775798
}
776799
}
777800
}
@@ -785,10 +808,12 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
785808
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Pion_Pos"), eta, momentum);
786809
ue.fill(HIST("ParticleFractions/hFractionVsPt_Pion_Pos"), eta, track.pt());
787810
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_Pi_v0_Pos"), momentum, tpcSignal, eta);
811+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_Pi_Pos"), eta, track.pt(), momentum);
788812
} else {
789813
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Pion_Neg"), eta, momentum);
790814
ue.fill(HIST("ParticleFractions/hFractionVsPt_Pion_Neg"), eta, track.pt());
791815
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_Pi_v0_Neg"), momentum, tpcSignal, eta);
816+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_Pi_Neg"), eta, track.pt(), momentum);
792817
}
793818

794819
} else {
@@ -802,9 +827,11 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
802827
if (charge > 0) {
803828
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Kaon_Pos"), eta, momentum);
804829
ue.fill(HIST("ParticleFractions/hFractionVsPt_Kaon_Pos"), eta, track.pt());
830+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_K_Pos"), eta, track.pt(), momentum);
805831
} else {
806832
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Kaon_Neg"), eta, momentum);
807833
ue.fill(HIST("ParticleFractions/hFractionVsPt_Kaon_Neg"), eta, track.pt());
834+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_K_Neg"), eta, track.pt(), momentum);
808835
}
809836
} else {
810837
ue.fill(HIST("Inclusive/hPtSecReco"), track.pt());
@@ -818,10 +845,12 @@ void MultiplicityPt::processMC(TrackTableMC const& tracks,
818845
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Proton_Pos"), eta, momentum);
819846
ue.fill(HIST("ParticleFractions/hFractionVsPt_Proton_Pos"), eta, track.pt());
820847
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_Pr_v0_Pos"), momentum, tpcSignal, eta);
848+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_Pr_Pos"), eta, track.pt(), momentum);
821849
} else {
822850
ue.fill(HIST("ParticleFractions/hFractionVsMomentum_Proton_Neg"), eta, momentum);
823851
ue.fill(HIST("ParticleFractions/hFractionVsPt_Proton_Neg"), eta, track.pt());
824852
ue.fill(HIST("DedxVsMomentum/dEdx_vs_Momentum_Pr_v0_Neg"), momentum, tpcSignal, eta);
853+
ue.fill(HIST("ResponseMatrix/heta_vs_pt_vs_p_Pr_Neg"), eta, track.pt(), momentum);
825854
}
826855
} else {
827856
ue.fill(HIST("Inclusive/hPtSecReco"), track.pt());

0 commit comments

Comments
 (0)