Averages_3Invpb Namespace Reference

Functions
def	vbtfRatioAverage
def	vbtfRatioCheck
def	vbtfXSectionAverage
def	vbtfXSectionCheck
	Utility functions. More...

Variables
tuple	absCorrRatioenu = Ratioenu*sqrt(0.017**2+0.0127**2)
tuple	absCorrRatiomunu = Ratiomunu*sqrt(0.021**2+0.0119**2)
tuple	absCorrWenu = Wenu*sqrt(0.008**2+0.0116**2)
tuple	absCorrWminusenu = Wminusenu*sqrt(0.015**2+0.0090**2)
tuple	absCorrWminusmunu = Wminusmunu*sqrt(0.019**2+0.0126**2)
tuple	absCorrWmunu = Wmunu*sqrt(0.011**2+0.0136**2)
tuple	absCorrWplusenu = Wplusenu*sqrt(0.009**2+0.0133**2)
tuple	absCorrWplusmunu = Wplusmunu*sqrt(0.013**2+0.0142**2)
tuple	absCorrWZe = WZe*sqrt(0.009**2+0.0103**2)
tuple	absCorrWZmu = WZmu*sqrt(0.011**2+0.0135**2)
tuple	absCorrZee = Zee*sqrt(0.011**2+0.0134**2)
tuple	absCorrZmumu = Zmumu*sqrt(0.012**2+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 = Ratiomunu*sqrt(relUncEff**2+relUncMomRes**2+relUncMCStat**2+relUncBkg**2)
float	absUncWenu = 0.5031
float	absUncWminusenu = 0.2432
tuple	absUncWminusmunu = Wminusmunu*sqrt(relUncEff**2+relUncMomRes**2+relUncRecoil**2+relUncMCStat**2+relUncBkg**2)
tuple	absUncWmunu = Wmunu*sqrt(relUncEff**2+relUncMomRes**2+relUncRecoil**2+relUncMCStat**2+relUncBkg**2)
float	absUncWplusenu = 0.3446
tuple	absUncWplusmunu = Wplusmunu*sqrt(relUncEff**2+relUncMomRes**2+relUncRecoil**2+relUncMCStat**2+relUncBkg**2)
float	absUncWZe = 0.524
tuple	absUncWZmu = WZmu*sqrt((absUncWmunu/Wmunu)**2 + (absUncZmumu/Zmumu)**2 - relSysSubtract**2 + relSysAdd**2)
float	absUncZee = 0.0565
tuple	absUncZmumu = Zmumu*sqrt(relUncEff**2+relUncMomRes**2+relUncBkg**2+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) More...
tuple	relSysSubtract = sqrt(relUncEffW**2+relUncMomResW**2+relUncMomResZ**2)
tuple	relUncBkg = sqrt(2.0e-2**2+0.2e-2**2)
tuple	relUncEff = sqrt(relUncFit**2+relUncPreTrig**2+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.

References mathSSE.sqrt(), vbtfRatioCheck(), vbtfXSectionAverage(), and vbtfXSectionCheck().

def vbtfRatioAverage(title, ratio1, ratio2, eratio1, eratio2, sysUnc1, sysUnc2, sysCor1, sysCor2):
      if OPTION== "StatisticalAverage":
            V11 = eratio1**2
            V22 = eratio2**2
            V12 = 0
      else:
            V11 = eratio1**2+(sysUnc1**2+sysCor1**2)
            V22 = eratio2**2+(sysUnc2**2+sysCor2**2)
            V12 = sysCor1*sysCor2

      a1 = (V22-V12)/(V11+V22-2*V12)
      a2 = (V11-V12)/(V11+V22-2*V12)
      average = a1*ratio1 + a2*ratio2
      errstat = sqrt(a1**2*eratio1**2+a2**2*eratio2**2)
      errunco = sqrt(a1**2*sysUnc1**2 + a2**2*sysUnc2**2)
      errtheo = sqrt(a1**2*sysCor1**2 + a2**2*sysCor2**2 + 2*a1*a2*sysCor1*sysCor2)
      errsyst = sqrt(errunco**2+errtheo**2)

      print "VBTF average: %s = %.4f +- %.4f (stat.)" % (title, average, errstat) 

      print "\tVBTF systematics (1): +- %.4f (exp) +- %.4f (the)" % (errunco, errtheo) 
      print "\tVBTF systematics (2): +- %.4f (exp+the)" % (errsyst)