From 78dda9ef2d52254c659a36b7efcb0de3ab3cd6b5 Mon Sep 17 00:00:00 2001 From: Adam Matyja Date: Thu, 5 Nov 2026 15:46:16 +0100 Subject: [PATCH 1/2] Control histograms added --- .../tauThreeProngEventTableProducer.cxx | 57 ++++++++++++++++--- 1 file changed, 50 insertions(+), 7 deletions(-) diff --git a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx index 37583d9a938..d99ab51f4c8 100644 --- a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx +++ b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx @@ -260,9 +260,6 @@ struct TauThreeProngEventTableProducer { registrySkim.add("skim/gapSide", ";Gap;events", {HistType::kTH1D, {{10, -1., 9.}}}); registrySkim.add("skim/trueGapSide", ";TrueGap;events", {HistType::kTH1D, {{10, -1., 9.}}}); - registrySkim.add("skim/etaTrk", ";#eta^{trk};events", {HistType::kTH1F, {{100, -1.5, 1.5}}}); - registrySkim.add("skim/ptTrk", ";p_{T}^{trk};events", {HistType::kTH1F, {{100, 0., 5.}}}); - registrySkim.add("skim/phiTrk", ";#phi^{trk};events", {HistType::kTH1F, {{128, -3.2, 3.2}}}); registrySkim.add("skim/nTof", ";N_{TOFtrk};events", {HistType::kTH1F, {{10, -1., 9.}}}); } if (doprocessMonteCarlo) { @@ -298,6 +295,23 @@ struct TauThreeProngEventTableProducer { registrySkim.add("skim/daughterEtaMC", ";#eta^{daughter};events", {HistType::kTH1F, {{100, -4., 4.}}}); registrySkim.add("skim/daughterPtMC", ";p_{T}^{daughter};events", {HistType::kTH1F, {{100, 0, 5.0}}}); } + // common histos for doprocessDoSkim and doprocessMonteCarlo + if (doprocessDoSkim || doprocessMonteCarlo) { + registrySkim.add("rec/zvertex_SGP", ";#it{V}_{z};events", {HistType::kTH1D, {{100, -15., 15.}}}); + registrySkim.add("rec/nPVtrack_SGP", ";#it{N}_{PVtrack};events", {HistType::kTH1D, {{100, 0., 100.}}}); + registrySkim.add("rec/qtot_SGP", ";#it{Q}_{tot};events", {HistType::kTH1D, {{30, -15., 15.}}}); + // registrySkim.add("skim/etaTrk", ";#eta^{trk};events", {HistType::kTH1F, {{100, -1.5, 1.5}}}); + // registrySkim.add("skim/ptTrk", ";p_{T}^{trk};events", {HistType::kTH1F, {{100, 0., 5.}}}); + // registrySkim.add("skim/phiTrk", ";#phi^{trk};events", {HistType::kTH1F, {{128, -3.2, 3.2}}}); + registrySkim.add("rec/phiPVtrack_SGP", ";#phi^{PVtrack};events", {HistType::kTH1F, {{128, -3.2, 3.2}}}); + registrySkim.add("rec/etaPVtrack_SGP", ";#eta^{PVtrack};events", {HistType::kTH1F, {{100, -1.5, 1.5}}}); + registrySkim.add("rec/ptPVtrack_SGP", ";#it{p}_{T}^{PVtrack};events", {HistType::kTH1F, {{100, 0, 5.0}}}); + registrySkim.add("rec/nTOFtrack_SGP", ";#it{N}_{TOFtrack};events", {HistType::kTH1F, {{50, 0, 50}}}); + registrySkim.add("rec/FT0A_SGP", ";#it{E}_{FT0A};events", {HistType::kTH1F, {{100, 0, 500}}}); + registrySkim.add("rec/FT0C_SGP", ";#it{E}_{FT0C};events", {HistType::kTH1F, {{100, 0, 500}}}); + registrySkim.add("rec/FV0A_SGP", ";#it{E}_{FV0A};events", {HistType::kTH1F, {{100, 0, 500}}}); + } + if (doprocessGenerated) { registrySkim.add("gen/efficiencyMC", ";efficiency;events", {HistType::kTH1D, {{16, -1., 15.}}}); registrySkim.get(HIST("gen/efficiencyMC"))->GetXaxis()->SetBinLabel(1, "All"); @@ -762,14 +776,15 @@ struct TauThreeProngEventTableProducer { int nTofTrk = 0; int nEtaIn15 = 0; int npT100 = 0; + int qtot = 0; ROOT::Math::LorentzVector> p; for (const auto& trk : PVContributors) { - // qtot += trk.sign(); + qtot += trk.sign(); // p.SetXYZM(trk.px(), trk.py(), trk.pz(), MassPiPlus); p.SetXYZT(trk.px(), trk.py(), trk.pz(), RecoDecay::e(trk.px(), trk.py(), trk.pz(), MassPiPlus)); - registrySkim.get(HIST("skim/etaTrk"))->Fill(p.Eta()); - registrySkim.get(HIST("skim/ptTrk"))->Fill(trk.pt()); - registrySkim.get(HIST("skim/phiTrk"))->Fill(p.Phi()); + registrySkim.get(HIST("rec/etaPVtrack_SGP"))->Fill(p.Eta()); + registrySkim.get(HIST("rec/ptPVtrack_SGP"))->Fill(trk.pt()); + registrySkim.get(HIST("rec/phiPVtrack_SGP"))->Fill(p.Phi()); if (std::abs(p.Eta()) < trkEtacut) nEtaIn15++; // 0.9 is a default @@ -779,6 +794,20 @@ struct TauThreeProngEventTableProducer { nTofTrk++; } // end of loop over PV tracks + // + // plots after Single Gap producer without any further cuts + // + registrySkim.get(HIST("rec/zvertex_SGP"))->Fill(dgcand.posZ()); + registrySkim.get(HIST("rec/nPVtrack_SGP"))->Fill(PVContributors.size()); + registrySkim.get(HIST("rec/qtot_SGP"))->Fill(qtot); + registrySkim.get(HIST("rec/nTOFtrack_SGP"))->Fill(nTofTrk); + registrySkim.get(HIST("rec/FT0A_SG"))->Fill(dgcand.totalFT0AmplitudeA()); + registrySkim.get(HIST("rec/FT0C_SGP"))->Fill(dgcand.totalFT0AmplitudeC()); + registrySkim.get(HIST("rec/FV0A_SGP"))->Fill(dgcand.totalFV0AmplitudeA()); + // + // end of plots after Single Gap producer + // + // check number of PV tracks if (PVContributors.size() != nPVtrackscut) // 4 or 6 return; @@ -1596,6 +1625,20 @@ struct TauThreeProngEventTableProducer { // timesFIT[1] = collFromMcColl.timeFT0C(); // timesFIT[2] = collFromMcColl.timeFV0A(); + // + // plots after Single Gap producer without any further cuts + // + registrySkim.get(HIST("rec/zvertex_SGP"))->Fill(zVertex); + registrySkim.get(HIST("rec/nPVtrack_SGP"))->Fill(nPVcontrib); + // registrySkim.get(HIST("rec/qtot_SGP"))->Fill(qtot); + // registrySkim.get(HIST("rec/nTOFtrack_SGP"))->Fill(nTofTrk); + registrySkim.get(HIST("rec/FT0A_SG"))->Fill(amplitudesFIT[0]); + registrySkim.get(HIST("rec/FT0C_SGP"))->Fill(amplitudesFIT[1]); + registrySkim.get(HIST("rec/FV0A_SGP"))->Fill(amplitudesFIT[2]); + // + // end of plots after Single Gap producer + // + // get particles associated to generated collision auto const& partsFromMcColl = mcParticles.sliceBy(partPerMcCollision, mccoll.globalIndex()); if (verbose) From 8832be3a0c1aa23c135e32bd7625eccd39654569 Mon Sep 17 00:00:00 2001 From: Adam Matyja Date: Thu, 5 Nov 2026 15:58:55 +0100 Subject: [PATCH 2/2] Fix typo in histogram name --- PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx index d99ab51f4c8..8afc2b486c8 100644 --- a/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx +++ b/PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx @@ -801,7 +801,7 @@ struct TauThreeProngEventTableProducer { registrySkim.get(HIST("rec/nPVtrack_SGP"))->Fill(PVContributors.size()); registrySkim.get(HIST("rec/qtot_SGP"))->Fill(qtot); registrySkim.get(HIST("rec/nTOFtrack_SGP"))->Fill(nTofTrk); - registrySkim.get(HIST("rec/FT0A_SG"))->Fill(dgcand.totalFT0AmplitudeA()); + registrySkim.get(HIST("rec/FT0A_SGP"))->Fill(dgcand.totalFT0AmplitudeA()); registrySkim.get(HIST("rec/FT0C_SGP"))->Fill(dgcand.totalFT0AmplitudeC()); registrySkim.get(HIST("rec/FV0A_SGP"))->Fill(dgcand.totalFV0AmplitudeA()); //