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 ) ]
) :
