Skip to content

Commit 2f727c3

Browse files
FIT: add output histos to FT0 QA task
Add the following to ft0Qa.cxx: - FT0 time res vs. Ncontributors - FT0 time res vs. FT0AC multiplicity - FT0 vertex - PV vs. Ncontributors - FT0 vertex - PV vs. FT0AC multiplicity Add documentation comment for FT0 time resolution
1 parent 0291c81 commit 2f727c3

File tree

2 files changed

+14
-1
lines changed

2 files changed

+14
-1
lines changed

Common/DataModel/FT0Corrected.h

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -25,7 +25,7 @@ DECLARE_SOA_DYNAMIC_COLUMN(T0ACorrectedValid, t0ACorrectedValid, //! Was T0ACorr
2525
[](float t0) -> bool { return t0 < 1e9; });
2626
DECLARE_SOA_DYNAMIC_COLUMN(T0CCorrectedValid, t0CCorrectedValid, //! Was T0CCorrected computable?
2727
[](float t0) -> bool { return t0 < 1e9; });
28-
DECLARE_SOA_DYNAMIC_COLUMN(T0resolution, t0resolution, //! Was T0CCorrected computable?
28+
DECLARE_SOA_DYNAMIC_COLUMN(T0resolution, t0resolution, //! (T0ACorrected - T0CCorrected) / 2. FT0 time resoluition = std. dev. of this quantity's distribution. On event level, handle negative values as needed.
2929
[](float t0A, float t0C) -> float { return 0.5f * (t0A - t0C); });
3030
DECLARE_SOA_DYNAMIC_COLUMN(T0ACValid, t0ACValid, //! Was T0AC computable?
3131
[](float t0a, float t0c) -> bool { return (t0a < 1e9) && (t0c < 1e9); });

Common/Tasks/ft0Qa.cxx

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -78,6 +78,10 @@ struct ft0QaTask {
7878
histos.add("hT0AC", "T0AC;T0AC time (ns);counts", kTH1F, {axisTime});
7979
histos.add("hT0res", "FT0 resolution", kTH1F, {axisColTimeRes});
8080
histos.add("hColTime", "", kTH1F, {axisTime});
81+
histos.add("hT0res_nContrib", "FT0 resolution vs. Ncontributors", kTH2F,
82+
{axisColTimeRes, axisNcontrib});
83+
histos.add("hT0res_MultT0AC", "FT0 resolution vs. T0AC multiplicity", kTH2F,
84+
{axisColTimeRes, axisMultT0AC});
8185

8286
// FT0 vertex
8387
histos.add("hT0vertex", "FT0 vertex;FT0 vertex (cm);counts", kTH1F,
@@ -103,6 +107,11 @@ struct ft0QaTask {
103107
histos.add("hPV_nContrib",
104108
"PV vs. Ncontributers;primary vertex (cm);(# contrubutors)",
105109
kTH2F, {axisVertex, axisNcontrib});
110+
histos.add("hT0vertexDiff_vs_nContrib", "FT0V - PV vs. Ncontributors;FT0 vertex - PV (cm);# contrubutors",
111+
kTH2F, {axisVertex, axisNcontrib});
112+
histos.add("hT0vertexDiff_vs_MultT0AC",
113+
"FT0V - PV vs. T0AC multiplicity;FT0 vertex - PV (cm);T0AC multiplicity (# ADC channels)",
114+
kTH2F, {axisVertex, axisMultT0AC});
106115

107116
// FT0 amplitude and multiplicity
108117
histos.add("hAmpT0A", "amplitude T0A;#ADC channels;counts", kTH1F,
@@ -363,7 +372,11 @@ struct ft0QaTask {
363372
histos.fill(HIST("hVertex_T0_PV"), ft0.posZ(), collision.posZ());
364373
histos.fill(HIST("hPV"), collision.posZ());
365374
histos.fill(HIST("hT0res"), collision.t0resolution());
375+
histos.fill(HIST("hT0res_nContrib"), collision.t0resolution(), nContrib);
376+
histos.fill(HIST("hT0res_MultT0AC"), collision.t0resolution(), multFT0M);
366377
histos.fill(HIST("hT0vertexDiff"), ft0.posZ() - collision.posZ());
378+
histos.fill(HIST("hT0vertexDiff_vs_nContrib"),ft0.posZ() - collision.posZ(), nContrib);
379+
histos.fill(HIST("hT0vertexDiff_vs_MultT0AC"), ft0.posZ() - collision.posZ(), multFT0M);
367380

368381
if (nContrib > 20) {
369382

0 commit comments

Comments
 (0)