Skip to content

Commit 78dda9e

Browse files
amatyjaamatyja
committed
Control histograms added
1 parent cc0fcac commit 78dda9e

1 file changed

Lines changed: 50 additions & 7 deletions

File tree

‎PWGUD/TableProducer/tauThreeProngEventTableProducer.cxx‎

Lines changed: 50 additions & 7 deletions
Original file line numberDiff line numberDiff line change
@@ -260,9 +260,6 @@ struct TauThreeProngEventTableProducer {
260260

261261
registrySkim.add("skim/gapSide", ";Gap;events", {HistType::kTH1D, {{10, -1., 9.}}});
262262
registrySkim.add("skim/trueGapSide", ";TrueGap;events", {HistType::kTH1D, {{10, -1., 9.}}});
263-
registrySkim.add("skim/etaTrk", ";#eta^{trk};events", {HistType::kTH1F, {{100, -1.5, 1.5}}});
264-
registrySkim.add("skim/ptTrk", ";p_{T}^{trk};events", {HistType::kTH1F, {{100, 0., 5.}}});
265-
registrySkim.add("skim/phiTrk", ";#phi^{trk};events", {HistType::kTH1F, {{128, -3.2, 3.2}}});
266263
registrySkim.add("skim/nTof", ";N_{TOFtrk};events", {HistType::kTH1F, {{10, -1., 9.}}});
267264
}
268265
if (doprocessMonteCarlo) {
@@ -298,6 +295,23 @@ struct TauThreeProngEventTableProducer {
298295
registrySkim.add("skim/daughterEtaMC", ";#eta^{daughter};events", {HistType::kTH1F, {{100, -4., 4.}}});
299296
registrySkim.add("skim/daughterPtMC", ";p_{T}^{daughter};events", {HistType::kTH1F, {{100, 0, 5.0}}});
300297
}
298+
// common histos for doprocessDoSkim and doprocessMonteCarlo
299+
if (doprocessDoSkim || doprocessMonteCarlo) {
300+
registrySkim.add("rec/zvertex_SGP", ";#it{V}_{z};events", {HistType::kTH1D, {{100, -15., 15.}}});
301+
registrySkim.add("rec/nPVtrack_SGP", ";#it{N}_{PVtrack};events", {HistType::kTH1D, {{100, 0., 100.}}});
302+
registrySkim.add("rec/qtot_SGP", ";#it{Q}_{tot};events", {HistType::kTH1D, {{30, -15., 15.}}});
303+
// registrySkim.add("skim/etaTrk", ";#eta^{trk};events", {HistType::kTH1F, {{100, -1.5, 1.5}}});
304+
// registrySkim.add("skim/ptTrk", ";p_{T}^{trk};events", {HistType::kTH1F, {{100, 0., 5.}}});
305+
// registrySkim.add("skim/phiTrk", ";#phi^{trk};events", {HistType::kTH1F, {{128, -3.2, 3.2}}});
306+
registrySkim.add("rec/phiPVtrack_SGP", ";#phi^{PVtrack};events", {HistType::kTH1F, {{128, -3.2, 3.2}}});
307+
registrySkim.add("rec/etaPVtrack_SGP", ";#eta^{PVtrack};events", {HistType::kTH1F, {{100, -1.5, 1.5}}});
308+
registrySkim.add("rec/ptPVtrack_SGP", ";#it{p}_{T}^{PVtrack};events", {HistType::kTH1F, {{100, 0, 5.0}}});
309+
registrySkim.add("rec/nTOFtrack_SGP", ";#it{N}_{TOFtrack};events", {HistType::kTH1F, {{50, 0, 50}}});
310+
registrySkim.add("rec/FT0A_SGP", ";#it{E}_{FT0A};events", {HistType::kTH1F, {{100, 0, 500}}});
311+
registrySkim.add("rec/FT0C_SGP", ";#it{E}_{FT0C};events", {HistType::kTH1F, {{100, 0, 500}}});
312+
registrySkim.add("rec/FV0A_SGP", ";#it{E}_{FV0A};events", {HistType::kTH1F, {{100, 0, 500}}});
313+
}
314+
301315
if (doprocessGenerated) {
302316
registrySkim.add("gen/efficiencyMC", ";efficiency;events", {HistType::kTH1D, {{16, -1., 15.}}});
303317
registrySkim.get<TH1>(HIST("gen/efficiencyMC"))->GetXaxis()->SetBinLabel(1, "All");
@@ -762,14 +776,15 @@ struct TauThreeProngEventTableProducer {
762776
int nTofTrk = 0;
763777
int nEtaIn15 = 0;
764778
int npT100 = 0;
779+
int qtot = 0;
765780
ROOT::Math::LorentzVector<ROOT::Math::PxPyPzE4D<double>> p;
766781
for (const auto& trk : PVContributors) {
767-
// qtot += trk.sign();
782+
qtot += trk.sign();
768783
// p.SetXYZM(trk.px(), trk.py(), trk.pz(), MassPiPlus);
769784
p.SetXYZT(trk.px(), trk.py(), trk.pz(), RecoDecay::e(trk.px(), trk.py(), trk.pz(), MassPiPlus));
770-
registrySkim.get<TH1>(HIST("skim/etaTrk"))->Fill(p.Eta());
771-
registrySkim.get<TH1>(HIST("skim/ptTrk"))->Fill(trk.pt());
772-
registrySkim.get<TH1>(HIST("skim/phiTrk"))->Fill(p.Phi());
785+
registrySkim.get<TH1>(HIST("rec/etaPVtrack_SGP"))->Fill(p.Eta());
786+
registrySkim.get<TH1>(HIST("rec/ptPVtrack_SGP"))->Fill(trk.pt());
787+
registrySkim.get<TH1>(HIST("rec/phiPVtrack_SGP"))->Fill(p.Phi());
773788

774789
if (std::abs(p.Eta()) < trkEtacut)
775790
nEtaIn15++; // 0.9 is a default
@@ -779,6 +794,20 @@ struct TauThreeProngEventTableProducer {
779794
nTofTrk++;
780795
} // end of loop over PV tracks
781796

797+
//
798+
// plots after Single Gap producer without any further cuts
799+
//
800+
registrySkim.get<TH1>(HIST("rec/zvertex_SGP"))->Fill(dgcand.posZ());
801+
registrySkim.get<TH1>(HIST("rec/nPVtrack_SGP"))->Fill(PVContributors.size());
802+
registrySkim.get<TH1>(HIST("rec/qtot_SGP"))->Fill(qtot);
803+
registrySkim.get<TH1>(HIST("rec/nTOFtrack_SGP"))->Fill(nTofTrk);
804+
registrySkim.get<TH1>(HIST("rec/FT0A_SG"))->Fill(dgcand.totalFT0AmplitudeA());
805+
registrySkim.get<TH1>(HIST("rec/FT0C_SGP"))->Fill(dgcand.totalFT0AmplitudeC());
806+
registrySkim.get<TH1>(HIST("rec/FV0A_SGP"))->Fill(dgcand.totalFV0AmplitudeA());
807+
//
808+
// end of plots after Single Gap producer
809+
//
810+
782811
// check number of PV tracks
783812
if (PVContributors.size() != nPVtrackscut) // 4 or 6
784813
return;
@@ -1596,6 +1625,20 @@ struct TauThreeProngEventTableProducer {
15961625
// timesFIT[1] = collFromMcColl.timeFT0C();
15971626
// timesFIT[2] = collFromMcColl.timeFV0A();
15981627

1628+
//
1629+
// plots after Single Gap producer without any further cuts
1630+
//
1631+
registrySkim.get<TH1>(HIST("rec/zvertex_SGP"))->Fill(zVertex);
1632+
registrySkim.get<TH1>(HIST("rec/nPVtrack_SGP"))->Fill(nPVcontrib);
1633+
// registrySkim.get<TH1>(HIST("rec/qtot_SGP"))->Fill(qtot);
1634+
// registrySkim.get<TH1>(HIST("rec/nTOFtrack_SGP"))->Fill(nTofTrk);
1635+
registrySkim.get<TH1>(HIST("rec/FT0A_SG"))->Fill(amplitudesFIT[0]);
1636+
registrySkim.get<TH1>(HIST("rec/FT0C_SGP"))->Fill(amplitudesFIT[1]);
1637+
registrySkim.get<TH1>(HIST("rec/FV0A_SGP"))->Fill(amplitudesFIT[2]);
1638+
//
1639+
// end of plots after Single Gap producer
1640+
//
1641+
15991642
// get particles associated to generated collision
16001643
auto const& partsFromMcColl = mcParticles.sliceBy(partPerMcCollision, mccoll.globalIndex());
16011644
if (verbose)

0 commit comments

Comments
 (0)