Functions
def	vbtfRatioAverage
def	vbtfRatioCheck
def	vbtfXSectionAverage
def	vbtfXSectionCheck
Variables
tuple	absCorrRatioenu = Ratioenusqrt(0.0172+0.0127*2)
tuple	absCorrRatiomunu = Ratiomunusqrt(0.0212+0.0119*2)
tuple	absCorrWenu = Wenusqrt(0.0082+0.0116*2)
tuple	absCorrWminusenu = Wminusenusqrt(0.0152+0.0090*2)
tuple	absCorrWminusmunu = Wminusmunusqrt(0.0192+0.0126*2)
tuple	absCorrWmunu = Wmunusqrt(0.0112+0.0136*2)
tuple	absCorrWplusenu = Wplusenusqrt(0.0092+0.0133*2)
tuple	absCorrWplusmunu = Wplusmunusqrt(0.0132+0.0142*2)
tuple	absCorrWZe = WZesqrt(0.0092+0.0103*2)
tuple	absCorrWZmu = WZmusqrt(0.0112+0.0135*2)
tuple	absCorrZee = Zeesqrt(0.0112+0.0134*2)
tuple	absCorrZmumu = Zmumusqrt(0.0122+0.0158*2)
float	absStatRatioenu = 0.02883
float	absStatRatiomunu = 0.026
float	absStatWenu = 0.1123
float	absStatWminusenu = 0.0690
float	absStatWminusmunu = 0.057
float	absStatWmunu = 0.090
float	absStatWplusenu = 0.0793
float	absStatWplusmunu = 0.069
float	absStatWZe = 0.425
float	absStatWZmu = 0.3675
float	absStatZee = 0.0370
float	absStatZmumu = 0.031
float	absUncRatioenu = 0.076
tuple	absUncRatiomunu = Ratiomunusqrt(relUncEff2+relUncMomRes2+relUncMCStat2+relUncBkg*2)
float	absUncWenu = 0.5031
float	absUncWminusenu = 0.2432
tuple	absUncWminusmunu = Wminusmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
tuple	absUncWmunu = Wmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
float	absUncWplusenu = 0.3446
tuple	absUncWplusmunu = Wplusmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
float	absUncWZe = 0.524
tuple	absUncWZmu = WZmusqrt((absUncWmunu/Wmunu)2 + (absUncZmumu/Zmumu)2 - relSysSubtract2 + relSysAdd*2)
float	absUncZee = 0.0565
tuple	absUncZmumu = Zmumusqrt(relUncEff2+relUncMomRes2+relUncBkg2+relUncMCStat*2)
string	OPTION = "StatisticalPlusSystematicAverage"
	Ratioenu = Wplusenu/Wminusenu
	Ratiomunu = Wplusmunu/Wminusmunu
tuple	relSysAdd = abs(relUncMomResW-relUncMomResZ)
int	relSysLumi = 11
	INPUTS FOLLOW (they are not the final version)
tuple	relSysSubtract = sqrt(relUncEffW2+relUncMomResW2+relUncMomResZ**2)
tuple	relUncBkg = sqrt(2.0e-22+0.2e-22)
tuple	relUncEff = sqrt(relUncFit2+relUncPreTrig2+relUncSysEff**2)
float	relUncEffW = 1.3e-2
float	relUncFit = 1.3e-2
float	relUncMCStat = 1.4e-3
float	relUncMomRes = 0.3e-2
float	relUncMomResW = 0.3e-2
float	relUncMomResZ = 0.2e-2
float	relUncPreTrig = 0.5e-2
float	relUncRecoil = 0.4e-2
float	relUncSysEff = 0.5e-2
float	Wenu = 10.1447
float	Wminusenu = 4.1401
float	Wminusmunu = 4.078
float	Wmunu = 9.922
float	Wplusenu = 5.9349
float	Wplusmunu = 5.844
	WZe = Wenu/Zee
	WZmu = Wmunu/Zmumu
float	Zee = 0.9595
float	Zmumu = 0.924

Function Documentation

def Averages_3Invpb::vbtfRatioAverage	(	title,
		ratio1,
		ratio2,
		eratio1,
		eratio2,
		sysUnc1,
		sysUnc2,
		sysCor1,
		sysCor2
	)

Definition at line 269 of file Averages_3Invpb.py.

00270                                                                                                  :
00271       if OPTION== "StatisticalAverage":
00272             V11 = eratio1**2
00273             V22 = eratio2**2
00274             V12 = 0
00275       else:
00276             V11 = eratio1**2+(sysUnc1**2+sysCor1**2)
00277             V22 = eratio2**2+(sysUnc2**2+sysCor2**2)
00278             V12 = sysCor1*sysCor2
00279 
00280       a1 = (V22-V12)/(V11+V22-2*V12)
00281       a2 = (V11-V12)/(V11+V22-2*V12)
00282       average = a1*ratio1 + a2*ratio2
00283       errstat = sqrt(a1**2*eratio1**2+a2**2*eratio2**2)
00284       errunco = sqrt(a1**2*sysUnc1**2 + a2**2*sysUnc2**2)
00285       errtheo = sqrt(a1**2*sysCor1**2 + a2**2*sysCor2**2 + 2*a1*a2*sysCor1*sysCor2)
00286       errsyst = sqrt(errunco**2+errtheo**2)
00287 
00288       print "VBTF average: %s = %.4f +- %.4f (stat.)" % (title, average, errstat) 
00289 
00290       print "\tVBTF systematics (1): +- %.4f (exp) +- %.4f (the)" % (errunco, errtheo) 
00291       print "\tVBTF systematics (2): +- %.4f (exp+the)" % (errsyst)

def Averages_3Invpb::vbtfRatioCheck	(	title,
		ratio,
		absStat,
		sysUnc,
		sysCor
	)

Definition at line 265 of file Averages_3Invpb.py.

00266                                                          :
00267       print "VBTF inputs: %s = %.4f +- %.4f (stat.) +- %.4f (exp.) +- %.4f (the.)" % (title, ratio, absStat, sysUnc, sysCor)

def Averages_3Invpb::vbtfXSectionAverage	(	title,
		xsection1,
		xsection2,
		exsection1,
		exsection2,
		sysUnc1,
		sysUnc2,
		sysCor1,
		sysCor2,
		relSysLumi
	)

Definition at line 240 of file Averages_3Invpb.py.

00241                                                                                                                             :
00242       if OPTION== "StatisticalAverage":
00243             V11 = exsection1**2
00244             V22 = exsection2**2
00245             V12 = 0
00246       else:
00247             V11 = exsection1**2+sysUnc1**2+sysCor1**2
00248             V22 = exsection2**2+sysUnc2**2+sysCor2**2
00249             V12 = sysCor1*sysCor2
00250 
00251       a1 = (V22-V12)/(V11+V22-2*V12)
00252       a2 = (V11-V12)/(V11+V22-2*V12)
00253       average = a1*xsection1 + a2*xsection2
00254       errstat = sqrt(a1**2*exsection1**2+a2**2*exsection2**2)
00255       errunco = sqrt(a1**2*sysUnc1**2 + a2**2*sysUnc2**2)
00256       errtheo = sqrt(a1**2*sysCor1**2 + a2**2*sysCor2**2 + 2*a1*a2*sysCor1*sysCor2)
00257       errsyst = sqrt(errunco**2+errtheo**2)
00258 
00259       print "VBTF average: %s = %.4f +- %.4f (stat.) [nb]" % (title, average, errstat) 
00260 
00261       absSysLumi = average*relSysLumi
00262       print "\tVBTF systematics (1): +- %.4f (exp) +- %.4f (the) +- %.4f (lumi) [nb]" % (errunco, errtheo, absSysLumi) 
00263       print "\tVBTF systematics (2): +- %.4f (exp+the) +- %.4f (lumi) [nb]" % (errsyst, absSysLumi)

def Averages_3Invpb::vbtfXSectionCheck	(	title,
		xsection,
		exsection,
		sysUnc,
		sysCor,
		relSysLumi
	)

Definition at line 235 of file Averages_3Invpb.py.

00236                                                                              :
00237       absSysLumi = xsection*relSysLumi
00238       print "VBTF inputs: %s = %.4f +- %.4f (stat.) +- %.4f (exp.) +- %.4f (the.) +- %.4f (lumi) [nb]" % (title, xsection, exsection, sysUnc, sysCor, absSysLumi)