@@ -50,6 +50,7 @@ struct wMuonFwdEfficiency {
5050 const AxisSpec axisPt{20 , 0.0 , +80.0 , " p_{T} (GeV/c)" };
5151 const AxisSpec axisDeltaPt{40 , 0.0 , +20.0 , " |p_{T}^{true} - p_{T}^{reco}|(GeV/c)" };
5252 const AxisSpec axisChi2{50 , 0.0 , +10.0 , " #chi^{2}" };
53+ const AxisSpec axisChi2Global{50 , 0.0 , +100.0 , " #chi^{2}" };
5354 const AxisSpec axisDCA{100 , 0.0 , +1000.0 , " pDCA" };
5455
5556 // create histograms
@@ -61,7 +62,9 @@ struct wMuonFwdEfficiency {
6162 histos.add (" PtTruthHist" , " PtTruthHist" , kTH1F , {axisPt});
6263 histos.add (" PtResolution" , " PtResolution" , kTH1F , {axisDeltaPt});
6364 histos.add (" chi2" , " chi2" , kTH1F , {axisChi2});
65+ histos.add (" chi2_global_track" , " chi2_global_track" , kTH1F , {axisChi2Global});
6466 histos.add (" chi2MatchMCHMID" , " chi2MatchMCHMID" , kTH1F , {axisChi2});
67+ histos.add (" chi2MatchMCHMFT" , " chi2MatchMCHMFT" , kTH1F , {axisChi2});
6568 histos.add (" pDCA" , " pDCA" , kTH1F , {axisDCA});
6669 }
6770
@@ -141,6 +144,12 @@ struct wMuonFwdEfficiency {
141144 histos.fill (HIST (" chi2MatchMCHMID" ), muChi2MatchMCHMID);
142145 histos.fill (HIST (" pDCA" ), muDca);
143146 }
147+
148+ // global muon tracks
149+ if (muTrackType == 0 && muChi2MatchMCHMFT > 0 ) {
150+ histos.fill (HIST (" chi2_global_track" ), muChi2);
151+ histos.fill (HIST (" chi2MatchMCHMFT" ), muMatchScoreMCHMFT);
152+ }
144153 }
145154 }
146155 }
@@ -174,7 +183,7 @@ struct wMuonFwdEfficiency {
174183 trackInfo.trackID , trackInfo.trackType , trackInfo.chi2 , trackInfo.chi2MatchMCHMID , trackInfo.chi2MatchMCHMFT , trackInfo.matchScoreMCHMFT , trackInfo.matchMFTTrackId );
175184 }
176185 }
177- LOGF (info, " Chosen track has index %d with track ID %ld" , chosenIndex, recoTracks[chosenIndex].trackID );
186+ // LOGF(info, "Chosen track has index %d with track ID %ld", chosenIndex, recoTracks[chosenIndex].trackID);
178187 }
179188 }
180189 }
0 commit comments