@@ -811,48 +811,48 @@ void MatchITSTPCQC::finalize()
811811 // first we use denominators and nominators to set the TEfficiency; later they are scaled
812812
813813 // some checks
814- for (int i = 0 ; i < matchType::SIZE ; ++i ) {
815- for (int i = 0 ; i < mPtDen [i ]->GetNbinsX (); ++i) {
816- if (mPtDen [i ]->GetBinContent (i + 1 ) < mPtNum [i ]->GetBinContent (i + 1 )) {
817- LOG (error) << title[i ] << " : bin " << i + 1 << " in [" << mPtNum [i ]->GetBinLowEdge (i + 1 ) << " , " << mPtNum [i ]->GetBinLowEdge (i + 1 ) + mPtNum [i ]->GetBinWidth (i + 1 ) << " ]: mPtDen[i] = " << mPtDen [i ]->GetBinContent (i + 1 ) << " , mPtNum[i] = " << mPtNum [i ]->GetBinContent (i + 1 );
814+ for (int ti = 0 ; ti < matchType::SIZE ; ++ti ) {
815+ for (int i = 0 ; i < mPtDen [ti ]->GetNbinsX (); ++i) {
816+ if (mPtDen [ti ]->GetBinContent (i + 1 ) < mPtNum [ti ]->GetBinContent (i + 1 )) {
817+ LOG (error) << title[ti ] << " : bin " << i + 1 << " in [" << mPtNum [ti ]->GetBinLowEdge (i + 1 ) << " , " << mPtNum [ti ]->GetBinLowEdge (i + 1 ) + mPtNum [ti ]->GetBinWidth (i + 1 ) << " ]: mPtDen[i] = " << mPtDen [ti ]->GetBinContent (i + 1 ) << " , mPtNum[i] = " << mPtNum [ti ]->GetBinContent (i + 1 );
818818 }
819819 }
820- for (int i = 0 ; i < mPtDen_noEta0 [i ]->GetNbinsX (); ++i) {
821- if (mPtDen_noEta0 [i ]->GetBinContent (i + 1 ) < mPtNum_noEta0 [i ]->GetBinContent (i + 1 )) {
822- LOG (error) << title[i ] << " : bin " << i + 1 << " in [" << mPtNum_noEta0 [i ]->GetBinLowEdge (i + 1 ) << " , " << mPtNum_noEta0 [i ]->GetBinLowEdge (i + 1 ) + mPtNum_noEta0 [i ]->GetBinWidth (i + 1 ) << " ]: mPtDen_noEta0[i] = " << mPtDen_noEta0 [i ]->GetBinContent (i + 1 ) << " , mPtNum_noEta0[i] = " << mPtNum_noEta0 [i ]->GetBinContent (i + 1 );
820+ for (int i = 0 ; i < mPtDen_noEta0 [ti ]->GetNbinsX (); ++i) {
821+ if (mPtDen_noEta0 [ti ]->GetBinContent (i + 1 ) < mPtNum_noEta0 [ti ]->GetBinContent (i + 1 )) {
822+ LOG (error) << title[ti ] << " : bin " << i + 1 << " in [" << mPtNum_noEta0 [ti ]->GetBinLowEdge (i + 1 ) << " , " << mPtNum_noEta0 [ti ]->GetBinLowEdge (i + 1 ) + mPtNum_noEta0 [ti ]->GetBinWidth (i + 1 ) << " ]: mPtDen_noEta0[i] = " << mPtDen_noEta0 [ti ]->GetBinContent (i + 1 ) << " , mPtNum_noEta0[i] = " << mPtNum_noEta0 [ti ]->GetBinContent (i + 1 );
823823 }
824824 }
825- for (int i = 0 ; i < mPhiDen [i ]->GetNbinsX (); ++i) {
826- if (mPhiDen [i ]->GetBinContent (i + 1 ) < mPhiNum [i ]->GetBinContent (i + 1 )) {
827- LOG (error) << title[i ] << " : bin " << i + 1 << " in [" << mPhiNum [i ]->GetBinLowEdge (i + 1 ) << " , " << mPhiNum [i ]->GetBinLowEdge (i + 1 ) + mPhiNum [i ]->GetBinWidth (i + 1 ) << " ]: mPhiDen[i] = " << mPhiDen [i ]->GetBinContent (i + 1 ) << " , mPhiNum[i] = " << mPhiNum [i ]->GetBinContent (i + 1 );
825+ for (int i = 0 ; i < mPhiDen [ti ]->GetNbinsX (); ++i) {
826+ if (mPhiDen [ti ]->GetBinContent (i + 1 ) < mPhiNum [ti ]->GetBinContent (i + 1 )) {
827+ LOG (error) << title[ti ] << " : bin " << i + 1 << " in [" << mPhiNum [ti ]->GetBinLowEdge (i + 1 ) << " , " << mPhiNum [ti ]->GetBinLowEdge (i + 1 ) + mPhiNum [ti ]->GetBinWidth (i + 1 ) << " ]: mPhiDen[i] = " << mPhiDen [ti ]->GetBinContent (i + 1 ) << " , mPhiNum[i] = " << mPhiNum [ti ]->GetBinContent (i + 1 );
828828 }
829829 }
830- for (int i = 0 ; i < mEtaDen [i ]->GetNbinsX (); ++i) {
831- if (mEtaDen [i ]->GetBinContent (i + 1 ) < mEtaNum [i ]->GetBinContent (i + 1 )) {
832- LOG (error) << title[i ] << " : bin " << i + 1 << " in [" << mEtaNum [i ]->GetBinLowEdge (i + 1 ) << " , " << mEtaNum [i ]->GetBinLowEdge (i + 1 ) + mEtaNum [i ]->GetBinWidth (i + 1 ) << " ]: mEtaDen[i] = " << mEtaDen [i ]->GetBinContent (i + 1 ) << " , mEtaNum[i] = " << mEtaNum [i ]->GetBinContent (i + 1 );
830+ for (int i = 0 ; i < mEtaDen [ti ]->GetNbinsX (); ++i) {
831+ if (mEtaDen [ti ]->GetBinContent (i + 1 ) < mEtaNum [ti ]->GetBinContent (i + 1 )) {
832+ LOG (error) << title[ti ] << " : bin " << i + 1 << " in [" << mEtaNum [ti ]->GetBinLowEdge (i + 1 ) << " , " << mEtaNum [ti ]->GetBinLowEdge (i + 1 ) + mEtaNum [ti ]->GetBinWidth (i + 1 ) << " ]: mEtaDen[i] = " << mEtaDen [ti ]->GetBinContent (i + 1 ) << " , mEtaNum[i] = " << mEtaNum [ti ]->GetBinContent (i + 1 );
833833 }
834834 }
835835
836836 // filling the efficiency
837- setEfficiency (mFractionITSTPCmatch [i ], mPtNum [i ], mPtDen [i ]);
838- setEfficiency (mFractionITSTPCmatch_noEta0 [i ], mPtNum_noEta0 [i ], mPtDen_noEta0 [i ]);
839- setEfficiency (mFractionITSTPCmatchPhi [i ], mPhiNum [i ], mPhiDen [i ]);
840- setEfficiency (mFractionITSTPCmatchEta [i ], mEtaNum [i ], mEtaDen [i ]);
841- setEfficiency (mFractionITSTPCmatchPhiVsPt [i ], mPhiVsPtNum [i ], mPhiVsPtDen [i ], true );
842- setEfficiency (mFractionITSTPCmatchEtaVsPt [i ], mEtaVsPtNum [i ], mEtaVsPtDen [i ], true );
843- setEfficiency (mFractionITSTPCmatch1OverPt [i ], m1OverPtNum[i ], m1OverPtDen[i ]);
837+ setEfficiency (mFractionITSTPCmatch [ti ], mPtNum [ti ], mPtDen [ti ]);
838+ setEfficiency (mFractionITSTPCmatch_noEta0 [ti ], mPtNum_noEta0 [ti ], mPtDen_noEta0 [ti ]);
839+ setEfficiency (mFractionITSTPCmatchPhi [ti ], mPhiNum [ti ], mPhiDen [ti ]);
840+ setEfficiency (mFractionITSTPCmatchEta [ti ], mEtaNum [ti ], mEtaDen [ti ]);
841+ setEfficiency (mFractionITSTPCmatchPhiVsPt [ti ], mPhiVsPtNum [ti ], mPhiVsPtDen [ti ], true );
842+ setEfficiency (mFractionITSTPCmatchEtaVsPt [ti ], mEtaVsPtNum [ti ], mEtaVsPtDen [ti ], true );
843+ setEfficiency (mFractionITSTPCmatch1OverPt [ti ], m1OverPtNum[ti ], m1OverPtDen[ti ]);
844844 if (mUseTrkPID ) { // Vs Tracking PID hypothesis
845845 for (int j = 0 ; j < o2::track::PID ::NIDs; ++j) {
846- setEfficiency (mFractionITSTPCmatchPtVsTrkPID [i ][j], mPtNumVsTrkPID [i ][j], mPtDenVsTrkPID [i ][j]);
847- setEfficiency (mFractionITSTPCmatchPhiVsTrkPID [i ][j], mPhiNumVsTrkPID [i ][j], mPhiDenVsTrkPID [i ][j]);
848- setEfficiency (mFractionITSTPCmatchEtaVsTrkPID [i ][j], mEtaNumVsTrkPID [i ][j], mEtaDenVsTrkPID [i ][j]);
846+ setEfficiency (mFractionITSTPCmatchPtVsTrkPID [ti ][j], mPtNumVsTrkPID [ti ][j], mPtDenVsTrkPID [ti ][j]);
847+ setEfficiency (mFractionITSTPCmatchPhiVsTrkPID [ti ][j], mPhiNumVsTrkPID [ti ][j], mPhiDenVsTrkPID [ti ][j]);
848+ setEfficiency (mFractionITSTPCmatchEtaVsTrkPID [ti ][j], mEtaNumVsTrkPID [ti ][j], mEtaDenVsTrkPID [ti ][j]);
849849 }
850850 }
851851 if (mUseMC ) {
852- setEfficiency (mFractionITSTPCmatchPhysPrim [i ], mPtPhysPrimNum [i ], mPtPhysPrimDen [i ]);
853- setEfficiency (mFractionITSTPCmatchPhiPhysPrim [i ], mPhiPhysPrimNum [i ], mPhiPhysPrimDen [i ]);
854- setEfficiency (mFractionITSTPCmatchEtaPhysPrim [i ], mEtaPhysPrimNum [i ], mEtaPhysPrimDen [i ]);
855- setEfficiency (mFractionITSTPCmatchPhysPrim1OverPt [i ], m1OverPtPhysPrimNum[i ], m1OverPtPhysPrimDen[i ]);
852+ setEfficiency (mFractionITSTPCmatchPhysPrim [ti ], mPtPhysPrimNum [ti ], mPtPhysPrimDen [ti ]);
853+ setEfficiency (mFractionITSTPCmatchPhiPhysPrim [ti ], mPhiPhysPrimNum [ti ], mPhiPhysPrimDen [ti ]);
854+ setEfficiency (mFractionITSTPCmatchEtaPhysPrim [ti ], mEtaPhysPrimNum [ti ], mEtaPhysPrimDen [ti ]);
855+ setEfficiency (mFractionITSTPCmatchPhysPrim1OverPt [ti ], m1OverPtPhysPrimNum[ti ], m1OverPtPhysPrimDen[ti ]);
856856 }
857857 }
858858 /*
0 commit comments