integ.GetZaxis().SetTitleOffset(1.2) integ.Draw('SURF4') canv.SetTheta(30) canv.SetPhi(55) canv.Update() if nTupleFilePath : # plot 2D efficiency binned function integrals from ROOT import TH2D integrals += [ TH2D( 'cthetaPhiEffHist', 'Angular Efficiency', 10, -1., +1., 10, -pi, +pi ) , TH2D( 'cpsiPhiEffHist', 'Angular Efficiency', 10, -1., +1., 10, -pi, +pi ) , TH2D( 'cpsiCthetaEffHist', 'Angular Efficiency', 10, -1., +1., 10, -1., +1. ) ] data.fillHistogram( integrals[-3], RooArgList( angles[1], angles[2] ) ) data.fillHistogram( integrals[-2], RooArgList( angles[0], angles[2] ) ) data.fillHistogram( integrals[-1], RooArgList( angles[0], angles[1] ) ) integrals[-3].Scale( 100. / 8. / pi ) integrals[-2].Scale( 100. / 8. / pi ) integrals[-1].Scale( 100. / 8. / pi ) canvs += [ TCanvas( 'cthetaPhiDataCanv', 'Angular Efficiency' ) , TCanvas( 'cpsiPhiDataCanv', 'Angular Efficiency' ) , TCanvas( 'cpsiCthetaDataCanv', 'Angular Efficiency' ) ] for hist, canv, xTitle, yTitle, zTitle, zScale\ in zip( integrals[ -3 : ], canvs[ -3 : ], [ xLabels[1], xLabels[0], xLabels[0] ], [ xLabels[2], xLabels[2], xLabels[1] ] , [ yLabels[2][0], yLabels[2][1], yLabels[2][2] ] , [ ( 0.94, 1.17 ), ( 0.94, 1.17 ), ( 0.94, 1.17 ) ] ) :