Functions
def	vbtfRatioAverage
def	vbtfRatioCheck
def	vbtfXSectionAverage
def	vbtfXSectionCheck
Variables
float	absCorrRatioenu = 0.0306
tuple	absCorrRatiomunu = Ratiomunusqrt(0.0212+0.0129*2)
float	absCorrWenu = 0.144
float	absCorrWminusenu = 0.073
tuple	absCorrWminusmunu = Wminusmunusqrt(0.0192+0.013*2)
tuple	absCorrWmunu = Wmunusqrt(0.0182+0.015*2)
float	absCorrWplusenu = 0.097
tuple	absCorrWplusmunu = Wplusmunusqrt(0.0132+0.014*2)
float	absCorrWZe = 0.1413
tuple	absCorrWZmu = WZmusqrt(0.0112+0.0135*2)
float	absCorrZee = 0.01715
tuple	absCorrZmumu = Zmumusqrt(0.0122+0.016*2)
float	absStatRatioenu = 0.0097
float	absStatRatiomunu = 0.0078
float	absStatWenu = 0.034
float	absStatWminusenu = 0.022
float	absStatWminusmunu = 0.01662
float	absStatWmunu = 0.02690
float	absStatWplusenu = 0.026
float	absStatWplusmunu = 0.02033
float	absStatWZe = 0.1190
tuple	absStatWZmu = WZmusqrt((absStatWmunu/Wmunu)2+(absStatZmumu/Zmumu)*2)
float	absStatZee = 0.0109
float	absStatZmumu = 0.008
float	absUncRatioenu = 0.0453
tuple	absUncRatiomunu = Ratiomunusqrt(relUncEff2+relUncMomRes2+relUncMCStat2+relUncBkg*2)
float	absUncWenu = 0.309
float	absUncWminusenu = 0.130
tuple	absUncWminusmunu = Wminusmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
tuple	absUncWmunu = Wmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
float	absUncWplusenu = 0.187
tuple	absUncWplusmunu = Wplusmunusqrt(relUncEff2+relUncMomRes2+relUncRecoil2+relUncMCStat2+relUncBkg*2)
float	absUncWZe = 0.2329
tuple	absUncWZmu = WZmusqrt((absUncWmunu/Wmunu)2 + (absUncZmumu/Zmumu)2 - relSysSubtract2 + relSysAdd*2)
float	absUncZee = 0.0446
tuple	absUncZmumu = Zmumusqrt(relUncEff2+relUncMomRes2+relUncBkg2+relUncTrigChanges*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)
float	relUncBkg = 1.5e-2
tuple	relUncEff = sqrt(relUncFit2+relUncPreTrig2+relUncSysEff**2)
float	relUncEffW = 1.3e-2
float	relUncFit = 0.4e-2
float	relUncMCStat = 1.4e-3
float	relUncMomRes = 0.1e-2
float	relUncMomResW = 0.3e-2
float	relUncMomResZ = 0.2e-2
float	relUncPreTrig = 0.5e-2
float	relUncRecoil = 0.4e-2
float	relUncSysEff = 0.4e-2
float	relUncTrigChanges = 0.1e-2
float	Wenu = 10.221
float	Wminusenu = 4.196
float	Wminusmunu = 4.09297
float	Wmunu = 10.03115
float	Wplusenu = 6.045
float	Wplusmunu = 5.93821
	WZe = Wenu/Zee
	WZmu = Wmunu/Zmumu
float	Zee = 0.9892
float	Zmumu = 0.961

Function Documentation

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

Definition at line 270 of file Averages_36Invpb.py.

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

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

Definition at line 266 of file Averages_36Invpb.py.

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

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

Definition at line 241 of file Averages_36Invpb.py.

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

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

Definition at line 236 of file Averages_36Invpb.py.

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