Skip to content

Commit 55c977e

Browse files
fgrosafgrosa
committed
Minor fixes
1 parent 902fdee commit 55c977e

2 files changed

Lines changed: 14 additions & 16 deletions

File tree

‎PWGHF/D2H/TableProducer/candidateCreatorBToJpsiReduced.cxx‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -275,7 +275,7 @@ struct HfCandidateCreatorBToJpsiReduced {
275275
} else if constexpr (DecChannel == DecayChannel::B0ToJpsiK0Star) {
276276
for (const auto& trackLf1 : tracksLfDau1ThisCollision) {
277277
// this track is among daughters
278-
if (trackLf1.trackId() == candJpsi.prongPosId() || trackLf1.trackId() == candJpsi.prongNegId()) {
278+
if (trackLf1.trackId() == candJpsi.prongPosId() || trackLf1.trackId() == candJpsi.prongNegId() || trackLf1.signed1Pt() * trackLf0.signed1Pt() > 0) {
279279
continue;
280280
}
281281
auto trackParCovLf1 = getTrackParCov(trackLf1);
@@ -352,7 +352,7 @@ struct HfCandidateCreatorBToJpsiReduced {
352352
} else if constexpr (DecChannel == DecayChannel::BsToJpsiPhi) {
353353
for (const auto& trackLf1 : tracksLfDau1ThisCollision) {
354354
// this track is among daughters
355-
if (trackLf1.trackId() == candJpsi.prongPosId() || trackLf1.trackId() == candJpsi.prongNegId()) {
355+
if (trackLf1.trackId() == candJpsi.prongPosId() || trackLf1.trackId() == candJpsi.prongNegId() || trackLf1.signed1Pt() * trackLf0.signed1Pt() > 0) {
356356
continue;
357357
}
358358
auto trackParCovLf1 = getTrackParCov(trackLf1);

‎PWGHF/D2H/TableProducer/dataCreatorJpsiHadReduced.cxx‎

Lines changed: 12 additions & 14 deletions
Original file line numberDiff line numberDiff line change
@@ -230,14 +230,14 @@ struct HfDataCreatorJpsiHadReduced {
230230
labelsEvents[Event::NoCharmHadPiSelected] = "without CharmHad-Pi pairs";
231231
labelsEvents[Event::CharmHadPiSelected] = "with CharmHad-Pi pairs";
232232
static const AxisSpec axisEvents = {NumBinsEvents, 0.5, NumBinsEvents + 0.5, ""};
233-
registry.add("hEvents", "Events;;entries", HistType::kTH1F, {axisEvents});
233+
registry.add("hEvents", "Events;;entries", HistType::kTH1D, {axisEvents});
234234
for (int iBin = 0; iBin < NumBinsEvents; iBin++) {
235235
registry.get<TH1>(HIST("hEvents"))->GetXaxis()->SetBinLabel(iBin + 1, labelsEvents[iBin].data());
236236
}
237237

238-
registry.add("hMassJpsi", "J/Psi mass;#it{M}_{#mu#mu} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{600, 2.8, 3.4, "#it{p}_{T} (GeV/#it{c})"}}});
239-
registry.add("hPtJpsi", "J/Psi #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1F, {{(std::vector<double>)binsPt, "#it{p}_{T} (GeV/#it{c})"}}});
240-
registry.add("hCpaJpsi", "J/Psi cos#theta_{p};J/Psi cos#theta_{p};Counts", {HistType::kTH1F, {{200, -1., 1, "J/Psi cos#theta_{p}"}}});
238+
registry.add("hMassJpsi", "J/Psi mass;#it{M}_{#mu#mu} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{600, 2.8, 3.4, "#it{p}_{T} (GeV/#it{c})"}}});
239+
registry.add("hPtJpsi", "J/Psi #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1D, {{(std::vector<double>)binsPt, "#it{p}_{T} (GeV/#it{c})"}}});
240+
registry.add("hCpaJpsi", "J/Psi cos#theta_{p};J/Psi cos#theta_{p};Counts", {HistType::kTH1D, {{200, -1., 1, "J/Psi cos#theta_{p}"}}});
241241
std::shared_ptr<TH1> hFitCandidatesJpsi = registry.add<TH1>("hFitCandidatesJpsi", "Jpsi candidate counter", {HistType::kTH1D, {axisCands}});
242242
std::shared_ptr<TH1> hFitCandidatesBPlus = registry.add<TH1>("hFitCandidatesBPlus", "hFitCandidatesBPlus candidate counter", {HistType::kTH1D, {axisCands}});
243243
std::shared_ptr<TH1> hFitCandidatesB0 = registry.add<TH1>("hFitCandidatesB0", "hFitCandidatesB0 candidate counter", {HistType::kTH1D, {axisCands}});
@@ -247,18 +247,18 @@ struct HfDataCreatorJpsiHadReduced {
247247
setLabelHistoCands(hFitCandidatesB0);
248248
setLabelHistoCands(hFitCandidatesBS);
249249
if (doprocessJpsiKData || doprocessJpsiKMc) {
250-
registry.add("hPtKaon", "Kaon #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1F, {{100, 0., 10.}}});
251-
registry.add("hMassJpsiKaon", "J/Psi Kaon mass;#it{M}_{J/#PsiK} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{800, 4.9, 5.7}}});
250+
registry.add("hPtKaon", "Kaon #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1D, {{100, 0., 10.}}});
251+
registry.add("hMassJpsiKaon", "J/Psi Kaon mass;#it{M}_{J/#PsiK} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{800, 4.9, 5.7}}});
252252
} else if (doprocessJpsiK0StarData || doprocessJpsiK0StarMc) {
253-
registry.add("hPtK0Star", "K*0 #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1F, {{100, 0., 10.}}});
254-
registry.add("hMassK0Star", "K*0 mass;#it{M}_{#piK} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{200, 0.9, 1.2}}});
255-
registry.add("hMassJpsiK0Star", "J/Psi K*0 mass;#it{M}_{J/#PsiK*0} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{800, 4.9, 5.7}}});
253+
registry.add("hPtK0Star", "K*0 #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1D, {{100, 0., 10.}}});
254+
registry.add("hMassK0Star", "K*0 mass;#it{M}_{#piK} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{200, 0.7, 1.1}}});
255+
registry.add("hMassJpsiK0Star", "J/Psi K*0 mass;#it{M}_{J/#PsiK*0} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{800, 4.9, 5.7}}});
256256
std::shared_ptr<TH1> hFitCandidatesK0Star = registry.add<TH1>("hFitCandidatesK0Star", "K*0 candidate counter", {HistType::kTH1D, {axisCands}});
257257
setLabelHistoCands(hFitCandidatesK0Star);
258258
} else if (doprocessJpsiPhiData || doprocessJpsiPhiMc) {
259-
registry.add("hPtPhi", "Phi #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1F, {{100, 0., 10.}}});
260-
registry.add("hMassPhi", "Phi mass;#it{M}_{KK} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{200, 0.9, 1.2}}});
261-
registry.add("hMassJpsiPhi", "J/Psi Phi mass;#it{M}_{J/#Psi#phi} (GeV/#it{c}^{2});Counts", {HistType::kTH1F, {{800, 4.9, 5.7}}});
259+
registry.add("hPtPhi", "Phi #it{p}_{T};#it{p}_{T} (GeV/#it{c});Counts", {HistType::kTH1D, {{100, 0., 10.}}});
260+
registry.add("hMassPhi", "Phi mass;#it{M}_{KK} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{200, 0.9, 1.2}}});
261+
registry.add("hMassJpsiPhi", "J/Psi Phi mass;#it{M}_{J/#Psi#phi} (GeV/#it{c}^{2});Counts", {HistType::kTH1D, {{800, 4.9, 5.7}}});
262262
std::shared_ptr<TH1> hFitCandidatesPhi = registry.add<TH1>("hFitCandidatesPhi", "Phi candidate counter", {HistType::kTH1D, {axisCands}});
263263
setLabelHistoCands(hFitCandidatesPhi);
264264
}
@@ -1187,8 +1187,6 @@ struct HfDataCreatorJpsiHadReduced {
11871187
}
11881188
} // kaon loop
11891189
if (fillHfCandJpsi) { // fill Jpsi table only once per Jpsi candidate
1190-
double invMassJpsi{0.};
1191-
invMassJpsi = HfHelper::invMassJpsiToMuMu(candidate);
11921190
hfJpsi(trackPos.globalIndex(), trackNeg.globalIndex(),
11931191
indexHfReducedCollision,
11941192
candidate.xSecondaryVertex(), candidate.ySecondaryVertex(), candidate.zSecondaryVertex(),

0 commit comments

Comments
 (0)