@@ -811,48 +811,48 @@ void MatchITSTPCQC::finalize()
811
811
// first we use denominators and nominators to set the TEfficiency; later they are scaled
812
812
813
813
// some checks
814
- for (int i = 0 ; i < matchType::SIZE; ++i ) {
815
- for (int i = 0 ; i < mPtDen [i ]->GetNbinsX (); ++i) {
814
+ for (int ti = 0 ; ti < matchType::SIZE; ++ti ) {
815
+ for (int i = 0 ; i < mPtDen [ti ]->GetNbinsX (); ++i) {
816
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 );
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 );
818
818
}
819
819
}
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 );
823
823
}
824
824
}
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 );
828
828
}
829
829
}
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 );
833
833
}
834
834
}
835
835
836
836
// 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 ]);
844
844
if (mUseTrkPID ) { // Vs Tracking PID hypothesis
845
845
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]);
849
849
}
850
850
}
851
851
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 ]);
856
856
}
857
857
}
858
858
/*
0 commit comments