Averages_36Invpb 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_36Invpb.vbtfRatioAverage	(		title,
			ratio1,
			ratio2,
			eratio1,
			eratio2,
			sysUnc1,
			sysUnc2,
			sysCor1,
			sysCor2
	)

Definition at line 271 of file Averages_36Invpb.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_36Invpb.vbtfRatioCheck	(		title,
			ratio,
			absStat,
			sysUnc,
			sysCor
	)

Definition at line 267 of file Averages_36Invpb.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_36Invpb.vbtfXSectionAverage	(		title,
			xsection1,
			xsection2,
			exsection1,
			exsection2,
			sysUnc1,
			sysUnc2,
			sysCor1,
			sysCor2,
			relSysLumi
	)

Definition at line 242 of file Averages_36Invpb.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_36Invpb.vbtfXSectionCheck	(		title,
			xsection,
			exsection,
			sysUnc,
			sysCor,
			relSysLumi
	)

INPUTS FOLLOW (they are not the final version)

Utility functions

Definition at line 237 of file Averages_36Invpb.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))