Averages_36Invpb Namespace Reference

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

Variables
float	absCorrRatioenu = 0.0306
tuple	absCorrRatiomunu = Ratiomunu*sqrt(0.021**2+0.0129**2)
float	absCorrWenu = 0.144
float	absCorrWminusenu = 0.073
tuple	absCorrWminusmunu = Wminusmunu*sqrt(0.019**2+0.013**2)
tuple	absCorrWmunu = Wmunu*sqrt(0.018**2+0.015**2)
float	absCorrWplusenu = 0.097
tuple	absCorrWplusmunu = Wplusmunu*sqrt(0.013**2+0.014**2)
float	absCorrWZe = 0.1413
tuple	absCorrWZmu = WZmu*sqrt(0.011**2+0.0135**2)
float	absCorrZee = 0.01715
tuple	absCorrZmumu = Zmumu*sqrt(0.012**2+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 = WZmu*sqrt((absStatWmunu/Wmunu)**2+(absStatZmumu/Zmumu)**2)
float	absStatZee = 0.0109
float	absStatZmumu = 0.008
float	absUncRatioenu = 0.0453
tuple	absUncRatiomunu = Ratiomunu*sqrt(relUncEff**2+relUncMomRes**2+relUncMCStat**2+relUncBkg**2)
float	absUncWenu = 0.309
float	absUncWminusenu = 0.130
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.187
tuple	absUncWplusmunu = Wplusmunu*sqrt(relUncEff**2+relUncMomRes**2+relUncRecoil**2+relUncMCStat**2+relUncBkg**2)
float	absUncWZe = 0.2329
tuple	absUncWZmu = WZmu*sqrt((absUncWmunu/Wmunu)**2 + (absUncZmumu/Zmumu)**2 - relSysSubtract**2 + relSysAdd**2)
float	absUncZee = 0.0446
tuple	absUncZmumu = Zmumu*sqrt(relUncEff**2+relUncMomRes**2+relUncBkg**2+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) More...
tuple	relSysSubtract = sqrt(relUncEffW**2+relUncMomResW**2+relUncMomResZ**2)
float	relUncBkg = 1.5e-2
tuple	relUncEff = sqrt(relUncFit**2+relUncPreTrig**2+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.

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)