Skip to content

Commit 7020250

Browse files
andreasmolanderandreasmolander
committed
Add time res vs tot ampl
1 parent cb4fb46 commit 7020250

1 file changed

Lines changed: 26 additions & 9 deletions

File tree

‎DPG/Tasks/FIT/qaFIT.cxx‎

Lines changed: 26 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -31,6 +31,7 @@
3131
#include <TH1F.h>
3232
#include <TH2.h>
3333
#include <TH2F.h>
34+
#include <TH3.h>
3435
#include <TH3F.h>
3536

3637
#include <array>
@@ -216,20 +217,24 @@ struct fitQa {
216217
/* 2D */
217218

218219
// FT0
219-
OutputObj<TH3F> ooFT0TimeVsFT0Vtx{"FT0TimeVsFT0Vtx"}; ///< FT0 collision time vs FT0 vertex
220-
OutputObj<TH3F> ooPVvsFT0Vtx{"PVvsFT0Vtx"}; ///< PV vs FT0 vertex
221-
OutputObj<TH3F> ooFT0TimeResVsNContrib{"FT0TimeResVsNContrib"}; ///< FT0 time resolution vs number of contributors
220+
OutputObj<TH3F> ooFT0TimeVsFT0Vtx{"FT0TimeVsFT0Vtx"}; ///< FT0 collision time vs FT0 vertex
221+
OutputObj<TH3F> ooPVvsFT0Vtx{"PVvsFT0Vtx"}; ///< PV vs FT0 vertex
222+
OutputObj<TH3F> ooFT0TimeResVsNContrib{"FT0TimeResVsNContrib"}; ///< FT0 time resolution vs number of contributors
223+
OutputObj<TH3F> ooFT0TimeResVsFT0TotAmpl{"FT0TimeResVsFT0TotAmpl"}; ///< FT0 time resolution vs FT0 total amplitude
222224

223225
// FV0
224-
OutputObj<TH3F> ooPVFV0FT0CVtxDiffNSVsNContrib{"PVFV0FT0CVtxDiffNSVsNContrib"}; ///< PV - FV0-FT0C vertex vs number of contributors
226+
OutputObj<TH3F> ooPVFV0FT0CVtxDiffNSVsNContrib{"PVFV0FT0CVtxDiffNSvsNContrib"}; ///< PV - FV0-FT0C vertex vs number of contributors
227+
OutputObj<TH3F> ooPVFV0FT0CVtxDiffNSVsFV0TotAmpl{"PVFV0FT0CVtxDiffNSvsFV0TotAmpl"}; ///< PV - FV0-FT0C vertex vs FV0 total amplitude
225228

226229
// FDD
227-
OutputObj<TH3F> ooFDDTimeVsFDDVtx{"FDDTimeVsFDDVtx"}; ///< FDD collision time vs FDD vertex
228-
OutputObj<TH3F> ooPVvsFDDVtx{"PVvsFDDVtx"}; ///< PV vs FDD vertex
229-
OutputObj<TH3F> ooPVFDDVtxDiffNSVsNContrib{"PVFDDVtxDiffNSVsNContrib"}; ///< PV - FDD vertex vs number of contributors
230+
OutputObj<TH3F> ooFDDTimeVsFDDVtx{"FDDTimeVsFDDVtx"}; ///< FDD collision time vs FDD vertex
231+
OutputObj<TH3F> ooPVvsFDDVtx{"PVvsFDDVtx"}; ///< PV vs FDD vertex
232+
OutputObj<TH3F> ooPVFDDVtxDiffNSVsNContrib{"PVFDDVtxDiffNSvsNContrib"}; ///< PV - FDD vertex vs number of contributors
233+
OutputObj<TH3F> ooPVFDDVtxDiffNSVsFDDTotAmpl{"PVFDDVtxDiffNSvsFDDTotAmpl"}; ///< PV - FDD vertex vs FDD total amplitude
230234

231235
// FT0, FV0
232-
OutputObj<TH3F> ooFT0TimeAFV0TimeDiffVsNContrib{"FT0TimeAFV0TimeDiffVsNContrib"}; ///< FT0A average time - FV0 average time vs number of contributors
236+
OutputObj<TH3F> ooFT0TimeAFV0TimeDiffVsNContrib{"FT0TimeAFV0TimeDiffVsNContrib"}; ///< FT0A average time - FV0 average time vs number of contributors
237+
OutputObj<TH3F> ooFT0TimeAFV0TimeDiffVsFV0TotAmpl{"FT0TimeAFV0TimeDiffVsFV0TotAmpl"}; ///< FT0A average time - FV0 average time vs FV0 total amplitude
233238

234239
/* 2D Quantities per X bin */
235240

@@ -435,15 +440,21 @@ struct fitQa {
435440
ooFT0TimeVsFT0Vtx.setObject(new TH3F(ooFT0TimeVsFT0Vtx.label.c_str(), "FT0 time vs FT0 vertex;$(\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FT0A}} \\rangle)/2 \\text{ (cm)}$;$(\\langle t_{\\text{FT0A}} \\rangle + \\langle t_{\\text{FT0C}} \\rangle)/2 \\text{ (ns)}$", nBinsVtx, vtxMin, vtxMax, nBinsT, tMin, tMax, conditions.size(), 0, conditions.size()));
436441
objs2D[&ooFT0TimeVsFT0Vtx] = [&]() { return std::make_pair(objs[&ooFT0Vtx](), objs[&ooFT0Time]()); };
437442

438-
ooPVvsFT0Vtx.setObject(new TH3F(ooPVvsFT0Vtx.label.c_str(), "FT0 vertex vs PV;$(\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FT0A}} \\rangle)/2 \\text{ (cm)}$;Primary vertex z position (cm)", nBinsVtx, vtxMin, vtxMax, nBinsVtx, vtxMin, vtxMax, conditions.size(), 0, conditions.size()));
443+
ooPVvsFT0Vtx.setObject(new TH3F(ooPVvsFT0Vtx.label.c_str(), "PV vs FT0 vertex;$(\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FT0A}} \\rangle)/2 \\text{ (cm)}$;Primary vertex z position (cm)", nBinsVtx, vtxMin, vtxMax, nBinsVtx, vtxMin, vtxMax, conditions.size(), 0, conditions.size()));
439444
objs2D[&ooPVvsFT0Vtx] = [&]() { return std::make_pair(objs[&ooFT0Vtx](), objs[&ooPV]()); };
440445

441446
ooFT0TimeResVsNContrib.setObject(new TH3F(ooFT0TimeResVsNContrib.label.c_str(), "FT0 time resolution vs number of contributors;Number of contributors to primary vertex;$\\text{PV} - (\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FT0A}} \\rangle)/2 \\text{ (ns)}$", nBinsNContrib / 10, nContribMin, nContribMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
442447
objs2D[&ooFT0TimeResVsNContrib] = [&]() { return std::make_pair(objs[&ooNcontrib](), objs[&ooFT0TimeRes]()); };
443448

449+
ooFT0TimeResVsFT0TotAmpl.setObject(new TH3F(ooFT0TimeResVsFT0TotAmpl.label.c_str(), "FT0 time resolution vs FT0 total amplitude;FT0 total amplitude (ADC);$\\text{PV} - (\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FT0A}} \\rangle)/2 \\text{ (ns)}$", nBinsTotAmpl / 10, totAmplMin, totAmplMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
450+
objs2D[&ooFT0TimeResVsFT0TotAmpl] = [&]() { return std::make_pair(objs[&ooFT0TotAmpl](), objs[&ooFT0TimeRes]()); };
451+
444452
ooPVFV0FT0CVtxDiffNSVsNContrib.setObject(new TH3F(ooPVFV0FT0CVtxDiffNSVsNContrib.label.c_str(), "PV - FV0-FT0C vertex vs number of contributors;Number of contributors to primary vertex;$\\text{PV} - (\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FV0}} \\rangle)/2 \\text{ (ns)}$", nBinsNContrib / 10, nContribMin, nContribMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
445453
objs2D[&ooPVFV0FT0CVtxDiffNSVsNContrib] = [&]() { return std::make_pair(objs[&ooNcontrib](), objs[&ooPVFV0FT0CVtxDiffNS]()); };
446454

455+
ooPVFV0FT0CVtxDiffNSVsFV0TotAmpl.setObject(new TH3F(ooPVFV0FT0CVtxDiffNSVsFV0TotAmpl.label.c_str(), "PV - FV0-FT0C vertex vs FV0 total amplitude;FV0 total amplitude (ADC);$\\text{PV} - (\\langle t_{\\text{FT0C}} \\rangle - \\langle t_{\\text{FV0}} \\rangle)/2 \\text{ (ns)}$", nBinsTotAmpl / 10, totAmplMin, totAmplMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
456+
objs2D[&ooPVFV0FT0CVtxDiffNSVsFV0TotAmpl] = [&]() { return std::make_pair(objs[&ooFV0TotAmpl](), objs[&ooPVFV0FT0CVtxDiffNS]()); };
457+
447458
ooFDDTimeVsFDDVtx.setObject(new TH3F(ooFDDTimeVsFDDVtx.label.c_str(), "FDD collision time vs FDD vertex;$(\\langle t_{\\text{FDDC}} \\rangle - \\langle t_{\\text{FDDA}} \\rangle)/2 \\text{ (cm)}$;$(\\langle t_{\\text{FDDA}} \\rangle + \\langle t_{\\text{FDDC}} \\rangle)/2 \\text{ (ns)}$", nBinsVtxFDD, vtxMinFDD, vtxMaxFDD, nBinsTFDD, tMinFDD, tMaxFDD, conditions.size(), 0, conditions.size()));
448459
objs2D[&ooFDDTimeVsFDDVtx] = [&]() { return std::make_pair(objs[&ooFDDVtx](), objs[&ooFDDTime]()); };
449460

@@ -453,9 +464,15 @@ struct fitQa {
453464
ooPVFDDVtxDiffNSVsNContrib.setObject(new TH3F(ooPVFDDVtxDiffNSVsNContrib.label.c_str(), "PV - FDD vertex vs number of contributors;Number of contributors to primary vertex;$\\text{PV} - (\\langle t_{\\text{FDDC}} \\rangle - \\langle t_{\\text{FDDA}} \\rangle)/2 \\text{ (ns)}$", nBinsNContrib / 10, nContribMin, nContribMax, nBinsTResFDD, tResMinFDD, tResMaxFDD, conditions.size(), 0, conditions.size()));
454465
objs2D[&ooPVFDDVtxDiffNSVsNContrib] = [&]() { return std::make_pair(objs[&ooNcontrib](), objs[&ooPVFDDVtxDiffNS]()); };
455466

467+
// ooPVFDDVtxDiffNSVsFDDTotAmpl.setObject(new TH3F(ooPVFDDVtxDiffNSVsFDDTotAmpl.label.c_str(), "PV - FDD vertex vs FDD total amplitude;FDD total amplitude (ADC);$\\text{PV} - (\\langle t_{\\text{FDDC}} \\rangle - \\langle t_{\\text{FDDA}} \\rangle)/2 \\text{ (ns)}$", nBinsTotAmpl / 10, totAmplMin, totAmplMax, nBinsTResFDD, tResMinFDD, tResMaxFDD, conditions.size(), 0, conditions.size()));
468+
// objs2D[&ooPVFDDVtxDiffNSVsFDDTotAmpl] = [&]() { return std::make_pair(objs[&ooFDDTotAmpl](), objs[&ooPVFDDVtxDiffNS]()); };
469+
456470
ooFT0TimeAFV0TimeDiffVsNContrib.setObject(new TH3F(ooFT0TimeAFV0TimeDiffVsNContrib.label.c_str(), "FT0A time - FV0 time vs number of contributors;Number of contributors to primary vertex;$\\langle t_{\\text{FT0A}} \\rangle - \\langle t_{\\text{FV0}} \\rangle \\text{ (ns)}$", nBinsNContrib / 10, nContribMin, nContribMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
457471
objs2D[&ooFT0TimeAFV0TimeDiffVsNContrib] = [&]() { return std::make_pair(objs[&ooNcontrib](), objs[&ooFT0TimeAFV0TimeDiff]()); };
458472

473+
ooFT0TimeAFV0TimeDiffVsFV0TotAmpl.setObject(new TH3F(ooFT0TimeAFV0TimeDiffVsFV0TotAmpl.label.c_str(), "FT0A time - FV0 time vs FV0 total amplitude;FV0 total amplitude (ADC);$\\langle t_{\\text{FT0A}} \\rangle - \\langle t_{\\text{FV0}} \\rangle \\text{ (ns)}$", nBinsTotAmpl / 10, totAmplMin, totAmplMax, nBinsTRes, tResMin, tResMax, conditions.size(), 0, conditions.size()));
474+
objs2D[&ooFT0TimeAFV0TimeDiffVsFV0TotAmpl] = [&]() { return std::make_pair(objs[&ooFV0TotAmpl](), objs[&ooFT0TimeAFV0TimeDiff]()); };
475+
459476
// 2D quantities per X bin
460477

461478
ooFT0AmplPerCh.setObject(new TH3F(ooFT0AmplPerCh.label.c_str(), "FT0 channel amplitudes;Channel ID;Amplitude (ADC)", nChFT0, 0, static_cast<float>(nChFT0), nADC, 0.f, static_cast<float>(nADC), conditions.size(), 0, conditions.size()));

0 commit comments

Comments
 (0)