@@ -110,6 +110,8 @@ std::array<std::shared_ptr<TH1>, nParticles> hPtGeneratedRecoEv;
110110std::array<std::shared_ptr<TH1 >, nParticles> hPtItsPrm;
111111std::array<std::shared_ptr<TH1 >, nParticles> hPtItsTpcPrm;
112112std::array<std::shared_ptr<TH1 >, nParticles> hPtTrkItsTpcPrm;
113+ std::array<std::shared_ptr<TH2 >, nParticles> hDeltaPtVsPtTrkItsTpcPrm;
114+ std::array<std::shared_ptr<TH2 >, nParticles> hPtGenVsPtTrkItsTpcPrm;
113115std::array<std::shared_ptr<TH1 >, nParticles> hPtItsTpcTofPrm;
114116std::array<std::shared_ptr<TH1 >, nParticles> hPtTrkItsTpcTofPrm;
115117std::array<std::shared_ptr<TH1 >, nParticles> hPtGeneratedPrm;
@@ -376,6 +378,8 @@ struct QaEfficiency {
376378 hPtItsPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/its" , PDGs[histogramIndex]), " ITS tracks (primaries) " + tagPt, kTH1D , {axisPt});
377379 hPtItsTpcPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/its_tpc" , PDGs[histogramIndex]), " ITS-TPC tracks (primaries) " + tagPt, kTH1D , {axisPt});
378380 hPtTrkItsTpcPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/trk/its_tpc" , PDGs[histogramIndex]), " ITS-TPC tracks (reco primaries) " + tagPt, kTH1D , {axisPt});
381+ hDeltaPtVsPtTrkItsTpcPrm[histogramIndex] = histos.add <TH2 >(Form (" MC/pdg%i/pt/prm/generated_vs_reco_delta" , PDGs[histogramIndex]), " Abs(Gen - Reco) pT vs Gen pT (primaries) " + tagPt, kTH2D , {axisPt, axisPt});
382+ hPtGenVsPtTrkItsTpcPrm[histogramIndex] = histos.add <TH2 >(Form (" MC/pdg%i/pt/prm/generated_vs_reco" , PDGs[histogramIndex]), " Reco pT vs Gen pT (primaries) " + tagPt, kTH2D , {axisPt, axisPt});
379383 hPtItsTpcTofPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/its_tpc_tof" , PDGs[histogramIndex]), " ITS-TPC-TOF tracks (primaries) " + tagPt, kTH1D , {axisPt});
380384 hPtTrkItsTpcTofPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/trk/its_tpc_tof" , PDGs[histogramIndex]), " ITS-TPC-TOF tracks (reco primaries) " + tagPt, kTH1D , {axisPt});
381385 hPtGeneratedPrm[histogramIndex] = histos.add <TH1 >(Form (" MC/pdg%i/pt/prm/generated" , PDGs[histogramIndex]), " Generated (primaries) " + tagPt, kTH1D , {axisPt});
@@ -1156,6 +1160,8 @@ struct QaEfficiency {
11561160 if (passedITS && passedTPC) {
11571161 hPtItsTpcPrm[histogramIndex]->Fill (mcParticle.pt ());
11581162 hPtTrkItsTpcPrm[histogramIndex]->Fill (track.pt ());
1163+ hDeltaPtVsPtTrkItsTpcPrm[histogramIndex]->Fill (mcParticle.pt (), abs (track.pt () - mcParticle.pt ()));
1164+ hPtGenVsPtTrkItsTpcPrm[histogramIndex]->Fill (mcParticle.pt (), track.pt ());
11591165 hEtaItsTpcPrm[histogramIndex]->Fill (mcParticle.eta ());
11601166 hEtaTrkItsTpcPrm[histogramIndex]->Fill (track.eta ());
11611167 hPhiItsTpcPrm[histogramIndex]->Fill (mcParticle.phi ());
0 commit comments