Averages_3Invpb Namespace Reference

Functions
def	vbtfRatioAverage (title, ratio1, ratio2, eratio1, eratio2, sysUnc1, sysUnc2, sysCor1, sysCor2)
def	vbtfRatioCheck (title, ratio, absStat, sysUnc, sysCor)
def	vbtfXSectionAverage (title, xsection1, xsection2, exsection1, exsection2, sysUnc1, sysUnc2, sysCor1, sysCor2, relSysLumi)
def	vbtfXSectionCheck (title, xsection, exsection, sysUnc, sysCor, relSysLumi)
	INPUTS FOLLOW (they are not the final version) More...

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) 

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

Definition at line 265 of file Averages_3Invpb.py.

Referenced by vbtfRatioAverage().

def vbtfRatioCheck(title, ratio, absStat, sysUnc, sysCor):
      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.

References mathSSE.sqrt().

Referenced by vbtfRatioAverage().

def vbtfXSectionAverage(title, xsection1, xsection2, exsection1, exsection2, sysUnc1, sysUnc2, sysCor1, sysCor2, relSysLumi):
      if OPTION== "StatisticalAverage":
            V11 = exsection1**2
            V22 = exsection2**2
            V12 = 0
      else:
            V11 = exsection1**2+sysUnc1**2+sysCor1**2
            V22 = exsection2**2+sysUnc2**2+sysCor2**2
            V12 = sysCor1*sysCor2

      a1 = (V22-V12)/(V11+V22-2*V12)
      a2 = (V11-V12)/(V11+V22-2*V12)
      average = a1*xsection1 + a2*xsection2
      errstat = sqrt(a1**2*exsection1**2+a2**2*exsection2**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.) [nb]" % (title, average, errstat) 

      absSysLumi = average*relSysLumi
      print "\tVBTF systematics (1): +- %.4f (exp) +- %.4f (the) +- %.4f (lumi) [nb]" % (errunco, errtheo, absSysLumi) 
      print "\tVBTF systematics (2): +- %.4f (exp+the) +- %.4f (lumi) [nb]" % (errsyst, absSysLumi) 

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

INPUTS FOLLOW (they are not the final version)

Utility functions

Definition at line 235 of file Averages_3Invpb.py.

Referenced by vbtfRatioAverage().

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