Skip to content

Commit a47126a

Browse files
f3schchiarazampolli
authored andcommitted
GLOQC: Add Cluster/Chi2 Efficiency Plots and DCAr_TPC vs P_T
Signed-off-by: Felix Schlepper <felix.schlepper@cern.ch>
1 parent 83cc85a commit a47126a

3 files changed

Lines changed: 190 additions & 22 deletions

File tree

Detectors/GlobalTracking/include/GlobalTracking/MatchITSTPCQC.h

Lines changed: 35 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -59,7 +59,7 @@ class MatchITSTPCQC
5959
bool init();
6060
void initDataRequest();
6161
void run(o2::framework::ProcessingContext& ctx);
62-
void setDataRequest(std::shared_ptr<o2::globaltracking::DataRequest> dr) { mDataRequest = dr; }
62+
void setDataRequest(const std::shared_ptr<o2::globaltracking::DataRequest>& dr) { mDataRequest = dr; }
6363
void finalize();
6464
void reset();
6565

@@ -87,6 +87,14 @@ class MatchITSTPCQC
8787
TH2F* getHistoEtaVsPtDen(matchType m) const { return mEtaVsPtDen[m]; }
8888
TEfficiency* getFractionITSTPCmatchEtaVsPt(matchType m) const { return mFractionITSTPCmatchEtaVsPt[m]; }
8989

90+
TH2F* getHistoClsVsPtNum(matchType m) const { return mClsVsPtNum[m]; }
91+
TH2F* getHistoClsVsPtDen(matchType m) const { return mClsVsPtDen[m]; }
92+
TEfficiency* getFractionITSTPCmatchClsVsPt(matchType m) const { return mFractionITSTPCmatchClsVsPt[m]; }
93+
94+
TH2F* getHistoChi2VsPtNum(matchType m) const { return mChi2VsPtNum[m]; }
95+
TH2F* getHistoChi2VsPtDen(matchType m) const { return mChi2VsPtDen[m]; }
96+
TEfficiency* getFractionITSTPCmatchChi2VsPt(matchType m) const { return mFractionITSTPCmatchChi2VsPt[m]; }
97+
9098
TH1F* getHistoPtPhysPrimNum(matchType m) const { return mPtPhysPrimNum[m]; }
9199
TH1F* getHistoPtPhysPrimDen(matchType m) const { return mPtPhysPrimDen[m]; }
92100
TEfficiency* getFractionITSTPCmatchPhysPrim(matchType m) const { return mFractionITSTPCmatchPhysPrim[m]; }
@@ -107,6 +115,9 @@ class MatchITSTPCQC
107115
TH1F* getHistoChi2Refit() const { return mChi2Refit; }
108116
TH2F* getHistoTimeResVsPt() const { return mTimeResVsPt; }
109117
TH1F* getHistoDCAr() const { return mDCAr; }
118+
TH2F* getHistoDCArVsPtNum() const { return mDCArVsPtNum; }
119+
TH2F* getHistoDCArVsPtDen() const { return mDCArVsPtDen; }
120+
TEfficiency* getFractionITSTPCmatchDCArVsPt() const { return mFractionITSTPCmatchDCArVsPt; }
110121

111122
TH1D* getHisto1OverPtNum(matchType m) const { return m1OverPtNum[m]; }
112123
TH1D* getHisto1OverPtDen(matchType m) const { return m1OverPtDen[m]; }
@@ -122,7 +133,7 @@ class MatchITSTPCQC
122133
/// \tparam T type of the publisher
123134
/// \param publisher the publisher e.g. getObjectsManager()
124135
template <typename T>
125-
void publishHistograms(std::shared_ptr<T> publisher)
136+
void publishHistograms(const std::shared_ptr<T>& publisher)
126137
{
127138
for (int i = 0; i < matchType::SIZE; ++i) {
128139
publisher->startPublishing(mPtNum[i]);
@@ -178,6 +189,14 @@ class MatchITSTPCQC
178189
publisher->startPublishing(mEtaVsPtDen[i]);
179190
publisher->startPublishing(mFractionITSTPCmatchEtaVsPt[i]);
180191

192+
publisher->startPublishing(mClsVsPtNum[i]);
193+
publisher->startPublishing(mClsVsPtDen[i]);
194+
publisher->startPublishing(mFractionITSTPCmatchClsVsPt[i]);
195+
196+
publisher->startPublishing(mChi2VsPtNum[i]);
197+
publisher->startPublishing(mChi2VsPtDen[i]);
198+
publisher->startPublishing(mFractionITSTPCmatchChi2VsPt[i]);
199+
181200
publisher->startPublishing(m1OverPtNum[i]);
182201
publisher->startPublishing(m1OverPtDen[i]);
183202
publisher->startPublishing(mFractionITSTPCmatch1OverPt[i]);
@@ -193,6 +212,9 @@ class MatchITSTPCQC
193212
publisher->startPublishing(mResidualPhi);
194213
publisher->startPublishing(mResidualEta);
195214
publisher->startPublishing(mDCAr);
215+
publisher->startPublishing(mDCArVsPtNum);
216+
publisher->startPublishing(mDCArVsPtDen);
217+
publisher->startPublishing(mFractionITSTPCmatchDCArVsPt);
196218
}
197219

198220
void setSources(GID::mask_t src) { mSrc = src; }
@@ -288,6 +310,14 @@ class MatchITSTPCQC
288310
TH2F* mEtaVsPtNum[matchType::SIZE] = {};
289311
TH2F* mEtaVsPtDen[matchType::SIZE] = {};
290312
TEfficiency* mFractionITSTPCmatchEtaVsPt[matchType::SIZE] = {};
313+
// Clusters
314+
TH2F* mClsVsPtNum[matchType::SIZE] = {};
315+
TH2F* mClsVsPtDen[matchType::SIZE] = {};
316+
TEfficiency* mFractionITSTPCmatchClsVsPt[matchType::SIZE] = {};
317+
// Chi2
318+
TH2F* mChi2VsPtNum[matchType::SIZE] = {};
319+
TH2F* mChi2VsPtDen[matchType::SIZE] = {};
320+
TEfficiency* mFractionITSTPCmatchChi2VsPt[matchType::SIZE] = {};
291321
// Eta split per PID hypothesis in tracking
292322
TH1D* mEtaNumVsTrkPID[matchType::SIZE][track::PID::NIDs] = {};
293323
TH1D* mEtaDenVsTrkPID[matchType::SIZE][track::PID::NIDs] = {};
@@ -301,6 +331,9 @@ class MatchITSTPCQC
301331
TH1F* mChi2Refit = nullptr;
302332
TH2F* mTimeResVsPt = nullptr;
303333
TH1F* mDCAr = nullptr;
334+
TH2F* mDCArVsPtNum = nullptr;
335+
TH2F* mDCArVsPtDen = nullptr;
336+
TEfficiency* mFractionITSTPCmatchDCArVsPt = nullptr;
304337
// 1/Pt
305338
TH1D* m1OverPtNum[matchType::SIZE] = {};
306339
TH1D* m1OverPtDen[matchType::SIZE] = {};

0 commit comments

Comments
 (0)