Ejemplo n.º 1
0
228.0 	0.00E+00 	4.59E-05 	1.35E-04 	7.10E-01 	2.89E-01 	2.62E+00
230.0 	0.00E+00 	4.45E-05 	1.31E-04 	7.09E-01 	2.90E-01 	2.73E+00
232.0 	0.00E+00 	4.32E-05 	1.28E-04 	7.09E-01 	2.91E-01 	2.84E+00
234.0 	0.00E+00 	4.19E-05 	1.25E-04 	7.08E-01 	2.92E-01 	2.95E+00
236.0 	0.00E+00 	4.07E-05 	1.22E-04 	7.07E-01 	2.93E-01 	3.06E+00
238.0 	0.00E+00 	3.95E-05 	1.20E-04 	7.06E-01 	2.94E-01 	3.18E+00
240.0 	0.00E+00 	3.84E-05 	1.17E-04 	7.06E-01 	2.94E-01 	3.30E+00
242.0 	0.00E+00 	3.73E-05 	1.14E-04 	7.05E-01 	2.95E-01 	3.42E+00
244.0 	0.00E+00 	3.63E-05 	1.12E-04 	7.04E-01 	2.96E-01 	3.54E+00
246.0 	0.00E+00 	3.53E-05 	1.09E-04 	7.04E-01 	2.96E-01 	3.67E+00
248.0 	0.00E+00 	3.44E-05 	1.07E-04 	7.03E-01 	2.97E-01 	3.80E+00
250.0 	0.00E+00 	3.34E-05 	1.04E-04 	7.02E-01 	2.97E-01 	3.93E+00
"""

fermioPhobicBR = _convertFromText(
    _fermioPhobicBRData,
    ['mass', 'gluonGluon', 'gammaGamma', "Zgamma", "WW", "ZZ", "totalWidth"])


def makeBRPlot():
    import pylab
    pylab.ioff()

    # don't plot gluon gluon (which is zero everywhere) in a log plot
    pylab.semilogy(fermioPhobicBR['mass'],
                   fermioPhobicBR['gammaGamma'],
                   label='gamma gamma',
                   linewidth=3)
    pylab.semilogy(fermioPhobicBR['mass'],
                   fermioPhobicBR['Zgamma'],
                   label='Z gamma',
Ejemplo n.º 2
0
820.0 	2.00E-04 	2.69E-07 	1.22E-06 	5.97E-01 	2.97E-01 	3.30E+02
830.0 	1.94E-04 	2.82E-07 	1.16E-06 	5.99E-01 	2.98E-01 	3.43E+02
840.0 	1.88E-04 	2.93E-07 	1.12E-06 	6.01E-01 	2.99E-01 	3.57E+02
850.0 	1.82E-04 	3.04E-07 	1.07E-06 	6.02E-01 	3.00E-01 	3.71E+02
860.0 	1.76E-04 	3.15E-07 	1.02E-06 	6.04E-01 	3.01E-01 	3.86E+02
870.0 	1.71E-04 	3.27E-07 	9.83E-07 	6.05E-01 	3.02E-01 	4.01E+02
880.0 	1.65E-04 	3.37E-07 	9.44E-07 	6.06E-01 	3.03E-01 	4.16E+02
890.0 	1.60E-04 	3.48E-07 	9.07E-07 	6.08E-01 	3.03E-01 	4.32E+02
900.0 	1.55E-04 	3.58E-07 	8.72E-07 	6.09E-01 	3.04E-01 	4.49E+02
910.0 	1.50E-04 	3.67E-07 	8.38E-07 	6.10E-01 	3.05E-01 	4.66E+02
920.0 	1.46E-04 	3.77E-07 	8.07E-07 	6.12E-01 	3.06E-01 	4.84E+02
930.0 	1.41E-04 	3.86E-07 	7.77E-07 	6.13E-01 	3.06E-01 	5.02E+02
940.0 	1.37E-04 	3.93E-07 	7.49E-07 	6.14E-01 	3.07E-01 	5.21E+02
950.0 	1.33E-04 	4.00E-07 	7.22E-07 	6.16E-01 	3.08E-01 	5.40E+02
960.0 	1.29E-04 	4.08E-07 	6.97E-07 	6.17E-01 	3.08E-01 	5.60E+02
970.0 	1.25E-04 	4.17E-07 	6.73E-07 	6.18E-01 	3.09E-01 	5.81E+02
980.0 	1.21E-04 	4.26E-07 	6.50E-07 	6.19E-01 	3.10E-01 	6.02E+02
990.0 	1.17E-04 	4.33E-07 	6.29E-07 	6.20E-01 	3.10E-01 	6.24E+02
1000.0 	1.14E-04 	4.39E-07 	6.08E-07 	6.21E-01 	3.11E-01 	6.47E+02
"""

# mH (GeV) 	H->gg 	H->gamma gamma 	H->Zgamma 	H-> WW 	H-> ZZ 	Total width (GeV)
gaugeBosonBR = _convertFromText(
    _gaugeBosonDecayData,
    ['mass', 'gluonGluon', 'gammaGamma', "Zgamma", "WW", "ZZ", "totalWidth"])


#----------------------------------------------------------------------
def getBR(mass, channel):
    return gaugeBosonBR[gaugeBosonBR['mass'] == mass][channel][0]
Ejemplo n.º 3
0
228.0 	0.00E+00 	4.59E-05 	1.35E-04 	7.10E-01 	2.89E-01 	2.62E+00
230.0 	0.00E+00 	4.45E-05 	1.31E-04 	7.09E-01 	2.90E-01 	2.73E+00
232.0 	0.00E+00 	4.32E-05 	1.28E-04 	7.09E-01 	2.91E-01 	2.84E+00
234.0 	0.00E+00 	4.19E-05 	1.25E-04 	7.08E-01 	2.92E-01 	2.95E+00
236.0 	0.00E+00 	4.07E-05 	1.22E-04 	7.07E-01 	2.93E-01 	3.06E+00
238.0 	0.00E+00 	3.95E-05 	1.20E-04 	7.06E-01 	2.94E-01 	3.18E+00
240.0 	0.00E+00 	3.84E-05 	1.17E-04 	7.06E-01 	2.94E-01 	3.30E+00
242.0 	0.00E+00 	3.73E-05 	1.14E-04 	7.05E-01 	2.95E-01 	3.42E+00
244.0 	0.00E+00 	3.63E-05 	1.12E-04 	7.04E-01 	2.96E-01 	3.54E+00
246.0 	0.00E+00 	3.53E-05 	1.09E-04 	7.04E-01 	2.96E-01 	3.67E+00
248.0 	0.00E+00 	3.44E-05 	1.07E-04 	7.03E-01 	2.97E-01 	3.80E+00
250.0 	0.00E+00 	3.34E-05 	1.04E-04 	7.02E-01 	2.97E-01 	3.93E+00
"""


fermioPhobicBR = _convertFromText(_fermioPhobicBRData,
                            [ 'mass', 'gluonGluon', 'gammaGamma', "Zgamma","WW", "ZZ", "totalWidth" ])


def makeBRPlot():
    import pylab
    pylab.ioff()

    # don't plot gluon gluon (which is zero everywhere) in a log plot
    pylab.semilogy(fermioPhobicBR['mass'], fermioPhobicBR['gammaGamma'], label='gamma gamma', linewidth = 3)
    pylab.semilogy(fermioPhobicBR['mass'], fermioPhobicBR['Zgamma'], label='Z gamma', linewidth = 3)
    pylab.semilogy(fermioPhobicBR['mass'], fermioPhobicBR['WW'], label='WW', linewidth = 3)
    pylab.semilogy(fermioPhobicBR['mass'], fermioPhobicBR['ZZ'], label='ZZ', linewidth = 3)

    # pylab.semilogy(whXsect['mass'],  whXsect['xsect'], label='WH')
    # pylab.semilogy(zhXsect['mass'],  zhXsect['xsect'], label='ZH')
    # pylab.semilogy(tthXsect['mass'], tthXsect['xsect'], label='ttH')