Skip to content

Commit ea71f25

Browse files
authored
Minor changes to debug histograms (#4394)
1 parent 278fa62 commit ea71f25

File tree

2 files changed

+9
-7
lines changed

2 files changed

+9
-7
lines changed

Analysis/Tasks/PWGHF/HFCandidateCreator2Prong.cxx

Lines changed: 0 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -41,9 +41,6 @@ struct HFCandidateCreator2Prong {
4141
Configurable<double> d_minparamchange{"d_minparamchange", 1.e-3, "stop iterations if largest change of any X is smaller than this"};
4242
Configurable<double> d_minrelchi2change{"d_minrelchi2change", 0.9, "stop iterations is chi2/chi2old > this"};
4343
Configurable<bool> b_dovalplots{"b_dovalplots", true, "do validation plots"};
44-
OutputObj<TH1F> hvtx_x_out{TH1F("hvtx_x", "2-track vtx", 100, -0.1, 0.1)};
45-
OutputObj<TH1F> hvtx_y_out{TH1F("hvtx_y", "2-track vtx", 100, -0.1, 0.1)};
46-
OutputObj<TH1F> hvtx_z_out{TH1F("hvtx_z", "2-track vtx", 100, -0.1, 0.1)};
4744
OutputObj<TH1F> hmass2{TH1F("hmass2", "2-track inv mass", 500, 0., 5.0)};
4845

4946
double massPi = RecoDecay::getMassPDG(kPiPlus);
@@ -119,9 +116,6 @@ struct HFCandidateCreator2Prong {
119116

120117
// fill histograms
121118
if (b_dovalplots) {
122-
hvtx_x_out->Fill(secondaryVertex[0]);
123-
hvtx_y_out->Fill(secondaryVertex[1]);
124-
hvtx_z_out->Fill(secondaryVertex[2]);
125119
hmass2->Fill(massPiK);
126120
hmass2->Fill(massKPi);
127121
}

Analysis/Tasks/PWGHF/HFTrackIndexSkimsCreator.cxx

Lines changed: 9 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -85,8 +85,12 @@ struct HFTrackIndexSkimsCreator {
8585
Configurable<double> d_minmassDp{"d_minmassDp", 1.5, "min mass dplus presel"};
8686
Configurable<double> d_maxmassDp{"d_maxmassDp", 2.1, "max mass dplus presel"};
8787
Configurable<bool> b_dovalplots{"b_dovalplots", true, "do validation plots"};
88+
89+
OutputObj<TH1F> hvtx_x_out{TH1F("hvtx_x", "2-track vtx", 1000, -2.0, 2.0)};
90+
OutputObj<TH1F> hvtx_y_out{TH1F("hvtx_y", "2-track vtx", 1000, -2.0, 2.0)};
91+
OutputObj<TH1F> hvtx_z_out{TH1F("hvtx_z", "2-track vtx", 1000, -10.0, 10.0)};
8892
OutputObj<TH1F> hmass2{TH1F("hmass2", "; Inv Mass (GeV/c^{2})", 500, 0, 5.0)};
89-
OutputObj<TH1F> hmass3{TH1F("hmass3", "; Inv Mass (GeV/c^{2})", 500, 0, 5.0)};
93+
OutputObj<TH1F> hmass3{TH1F("hmass3", "; Inv Mass (GeV/c^{2})", 500, 1.6, 2.1)};
9094

9195
Filter filterSelectTracks = aod::seltrack::issel == 1;
9296
using SelectedTracks = soa::Filtered<soa::Join<aod::Tracks, aod::TracksCov, aod::TracksExtra, aod::SelTrack>>;
@@ -160,6 +164,7 @@ struct HFTrackIndexSkimsCreator {
160164
array<float, 3> pvec1;
161165
df.getTrack(0).getPxPyPzGlo(pvec0);
162166
df.getTrack(1).getPxPyPzGlo(pvec1);
167+
const auto& secondaryVertex = df.getPCACandidate();
163168

164169
// calculate invariant masses
165170
auto arrMom = array{pvec0, pvec1};
@@ -169,6 +174,9 @@ struct HFTrackIndexSkimsCreator {
169174
if (b_dovalplots) {
170175
hmass2->Fill(mass2PiK);
171176
hmass2->Fill(mass2KPi);
177+
hvtx_x_out->Fill(secondaryVertex[0]);
178+
hvtx_y_out->Fill(secondaryVertex[1]);
179+
hvtx_z_out->Fill(secondaryVertex[2]);
172180
}
173181

174182
// fill table row

0 commit comments

Comments
 (0)