(************** Content-type: application/mathematica **************

                    Mathematica-Compatible Notebook

This notebook can be used with any Mathematica-compatible
application, such as Mathematica, MathReader or Publicon. The data
for the notebook starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do
one of the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the
  application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing
the word CacheID, otherwise Mathematica-compatible applications may
try to use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: info@wolfram.com
  phone: +1-217-398-0700 (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
*******************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[     81913,       2234]*)
(*NotebookOutlinePosition[     82552,       2256]*)
(*  CellTagsIndexPosition[     82508,       2252]*)
(*WindowFrame->Normal*)



Notebook[{

Cell[CellGroupData[{
Cell["Phenomenology of pp\[Rule]H+X at NLO", "Subtitle"],

Cell[CellGroupData[{

Cell["Introduction", "Section"],

Cell[TextData[{
  StyleBox["In this notebook we calculate the inclusive cross section for \
Higgs production at hadron colliders, at NLO in the strong coupling. We use \
the analytic results obtained in a previous notebook corresponding to the \
original calculation by Sally Dawson (Nuclear Physics B (1991) 283). To get \
useful numbers we use a modern set of PDF, i.e. the MRST as implemented in ",
    FontSize->16],
  StyleBox["Mathematica",
    FontSize->16,
    FontSlant->"Italic"],
  StyleBox[" by J.Andersen (many thanks!).  A description of the calculation, \
including the formulas used here for the numerical results, is given in the \
notes.",
    FontSize->16]
}], "Text"]
}, Closed]],

Cell[CellGroupData[{

Cell["Preliminaries", "Section"],

Cell["\<\
Off[General::spell];
Off[General::spell1];
Clear[\"Global`*\"];\
\>", "Input"],

Cell[CellGroupData[{

Cell["Install Vegas", "Subsection"],

Cell["\<\
I use the  Vegas package of CUBA library of T. Hahn \
(hep-ph/0404043). If you do not want to use it, use NIntegrate instead of \
Vegas\
\>", "Text"],

Cell[CellGroupData[{

Cell[BoxData[
    \(Install["\</Users/fabiomaltoni/Physics/Codes/CUBA/VegasX\>"] 
      Install["\</Users/fabiomaltoni/Physics/Codes/CUBA/SuaveX\>"]\)], "Input"],

Cell[BoxData[
    \(TraditionalForm\`LinkObject[
        "/Users/fabiomaltoni/Physics/Codes/CUBA/SuaveX", 3, 3]\ LinkObject[
        "/Users/fabiomaltoni/Physics/Codes/CUBA/VegasX", 2, 2]\)], "Output"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{

Cell["Call special graphics routines", "Subsection"],

Cell["\<\
<< Graphics`Colors`;
<< Graphics`Graphics`;\
\>", "Input"]
}, Closed]],

Cell[CellGroupData[{

Cell["Install PDF'S: for help read the PDF-HOWTO document", "Subsection"],

Cell[CellGroupData[{

Cell[BoxData[{
    \(SetDirectory["\<~/Physics/Teaching/pp>Higgs\ at\ NLO/Jeppe\>"]\), "\
\[IndentingNewLine]", 
    \(Install["\<pdf.exe\>"]\)}], "Input"],

Cell[BoxData[
    \(TraditionalForm\`"/Users/fabiomaltoni/Physics/Teaching/pp>Higgs at \
NLO/Jeppe"\)], "Output"],

Cell[BoxData[
    \(TraditionalForm\`LinkObject["./pdf.exe", 4, 4]\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\

<<\"loadCTEQ5.m\";\
\>", "Input"],

Cell[BoxData[
    \(TraditionalForm\`"Loading Package: PDF"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"PDF's from CTEQ5"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"Eur.Phys.J.C12:375-392,2000 "\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"*** Warning *** Unofficial release. "\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"*** Cross-check with official Fortran code."\)], \
"Print"],

Cell[BoxData[
    \(TraditionalForm\`" Version 1.0: Written by Tamara Trout & Fred Olness, \
October 1, 2000"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`" "\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"In case of problems, contact:"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"Fred Olness: olness@mail.physics.smu.edu"\)], "Print"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Wrapper to call the pdf. Notation is self-explanatory. Just notice \
that all parton distribution codes usually return x f(x). To avoid confusion, \
at expense of a couple more floating point operations, I divide all the \
values by the corresponding x.  pdfcall calculates the parton-parton \
luminoties, gg qg, qq~\
\>", "Subsubsection"],

Cell["\<\
(* MRST *)

pdfMRST[X1_,X2_,q_]:=Module[
{Q,xgDF,xdDF,xdbDF,xuDF,xubDF,xsDF,xcDF,xbDF,pd1,pd2,xgg,xqg,xqq},
Q=q*1.;
xgDF =f[3,X1,Q]/X1;
xdbDF=f[8,X1,Q]/X1;
xdDF =f[2,X1,Q]/X1+xdbDF;
xubDF=f[4,X1,Q]/X1;
xuDF =f[1,X1,Q]/X1+xubDF;
xsDF =f[6,X1,Q]/X1;
xcDF =f[5,X1,Q]/X1;
xbDF =f[7,X1,Q]/X1;

pd1={xgDF,xdDF,xdbDF,xuDF,xubDF,xsDF,xcDF,xbDF};
xgDF =f[3,X2,Q]/X2;
xdbDF=f[8,X2,Q]/X2;
xdDF =f[2,X2,Q]/X2+xdbDF;
xubDF=f[4,X2,Q]/X2;
xuDF =f[1,X2,Q]/X2+xubDF;
xsDF =f[6,X2,Q]/X2;
xcDF =f[5,X2,Q]/X2;
xbDF =f[7,X2,Q]/X2;
pd2={xgDF,xdDF,xdbDF,xuDF,xubDF,xsDF,xcDF,xbDF};


xgg=pd1[[1]]*pd2[[1]];
xqg=pd1[[1]]*(pd2[[2]]+pd2[[3]]+pd2[[4]]+pd2[[5]]+2(pd2[[6]]+pd2[[7]]+pd2[[8]]\
))+
   pd2[[1]]*(pd1[[2]]+pd1[[3]]+pd1[[4]]+pd1[[5]]+2(pd1[[6]]+pd1[[7]]+pd1[[8]])\
);
xqq=pd1[[2]]*pd2[[3]]+pd1[[4]]*pd2[[5]]+pd1[[3]]*pd2[[2]]+pd1[[5]]*pd2[[4]]+
   2(pd1[[6]] pd2[[6]]+pd1[[7]] pd2[[7]]+pd1[[8]] pd2[[8]]);
  
  Return[{xgg,xqg,xqq}]; 
];\
\>", "Input"],

Cell["\<\

(* cteq5 *)\
\>", "Input"],

Cell["\<\
pdfCTEQ[X1_,X2_,q_]:=Module[
{Q,pd1,pd2,xgg,xqg,xqq},
Q=q*1.;

pd1={cteq5pdf[1,0,X1,Q],
     cteq5pdf[1,2,X1,Q],
     cteq5pdf[1,-2,X1,Q],
     cteq5pdf[1,1,X1,Q],
     cteq5pdf[1,-1,X1,Q],
     cteq5pdf[1,3,X1,Q],
     cteq5pdf[1,4,X1,Q],
     cteq5pdf[1,5,X1,Q]};
pd2={cteq5pdf[1,0,X2,Q],
     cteq5pdf[1,2,X2,Q],
     cteq5pdf[1,-2,X2,Q],
     cteq5pdf[1,1,X2,Q],
     cteq5pdf[1,-1,X2,Q],
     cteq5pdf[1,3,X2,Q],
     cteq5pdf[1,4,X2,Q],
     cteq5pdf[1,5,X2,Q]};
     
xgg=pd1[[1]]*pd2[[1]];
xqg=pd1[[1]]*(pd2[[2]]+pd2[[3]]+pd2[[4]]+pd2[[5]]+2(pd2[[6]]+pd2[[7]]+pd2[[8]]\
))+
   pd2[[1]]*(pd1[[2]]+pd1[[3]]+pd1[[4]]+pd1[[5]]+2(pd1[[6]]+pd1[[7]]+pd1[[8]])\
);
xqq=pd1[[2]]*pd2[[3]]+pd1[[4]]*pd2[[5]]+pd1[[3]]*pd2[[2]]+pd1[[5]]*pd2[[4]]+
   2(pd1[[6]] pd2[[6]]+pd1[[7]] pd2[[7]]+pd1[[8]] pd2[[8]]);
  
  Return[{xgg,xqg,xqq}]; 
];

pdfCTEQLO[X1_,X2_,q_]:=Module[
{Q,xgDF,xdDF,xdbDF,xuDF,xubDF,xsDF,xcDF,xbDF,pd1,pd2,xgg,xqg,xqq},
Q=q*1.;

pd1={cteq5pdf[3,0,X1,Q],
     cteq5pdf[3,2,X1,Q],
     cteq5pdf[3,-2,X1,Q],
     cteq5pdf[3,1,X1,Q],
     cteq5pdf[3,-1,X1,Q],
     cteq5pdf[3,3,X1,Q],
     cteq5pdf[3,4,X1,Q],
     cteq5pdf[3,5,X1,Q]};
pd2={cteq5pdf[3,0,X2,Q],
     cteq5pdf[3,2,X2,Q],
     cteq5pdf[3,-2,X2,Q],
     cteq5pdf[3,1,X2,Q],
     cteq5pdf[3,-1,X2,Q],
     cteq5pdf[3,3,X2,Q],
     cteq5pdf[3,4,X2,Q],
     cteq5pdf[3,5,X2,Q]};
     
xgg=pd1[[1]]*pd2[[1]];
xqg=pd1[[1]]*(pd2[[2]]+pd2[[3]]+pd2[[4]]+pd2[[5]]+2(pd2[[6]]+pd2[[7]]+pd2[[8]]\
))+
   pd2[[1]]*(pd1[[2]]+pd1[[3]]+pd1[[4]]+pd1[[5]]+2(pd1[[6]]+pd1[[7]]+pd1[[8]])\
);
xqq=pd1[[2]]*pd2[[3]]+pd1[[4]]*pd2[[5]]+pd1[[3]]*pd2[[2]]+pd1[[5]]*pd2[[4]]+
   2(pd1[[6]] pd2[[6]]+pd1[[7]] pd2[[7]]+pd1[[8]] pd2[[8]]);
  
  Return[{xgg,xqg,xqq}]; 
];\
\>", "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Decide the pdf family to be used ", "Subsubsection"],

Cell["pdfcall[x__] = pdfCTEQ[x];", "Input"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Alpha_S: Very basic implementation of Alpha_S. 
Check that the value of Lamda_4 or Lambda_5 is consistent with that of the \
PDF. 
One simple, but indirect way to do it is to compare the value of alphas(MZ) \
with the one quoted by MRST.\
\>", "Subsection"],

Cell[BoxData[{
    \(\(b[nf_] := \(\((33 - 2\ nf)\)/12\)/
          Pi;\)\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \
\ \ \ \), "\[IndentingNewLine]", 
    \(\(bp[
          nf_] := \(\(\((153 - 19\ nf)\)/2\)/
            Pi\)/\((33 - 2\ nf)\);\)\), "\[IndentingNewLine]", 
    \(\(t[nf] := 
        Log[q^2/\[CapitalLambda][nf]^2];\)\), "\[IndentingNewLine]", 
    \(\(\(mb = 4.75;\)\(\[IndentingNewLine]\)
    \)\), "\[IndentingNewLine]", 
    \(\(\[CapitalLambda][5] = 0.146;\)\), "\[IndentingNewLine]", 
    \(\(asLO[q_, nf_] = 1/\((b[nf]\ t[nf])\);\)\), "\n", 
    \(\(\[CapitalLambda][5] = 0.226;\)\), "\[IndentingNewLine]", 
    \(\(asNLO[q_, nf_] = 
        1/\((b[nf]\ t[nf])\) \((1 - 
              bp[nf]/b[nf]\ Log[t[nf]]/t[nf])\);\)\)}], "Input"],

Cell[CellGroupData[{

Cell["\<\
Check the values of alpha_S at the scale MZ\
\>", "Subsubsection"],

Cell[CellGroupData[{

Cell["asNLO[91.118,5]", "Input"],

Cell[BoxData[
    \(TraditionalForm\`0.11799518835498905`\)], "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Color factors ", "Subsection"],

Cell["\<\
CF=4/3;
CA=3;
nf=5;
TF=1/2;
b0=11/6 CA -2 nf TF/3;
\
\>", "Input"]
}, Closed]]
}, Open  ]],

Cell[CellGroupData[{

Cell["LO cross section", "Section"],

Cell[CellGroupData[{

Cell["\<\
Integral of the LO loop. The sign of the Imaginary part is given by \
the usual prescription mt^2-i eps which is equivalent to t+i eps in the \
notation below.
t=mh^2/4/mt^2;\
\>", "Subsubsection"],

Cell[" ", "Input"],

Cell[CellGroupData[{

Cell["\<\
eps=0.00000001;
inte[t_]=3*Integrate[(1-4 x y)/(1-4 t x y),{x,0,1},{y,0,1-x}]//Simplify;
Plot[{Re[inte[t+I eps]],Im[inte[t+I \
eps]]},{t,0,10},PlotStyle\[Rule]{{Blue,Thickness[0.007]},{Red,Thickness[0.007]\
}}];\
\>", "Input"],

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .61803 
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10  scalefont  setfont
% Scaling calculations
0.0238095 0.0952381 0.0147154 0.338579 [
[.21429 .00222 -3 -9 ]
[.21429 .00222 3 0 ]
[.40476 .00222 -3 -9 ]
[.40476 .00222 3 0 ]
[.59524 .00222 -3 -9 ]
[.59524 .00222 3 0 ]
[.78571 .00222 -3 -9 ]
[.78571 .00222 3 0 ]
[.97619 .00222 -6 -9 ]
[.97619 .00222 6 0 ]
[.01131 .09936 -24 -4.5 ]
[.01131 .09936 0 4.5 ]
[.01131 .18401 -18 -4.5 ]
[.01131 .18401 0 4.5 ]
[.01131 .26865 -24 -4.5 ]
[.01131 .26865 0 4.5 ]
[.01131 .35329 -6 -4.5 ]
[.01131 .35329 0 4.5 ]
[.01131 .43794 -24 -4.5 ]
[.01131 .43794 0 4.5 ]
[.01131 .52258 -18 -4.5 ]
[.01131 .52258 0 4.5 ]
[.01131 .60723 -24 -4.5 ]
[.01131 .60723 0 4.5 ]
[ 0 0 0 0 ]
[ 1 .61803 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 g
.25 Mabswid
[ ] 0 setdash
.21429 .01472 m
.21429 .02097 L
s
[(2)] .21429 .00222 0 1 Mshowa
.40476 .01472 m
.40476 .02097 L
s
[(4)] .40476 .00222 0 1 Mshowa
.59524 .01472 m
.59524 .02097 L
s
[(6)] .59524 .00222 0 1 Mshowa
.78571 .01472 m
.78571 .02097 L
s
[(8)] .78571 .00222 0 1 Mshowa
.97619 .01472 m
.97619 .02097 L
s
[(10)] .97619 .00222 0 1 Mshowa
.125 Mabswid
.07143 .01472 m
.07143 .01847 L
s
.11905 .01472 m
.11905 .01847 L
s
.16667 .01472 m
.16667 .01847 L
s
.2619 .01472 m
.2619 .01847 L
s
.30952 .01472 m
.30952 .01847 L
s
.35714 .01472 m
.35714 .01847 L
s
.45238 .01472 m
.45238 .01847 L
s
.5 .01472 m
.5 .01847 L
s
.54762 .01472 m
.54762 .01847 L
s
.64286 .01472 m
.64286 .01847 L
s
.69048 .01472 m
.69048 .01847 L
s
.7381 .01472 m
.7381 .01847 L
s
.83333 .01472 m
.83333 .01847 L
s
.88095 .01472 m
.88095 .01847 L
s
.92857 .01472 m
.92857 .01847 L
s
.25 Mabswid
0 .01472 m
1 .01472 L
s
.02381 .09936 m
.03006 .09936 L
s
[(0.25)] .01131 .09936 1 0 Mshowa
.02381 .18401 m
.03006 .18401 L
s
[(0.5)] .01131 .18401 1 0 Mshowa
.02381 .26865 m
.03006 .26865 L
s
[(0.75)] .01131 .26865 1 0 Mshowa
.02381 .35329 m
.03006 .35329 L
s
[(1)] .01131 .35329 1 0 Mshowa
.02381 .43794 m
.03006 .43794 L
s
[(1.25)] .01131 .43794 1 0 Mshowa
.02381 .52258 m
.03006 .52258 L
s
[(1.5)] .01131 .52258 1 0 Mshowa
.02381 .60723 m
.03006 .60723 L
s
[(1.75)] .01131 .60723 1 0 Mshowa
.125 Mabswid
.02381 .03164 m
.02756 .03164 L
s
.02381 .04857 m
.02756 .04857 L
s
.02381 .0655 m
.02756 .0655 L
s
.02381 .08243 m
.02756 .08243 L
s
.02381 .11629 m
.02756 .11629 L
s
.02381 .13322 m
.02756 .13322 L
s
.02381 .15015 m
.02756 .15015 L
s
.02381 .16708 m
.02756 .16708 L
s
.02381 .20093 m
.02756 .20093 L
s
.02381 .21786 m
.02756 .21786 L
s
.02381 .23479 m
.02756 .23479 L
s
.02381 .25172 m
.02756 .25172 L
s
.02381 .28558 m
.02756 .28558 L
s
.02381 .30251 m
.02756 .30251 L
s
.02381 .31944 m
.02756 .31944 L
s
.02381 .33637 m
.02756 .33637 L
s
.02381 .37022 m
.02756 .37022 L
s
.02381 .38715 m
.02756 .38715 L
s
.02381 .40408 m
.02756 .40408 L
s
.02381 .42101 m
.02756 .42101 L
s
.02381 .45487 m
.02756 .45487 L
s
.02381 .4718 m
.02756 .4718 L
s
.02381 .48873 m
.02756 .48873 L
s
.02381 .50566 m
.02756 .50566 L
s
.02381 .53951 m
.02756 .53951 L
s
.02381 .55644 m
.02756 .55644 L
s
.02381 .57337 m
.02756 .57337 L
s
.02381 .5903 m
.02756 .5903 L
s
.25 Mabswid
.02381 0 m
.02381 .61803 L
s
0 0 m
1 0 L
1 .61803 L
0 .61803 L
closepath
clip
newpath
0 0 1 r
.007 w
.02381 .35329 m
.04262 .3703 L
.06244 .39214 L
.08255 .42072 L
.09396 .44167 L
.09905 .45278 L
.10458 .46661 L
.11004 .48287 L
.11243 .49122 L
.11508 .5017 L
.11756 .51355 L
.11989 .52899 L
.12251 .54656 L
.125 .56035 L
.12723 .57062 L
.12967 .57981 L
.13191 .58664 L
.13395 .59162 L
.13636 .59622 L
.1376 .59807 L
.13894 .59973 L
.14018 .60096 L
.14132 .60184 L
.14242 .60249 L
.14357 .60296 L
.14479 .60325 L
.14548 .60331 L
.14612 .60332 L
.14729 .60319 L
.14853 .60287 L
.14964 .60244 L
.15068 .60191 L
.15305 .60029 L
.1556 .59799 L
.15794 .59542 L
.16321 .58831 L
.18154 .55452 L
.22117 .46854 L
.25928 .39356 L
.29984 .32787 L
.33888 .27688 L
.38038 .23315 L
.42036 .19894 L
.45883 .17172 L
.49975 .14757 L
.53915 .12804 L
.57704 .11202 L
.61738 .09739 L
.6562 .08524 L
.69748 .07404 L
Mistroke
.73724 .06464 L
.77548 .05669 L
.81618 .04921 L
.85536 .04283 L
.897 .03681 L
.93712 .03165 L
.97572 .02719 L
.97619 .02714 L
Mfstroke
1 0 0 r
.02381 .01472 m
.02499 .01472 L
.02605 .01472 L
.02846 .01472 L
.03279 .01472 L
.04262 .01472 L
.05293 .01472 L
.06244 .01472 L
.07283 .01472 L
.0782 .01472 L
.0841 .01472 L
.08889 .01472 L
.09404 .01472 L
.09948 .01472 L
.10458 .01472 L
.10712 .01472 L
.10982 .01472 L
.11237 .01472 L
.1135 .01472 L
.11468 .01472 L
.11536 .01472 L
.11611 .01472 L
.11741 .01472 L
.11815 .01472 L
.11894 .01472 L
.11968 .01557 L
.12036 .01723 L
.12168 .02163 L
.12311 .02756 L
.12571 .04023 L
.1458 .15347 L
.16558 .24403 L
.17556 .27924 L
.18659 .31111 L
.20593 .35277 L
.21681 .36992 L
.22679 .38254 L
.23594 .39192 L
.24571 .4 L
.25644 .40693 L
.26641 .41184 L
.27666 .4156 L
.28243 .41723 L
.28776 .41845 L
.29282 .41938 L
.29755 .42007 L
.30248 .42062 L
.3051 .42084 L
.30791 .42103 L
.31046 .42116 L
Mistroke
.31188 .42122 L
.31317 .42126 L
.31438 .42128 L
.31568 .42131 L
.31691 .42132 L
.31804 .42132 L
.31868 .42132 L
.31937 .42132 L
.32061 .42131 L
.32136 .4213 L
.32208 .42129 L
.32342 .42126 L
.32483 .42122 L
.32613 .42117 L
.32909 .42104 L
.33373 .42076 L
.33884 .42034 L
.34919 .41922 L
.36857 .41627 L
.38944 .41216 L
.43064 .40226 L
.47031 .39155 L
.50847 .38087 L
.54909 .36953 L
.58818 .35886 L
.62973 .34793 L
.66977 .33786 L
.70829 .32862 L
.74926 .31928 L
.78871 .31076 L
.82665 .30297 L
.86704 .2951 L
.90592 .28792 L
.94725 .28068 L
.97619 .27584 L
Mfstroke
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
  ImageSize->{439.375, 271.562},
  ImageMargins->{{43, 0}, {0, 0}},
  ImageRegion->{{0, 1}, {0, 1}},
  ImageCache->GraphicsData["Bitmap", "\<\
CF5dJ6E]HGAYHf4PAg9QL6QYHg<PAVmbKF5d0`40006g00013b000`400?l00000o`00003ooooooonh
ooooo`00oooooonhooooo`00oooooonhooooo`00J_ooool5o`0004cooooo0ol0001<ooooo`?o0000
COooool3o`0004[ooooo0ol00004ooooo`;o00001oooool006_ooooo00Co0000oooooooooooo0000
COooool00ol0003oooooooooo`1:ooooo`05o`000?ooooooooooooooool00000Boooool01Ol0003o
oooooooooooooooo000004[ooooo00?o0000ooooooooool00_ooool01?l0003oooooooooool00006
ooooo`00K?ooool00ol0003oooooooooo`1:ooooo`Go0000Boooool01Ol0003ooooooooooooooooo
000004_ooooo00Go0000ooooooooooooooooo`00001:ooooo`03o`000?oooooooooo00;ooooo00Co
0000oooooooooooo00001_ooool006gooooo00?o0000ooooooooool0BOooool01?l0003ooooooooo
ool0001<ooooo`Co0000COooool3o`0004_ooooo00?o0000ooooooooool00_ooool01?l0003ooooo
ooooool00006ooooo`00J_ooool01Ol0003ooooooooooooooooo000004_ooooo00?o0000ooooool0
0000C?ooool00ol0003oooooooooo`1=ooooo`05o`000?ooooooooooooooool00000B_ooool00ol0
003oooooooooo`02ooooo`04o`000?ooooooooooo`0000Kooooo001Zooooo`05o`000?oooooooooo
ooooool00000C?ooool2o`0004gooooo00?o0000ooooooooool0C?ooool01Ol0003ooooooooooooo
oooo000004Wooooo0_l00004ooooo`04o`000?ooooooooooo`0000Kooooo001[ooooo`?o0000C_oo
ool00ol0003oooooooooo`1<ooooo`;o0000COooool3o`0004_ooooo00?o0000ooooooooool00ooo
ool2o`0000Oooooo000Mooooo`03o`000?oooooooooo0?ooooooV?ooool001gooooo00?o0000oooo
ooooool0oooooonHooooo`007Oooool00ol0003oooooooooo`3ooooooiSooooo000Mooooo`03o`00
0?oooooooooo0?ooooooV?ooool001gooooo00?o0000ooooooooool0oooooonHooooo`007Oooool0
0ol0003oooooooooo`3ooooooiSooooo000Cooooo`Wo0000:_oo003oo`0007;o0000000Looooob[o
o`004oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool04?ooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool04Oooool00ol0003oooooooooo`0Aooooo`03o`000?oooooooooo017ooooo00?o0000oooooooo
ool02?ooool001gooooo00?o0000ooooooooool09?ooool2ool002Oooooo00?o0000ooooooooool0
C?ooool00ol0003oooooooooo`1=ooooo`03o`000?oooooooooo04gooooo00?o0000ooooooooool0
COooool00ol0003oooooooooo`08ooooo`007Oooool00ol0003oooooooooo`0Tooooo`?oo`009_oo
ool00ol0003oooooooooo`1<ooooo`03o`000?oooooooooo04gooooo00?o0000ooooooooool0COoo
ool00ol0003oooooooooo`1=ooooo`03o`000?oooooooooo00Sooooo000Mooooo`03o`000?oooooo
oooo02Gooooo0_oo003oooooog7ooooo000Mooooo`03o`000?oooooooooo02Gooooo0_oo003ooooo
of7ooooo1_l00?l:ooooo`007Oooool00ol0003oooooooooo`0Uooooo`;oo`00oooooomIooooo`ko
003o2_ooool001gooooo0_l0000Vooooo`;oo`00oooooomAoooooa7o003o3oooool001gooooo00?o
0000ooooooooool09Oooool2ool00?ooooooB?oooolBo`00oaOooooo000Mooooo`03o`000?oooooo
oooo02Gooooo0ooo003ooooood3ooooo4Ol00?lOooooo`007Oooool00ol0003oooooooooo`0Voooo
o`;oo`00ooooooljooooo`ko003o:?ooool001gooooo00?o0000ooooooooool09_ooool2ool00?oo
oooo=?ooool=o`00oboooooo000Mooooo`03o`000?oooooooooo02Kooooo0_oo003oooooobgooooo
3_l00?leooooo`007Oooool00ol0003oooooooooo`0Vooooo`;oo`00oooooolUoooooa3o003o>ooo
ool001gooooo0_l0000Wooooo`;oo`00oooooolNoooooa3o003o@_ooool001gooooo00?o0000oooo
ooooool09_ooool3ool00?oooooo5oooool>o`00od[ooooo000Mooooo`03o`000?oooooooooo02Oo
oooo0_oo003ooooooa;ooooo3?l00?mAooooo`007Oooool00ol0003oooooooooo`0Wooooo`;oo`00
ooooool=ooooo`_o003oEoooool001gooooo00?o0000ooooooooool09oooool2ool00?oooooo2Ooo
ool:o`00oecooooo000Mooooo`03o`000?oooooooooo02Oooooo0_oo003oooooo`Gooooo2Ol00?mQ
ooooo`007Oooool00ol0003oooooooooo`0Wooooo`?oo`00ooooool9o`00ofGooooo000Mooooo`;o
0000:Oooool2ool00?_ooooo2Ol00?mYooooo`007Oooool00ol0003oooooooooo`0Xooooo`;oo`00
moooool9o`00ofgooooo000Mooooo`03o`000?oooooooooo02Sooooo0_oo003booooo`[o003oLOoo
ool001gooooo00?o0000ooooooooool0:?ooool2ool00>kooooo2_l00?meooooo`007Oooool00ol0
003oooooooooo`0Xooooo`?oo`00jOooool9o`00og[ooooo000Mooooo`03o`000?oooooooooo02Wo
oooo0_oo003Uooooo`Wo003oO_ooool001gooooo00?o0000ooooooooool0:Oooool2ool00>3ooooo
2_l00?n2ooooo`007Oooool2o`0002[ooooo0_oo003Looooo`[o003oQ_ooool001gooooo00?o0000
ooooooooool0:Oooool2ool00=Wooooo2?l00?n;ooooo`007Oooool00ol0003oooooooooo`0Yoooo
o`?oo`00eOooool7o`00ohoooooo000Mooooo`03o`000?oooooooooo02[ooooo0_oo003Booooo`Oo
003oT_ooool001gooooo00?o0000ooooooooool0:_ooool2ool00=3ooooo1_l00?nEooooo`000_oo
ool2o`0000Cooooo00Co0000oooooooooooooooo1Ol00002ooooo`?o00001oooool00ol0003ooooo
ooooo`0Zooooo`;oo`00cOooool6o`00oiSooooo00001Ooooooo0000oooooooooooo000000Sooooo
00Ko0000oooooooooooo0000ooooool00003ooooo`03o`000?oooooooooo00Cooooo00?o0000oooo
ooooool0:_ooool2ool00<[ooooo1ol00?nJooooo`0000Goooooo`000?ooooooooooo`000009oooo
o`03o`000?oooooooooo00Gooooo00?o0000ooooooooool01?ooool3o`0002[ooooo0ooo0037oooo
o`Ko003oWOooool00005ooooool0003oooooooooool000002_ooool00ol0003oooooooooo`04oooo
o`03o`000?oooooooooo00Cooooo00?o0000ooooooooool0:oooool2ool00<Cooooo1_l00?nPoooo
o`0000Goooooo`000?ooooooooooo`000007ooooo`06o`000?ooooooooooooooool0003ooooo1?l0
0007ooooo`03o`000?oooooooooo02_ooooo0_oo0031ooooo`Oo003oX_ooool00005ooooool0003o
ooooooooool000001oooool01ol0003ooooooooooooooooo0000ooooool000002_ooool00ol0003o
ooooooooo`0[ooooo`;oo`00__ooool7o`00ojGooooo0002ooooo`;o00002Oooool3o`0000;ooooo
1Ol00006ooooo`03o`000?oooooooooo02_ooooo0_oo002kooooo`Oo003oZ?ooool001gooooo00?o
0000ooooooooool0:oooool2ool00;Sooooo1ol00?n[ooooo`007Oooool00ol0003oooooooooo`0[
ooooo`?oo`00]Oooool6o`00ojkooooo000Mooooo`;o0000;Oooool2ool00;;ooooo1_l00?naoooo
o`007Oooool00ol0003oooooooooo`0/ooooo`;oo`00[oooool7o`00ok?ooooo000Mooooo`03o`00
0?oooooooooo02cooooo0_oo002]ooooo`Ko003o]_ooool001gooooo00?o0000ooooooooool0;?oo
ool2ool00:[ooooo1_l00?niooooo`007Oooool00ol0003oooooooooo`0/ooooo`?oo`00Y_ooool7
o`00ok_ooooo000Mooooo`03o`000?oooooooooo02gooooo0_oo002Tooooo`Ko003o__ooool001go
oooo00?o0000ooooooooool0;Oooool2ool00:;ooooo1Ol00?o1ooooo`007Oooool00ol0003ooooo
ooooo`0]ooooo`;oo`00X?ooool5o`00ol?ooooo000Mooooo`;o0000;_ooool2ool009kooooo1Ol0
0?o5ooooo`007Oooool00ol0003oooooooooo`0]ooooo`?oo`00Voooool5o`00olOooooo000Moooo
o`03o`000?oooooooooo02kooooo0_oo002Jooooo`Co003obOooool001gooooo00?o0000oooooooo
ool0;_ooool2ool009Sooooo1?l00?o;ooooo`007Oooool00ol0003oooooooooo`0^ooooo`;oo`00
U_ooool5o`00olcooooo000Mooooo`03o`000?oooooooooo02kooooo0_oo002Dooooo`Go003oc_oo
ool001gooooo00?o0000ooooooooool0;_ooool3ool0097ooooo1Ol00?o@ooooo`007Oooool2o`00
033ooooo0_oo002?ooooo`Go003od_ooool001gooooo00?o0000ooooooooool0;oooool2ool008ko
oooo1?l00?oDooooo`007Oooool00ol0003oooooooooo`0_ooooo`;oo`00S?ooool4o`00omKooooo
000Mooooo`03o`000?oooooooooo02oooooo0_oo002:ooooo`Go003oeoooool001gooooo00?o0000
ooooooooool0;oooool3ool008Sooooo1?l00?oIooooo`007Oooool00ol0003oooooooooo`0`oooo
o`;oo`00Q_ooool4o`00om_ooooo000Mooooo`03o`000?oooooooooo033ooooo0_oo0024ooooo`Go
003og?ooool001gooooo0_l0000aooooo`;oo`00Poooool4o`00omkooooo000Mooooo`03o`000?oo
oooooooo033ooooo0_oo0021ooooo`Co003oh?ooool001gooooo00?o0000ooooooooool0<?ooool3
ool007kooooo1Ol00?oQooooo`007Oooool00ol0003oooooooooo`0aooooo`;oo`00OOooool4o`00
on?ooooo0008ooooo`;o00001?ooool00ol0003oooooooooo`02ooooo`?o00001oooool00ol0003o
ooooooooo`0aooooo`;oo`00Noooool4o`00onGooooo0007ooooo`04o`000?ooooooooooo`0000Oo
oooo00Go0000ooooooooooooooooo`000006ooooo`03o`000?oooooooooo037ooooo0_oo001joooo
o`Co003oi_ooool000Oooooo00Co0000oooooooooooo00002oooool00ol0003oooooooooo`04oooo
o`03o`000?oooooooooo037ooooo0_oo001hooooo`Co003oj?ooool000Oooooo00Co0000oooooooo
oooo00002oooool00ol0003oooooooooo`04ooooo`?o0000<Oooool3ool007Kooooo1?l00?oYoooo
o`001oooool01?l0003oooooooooool00007ooooo`Co00001oooool00ol0003oooooooooo`0boooo
o`;oo`00MOooool3o`00on_ooooo0007ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000oooo
ooooool02?ooool00ol0003oooooooooo`0booooo`;oo`00Loooool4o`00oncooooo0008ooooo`;o
00002?ooool5o`0000Kooooo00?o0000ooooooooool0<_ooool2ool007;ooooo1?l00?o]ooooo`00
7Oooool00ol0003oooooooooo`0booooo`;oo`00L?ooool4o`00onoooooo000Mooooo`03o`000?oo
oooooooo03;ooooo0ooo001^ooooo`Co003ol?ooool001gooooo00?o0000ooooooooool0<oooool2
ool006cooooo1?l00?obooooo`007Oooool2o`0003Cooooo0_oo001[ooooo`Co003oloooool001go
oooo00?o0000ooooooooool0<oooool2ool006[ooooo0ol00?oeooooo`007Oooool00ol0003ooooo
ooooo`0cooooo`?oo`00J?ooool3o`00ooKooooo000Mooooo`03o`000?oooooooooo03Cooooo0_oo
001Wooooo`?o003omoooool001gooooo00?o0000ooooooooool0=?ooool2ool006Gooooo1?l00?oh
ooooo`007Oooool00ol0003oooooooooo`0dooooo`;oo`00I?ooool4o`00ooWooooo000Mooooo`03
o`000?oooooooooo03Cooooo0_oo001Sooooo`?o003onoooool001gooooo0_l0000eooooo`?oo`00
HOooool3o`00oocooooo000Mooooo`03o`000?oooooooooo03Gooooo0_oo001Oooooo`Co003ooOoo
ool001gooooo00?o0000ooooooooool0=Oooool2ool005kooooo1?l00?onooooo`007Oooool00ol0
003oooooooooo`0eooooo`;oo`00GOooool3o`00oooooooo0Oooool001gooooo00?o0000oooooooo
ool0=Oooool2ool005cooooo0ol00?ooooooo`;ooooo000Mooooo`03o`000?oooooooooo03Gooooo
0ooo001Jooooo`?o003oooooool3ooooo`007Oooool00ol0003oooooooooo`0fooooo`;oo`00FOoo
ool3o`00oooooooo1?ooool001gooooo0_l0000gooooo`;oo`00FOooool2o`00oooooooo1Oooool0
01gooooo00?o0000ooooooooool0=_ooool2ool005Sooooo0ol00?ooooooo`Gooooo000Mooooo`03
o`000?oooooooooo03Kooooo0_oo001Gooooo`?o003oooooool6ooooo`007Oooool00ol0003ooooo
ooooo`0fooooo`?oo`00EOooool3o`00oooooooo1oooool001gooooo00?o0000ooooooooool0=ooo
ool2ool005Cooooo0ol00?ooooooo`Sooooo000Mooooo`03o`000?oooooooooo03Oooooo0_oo001C
ooooo`?o003oooooool9ooooo`007Oooool00ol0003oooooooooo`0gooooo`;oo`00D_ooool3o`00
oooooooo2_ooool001gooooo0_l0000hooooo`?oo`00D?ooool3o`00oooooooo2oooool001gooooo
00?o0000ooooooooool0>?ooool2ool004oooooo0ol00?ooooooo`cooooo000Mooooo`03o`000?oo
oooooooo03Sooooo0_oo001>ooooo`?o003oooooool=ooooo`007Oooool00ol0003oooooooooo`0h
ooooo`;oo`00C_ooool2o`00oooooooo3_ooool001gooooo00?o0000ooooooooool0>?ooool3ool0
04cooooo0ol00?ooooooo`kooooo0002ooooo`;o00001?ooool00ol0003oooooooooo`03ooooo`03
o`000?oooooooooo00;ooooo0ol00007ooooo`03o`000?oooooooooo03Wooooo0_oo001;ooooo`?o
003oooooool?ooooo`0000Goooooo`000?ooooooooooo`000009ooooo`05o`000?oooooooooooooo
ool000000oooool00ol0003oooooooooo`04ooooo`03o`000?oooooooooo03Wooooo0_oo001:oooo
o`?o003oooooool@ooooo`0000Goooooo`000?ooooooooooo`00000:ooooo`03o`000?oooooooooo
00Cooooo00?o0000ooooooooool01?ooool3o`0003Wooooo0ooo0018ooooo`?o003ooooooolAoooo
o`0000Goooooo`000?ooooooooooo`00000:ooooo`03o`000?oooooooooo00Cooooo00?o0000oooo
ooooool01?ooool00ol0003oooooooooo`0jooooo`;oo`00Aoooool3o`00oooooooo4_ooool00005
ooooool0003oooooooooool000002oooool00ol0003oooooo`000003o`0000Oooooo00?o0000oooo
ooooool0>_ooool2ool004Kooooo0ol00?oooooooa?ooooo00001Ooooooo0000oooooooooooo0000
00Oooooo00Oo0000ooooooooooooooooo`000?oooooo000000[ooooo00?o0000ooooooooool0>_oo
ool2ool004Gooooo0ol00?ooooooo`Gooooo1Ooo000:ooooo`000_ooool2o`0000Sooooo1Ol00000
0ooooooo0000o`000003o`0000Kooooo00?o0000ooooooooool0>_ooool3ool004?ooooo0ol00?oo
ooooo`_oo`002_ooool001gooooo00?o0000ooooooooool0>oooool2ool004?ooooo0_l00?ojoooo
o`goo`003_ooool001gooooo00?o0000ooooooooool0>oooool2ool004;ooooo0ol00?odooooo`go
o`005?ooool001gooooo0_l0000looooo`?oo`00@?ooool3o`00onoooooo3Ooo000Jooooo`007Ooo
ool00ol0003oooooooooo`0looooo`;oo`00?oooool3o`00on[ooooo3Ooo000Pooooo`007Oooool0
0ol0003oooooooooo`0looooo`;oo`00?_ooool3o`00onKooooo3?oo000Vooooo`007Oooool00ol0
003oooooooooo`0looooo`?oo`00?Oooool2o`00on;ooooo2ooo000/ooooo`007Oooool00ol0003o
ooooooooo`0mooooo`;oo`00??ooool3o`00omgooooo2ooo000aooooo`007Oooool00ol0003ooooo
ooooo`0mooooo`?oo`00>_ooool3o`00om[ooooo2_oo000fooooo`007Oooool00ol0003ooooooooo
o`0nooooo`;oo`00>Oooool3o`00omKooooo2_oo000kooooo`007Oooool00ol0003oooooooooo`0n
ooooo`;oo`00>?ooool3o`00om?ooooo2_oo000oooooo`007Oooool2o`0003oooooo0ooo000foooo
o`?o003ocoooool:ool004Cooooo000Mooooo`03o`000?oooooooooo03oooooo0_oo000fooooo`;o
003oboooool;ool004Sooooo000Mooooo`03o`000?oooooooooo03oooooo0_oo000eooooo`?o003o
a_ooool;ool004gooooo000Mooooo`03o`000?oooooooooo03oooooo0ooo000cooooo`?o003o`_oo
ool;ool005;ooooo000Mooooo`03o`000?oooooooooo043ooooo0_oo000booooo`?o003o_oooool:
ool005Oooooo000Mooooo`03o`000?oooooooooo043ooooo0ooo000`ooooo`?o003o_?ooool9ool0
05cooooo000Mooooo`03o`000?oooooooooo047ooooo0_oo000`ooooo`;o003o^Oooool9ool0063o
oooo000Mooooo`;o0000@_ooool3ool002kooooo0ol00?ndooooo`[oo`00I?ooool001gooooo00?o
0000ooooooooool0@_ooool2ool002gooooo0ol00?n_ooooo`coo`00J?ooool001gooooo00?o0000
ooooooooool0@_ooool2ool002cooooo0ol00?n[ooooo`coo`00KOooool001gooooo00?o0000oooo
ooooool0@_ooool3ool002[ooooo0ol00?nWooooo`_oo`00Loooool001gooooo00?o0000oooooooo
ool0@oooool2ool002[ooooo0_l00?nTooooo`[oo`00N?ooool001gooooo00?o0000ooooooooool0
@oooool3ool002Sooooo0ol00?nOooooo`[oo`00OOooool001gooooo00?o0000ooooooooool0A?oo
ool2ool002Oooooo0ol00?nLooooo`[oo`00POooool001gooooo0_l00015ooooo`?oo`009_ooool2
o`00oiWooooo2Ooo0026ooooo`007Oooool00ol0003oooooooooo`15ooooo`;oo`009Oooool3o`00
oiKooooo2?oo002:ooooo`007Oooool00ol0003oooooooooo`15ooooo`?oo`009?ooool2o`00oiCo
oooo1ooo002>ooooo`007Oooool00ol0003oooooooooo`16ooooo`;oo`008oooool3o`00oi7ooooo
1ooo002Aooooo`004oooool3o`0000Oooooo00?o0000ooooooooool0A_ooool2ool002;ooooo0ol0
0?n?ooooo`Ooo`00U?ooool001Cooooo00?o0000ooooooooool01_ooool00ol0003oooooooooo`16
ooooo`?oo`008Oooool2o`00ohcooooo2?oo002Gooooo`005?ooool00ol0003oooooooooo`06oooo
o`03o`000?oooooooooo04Oooooo0_oo000Pooooo`?o003oQoooool:ool009[ooooo000Dooooo`03
o`000?oooooooooo00Gooooo0ol00?l00ol0003oooooooooo`15ooooo`?oo`007oooool2o`00ohCo
oooo2_oo002Nooooo`005?ooool00ol0003oooooooooo`05ooooo`Co003oB?ooool2ool001kooooo
0ol00?moooooo`[oo`00Xoooool001?ooooo0_l00008ooooo`04o`000?l00?oo003oo`00odOooooo
0ooo000Looooo`?o003oO?ooool:ool00:Oooooo000Dooooo`03o`000?oooooooooo00Kooooo00?o
0000ooooool00?l00_l00?m7ooooo`?oo`006oooool2o`00og[ooooo2?oo002/ooooo`007Oooool0
0ol0003oooooooooo`03o`00odOooooo0ooo000Iooooo`?o003oMoooool7ool00;3ooooo000Moooo
o`04o`000?ooooooooooooooo`?o003oAoooool2ool001Wooooo0_l00?meooooo`Ooo`00/oooool0
01gooooo00?o0000ooooooooool00_ooool3o`00odKooooo0ooo000Gooooo`?o003oL_ooool7ool0
0;Kooooo000Mooooo`;o00001?ooool3o`00odKooooo0ooo000Eooooo`?o003oKoooool8ool00;Wo
oooo000Mooooo`03o`000?oooooooooo00Cooooo0_l00?m7ooooo`;oo`005Oooool2o`00ofcooooo
2Ooo002looooo`007Oooool00ol0003oooooooooo`04ooooo`?o003oA_ooool3ool001?ooooo0ol0
0?mWooooo`[oo`00`?ooool001gooooo00?o0000ooooooooool01Oooool3o`00odKooooo0ooo000B
ooooo`;o003oI?ooool:ool00<Cooooo000Mooooo`03o`000?oooooooooo00Kooooo0ol00?m6oooo
o`?oo`004?ooool3o`00of3ooooo2Ooo0039ooooo`007Oooool00ol0003oooooooooo`07ooooo`?o
003oA_ooool3ool000kooooo0ol00?mNooooo`Soo`00cOooool001gooooo00?o0000ooooooooool0
2?ooool3o`00odKooooo0ooo000=ooooo`;o003oG?ooool7ool00=7ooooo000Mooooo`;o00002_oo
ool3o`00odKooooo0ooo000;ooooo`?o003oFOooool7ool00=Cooooo000Mooooo`03o`000?oooooo
oooo00[ooooo0ol00?m6ooooo`?oo`002_ooool2o`00oeOooooo1ooo003Gooooo`007Oooool00ol0
003oooooooooo`0;ooooo`?o003oA_ooool4ool000Oooooo0ol00?mCooooo`Soo`00f_ooool001go
oooo00?o0000ooooooooool03?ooool2o`00odOooooo1?oo0005ooooo`?o003oD?ooool9ool00=go
oooo000Mooooo`03o`000?oooooooooo00cooooo0ol00?m8ooooo`?oo`001?ooool2o`00odcooooo
2_oo003Qooooo`007Oooool00ol0003oooooooooo`0=ooooo`?o003oB?ooool4ool00004ooooool0
0?oo003oo`00odSooooo2_oo003Uooooo`007Oooool00ol0003oooooooooo`0>ooooo`;o003oBOoo
ool4ool000;o003oAOooool9ool00>[ooooo000Mooooo`;o00003oooool3o`00od[ooooo0ooo0000
0ol00?ooooooooooo`0oooooo`Woo`00k_ooool001gooooo00?o0000ooooooooool03oooool3o`00
od[ooooo1?oo000kooooo`[oo`00l_ooool001gooooo00?o0000ooooooooool04?ooool2o`00odWo
oooo0_l00?l5ool003Gooooo2_oo003fooooo`007Oooool00ol0003oooooooooo`0@ooooo`?o003o
B?ooool2o`00o`;ooooo1Ooo000_ooooo`Woo`00noooool001gooooo00?o0000ooooooooool04Ooo
ool3o`00odKooooo0ol00?l4ooooo`Koo`009oooool:ool00?oooooo000Mooooo`03o`000?oooooo
oooo01;ooooo0_l00?m6ooooo`;o003o1oooool6ool001oooooo3?oo003oooooo`Cooooo000Moooo
o`03o`000?oooooooooo01;ooooo0ol00?m4ooooo`?o003o2_ooool7ool001Gooooo3Ooo003ooooo
o`Wooooo000Mooooo`;o00005?ooool3o`00od?ooooo0_l00?l=ooooob7oo`00ooooool?ooooo`00
7Oooool00ol0003oooooooooo`0Dooooo`;o003o@_ooool3o`00oa7ooooo5ooo003ooooooaGooooo
000Mooooo`03o`000?oooooooooo01Cooooo0ol00?m1ooooo`;o003ooooooolnooooo`007Oooool0
0ol0003oooooooooo`0Eooooo`;o003o@?ooool3o`00oooooooo?_ooool001gooooo00?o0000oooo
ooooool05Oooool3o`00ocoooooo0_l00?ooooooocoooooo00001?oooooo0000o`000?l00004oooo
o`04o`000?ooooooooooooooo`Go00000_ooool3o`0000Oooooo00?o0000ooooooooool05_ooool3
o`00ocgooooo0ol00?ooooooocoooooo0002ooooo`03o`000?oooooooooo00Sooooo00Ko0000oooo
oooooooo0000ooooool00003ooooo`03o`000?oooooooooo00Cooooo00?o0000ooooooooool05ooo
ool2o`00occooooo0ol00?oooooood3ooooo0002ooooo`03o`000?oooooooooo00Wooooo00?o0000
ooooooooool01Oooool00ol0003oooooooooo`04ooooo`?o00005oooool3o`00oc_ooooo0_l00?oo
oooood7ooooo0002ooooo`03o`000?oooooooooo00[ooooo00?o0000ooooooooool01?ooool00ol0
003oooooooooo`04ooooo`03o`000?oooooooooo01Sooooo0_l00?ljooooo`?o003ooooooom1oooo
o`000_ooool00ol0003oooooooooo`07ooooo`06o`000?ooooooooooooooool0003ooooo1?l00007
ooooo`03o`000?oooooooooo01Sooooo0ol00?liooooo`;o003ooooooom2ooooo`0000?oooooo`00
0?l000002Oooool01ol0003ooooooooooooooooo0000ooooool000002_ooool00ol0003ooooooooo
o`0Iooooo`;o003o>?ooool3o`00oooooooo@_ooool000;ooooo00?o0000ooooooooool02?ooool3
o`0000;ooooo1Ol00006ooooo`03o`000?oooooooooo01Wooooo0ol00?lgooooo`;o003ooooooom3
ooooo`007Oooool00ol0003oooooooooo`0Jooooo`;o003o=_ooool3o`00oooooooo@oooool001go
oooo00?o0000ooooooooool06_ooool2o`00ocKooooo0_l00?ooooooodCooooo000Mooooo`;o0000
6oooool3o`00ocCooooo0ol00?ooooooodCooooo000Mooooo`03o`000?oooooooooo01_ooooo0_l0
0?ldooooo`;o003ooooooom5ooooo`007Oooool00ol0003oooooooooo`0Kooooo`?o003o<_ooool3
o`00ooooooooAOooool001gooooo00?o0000ooooooooool07?ooool3o`00oc7ooooo0_l00?oooooo
odKooooo000Mooooo`03o`000?oooooooooo01gooooo0_l00?l`ooooo`?o003ooooooom6ooooo`00
7Oooool00ol0003oooooooooo`0Mooooo`?o003o;oooool2o`00ooooooooAoooool001gooooo00?o
0000ooooooooool07_ooool2o`00obkooooo0ol00?ooooooodOooooo000Mooooo`;o00007oooool2
o`00obkooooo0_l00?ooooooodSooooo000Mooooo`03o`000?oooooooooo01kooooo0ol00?l/oooo
o`?o003ooooooom8ooooo`007Oooool00ol0003oooooooooo`0Oooooo`;o003o;?ooool2o`00oooo
ooooBOooool001gooooo00?o0000ooooooooool07oooool2o`00ob_ooooo0ol00?ooooooodWooooo
000Mooooo`03o`000?oooooooooo01oooooo0ol00?lZooooo`;o003ooooooom:ooooo`007Oooool0
0ol0003oooooooooo`0Pooooo`;o003o:_ooool2o`00ooooooooB_ooool001gooooo00?o0000oooo
ooooool08?ooool2o`00obWooooo0ol00?oooooood[ooooo000Mooooo`03o`000?oooooooooo023o
oooo0_l00?lYooooo`;o003ooooooom;ooooo`007Oooool2o`00027ooooo0ol00?lWooooo`?o003o
oooooom;ooooo`007Oooool00ol0003oooooooooo`0Qooooo`;o003o9oooool2o`00ooooooooC?oo
ool001gooooo00?o0000ooooooooool08Oooool2o`00obKooooo0ol00?ooooooodcooooo000Moooo
o`03o`000?oooooooooo027ooooo0_l00?lVooooo`;o003ooooooom=ooooo`007Oooool00ol0003o
ooooooooo`0Qooooo`?o003o9?ooool3o`00ooooooooCOooool001gooooo00?o0000ooooooooool0
8_ooool2o`00obCooooo0_l00?ooooooodkooooo000Mooooo`03o`000?oooooooooo02;ooooo0_l0
0?lTooooo`;o003ooooooom>ooooo`007Oooool2o`0002?ooooo0_l00?lSooooo`?o003ooooooom>
ooooo`007Oooool00ol0003oooooooooo`0Rooooo`?o003o8_ooool2o`00ooooooooCoooool001go
oooo00?o0000ooooooooool08oooool2o`00ob7ooooo0ol00?ooooooodoooooo000Mooooo`03o`00
0?oooooooooo02?ooooo0_l00?lQooooo`;o003ooooooom@ooooo`001oooool3o`0000Cooooo00?o
0000ooooooooool00_ooool3o`0000Oooooo00?o0000ooooooooool08oooool2o`00ob3ooooo0ol0
0?oooooooe3ooooo0008ooooo`03o`000?oooooooooo00Oooooo00Go0000ooooooooooooooooo`00
0006ooooo`03o`000?oooooooooo02?ooooo0_l00?lPooooo`;o003ooooooomAooooo`002?ooool0
0ol0003oooooooooo`0;ooooo`03o`000?oooooooooo00Cooooo00?o0000ooooooooool08oooool2
o`00oaoooooo0ol00?oooooooe7ooooo0008ooooo`03o`000?oooooooooo00_ooooo00?o0000oooo
ooooool01?ooool3o`0002?ooooo0ol00?lNooooo`;o003ooooooomBooooo`002?ooool00ol0003o
ooooooooo`07ooooo`Co00001oooool00ol0003oooooooooo`0Tooooo`;o003o7_ooool2o`00oooo
ooooD_ooool000Oooooo0_l00009ooooo`03o`000?oooooooooo00Sooooo00?o0000ooooooooool0
9?ooool2o`00oagooooo0ol00?oooooooe;ooooo0008ooooo`03o`000?oooooooooo00Oooooo1Ol0
0006ooooo`03o`000?oooooooooo02Cooooo0_l00?lMooooo`;o003ooooooomCooooo`007Oooool0
0ol0003oooooooooo`0Tooooo`;o003o7?ooool3o`00ooooooooDoooool001gooooo00?o0000oooo
ooooool09?ooool2o`00oacooooo0_l00?oooooooeCooooo000Mooooo`03o`000?oooooooooo02Co
oooo0_l00?lKooooo`?o003ooooooomDooooo`007Oooool2o`0002Gooooo0ol00?lJooooo`;o003o
oooooomEooooo`007Oooool00ol0003oooooooooo`0Uooooo`;o003o6Oooool3o`00ooooooooEOoo
ool001gooooo00?o0000ooooooooool09Oooool2o`00oaWooooo0_l00?oooooooeKooooo000Moooo
o`03o`000?oooooooooo02Gooooo0_l00?lIooooo`;o003ooooooomFooooo`007Oooool00ol0003o
ooooooooo`0Uooooo`;o003o6?ooool3o`00ooooooooE_ooool001gooooo00?o0000ooooooooool0
9Oooool2o`00oaSooooo0_l00?oooooooeOooooo000Mooooo`03o`000?oooooooooo02Gooooo0_l0
0?lGooooo`?o003ooooooomGooooo`007Oooool2o`0002Kooooo0ol00?lFooooo`;o003ooooooomH
ooooo`007Oooool00ol0003oooooooooo`0Vooooo`;o003o5Oooool3o`00ooooooooF?ooool001go
oooo00?o0000ooooooooool09_ooool2o`00oaGooooo0_l00?oooooooeWooooo000Mooooo`03o`00
0?oooooooooo02Kooooo0_l00?lDooooo`?o003ooooooomIooooo`007Oooool00ol0003ooooooooo
o`0Vooooo`;o003o4oooool3o`00ooooooooF_ooool001gooooo00?o0000ooooooooool09_ooool3
o`00oa;ooooo0_l00?oooooooe_ooooo000Mooooo`03o`000?oooooooooo02Oooooo0_l00?lAoooo
o`?o003ooooooomKooooo`007Oooool2o`0002Sooooo0_l00?lAooooo`;o003ooooooomLooooo`00
7Oooool00ol0003oooooooooo`0Wooooo`;o003o4?ooool3o`00ooooooooG?ooool001gooooo00?o
0000ooooooooool09oooool3o`00o`oooooo0_l00?oooooooegooooo000Mooooo`03o`000?oooooo
oooo02Sooooo0_l00?l>ooooo`?o003ooooooomMooooo`007Oooool00ol0003oooooooooo`0Xoooo
o`;o003o3Oooool3o`00ooooooooG_ooool001gooooo00?o0000ooooooooool0:?ooool3o`00o`co
oooo0_l00?oooooooeoooooo000Mooooo`03o`000?oooooooooo02Wooooo0_l00?l;ooooo`?o003o
oooooomOooooo`007Oooool2o`0002[ooooo0_l00?l;ooooo`;o003ooooooomPooooo`007Oooool0
0ol0003oooooooooo`0Yooooo`?o003o2Oooool3o`00ooooooooH?ooool001gooooo00?o0000oooo
ooooool0:_ooool2o`00o`Sooooo0ol00?ooooooof7ooooo000Mooooo`03o`000?oooooooooo02[o
oooo0ol00?l6ooooo`?o003ooooooomRooooo`007Oooool00ol0003oooooooooo`0[ooooo`?o003o
1?ooool3o`00ooooooooHoooool00004ooooool0003o0000o`0000Cooooo00?o0000ooooooooool0
0oooool00ol0003oooooooooo`02ooooo`?o00001oooool00ol0003oooooooooo`0/ooooo`So003o
oooooomTooooo`000_ooool00ol0003oooooooooo`09ooooo`05o`000?ooooooooooooooool00000
0oooool00ol0003oooooooooo`04ooooo`03o`000?oooooooooo02gooooo1_l00?ooooooofGooooo
0002ooooo`03o`000?oooooooooo00[ooooo00?o0000ooooooooool01?ooool00ol0003ooooooooo
o`04ooooo`?o0000;_ooool3o`00ooooooooIoooool000;ooooo00?o0000ooooooooool02_ooool0
0ol0003oooooooooo`04ooooo`03o`000?oooooooooo00Cooooo00?o0000ooooooooool0oooooonH
ooooo`000_ooool00ol0003oooooooooo`0;ooooo`03o`000?oooooo000000?o00001oooool00ol0
003oooooooooo`3ooooooiSooooo00000ooooooo0000o`000009ooooo`07o`000?oooooooooooooo
ool0003oooooo`00000:ooooo`03o`000?oooooooooo0?ooooooV?ooool000;ooooo00?o0000oooo
ooooool01oooool5o`000003ooooool0003o000000?o00001_ooool00ol0003oooooooooo`3ooooo
oiSooooo000Mooooo`03o`000?oooooooooo0?ooooooV?ooool00?oooooo^?ooool00001\
\>"],
  ImageRangeCache->{{{0, 438.375}, {270.562, 0}} -> {-0.732569, -0.112555, \
0.025053, 0.00704711}}]
}, Open  ]],

Cell["\<\
Plot of the real and imaginary part of 3*I(a) as defined in the \
notes (I use 3*I(a) so that the function goes to 1 as a->0.)\
\>", "Text"]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
  "Function to be integrated to get the LO cross section. In general, we use \
the convention that all functions to be numerically integrated in the \
hypercubes ",
  Cell[BoxData[
      \(TraditionalForm\`\([0, 1]\)\^d\)]],
  ", where d is the dimension of the integration space. dsigma depends on the \
Higgs mass (mh), the total collider energy (sqrtS) and from the arbitrary \
normalization and factorization scales (mur,muf)."
}], "Subsubsection"],

Cell["\<\
dsigmaLO[zz_,mh_,sqrtS_,muf_,mur_]:=Module[
{y,x1,x2,gg0,qg0,qq0,s0,ymin,ymax,JAC,v,S,tau0},
Muf=muf*1.;
v=246.;
S=sqrtS^2;
tau0=mh^2/S;
ymax=-Log[Sqrt[tau0]];
ymin=-ymax;
y=ymin+(ymax-ymin)*zz;
JAC=ymax-ymin;
x1=Sqrt[tau0] Exp[y];
x2=Sqrt[tau0] Exp[-y];
{gg0,qg0,qq0}=pdfCTEQLO[x1,x2,Muf];
s0=asLO[mur,5]^2/576/Pi/v^2*tau0;
s0=s0*gg0;
s0=s0*389379660; (*to picobarns*)
s0=s0*JAC;
Return[s0];
];\
\>", "Input"],

Cell[CellGroupData[{

Cell["sLO=NIntegrate[dsigmaLO[xx,100,14000,100,100],{xx,0,1}]", "Input"],

Cell[BoxData[
    \(TraditionalForm\`28.056643821030747`\)], "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Get the cross section for various Higgs masses and including the \
form factor of the loop.
I first build a table with the results of the cross section in picobarns, and \
then plot it.\
\>", "Subsection"],

Cell["\<\
resEFT=Table[{i,NIntegrate[dsigmaLO[x,i*1.,14000.,i*1.,i*1.],{x,0,1}\
]},{i,20,600,10}];\
\>", "Input"],

Cell["\<\
resFULL=Table[{resEFT[[i]][[1]],Abs[inte[resEFT[[i]][[1]]^2/175^2/4]\
]^2*resEFT[[i]][[2]]},{i,1,Length[resEFT]}];\
\>", "Input"],

Cell["resFULL", "Input"],

Cell[CellGroupData[{

Cell[TextData[{
  "Show[LogListPlot[resEFT, ",
  StyleBox["PlotJoined", "MR"],
  " ",
  StyleBox["->", "MR"],
  " ",
  StyleBox["True,", "MR"],
  "PlotStyle\[Rule]{Red}],\nplotFULL=LogListPlot[resFULL, ",
  StyleBox["PlotJoined", "MR"],
  " ",
  StyleBox["->True,", "MR"],
  "PlotStyle\[Rule]{Blue}]];\n\n"
}], "Input"],

Cell[GraphicsData["PostScript", "\<\
%!
%%Creator: Mathematica
%%AspectRatio: .61803 
MathPictureStart
/Mabs {
Mgmatrix idtransform
Mtmatrix dtransform
} bind def
/Mabsadd { Mabs
3 -1 roll add
3 1 roll add
exch } bind def
%% Graphics
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10  scalefont  setfont
% Scaling calculations
0.0238095 0.0015873 0.078927 0.221776 [
[.02381 -0.0125 -3 -9 ]
[.02381 -0.0125 3 0 ]
[.18254 -0.0125 -9 -9 ]
[.18254 -0.0125 9 0 ]
[.34127 -0.0125 -9 -9 ]
[.34127 -0.0125 9 0 ]
[.5 -0.0125 -9 -9 ]
[.5 -0.0125 9 0 ]
[.65873 -0.0125 -9 -9 ]
[.65873 -0.0125 9 0 ]
[.81746 -0.0125 -9 -9 ]
[.81746 -0.0125 9 0 ]
[.97619 -0.0125 -9 -9 ]
[.97619 -0.0125 9 0 ]
[-0.0125 .07893 -10 -6 ]
[-0.0125 .07893 0 6 ]
[-0.0125 .14569 -10 -6 ]
[-0.0125 .14569 0 6 ]
[-0.0125 .23394 -10 -6 ]
[-0.0125 .23394 0 6 ]
[-0.0125 .3007 -16 -6 ]
[-0.0125 .3007 0 6 ]
[-0.0125 .36746 -16 -6 ]
[-0.0125 .36746 0 6 ]
[-0.0125 .45572 -16 -6 ]
[-0.0125 .45572 0 6 ]
[-0.0125 .52248 -22 -6 ]
[-0.0125 .52248 0 6 ]
[-0.0125 .58924 -22 -6 ]
[-0.0125 .58924 0 6 ]
[ -0.0005 0 0 0 ]
[ 1 .61803 0 0 ]
] MathScale
% Start of Graphics
1 setlinecap
1 setlinejoin
newpath
0 g
.25 Mabswid
[ ] 0 setdash
.02381 0 m
.02381 .00625 L
s
[(0)] .02381 -0.0125 0 1 Mshowa
.18254 0 m
.18254 .00625 L
s
[(100)] .18254 -0.0125 0 1 Mshowa
.34127 0 m
.34127 .00625 L
s
[(200)] .34127 -0.0125 0 1 Mshowa
.5 0 m
.5 .00625 L
s
[(300)] .5 -0.0125 0 1 Mshowa
.65873 0 m
.65873 .00625 L
s
[(400)] .65873 -0.0125 0 1 Mshowa
.81746 0 m
.81746 .00625 L
s
[(500)] .81746 -0.0125 0 1 Mshowa
.97619 0 m
.97619 .00625 L
s
[(600)] .97619 -0.0125 0 1 Mshowa
.125 Mabswid
.05556 0 m
.05556 .00375 L
s
.0873 0 m
.0873 .00375 L
s
.11905 0 m
.11905 .00375 L
s
.15079 0 m
.15079 .00375 L
s
.21429 0 m
.21429 .00375 L
s
.24603 0 m
.24603 .00375 L
s
.27778 0 m
.27778 .00375 L
s
.30952 0 m
.30952 .00375 L
s
.37302 0 m
.37302 .00375 L
s
.40476 0 m
.40476 .00375 L
s
.43651 0 m
.43651 .00375 L
s
.46825 0 m
.46825 .00375 L
s
.53175 0 m
.53175 .00375 L
s
.56349 0 m
.56349 .00375 L
s
.59524 0 m
.59524 .00375 L
s
.62698 0 m
.62698 .00375 L
s
.69048 0 m
.69048 .00375 L
s
.72222 0 m
.72222 .00375 L
s
.75397 0 m
.75397 .00375 L
s
.78571 0 m
.78571 .00375 L
s
.84921 0 m
.84921 .00375 L
s
.88095 0 m
.88095 .00375 L
s
.9127 0 m
.9127 .00375 L
s
.94444 0 m
.94444 .00375 L
s
.25 Mabswid
0 0 m
1 0 L
s
0 .07893 m
.00625 .07893 L
s
gsave
-0.0125 .07893 -71 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(1) show
69.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .14569 m
.00625 .14569 L
s
gsave
-0.0125 .14569 -71 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(2) show
69.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .23394 m
.00625 .23394 L
s
gsave
-0.0125 .23394 -71 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(5) show
69.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .3007 m
.00625 .3007 L
s
gsave
-0.0125 .3007 -77 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(10) show
75.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .36746 m
.00625 .36746 L
s
gsave
-0.0125 .36746 -77 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(20) show
75.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .45572 m
.00625 .45572 L
s
gsave
-0.0125 .45572 -77 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(50) show
75.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .52248 m
.00625 .52248 L
s
gsave
-0.0125 .52248 -83 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(100) show
81.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
0 .58924 m
.00625 .58924 L
s
gsave
-0.0125 .58924 -83 -10 Mabsadd m
1 1 Mabs scale
currentpoint translate
0 20 translate 1 -1 scale
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
63.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
(200) show
81.000 13.000 moveto
%%IncludeResource: font Courier
%%IncludeFont: Courier
/Courier findfont 10.000 scalefont
[1 0 0 -1 0 0 ] makefont setfont
0.000 0.000 0.000 setrgbcolor
0.000 0.000 rmoveto
1.000 setlinewidth
grestore
.001 w
0 .18474 m
.00375 .18474 L
s
0 .21245 m
.00375 .21245 L
s
0 .2515 m
.00375 .2515 L
s
0 .26635 m
.00375 .26635 L
s
0 .27921 m
.00375 .27921 L
s
0 .29056 m
.00375 .29056 L
s
0 .40652 m
.00375 .40652 L
s
0 .43423 m
.00375 .43423 L
s
0 .47328 m
.00375 .47328 L
s
0 .48813 m
.00375 .48813 L
s
0 .50099 m
.00375 .50099 L
s
0 .51233 m
.00375 .51233 L
s
.25 Mabswid
0 0 m
0 .61803 L
s
0 0 m
1 0 L
1 .61803 L
0 .61803 L
closepath
clip
newpath
1 0 0 r
.5 Mabswid
.05556 .60332 m
.07143 .56246 L
.0873 .52864 L
.10317 .49984 L
.11905 .47468 L
.13492 .4523 L
.15079 .43206 L
.16667 .41355 L
.18254 .39647 L
.19841 .38059 L
.21429 .36574 L
.23016 .35176 L
.24603 .33854 L
.2619 .326 L
.27778 .31405 L
.29365 .30264 L
.30952 .29171 L
.3254 .28121 L
.34127 .2711 L
.35714 .26135 L
.37302 .25194 L
.38889 .24283 L
.40476 .23401 L
.42063 .22545 L
.43651 .21714 L
.45238 .20905 L
.46825 .20117 L
.48413 .1935 L
.5 .18601 L
.51587 .17869 L
.53175 .17154 L
.54762 .16454 L
.56349 .1577 L
.57937 .15099 L
.59524 .14441 L
.61111 .13797 L
.62698 .13164 L
.64286 .12543 L
.65873 .11933 L
.6746 .11334 L
.69048 .10744 L
.70635 .10165 L
.72222 .09594 L
.7381 .09033 L
.75397 .0848 L
.76984 .07935 L
.78571 .07398 L
.80159 .06868 L
.81746 .06346 L
.83333 .05831 L
Mistroke
.84921 .05322 L
.86508 .0482 L
.88095 .04325 L
.89683 .03835 L
.9127 .03351 L
.92857 .02873 L
.94444 .02401 L
.96032 .01934 L
.97619 .01472 L
Mfstroke
0 0 1 r
.05556 .60347 m
.07143 .5628 L
.0873 .52923 L
.10317 .50076 L
.11905 .47602 L
.13492 .45412 L
.15079 .43445 L
.16667 .41658 L
.18254 .40023 L
.19841 .38517 L
.21429 .37121 L
.23016 .35822 L
.24603 .3461 L
.2619 .33474 L
.27778 .32408 L
.29365 .31407 L
.30952 .30464 L
.3254 .29578 L
.34127 .28743 L
.35714 .27959 L
.37302 .27223 L
.38889 .26534 L
.40476 .25892 L
.42063 .25295 L
.43651 .24745 L
.45238 .24243 L
.46825 .23791 L
.48413 .23392 L
.5 .23052 L
.51587 .22778 L
.53175 .22583 L
.54762 .2249 L
.56349 .22544 L
.57937 .22909 L
.59524 .23671 L
.61111 .241 L
.62698 .24264 L
.64286 .24222 L
.65873 .24019 L
.6746 .23691 L
.69048 .23265 L
.70635 .2276 L
.72222 .22194 L
.7381 .21579 L
.75397 .20926 L
.76984 .20242 L
.78571 .19533 L
.80159 .18806 L
.81746 .18065 L
.83333 .17312 L
Mistroke
.84921 .16551 L
.86508 .15785 L
.88095 .15015 L
.89683 .14243 L
.9127 .1347 L
.92857 .12697 L
.94444 .11926 L
.96032 .11157 L
.97619 .10391 L
Mfstroke
% End of Graphics
MathPictureEnd
\
\>"], "Graphics",
  ImageSize->{413, 255.125},
  ImageMargins->{{43, 0}, {0, 0}},
  ImageRegion->{{0, 1}, {0, 1}},
  ImageCache->GraphicsData["Bitmap", "\<\
CF5dJ6E]HGAYHf4PAg9QL6QYHg<PAVmbKF5d0`40006M0000ob000`400?l00000o`00003ooooooonN
ooooo`00oooooonNooooo`00oooooonNooooo`008_ooool2o`0003Gooooo0ol00004ooooo`;o0000
1?ooool2o`0002gooooo1Ol00003ooooo`;o00001?ooool2o`0002kooooo0ol00004ooooo`;o0000
1?ooool2o`0002oooooo0ol00003ooooo`;o00001?ooool2o`0002kooooo0ol00004ooooo`;o0000
1?ooool2o`0002kooooo0ol00004ooooo`;o00001?ooool2o`0000Cooooo000Qooooo`04o`000?oo
ooooooooo`0003Gooooo00?o0000ooooooooool00_ooool01?l0003oooooooooool00002ooooo`04
o`000?ooooooooooo`0002gooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool0
0002ooooo`04o`000?ooooooooooo`0002cooooo00Go0000ooooooooooooooooo`000002ooooo`04
o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000;oooool01Ol0003ooooooooooooo
oooo000000;ooooo00Co0000oooooooooooo00000_ooool00ol0003oooooooooo`0Zooooo`05o`00
0?ooooooooooooooool000000_ooool01?l0003oooooooooool00002ooooo`04o`000?oooooooooo
o`0002cooooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo
00Co0000oooooooooooo00000oooool0027ooooo00Co0000oooooooooooo0000=Oooool00ol0003o
ooooooooo`02ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000;_ooool0
0ol0003oooooooooo`02ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000
<?ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00?o0000oooo
ooooool0:_ooool5o`0000;ooooo00Co0000oooooooooooo00000_ooool01?l0003oooooooooool0
000`ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo00000_ooool00ol0003o
ooooooooo`0Zooooo`05o`000?ooooooooooooooool000000_ooool01?l0003oooooooooool00002
ooooo`04o`000?ooooooooooo`0000?ooooo000Qooooo`04o`000?ooooooooooo`0003Gooooo00?o
0000ooooooooool00_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0002oo
oooo00Go0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00?o0000
ooooooooool0;?ooool2o`0000?ooooo00Co0000oooooooooooo00000_ooool01?l0003ooooooooo
ool0000/ooooo`04o`000?ooooooooooo`0000?ooooo00Co0000oooooooooooo00000_ooool01?l0
003oooooooooool0000`ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000
0_ooool00ol0003oooooooooo`0Zooooo`Co00000oooool01?l0003oooooooooool00002ooooo`04
o`000?ooooooooooo`0000?ooooo000Qooooo`04o`000?ooooooooooo`0003Gooooo00?o0000oooo
ooooool00_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0002cooooo00Go
0000ooooooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooo
oooo0000<?ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000;ooooo00?o
0000ooooooooool0:oooool00ol0003oooooo`000003ooooo`04o`000?ooooooooooo`0000;ooooo
00Co0000oooooooooooo0000;?ooool4o`0000?ooooo00Co0000oooooooooooo00000_ooool01?l0
003oooooooooool0000/ooooo`03o`000?oooooooooo00Cooooo00Co0000oooooooooooo00000_oo
ool01?l0003oooooooooool00003ooooo`008Oooool01?l0003oooooooooool0000dooooo`;o0000
1?ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0002cooooo00Go0000oooo
ooooooooooooo`000002ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000oooooooooooo0000
;?ooool01Ol0003ooooooooooooooooo000000;ooooo00Co0000oooooooooooo00000_ooool01?l0
003oooooooooool0000^ooooo`;o00000oooool01?l0003oooooooooool00002ooooo`04o`000?oo
ooooooooo`0002cooooo00?o0000ooooooooool01?ooool01?l0003oooooooooool00002ooooo`04
o`000?ooooooooooo`0002gooooo00?o0000ooooooooool00oooool01?l0003oooooooooool00002
ooooo`04o`000?ooooooooooo`0000?ooooo000Rooooo`;o0000=_ooool00ol0003oooooooooo`03
ooooo`;o00001?ooool2o`0002kooooo0ol00004ooooo`;o00001?ooool2o`0002kooooo0ol00004
ooooo`;o00001?ooool2o`00033ooooo00?o0000ooooooooool00_ooool2o`0000Cooooo0_l0000]
ooooo`Go00000oooool2o`0000Cooooo0_l0000_ooooo`;o00001?ooool2o`0000Cooooo0_l00004
ooooo`00oooooonNooooo`00oooooonNooooo`00oooooonNooooo`00oooooonNooooo`00oooooonN
ooooo`006_ooooooo`0008Co0000000Jooooo`03o`000?oooooooooo00Kooooo00?o0000oooooooo
ool02_ooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`0:ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00[ooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`0:ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00[ooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02_ooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00[ooooo00?o0000oooooooo
ool02Oooool00ol0003oooooooooo`09ooooo`03o`000?oooooooooo00Wooooo00?o0000oooooooo
ool02?ooool001[ooooo00?o0000ooooooooool01_ooool00ol0003oooooooooo`0:ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00[ooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02_ooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00[ooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02_ooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`0:ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02_ooool00ol0003oooooooooo`09ooooo`03o`00
0?oooooooooo00Wooooo00?o0000ooooooooool02Oooool00ol0003oooooooooo`08ooooo`006_oo
ool00ol0003oooooooooo`06ooooo`03o`000?oooooooooo03_ooooo00?o0000ooooooooool0>_oo
ool00ol0003oooooooooo`0jooooo`03o`000?oooooooooo03[ooooo00?o0000ooooooooool0>_oo
ool00ol0003oooooooooo`0jooooo`03o`000?oooooooooo00Sooooo000Jooooo`03o`000?oooooo
oooo0?ooooooPOooool001[ooooo00?o0000ooooooooool0oooooon1ooooo`006_ooool00ol0003o
ooooooooo`3oooooogGooooo0_oo000:ooooo`006_ooool00ol0003oooooooooo`3oooooog;ooooo
0ooo000<ooooo`006_ooool00ol0003oooooooooo`3oooooofoooooo0ooo000?ooooo`006_ooool0
0ol0003oooooooooo`3oooooofcooooo0ooo000Booooo`006_ooool00ol0003oooooooooo`3ooooo
ofOooooo1Ooo000Eooooo`006_ooool00ol0003oooooooooo`3oooooof?ooooo1?oo000Jooooo`00
6_ooool00ol0003oooooooooo`3oooooof3ooooo0ooo000Nooooo`006_ooool00ol0003ooooooooo
o`3ooooooegooooo0ooo000Qooooo`006_ooool00ol0003oooooooooo`3ooooooe[ooooo0ooo000T
ooooo`006_ooool00ol0003oooooooooo`3ooooooeOooooo0ooo000Wooooo`006_ooool00ol0003o
ooooooooo`3ooooooeCooooo0ooo000Zooooo`006_ooool00ol0003oooooooooo`3ooooooe7ooooo
0ooo000]ooooo`006_ooool00ol0003oooooooooo`3oooooodkooooo0ooo000`ooooo`006_ooool0
0ol0003oooooooooo`3ooooood[ooooo1?oo000cooooo`006_ooool00ol0003oooooooooo`3ooooo
odOooooo0ooo000gooooo`006_ooool00ol0003oooooooooo`3oooooodCooooo0ooo000jooooo`00
6_ooool00ol0003oooooooooo`3ooooood7ooooo0ooo000mooooo`006_ooool00ol0003ooooooooo
o`3oooooockooooo0ooo0010ooooo`006_ooool00ol0003oooooooooo`3ooooooc_ooooo0ooo0013
ooooo`006_ooool00ol0003oooooooooo`3oooooocSooooo0ooo0016ooooo`006_ooool00ol0003o
ooooooooo`3oooooocGooooo0ooo0019ooooo`006_ooool00ol0003oooooooooo`3ooooooc;ooooo
0ooo001<ooooo`006_ooool00ol0003oooooooooo`3ooooooboooooo0ooo001?ooooo`003oooool3
o`0000Sooooo00?o0000ooooooooool0ooooool/ooooo`?oo`00D_ooool0013ooooo00?o0000oooo
ooooool01oooool00ol0003oooooooooo`3oooooobWooooo0ooo001Eooooo`004?ooool00ol0003o
ooooooooo`07ooooo`Co0000oooooolUooooo`?oo`00F?ooool0013ooooo00?o0000ooooooooool0
1oooool00ol0003oooooooooo`3oooooob?ooooo0ooo001Kooooo`004?ooool00ol0003ooooooooo
o`07ooooo`03o`000?oooooooooo0?oooooo8?ooool3ool005kooooo000?ooooo`;o00002Oooool0
0ol0003oooooooooo`3ooooooagooooo0ooo001Qooooo`004?ooool00ol0003oooooooooo`07oooo
o`03o`000?oooooooooo0?oooooo6_ooool3ool006Cooooo000Jooooo`03o`000?oooooooooo0?oo
oooo5oooool3ool006Oooooo000Jooooo`03o`000?oooooooooo0?oooooo4oooool4ool006[ooooo
000Jooooo`03o`000?oooooooooo0?oooooo4?ooool3ool006kooooo000Jooooo`03o`000?oooooo
oooo0?oooooo3Oooool3ool0077ooooo000Jooooo`03o`000?oooooooooo0?oooooo2oooool2ool0
06Sooooo0_l00?l:ooooo`006_ooool00ol0003oooooooooo`3oooooo`Wooooo0_oo001Xooooo`;o
003o3?ooool001[ooooo00?o0000ooooooooool0ooooool7ooooo`;oo`00J?ooool2o`00o`kooooo
000Jooooo`03o`000?oooooooooo0?oooooo1?ooool3ool006Sooooo0_l00?l@ooooo`006_ooool0
0ol0003oooooooooo`3oooooo`7ooooo0ooo001Yooooo`;o003o4_ooool001[ooooo00?o0000oooo
ooooool0oOooool3ool006[ooooo0_l00?lDooooo`006_ooool00ol0003oooooooooo`3jooooo`?o
o`00Joooool2o`00oaKooooo000Jooooo`03o`000?oooooooooo0?Sooooo0_oo001/ooooo`;o003o
6?ooool001[ooooo00?o0000ooooooooool0m_ooool2ool006cooooo0_l00?lJooooo`006_ooool0
0ol0003oooooooooo`3dooooo`;oo`00K?ooool2o`00oacooooo000Jooooo`03o`000?oooooooooo
0?7ooooo0ooo001/ooooo`;o003o7_ooool001[ooooo00?o0000ooooooooool0k_ooool3ool006go
oooo0_l00?lPooooo`006_ooool00ol0003oooooooooo`3[ooooo`?oo`00K_ooool2o`00ob;ooooo
000Jooooo`03o`000?oooooooooo0>Sooooo0ooo001_ooooo`;o003o9?ooool000kooooo1Ol00007
ooooo`03o`000?oooooooooo0>Kooooo0_oo001`ooooo`;o003o9_ooool000oooooo00Co0000oooo
oooooooo00001oooool00ol0003oooooooooo`3Tooooo`;oo`00L?ooool2o`00obSooooo000@oooo
o`03o`000?oooooooooo00Oooooo1?l0003Pooooo`?oo`00L?ooool2o`00ob[ooooo000Aooooo`03
o`000?oooooooooo00Kooooo00?o0000ooooooooool0g_ooool3ool0077ooooo0_l00?l/ooooo`00
3_ooool01Ol0003ooooooooooooooooo000000Oooooo00?o0000ooooooooool0foooool3ool007;o
oooo0_l00?l^ooooo`003_ooool01Ol0003ooooooooooooooooo000000Oooooo00?o0000oooooooo
ool0fOooool2ool007?ooooo0_l00?l`ooooo`003oooool3o`0000Sooooo00?o0000ooooooooool0
eoooool2ool007?ooooo0_l00?lbooooo`006_ooool00ol0003oooooooooo`3Eooooo`;oo`00Looo
ool2o`00ocCooooo000Jooooo`03o`000?oooooooooo0=?ooooo0_oo001booooo`?o003o=_ooool0
01[ooooo00?o0000ooooooooool0dOooool2ool007;ooooo0_l00?liooooo`006_ooool00ol0003o
ooooooooo`3?ooooo`;oo`00L_ooool2o`00oc_ooooo000Jooooo`03o`000?oooooooooo0<cooooo
0ooo001booooo`;o003o?Oooool001[ooooo00?o0000ooooooooool0bOooool3ool007?ooooo0_l0
0?loooooo`006_ooool00ol0003oooooooooo`37ooooo`;oo`00M?ooool2o`00od7ooooo000Joooo
o`03o`000?oooooooooo0<Gooooo0_oo001dooooo`;o003o@oooool001[ooooo00?o0000oooooooo
ool0`oooool2ool007Cooooo0_l00?m5ooooo`006_ooool00ol0003oooooooooo`31ooooo`;oo`00
M?ooool2o`00odOooooo000Jooooo`?o0000_oooool2ool007Cooooo0_l00?m9ooooo`006_ooool0
0ol0003oooooooooo`2mooooo`;oo`00M?ooool2o`00od_ooooo000Jooooo`03o`000?oooooooooo
0;_ooooo0_oo001dooooo`;o003oCOooool001[ooooo00?o0000ooooooooool0^Oooool2ool007Co
oooo0_l00?m?ooooo`006_ooool00ol0003oooooooooo`2gooooo`;oo`00M?ooool2o`00oe7ooooo
000Jooooo`03o`000?oooooooooo0;Gooooo0_oo001dooooo`;o003oDoooool001[ooooo00?o0000
ooooooooool0/oooool2ool007?ooooo0ol00?mEooooo`006_ooool00ol0003oooooooooo`2aoooo
o`;oo`00L_ooool3o`00oeSooooo000Jooooo`03o`000?oooooooooo0:oooooo0_oo001booooo`;o
003oFoooool001[ooooo00?o0000ooooooooool0[Oooool2ool007;ooooo0_l00?mMooooo`006_oo
ool00ol0003oooooooooo`2[ooooo`;oo`00L_ooool2o`00oeoooooo000Jooooo`?o0000Z?ooool3
ool0077ooooo0ol00?mQooooo`006_ooool00ol0003oooooooooo`2Vooooo`;oo`00LOooool3o`00
ofCooooo000Jooooo`03o`000?oooooooooo0:Cooooo0_oo001aooooo`;o003oIoooool001[ooooo
00?o0000ooooooooool0X_ooool2ool0077ooooo0_l00?mYooooo`006_ooool00ol0003ooooooooo
o`2Pooooo`;oo`00L?ooool3o`00of_ooooo000Jooooo`03o`000?oooooooooo09kooooo0_oo000W
oooooa;o003o=_ooool3o`00ofkooooo000?ooooo`?o00002?ooool00ol0003oooooooooo`2Loooo
o`;oo`008oooool6o`00oa;ooooo0ol00?l`ooooo`?o003oLOooool000kooooo00Go0000oooooooo
ooooooooo`000007ooooo`03o`000?oooooooooo09[ooooo0_oo000Oooooo`Ko003o6oooool3o`00
ob[ooooo0ol00?mdooooo`004_ooool00ol0003oooooooooo`05ooooo`Co0000Uoooool2ool001go
oooo1?l00?lTooooo`?o003o9?ooool3o`00ogOooooo000Booooo`03o`000?oooooooooo00Gooooo
00?o0000ooooooooool0Uoooool00ooo003oooooooooo`0Jooooo`?o003o:oooool4o`00oagooooo
0ol00?mjooooo`003_ooool4o`0000Sooooo00?o0000ooooooooool0UOooool2ool001[ooooo0ol0
0?lbooooo`?o003o5Oooool5o`00oggooooo000>ooooo`03o`000?oooooooooo00Wooooo00?o0000
ooooooooool0U?ooool00ooo003oooooooooo`0Gooooo`?o003o>?ooool4o`00o`_ooooo1_l00?n2
ooooo`003_ooool5o`0000Oooooo00?o0000ooooooooool0T_ooool2ool001Kooooo1?l00?looooo
o`_o003oR?ooool001[ooooo00?o0000ooooooooool0T?ooool2ool001Gooooo0ol00?oFooooo`00
6_ooool00ol0003oooooooooo`2>ooooo`;oo`005?ooool3o`00omWooooo000Jooooo`?o0000S?oo
ool2ool001?ooooo0ol00?oLooooo`006_ooool00ol0003oooooooooo`2;ooooo`03ool00?oooooo
oooo013ooooo0ol00?oOooooo`006_ooool00ol0003oooooooooo`29ooooo`;oo`004?ooool3o`00
on;ooooo000Jooooo`03o`000?oooooooooo08Sooooo00?oo`00ooooooooool03Oooool3o`00onGo
oooo000Jooooo`03o`000?oooooooooo08Kooooo0_oo000>ooooo`;o003oj?ooool001[ooooo00?o
0000ooooooooool0QOooool00ooo003oooooooooo`0<ooooo`;o003oj_ooool001[ooooo0ol00023
ooooo`;oo`003Oooool2o`00oncooooo000Jooooo`03o`000?oooooooooo08;ooooo00?oo`00oooo
ooooool02_ooool3o`00onkooooo000Jooooo`03o`000?oooooooooo083ooooo0_oo000:ooooo`?o
003olOooool001[ooooo00?o0000ooooooooool0Ooooool00ooo003oooooooooo`08ooooo`;o003o
m?ooool001[ooooo00?o0000ooooooooool0OOooool2ool000Wooooo0_l00?ofooooo`006_ooool3
o`0007cooooo00?oo`00ooooooooool01oooool2o`00ooSooooo000Jooooo`03o`000?oooooooooo
07[ooooo0_oo0008ooooo`;o003on_ooool001[ooooo00?o0000ooooooooool0NOooool00ooo003o
ooooooooo`06ooooo`;o003oo?ooool001[ooooo00?o0000ooooooooool0Moooool2ool000Oooooo
0_l00?onooooo`006_ooool3o`0007Kooooo00?oo`00ooooooooool01Oooool2o`00oooooooo0Ooo
ool001[ooooo00?o0000ooooooooool0M?ooool2ool000Kooooo0_l00?ooooooo`?ooooo0009oooo
o`?o00001?ooool2o`0000Sooooo00?o0000ooooooooool0L_ooool2ool000Kooooo0_l00?oooooo
o`Gooooo000:ooooo`03o`000?oooooooooo00;ooooo00Co0000oooooooooooo00001oooool00ol0
003oooooooooo`1aooooo`03ool00?oooooooooo00Gooooo00?o003oooooooooool0ooooool5oooo
o`002_ooool00ol0003oooooooooo`02ooooo`04o`000?ooooooooooo`0000Oooooo1?l0001^oooo
o`;oo`001_ooool2o`00oooooooo2?ooool000[ooooo00?o0000ooooooooool00_ooool01?l0003o
ooooooooool00007ooooo`03o`000?oooooooooo06gooooo0_oo0007ooooo`03o`00oooooooooooo
0?oooooo2?ooool000[ooooo00?o0000ooooooooool00_ooool01?l0003oooooooooool00007oooo
o`03o`000?oooooooooo06cooooo00?oo`00ooooooooool01Oooool2o`00oooooooo2oooool000Wo
oooo0_l00004ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000ooooooooool0J_ooool2ool0
00Gooooo0ol00?ooooooo`gooooo000:ooooo`03o`000?oooooooooo00?ooooo0_l00008ooooo`03
o`000?oooooooooo06Wooooo00?oo`00ooooooooool00oooool2o`00oooooooo4?ooool001[ooooo
00?o0000ooooooooool0Ioooool2ool000Cooooo0_l00?oooooooa;ooooo000Jooooo`03o`000?oo
oooooooo06Kooooo00?oo`00ooooooooool00oooool00ol00?ooooooooooo`3ooooooa;ooooo000J
ooooo`03o`000?oooooooooo06Gooooo00?oo`00ooooooooool00_ooool2o`00oooooooo5Oooool0
01[ooooo00?o0000ooooooooool0I?ooool00ooo003oooooooooo`02ooooo`03o`00oooooooooooo
0?oooooo5Oooool001[ooooo00?o0000ooooooooool0Hoooool01?oo003oooooooooooooool2o`00
oooooooo6?ooool001[ooooo00?o0000ooooooooool0HOooool2ool000?ooooo00?o003ooooooooo
ool0oooooolHooooo`006_ooool00ol0003oooooooooo`1Pooooo`04ool00?ooooooooooooooo`;o
003ooooooolKooooo`006_ooool00ol0003oooooooooo`1Oooooo`05ool00?ooooooooooooooool0
0?l0oooooolMooooo`006_ooool00ol0003oooooooooo`1Nooooo`03ool00?oooooooooo00;o003o
oooooolNooooo`006_ooool00ol0003oooooooooo`1Mooooo`04ool00?ooooooooooo`00oooooooo
8?ooool001[ooooo00?o0000ooooooooool0Foooool2ool000;ooooo00?o003oooooooooool0oooo
oolOooooo`006_ooool00ol0003oooooooooo`1Jooooo`05ool00?ooooooooooooooool00?l0oooo
oolRooooo`006_ooool00ol0003oooooooooo`1Iooooo`05ool00?ooooooooooooooool00?l0oooo
oolSooooo`006_ooool00ol0003oooooooooo`1Hooooo`03ool00?oooooooooo00;o003ooooooolT
ooooo`006_ooool00ol0003oooooooooo`1Gooooo`04ool00?ooooooooooo`00oooooooo9_ooool0
01[ooooo00?o0000ooooooooool0EOooool2ool00003ooooool00?oo003o0?oooooo9oooool001[o
oooo00?o0000ooooooooool0E?ooool01?oo003oooooooooool00?ooooooobWooooo000Jooooo`03
o`000?oooooooooo05?ooooo00Coo`00ooooool00?oo003ooooooolZooooo`006_ooool00ol0003o
ooooooooo`1Booooo`03ool00?oooooo003o0?oooooo;?ooool000Sooooo1Ol00003ooooo`;o0000
2?ooool00ol0003oooooooooo`1Aooooo`03ool00?oooooo003o0?oooooo;Oooool000Wooooo00Co
0000oooooooooooo00000_ooool01?l0003oooooooooool00007ooooo`03o`000?oooooooooo053o
oooo00?oo`00ooooool00?l0ooooool^ooooo`002_ooool00ol0003oooooooooo`02ooooo`04o`00
0?ooooooooooo`0000Oooooo1?l0001>ooooo`03ool00?oooooo003o0?oooooo;oooool000_ooooo
00Go0000ooooooooooooooooo`000002ooooo`03o`000?oooooooooo00Gooooo00?o0000oooooooo
ool0C_ooool00ooo003oooooo`00o`3ooooooc3ooooo0008ooooo`05o`000?ooooooooooooooool0
00000_ooool01?l0003oooooooooool00007ooooo`03o`000?oooooooooo04gooooo00?oo`00oooo
ool00?l0oooooolaooooo`002?ooool01Ol0003ooooooooooooooooo000000;ooooo00Co0000oooo
oooooooo00001oooool00ol0003oooooooooo`1<ooooo`03ool00?oooooo003o0?oooooo<_ooool0
00Wooooo0ol00004ooooo`;o00002?ooool00ol0003oooooooooo`1;ooooo`03ool00?oooooo003o
0?oooooo<oooool001[ooooo00?o0000ooooooooool0B_ooool00ooo003oooooo`00o`3oooooocCo
oooo000Jooooo`03o`000?oooooooooo04Wooooo00?oo`00ooooool00?l0ooooooleooooo`006_oo
ool00ol0003oooooooooo`18ooooo`03ool00?l00?oo003o0?oooooo=_ooool001[ooooo00?o0000
ooooooooool0Aoooool00ooo003o003oooooo`3oooooocOooooo000Jooooo`03o`000?oooooooooo
04Kooooo00?oo`00o`00ooooool0oooooolhooooo`006_ooool00ol0003oooooooooo`15ooooo`03
ool00?l00?oooooo0?oooooo>Oooool001[ooooo00?o0000ooooooooool0A?ooool00ooo003o003o
ooooo`3ooooooc[ooooo000Jooooo`03o`000?oooooooooo04Cooooo00?o003oooooooooool0oooo
ooljooooo`006_ooool00ol0003oooooooooo`13ooooo`03ool00?l00?oooooo0?oooooo>oooool0
01[ooooo00?o0000ooooooooool0@_ooool00ooo003o003oooooo`3ooooooccooooo000Jooooo`?o
0000@Oooool00ooo003o003oooooo`3oooooocgooooo000Jooooo`03o`000?oooooooooo043ooooo
00?oo`00o`00ooooool0oooooolnooooo`006_ooool00ol0003oooooooooo`0oooooo`03ool00?l0
0?oooooo0?oooooo?oooool001[ooooo00?o0000ooooooooool0?_ooool00ooo003o003oooooo`3o
oooood3ooooo000Jooooo`03o`000?oooooooooo03gooooo00?oo`00o`00ooooool0oooooom1oooo
o`006_ooool00ol0003oooooooooo`0looooo`03ool00?l00?oooooo0?oooooo@_ooool001[ooooo
00?o0000ooooooooool0>oooool00ooo003o003oooooo`3ooooood?ooooo000Jooooo`03o`000?oo
oooooooo03[ooooo00?oo`00o`00ooooool0oooooom4ooooo`006_ooool00ol0003oooooooooo`0i
ooooo`03ool00?l00?oooooo0?ooooooAOooool001[ooooo00?o0000ooooooooool0>?ooool00ooo
003o003oooooo`3oooooodKooooo000Jooooo`03o`000?oooooooooo03Oooooo00?oo`00o`00oooo
ool0oooooom7ooooo`006_ooool3o`0003Oooooo00?o003oooooooooool0oooooom7ooooo`006_oo
ool00ol0003oooooooooo`0fooooo`03ool00?l00?oooooo0?ooooooB?ooool001[ooooo00?o0000
ooooooooool0=Oooool00ooo003o003oooooo`3oooooodWooooo000Jooooo`03o`000?oooooooooo
03Cooooo00?oo`00o`00ooooool0oooooom:ooooo`006_ooool00ol0003oooooooooo`0dooooo`03
o`00oooooooooooo0?ooooooB_ooool001[ooooo00?o0000ooooooooool0<oooool00ol00?oooooo
ooooo`3ooooood_ooooo0009ooooo`?o00001?ooool2o`0000Sooooo00?o0000ooooooooool0<_oo
ool00ol00?ooooooooooo`3oooooodcooooo0008ooooo`05o`000?ooooooooooooooool000000_oo
ool01?l0003oooooooooool00007ooooo`03o`000?oooooooooo037ooooo00?o003oooooooooool0
oooooom=ooooo`003?ooool01?l0003oooooooooool00002ooooo`03o`000?oooooooooo00Gooooo
1?l0000`ooooo`03o`00oooooooooooo0?ooooooCOooool000cooooo00Co0000oooooooooooo0000
0_ooool00ol0003oooooooooo`05ooooo`03o`000?oooooooooo033ooooo00?o003oooooooooool0
oooooom>ooooo`002?ooool4o`0000?ooooo00Co0000oooooooooooo00001oooool00ol0003ooooo
ooooo`0_ooooo`03o`00oooooooooooo0?ooooooCoooool000Sooooo00?o0000ooooooooool01?oo
ool01?l0003oooooooooool00007ooooo`03o`000?oooooooooo02kooooo00?oo`00o`00ooooool0
oooooom@ooooo`002?ooool5o`0000?ooooo0_l00008ooooo`03o`000?oooooooooo02kooooo00?o
003oooooooooool0oooooom@ooooo`006_ooool00ol0003oooooooooo`0]ooooo`03o`00oooooooo
oooo0?ooooooDOooool001[ooooo00?o0000ooooooooool0;?ooool00ooo003o003oooooo`3ooooo
oe;ooooo000Jooooo`?o0000:oooool00ooo003o003oooooo`3ooooooe?ooooo000Jooooo`03o`00
0?oooooooooo02_ooooo00?o003oooooooooool0oooooomCooooo`006_ooool00ol0003ooooooooo
o`0Zooooo`03ool00?l00?oooooo0?ooooooE?ooool001[ooooo00?o0000ooooooooool0:Oooool0
0ooo003o003oooooo`3ooooooeGooooo000Jooooo`03o`000?oooooooooo02Wooooo00?o003ooooo
ooooool0oooooomEooooo`006_ooool3o`0002Sooooo00?oo`00o`00ooooool0oooooomFooooo`00
6_ooool00ol0003oooooooooo`0Xooooo`03o`00oooooooooooo0?ooooooE_ooool001[ooooo00?o
0000ooooooooool09oooool00ol00?ooooooooooo`3ooooooeOooooo000Jooooo`03o`000?oooooo
oooo02Kooooo00?oo`00o`00ooooool0oooooomHooooo`006_ooool00ol0003oooooooooo`0Voooo
o`03o`00oooooooooooo0?ooooooF?ooool001[ooooo0ol0000Uooooo`03o`00oooooooooooo0?oo
ooooFOooool001[ooooo00?o0000ooooooooool09Oooool00ol00?ooooooooooo`3ooooooeWooooo
000Jooooo`03o`000?oooooooooo02Cooooo00?o003oooooooooool0oooooomJooooo`006_ooool0
0ol0003oooooooooo`0Tooooo`03o`00oooooooooooo0?ooooooF_ooool001[ooooo00?o0000oooo
ooooool08oooool00ol00?ooooooooooo`3ooooooe_ooooo000Jooooo`?o00008oooool00ol00?oo
ooooooooo`3ooooooe_ooooo000Jooooo`03o`000?oooooooooo02;ooooo00?o003oooooooooool0
oooooomLooooo`000oooool3o`0000Cooooo0_l00004ooooo`;o00002?ooool00ol0003ooooooooo
o`0Qooooo`03o`00oooooooooooo0?ooooooGOooool000Cooooo00?o0000ooooooooool00_ooool0
1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000ooooooooool0
8Oooool00ol00?ooooooooooo`3ooooooegooooo0004ooooo`03o`000?oooooooooo00;ooooo00Co
0000oooooooooooo00000_ooool01?l0003oooooooooool00007ooooo`Co00007oooool00ol00?oo
ooooooooo`3ooooooekooooo0004ooooo`03o`000?oooooooooo00;ooooo00Co0000oooooooooooo
00000_ooool01?l0003oooooooooool00007ooooo`03o`000?oooooooooo023ooooo00?o003ooooo
ooooool0oooooomNooooo`001?ooool00ol0003oooooooooo`02ooooo`04o`000?ooooooooooo`00
00;ooooo00Co0000oooooooooooo00001oooool00ol0003oooooooooo`0Oooooo`03o`00oooooooo
oooo0?ooooooGoooool000?ooooo0_l00004ooooo`04o`000?ooooooooooo`0000;ooooo00Co0000
oooooooooooo00001oooool00ol0003oooooooooo`0Oooooo`03o`00oooooooooooo0?ooooooGooo
ool000Cooooo00?o0000ooooooooool00oooool2o`0000Cooooo0_l00008ooooo`03o`000?oooooo
oooo01kooooo00?o003oooooooooool0oooooomPooooo`006_ooool00ol0003oooooooooo`0Noooo
o`03o`00oooooooooooo0?ooooooH?ooool001[ooooo00?o0000ooooooooool07Oooool00ol00?oo
ooooooooo`3oooooof7ooooo000Jooooo`03o`000?oooooooooo01gooooo00?o003oooooooooool0
oooooomQooooo`006_ooool00ol0003oooooooooo`0Looooo`03o`00oooooooooooo0?ooooooH_oo
ool001[ooooo00?o0000ooooooooool07?ooool00ol00?ooooooooooo`3oooooof;ooooo000Joooo
o`03o`000?oooooooooo01cooooo00?o003oooooooooool0oooooomRooooo`006_ooool00ol0003o
ooooooooo`0Kooooo`03o`00oooooooooooo0?ooooooHoooool001[ooooo00?o0000ooooooooool0
6oooool00ol00?ooooooooooo`3oooooof?ooooo000Jooooo`03o`000?oooooooooo01[ooooo00?o
003oooooooooool0oooooomTooooo`006_ooool00ol0003oooooooooo`0Jooooo`03o`00oooooooo
oooo0?ooooooI?ooool001[ooooo00?o0000ooooooooool06Oooool00ol00?ooooooooooo`3ooooo
ofGooooo000Jooooo`03o`000?oooooooooo01Wooooo00?o003oooooooooool0oooooomUooooo`00
6_ooool00ol0003oooooooooo`0Hooooo`03o`00oooooooooooo0?ooooooI_ooool001[ooooo00?o
0000ooooooooool06?ooool00ol00?ooooooooooo`3oooooofKooooo000Jooooo`03o`000?oooooo
oooo01Sooooo00?o003oooooooooool0oooooomVooooo`006_ooool00ol0003oooooooooo`0Goooo
o`03o`00oooooooooooo0?ooooooIoooool001[ooooo00?o0000ooooooooool05oooool00ol00?oo
ooooooooo`3oooooofOooooo000Jooooo`03o`000?oooooooooo01Oooooo00?o003oooooooooool0
oooooomWooooo`000_ooool5o`0000?ooooo0_l00004ooooo`;o00002?ooool00ol0003ooooooooo
o`0Fooooo`03o`00oooooooooooo0?ooooooJ?ooool000?ooooo00Co0000oooooooooooo00000_oo
ool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000oooooooo
ool05_ooool00ol00?ooooooooooo`3oooooofSooooo0004ooooo`03o`000?oooooooooo00;ooooo
00Co0000oooooooooooo00000_ooool01?l0003oooooooooool00007ooooo`Co00005?ooool00ol0
0?ooooooooooo`3oooooofWooooo0005ooooo`05o`000?ooooooooooooooool000000_ooool01?l0
003oooooooooool00002ooooo`03o`000?oooooooooo00Gooooo00?o0000ooooooooool05Oooool0
0ol00?ooooooooooo`3oooooofWooooo0002ooooo`05o`000?ooooooooooooooool000000_ooool0
1?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000ooooooooool0
5Oooool00ol00?ooooooooooo`3oooooofWooooo0002ooooo`05o`000?ooooooooooooooool00000
0_ooool01?l0003oooooooooool00002ooooo`04o`000?ooooooooooo`0000Oooooo00?o0000oooo
ooooool05?ooool00ol00?ooooooooooo`3oooooof[ooooo0003ooooo`?o00001?ooool2o`0000Co
oooo0_l00008ooooo`03o`000?oooooooooo01Cooooo00?o003oooooooooool0oooooomZooooo`00
6_ooool00ol0003oooooooooo`0Dooooo`03o`00oooooooooooo0?ooooooJ_ooool001[ooooo00?o
0000ooooooooool04oooool00ol00?ooooooooooo`3oooooof_ooooo000Jooooo`03o`000?oooooo
oooo0?ooooooPOooool001[ooooo00?o0000ooooooooool0oooooon1ooooo`006_ooool00ol0003o
ooooooooo`3ooooooh7ooooo000Jooooo`03o`000?oooooooooo0?ooooooPOooool001[ooooo00?o
0000ooooooooool0oooooon1ooooo`00oooooonNooooo`00\
\>"],
  ImageRangeCache->{{{0, 412}, {254.125, 0}} -> {-58.8613, -0.530814, \
1.63559, 0.0117063}}]
}, Open  ]],

Cell["\<\
Cross section (pb) as a function of the Higgs mass. This plot shows \
how well the EFT does. For a Higgs<200 GeV the approximation is very \
good.\
\>", "Text"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Various ways to do the integration over the x1 and x2.", "Subsection"],

Cell[CellGroupData[{

Cell["\<\

uno=Integrate[x1*x2,{x1,t0,1},{x2,t0/x1,1}]//Expand
due=Integrate[ta,{ta,t0,1},{y,Log[Sqrt[ta]],-Log[Sqrt[ta]]}]//PowerExpand//\
Expand
tre=Integrate[t0/z \
t0/z^2,{z,t0,1},{y,Log[Sqrt[t0/z]],-Log[Sqrt[t0/z]]}]//PowerExpand//Expand\
\>\
", "Input"],

Cell[BoxData[
    \(TraditionalForm\`1\/2\ \(log(t0)\)\ t0\^2 - t0\^2\/4 + 1\/4\)], "Output"],

Cell[BoxData[
    \(TraditionalForm\`1\/2\ \(log(t0)\)\ t0\^2 - t0\^2\/4 + 1\/4\)], "Output"],

Cell[BoxData[
    \(TraditionalForm\`1\/2\ \(log(t0)\)\ t0\^2 - t0\^2\/4 + 1\/4\)], "Output"]
}, Open  ]]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{

Cell["NLO cross section", "Section"],

Cell["\<\
We have integrate over the angular variables, so we are only left \
with two integrations. One is over z (=mh/S/x1/x2) and the other over the \
rapidity y of the partonic cms. For every point in the phase space we have to \
calculate an event and a corresponding counter-event with z=1 to implement \
the + distributions in the gg channel. Various contributions add to the final \
result:

virtual: gg>h  at 1-loop + corrections to the effective lagrangian (UV and IR \
divergent)
real: qq~ >h g (finite)
           qg > qh (collinear divergent)
            gg>gh  (soft and collinear divergent)
\
\>", "Text",
  FontSize->16],

Cell[CellGroupData[{

Cell["Born+Virtual ", "Subsection"],

Cell[TextData[{
  "dsigmaBV[yy_,mh_,sqrtS_,muf_,mur_]:=Module[\n",
  StyleBox["(* local variables *)",
    FontColor->RGBColor[0, 1, 0]],
  "\n{y0,x10,x20,sig0,ymax0,JAC0,v,S,tau0,beta,gg0,qg0,qq0},\nMuf=muf 1.;\n\
Mur=mur 1.;\nv=246.;\nS=sqrtS^2;\ntau0=mh^2/S;\nbeta=Sqrt[1-tau0];\n\n",
  StyleBox["(* calculate quantities for z=1 *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nymax0=-Log[Sqrt[tau0]];\ny0=-ymax0+2*ymax0*yy;\nJAC0=2*ymax0;\n\
x10=Sqrt[tau0] Exp[y0];\nx20=Sqrt[tau0] Exp[-y0];\n\
{gg0,qg0,qq0}=pdfcall[x10,x20,Muf];\n\n",
  StyleBox["(* sigma0 *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nsig0=asNLO[Mur,5]^2/576/Pi/v^2*tau0;\nsig0=sig0+sig0*asNLO[Mur,5]/2/Pi*\n\
(11/3 CA+ 2 Pi^2 - 2 b0 2 Log[Muf/Mur]+\n16 CA Log[beta] Log[mh/Muf]+16 CA \
Log[beta]^2);\nsig0=sig0*gg0;\nsig0=sig0*389379660; (*to picobarns*)\n\
sig0=sig0*JAC0;\n\n\
(*Print[{muf,mur,v,S,tau0,beta,ymax0,y0,JAC0,x10,x20,gg0,qg0,qq0,sig0}];*)\n\
Return[sig0];\n];"
}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Real contributions", "Subsection"],

Cell[TextData[{
  "dsigmaR[xx_,yy_,mh_,sqrtS_,muf_,mur_]:=Module[\n",
  StyleBox["(* local variables *)",
    FontColor->RGBColor[0, 1, 0]],
  "\n{v,S,tau0,\n y,y0,z,tau,\n x1,x2,x10,x20,\n ymax0,ymax,\n JAC,JAC0,\n \
gg,qg,qq,gg0,qg0,qq0,\n qqterm,qgterm,ggterm,ggterm0,\n sig,sig0,\n Muf},\n\n\
v=246.;\nS=sqrtS^2;\ntau0=mh^2/S;\nMuf=muf*1.0;\n\n",
  StyleBox["(* calculate quantities for an event *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nz=tau0+(1-tau0)*xx;\ntau=tau0/z;\nymax =-Log[Sqrt[tau]];\n\
y=-ymax+2*ymax*yy;\nJAC =2*ymax*(1-tau0)*tau0/z^2;\nx1=Sqrt[tau] Exp[y];\n\
x2=tau/x1;\n",
  StyleBox["(* call the pdf *)\n",
    FontColor->RGBColor[1, 0, 0]],
  "\n{gg,qg,qq}=pdfcall[x1,x2,Muf];\n\n",
  StyleBox["(* calculate quantities for counter-event *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nymax0=-Log[Sqrt[tau0]];\ny0=-ymax0+2*ymax0*yy;\n\
JAC0=2*ymax0*(1-tau0)*tau0;\nx10=Sqrt[tau0] Exp[y0];\nx20=tau0/x10;\n",
  StyleBox["(* call the pdf at z=1 *)",
    FontColor->RGBColor[1, 0, 0]],
  "\n{gg0,qg0,qq0}=pdfcall[x10,x20,Muf];\n\n",
  StyleBox["(* sigma0 *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nsig0=asNLO[mur,5]^2/576/Pi/v^2;\nsig0=sig0*asNLO[mur,5]/2/Pi;\n\n",
  StyleBox["(* qq channnel : no counter event *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nqqterm=64/27*(1-z)^3;\nqqterm=qqterm*sig0*JAC*qq;\n\n",
  StyleBox["(* qg channnel : no counter event *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nqgterm=CF*( (1+(1-z)^2)/z (2*Log[mh/muf]+2 Log[1-z]-Log[z])\n           \
+(z^2-3/2(1-z)^2)/z )*z;\nqgterm=qgterm*sig0*JAC*qg;\n\n",
  StyleBox["(* gg channnel *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nggterm=CA*(2 (2 (z/(1-z)+(1-z)/z+z (1-z) )) * (2*Log[mh/muf])-\n       \
11/3 (1-z)^3/z -\n       4 (1-z+z^2)^2/z/(1-z) Log[z]+\n       8 \
(1-z+z^2)^2/z Log[1-z]/(1-z) )*z;\nggterm=ggterm*sig0*JAC*gg;\n",
  StyleBox["(* gg counter-event *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nggterm0=CA*(-4/(1-z) 2*Log[mh/muf] - 8*Log[1-z]/(1-z) );\n\
ggterm0=ggterm0*sig0*JAC0*gg0;\n\n\n",
  StyleBox["(* total *)",
    FontColor->RGBColor[1, 0, 0]],
  "\nsig=0;\nsig=sig+qqterm;\nsig=sig+qgterm;\nsig=sig+ggterm+ggterm0;\n\
sig=sig*389379660; (*to picobarns*)\n\nReturn[sig];\n\n];\n"
}], "Input"],

Cell[CellGroupData[{

Cell["\<\
virt=NIntegrate[dsigmaBV[xvar,100,14000,100,100],{xvar,0,1}];
virt\
\>", "Input"],

Cell[BoxData[
    \(TraditionalForm\`33.56877602575538`\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
eps=0.0000000001;
real=Vegas[dsigmaR[xvar,yvar,100,14000,100,100],{xvar,eps,1-eps},{yvar,eps,1-\
eps},Compiled->False,NStart->1000,MaxPoints->10000]\
\>", "Input"],

Cell[BoxData[
    \(TraditionalForm\`"\nIteration 1:  1000 integrand evaluations so far\n\
[1] 16.6638 +- 0.729445  \tchisq -0 (0 df)"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`"\nIteration 2:  2500 integrand evaluations so far\n\
[1] 16.6521 +- 0.173781  \tchisq 0.000273976 (1 df)"\)], "Print"],

Cell[BoxData[
    \(TraditionalForm\`$Aborted\)], "Output"]
}, Open  ]],

Cell[CellGroupData[{

Cell["asNLO[100,5]", "Input"],

Cell[BoxData[
    \(TraditionalForm\`0.11636176382460801`\)], "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
FrontEndVersion->"4.1 for Macintosh",
ScreenRectangle->{{0, 1280}, {0, 832}},
WindowSize->{872, 573},
WindowMargins->{{131, Automatic}, {106, Automatic}}
]

(*******************************************************************
Cached data follows.  If you edit this Notebook file directly, not
using Mathematica, you must remove the line containing CacheID at
the top of  the file.  The cache data will then be recreated when
you save this file from within Mathematica.
*******************************************************************)

(*CellTagsOutline
CellTagsIndex->{}
*)

(*CellTagsIndex
CellTagsIndex->{}
*)

(*NotebookFileOutline
Notebook[{

Cell[CellGroupData[{
Cell[1727, 52, 56, 0, 65, "Subtitle"],

Cell[CellGroupData[{
Cell[1808, 56, 31, 0, 56, "Section"],
Cell[1842, 58, 687, 14, 114, "Text"]
}, Closed]],

Cell[CellGroupData[{
Cell[2566, 77, 32, 0, 36, "Section"],
Cell[2601, 79, 88, 4, 57, "Input"],

Cell[CellGroupData[{
Cell[2714, 87, 35, 0, 46, "Subsection"],
Cell[2752, 89, 159, 4, 32, "Text"],

Cell[CellGroupData[{
Cell[2936, 97, 161, 2, 43, "Input"],
Cell[3100, 101, 201, 3, 28, "Output"]
}, Open  ]]
}, Closed]],

Cell[CellGroupData[{
Cell[3350, 110, 52, 0, 30, "Subsection"],
Cell[3405, 112, 68, 3, 42, "Input"]
}, Closed]],

Cell[CellGroupData[{
Cell[3510, 120, 73, 0, 30, "Subsection"],

Cell[CellGroupData[{
Cell[3608, 124, 155, 3, 43, "Input"],
Cell[3766, 129, 113, 2, 28, "Output"],
Cell[3882, 133, 80, 1, 28, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[3999, 139, 44, 3, 42, "Input"],
Cell[4046, 144, 72, 1, 24, "Print"],
Cell[4121, 147, 68, 1, 24, "Print"],
Cell[4192, 150, 80, 1, 24, "Print"],
Cell[4275, 153, 88, 1, 24, "Print"],
Cell[4366, 156, 97, 2, 24, "Print"],
Cell[4466, 160, 122, 2, 24, "Print"],
Cell[4591, 164, 53, 1, 24, "Print"],
Cell[4647, 167, 81, 1, 24, "Print"],
Cell[4731, 170, 92, 1, 24, "Print"]
}, Open  ]],

Cell[CellGroupData[{
Cell[4860, 176, 349, 6, 58, "Subsubsection"],
Cell[5212, 184, 945, 37, 522, "Input"],
Cell[6160, 223, 37, 3, 42, "Input"],
Cell[6200, 228, 1739, 64, 897, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[7976, 297, 58, 0, 42, "Subsubsection"],
Cell[8037, 299, 43, 0, 27, "Input"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[8129, 305, 267, 6, 82, "Subsection"],
Cell[8399, 313, 778, 16, 155, "Input"],

Cell[CellGroupData[{
Cell[9202, 333, 76, 2, 42, "Subsubsection"],

Cell[CellGroupData[{
Cell[9303, 339, 32, 0, 27, "Input"],
Cell[9338, 341, 71, 1, 28, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[9470, 349, 36, 0, 46, "Subsection"],
Cell[9509, 351, 76, 7, 102, "Input"]
}, Closed]]
}, Open  ]],

Cell[CellGroupData[{
Cell[9634, 364, 35, 0, 56, "Section"],

Cell[CellGroupData[{
Cell[9694, 368, 207, 5, 58, "Subsubsection"],
Cell[9904, 375, 18, 0, 27, "Input"],

Cell[CellGroupData[{
Cell[9947, 379, 236, 6, 72, "Input"],
Cell[10186, 387, 28784, 687, 280, 5861, 400, "GraphicsData", "PostScript", \
"Graphics"]
}, Open  ]],
Cell[38985, 1077, 150, 3, 32, "Text"]
}, Open  ]],

Cell[CellGroupData[{
Cell[39172, 1085, 468, 9, 74, "Subsubsection"],
Cell[39643, 1096, 420, 20, 297, "Input"],

Cell[CellGroupData[{
Cell[40088, 1120, 72, 0, 27, "Input"],
Cell[40163, 1122, 70, 1, 28, "Output"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{
Cell[40282, 1129, 215, 5, 64, "Subsection"],
Cell[40500, 1136, 113, 3, 27, "Input"],
Cell[40616, 1141, 139, 3, 42, "Input"],
Cell[40758, 1146, 24, 0, 27, "Input"],

Cell[CellGroupData[{
Cell[40807, 1150, 319, 12, 72, "Input"],
Cell[41129, 1164, 34897, 892, 264, 10632, 588, "GraphicsData", "PostScript", \
"Graphics"]
}, Open  ]],
Cell[76041, 2059, 170, 4, 30, "Text"]
}, Open  ]],

Cell[CellGroupData[{
Cell[76248, 2068, 76, 0, 46, "Subsection"],

Cell[CellGroupData[{
Cell[76349, 2072, 259, 8, 72, "Input"],
Cell[76611, 2082, 93, 1, 46, "Output"],
Cell[76707, 2085, 93, 1, 46, "Output"],
Cell[76803, 2088, 93, 1, 46, "Output"]
}, Open  ]]
}, Closed]]
}, Closed]],

Cell[CellGroupData[{
Cell[76957, 2096, 36, 0, 36, "Section"],
Cell[76996, 2098, 636, 15, 194, "Text"],

Cell[CellGroupData[{
Cell[77657, 2117, 35, 0, 46, "Subsection"],
Cell[77695, 2119, 958, 19, 462, "Input"]
}, Open  ]],

Cell[CellGroupData[{
Cell[78690, 2143, 40, 0, 46, "Subsection"],
Cell[78733, 2145, 2219, 45, 1107, "Input"],

Cell[CellGroupData[{
Cell[80977, 2194, 91, 3, 42, "Input"],
Cell[81071, 2199, 69, 1, 28, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[81177, 2205, 173, 4, 57, "Input"],
Cell[81353, 2211, 147, 2, 60, "Print"],
Cell[81503, 2215, 156, 2, 60, "Print"],
Cell[81662, 2219, 59, 1, 28, "Output"]
}, Open  ]],

Cell[CellGroupData[{
Cell[81758, 2225, 29, 0, 27, "Input"],
Cell[81790, 2227, 71, 1, 28, "Output"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)



(*******************************************************************
End of Mathematica Notebook file.
*******************************************************************)

