Skip to content

Commit a3fa289

Browse files
monamelopmonamelop
committed
[PWGJE] improve QA histograms in jetDsSpecSubs
1 parent 94f2054 commit a3fa289

1 file changed

Lines changed: 83 additions & 87 deletions

File tree

‎PWGJE/Tasks/jetDsSpecSubs.cxx‎

Lines changed: 83 additions & 87 deletions
Original file line numberDiff line numberDiff line change
@@ -111,34 +111,29 @@ struct JetDsSpecSubs {
111111
"registry",
112112
{
113113
{"h_collisions", "event status;event status;entries", {HistType::kTH1F, {{10, 0.0, 10.0}}}},
114-
{"h_collision_counter_data", ";event counter;entries", {HistType::kTH1F, {{10, 0., 10.}}}},
114+
{"h_event_counter_data", ";Selection step;Events", {HistType::kTH1F, {{3,0.5,3.5}}}},
115115

116116
{"h_track_pt", ";#it{p}_{T,track};entries", {HistType::kTH1F, {{200, 0., 200.}}}},
117117
{"h_track_eta", ";#eta_{track};entries", {HistType::kTH1F, {{100, -1., 1.}}}},
118118
{"h_track_phi", ";#varphi_{track};entries", {HistType::kTH1F, {{80, -1., 7.}}}},
119119

120120
// Data histograms
121-
{"h_dsjet_counter_data", ";type;counts", {HistType::kTH1F, {{3, 0., 3.}}}},
122-
123121
{"h_jet_pt_data", "jet pT;#it{p}_{T,jet} (GeV/#it{c});entries", {HistType::kTH1F, {{200, 0., 200.}}}},
124122
{"h_jet_eta_data", "jet #eta;#eta_{jet};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}}},
125123
{"h_jet_phi_data", "jet #phi;#phi_{jet};entries", {HistType::kTH1F, {{80, -1.0, 7.}}}},
126124

127-
{"h_ds_mass_data", ";m_{D_{S}} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{1000, 0., 6.}}}},
128-
{"h_ds_pt_data", ";#it{p}_{T,D_{S}} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
129-
{"h_ds_eta_data", ";#eta_{D_{S}};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
130-
{"h_ds_phi_data", ";#phi_{D_{S}};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
131-
132-
{"h_ds_jet_pt_data", ";#it{p}_{T,D_{S} jet} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
133-
{"h_ds_jet_eta_data", ";#eta_{D_{S} jet};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
134-
{"h_ds_jet_phi_data", ";#phi_{D_{S} jet};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
125+
{"h_ds_mass_data", ";m_{D_{S}} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{300, 1.7, 2.15}}}},
126+
{"h_ds_pt_data", ";#it{p}_{T,D_{S}} (GeV/#it{c});entries", {HistType::kTH1F, {{250, 0., 100.}}}},
127+
{"h_ds_eta_data", ";#eta_{D_{S}};entries", {HistType::kTH1F, {{100, -1., 1.}}}},
128+
{"h_ds_phi_data", ";#phi_{D_{S}};entries", {HistType::kTH1F, {{80, -1., 7.}}}},
135129

136-
{"h_ds_jet_projection_data", ";z^{D_{S},jet}_{||};entries", {HistType::kTH1F, {{1000, 0., 2.}}}},
137-
{"h_ds_jet_distance_data", ";#DeltaR_{D_{S},jet};entries", {HistType::kTH1F, {{1000, 0., 1.}}}},
138-
{"h_ds_jet_mass_data", ";m_{jet}^{ch} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{200, 0., 50.}}}},
139-
{"h_ds_jet_lambda11_data", ";#lambda_{1}^{1};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
140-
{"h_ds_jet_lambda12_data", ";#lambda_{2}^{1};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
141-
{"hSparse_ds_data", ";m_{D_{S}};#it{p}_{T,D_{S}};#it{p}_{T,jet};z^{D_{S},jet}_{||};#DeltaR_{D_{S},jet};Origin(D_{S})", {HistType::kTHnSparseF, {{60, 1.7, 2.1}, {60, 0., 100.}, {60, 0., 100.}, {60, 0., 2.}, {60, 0., 1.0}}}},
130+
{"h_ds_jet_projection_data", ";z^{D_{S},jet}_{||};entries", {HistType::kTH1F, {{200, 0., 1.2}}}},
131+
{"h_ds_jet_distance_data", ";#DeltaR_{D_{S},jet};entries", {HistType::kTH1F, {{200, 0., 1.}}}},
132+
{"h_ds_jet_mass_data", ";m_{jet}^{ch} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{300, 0., 25.}}}},
133+
{"h_ds_jet_lambda11_data", ";#lambda_{1}^{1};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
134+
{"h_ds_jet_lambda12_data", ";#lambda_{2}^{1};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
135+
{"hSparse_ds_data", ";m_{D_{S}};#it{p}_{T,D_{S}};#it{p}_{T,jet};z^{D_{S},jet}_{||};#DeltaR_{D_{S},jet}",
136+
{HistType::kTHnSparseF, {{60, 1.7, 2.15}, {60, 0., 80.}, {60, 0., 100.}, {20, 0., 1.2}, {20, 0., 1.0}}}},
142137

143138
// MC general histograms
144139
{"McEffJet", "N_{jet};", {HistType::kTH1F, {{4, 0., 4.0}}}},
@@ -149,38 +144,39 @@ struct JetDsSpecSubs {
149144
{"h_jet_eta_mcd", "detector-level jet #eta;#eta_{jet}^{det};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}}},
150145
{"h_jet_phi_mcd", "detector-level jet #phi;#phi_{jet}^{det};entries", {HistType::kTH1F, {{80, -1.0, 7.}}}},
151146

152-
{"h_ds_mass_mcd", ";m_{D_{S}}^{rec} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{1000, 0., 6.}}}},
153-
{"h_ds_pt_mcd", ";#it{p}_{T,D_{S}}^{det} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
154-
{"h_ds_eta_mcd", ";#eta_{D_{S}}^{det};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
155-
{"h_ds_phi_mcd", ";#phi_{D_{S}}^{det};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
156-
157-
{"h_ds_jet_pt_mcd", ";#it{p}_{T,D_{S} jet}^{det} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
158-
{"h_ds_jet_eta_mcd", ";#eta_{D_{S} jet}^{det};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
159-
{"h_ds_jet_phi_mcd", ";#phi_{D_{S} jet}^{det};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
160-
{"h_ds_jet_projection_mcd", ";z^{D_{S},jet}_{||, det};entries", {HistType::kTH1F, {{1000, 0., 2.}}}},
161-
{"h_ds_jet_distance_mcd", ";#DeltaR_{D_{S},jet}^{det};entries", {HistType::kTH1F, {{1000, 0., 1.}}}},
162-
{"h_ds_jet_mass_mcd", ";m_{jet}^{ch, det} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{200, 0., 50.}}}},
163-
{"h_ds_jet_lambda11_mcd", ";#lambda_{1}^{1, det};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
164-
{"h_ds_jet_lambda12_mcd", ";#lambda_{2}^{1, det};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
165-
{"hSparse_ds_mcd", ";m_{D_{S}}^{rec};#it{p}_{T,D_{S}}^{det};#it{p}_{T,jet}^{det};z^{D_{S},jet}_{||, det};#DeltaR_{D_{S},jet}^{det}", {HistType::kTHnSparseF, {{60, 1.7, 2.1}, {60, 0., 100.}, {60, 0., 100.}, {60, 0., 2.}, {60, 0., 1.0}, {2, -0.5, 1.5}}}},
147+
{"h_ds_pt_mcd", ";#it{p}_{T,D_{S} jet}^{det} (GeV/#it{c});entries", {HistType::kTH1F, {{250, 0., 100.}}}},
148+
{"h_ds_eta_mcd", ";#eta_{D_{S} jet}^{det};entries", {HistType::kTH1F, {{100, -1., 1.}}}},
149+
{"h_ds_phi_mcd", ";#phi_{D_{S} jet}^{det};entries", {HistType::kTH1F, {{80, -1., 7.}}}},
150+
{"h_ds_mass_mcd", ";m_{D_{S}}^{det} (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{300, 1.7, 2.15}}}},
151+
152+
{"h_ds_jet_lambda11_mcd", ";#lambda_{1}^{1, det};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
153+
{"h_ds_jet_lambda12_mcd", ";#lambda_{2}^{1, det};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
154+
155+
// MCD - Sparse 1: mass, p_{T,Ds}, p_{T,jet}, z|| and prompt/non-prompt
156+
{"hSparse_ds_mcd1", ";m_{D_{S}}^{rec};#it{p}_{T,D_{S}}^{det};#it{p}_{T,jet}^{det};z^{D_{S},jet}_{||, det};Origin(D_{S})",
157+
{HistType::kTHnSparseF, {{60, 1.7, 2.15}, {60, 0., 80.}, {60, 0., 100.}, {20, 0., 1.2}, {2, -0.5, 1.5}}}},
158+
// MCD - Sparse 2: p_{T,Ds}, p_{T,jet}, and DeltaR
159+
{"hSparse_ds_mcd2", ";#it{p}_{T,D_{S}}^{det};#it{p}_{T,jet}^{det};#DeltaR_{D_{S},jet}^{det}",
160+
{HistType::kTHnSparseF, {{60, 0., 80.}, {60, 0., 100.}, {20, 0., 1.0}}}},
161+
// MCD - Sparse 3: p_{T,jet}, z|| and DeltaR
162+
{"hSparse_ds_mcd3", ";#it{p}_{T,jet}^{det};z^{D_{S},jet}_{||, det};#DeltaR_{D_{S},jet}^{det}",
163+
{HistType::kTHnSparseF, {{60, 0., 100.}, {20, 0., 1.2}, {20, 0., 1.0}}}},
166164

167165
// MC particle-level histograms
168166
{"h_jet_pt_mcp", "particle-level jet pT;#it{p}_{T,jet}^{part} (GeV/#it{c});entries", {HistType::kTH1F, {{200, 0., 200.}}}},
169167
{"h_jet_eta_mcp", "particle-level jet #eta;#eta_{jet}^{part};entries", {HistType::kTH1F, {{100, -1.0, 1.0}}}},
170168
{"h_jet_phi_mcp", "particle-level jet #phi;#phi_{jet}^{part};entries", {HistType::kTH1F, {{80, -1.0, 7.}}}},
171169

172-
{"h_ds_pt_mcp", ";#it{p}_{T,D_{S}}^{part} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
173-
{"h_ds_eta_mcp", ";#eta_{D_{S}}^{part};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
174-
{"h_ds_phi_mcp", ";#phi_{D_{S}}^{part};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
175-
176-
{"h_ds_jet_pt_mcp", ";#it{p}_{T,D_{S} jet}^{part} (GeV/#it{c});entries", {HistType::kTH1F, {{1000, 0., 100.}}}},
177-
{"h_ds_jet_eta_mcp", ";#eta_{D_{S} jet}^{part};entries", {HistType::kTH1F, {{250, -1., 1.}}}},
178-
{"h_ds_jet_phi_mcp", ";#phi_{D_{S} jet}^{part};entries", {HistType::kTH1F, {{250, -1., 7.}}}},
179-
{"h_ds_jet_projection_mcp", ";z^{D_{S},jet}_{||, part};entries", {HistType::kTH1F, {{1000, 0., 2.}}}},
180-
{"h_ds_jet_distance_mcp", ";#DeltaR_{D_{S},jet}^{part};entries", {HistType::kTH1F, {{1000, 0., 1.}}}},
181-
{"h_ds_jet_lambda11_mcp", ";#lambda_{1}^{1, part};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
182-
{"h_ds_jet_lambda12_mcp", ";#lambda_{2}^{1, part};entries", {HistType::kTH1F, {{200, 0., 1.0}}}},
183-
{"hSparse_ds_mcp", ";#it{p}_{T,D_{S}}^{part};#it{p}_{T,jet}^{part};z^{D_{S},jet}_{||, part};#DeltaR_{D_{S},jet}^{part}", {HistType::kTHnSparseF, {{60, 0., 100.}, {60, 0., 100.}, {60, 0., 2.}, {60, 0., 1.0}}}},
170+
{"h_ds_pt_mcp", ";#it{p}_{T,D_{S} jet}^{part} (GeV/#it{c});entries", {HistType::kTH1F, {{250, 0., 100.}}}},
171+
{"h_ds_eta_mcp", ";#eta_{D_{S} jet}^{part};entries", {HistType::kTH1F, {{100, -1., 1.}}}},
172+
{"h_ds_phi_mcp", ";#phi_{D_{S} jet}^{part};entries", {HistType::kTH1F, {{80, -1., 7.}}}},
173+
174+
{"h_ds_jet_lambda11_mcp", ";#lambda_{1}^{1, part};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
175+
{"h_ds_jet_lambda12_mcp", ";#lambda_{2}^{1, part};entries", {HistType::kTH1F, {{100, 0., 1.0}}}},
176+
177+
// MCP - Sparse: p_{T,Ds}, p_{T,jet}, z|| and DeltaR
178+
{"hSparse_ds_mcp", ";#it{p}_{T,D_{S}}^{part};#it{p}_{T,jet}^{part};z^{D_{S},jet}_{||, part};#DeltaR_{D_{S},jet}^{part}",
179+
{HistType::kTHnSparseF, {{60, 0., 80.}, {60, 0., 100.}, {20, 0., 1.2}, {20, 0., 1.0}}}},
184180
}};
185181
//========
186182
// INIT
@@ -191,10 +187,10 @@ struct JetDsSpecSubs {
191187
eventSelectionBits = jetderiveddatautilities::initialiseEventSelectionBits(static_cast<std::string>(eventSelections));
192188
trackSelection = jetderiveddatautilities::initialiseTrackSelection(static_cast<std::string>(trackSelections));
193189

194-
auto hData = registry.get<TH1>(HIST("h_dsjet_counter_data"));
195-
hData->GetXaxis()->SetBinLabel(1, "Ds-jet entries");
196-
hData->GetXaxis()->SetBinLabel(2, "Ds candidates");
197-
hData->GetXaxis()->SetBinLabel(3, "Ds jets with >=1 cand.");
190+
auto hEvt = registry.get<TH1>(HIST("h_event_counter_data"));
191+
hEvt->GetXaxis()->SetBinLabel(1,"Input collisions");
192+
hEvt->GetXaxis()->SetBinLabel(2,"Event selection");
193+
hEvt->GetXaxis()->SetBinLabel(3,"|z| < 10 cm");
198194

199195
// Labels
200196
auto mcCollisionCounter = registry.get<TH1>(HIST("McEffCol"));
@@ -209,8 +205,8 @@ struct JetDsSpecSubs {
209205
jetCounter->GetXaxis()->SetBinLabel(BinMCJetCntr::DetectorLevelJetWithMatchedCandidate, "particle matched jets");
210206
jetCounter->GetXaxis()->SetBinLabel(BinMCJetCntr::ParticleLevelJetWithMatchedCandidate, "detector matched jets");
211207

212-
auto hSparse_ds_mcd = registry.get<THnSparse>(HIST("hSparse_ds_mcd"));
213-
auto* axisOrigin = hSparse_ds_mcd->GetAxis(5);
208+
auto hSparse_ds_mcd1 = registry.get<THnSparse>(HIST("hSparse_ds_mcd1"));
209+
auto* axisOrigin = hSparse_ds_mcd1->GetAxis(4);
214210
axisOrigin->SetBinLabel(1, "Prompt");
215211
axisOrigin->SetBinLabel(2, "Non-prompt");
216212
}
@@ -302,17 +298,19 @@ struct JetDsSpecSubs {
302298
void processDataChargedSubstructure(aod::JetCollision const& collision, FilteredDsDataJets const& jets,
303299
aod::CandidatesDsData const&, aod::JetTracks const&)
304300
{
305-
registry.fill(HIST("h_collision_counter_data"), 2.0);
301+
registry.fill(HIST("h_event_counter_data"), 1);
306302

307-
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits) ||
308-
!(std::abs(collision.posZ()) < vertexZCut)) {
303+
if (!jetderiveddatautilities::selectCollision(collision, eventSelectionBits)) {
309304
return;
310305
}
306+
registry.fill(HIST("h_event_counter_data"), 2);
311307

312-
registry.fill(HIST("h_collision_counter_data"), 3.0);
308+
if (std::abs(collision.posZ()) >= vertexZCut) {
309+
return;
310+
}
311+
registry.fill(HIST("h_event_counter_data"), 3);
313312

314313
for (const auto& jet : jets) {
315-
registry.fill(HIST("h_dsjet_counter_data"), 0.5); // DsChargedJets entries
316314

317315
registry.fill(HIST("h_jet_pt_data"), jet.pt());
318316
registry.fill(HIST("h_jet_eta_data"), jet.eta());
@@ -327,12 +325,8 @@ struct JetDsSpecSubs {
327325

328326
TVector3 jetVector(jet.px(), jet.py(), jet.pz());
329327

330-
int nDsInJet = 0;
331-
332328
// Loop over Ds candidates (particle level)
333329
for (const auto& dsCandidate : jet.candidates_as<aod::CandidatesDsData>()) {
334-
++nDsInJet;
335-
registry.fill(HIST("h_dsjet_counter_data"), 1.5); // Ds candidates associated with the jet
336330

337331
TVector3 dsVector(dsCandidate.px(), dsCandidate.py(), dsCandidate.pz());
338332

@@ -351,6 +345,7 @@ struct JetDsSpecSubs {
351345

352346
registry.fill(HIST("h_ds_jet_distance_data"), deltaR);
353347

348+
354349
// Main THnSparse: invariant mass, pT, z, and DeltaR
355350
registry.fill(HIST("hSparse_ds_data"),
356351
dsCandidate.m(),
@@ -360,25 +355,17 @@ struct JetDsSpecSubs {
360355
deltaR);
361356
}
362357

363-
// Jet-level quantities (filled once per jet containing at least one Ds)
364-
if (nDsInJet > 0) {
365-
366-
registry.fill(HIST("h_dsjet_counter_data"), 2.5); // Ds jets with at least one associated candidate
358+
if (!jet.candidates_as<aod::CandidatesDsData>().empty()) {
359+
// Jet mass
360+
registry.fill(HIST("h_ds_jet_mass_data"),mjet);
367361

368-
// Jet properties
369-
registry.fill(HIST("h_ds_jet_pt_data"), jet.pt());
370-
registry.fill(HIST("h_ds_jet_eta_data"), jet.eta());
371-
registry.fill(HIST("h_ds_jet_phi_data"), jet.phi());
372-
// Jet Mass
373-
registry.fill(HIST("h_ds_jet_mass_data"), mjet);
374-
375-
// Jet substructure observables
376-
if (lambda11 >= 0.f) {
377-
registry.fill(HIST("h_ds_jet_lambda11_data"), lambda11);
378-
}
379-
if (lambda12 >= 0.f) {
380-
registry.fill(HIST("h_ds_jet_lambda12_data"), lambda12);
381-
}
362+
// Jet angularity
363+
if (lambda11 >= 0.f) {
364+
registry.fill(HIST("h_ds_jet_lambda11_data"), lambda11);
365+
}
366+
if (lambda12 >= 0.f) {
367+
registry.fill(HIST("h_ds_jet_lambda12_data"), lambda12);
368+
}
382369
}
383370
}
384371
}
@@ -464,19 +451,28 @@ struct JetDsSpecSubs {
464451
registry.fill(HIST("h_ds_jet_lambda11_mcd"), mcd_lambda11);
465452
registry.fill(HIST("h_ds_jet_lambda12_mcd"), mcd_lambda12);
466453
// Detector-level Ds Histgrams
467-
registry.fill(HIST("h_ds_jet_pt_mcd"), mcdDscand.pt());
468-
registry.fill(HIST("h_ds_jet_mass_mcd"), mcdDscand.m());
469-
registry.fill(HIST("h_ds_jet_eta_mcd"), mcdDscand.eta());
470-
registry.fill(HIST("h_ds_jet_phi_mcd"), mcdDscand.phi());
454+
registry.fill(HIST("h_ds_pt_mcd"), mcdDscand.pt());
455+
registry.fill(HIST("h_ds_mass_mcd"), mcdDscand.m());
456+
registry.fill(HIST("h_ds_eta_mcd"), mcdDscand.eta());
457+
registry.fill(HIST("h_ds_phi_mcd"), mcdDscand.phi());
471458

472-
// Main THnSparse: invariant mass, pT, z, and DeltaR
473-
registry.fill(HIST("hSparse_ds_mcd"),
459+
// MCD THnSparse1: invariant mass, p{T,Ds}, pT, z, and origin (prompt/non-prompt)
460+
registry.fill(HIST("hSparse_ds_mcd1"),
474461
mcdDscand.m(),
475462
mcdDscand.pt(),
476463
mcdjet.pt(),
477464
mcd_zParallel,
478-
mcd_deltaR,
479465
origin);
466+
// MCD THnSparse2: invariant p{T,Ds}, pT and DeltaR
467+
registry.fill(HIST("hSparse_ds_mcd2"),
468+
mcdDscand.pt(),
469+
mcdjet.pt(),
470+
mcd_deltaR);
471+
// MCD THnSparse3: invariant pT z and DeltaR
472+
registry.fill(HIST("hSparse_ds_mcd3"),
473+
mcdjet.pt(),
474+
mcd_zParallel,
475+
mcd_deltaR);
480476
}
481477
}
482478
// Particle level
@@ -505,13 +501,13 @@ struct JetDsSpecSubs {
505501
registry.fill(HIST("h_jet_pt_mcp"), mcpjet.pt());
506502
registry.fill(HIST("h_jet_eta_mcp"), mcpjet.eta());
507503
registry.fill(HIST("h_jet_phi_mcp"), mcpjet.phi());
508-
registry.fill(HIST("h_ds_jet_projection_mcp"), mcp_zParallel);
504+
509505
registry.fill(HIST("h_ds_jet_lambda11_mcp"), mcp_lambda11);
510506
registry.fill(HIST("h_ds_jet_lambda12_mcp"), mcp_lambda12);
511507
// Particle-level Ds Histgrams
512-
registry.fill(HIST("h_ds_jet_pt_mcp"), mcpDscand.pt());
513-
registry.fill(HIST("h_ds_jet_eta_mcp"), mcpDscand.eta());
514-
registry.fill(HIST("h_ds_jet_phi_mcp"), mcpDscand.phi());
508+
registry.fill(HIST("h_ds_pt_mcp"), mcpDscand.pt());
509+
registry.fill(HIST("h_ds_eta_mcp"), mcpDscand.eta());
510+
registry.fill(HIST("h_ds_phi_mcp"), mcpDscand.phi());
515511

516512
// Main THnSparse: invariant mass, pT, z, and DeltaR
517513
registry.fill(HIST("hSparse_ds_mcp"),

0 commit comments

Comments
 (0)