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 270 of file Averages_3Invpb.py.

References edm.print(), 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 266 of file Averages_3Invpb.py.

References edm.print().

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 241 of file Averages_3Invpb.py.

References edm.print(), and 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 236 of file Averages_3Invpb.py.

References edm.print().

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))