highT.nb

(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 12.2' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    122319,       2714]
NotebookOptionsPosition[    116423,       2611]
NotebookOutlinePosition[    116828,       2627]
CellTagsIndexPosition[    116785,       2624]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Phase Transition under High-T Expansion", "Title",ExpressionUUID->"6cd35c81-39c0-48fd-830b-7d456a3fdc89"],

Cell[BoxData[
 RowBox[{"ClearAll", "[", "\"\<Global`*\>\"", "]"}]], "Input",
 CellLabel->
  "In[826]:=",ExpressionUUID->"33b22c75-f5aa-47cf-a757-aba03fa7960e"],

Cell[CellGroupData[{

Cell["Analytical Calculation", "Section",ExpressionUUID->"09784b56-5c1a-48b1-ae64-9874125f0427"],

Cell[CellGroupData[{

Cell["Expression for critical temperature", "Subsection",ExpressionUUID->"eeecd131-c144-47e2-b47a-2c75fcfa74c5"],

Cell["\<\
This section calculates the formula for critical temperature.\
\>", "TextIndent",ExpressionUUID->"e6d25092-ca41-45b6-8ad4-525823e653cd"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"$Assumptions", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "Dterm", ",", "h", ",", "Eterm", ",", "T", ",", "\[Lambda]term", ",", 
       "T0"}], "}"}], "\[Element]", "PositiveReals"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Vana", "=", 
    RowBox[{
     RowBox[{"Dterm", " ", 
      RowBox[{"(", 
       RowBox[{
        SuperscriptBox["T", "2"], "-", 
        SuperscriptBox["T0", "2"]}], ")"}], 
      SuperscriptBox["h", "2"]}], "-", 
     RowBox[{"Eterm", " ", "T", " ", 
      SuperscriptBox["h", "3"]}], "+", 
     RowBox[{
      FractionBox["1", "4"], "\[Lambda]term", " ", 
      SuperscriptBox["h", "4"]}]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"General", " ", "form", " ", "of", " ", "high"}], "-", 
    RowBox[{"T", " ", "potential"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"hmin", "=", 
   RowBox[{
    RowBox[{"h", "/.", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{
        RowBox[{"Solve", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"D", "[", 
            RowBox[{"Vana", ",", "h"}], "]"}], "\[Equal]", "0"}], ",", "h"}], 
         "]"}], "//", "Normal"}], ")"}], "[", 
      RowBox[{"[", "2", "]"}], "]"}]}], "//", "Simplify"}]}], 
  RowBox[{"(*", 
   RowBox[{"Non", "-", 
    RowBox[{"zero", " ", "local", " ", "minimum"}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Tcana", "=", 
   RowBox[{
    RowBox[{
     RowBox[{"T", "/.", 
      RowBox[{"Solve", "[", 
       RowBox[{
        RowBox[{
         RowBox[{"(", 
          RowBox[{
           RowBox[{"Vana", "/.", 
            RowBox[{"{", 
             RowBox[{"h", "\[Rule]", "hmin"}], "}"}]}], "//", "Simplify"}], 
          ")"}], "\[Equal]", "0"}], ",", "T"}], "]"}]}], "//", "Normal"}], "//",
     "Simplify"}]}], 
  RowBox[{"(*", 
   RowBox[{"Critical", " ", "temperature"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"hoverT", "=", 
  RowBox[{
   RowBox[{
    FractionBox["hmin", "Tc"], "/.", 
    RowBox[{"{", 
     RowBox[{"T", "\[Rule]", "Tc"}], "}"}]}], "//", "Simplify"}]}]}], "Input",\

 CellLabel->
  "In[541]:=",ExpressionUUID->"69d3d85d-3aa9-4f5b-8164-a8603914bbae"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"3", " ", "Eterm", " ", "T"}], "+", 
   SqrtBox[
    RowBox[{
     RowBox[{"9", " ", 
      SuperscriptBox["Eterm", "2"], " ", 
      SuperscriptBox["T", "2"]}], "+", 
     RowBox[{"8", " ", "Dterm", " ", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"-", 
         SuperscriptBox["T", "2"]}], "+", 
        SuperscriptBox["T0", "2"]}], ")"}], " ", "\[Lambda]term"}]}]]}], 
  RowBox[{"2", " ", "\[Lambda]term"}]]], "Output",
 CellLabel->
  "Out[543]=",ExpressionUUID->"f615bf16-5f90-47aa-87b7-d7fadc145000"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"T0", " ", 
   SqrtBox[
    FractionBox[
     RowBox[{"Dterm", " ", "\[Lambda]term"}], 
     RowBox[{
      RowBox[{"-", 
       SuperscriptBox["Eterm", "2"]}], "+", 
      RowBox[{"Dterm", " ", "\[Lambda]term"}]}]]]}], "}"}]], "Output",
 CellLabel->
  "Out[544]=",ExpressionUUID->"a238d36f-4ce5-4cef-a3e9-3239616a2953"],

Cell[BoxData[
 FractionBox[
  RowBox[{
   RowBox[{"3", " ", "Eterm", " ", "Tc"}], "+", 
   SqrtBox[
    RowBox[{
     RowBox[{"9", " ", 
      SuperscriptBox["Eterm", "2"], " ", 
      SuperscriptBox["Tc", "2"]}], "+", 
     RowBox[{"8", " ", "Dterm", " ", 
      RowBox[{"(", 
       RowBox[{
        SuperscriptBox["T0", "2"], "-", 
        SuperscriptBox["Tc", "2"]}], ")"}], " ", "\[Lambda]term"}]}]]}], 
  RowBox[{"2", " ", "Tc", " ", "\[Lambda]term"}]]], "Output",
 CellLabel->
  "Out[545]=",ExpressionUUID->"1da25cff-277e-4fbd-8bc6-8a6dec03a7e9"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Expression for eliminating ",
 Cell[BoxData[
  FormBox["S", TraditionalForm]], "None",ExpressionUUID->
  "bfc58de3-3664-429d-a682-b5708f98020e"]
}], "Subsection",ExpressionUUID->"9c08ab25-f758-456f-83bb-d617fea16fa1"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{
  RowBox[{"$Assumptions", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
      "Dterm", ",", "h", ",", "Eterm", ",", "T", ",", "\[Lambda]term", ",", 
       "T0", ",", "\[Mu]H", ",", "A", ",", "\[Mu]S", ",", "cs", ",", "S", ",",
        "v"}], "}"}], "\[Element]", "PositiveReals"}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"VwithS", "=", 
   RowBox[{
    RowBox[{"Dterm", " ", 
     SuperscriptBox["T", "2"], 
     SuperscriptBox["h", "2"]}], "-", 
    RowBox[{
     FractionBox["1", "2"], 
     SuperscriptBox["\[Mu]H", "2"], 
     SuperscriptBox["h", "2"]}], "-", 
    RowBox[{"Eterm", " ", "T", " ", 
     SuperscriptBox["h", "3"]}], "+", 
    RowBox[{
     FractionBox["1", "4"], "\[Lambda]term", " ", 
     SuperscriptBox["h", "4"]}], "+", 
    RowBox[{
     FractionBox["1", "2"], 
     SuperscriptBox["\[Mu]S", "2"], 
     SuperscriptBox["S", "2"]}], "-", 
    RowBox[{
     FractionBox["1", "2"], "A", 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["h", "2"], "+", 
       RowBox[{"cs", " ", 
        SuperscriptBox["T", "2"]}], "-", 
       RowBox[{"2", 
        SuperscriptBox["v", "2"]}]}], ")"}], "S"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Svev", "=", 
  RowBox[{
   RowBox[{
    RowBox[{"S", "/.", 
     RowBox[{"Solve", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"D", "[", 
         RowBox[{"VwithS", ",", "S"}], "]"}], "\[Equal]", "0"}], ",", "S"}], 
      "]"}]}], "//", "Normal"}], "//", 
   "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"Dprime", "=", 
  RowBox[{
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"VwithS", "/.", 
      RowBox[{"{", 
       RowBox[{"S", "->", "Svev"}], "}"}]}], ",", "h", ",", "2"}], "]"}], "//",
    "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"Eprime", "=", 
  RowBox[{
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"VwithS", "/.", 
      RowBox[{"{", 
       RowBox[{"S", "->", "Svev"}], "}"}]}], ",", "h", ",", "3"}], "]"}], "//",
    "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{"\[Lambda]prime", "=", 
  RowBox[{
   RowBox[{"Coefficient", "[", 
    RowBox[{
     RowBox[{"VwithS", "/.", 
      RowBox[{"{", 
       RowBox[{"S", "->", "Svev"}], "}"}]}], ",", "h", ",", "4"}], "]"}], "//",
    "Simplify"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"VwithS", "/.", 
   RowBox[{"{", 
    RowBox[{"S", "->", "Svev"}], "}"}]}], "//", "FullSimplify"}]}], "Input",
 CellLabel->
  "In[546]:=",ExpressionUUID->"64989a1f-b188-48e1-ae1b-76b6648eea81"],

Cell[BoxData[
 RowBox[{"{", 
  FractionBox[
   RowBox[{"A", " ", 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["h", "2"], "+", 
      RowBox[{"cs", " ", 
       SuperscriptBox["T", "2"]}], "-", 
      RowBox[{"2", " ", 
       SuperscriptBox["v", "2"]}]}], ")"}]}], 
   RowBox[{"2", " ", 
    SuperscriptBox["\[Mu]S", "2"]}]], "}"}]], "Output",
 CellLabel->
  "Out[548]=",ExpressionUUID->"d2ce8b13-ffa6-4811-b9f0-3dea23d3a418"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["1", "4"], " ", 
   RowBox[{"(", 
    RowBox[{
     RowBox[{"4", " ", "Dterm", " ", 
      SuperscriptBox["T", "2"]}], "-", 
     RowBox[{"2", " ", 
      SuperscriptBox["\[Mu]H", "2"]}], "+", 
     FractionBox[
      RowBox[{
       SuperscriptBox["A", "2"], " ", 
       RowBox[{"(", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "cs"}], " ", 
          SuperscriptBox["T", "2"]}], "+", 
         RowBox[{"2", " ", 
          SuperscriptBox["v", "2"]}]}], ")"}]}], 
      SuperscriptBox["\[Mu]S", "2"]]}], ")"}]}], "}"}]], "Output",
 CellLabel->
  "Out[549]=",ExpressionUUID->"395fa0e5-48ae-4841-a4cd-91bbecadee05"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "Eterm"}], " ", "T"}], "}"}]], "Output",
 CellLabel->
  "Out[550]=",ExpressionUUID->"38497e21-7299-4ba5-8c4b-6665c41f4c95"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   FractionBox["\[Lambda]term", "4"], "-", 
   FractionBox[
    SuperscriptBox["A", "2"], 
    RowBox[{"8", " ", 
     SuperscriptBox["\[Mu]S", "2"]}]]}], "}"}]], "Output",
 CellLabel->
  "Out[551]=",ExpressionUUID->"7beb5017-34e2-405b-9181-8b6cb65c0cdc"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{
    FractionBox["1", "4"], " ", 
    SuperscriptBox["h", "2"], " ", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       RowBox[{"-", "4"}], " ", "Eterm", " ", "h", " ", "T"}], "+", 
      RowBox[{"4", " ", "Dterm", " ", 
       SuperscriptBox["T", "2"]}], "+", 
      RowBox[{
       SuperscriptBox["h", "2"], " ", "\[Lambda]term"}], "-", 
      RowBox[{"2", " ", 
       SuperscriptBox["\[Mu]H", "2"]}]}], ")"}]}], "-", 
   FractionBox[
    RowBox[{
     SuperscriptBox["A", "2"], " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{
        SuperscriptBox["h", "2"], "+", 
        RowBox[{"cs", " ", 
         SuperscriptBox["T", "2"]}], "-", 
        RowBox[{"2", " ", 
         SuperscriptBox["v", "2"]}]}], ")"}], "2"]}], 
    RowBox[{"8", " ", 
     SuperscriptBox["\[Mu]S", "2"]}]]}], "}"}]], "Output",
 CellLabel->
  "Out[552]=",ExpressionUUID->"4ab905f3-58b1-46ed-8624-8243c274cfc7"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(", 
   RowBox[{
    RowBox[{"Dprime", " ", 
     SuperscriptBox["h", "2"]}], "+", 
    RowBox[{"Eprime", " ", 
     SuperscriptBox["h", "3"]}], "+", 
    RowBox[{"\[Lambda]prime", " ", 
     SuperscriptBox["h", "4"]}], "-", 
    RowBox[{"(", 
     RowBox[{"VwithS", "/.", 
      RowBox[{"{", 
       RowBox[{"S", "->", "Svev"}], "}"}]}], ")"}]}], ")"}], "//", 
  "FullSimplify", 
  RowBox[{"(*", 
   RowBox[{
   "Check", " ", "whether", " ", "we", " ", "have", " ", "some", " ", 
    "remaining", " ", "terms"}], "*)"}]}]], "Input",
 CellLabel->
  "In[553]:=",ExpressionUUID->"6bfa97e4-7db2-4078-9334-574123f43310"],

Cell[BoxData[
 RowBox[{"{", 
  FractionBox[
   RowBox[{
    SuperscriptBox["A", "2"], " ", 
    SuperscriptBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{"cs", " ", 
        SuperscriptBox["T", "2"]}], "-", 
       RowBox[{"2", " ", 
        SuperscriptBox["v", "2"]}]}], ")"}], "2"]}], 
   RowBox[{"8", " ", 
    SuperscriptBox["\[Mu]S", "2"]}]], "}"}]], "Output",
 CellLabel->
  "Out[553]=",ExpressionUUID->"dbbca74e-433b-4ee5-8080-dd2298de43fb"]
}, Open  ]],

Cell["This is field-independent, and can be ignored.", "TextIndent",ExpressionUUID->"39cec182-a215-4bf2-81e2-baa4eb199a6e"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Define ",
 Cell[BoxData[
  FormBox[
   SubscriptBox["V", "eff"], TraditionalForm]], "None",ExpressionUUID->
  "aee47986-e905-4a97-a7c2-b2da209cc48e"],
 " in high-T with physical parameters"
}], "Section",ExpressionUUID->"9661d2b9-84dc-431f-850a-75d73e467719"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"g", "=", "0.65"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"gY", "=", "0.36"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"yt", "=", "0.9945"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"DSM", "=", 
   RowBox[{
    FractionBox["1", "16"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"3", 
       SuperscriptBox["g", "2"]}], "+", 
      SuperscriptBox["gY", "2"], "+", 
      RowBox[{"4", 
       SuperscriptBox["yt", "2"]}]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ESM", "=", 
   RowBox[{
    FractionBox["1", 
     RowBox[{"48", "\[Pi]"}]], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{"2", 
       SuperscriptBox["g", "3"]}], "+", 
      SuperscriptBox[
       RowBox[{"(", 
        RowBox[{
         SuperscriptBox["g", "2"], "+", 
         SuperscriptBox["gY", "2"]}], ")"}], 
       RowBox[{"3", "/", "2"}]]}], ")"}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cs", "=", 
   FractionBox["1", "3"]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"mh", "=", "125.13"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"GF", " ", "=", " ", 
   RowBox[{"1.16637", " ", "*", 
    SuperscriptBox["10", 
     RowBox[{"-", "5"}]]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"v", "=", 
   FractionBox["1", 
    SqrtBox[
     RowBox[{"2", 
      SqrtBox["2"], "*", "GF"}]]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Mu]H2", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   RowBox[{
    FractionBox["1", "2"], 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       SuperscriptBox["mh", "2"], 
       RowBox[{"(", 
        RowBox[{"1", "-", 
         SuperscriptBox["sin\[Theta]", "2"]}], ")"}]}], "+", 
      RowBox[{
       SuperscriptBox["mS", "2"], 
       SuperscriptBox["sin\[Theta]", "2"]}]}], ")"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Mu]S2", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{
     SuperscriptBox["sin\[Theta]", "2"], 
     SuperscriptBox["mh", "2"]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{"1", "-", 
       SuperscriptBox["sin\[Theta]", "2"]}], ")"}], 
     SuperscriptBox["mS", "2"]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"A", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   FractionBox[
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["mh", "2"], "-", 
       SuperscriptBox["mS", "2"]}], ")"}], "sin\[Theta]", 
     SqrtBox[
      RowBox[{"(", 
       RowBox[{"1", "-", 
        SuperscriptBox["sin\[Theta]", "2"]}], ")"}]]}], 
    RowBox[{
     SqrtBox["2"], "v"}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Lambda]", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   FractionBox[
    RowBox[{
     RowBox[{
      RowBox[{"(", 
       RowBox[{"1", "-", 
        SuperscriptBox["sin\[Theta]", "2"]}], ")"}], 
      SuperscriptBox["mh", "2"]}], "+", 
     RowBox[{
      SuperscriptBox["sin\[Theta]", "2"], 
      SuperscriptBox["mS", "2"]}]}], 
    RowBox[{"4", 
     SuperscriptBox["v", "2"]}]]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Spath", "[", 
   RowBox[{"mS_", ",", "sin\[Theta]_", ",", "h_", ",", "T_"}], "]"}], ":=", 
  FractionBox[
   RowBox[{
    RowBox[{"A", " ", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
    RowBox[{"(", 
     RowBox[{
      SuperscriptBox["h", "2"], "+", 
      RowBox[{"cs", " ", 
       SuperscriptBox["T", "2"]}], "-", 
      RowBox[{"2", " ", 
       SuperscriptBox["v", "2"]}]}], ")"}]}], 
   RowBox[{"2", " ", 
    RowBox[{"\[Mu]S2", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"VhighT2d", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_", ",", "h_", ",", "S_", ",", "T_"}], 
    "]"}], "=", 
   RowBox[{
    RowBox[{"DSM", " ", 
     SuperscriptBox["T", "2"], 
     SuperscriptBox["h", "2"]}], "-", 
    RowBox[{
     FractionBox["1", "2"], 
     RowBox[{"\[Mu]H2", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
     SuperscriptBox["h", "2"]}], "-", 
    RowBox[{"ESM", " ", "T", " ", 
     SuperscriptBox["h", "3"]}], "+", 
    RowBox[{
     FractionBox["1", "4"], 
     RowBox[{"\[Lambda]", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
     SuperscriptBox["h", "4"]}], "+", 
    RowBox[{
     FractionBox["1", "2"], 
     RowBox[{"\[Mu]S2", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
     SuperscriptBox["S", "2"]}], "-", 
    RowBox[{
     FractionBox["1", "2"], 
     RowBox[{"A", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
     RowBox[{"(", 
      RowBox[{
       SuperscriptBox["h", "2"], "+", 
       RowBox[{"cs", " ", 
        SuperscriptBox["T", "2"]}], "-", 
       RowBox[{"2", 
        SuperscriptBox["v", "2"]}]}], ")"}], "S"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Deff", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   RowBox[{"DSM", "-", " ", 
    FractionBox[
     RowBox[{"cs", " ", 
      SuperscriptBox[
       RowBox[{"A", "[", 
        RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "2"]}], 
     RowBox[{"4", 
      RowBox[{"\[Mu]S2", "[", 
       RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Lambda]eff", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"\[Lambda]", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "-", 
    FractionBox[
     SuperscriptBox[
      RowBox[{"A", "[", 
       RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "2"], 
     RowBox[{"2", " ", 
      RowBox[{"\[Mu]S2", "[", 
       RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"T0", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], ":=", 
   SqrtBox[
    FractionBox[
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        FractionBox["1", "2"], 
        RowBox[{"\[Mu]H2", "[", 
         RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}], "-", 
       FractionBox[
        RowBox[{
         SuperscriptBox["v", "2"], 
         SuperscriptBox[
          RowBox[{"A", "[", 
           RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "2"]}], 
        RowBox[{"2", 
         RowBox[{"\[Mu]S2", "[", 
          RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}], ")"}], 
     RowBox[{"(", 
      RowBox[{"DSM", "-", " ", 
       FractionBox[
        RowBox[{"cs", " ", 
         SuperscriptBox[
          RowBox[{"A", "[", 
           RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "2"]}], 
        RowBox[{"4", 
         RowBox[{"\[Mu]S2", "[", 
          RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}], ")"}]]]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Tc", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"T0", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
    SqrtBox[
     FractionBox[
      RowBox[{
       RowBox[{"Deff", "[", 
        RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
       RowBox[{"\[Lambda]eff", "[", 
        RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}], 
      RowBox[{
       RowBox[{"-", 
        SuperscriptBox["ESM", "2"]}], "+", 
       RowBox[{
        RowBox[{"Deff", "[", 
         RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
        RowBox[{"\[Lambda]eff", "[", 
         RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]}]]]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"strength", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], "=", 
   FractionBox[
    RowBox[{"2", "ESM"}], 
    RowBox[{"\[Lambda]eff", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"VhighT", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_", ",", "h_", ",", "T_"}], "]"}], "=", 
   RowBox[{
    RowBox[{
     RowBox[{"Deff", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
     SuperscriptBox["T", "2"], 
     SuperscriptBox["h", "2"]}], "+", 
    RowBox[{
     RowBox[{"(", 
      RowBox[{
       RowBox[{
        RowBox[{"-", 
         FractionBox["1", "2"]}], 
        RowBox[{"\[Mu]H2", "[", 
         RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}], "+", 
       FractionBox[
        RowBox[{" ", 
         RowBox[{
          SuperscriptBox["v", "2"], 
          SuperscriptBox[
           RowBox[{"A", "[", 
            RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "2"]}]}], 
        RowBox[{"2", 
         RowBox[{"\[Mu]S2", "[", 
          RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]]}], ")"}], 
     SuperscriptBox["h", "2"]}], "-", 
    RowBox[{"ESM", " ", "T", " ", 
     SuperscriptBox["h", "3"]}], "+", 
    RowBox[{
     FractionBox["1", "4"], 
     RowBox[{"\[Lambda]eff", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], 
     SuperscriptBox["h", "4"]}]}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"Positive\[ReversePrime]factor", "[", 
    RowBox[{"mS_", ",", "sin\[Theta]_"}], "]"}], "=", 
   RowBox[{
    RowBox[{"-", 
     SuperscriptBox["ESM", "2"]}], "+", 
    RowBox[{
     RowBox[{"Deff", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
     RowBox[{"\[Lambda]eff", "[", 
      RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]}]}], ";"}]}], "Input",
 CellLabel->
  "In[827]:=",ExpressionUUID->"0a8b36df-21ed-4130-9d65-d0a68edb84ce"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"sin", "[", 
    RowBox[{"mS_", ",", "ft_"}], "]"}], ":=", 
   SqrtBox[
    RowBox[{
     FractionBox[
      SuperscriptBox["mS", "2"], 
      RowBox[{
       SuperscriptBox["mh", "2"], "-", 
       SuperscriptBox["mS", "2"]}]], 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"1", "/", "ft"}], "-", "1"}], ")"}]}]]}], ";"}]], "Input",
 CellLabel->
  "In[850]:=",ExpressionUUID->"5499e233-fd83-4c7d-9b85-367412aebf6a"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Plot", "[", 
  RowBox[{
   RowBox[{"Tc", "[", 
    RowBox[{
     SuperscriptBox["10", "logms"], ",", 
     RowBox[{"sin", "[", 
      RowBox[{
       SuperscriptBox["10", "logms"], ",", "0.1"}], "]"}]}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"logms", ",", 
     RowBox[{"-", "1.5"}], ",", "1"}], "}"}], ",", 
   RowBox[{"PlotRange", "\[Rule]", "Full"}]}], "]"}]], "Input",
 CellLabel->"In[32]:=",ExpressionUUID->"3db30cd6-516d-49ae-b089-bcfff54becc1"],

Cell[BoxData[
 GraphicsBox[{{{}, {}, 
    TagBox[
     {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[4], 
      AbsoluteThickness[1.6], Opacity[1.], LineBox[CompressedData["
1:eJwd13c8lV8YAHCaVhnZUVYUISEyOhRKRmlJJSVSRhlFvwZJSUV2pTKipKjM
ouEpZGRzyXbteV17VPid8/51P9/P+95znvM85znn80raXjxgv4SFhaWGlYWF
/N7oMxtfXJwBDmaph1mWvh6b0u9138tnQGX/JcoRl06Z3Hw2AzFRQpSbhF6d
PnVrBpi6fe7EEl8GvPTOz8BjnyzKyaweL1k0ZqA2zYyy4EPFj/lrZ+Aig4/y
LdH+In+WGdhwuN6N+Kiq9RDXr2m4dsGScj4IzVe+n4ZwVgHKSqY1q8PDp6HR
r8qVeLn97q0i1tPwfo0+5bOqj63e7pqGWo6/F4mLWftuastPQ9dcOuXAGP9K
65kpYPcTo7ym4adzfOgUfPOmuRBfvRxnknFlChrVvSl38l2TL7CZguErGyin
m23p71GcAu3jzs7UfM+jLS9ITkHJBR3KokMcRdP8U7C5jpPyrYCeVyv/TcKT
kCQn4gP5T89s+jUJ1641OhJHrcjfn/B1EqRMYinTjYd0xT9MQtA2O8ouVdrC
vBGTkHpg6DxxQGtL2aw1ttbIOeJvM+u2FY9NgP7F/LPEC828Spt6JoBb3Zuy
3vdlsvcbJuD+vAblgoBBAVOYAI2AJHviMpGPU5WBE3Br7Q074tULSYwtNyeg
5qEy5f2dz3pCPSbgVC/9DDHtrW/dwWMTILRVj3KLtlnWb7kJCCwcPU2cxv9X
Xl50AmbCwij7M5LirnNNwPuBrZSVY5cFSo6Nw/N1F09R+Vny1dbp8ziwLdJO
Eg+PhEVHZIzDBmYs5SPN5xu+pYyDg6cjZYVMIXOe2HGQ2b1gTcVj77E9y28c
LJskKMuVyPMsmo7DHLvFcSrf57s19xiNw4KEIGVv9pjToWgc8n82HSOuN+bN
kN46DmOZtpT9f80c3iM0Dj68LlbEPaUFz0I7xsB41tGSOKHipJyM5xhI/lx2
kNgUbtTRL46Bb03uAeLpD8/9os+PwYnTVygbhza1C1qPQdfigAUx8+CRJ+wG
Y1AoW7SfWLvRjJ3JOwZ1e8+ZE3eXOH9K5hyDryyilB9+fmB/bvkYcG4qNSPu
eFbyo2N2FPLLFCjfPWl4ldY+ChsNek2Ia7t0hnLejYIEOzKm4hP58uLQ21Fw
ecpNudB8+1Fm4ijsm2/fQ5ydo/5TJm4UWDRuUo4JVooNDhuFzZdzdxM7aksc
tPcahe4JFSNi1vClX3j0R+FT8sROql+Kfd2SdfH8a3IpT84vyhlpjYK4UQDl
3nP/Iq6pjkLBKjHKZTumLvRuGAXPdfr6xE8Ge6W+cozCViFfRHzydoRAZRIT
2IfqtYkr0PTKK/FMuH8mhjL6a/lH4jkTHJLsKUu4r213D2aC1aMJLapfbeKT
BD2ZsHKBi/JZ7VTtk7uYcLZGXZP4wnjpaUbrCCirnVSj9sOppe8514xAjOOC
AnFq5e6HNK4ReFaeTblzR+CF6BUjoDPuTtlIXEBZ+Q8Drvj0ylP90iSbeqCD
AXkJvzZR+Tu4Nz3qAwN+dfjLEYNR6Ec5UwbYR5VLEY9n1T0aNWSAyncfyjIb
RL1yEAMMXqhQDliaoGGiyoCNoRGSVP99z8y+IMoAfklLCep82N7wOat/GN6V
VIgTsyiuh113hsHUw1+Y+Lqg/sk73sPg9keO8syi7Xyh1zB83lAiRDxSk6i9
13EYrBw4KDdfUcrev28YmgMDBIg/FuimWYsMwx8eDz7q/DlxIsHr3RC4Oi1w
Uuf9g6i7yXWDwP2iYd6U7LefEimnKgch2/gc5bLF11UCJYNwK3DmH3Gux0eR
m18HoayNn3L88drkQwmD8CBl7x9iR4VVVfOug1AZFT9N/PfXTeH9XINwiZ+V
STy5jE13+YpBsDTyGyEe2RF8+vPiAKxbsZwyPT36rczEACQprmQQFzz9rDPX
OADHeVcMEQc6Tp6Kfz0AHyoYvcRiHOfeTO4cANho20Z8MrcsUEpnAFKdS1qJ
X7iruO5XH4DLolsoyzb/2ZYiNwDsaX+biZVTggrOcA2A+OcHjcQ7zTPaa+r6
ofzTQxrx+bAFgbRz/bBaQuwXcYrRmbm20/3wy+NiCRX/n6IWruP9kCGeV0zs
fiY04ZxZP3i52RcRX1eTUZHY2g8mK14WEIfUG5sG/+2D9H0sQPxRNPKWy8M+
aAzdkkH8+VbLwS93+yB0xDWdyveg9AZ23z7YF5aaRlz4Ob3opUcfMJIVU4nr
j1VxtRztA37+9e+Ip55yPt4r3QdbhNsSidXW3kqWy+6FgdKiKOL0tRdpHR09
cOvYuhvUc9pahYTmHrg6p3WdOCuw2Neurgeg9Mg14ux5SeW+4h7YtTLwPyqe
dtq94Q89wFY4fJm4NEELzdzogZxv/heIuxVWvOUS7QFhMyEb4rPd6SwVa3rg
eZfwSeK+5zaWwat6QDxSxJp4cFXOMj7WHsg1FD5OPMp0shHq74a2Qg5L4vmM
an6pj93gFJJvTiyoE+2jcaAb+Cy9dYk7LshFupl0w64gYx2qPi/S3iYbdAMj
co02Vc+VhXXrNbphWOSlJvGFmhEFNrFusCr7qEr885xeQ0NPF4jXvthEfCmy
a8vV/7qgiXGHn7hmZFPn17hO6L6xo9cEW9Rd5Ibj007YWFrRTXxmik1YOKIT
DEqsu4gn//aZedztBJtaLzqxIHtizqYLnXB4/lkT8XFpqdBH2p1w9HFCBXG3
paiea30HcEi8ziLe3MzeLF7VAQvSnJnEl0/OXS4t6YBbqS7pxCvsG5Jlv3XA
O2ulD8Qb3R8JtSZ0AHdXdBKxSyAf09itAzTPqT0jnv7OESPN1QHa99V9iANy
f88YL++AgD36N6j1fH1p4bpAh2xhk2vEOz7tWP5tlA78fUe9iP3fuTsfqaOD
xArri8QCT5u07sXQweTJchtiNY+3v0e20CFvvkSHuNDVS0VAng5WkhFaxEcv
GDzQlqbDosAJTeLr59t2BAjQYZljjyo1PsiY7ppsh0SfVnniV08m+5P820HM
47QwcefeCKlLyW2QFJ40thf73ruJNJaYNtiorMYk3sJzUD8opA10zn8bJvau
4z2VeLkNUGpxH7HoqeDoBtQGNbvzWoktLt8X1qW1QnwgKiH+Hntz9Yr5FnBw
soghjp50/vPYvBlkLu7UI9abfRTBubMZjldb6hJ3/v2u5KPeDIYhjlrEsksF
zpwVa4Zg0QA14ve830rVBpsgTjxhIzW+0qro6jtN0PDlCA/1/3PvEOe3Rrhb
G99qjH3b+Xejd2ojNLLYNBHLurJeGk9ohE1awr+JHb0OJzXeb4Rlhb5VxBN3
5nmSjjbCoZkt+cTLEsw6DSYb4KT1/GtiuVbGbW+FBvgbpXWR+ERdzfMu8QZw
nd3lRBxS/ilzD08DvPcwdiD+8823m2/qNzim7bEhrogRMHid+xt6W4X2EXue
0l1aZfEbX3tTisQFXYG+klfqodYwrH8P9lyzW5S/Yz1UhAh1EyvRjqQNnagH
Lq8n7cRPCiQ6svTrQXRZUD2x06tMvb2c9dCeqFdAzOfQsugeUwdCqhUxxKcH
FbwLCmiwMt52P3GYX3McTzYN+ouSTIjzxB7kn0imwVzUoBGx1L5BtqlQGvze
elqXuDMjKXyDDQ2+F4nIE9vd2PDGf64WDk9kLiE+y7OetkepFmLUPVN3k/xu
45Mvf1wNmU+V5o2wxTXubrR6UA13+U7MEFdp/JPt8a6G8/tvjxFrbO+Vnrev
hnucP3uIl+l8FldUq4bxPysqiGP0bXmDqqrgX2BANHG1afqsKXsVHOI+pUnM
9sx8JG2uElBbiwrxxeeaQyX0SnhPO6hArBPL1TP3vhKYUlvWEde/zGywMqkE
2bQXS4nZU5d9F7ldAbldZ8oNsV0LXz6MmioDm3U7ThD/S112N6KpDIzKjhwi
Dnhm5xMMZfC4yMGUOM5VxvXOvTLwDHDXIa4SfWnhJl4Gx23kxYiVLybwGxuW
QmGyWJMBud+E4p/ORpQA3Ye2j9j5fGyipVohsM1+ldhFnkdGZEqsLYTetMNC
xG559/IGWAtBmt6zithz7aW2a5U/IeXq8J+d2D7lewTjHH/Cj7JQGnGYyvid
/rgC2Dyp7k+cM2dg/9+qfDDck9ypT/Ihq+WxczIPjk1LNRDnHlD25WzOA2m+
0HLivGSR6OikPGB13Z9NXGY9QvuxKw/uVx4LIqb/eGzAcfUHpFk+UidexV+3
e1s5wPvkc956pD9Q8rGHjwBS/8W4Eac4+rr02gD8yyy2Izb9oRj+aDwXfDZN
7yUOvODfNi2UC/wFnwSJOUo0Ln+y/Qrj8b+SEamv99N4zdls2H2vtnQHOf/e
uGaF/MiGe/nW34iTaEbF/fezwRAa3xPvVphgPBHH7+cnhhL7N5hsn9v5CRpd
w48QL1edr8wJyoINNW9bdbGX9NvMa0mng4SmI10H+2SaiU8YIw3ej56oIv58
VYN18FMa2ErrfSd2X8W9PMo0DQ5++R1L3LU1l2vmcioIrdp7kjjgr6bm2dB3
oFii2KCNzSWXfvlISyKUc/uBFrb8vz0dLt1hMGb29agmtrlKfKSOWBhEir83
JPY4+9eY81AotH97uJXY68Vxe4u7waC6VYyL+OneDWxFhvfgsdrZXA1Sv+1u
z+zu/AeThUnriTff0Dcf8b+CNp3wpKljP0tquZ7x+hqqlJIGYg6aV/KVYm/E
//7rG+JLApa3We/cQnLlWd7EFs5BHQ8KA1D6rlw54hWr5s2qXoah2jBhTzVy
/n4+z9YjF4t4HDyWqmJrxmzR8WGJQxKLIQNbsXlsNk9yusShJMPISuLvdGm7
DUYvEO2qxTPi9d18Bkdn45Hyr0gV4rbB0aW5JxLR22oNKxVyn8+m+N3b8A7Z
/ZWPUMZOvhemH3I8C7Xn5qUrkPz12rtrPcxCv5WGw4mTdm5P6P6ehS57zV0i
Tvzbvmy77EeUqwDbiONdFIs7Rj8iFZPJbHmSP4tic1X/bLTReeTTJuz7IgvH
61O/oEJn+1g57CyamvbqpT/Qh6cyq6VJP2ased619gfySVTsl8LeGD7+L1vt
B3qazZ8nRY2Xmnvm7A8k4nn9MvHNCoWdOSU/kNS/902S2HuLpYztQvJQNfud
GAnSP194LD+LF6B+HkHOddgiz5ifgrcVILaueLo49oGrFcL2+wrQ0mz2j8Q/
NQObuH0L0I1Ou1PEb7NWnrTvLkDWV80yxLAvf1iw40n+if4eLrFYS+qbMOx+
VrMI+e9xdxXGXhd1uOHowSLko8nUId4anKtrcqEIGUlqsRNbXQth2/KyCMWI
bn8hRPJ5UC1mjrsYNe+vrxAk/bTseklgXzEaNreVECD3qz2nRMajX+hFzLFn
vNiLcvLli5PlyGjzZDc7tljLvJjM6gokfy37NbFWSLXzHrkKtGK/oSPxpdkr
XKFWFcgp79MIG/ZgcZGJVG4FeoH+TK7Epp23/2VwtxJ1phVPL8euO1zVVRBY
hdS9RCtYsbuVP3NaRlWhypj794jH2V+qDryqQrNJXw2IeXK9/FZBFfL2s//M
QvIvu176yFgVqhIviF/MxPtx2sW273A1moq+5TCP/fYxZwf7+hq0N9kBZrFl
b2g9eLu5BjH4NjsTx9ueVzfRqkG6OTHCxE8Vi+4FHq5Bh24Guc5g38+/tZU7
sAb9FQsQm8Y+z5zz45+rQTnv3thNYG/c3S+7vrYWbTIyz2Bgv9wsVA30WnS4
RvMwsQSf0bVTI7XIcKpwehhbuDWh8gU7DWlElGkSs3nYXJHRo6H9qh7Zg9h9
MfUlm97R0NmPLa/6sBOnC5zV7tYhq0q/bR3YctqbuEbC65BbtmYtHTvJJyj5
dVwdehIedpH4DduRQdHPdYherfq6HTtFtM+BlVGHBLSk+dqw03awn6k4UI+6
hZ3qm7C/+ptZnVv3Gxn93C5Iw64RrDd8ntWA1ss8OfkTu2riErtgXgMK2a/Z
XYBdXrWmPLiiAeVddTxPXHzf4pBfXwP6puPnno8NC2W2jiKNiCX0ms8PMn9v
vo/mjUbUJK18+xu2/8e0HJpBExKwKlPPwr4dvv+GuUUTiqnqScnE9nVl6hVb
NyErzdvSxNflFYu/eDahYHFz7gxsj5ik+vikJuSwmqMjFfuUf+yEK1czyssB
5xTs7UeCFFfRmlH942nJBFKPHW4ekR3NaGp19bV47NVyh3PEmc3I0Ze97gV2
94yYoRJHCzpAD70Thx0SlXJin14LmhNhtkWTerSUPghNaUF3xc3dn2BHnuEY
FLjTigxi+ISCsMdc7yau29aO3Mq2bL1E4smO64Cd7ajK+KKzB8nfYo7Y6X3t
KOX33Ct37NSHw2EJ59rRjchRQTfsCykWPhuftaOA4AOjLtjTE6s2i7DQUYb2
kZsOZL3KDcuk7elIruQtzzFs+Vq7eA8nOmL7IaVqhd3pOYYK3OiIJl996Cj2
gVyOa/bedGTkWRF5BHurme7Y68d0ZJq9bfVB7FHH+NbNv+ioQ1S0wQTbJdH5
4zblDsTrO/hLB7ulcm3k6m0d6JrexDdt7L1zvzx6dTqQUen4By3SP6abVB7t
7UCdm2tCNLHpYz3JU/YdiMV0ylAN+5CuzYvM5x0ocaHzugK2Ds0icCtHJ7rt
XScgQvpxnsWJg7cTbTr5pkcIW0Qu1bhTqBN5WllkCpL1/7d6ZdiGTmS1TMSc
n+RrfanvmF4nYoad9eDGlnYy8Er16kRlnON+y7E5lmjYKvV0oqcb92SOZujr
ffJaerB3qBMNeUbYMrHtGJW7YsY7kT0zinsEO7fh3IbVLF3ogEup/RC224en
/QzRLiT049WSXuyGEwsX3u3vQutGxOabsBM//by2+WsXUr85P5SPvcvlUKR8
eDfaWnfHLhTbPewLy8WobhTU7VcbjB3/Sco5I7Ybqa4S1n+IzbJkbKdOSjdS
D5kXekDmfxQ0avazG1mfjcq8ja35vdDEbaYb8ZaOvvbC3rxm+5Ls4z3oXnW+
ywns45pxLv9O9yDeLdMpx7DvW69s1DvXg96Fhw8exR54Xfeh5FIPYs1iP32Y
xKvjbt0c1IOMZL5omWNLnk3OXoAe9Ch0OBZhC+aIXzSU6UXthXxyktj/zrxe
y6vQiwqenFixHrtztUpxi0ovar0c2iWG/d7OQPIS6kVNmzMfC2Mb8jjVxh/r
RbLiL3u5sS85ZGuwhPaiw73HxRfS9fVqBQ6yfJnvReXtiu0N2CHuAaFbfveh
5ynH1O9iz5alZAe09qGZN2Ypt7FPyVW307v6kKPQqMQtbJUWEaXQ0T5UVKW1
9AZ2jWHyrzH2fhShlpLsjs0vUrE0Q6cfWawVu2+NHQVrPNXj+5EHM91FBTue
K/aYlssAkr7Daladpq9nnEEPd/cYQHEms4oV2KNWUuVv/xtAvHtkuEqxdV+/
Qmv9B9BAP+/3AuwG/RSZfzED6MqV+eXZ2Ku9cka+VQ0gE3rzyhjsqx21t/TV
BhE3w6XaAVsyQODrf1qDKPo/zkN22MVKllNpeoMoY9X92lPYAtebHKTMBpGp
3PoSK+xUwQ7TpWcHkVTosYcm2L0mI4I/Hw+iNSFR3krYBz6yJe/5O4hOW5Z1
jKXi9zudN7AtGUIbE8dfMbBXcVfHFrENIdWr8/YD2MUOURFGgkOorSS1kY6t
KyLvbaAyhBR+iwdXYm+8bmKh5zCEOFz8lFKwF/QfzmjUDqFEGS0pW+y35WsM
5JKH0QHl6B74gPfXG77tHmnD6An7ndbP2I53eJXg0zC63N9elYW9SZdb2LJg
GCkbt71Jxk5K4Rj2bx1G59vvyTzGTnzAGtHLzUD9TuOaLtjxxsyuV5cYyLdv
DASwb20YaRi7ykCsrxfFebBtWRnlur4MdLhxypMDWzJn8FNdEAOp6twTWXyP
v7829j5Y/pqB3kb/kujHjlnZqmrfyEBnjWJVc7Cf/vzlJ7NjBL0zYOyxxP7O
bj6mYDiCGKWwYz92j1n1SVXTESS/+YqiMfaW+t/bdx4bQTkDN6a0sAt7u5k2
niMobmCN8jrs8ZWLx5+9H0EvG1bs7HqHz7+9qmp865lIqd0p3QHb9eHHFyKy
TNQ89vSjDXZkzXZuSUUmCrD5lm6JTT+mN6iszURf5r6E7cb2dDKLM7NkojrZ
r6Oy2PGBDlz3gpnoanSQeHcKrl/V4H8hj5iINiIU3YzN4L/Q9ziaiSpiefhr
sTWjL+UlJjORQq58/w/sine+VwqKmMjLhntjLPbk2LLe0gommj0mavsIW2Rb
wMHaOiZ6y2ceHoRtl/tQqbOLiQ7qcTVfw76/hO/5wCATyffcHnPH/mAUyT42
xkTDCUsWz2PT7gt7zc4y0b8b15eexv5T8ax7cZGJvjoNLFhi/w9X/XHJ
       "]]},
     Annotation[#, "Charting`Private`Tag$8048#1"]& ]}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 35.34643541715155},
  AxesStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  FrameTicksStyle->Directive[
    GrayLevel[0], FontSize -> 12],
  GridLines->{None, None},
  GridLinesStyle->Directive[
    AbsoluteThickness[0.5], 
    Opacity[0.5]],
  ImagePadding->All,
  ImageSize->{329.79461669921875`, Automatic},
  ImageSizeRaw->{{180}, {180}},
  LabelStyle->Directive[
    GrayLevel[0], FontSize -> 12],
  Method->{"DefaultBoundaryStyle" -> Directive[
      GrayLevel[0]], 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None, 
    "CoordinatesToolOptions" -> {"DisplayFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& ), "CopiedValueFunction" -> ({
        (Identity[#]& )[
         Part[#, 1]], 
        (Identity[#]& )[
         Part[#, 2]]}& )}},
  PlotRange->{{-1.5, 1}, {35.34643541715147, 35.416426157547974`}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic},
  TicksStyle->Directive[
    GrayLevel[0], FontSize -> 12]]], "Output",
 CellLabel->"Out[32]=",ExpressionUUID->"0295452a-14a4-4bd4-8843-cc9ab896f769"]
}, Open  ]],

Cell[BoxData[
 RowBox[{"(*", 
  RowBox[{
   RowBox[{"Tc", " ", "drops", " ", "with", " ", "fixed", " ", "fine"}], "-", 
   
   RowBox[{"tuning", "."}]}], "*)"}]], "Input",ExpressionUUID->"ae3b750b-14d5-\
4e52-89d2-75e3af983480"],

Cell["\<\
Now we should fix the same Tc to scan over the mass, so that we can finish it.\
\>", "TextIndent",ExpressionUUID->"46a938ad-3f76-4d1d-ab00-60484f047b0f"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"sin", "[", 
  RowBox[{"6.73958", ",", "0.1"}], "]"}]], "Input",
 CellLabel->
  "In[849]:=",ExpressionUUID->"c6c6b4a4-95e3-4309-bdab-5f694857b383"],

Cell[BoxData[
 RowBox[{"sin", "[", 
  RowBox[{"6.73958`", ",", "0.1`"}], "]"}]], "Output",
 CellLabel->
  "Out[849]=",ExpressionUUID->"11eea865-aa48-4f79-a17e-dbf42ee47b37"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tc", "[", 
  RowBox[{"1", ",", 
   RowBox[{"sin", "[", 
    RowBox[{"1", ",", "0.1"}], "]"}]}], "]"}]], "Input",
 CellLabel->
  "In[851]:=",ExpressionUUID->"68bd2c68-b2ae-4b89-88b2-c94ef87b250b"],

Cell[BoxData["35.41572488554357`"], "Output",
 CellLabel->
  "Out[851]=",ExpressionUUID->"d908dacc-ed12-4cca-a23a-f186f6b5918a"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Tcfixlist", "=", 
  RowBox[{"Table", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      SuperscriptBox["10", "logms"], ",", 
      RowBox[{"x", "/.", 
       RowBox[{"FindRoot", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"Tc", "[", 
           RowBox[{
            SuperscriptBox["10", "logms"], ",", "x"}], "]"}], "==", 
          RowBox[{"Tc", "[", 
           RowBox[{"1", ",", 
            RowBox[{"sin", "[", 
             RowBox[{"1", ",", "0.1"}], "]"}]}], "]"}]}], ",", 
         RowBox[{"{", 
          RowBox[{"x", ",", 
           RowBox[{"sin", "[", 
            RowBox[{
             SuperscriptBox["10", "logms"], ",", "0.1"}], "]"}]}], "}"}]}], 
        "]"}]}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"logms", ",", 
      RowBox[{"Log10", "[", "0.01", "]"}], ",", 
      RowBox[{"Log10", "[", "12", "]"}], ",", 
      FractionBox[
       RowBox[{
        RowBox[{"Log10", "[", "12", "]"}], "-", 
        RowBox[{"Log10", "[", "0.01", "]"}]}], "24"]}], "}"}]}], 
   "]"}]}]], "Input",
 CellLabel->
  "In[109]:=",ExpressionUUID->"58613e4b-6b4e-4947-bd69-0984f9e02323"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{"0.01`", ",", "0.0002397560306391787`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.013436904130815223`", ",", "0.0003221578800956349`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.01805503926207193`", ",", "0.0004328804555831094`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.024260383164256546`", ",", "0.000581657319632857`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.03259844427549589`", ",", "0.0007815673676203946`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.04380221705435605`", ",", "0.001050184587615302`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.05885661912765422`", ",", "0.0014111229831271217`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.07908507486821956`", ",", "0.001896112474526031`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.10626585691826106`", ",", "0.0025477882764638753`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.14278841317896016`", ",", "0.003423438978397592`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.19186342188769212`", ",", "0.004600042854977942`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.25780504061150755`", ",", "0.00618103523061008`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"0.34641016151377524`", ",", "0.008305402019729385`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.46546801302008145`", ",", "0.011159899351290603`"}], "}"}], 
   ",", 
   RowBox[{"{", 
    RowBox[{"0.625444906691189`", ",", "0.014995474712154629`"}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"0.8404043250316183`", ",", "0.02014933613039721`"}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"1.1292432346572332`", ",", "0.02707461664614097`"}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"1.517353308446092`", ",", "0.0363802592557372`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2.038853093816544`", ",", "0.04888467094901649`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"2.7395873558428945`", ",", "0.06568796593801633`"}], "}"}], ",", 
   
   RowBox[{"{", 
    RowBox[{"3.681157265845454`", ",", "0.08826940261496133`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"4.946335727161925`", ",", "0.11861920763219742`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6.646343896470106`", ",", "0.15941795260151267`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"8.930628575729772`", ",", "0.21428331342528117`"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"11.999999999999984`", ",", "0.288116782395531`"}], "}"}]}], 
  "}"}]], "Output",
 CellLabel->
  "Out[109]=",ExpressionUUID->"ddd7316e-dd4a-4a7c-a0d3-f5642bffca19"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Export", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"NotebookDirectory", "[", "]"}], "<>", 
    "\"\<Tcfix_input.csv\>\""}], ",", "Tcfixlist"}], "]"}]], "Input",
 CellLabel->
  "In[110]:=",ExpressionUUID->"b4347670-6e0f-46f7-a474-bcb023deec4e"],

Cell[BoxData["\<\"/Users/isaac/Dropbox/Research_Project/2022-1-SFOPT_light_\
singlet/02-Analysis/phase_transition/Tcfix_input.csv\"\>"], "Output",
 CellLabel->
  "Out[110]=",ExpressionUUID->"0d2e608e-67f9-4a9a-9d4d-0f37d19d8a32"]
}, Open  ]],

Cell["Check whether the minimum is successfully found.", "TextIndent",ExpressionUUID->"04790823-db8b-4fcc-bfe9-f13cf98b72ad"],

Cell[CellGroupData[{

Cell[BoxData[{
 RowBox[{"{", 
  RowBox[{
   RowBox[{"\[Mu]H2", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}], ",", 
   RowBox[{"\[Mu]S2", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}], ",", 
   RowBox[{"A", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}], ",", 
   RowBox[{"\[Lambda]", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}]}], "}"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"strength", "[", 
   RowBox[{"5", ",", "0.17"}], "]"}], "*", 
  RowBox[{"Tc", "[", 
   RowBox[{"5", ",", "0.17"}], "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"FindRoot", "[", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"VhighT2d", "[", 
         RowBox[{"5", ",", "0.17", ",", "h", ",", "S", ",", 
          RowBox[{"Tc", "[", 
           RowBox[{"5", ",", "0.17"}], "]"}]}], "]"}], ",", "h"}], "]"}], "==",
       "0"}], ",", 
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"VhighT2d", "[", 
         RowBox[{"5", ",", "0.17", ",", "h", ",", "S", ",", 
          RowBox[{"Tc", "[", 
           RowBox[{"5", ",", "0.17"}], "]"}]}], "]"}], ",", "S"}], "]"}], "==",
       "0"}]}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"h", ",", "100"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"S", ",", "700"}], "}"}]}], "}"}]}], 
  "]"}], "\[IndentingNewLine]", 
 RowBox[{"hc", "=", 
  RowBox[{"h", "/.", 
   RowBox[{"FindRoot", "[", 
    RowBox[{
     RowBox[{
      RowBox[{"D", "[", 
       RowBox[{
        RowBox[{"VhighT", "[", 
         RowBox[{"5", ",", "0.17", ",", "h", ",", 
          RowBox[{"Tc", "[", 
           RowBox[{"5", ",", "0.17"}], "]"}]}], "]"}], ",", "h"}], "]"}], "==",
       "0"}], ",", 
     RowBox[{"{", 
      RowBox[{"h", ",", "100"}], "}"}]}], "]"}]}]}], "\[IndentingNewLine]", 
 RowBox[{"VhighT", "[", 
  RowBox[{"5", ",", "0.17", ",", "0", ",", 
   RowBox[{"Tc", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{"VhighT", "[", 
  RowBox[{"5", ",", "0.17", ",", "hc", ",", 
   RowBox[{"Tc", "[", 
    RowBox[{"5", ",", "0.17"}], "]"}]}], "]"}]}], "Input",
 CellLabel->
  "In[576]:=",ExpressionUUID->"1bec00ec-271c-4f3f-a8dd-06f7273ac3fd"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "7602.868580794999`", ",", "476.77973841000005`", ",", 
   "10.636175006840002`", ",", "0.12540903386192231`"}], "}"}]], "Output",
 CellLabel->
  "Out[576]=",ExpressionUUID->"0a556395-5d99-4195-9e07-f7df47f957ee"],

Cell[BoxData["48.45010900714836`"], "Output",
 CellLabel->
  "Out[577]=",ExpressionUUID->"064e7518-94eb-4a76-a1a1-2705c993e0a3"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"h", "\[Rule]", "48.45010900715298`"}], ",", 
   RowBox[{"S", "\[Rule]", 
    RowBox[{"-", "647.562844270429`"}]}]}], "}"}]], "Output",
 CellLabel->
  "Out[578]=",ExpressionUUID->"3f20aae7-1e7a-4d65-9b6c-80cd327afa35"],

Cell[BoxData["48.450109007148484`"], "Output",
 CellLabel->
  "Out[579]=",ExpressionUUID->"8280af7f-fe32-4453-b489-314afa9b6111"],

Cell[BoxData["0.`"], "Output",
 CellLabel->
  "Out[580]=",ExpressionUUID->"15e6467c-2604-47db-a6d0-5ef5c8d0956d"],

Cell[BoxData[
 RowBox[{"-", "5.820766091346741`*^-11"}]], "Output",
 CellLabel->
  "Out[581]=",ExpressionUUID->"2647fa55-9030-4f81-aef4-1fd313543705"]
}, Open  ]],

Cell["Success!", "TextIndent",ExpressionUUID->"2f0eac60-0248-401f-b530-dfc6cedeb210"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plots", "Section",ExpressionUUID->"db68c9f9-089a-4f7f-bb60-1fc7f5ccaaeb"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"RegionPlot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{
     RowBox[{"-", 
      SuperscriptBox["ESM", "2"]}], "+", 
     RowBox[{
      RowBox[{"Deff", "[", 
       RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], " ", 
      RowBox[{"\[Lambda]eff", "[", 
       RowBox[{"mS", ",", "sin\[Theta]"}], "]"}]}]}], "<=", "0"}], ",", 
   RowBox[{"{", 
    RowBox[{"mS", ",", "0", ",", "16"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"sin\[Theta]", ",", "0", ",", "0.3"}], "}"}]}], "]"}]], "Input",
 CellLabel->
  "In[582]:=",ExpressionUUID->"1bd87dd9-c963-4e10-b94d-c3399d7e51b8"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxll3k0lPsfx41JdSUULdbkpJSU5Yeprj5IFJdsV3JLqVS3FCKk5YdLVJdL
GiQke5Yac+0mmrE2mAUZJPtOlomUStf94/udc7pzzpw5r/OeZ57nO+/v+/N+
ns2n3WxchIWEhFhL738/V2Zniz81zNnfKbU+RlEtAhCfzvSSUFRLwtxq4tSc
Y5CJOY2fUq+gRsHcTY9VH9TJxywTGfZXjkExZttTgTOe5jTMtTFlzxTUXmE+
+73AMUCDgZlwnrJqUKcKcyIri3ZoXw3mfbppl3MM6jC3JTxRlDRhYr4mEsf2
NG+Aqvot1qH8MUCM9LWXH/rzrFhYR4x0h4u1sQpqHKwjRvq2RGMPtjoX64iR
/pFLPxyg0YR1xEivEgFlbe1mrCNGOnvMIX2Tmj+86Ny/7h6/A6ZrbybzrMox
R9eTAtnqdMxa2v++KjFz4tQHoknVmF2FVckL+rWYRS8qmzgdeI05nSs/Tz9U
jxn5hRj5hRj5hRj5hRj5hRj5hRj5hdaPzi+u0Jc7pDMAFXpdor/XCvy0on53
zDWIxRzv87vHoE4qZsn1tjpe5jkQtHg1IYzPhNxbFkM8q+dYN3JMvBmgQcW6
ZLjcr9raf2P95vaorkP7CrFudkKyaEG/COuFn0IMJE1KsR60U0SGcagMs9fW
zals9TzMz+fWLY8hFWAeqRK94HSgBK+v6H+unKmSasx28ecCbCtqMc8QnbWK
ql7j8xfEdlVtVnMXXP8P6/1xfT+uB10/+v0Roal5/aIqzMEXRjOf0mowK3P6
j4kwBH4Z79mQqWLZiP1CjPTepDXH7tqysY4Y6QbLlgd75wjyhBjp4qQ/Tqi+
EOQJsf5L7YWHvGmQHT8zfOOKgIWDk5tbprj4+HeXhHQ78gR5Q4y+P6bYW6Hu
IWBO8aacO3wBo/yh/wvNN8RoviFG8w0dj/KPGOXddvRBhjOxE+cp1MCQQyH0
4DydEytpoBD6cJ4mdE/62iz24zw1CPO8bBYHcZ4yWZYeNovDOE8an2flZxdG
/5Onsrq+sksG9TAQ6BchHzcI9GQ98YyWatglfmJX+6thWLYyZcCBWAd+j6CB
PDwKJpfFy0S1mLDAT7tAXp0EaP5E0Lds8m9sBM+5nYUiss+hT/pxUKM5F0yy
1Ry0ffLAkXbkPDGpCXpeeEebc6jAadBckxb/BjaSffN4LQUwsWF1hcmxVhAb
bLV+GFUI0lvjJjum24GUuK+IplUKtz2Z6++tXZpfsnqr39kv7U92yJXG+ndw
/EhMnYpsOdiQlc+EMnrgAOOLYdHkK5A9QJDY8aQPWjuzn/zGYcAXP0qYq8EM
2F036Hbw5sLe+C6GUOYMHMonb9Hu4cLcFClfzoYPNvL7vCLDmsBT2cU9+SEf
SNzMZyrcJsiNLXNLnuVj/69a//aQSuBg/3eNmoVRCS3Y/wn/vXeohDbsP9p/
z2Ujk9eRykCxW0NudkEwz260HnZo+yTYHxXj8ru2z8/g/RHt3mRE4eXCQvtz
04My4VDuaVT0weopnG2K1In2CYC9oYYHU8TGQM5bQcbyKBMcXIrqyMPjYJ+0
snHtngbwMnV5Uf9lDEoznLmeh+vBzWFNLHN5MaB5+SDdyX+Q8xKSDB0iq2US
IVLm6+q28RLwi/j6/+NacUB5O6zvQqwCNE8vi2a5fXasgc33XcErKh2sK1j0
o4nVcMsxyPxTSxq0ezy5WWnCghJX0aOMiOfQ5j5ktserGXYpNQWJpP8NikxL
Vbc6HlhuJ98bIxZDyTPWrKpZHrTSTjdK+UQCPcdgd/LRQhhZF+RZEBUNZ1eQ
dPwt6PDMXY9v6/MUtKzKL0inMSBbN5o3xk6G/T4ulrunXoOFMTc1XisL7Ffs
2fh4AxO+Jdlb3rTPgp+72/slqzmQGXy/784kBVh7Y/bZ3W4BTul4AjsvH05H
20XF6LXBpynzYxuvFYNGgWiqM7EL96lQl8YNqmEXiIaGSjpHlEOQBdtIMasb
zBzFOM++VID31lssCqEX9+0nZj/x4OFemCoeyxf2oUPPH+9SCZReOHV05vGj
bjrY82eZFEI/7uPJrbuZFI9+4KlH9rakVUKmJs0q5c9+mPN4+lmfXgmd8o+9
bRYHcF/XiGu1kM8Pgkb8n7nCKdVQOrLW02ZxCPd3OE3SX894GCoLf16VEVQL
sfn33GwWR3CfZz26Q9yxaRSEerdVXDr3GgY/7zS0ZSaAtMjqjWayt3GfqGhk
MXfulARiSNj7HXIRgPpFdVXtbhdiw3/6JVE64W2kbjZYKTRHUAlNuG/EChh3
qYRW3Dd1v1ADqYQO3DeigYP5JjKdoCTV1x/uTIOv8uVKEjo9cP12TOKjk6+A
7/vTsN+JPjjzQb/+UjgDSEo+D+T9BkBqRYytErUKxlrscoeCB8D8wsC8/ssq
EGL9ur/cfggyNpvmVUXVgNYpPfOUnCEIbzzCmhqvAfnsejpFcwR6zLVc7l6r
gx20wDUSb0bA1+jjAl32NaT/xKshz08I+myL8kUb8QkIaboeSTvVAEO2V0Ul
3kzANZ3FrXPbG8HLSP2XUMZ7iGR+8OfJsqC8ydQiVPk9fBNrzFDxbYSgN5aV
5PlJ3HfilGUrJXQmoV5v8PyVv1ggrDkSRqBMQt7J0Y6SVWwYstybUf94Cq4t
k71P+8oGpzA6g2w6BbGnv20rSmfDcn+/Ijm/adyPR3Uq9ourT0N+jejsdycO
lCX6eSUbTwOzIeQq25sDhPUO606WTYNDTOATgykOfE9RiLcunAb5zvV2qWMc
OES2CHP1ncF9+mELQdlYcwZKR+Pm6ee5cCSgw76tbgbGlu+Wfr/k670spTMv
vfm4PwOkaGdeKvKh7XaEdapDEyQ4xaw5+ZYPhZqmLToyzfBRVSuduZwC6Plh
9sAnYxfiK0DPD66ePF2XpZ5Bzw/Fx0nB8gvj+H58f0iXRvnH9/j+2z3w/TKJ
uSl8v4361M/veJC3GBf4M2vGO6bfgpJKjrQziQZ13l/vyC/N2z2rpDw1luZt
xznvICqhHd9v/QN7flrk
    "], {{
      {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[4], 
       AbsoluteThickness[1.6], Opacity[0.3], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJwBxAE7/iFib1JiAgAAAJEAAAADAAAAKikhPC4tTk0XSkkVREMPPzQzPTAv
PjIxJB0JIxwIJR4KRkURGwgHIhsHRxQSRRIQSEcTTRgWHgsKSRYUMzQHHAkI
OCYeKCohISkPHCQJKCENJyAMOiggGyMIDiEPNCIHRUYTk5QcJh8LOScfHSUK
NS8ELS4EICgNLzAFNzMGMTIGHiYLHycMNjEFQxAOXs0pHwwLIQ4NNj4xNz8z
EUUQUzYFNT0vJkEfJUAeQ0QRJ0IgEEMREkUTVUsVFEcVR0gVFUkUE0cSFkkX
S0oVSUoXSkwXF00WD0MOGE0ZTU4ZY4QuZYUwZ4YyMlQGTFYXUjUEVDcGHQoJ
lpcdLlIEMFMFTlAZn2IaVk8XT04XiH4rkJEbioEsUFEZbsFMcMRQIA0MWakk
j843mFolxYcBml84qqtAxokCur5It7VGrq9CybYRyrsTdcwiy3dVyLEPwG1L
w29PnmE6srBEzXYpzFcizmg3lVkkdn0pfoArfMgPklgjecUBnGA5x4sDjmY2
o6I8WKgjoKE7rK1BgscDgYMss8kRsLREYnoabJYdvLpIuMoTjWQ1V6ciapAb
a5Mcf8YCtblGElB7QA==
          "]], 
         PolygonBox[{{114, 141, 53, 82}, {115, 142, 54, 83}, {165, 103, 50, 
          62}, {166, 105, 52, 63}, {137, 138, 44, 2}, {109, 191, 74, 75}, {
          122, 121, 1, 26}, {135, 136, 43, 1}, {140, 163, 60, 45}, {169, 108, 
          29, 36}, {168, 107, 28, 35}, {148, 149, 36, 28}, {145, 146, 35, 
          27}, {131, 130, 3, 44}, {128, 127, 2, 43}, {139, 140, 45, 3}, {164, 
          101, 48, 61}, {162, 99, 46, 60}, {111, 194, 78, 79}, {93, 123, 42, 
          40}, {167, 106, 27, 34}, {151, 152, 37, 29}, {91, 172, 65, 38}, {
          105, 117, 34, 52}, {185, 184, 19, 70}, {182, 183, 70, 17}, {134, 
          116, 84, 50}, {133, 115, 83, 48}, {132, 114, 82, 46}, {104, 166, 63,
           55}, {193, 120, 86, 76}, {102, 165, 62, 54}, {171, 153, 30, 64}, {
          90, 170, 64, 37}, {175, 157, 32, 66}, {92, 174, 66, 39}, {180, 179, 
          17, 68}, {177, 178, 68, 15}, {100, 164, 61, 53}, {157, 158, 58, 
          32}, {187, 188, 72, 19}, {189, 203, 85, 21}, {123, 94, 41, 42}, {
          125, 124, 15, 41}, {120, 195, 79, 86}, {119, 192, 75, 85}, {116, 
          143, 55, 84}, {97, 93, 40, 58}, {190, 189, 21, 72}, {173, 155, 31, 
          65}, {96, 92, 39, 57}, {196, 113, 81, 80}, {194, 112, 80, 78}, {191,
           110, 76, 74}, {161, 159, 26, 59}, {155, 156, 57, 31}, {95, 91, 38, 
          56}, {153, 154, 56, 30}}]}]}, {}, {}, {}, {}}, 
     {GrayLevel[0], AbsoluteThickness[1.6], 
      LineBox[{26, 1, 43, 2, 44, 3, 45, 4, 47, 5, 49, 6, 51, 7, 8, 9, 10, 11, 
       12, 13, 14, 16, 18, 20, 22, 24, 25, 81, 113, 196, 112, 194, 111, 195, 
       120, 193, 110, 191, 109, 192, 119, 203, 189, 190, 186, 188, 187, 202, 
       184, 185, 181, 183, 182, 201, 179, 180, 176, 178, 177, 200, 124, 125, 
       118, 205, 94, 123, 93, 97, 158, 157, 175, 174, 92, 96, 156, 155, 173, 
       172, 91, 95, 154, 153, 171, 170, 90, 152, 151, 150, 108, 169, 89, 149, 
       148, 147, 107, 168, 88, 146, 145, 144, 106, 167, 87, 204, 117, 105, 
       166, 104, 206, 143, 116, 134, 103, 165, 102, 142, 115, 133, 101, 164, 
       100, 141, 114, 132, 99, 162, 163, 140, 139, 199, 130, 131, 129, 138, 
       137, 198, 127, 128, 126, 136, 135, 197, 121, 122, 98, 159, 161, 160, 
       59, 26}]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  AxesStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{None, None}, {None, None}},
  FrameStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  FrameTicksStyle->Directive[
    GrayLevel[0], FontSize -> 12],
  GridLines->{None, None},
  GridLinesStyle->Directive[
    AbsoluteThickness[0.5], 
    Opacity[0.5]],
  ImagePadding->All,
  ImageSize->{371.8105773925781, Automatic},
  ImageSizeRaw->{{180}, {180}},
  LabelStyle->Directive[
    GrayLevel[0], FontSize -> 12],
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 16}, {0, 0.3}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic},
  TicksStyle->Directive[
    GrayLevel[0], FontSize -> 12]]], "Output",
 CellLabel->
  "Out[582]=",ExpressionUUID->"34d94513-a59e-404b-81d4-70d7aa067f60"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"ContourPlot", "[", 
  RowBox[{
   RowBox[{"Tc", "[", 
    RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], ",", 
   RowBox[{"{", 
    RowBox[{"mS", ",", "0", ",", "16"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"sin\[Theta]", ",", "0", ",", "0.32"}], "}"}], ",", 
   RowBox[{"BaseStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontSize", "\[Rule]", "16"}], ",", 
      RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}]}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"BoundaryStyle", "\[Rule]", "None"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Style", "[", 
       RowBox[{"\"\<\!\(\*SubscriptBox[\(m\), \(S\)]\) [GeV]\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}], ",", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<sin\[Theta]\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}]}], "}"}]}], ",", 
   RowBox[{"LabelStyle", "\[Rule]", 
    RowBox[{"{", "Black", "}"}]}], ",", 
   RowBox[{"FrameTicksStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Directive", "[", 
      RowBox[{
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}], "}"}]}], ",", 
   RowBox[{"Contours", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"25", ",", "30", ",", "40", ",", "50"}], "}"}]}], ",", 
   RowBox[{"ContourShading", "\[Rule]", "None"}]}], "]"}]], "Input",
 CellLabel->
  "In[583]:=",ExpressionUUID->"8bd1f3f7-d846-4b3b-803c-80a05f3bab69"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJykfXlUTt8XdxIZ02SmVCKhTJm1L2UeSoZQSqYkpGTKHCGZUkoSGkgDzYMG
DZrH52lUkUqaNQuF3uer9nnWPb3v+v3xtr5rfdfHvs+9Z5+zz+fsvc8ks9dE
64CggIDAuuECAv/9f6fSYp/jUz6pOu6avuNztBY8mLyW918rCPT+9WBBBnFJ
6FHef0MZtlyc4H/iklEEF//7h/EE/+3+70+aYNt/v5clWG7tvzcSHPLvfVMI
Xv1PPI3gon/vn06w8b8CziS4528WhQUEalQjL3/Y+pjoG+A1r/KjoSOgvohR
X7ZcnGDUFzHqixj1RYz6IkZ9EaO+iFFfxKgvYtQXMerLxgICO84nap9P9KT0
9aT09aT09aT09aT09aT09aT09aT09aT09aT09aT09aT09aT09aT05eu38tJU
n693Aoh+iFE/tlycYNQPMeqHGPVDjPohRv0Qo36IUT/EqB9i1A8x6scuL7Zf
GNV+YVT7hVHtF0a1XxjVfmFU+4VR7RdGtV8Y1X5hVPuFUe0XRrVfGNV+YUS/
iXnLef9FE/0Qo35suTjBqB9i1A8x6ocY9UOM+iFG/RCjfohRP8SoH2LUj11e
tM84yj7jKPuMo+wzjrLPOMo+4yj7jKPsM46yzzjKPuMo+4yj7DOOss84yj7j
iH7nJN+sLrdIIPohRv3YcnGCUT/EqB9i1A8x6ocY9UOM+iFG/RCjfohRP8So
H7u82P+Sqf6XTPW/ZKr/JVP9L5nqf8lU/0um+l8y1f+Sqf6XTPW/ZKr/JVP9
L5nqf8lEv+p/42Ma0Q8x6seWixOM+iFG/RCjfohRP8SoH2LUDzHqhxj1Q4z6
IUb92OVFfsmk+CWT4pdMil8yKX7JpPglk+KXTIpfMil+yaT4JZPil0yKXzIp
fsmk+CWT6Of9X3OO5hD9EKN+bLk4wagfYtQPMeqHGPVDjPohRv0Qo36IUT/E
qB9i1A9xpcEe8WFBLXB9oRz3YD2X4k8uxZ9cij+5FH9yKf7kUvzJpfiTS/En
l+JPLsWfXIo/uRR/cin+5MKe6lVhL7itoPTz10iNEblEv7LUw2f3CeUS/RCj
fohRP8SoH2LUDzHqhxj1Q4z6IUb9EKN+iFE/xKgfu7w4HuRR40EeNR7kUeNB
HjUe5FHjQR41HuRR40EeNR7kUeNBHjUe5FHjQR41HmB5K3rHgxsE97SnK8FS
//rna4J7fh9McA8/RRDcY++xMCchvrMruQKev+qMX9YZ31tfXbCh/XTgrtIU
gi2t/27anp5OcLiUVf3m8CyCfxYNSHiZwyFYyeT2HreYXPI9jBPY+mC806q6
ptxC8g1PP/w96tvTPkoMYvZ4coOqH0+qfgLId0/3++8JXyLPritbEmQdTuQ9
5Q0ictqfouuD5nu6fth83wVsfuzq/V1Or34DGcQoZ/e3gQzWJ8rZ9qtE5IjZ
9t9F9EQ5YpSz+UyQlAcxvh/rH3HylyVHnbRaiT1gednfr4fHTv/9pWI5mEU/
/hkMwf+qa2IeqZee9jxNcE/5bAnuia+9Yc1//zCBN45sHWfWqOYPi3MNSmzS
csCh/5N1gzsDyff/a53xcim95esg5fkzpSV/a3AdMFLTZs3vl07k/5pLlI85
SpPmP7yeQfB/rSPxJ5PggTxl7p/IJvjGk64Gwb25RD/EKEd9v+WbXOFKPQRV
zVdjp245CTOG7hw4ei8X6HgMy4v6oD2inC5fw+b/WoRLMNY3Yrb9dQCb/7vI
84jp8rPtrQNoe6Pbk21vygy2Lz4/Sr9lpH5LItGP7l/YHlifbP4XoOxNmbKn
1j7v389rPRffJCJH/sLy9rzPo3c8qAcPHTnnhKuJRN/Nwic1hU9yCH7q67Bg
1Swu+b4A+4/B9sD+g+Vj89lpUj62v98FbP++g9gzlh8xfp9uDzb/txKM3+vh
r0iKr2Ipfy+WPN9TX7EUX70nzyPG53vs5T2p357xyoWy53CCe/jkHcE9v4un
4qMbVH3x648eb1C/PUIxl4VsvpDyYH3S/I7txebbYPJ9xPi9HvsKoZ5PIO9P
eZvN66GJhB/Z9iZIMNoL2hditn3V99ZvEGUfYZQ+4UTOzl90UPXbBez4v4PU
N8rRv2A/H0DxRxDRjx4PcPz799oHQkQ/dvuFEczOV9RT9YnjS3Dv+yQJRjna
d0/1iTJoj/g84iP/PiDPID/8G08mj2VwPO+RT+yDsb16npch4yG+j/Z/6PGM
bb9dwK6PLsIn2A49zxv2yksIRnlPfV2hnqe/x3+ezcf1ZHxlt68rsXfEWL+I
2fbnSeRsf1SQYLTnHnvwJ+2Dcvb7A4ic3f8EGXb/G9qn/en8JMrxfWgf2N4o
R8zmx6EEs+05ltgTbW+0HDHaDz6PmM2vQxk2/4oTjO/H/oHlpb9Hy9He0V7x
ecQox/djf8D3IcbysvvPRCLvGS/liBwxylEfur/R30M+Rjkby1C/lyFyjMdQ
3hO/Te2D8fnQ3vgPMepDfx/9EXb7pRH98HnEyBf3/ymkSOSIUY7vR77E97Ox
PMU/8kQu+0/fGUSOGOUYzyLG+qK/h/yMcuQ3rC98HjHKQ3vjCTbfKxE5vh95
COsH5YhRjvWDcsRsP3JWH/5FOdue+XyMcsQox/pCOWLkb3wf+lP4ezp+Z8dn
9VT8WkXxcwXFt/WUP9RF4uMe/ZWIv4iYHf91kPex/eMgwt80/2L8hL9HPmb7
f/6UPvz5S3o8YM+XCfSZD2N/X7QPv9PzRzQfs+UD+8zP0PyLcuQTer6G5l96
fofNj2P7zH/Q/MmerxAl+Xy2f5JA5cv4/EnPN9D8Sc9f0HyJ8QG7vvj5fpov
6fw+zY/0/AHyF5aP5js6f07zHcqxfJgfYdc/P39N8x2dr6b5js6P03xH54tp
PqPzwzSfoRzLS+dvaT6j87U0n9H5YOQf1B8x6kvzEzteL+nt33cI3/TYiz3l
/zmS59nyCgq3UnzSTL2/g+RfsL8jf6Gc7X+2EozvQz5k88cNwi/sfB/yiS31
/huUv8rnP/wevh/9XYw/UN5jH4MZjA9RjvyJcsQox3gP5ex4K723fNpUfsSD
/B4xWx9Pog8rv8V7P2L8PfIz1hfbnxYgfNtjTv2Yt//s+Q0c+/eBP73xu38v
Hk5+j+9Hfu75vRjRH59H/dn6BJPnsT7weazff+bePZr8vqf/jCPPo75s/3ww
w45XJcj78Hsox+9hfuJf9YaO7PN7lPf0twnkfT38IEWex/Li81hejKN69J9E
5D38IEPkWD6M/7E8iLE8mI/C8mB+AsvDzh/IEnkPf8gROZYX34ffw/wQlh8x
Po/5RdQHf4/64PM9/DaVyHv4SIHIsfz4PiwvYiwv5vN63j+NyIP/Pa9I5Fhe
zLdifSDG72H+tcde5BnMh/bw7RTyfM/3pxP5qn/PzyRy1B/fh/pjfI76oxz1
Rzn6ZyjvaV9lIsfys/OzsiS/g/WDcqwflGP9oBzLi/la1B8x6osY9cV8Lo4H
NL+y/bWy3vyOXe/7f/Xml14RPsfxg/17Fyo/60HxsyfhD5rvkI8QIx9h/6b5
meZj9Bex/6EcMcrxfTT/oD+Icpo/6P5M9zd6/QO2P9tf4+cn0N/D8rHzk6OI
HN/P7v/SxL9A+2LnjypIe2D7Xf3Hdz5kPECMchwvUF+6/pDfsDyI8fvsfFbJ
/yP+uULls05T+AbxX/B5HMd79DnaKw/slZuS8bdHfgEOWR2TvrQ7HYKGd9+9
dtWUypeZ9tZvQZ/fs/NlgWQ8x98jZvsztlS+0p7yX/j5Mjr/jf0D/SXEKGfH
Wx1Av4/2L+h8Gu0fYHui/4CY3V+Dif3T+TV8HuU4XrP9iVhiLzg+or0jRjk7
/yVBja+SBLPtMZH4u+ifYn0hZtv/HSoet6XiYxfSnrh+F59n4w7iT2N/QXlP
edspf1qQ8BX6Y2hXiNnzz4KED1GOfIv9ld0/u4D277C9sfzs8tb2zlc8JPXF
Xq8sSPgav99TfncqPxhC+j9i5Cv2fIhEH3+QvX5RlPArtjfyKY7/6B9gPIqY
7T+m9P6+kuiHGMcvxHS8g7hH3xbyPLaTAOuvmdQntjfWX4++Qxn2/MxQKj6Q
YNjzceJEjvWHcqxf9E8R0+tVsf7Y8XU08e9Qju9HzF5/xfdvaX8X8xvYXzEf
geVBzI6P44h/hv0Xv4eYvV6I7+/S/i+9nhPHX3b+PoHoy/aXRxGM30f+QP1o
e2PnU2SofM4k8j4sP2LMB7Dz5TJ9/Gl2fkWOWr8pS/nrspT/PZ7kP7C86M8i
ptfvsdfrKPfxT9n5ESmSn2Dnt/j+JvqfqD/bX5Wn8unyfeoTy4+/R4zfR32w
fjC/ivXBnm+TIflRNt/doPr7NeLP9vTXW+T5Hv2fk3wPvq/HfpUYBZ72Yhvy
KH/EFFa7XB23PjMd2PkEedLeyGeoH9ojYrRXjHcQo/5YH+x4SYb475gvYuer
5YkcMXs9RCDxb3rKk9RbP+a99cMB9K96+Ppbb/3cJ+MH8iN7vHlOreeNIP4I
m98G9uFD5AusL8wXIr+w48lRBKMc7Q3rj73edwKpP8T0fA7Ksb7o+QHMt2P8
hPWJ+ULEKGe3Ry2xN3b840jqjz1+dBHc0z9mEvvD+BIxtg89nrNxF7DzN3/o
cY34I1h/2D6IMV5hz8dhe4r3tsdbyt4jCR8hZvMXX472gu2D4xNiHM/Y82P8
/Aydr8fxC/sPYna8xJfj7/F7dH4exyN2/oo//rD7+6Q++XaUYz4CMeajaczO
r08l/I/5avZ6Ntk+85Xs9WBKffIZGB8inyJGOcaTKEfew+8h3yM/0/kIfJ49
H5xL8huYj0D+Yucr5Kj5PlEqfhpKMMbrdPzDzqf2Z9j5ClEiR3+UnV+VpNYv
iFP5VDEqXypJzXdJEP+PzU/h1Pw2rj+VJP0H64vtj8lS6wdkiP2y5xviSPug
f8We73xPrb/m5w/RvtnzmwmED9G/wd+jveHvkY/x9+z5HnnCr2i/7PF+CjVf
PqWP/aIc+RYx8i0+j/0N/R+0R8zvYX2hP4T1y54PlyH9BcvL3o8xlbwP53Pw
9zjfw56fV2DY60OVyPuQz+n3seebppF8CpYHxxesP8RYfzgeof0hf6Oc5g86
H8gej0pIvIP+E44vPeICkm/BfosYx3v2esChBGP/Q3tl8zeHYHq/CMqRr7F+
2HzN34+AGO0B9UX7odcrY74J6xvtGdsT7R3bmz1fOY3YO7Yntjfmb9Fe8f1Y
b/g+xPg+5Fd8H7u9Ffvsl2DrP53oi/XMzo/kEIzxLFveSjDyHztfi/ML4RR/
xxN7w/rA/kP7s2j/WL+IUX92/5hC+bcypH2wfdnrtwuIf47jAz3fR8fb9PwD
9l/EaG/o32F8gfbOzg8IkflR9I8Roz+G+R/E7HyRGCkvjkeI0f/qKW8kqV92
voRe7yTLsNeD88d3bD+sX7RXfB4xu/4HM+x8pgS1/1eO2t8mS83HTKX2E02l
7LWC8BHyT48+j4gfy+afeIKRvzBzw7YHbco+DEn7s/cT8v0XfD87fxTfZ76c
PX8tyvzfsSPhP6w/tj/hSs0XeJL6Zc+/S5L+x+4PEdR4GE36Gzv/I0eeZ7dX
NGkvfB7Hd7Y/Ik+tN5Qn/R/jJ7Q/xFi/yL/s/qpE2asSeR7HS3y+p3+0kPgK
+Qv5GzF7vUAbFf8rUf6tEmlv/B7Nn/i+ghZ1o7cWL6D0yjuhH35W1PqJK+R7
yA89OJ1gtj1do+av+OclsPPBYn3OF0C+wPal92fjd9Ce6f3M7HwWfz0Pxkf0
+hz2+MbfX4rtSa/HYfN53/U4NP/T63FQjuMje7znr8/D8rD30zXDqX/7vzJ6
8Idi+K92xCb05lNlayA0aral9cEUwjfo36BdsPeX5ZD50dKgitRTq/PIfgDM
F9Hfn6Z7L/67QgY8n/lLbKPGR7i74sJi10dRPeV/UA3/xHnJxF7nzf3vL5fg
4v/Y934v34q3g/x//749j8ivfPhvA1Quld/pzUeVVZP9IVge9D/w9/+NPnb5
OQSj/SNGf9LOS/7PtPntZL8Pytn+Yz3Zb2PpFTS2f78qWLZtgZ7VsiSojhAZ
sDGjCmZaH1l+YmgylU9zJfk42X8GmABv1nC2XZ1aTe2HqSf7S/D5nu/zcU95
Ewhmz6dWATufK0TGK+RXev6LXt/H9lflqfXP4/v4CxgPoL9Krydjzw9PIvyD
/I31i3zLzncOJPEv+g/0fDm9/pI9v86PJ5Hv6fWHWH84PtP5a9ofRoz9mZ0P
nkbxuWKf/Bu9/g/tFf1bxOiP0OsB2fY9qc96QGwf1IdeD4hy5H/2ekD+fm3U
h73emb8/G/VD/x3LT/tTbP9rep/8I3v/XQHJxyJm5xfL+sxvoT6NzOdElwNN
kJfQ71JBQjaR0+sF2PNXJcCen6wA9vxxBdlvdMrYdJ6fYhF8MNcVW238lsjR
f7Dbd35xgH45LHX7YedH8hVVvfsN4wn+9zm7RPL9nvxYKMHo73z8D1t9gub9
IVGzc6IJv+E4ixjtA/mq6o4Yr8h8vmOt3+Y29uzXUsom30NeRYz6Yv3h+ITf
+1fdz/hY7N8Peu1vQg55n1jO6ZkJYsVQ1et/6c0b/iWCx29m8l+Wb4zml5fa
L9u7X9kdkv0C849YJMLnF8Lyf/xeQI3RmYkvyuPBqkEk2UvUnZSPtX5JohnM
lg0pdFfnwH37XybvLJqhYmh8/qvXHFL/rPVpE2PhxXHeD9Sew6ATuiaORRHg
137NeO0EZ3h1ZfxesX22UDJl5OEpd3F9hRuZT5iw6HG2cOdtkO6dTzDR3HH/
mNJp4n9fnJfRYXz3Jgy95r9aYKsR0Q/t2/5kRUfpeH+efSt+56rdh8inw9QH
dD6BCb3+ffyiC8aSY12h9tiJueZq/PUhPfxoT/lrjnBB0MQnpzwQug+sPuIt
akfah5Uvymqi9muX9NZHCATX6VxrnOEFztJXarWNrSHcfoXbAeGXMEvrwEIB
1etUviQMzvyNPOnZ8bHXX4gmctZ6JYUP8F9tSR4MJ+MhvPg6e2vtYwj77/Gu
KGq9SDqJfzSG6nEtim6BnnIItI03hlN3lPzDb76GyPO2U9q5d0h79PCDBbw5
pjnG2yIEJH2MVX/68f0b5At8P/IJYrY/FE/ksQmSiTO1W4A1X7iiBQT/q4Az
XAiwXxqqtbuR7A/vcS9rYca7iunvDvfGp2dre/Y3x6XCJP3AvRdmpELVupwh
R4y9er5XlAbtSYOmdI73gU271noUDiiAo4sfrqtZFtrTfue/gf/PW34/czLh
Qmma82CnGrKf81/3KW4j/su/2h9SR/YXR6sbSaXJtrH3Nz9rBlb++08l6Q89
9lED/+ZfJFN79KmvBYv/ur9KGjh0XnZruZEJfpFKEkI5vfljiTaYPpzXgzf3
+lOH2uDifw7pq97vqbTBiLrnInXmuTAjhFnz/VEbuH5/9/lScS6E/yePbgPu
Zx4hd+dCsmSbedTmWrIffIKIRczKtErYHrCqKF8wscdepCt79N3We76QeiXx
d/7x5cgmMh6M1vwoqNjVTPi/hx9aevhibe95JC/LyH7ff+P7wq/AY4/jy8x6
93NxC+Df8D20dz1bYhlkXeAp2NK7/nZUGlnPbOZ7VuDz5gy41fztY/4j357y
vPhG+DOcE2u8bvg3GKdwX0piX2aP/3KoBdz/459nvfMLFvk99e4XAm77rKs+
WZaBq9zWGZve9O5/i8gn/euV8tbszq+FZLz4V98nCsh48s+P3utK4rUBCe1V
UV8zQNR87mclpdegVdfk56yRDENsR6tnW3uC+pzIRZLt78G482N0urwHeItq
ORRKBYOmhdDv95r28LRQdvmvaeHgsn3u+x8HHkHBmT/mJ7Mryf5b9/HVeldP
8+wVz0vKG31+1ammXv5Pg4I/HQPOOLX08mUmKO3yv56f0AArSvuPFEtJgo9p
4vZBO5pgqdno5ILaVDhlWfkl6k4LtPjt6vdBLhP0SloNxrfkAmf5H+ZSqgvI
JGR9WdCQDzVyBUdrvrrDjRgTiVlqxXB5cpnXA2cfUJt41TP4WgncWq4+2yP4
DXDdhUZrbi6DzcsUDNUVg8HHf4y4smMZnN5/ZplmbDA4fHs1VW3lF0j6bul/
ZUYEWA/WPtDQ+hVG+xQvSXgUA7syD0yfsuIDiZe2dN89y+n+xD/P4KC4js3W
xzA1b6S7a4MvzMpYV3Gmsry3PsJgRzazyeQijyft/5iE9Q+GI+e3L6wubyT9
mTH0vMd84Y/nH3WL/B2GthG+fjIvy/icYDvcmrl/3bf1OVCv5nTjbVgrKN80
uOpSzgVm6ZEDZ6SbIHuw1rXmjakgVT289bNsBpyNGZNckXkezIUdRyqd4cCx
Ga+O+Mk+gGXGqZXGOzngt7d9unoTb/y5fPaa8qQimCv8LO31Ui+4PX+X1/En
nyDo8vLvHhMDoKqq6V3w8ShwkdbbX2oeBPv+uEf5ZHwDYdt6rTutyaAzRWFL
/YFGWFXtm9agnwKjb7lcdr/VBE1xChtDfdMgQlmq2LmmFVwKxwm8HMKBe0Wv
XhzY3Q4RkxNDN2jmgUOltcm3ZU388ytM3qmZL2sh8dyckFDxpqg2GGesk27G
8zevj2oUXPirDT4VpC3nTsmBMx5bODc12yE47MTvYT45YGV9MPGyTjtY9MZz
X7VTFIyv8vwB9/aWd0PSgaNlYivymx/PxWwdE/XRrQHuJhi7V9glgX0DYz5B
tY3wpWSE2oY5kwvgsf7g7qwxL2C65lSxp0w9fDWydLo2JxGiB7gviGtphet1
BYO+TuaATMxb59qNbfBDt2ZuWjQHKkBvygy3NrjuNHvpkhdc+MWxihqk0wbd
LZIj7tdwwPmzwqBBsxpJPCsl4tdgteULOc+gMf+X5kuDL+AXVz4i4koE+FSZ
he2/UglQM+tBTVUU6NduMD+6uIqcZ3B9ROVNtTtVwC3ovJR0LBaSSqOvt6tV
Q5BidPX0SXHwVtJ5l5luDXy7nDmGc/k97LUVEx5zogY2/mgq9Pj9Hqw1HquM
S8yFrP7X8palPIH5jYNvRXRVwgqT7ZZa/tFwanJXyarGasgRsl74RzgeZJZ8
FHQaUUnOIxC6+ND5lkgNaJ1d6iSjEw9+P7Vfj/SthpuOrU4VPP9YmHPQ4Jl9
HeFv7mX9a0FSaXDONGbIpbhdYDp21KIjmuRcFYHkwuLMFV84ENjr7/xea/rL
qCUVLC9arVz/eAusevxRrIY3rqA/ETL+nXRXYCbhP6Mysdf3ZmWAsOfzghvv
LsD1iONawllpve8zhLVRBYckhqT19ldtCD7kPttyWkYvv10As4SdsdYSHDJf
4HPkSdobIw6s6s2PqpjmZjpPySb5M/vZ558qe2eB42cRvz9XreBYXN3l57Oz
YIP6W/Gdj67CypkSh/etawFlufJ2bb8MOOMSO/vbuBKSf/mp22Q2/WwZyb/o
emkfz/3dDDbbrghaCmSA7bmUb3MUW2B7wcg/Pjsy4MAGTauvu9ohWvJWXoJS
LpzrlFXS59n//B2rChwX5IKA+QrI02uHy/4+55rtcyF+SaWEOw8Lj/hc4X87
F36YlN3qz+uPA/WLj4/flQcHV3pvvaPfDlt3N0vsq8uFhBWZqy7zsKi/+AzL
1lwIenFfCba1g1j57JvHPuWA9bJ3D45v5fW/HftFBuTlQMnoqGvaO9phbuXt
WxkDcmHS27NnNXnfq2o+lH5qXh7YLOqf9Uwzv7f+XEF+4evqF651wN1scHee
TQLobEi5fQNq4IfBC9FXAu/Bn+t58r5GDTR2DnZYpv4emtbKPlO7WgcvrkdF
Sx9KgJMZ9rstD9WQfMg2E+mv8IYfH/rPynVsu9ZG/AHrAQ3d889+hrIP4S7G
/YLARdhJbX7hZ1g1PUpt1YEgqDS/xomVKifj7eSxu1OWe/DjwXED5EU8jn0l
8d+L81JrvvL616egcy/H8PqX8sb7O4Ska8A1ybex7HY8HJfasjBLqQYOfnQ+
VRcXD8yLaJVr08tBaH6btrh0KOwaa+tZGtsMZy+sznxwOx26c1+67fJqhvmn
7k98qpUOaQ/rcqXWtcPuQMlls0/kQEJx2q/hG9ph4ZalJ8osc2CRxlHnZRtL
4eH6rD8bRALBvFZTpHX4Z5CV3Bxz0C4Q5BwKR5t/byH+yNBDIv2WurVAQbC+
303dTNBSerx8TnML/LJYnVDYmAnTv41tSvrVAllHuQWJg7PAs+Fi+UCev4bx
fsfp8tMJ8m0k/vKtmjtMekQ5aAcMez8pJgSM0lvFwtfx49EPvnouGWO5UJfy
sdCm8TFc1Pi770UdB6R/ljgM+OAIix7XmiaLcWHEOqsOZb3HsCXt/pubvhw4
3+u/K1wcu13kAi9unajgunyQNRyeHxDh15INfvebouOSbMC6sjV/RwkHpr3b
8zWr1QFeDduh7xfKgW0aXs+8u+xB4lSesYU+Px7+9TTc/4seP/83/JnWofW7
+fk+fbm1HZv+PubxrverPd2+0GRyyNbwvRP8OORg7VbiC4U2TZ+Mee+Tuj6n
27ArFx7tymse+aMFGtLEq+L7ZcHRmBKbK7z+56u25trMxbnwtfDBmm+873XZ
qq62ccqF7/VGcw1533s0pNOibH8uDP34TcSK975EjYrQ/3iXjj/o8/tKpL4X
TONhbI/HKg/sS/XLYdixqyc/8fiVxoYKywVDJHj9AWSHczfm9DmfiY53C5WK
lx04Vg9xmZbjO3nxOh2/0/JR3fJjlW+3goHLvtQbKTkwc/to1ZlrK6BA5OP2
+rPxfc6noc9n01x+K+iYnxvM4FTKv/tk2Sf/TJ9/tf/t6Nv+x6rIeUc0pvOz
dL6Vztd+PHDx1xKev0jOC6Lql85n0Pkgun7o87f8xOPKLp6sgIW9+SC6vun8
AX2+OC2n8wV0e9L5g/+VX6Drn47/aUzPD9H5AHo+oE6mWnnkslbCB3T+hI73
6fzJ+bGdG5L2BpD5SBrT+RQa0/kDOt9Bx1t0Pk9BJ0Z53aNC2H7v3MJHeWGg
eFfkzwCjUsL/dxSjZ9tzGsDbMfb8cXlefHOdGT5MrRk0e/lxmFjT3YfyDXBP
U/Hu59MZUFlkP8L7xGcQ43bN6UiMgb9z5z82m9IIZceFcuwvZ4GTRHKE3uQ6
mHZp3NnZJmkw/O+42NyQr8RfXZA4bZn73jhQ7x0/vYcWjLMwqYSXxwyla14m
wJhjlUXK75ohQt3lFLRy4NK7jFrPj6Xwy0ppLDMzBsZazPlaPbYMBhqcHCyv
Hgvnk06rvDj5FZwPfhmk758IanJvjVwuthG+DNOqCmWefSDjn8Mcvc2n87+B
ySBll21KWfBA02fdq9JmaD2ffHWVBBfG3pg/3qahgeSfHlua1m0v+Uz8U/l7
sksCZxSSfNKGIvM7x6elwIlHjQkZB171yffR+S36eTo/QcfjdPxLx7t0vpGO
j+n8IB0f0/knOl6m4z86nqLjLzp+peNBOj6j4z86nqLjMzr+pOMrOt6k43M6
vqfjPTreouMrOn6n4yc6/qLjJTr+pfMJdD6CjlfpeIuOr+j4iY7H6HiLzofQ
8RIdH9HxGx0f0fEVHQ/R8Q+9Xp+Od+h4iI5/6P2EdDxEx2d0vEXnG+j4ic43
0PETHR/R8RMdL9HxOJ2fofNBdHxFx1N0PoKOr+h8BB3/0PEOHY/R/iYdD9Hx
F52/oOM1+vxGOp9B58v+V7xF5zPoeIvOD9D5Azr+ouMzOt6n4zM6P0THa3T+
hI7X6PweHa/R+T06fqPzfXQ8R+dP6HiOjlfp/AcdD9LxH84n4vm9dPxHx4d0
vpCOD+n4nY4H6fiXzmfS8SQdH9LxIx0v0ufn0vHj/8rf0fkEOt6k82t0vpjO
N9P5FDr/S8erdHxNx7d0foCe36LzZ3R8S8fDdH6Ljo/p/C0dH9PxMJ2vovNT
dDxNx8t0PE3nf+j4l84v0PEwnY+i42M6P0XHy3R8TsfPdHxPx9N0PE/n6+n4
PiV+cnHRmBbib81Sb4qoe8zHgi83xHCdefZ1emVTbC4XXv0deXfjixZoNlR4
degTF7iS9usXG5WAs5LRuhOikXBEQqO/1MbPEDNg6tnk2Ah4arCjn6JKK8SJ
rEr5YpIDT9tEnh9Y2Apu2gUeSrz+e9EywOeZZxtYnNk5adWGHBif/Wu1Z1Yb
OG5U/CHNw4U+6X8H9q+Dua/3JT48nwp37Cc0OT+rg2be4DrBOhVGqj15eE6x
DkQL3hmcDEqFAXdVcoVD62DxjA2m61+lwsPCWbkTLLPhmFXhpFtR3hBo7jBw
XR4H7q8a6Xj5vjc4dv+qa/tdA+uNPr6e1JIMBsyJR8MG1YHkoIWciwNT4Hui
+PBTPP6u6W9fZl6fA8t0lt/PjWqFO34rhcZ05kCdRr4VY8rz31I8pR+/zQHX
yTa+DLcN/jamX9wRlAO2WluP1Tjx5ysXOJjG6HbwsensGRkr1VuI/6fDkTtl
EczH9ke2KOlZJULnpGnRIavc4cTGA1/HNaWDU6ya84A0dzj49HHAnsgw2Hwz
bsSvE3bwdr1hVPPcaJgzdeLT7i92MHiE6vfIM+/AfknZDNmiKzAnYc9awUex
8HfguoNuYVdAT+9Q4/FsD2hzlQka0J8X39auPaY38gVIfXjl8XytFTgu6nyk
tcQfEr9liE6edwkEZFd0CEkFwL0ZAyM27r4IlouePV3Fa/+53btE10+JBD/l
kqKLvPZvtMs8dehIJFx0H3tmdEMdLHokfdpiZgqMGp/RckWuHl7dHX97KC+u
2VAxefwVnjxQ3uPS7/QUUHL5m9rGk9d+36dTmpsCSdrfBA147eOtOPzZed9k
eFCrtGW4cB2Uis39uvBeMpSYm05LWtAOLY9XDJuzIw8MRY4bDVjbDlExQ4Il
d+ZB7p9Rsmvf1hL/t63eK+fTiTqCEwJTO5bKtcOC67pVxoZ5cPMq53eHYjuc
XHDScbdJHiTOSGx4cqQauvzdb9q9ToYZVuHDDTWrYZvSeZuJtsnQb4f3mDFh
lWCrlh0mvTYBZmsn+tx0/Qr5p3T3681PgCBX1xuvy4sh99GSyDr5SAjRHi88
dHIJZCe/tD5hFAmcUp0Ei8f89WzXjn48sCmUj+eGzogx8mkh/oPkxnGNQY18
3Oi1YUHhgCrofyNyl5dlIqg/048/sLcK6s911one4tnNocUu2VOKYJz06CIo
CYUC0V/pf6SK4Wyp8MQDyaEgfahd4FR8GbQuajPdfjkaLHRqs5wOlcMT36U6
cQei4YTApneVPPn+4PzA1K3vYI2KleUMo3JYdW/p5TmX3sHOW9HpSSnNcFHO
e+fNjiwo/JaddP14C8zgLjQd9DEL5sUIrHJKaAZb/+bm+k3Z0DLxmqnO4RYw
trX7fkcvGxjFhk73xBbYsTxcfttfLozsEvw05l0LzBp17MfLDi58jT5/4WNG
MRh2p3I+8/jldmK9mcXlYjijtnvuhIQIXj/3Fem4kQgJ1zL3Rha4gWWBQkT3
j3S4NnLoia61bvBrhnT8KsVY8AgSfZE/zR32q5565yITBFnnhH3/HHSHcenl
+t39WuC3o2mG3WwOrKxQ+pnHG1+Z8UF/9k3hQNfY0446gi2gyAkWE+X59XNF
vDYa3m2BqD+yKsMOcMBwhdluYRNX0LaoXL308RVYHtzUXrjSFUZNHPLu5fEr
8HvqTbvjMWFQE/D84eqXD+CyhfKbh0uj4UxEeJS60gMQFLA0fRbhAxMWCU8e
ssMOftw4+db+0HPI2Aatd17ZwfHb2zjCzk1w6s/vmrrObNg/IWTZIqsmyLn5
MbmgMhukhwks9WzKh++BQzWfvAyBGkmDy/IrP8CT5VEvs0+FQEj/ZeGD5n+C
m/tnmqVeiQa7kwr3FSZ8AqO1Owt8DaNh5f2bflMOtMLaY0Yiic9z4Hi+ZGva
1lbgPt6zZPz9HKhZmfJ8s3gRGf9lb5iMaBAvJvibWHR0v6gi0Ch3lc5aGQ5/
7SUXnDxZDA8zGkb+nB8Ocf565tk8+YTf444rPA8Hb7lXF1xPFUO7cUVTcUQ4
rJx+0GnQ+UBYenXhyi9eD0BvsuASlcE8P+zcjEOtKg9g2gtu4IhAV9CdvsR3
c+kVMHoE4y/auoLYu0nW56N58Z2D4srM641w/v6Nh8+jM+Hv1o5Da6c1waTf
Dyuf+2eCzu7E+4aK/HhEf+EWr6rlfFw37EBVTX/+eqlnFnp/8/bw8fPlZf1V
l9TBslhH1c8dKfAyqlkt/l0d6HzOTdnZkAJ1E+JlFZXrYbH7nPOFGmlQXnIt
7vS6eni456nnQK00OJIuM07sdQucyXHyjF/KBaH4O4L3eP6ll+iwI7oLuNAS
K1Sx1bUFTL7NTtEp4sIP7dRb91+1wO25r33flnHh4LgxjVffFJP1DgOOCJTW
d/PxM+aefEV5JdT+fCvvMS4Rvg1t3q36oRLK4z3+1osmgnjwJou/1jw/+jtn
b0KWD6yPmN35qCQdXiS8rlaT8YXlI2Z63Qz8BiXhBec2DMkCxifKKz32GyjU
3s0WHJsF26Z0Kgzf0gRKvmV72wOyocyt32nXXU2wTM+mw+ttNhyf9S32yr58
ULnS/W2YUwi4Dvuyz6YpDwSutM9eODkENkVKRdhqfINai/6HRV9mwr6KEX6W
877ByYd24mG3MsHgu6Khxvg22JN/af/aXbmwc4yH1G6pNnB4cCTu+O5caHct
rr7xtwCsA6Q9xBaFgaeoTPYesUKIvqmqOmNDGCjJ/30oEFgDS0eOfmr9IwWM
Cx/oZnvWgOq9kDvLv6XA7D2nBkRJ1oHEAk3uUP006EodoyY1vg52zHO1cj+Q
Bpb32ndUCbfAeun0VxOXcSHyceKSIF5/L9GzsfjMax+dDcf/OG77CqYWhxLe
30uEc3esJN7s/QoGA9yKJFx48agfRBkLVMJru4/C3SoJcKRR/9QFzhdQ9BK9
7zsgAa4IFqwedaQWQjJ2T/zB80+2MT5qWw1q4bK1meNcr1SYnLf1QXtzK5RW
l4yQUs4F93iDfTMqW0H70sYHMrK5oBR7plnVvA2mRYXe+RGcC2kJe19anGkD
jYiVJlpvc6Hy6bWpj73bIK3LTGVfQy7sG7BW/KB7G8zZ//nZyK+5cEPtlNXk
6zVQG6/6dA5vfJ0yUHBHyZ0aGK90PU2IN75+B6cR3Hm1UBrrHPPzXCqcE069
+Yyphdu6mzeK3UyF8S8vhpcotIGMndY1xUO58KT+gInBlDaQO9nFeB7IhfQL
CxN8bZth3Od1W1MSODx/8pGs5s1mmJSm+Fw7kgMWH+V22XjHw5cm/4cDitxg
irruL5/J8WCd6lKWv9kNfpqMkhpkXAk/V9e9M3RLgJ8HxEM+HKkEv9dZKqoe
CbB+8eaT5aMKYLBAbq6URihcb8wPbxhRAJfibPJHrwmFX1/Czwg9aYTso3Oi
kn5kway0XR3Glo0w45uDY+mnLDBoPVUouTkNQiq9fa14/l738KKnG8TTQG3+
J8cZmt4wJCXFf7NQKZxXlZypp/UOhkec+nJsZinY1zSl/Tn3DiJDdNcb8/rn
gcCXj8x4/bOosDXHi9c/67QUjHJ5/VN4Y/rFgzotoJa75b7SAy4EqWavfr2z
BW6J7Ik5dp8LrsLxpd+9m8Fx/vLdl0o5EKOrvz07sBmuHVnvcK+KA0NC85Rl
pzZD4+W3Ihd3cMD8UuP9G3OaoVDbw3Lyfg4srulQbdDijY+/dUbes+bCpUTL
4V5bW2Bz69gLW2244GLtfbT4az5wno9qujwxFFRSH4bXtuZD+PPz9yIVQ6G9
fo3LgT0fQCmk6Znn03AY124xbITFB/j02HuJV3g4JJ1YfQMGf4BByzZ3260K
h/VqBaeLOguhxVWgfszCcDjWuk/lHI8f86N8xyYJJoJ0keRVhsePFV/uMC7N
CXCvtKDtoVI91OxxFJvDs+/h6Y7DC9fWQz0TOUue53/Pfu/4dbZgHZgdX77q
vEoqdB67tGj20zrIv3xUbYxoKrROXH9bh+f/BFR+eHZ8SgLsE58bMZPn/4hM
0be5KpsAu5ydFXaObIUfhndWOKlyIcBpSt276a1gcEn6uuFCLgTfcVfI29gI
SZUaQ1sqMmFIxqqxJYOaoE5gy8nSxExQkbn5UoPnH6UUFI4s2ZoIOyvExVwN
qqDlc9dC59WJ8MfC5PIw3vgmMkz8vt3VEKiXPzn/AW98q7ANHO29NwQWF0mE
1w5+D4viN8RuWeQOHs5FV2OM38OfE+/zpZ+4w7BNdiLXFgTBjcWRv2cZ2sFy
IyEB5UdBvPInHcpIsYNRT/yzsnn+acXcYac2bcmDC2OO3JDg+adN06ZUNGzO
g98OWXd+P2uD+asFlH7PzAGNe5Gqq5PaYGrJl4QaxRzgausWqb1vg/FlSV3x
arlweFBUv4kDeP7qR2Wji6q5YCiVNcSwuQ0mvP4tlOSSC6O35znbT+K939pm
Wz0vnlRRcZQYI8iLZ7mKzmutcuG1w01NE954WO0nkHr1Zi58tvj5a4dkOzyx
bNa0+pkL1lMWfHys0g5Odk7zPn7PBU/z7h9GJ3jv37nBeFdWDtw/bSuqldcG
Lp+dZNI5ObDo0pn8dp5+Fq8UHjquzIMQwdvha3n6DdppoFe8LA9kpylbDaxo
g/ZlKllaR3Lgl2J+ZqAI73t33cxW8uJBTVl9+7xjbdB2+ffMNS9yoPJnZR0n
sw0mhW88aemQA6+NHrmP3NtE4u2xn4NEF71qJlhlsnaYMU8fo8UFS2sH5cHG
+a6DGnjju/XOIdAgyBvvhz973L9fOyRb6Q6tMOLpe7bfzsnT2mHRyFmHQwxz
/6f/9v/b/vfKS/ynrzAErYFF7U+6N0NOrKbNRcvzAJfs9N7XacLJWTY/PWXv
wbwK0f2TLDcCIxC2bZDgI6hL7BwSobQWliX5PN/F80ep+7bg2pVFElX3HsHq
pfVFE+/0o+63FWSa1LUFD59+Bc9zd4j6aA4h+81mn/9x4QqIUecdDGUefBXm
fNgSCD/NvX5e2D24z/mG6pqvwwa6BoKQ91rJKGv++Ys21eo6Zv2H9DnP++6Z
AI9rR8JBeKGPz8u6kX3O5z54s+14udU7UAk0Vf2xjn8+1aOQI6fSVfj3fd46
9ON60Gz+fZ/nZWLER4ya3Oc8KXel5Mk7bfnna/8QPVPxcgz/vAOJCYVbVbsS
gLNYPEclmH//4vcLAaaV0/n3IRqsbghXkuPjxcUGq9+L8bH40cJ8bQE+xv1x
uhKp85puZoLc+m0JBaf45zFZFF3YO3sp/77AnlmNWQwwQ64763HA2qmsKmUR
/z66QKPlxRW6o5gJhnWnvySGgtRCbUOLieKMRpEv06SXANHVQ292XZ/CyGoO
e5rYlQRtsSv/jhsszax1bKvhuGdA5JoPXp/9ePGLTHDo9BoudJ+sGbmMh38W
N8++Xs2FG5xBjT94/uihTYcVllRxYcT0aUmBvi1QFPdFqrmSCyrvHaIz5WUZ
2e7aTWr78+DWzEeT+u+uhHUbnjSLxVwDz0OLuTbl5ZCqW29+uuEGzz674CG3
HBS8n/r+MoyAXzdkOpf6l4OAso6FxtEIMBt69Njto+XgtzBvZNfpCKi/G15R
olgOeis2Bry0jIDdeY8nCixthQz1K9LFp3LgSlhwQMqSVpjpv/mNKw+nXqgp
2snDjddb02fxsJj6BMH6xa3gt9uuNuZkDsz7JrXtqWcbJB1+JfNf/ujInMrp
Gx34+/d+Owblym2qB3ONyd/FHVMhlTPJxlOtHu6NElqw1CEV9q7o+vWGN34l
eLyPuGufCl1B+YdCx9bDrzlXVcsfpMLZ+DdT7DbXw6O6S1+LnFPhlpdz+eI1
9TDruvJqpyep8Pr0Se1bKvWwJ/rukDFPU4GzUiOzSLYe7rd/y3R4lgoN2dMG
VfN48H4/w64vVeYgPGfkrtdxmSAzLf6uUcI58Pn593R1MRf0OmSrf7SfBpMx
T8zFhXMg26g+/rC+N6zNc5P1aOfC/ijNgN17veHPev+LydlcuGcw/q7aIW8I
eB9dXOfNhYiBVYenHfOGS+v2mh/6UQJesbHvpfoHwIzI+uTKiyWwdfmZ6Te6
/MHqyd9jn+4Ug9eK6k7tBn8ovSA2quBJEfxV256m8MkfJN7u03o/ux5e7uP8
0ZBKgZsu886XSddDxIplfuUTUsDsYEH2uN91UC4wWnzAuBTQVT4jt+1DHQyK
aTZXGJ0CDi+XOFudsocchTrpxV7mIDLPyttRwx5GnNQPqb9iDn83jUmN7GcP
N4evK9s01xwajXyqPwfZQb+X84YGVp6AEQEa1YozuHBqXFP7Uu0AyA9SrzKQ
4YJo1ZElmzYEwJ63SpXZQ7mw8sKhlOPzA6D+3ZiKZR28fia5f5vdpAAY9zZ9
y4RLXNjVUHR4sE8APLX5uKbuGhfi1pae+uQfAJ6z/s65zvPL7glV37WKCgD/
fOmJsk+4sDum6eXMlADwaznnWpDRClo+qgGHhuTCjDMfdpelt8KyPX/eHx2c
C/JWtgVtqa0gkWZRe1o4F6QerNMQTmmF2rk/RSwH5MKizvlLlpm2wXdX2+3/
5QOFzZ7PSanmz2+5SLnMOLybfz9a9SODxYu0+djiTt22wnV8PNzyhOlJVT4+
c8Xzg81s/n1n3o8MNw5S5GOR5Or5/cbzsdl3w0mdQ/n4qoz6tWu3uHDWUKVs
r64HSPvr69zewgX5nd1pjls8YEemw+AnXzlwYand834rPMC2LjPM5zQH8qWn
njKe5QFT+6u06GRzQN+1cGp4hAs45whalm7hwKwHe82XVrjAgC9LVye/zYKc
E2eGr5z9FPImmd/ipmeA+ba7LzdefgqanOHy8y53wiLFV0X5Ra/hrnnuqOJA
IebhgiCpJWODAHTvjnpsNpCRm+cr+TM6GO494MiY75ZkHParrL/yjoe3OG/4
MVuYeTNFX+XUxAhQXX1jcomtKCNptPrIEM94iBz+0NLBUZJpTtKL/yvwHga+
iXkS8kyeyZhjLCMwIglKsvWt3KeNYa6uS7kUrZgOZ5vebj5YOoaJSVkB8DMd
zhm+TP3eOpHJ21E3+LNHOuRvT5Z5UTSRubcns/laMgdq1KbGTD8swwyojZ47
2i4X/lzenvPmujxjb1+qerqUFx991RVb/ymPzD8kX5Q5r29c1Oc83Bujr9/S
/n4NTLcL2h2XSYCFvwpeaA09AIP2tmzbPPEtOC9T0euIOgYpD6eP+v0oCkD8
3Yv1lkfA469GukJKHFyPzJ7XMewQ7EkpPKFemkDmwwtvj1UaeSyVOr8/HqS+
KnF2fDeHCzP76w4tqIPU4wbqPxptwXIqp1orqxO22TnMqWlzhTzDNa+imSay
n1f3vbOq2JA2gjcodWne2Purz3n3pXW/96818IKfk1RHnNcZyNT8LEx79PU1
qO+Nkd9YKMDcrq9OWjTeD5xk/a/ezP8DBXOXJOen+INwqpbBQ2sRJqFLt0ZS
xg/mp0W47NUYwdTMv/ihtp5/X9SJgnFa2+2C4dOlynt58/j3R+knit0cmM3f
/5v+WaWksmU4c1exO/TOi0gYvDvQfY75WKYrL+hGQ+NbUFm7oeCD9TByXoKR
WUfzQRlJZvWyY3O3f40A87FOu9cWijLCinJPVcfHgMXf6WJQMJyZqnJN1zkl
Fgb9eLb12S0ppjQ26kKsdAxIOyhwL7sMI+cpFGfLuW06zj8PV+baDZvrC0b0
2X/cX3ZzirtAPOgwRvLcx2MYpdG6d2xNYmH8s9shl79IMBO6262WGCTBnKFX
zcokJjIbfbyVAo+8h/Jzna6bR8oxQ3f7lt88+R4qGrsk4raOYZ41Z5bon0oC
h78al2feGckM61pyO2VZCpSst/o8SGUc82XYYq1O5RR49tR504qkScyYCaV7
Rh9OgjWtbzerTFNgdrdY2zp7pUDb26fD/Cby71Ny2dhxcrugHANWyh/SDVNh
whwrWfVAReZ+aeWVZHVe/1GumDBZVZ5xPnDpTvGBDLJ/2vbc5IzAgiyYaiQX
MX8L//6iqN9LVtz4MImxUW8I9RybDS33kmpjn/P3UyeeLTvvdnYq81Jjc9Sm
/dnw8Xt8meNFJSZB4f2sKcrZ0PB7W9MxTzmm6Hxw+dzJOXDxsp7bH64icyt0
pdG9Fi7Y+wWtHe0wiQm+/c12q0cOfHmSNrqyvu9+bO2FYxYYqufCiNEaGcuG
zWR+Xzug8NCUfz9xR8y3iQ8qc8EjdtZwz8v88wC35Fm+2tC/mjqPoh7W243/
mT/8NpnfzuVcDDduvw+vL127ID7kC1m/t3qHdaE385ngJQkbBLU+fYODJisu
hQc6Q5zZ4ls1S2rA+pVt2rH5zjDw06mLY/r/IP7/r+vDz+yfpsQI6h/h2M3O
g2XBEV27bAWYDyK1prqegaQ/es3apHX14xuwm9vecOHJb3KecHuT5sJbe2up
+0TqYfiqlqI1Bk6kv3kJLxId9cYNDpk969q+ql+f+4IsatNqX6/0hCPTvjvO
n99N39fLaBn0f2LSLxBuWWWGJMUP67N/vVkqLcZtvT/cGT3wemjwiD7nCbuN
6f666lY0fJh1/OWLt33PD7w5y/Lt/tXviP89c8rzs7s3xYKR5Jr1pxz73u9z
3NaWG3omDjRjlP8sOtf3fMDg9ye2zfvwHsraBizf9qrv/vnKjj+Jw7Lfw+o9
Mhbbk0f3ub+n4XyBsIsc//4eDqPwoHFeMvhVqQ+5cbTveYLPViVkeVcnkf7y
brhb5NDNqbCiM3Vndyv/PORnEgt16p70va9n6frhJ0x/pMIUkVnH6lL77sff
uKpa6faHNJDODDsZld73Pp7ocZHTo4T59/E8WZqyxGtABrw8l/m7U6Xv/Tzj
dh4ucJ2TScp/sH3V2/ce2eAUaXap+Vff+3gMKq9u2TWMA8JtckfPDet7v9hN
hYKxT+5mk/NNtMbYflpZxgGHhQuMvBfz7+O5IP1uoP0A/n08Avs+SZ+Zyt+v
v3F3rP9wVcU++/U9h+3td+khF7abHXrQpT6pz308+neKxhvt4vGDovzBU1v5
9/MoHnaT9z48o8/+/fs3458ck84FccXX8ldP8fxn4w3+g8vuwL281o+cJR9A
5eFOd/35NqAu+r31iVkV/HBsgo5oe0iw27r5S2AlHGgKGrh8rB10NlnkXp7X
976evO5jhX6vH8CQwXbq3qtaCD/U+t6Lcp39A6J/fP/hO9EbnH9Hec6aPoip
013uHTk9EroCdasUSwYx2oc90zU7I8F//3XtU45CzNp7WTsk+gVBU4n02gv5
Q5hTb2NmazUFQHvxkbDvn/oziofCXifvD4VNKgnC/ZSHMLtbjxoUzA2HtSlp
LjtT2uCzkcCiJ2YvYd+QsAyxox3w5YfeFcs4D/jx6G/Lns/tsCZk7b4IFXcY
Gq8aEKQrxMy7W+JyOsYDtM7IOtzdNZBZUnYcauJfws3NSkYnctqJv9DT/h1g
/qbp58KpLwjW+CkTWfbeC75s9J6davSD7De9ZHp7uNHt3wQXn3OcP8d8MOOi
oyPxrdEL+tsd3f4rUJA5+lF7g+K6AMJfZ/t7DRCp9ofDM5b94Zj2I+dlDPHZ
27AmpLvPfWr5xnlDde68hnGOspNtlP5C6tDVZ9t++sOBI7WSRsL9Sb4Ezxtp
3bv498HfgbCm9sOotxwxRrr96odBCoGgKHQ91HsO/37TMmMvV/tKMWZEjKns
w5owKPjh4KcnOoLxd/474fSrKLj1ap13Ze0o5uoXmLhzWxiovNp9XMNImPgn
aSUK28MvD2HuHt4JV4qiYWrt1Mh9s0UYq8W7n338FQunV9gdnqI3inG8P3OP
ckU0RC3XE9IOFWfMRFNyBNfEQejQzgPTu8czv86ZaL+4EQVP/ZSnHSgcx9xq
V28f2xgHynmMwzybScxnt50FoSWxcPtMSl4cdxJzIWL+qK9ycfADRK5s4Mgw
3zf759spxEPhx4mHjJeKM6Xbb38acDUBHtW8kIGToxgJTmXiqqgESFWc/J3R
lGZkOWK3deoTYH/wsvtndk1mJAbu9dKKT4Ds/YduG04eyYSfCPqzJzEZnA6r
OF0vHctEZma9/FyaDF7V/ZO2qoxmRCYf+7RpPC/eOfrTc2ulDGP96suA6ppk
kHszYXp05RSm6ul3sVadZBAMdPk4Q2IC83YPY9BWlwq+uzw9Vk2azATOaV2Y
XJUGcw6nHos4Oo158uHK2Xk8vnxjbeehd2Mc8/zWuUTIz4QZ4fGgMHs8o2o4
n7lvmwULxVLOHs6RZoa1TnwioJ8Jc1IyDwvNmsAIT3846WRgNqw3XZDwTWcK
I2WkA8JfM0FVLakqrWM6c/PFjA16DzLgqW36tKXbZjCySxYN8GvLgP4R3Qsy
F89k9tedePa0LBPaFmyMt0yRYZQsdFsijDgwP/r5/WvjZJmB6eHh2cEc2Hzb
S365owLz2nj1kiJTDvyI2ZCmHa/MrMvixn1S5cD0+kotAUtpZluD/PkFoVy4
+m5zZ0uYLKM+6FdW1youeLqsbo8Insakqau/cZvNBZN1+x/EGkxjbh/rep71
iwO7is3Ovnk+iZF8+cCt0y8HPkyONz3rOJn5rXNJV6xfLnyYWZ44rWk6455t
ZnsrIIecdzDO6bHHLNF7YGZg+jtcqJP0H/E1N7wdbpdBW9ZrMyPvB7DV8aPt
0gkVsCUhoUBOzx5kLEw6hZQ7YVWTU3Hu3FcwjfvHZoVmE0y28Luxq94ePuza
pWjr2woB4VME1Te5Qnm9lpTg1g4Yt3qEu2KTC1zXkrhiu6YbgkNdw1Y18++f
iV8iLWdh/AIK+tvbmkoJMcb717hL13rDyKPdQpcfDWJmGy+aOtyUf/9M98o1
TxMy/UHh7qxHJ8SHM+smT/V3PhoEXrkCL1M0BjHXR/nOv3s2lDz/Ysy3B5rz
w2CzdvrvyUmijMp4jZqlh0PJ+V2D3kXFrKsKBQndgAbOiGGMX6xu9Nhvb8nv
HwmMmSe3IQpcIk6rh3BHMKrSqx+a60cDpzBuQ+e0YYxO0QkPrdAYmJ4xYvPG
ERKM0tUt/rfNEmGbjr7sh3BJpnHML3NljWQYxYjK7ODw7wdEf0MhbK9VTWcw
bLK9tq2gln8/8upbk8dUHRRjRA+HbFB++B5sDh7UdHXm35f86Ljc04sBksyk
L51RumPew6AhyTfUro9kzhacro0sfQ+uupNnOTXy70dOP2fg47RIisn6OVX1
wm8OvE5N6hj0RoqJ2WxekqHKBUGLVZZlwnVkffSI+V8ak6f9hOyBdzLzFH1A
Ie6zZUTxT5g44NDbR50+MHX9Pa3alj+gx2VKbr3wgQ3qanUWot0wJVxUyfqV
L9wSit1eqF4Ke89bHxQ/ewkifnh05z7Phfz+ly5V21yH6TeLnWu68yF2vulR
fWErGOa1Mk/sXTE4PNm2VuMv/36aKaHrBh1Pt4Tlin+aRDeXk/XgCx78CtzZ
v5DEx2jf84KHzbsYfB6qlwo6NtmWU/59CcQbr/tWbmdNcGuStML+M5Yknq8S
fRpi3H4ZlN/EJv8dkUu+t6+fWe4Q/WKCXV7XT61ZlE2wTbDJ6cN708DdN3FN
VcMFsBzYeU/VOQVmhKx9YXbiHLTFFXVvO5kFH9qaNs+aaQI73ZWunY5Ng7Lj
1vsL6/j34XyqcrnnKsLHRoZtnI3xhyA6ff+8pyK8OHH9oob+d0zh3Zs2w1DN
NAgWPZITds0Ukp8/1KzeGkPOwxds2liudzAdCh7PuyvucwRyJ+SZ9l9vDbqK
hqMW3t4CgQpvpTmDn8Hn4l96FuFbQatw2xLndn+oWiLd/FdUGyZvtHmX+zwY
VG+57VQ/tBWuF60aUpzxEqQUpN/vLDWAiQ1vNrv/9gZ15Zlrfv3Wh5c3oXnc
mwDY+m1IxrULO0j5Z519rC+duA+Wut5rPfIjrM99PUtXnhs+Ve8+iV/6L2xs
PKr/AAxzpb7bhzf2ub8nxHKuyiG95+CRNXyASOZvYq9N6uvdlI43w3fO0duO
UW7EHxv3ofvh/IcesGdmh2fwPP59PgYx/Tr153X0ud/H/caLZ1Kiz0DK5uz1
hds7+9z3069xpON5z6ckPjtpbbfaQd8Nsp+3WqmF/yXv+1QpsPtX4wASL6G/
MUs76/rBqz4k/ooVqn94/N4b0N+kmVooPYj4NwPOZQhrLxFkHo0bFBN6NRhm
z9vVP9myP7nPxdw3vnVGoEjf89Eicy2+F/Hvfz+rp9FycmIIOP768OB5inCf
+4JGPpN4P9UkkuA8UelVAmGxIHc94PWFJB7/PM77PMCYf1/Ql+Xt73eUxYL5
UvhWvbjv/UGK+4IDd3nGwGMmyCbypQSz53ue/8vJ/PuDdl44vF7+dzycLHfL
FuLy7xfymTsvt9uSf58Qnv8WUH/CQLk4EdRz+j239a2ATWrFoWqbbOHwpPoa
ia1VMDNq5ieJpjtwwGv6RNPMhj73DWV+Dfuace8mHF8u/Hl++1ciTzbc8Hdt
0meCKx5t4Rhvr6PuX6sCAVmT6g2ed4k9zVRdlmjc/hzG7435vSqpidjfhxVl
SmPbvxPc7dSmeLW8BYSvzb5h2eAI5oOPDPSdWQ1TZ554pzLpMTiIyg96M6kO
TrufXtFv2SNiX2c15AOE1jhC6Vj/DnfzX2ACbmd1Gp1A4tQO1TbddoCoQyun
xNiDvkXmif47OuBjBSf9ZYwDnJ1w1kSoln/f8GvtN4lTxgsyGhFGw8x2B8PY
m9sqthzsx1xekfr59sMgkNM6ek30XX9ynuV22UTfsV/7McemphbUVXvCbKdj
5RHcbgirePrnuP5r8JcfmKU37Te87v9eITTUD07um3i0eVPf+44MYrarPw/z
h8ah4l90cipJfgLr06lQ6cXCqQ8Jlqlbv08owQkG6cW+ytCogfNed17YnHWA
jXfsfF8YVZPzZp+PKBY5cruB4HmHzb4byPHv71TRuBmrw+3H7Lmv4L5wqhvJ
P9xwjFYdvycEsgKz646Y8u9LmnRf03Ggx0jmztFfcWcHh5L5xgfP1JokV4dB
7ALpgZlCEoxJ3oq3Hvvfgt3Vudo3X/DvR/Ku1VVOnca/H6l4jtvCMrnxzIrd
2/c+GZ4KV32+LsmN7ns/UqHMs1d2K1IgfNtqHb9PX6D70oWQ9v0O0Lp4X/Rp
x0qQ3JF0QryfHXj9DXzuN4d/f+yCTRpDv8jx77MKcWvbe43bAva5I54LTX1A
zjFDvkZ+7X52ZdAS44cwaGX1QTepdjgh7yBTXvMYxJfPO2h2ZVif+5GeKh5o
+bTlLRwOSzri2TaEnI/vMoYbviKCf18S9s8/12IdUlvfQmrkg8v93fj3J1XK
Pp+SCiOZ0JMxz2L2RZHzfOGot8bp7ijyvc9PQlxOlbyDw1E2ujka/PuTRLJG
XLp4fjSzNjv+/Ys9MTDHZvyJNEH+/Ukz9RSlpN6OZexT5UpyleJI/c5dp38/
4nYc+d5zObOjW4YmgNCcyAEvv49nGr8rCe5/zb//aNE9hYormvEQ+9TkyYR9
/PuQ/CyjZkrOkWK6jXPdhWYkwL67p3zUffj3H+F5rA7rc78Mf5IAnXl7Ipw/
8u9DYhi7iQn7RzPP3p4LOeuUCEp5gWp6v/n3H+2aqGo0deVYxjZ94QPBLclw
NV0+YuJu/n3Ra1v7qc1dMIkZXDhaPMybf7+R6TbPQ7X6iZD8Jt1j9nv+fUdu
yx3Pbdwmw4R/cnU/Py0ZRPILpNdZ8e83wvNZO1d+ePPAJxl0Eu4GaJ+QZRqv
Hgza/op/v1H7y4Zbl01SwKPowhspMf59R/YqOTVHAscy9+03xg2pzgBw+eqe
OoB/vxHWf5DS9mDvu5mgOP6B7dlJysyEHXmuxwr49xvZXv6qkiLOv99IX2vh
KdkdXBj2zP2s3vRJzK3guN97Z+YQfY85KTGjf3Phcf8rGr8q+t5vdPyb2mcZ
hSTI2DXPulZKiuS7UsYk1O3+IsW8nTrw8eUrGfCoy655qJ0Myf+OvuVz5HWe
DCPQtVB8xhD+fUU6t4fbl6y+AZfi3er3aFfBhgTduIJ31+Hhmk0uWp8+gaea
yITwQBsY/rPtsc2pOjC5YrVKqMwFkt4kj3c9rsSkBXRwnA349xM5LXQYZpRu
BscuHBiip8EFjZvfTwbamMII0wH/p7izjIqqbeM9KaECIpJS0qGohIR4DQgC
goI0IgIiWISUIAoCSqpIKQJSUtIhXdIdMwwNonTXACKKcPas9T7Ms14/nHXW
++F8/K2Zgdn73nvf13Xv2b+/TmQIZr9+FYkskz+a0gb55zYmvTsfA9dsXV+j
XRvYPqt+w+/mADvn5xkKSpj3t8/whyHDLDUzqvLxm8fHqdpAV+unlL4V3/72
hb+gPtGWw4sKMhOkuMSEhTekpjUMtHyo3afFM3V2WFAYtPzq9bNg37dtOBel
lFVUus91L7ypvCNr/spf5PAeObyeY7nPHuZVNmMpt/bZDUVm/FTICPa+Z6mq
NddB2vpcZorDIwhykRimt2kE6SmFHq8aJ6Bkbui1/FYPW2+tX1RKOsLcRZ7+
4MBlOOM+RoXlfAPzeedDbwbuwuww/a2jXHHQ0JHfsexDjlK4OPv6gSIhv0iP
7arsWlkp/HZmVXY9TsgzEic1ieY5QMgz+ud6lSVpzh66UgMdHIzviswYUUJF
1OaTooQ8I8G5RkvzghZQ98v9fLr+73yjqkOkS5rWzTCXHThO8YCQb/RPHk2g
32LgxE4rXLjOoTVOwo5a5dpVp2huhTKqcGkv/+N/5SHx859m5LxJyD9KHeu9
SGXTBsbTo2oxOoQ8pK/3x7PPRfPur2euUk6oVUWdRHWqbL6Mx2AhUL9RN/jg
qb/yj2hTZpx2vLBgttVsUe00Bza064znlgMAk7ESzbb3Gya14XWp34f9etYK
ta4s+CMapEZHCrQmCXmK//hXr1gZ+09c6oG9sp2pcS1C/lHhuer50Cd7YPLl
G+mlkjw4VSLj+kyBkGekW2fUjzvFjJq/Hx6u+LEEpIvVeg3+MO7f7/tn/5Uc
CB1XcCzZX38X4qWQPBFbCjcaNTW5bJj/ym/5VjaFYn1cDseEp5UedxPyjiiV
Xy1+KGVFUdFMM5vJEfKOhCxWbVNWy+HJrrxrwWVC/tGmVNRTqRGW/fltP39P
amFirKkMsNM3vAaU2FFjukpmddjK/fE7uvdKW52tEsTu8n4eFCTkVbvu2XNN
sjLv37/4Z/sGZlbSRA8Q8pF+Jxn22LVUw9S4vPT6WUJe0ggu/fkTNCeq1Vjm
BEayGl7qbtI0enP9lZfkw8mL8nT4sv99yn5UuBoz1MHg7mMLlbC/85IinFZF
dvgboZxBQ/BGKPNfeUl95FjKZp5GaMWYDFtvE/Kn9xQVLxgp8qMS6fh1tQr/
lY/koXTfX7sJhse46LZ8BVBGHjc+tJAQ8pGuVNp+vvmsGRzIvPa0U//OS5pz
vizpLN6yvz/SdmxJGC+1Q2mSleshU0J+ktj9u7jiGkJ+koelpdBD1AkUhjET
Jz2AhjN6rgrC8YT8pNKlZJ0YMR5Ub6I/Y1IhIR/JzPBkMt8yBt5pnl5avP93
XtLt+9uH6eO7IetX/7AFsm8PmFfTZ4US8pEwgwv06jvd8JXfLuKrJCEvyZnJ
ONkOzfNXfsOvpM2+QOdueHZLyqpn/gCKeOy0Ar8uIW/5nW58W8DPJAhQZ70U
R0qKMmoo/cFvlQb6a6DrvXsQ9bnxYLW7fyZYSEk6fwykQ9HEL43Zj+SCpK2M
KgnTUdSR6xxMpo6fYXahoJUo/CiK07PxQB9PGZw7JYseOsmN4tgMLiL1rgFi
+jKqVgdkPn94uJ2hqxZORlsnjxXzoPgVk0/j2OqB+DbxDsMfXpR6Xzpj5t0G
2AiM6NYK4UH5F/jEXxJohe5c/pZzyHwiGI7qvILMU7eT7J8e7OBDfe8mSrLh
64AP0j6TU+9FUW8TiCxNZDrgj0Da4SYzLlRILeWQ/flOUGK8UGyvxY0qdG30
SP3RBd3jpujWeX7UuZX7r1OTOkFZyCXZiVIQpdss5757BA3W8xIen0dOorQG
Qs6O+XXCLYb6AQp+IdRMaXNOMQYNG00LdQ92kflOpPAiwxVCvo7pExuqs+xF
oNU6OM8uNw71Mg81hp+HA5pVmb3ZeBDeddSLn6RxAwUWicqQzG4oKhVOunjV
CwRwusbMun2g7EAiSrNCyM/JbXkb1fjGFYYaSy5f2Ojdz7MYWP3GYFjXte8L
OFHw4qnjU1oU7g+fwG5xDSi1hofYkrGjNg8oOHtqocGL02WGmIgbtZKpy9HK
hoV2uc32R22CKIEvM5b3tFogVGLy7B6tMCrGJIhlN6YVnu3e1fHXFUEJVpqN
hiy0QU3gNonpRWGUctLrqxFI3XBz3nqqyF0UtXVR+dvrI1iw1TNTxQ5VQ7Rk
wFei4tdAehU7xWBRDcQOXNevfn0Nxh4OCqdnqqB9/qrOnFAQ5JpdKJairIQI
C3eNF4+CgNk39VYvxdB/5fP2QRsfxxWMigvswmPXJzfYUJr6TJJGdZ3Qyef0
SXZhGi43U5sk6vhBv9hU7DszERSZtg3ny2vd8PEwZuHDg0S4LJwqRbLuCdRS
vlHdFK37+0+b5OkD35QacLq8Ff3VLBwO+HO7J7jUgHcbr4CJUTjYhlnH97PU
QFzHM/kbl8OhL7a47nBFNVR0DekYy4VD1423HEmuVaBdMqPNM+gFqkUFO5oW
yXCcyctX570PCNExnUu4mQwXngbwRIb5wHzdcMZFvWTwuhRlwhvoAxku8VzT
6slQT5cRkePpAyOlktjLcrlw6M1jWxGJZ/B4uXediphi//e4/6wHl8TdvGmZ
kLbfv/TvfZt5tJUBIrtvKyUnqP7Kh3DLjxawdMoDAfenXt9yCPkQMcw2V8lk
CfkQ/9TnBmsYywGa2v16vGlas+9tUi2EDb/Jv7rG/Vd+xDGDY+xrt+uhju10
IXEHIT8isuCwtlUBz1/5KS3s3LvO5a371zNV4wXf58/aYIDUQpoRRciTUJn/
ybowzvdf+Vj8KP40yxezbR1gQaLvYZAg+lfeROUbYZ7hlg4gN3woGQqEvC8i
EUmj73bc/+Xb50axvpDsxXs2/plPzB/N6nEFd8L4jDy9oy0hj0LrWVJHKx0h
j+Ifnz0vpoEzwJaQT9HCxuQ5b9YJPFY93gcM/s6rGJ4erqlCo6HDV1RS4IPg
vu9f05OSa+E0Ia/iZnH4I957w+DBtfZJlb8c/Ns+hxxbnIc0cqpY/POvh4hc
JJ4h3LzhJIt/3vWf9SZb7TOPK8LCgU7nF4u598S+r/klHTfzB7JuCI+T+E2W
lgX5vUsNsmsYCOWsVRVvzII9Uz9u11YMBHF+nXhDlA1X5rncCxMxEMj502NR
Nns/z2I0yK7yjGsPmImrjBrszIJ14Vyce2YTtGBpQxvObYCuuDsZ/vnZsQKb
CsV/PS/7lrakOkx8A46MU92VRj4/JP5QnRbhmHo76VyXHrj99LvZgbMbUBQg
Nxj3qAeW67Ue/Tm9AV1XKdzeOPeAHa/5ZIz1DEzYOniHZzXB64HenwMp0yBU
wJ77iLseTmXsbvFFTINm4hXJkhP18FBL68cXl2n44JUFcnz1kL+ZsGFkOA1z
pocvVwnUg8QQPTWZ7TBsc48LlxWUQ6f5nH4osr+fp2RaS+SWA7H9u9Vss2F4
l6OwKphaDhKeFwPbDIYhraTPkT22HJjIjsdo9xKer61QiQt3xhL4QFfy8/Iu
AmvpZzoQt/+dNzFRL/xYTaAYymPzfxlxMuyvp/7TL1AMzqYPtn2GPRU/owuK
hLwJ2vc+yoXP+f7r9yB8KNvQJB65kAZ4cnF2jP0XIX/in/ps2bmI/REtIW+C
S7yza/kzGlJUyIbyVAj5E6minR168HfehOqOtYI70q92npfKeNZGyJ/4Z/1W
Qr3/NIdaNohaKLk8jafdX0/+p97yc334YlumDkZ+X8hiNyXkSbh7l6fXrPL/
V56EAOrypXf0DHNNQLxwcytT8O98iXMyxqJqGS3QxjU3U3eUkDchbMzB8MlE
8L/yhgRRQ+2R7rbHWvbzGRRmWWzzFFphSPkYPZ3n3/kTTyxyCqMutoMSVyfW
0ICQP5HNJm2btfB3/kSvoordWaRv1l2Y1lhNIuRN/FPfHbmzeCjiXSdcIjqo
fJv67/wJh2Hms7d1uyBwWU/Awufv/AlfIQvBni1C/sRWL93CszuE/Amy6Diz
o+IYkG6ve6XaQsibqE80IIn49nfexM5CbfvqFmY/bwI7xqk8W9EN7q/Fw3I2
CXkT19MTcSkn++FXqOauNFcgeKoyuWRzDQGR9KPfRPL+YJ3BV1sgQ/DTDJ57
KikuSeACLaHjp0UJ/OZeH8lJHgK//8j1ILhlGlxMVBquhTaAUZx6sUndNNws
/dDbjRzPI+8TNH8WTYOIvcqPnuAGMA3dmglNn4afgjjG/jcNYJV6cTlWfBhk
47iqLuYVAaM7hSkRzzCYvvdie59fBI8E36DqiJHzN/BiH6qwCPqwjCd8vw/B
+lPykLniIhC8Vnj1Ye13OEWkbGzgWQlsSljqbwiv8EReaNatAvncGwdrmldB
gVhRMeBHJ5DFskaE1a9CHGPZHfzz6O+JT8gnNayB8Qj3Pfzz50dQ23GuOd2g
4fuk2Ws0H5qFfujSpHcDsfSSHNFgPkj+WamSieuGB9HoE3868iERPSdoGd4N
vbufqX7V5gOjlCg5g+HQ/rrGS8ESNwtqAjsfRluS2A3usyluRisBO7DPbDPz
bN/ahyDqbbDDOF05qGa4va70a4CD2DNS1X0foQAltGeVFQIdu4rGghmJcGEm
4v0TqbegM9uxTpuaCDoSzm1Jqh/AtHyC91ts4n/enwDWQdt62e8SQdTYN/wX
8RqYjyacCj+DBqMMVKs+yRp0Mzz/RIvMp7ct9Mq5C9EwF27St3U9AW5Re5ow
yqGBg+GF4Vv1BKjMTb+497oT/OgxN9CiCVCkodeOSW6HVToO84M0CRDlW1pm
HfIAyD8G62Lu28KNecqHToEPIeqMjei5FFt4Im3+dvXwI6jT4h8hYbf7z/uf
wKLdt5edYXaAuY3p7TucB5wmpRYLr+2hKHJLmcwmF2jCKq4u4+9PUvhrVOln
w8JuLe8Gwu+HuA47YTOg+X7Lr22EGV/132sciweqx91Prx52gpCF5Ahet3gw
+pR6q3DJEWj1fKZr9OJg+4D2qefZjrDTXRUjWxwDkZY721p2jiC7q5Rf9+s1
RNKfkjqkbAseEQdO259D6tw6pyNa3rZA3hog0p4bBMd4SZr7d2zh1c4hAQHh
N0Dv88bD1MUOxM+SfZOgTdm/X/nQ7sFAsmv6Pm9ufJiqS87b5z6iDs5Smpdw
IMyGrZwrGArEszMpzwcDy6Be0PsTwaAcJT9tlvwWLljyO+nwB//n/dFgsfbD
iEY4GEaaN04r1FSDM0VUlfeNEMBcO1o5erka1K2kH0dcCQGXpvh08Yoq2OJx
Xq4WCwGNBdbhtPEKSBzLs5g/EgK5iedO/JiKB/trx0/203rBEebVRkq7BHha
Pl4gH+UFZKoXrVBfiuEY7k++ekooyLQ/jc/W9oACjT13DZJwCHCb3I4+4w1u
s6EyszthYJ54+gqfoj9QcJXlcuLC/vP+IAg3GBMwmAmD+xSVbyijV4B6pzF9
4VcXfAo4Y8/9bQhqNVrV/RUKQY6+N5qjeAgyp8hKy5UK4RbP8TXmO0PgyeQW
dkK9EAIkLC4dZRoCXbUVm4BrhUBeeD2aSuorLOaYGLR6VUKsGjsIWOLg44uM
eLzPwOB7fMzqaUJeBw+2Ka2Gk8AVJcdqLf4QfAZ6sRZD5MMErlHfNf9dPgic
rEMJeL/B0UXzYy0Vg/D+SJOPQHwJUN3fflwmh9m/f4eR1fBNm0Dvs9nIWqyf
OIFHf5+yazvQtc8dsXxtjULlkKz8iJXMJhpkP5F9RMuUw+jFO11C7tFwtLxW
c9KoHK6h1GScY6JhscPzz9bjcqg/L5JYXRkNrgMQcOhpPshTj1yYSQsFl0Kv
JR6OPDjNrlbxXjUElrkFKu7R5UGgwaZNhWwIDP+yN9rcyIUL7VcwpJwh0Iyt
+OE5mAs4VIqEOmkIbI2sbEfvZIEwxU2NrUtB4IwTkthizYZ3gstLMZZB4MAe
+ynLIBvsHtCcWfgYBHaq9Oy3w7NBNfuUk/T3IEgJzEbT5SeA6+Lddu1RL3Aa
WoRsjgRgzC69yK/mBeuURJ+bfJdh0CGOIqGyA/T9A2qt/uVn2MtoK5r6l4+h
85Wvl4zcPKwNrEbh/Qvq6cfFhcQWICD7wD28b2GztrKFJmsN2Hb+9OH9CirV
pHtuyLz6Tz8VHtxY4SdI4FNnNKQT0tH73IzBfC4TJbCOrv9HvI+hliLAD+9j
GIrIEUx9hAH/LMl7zwNS4SudRqMbMt5TUSW826GpMEHiQ3yiAQ3xzlXRk8mp
MLdReb7lKhqMterp0SWpcNEojMfrX/6GAU8RN7yvwZKa6gne1+AAe0JBUu2w
ZMlpIu2eCTvyC3JMZ9qB+eQP5nyXTKCTqzMV4m8Hu/LE4GTrTOCVjn4ux9YO
jWrXKCPNM+GqzO1TeeglyCmL+3ZXqBNUvp3Lq+hagnRfuzRrwU6QEflWhO1Y
gngOWnkn/k4QcfGtmG9bgrdF2egnvJ1wH0N9msFiBXJORjYo1HbBcq+pQuit
FXA7/+pDZU0XbI5S3Q4yXwG66w6XC6u7YGfms2+g2QqkuBj+zPzSBWMXtLDU
Zi2w8KnQoD42DfJfNL29fr8FyI34H1OlpYEOXRHrL48WkCmlqAirTIPN6KS4
yLAWsGGZ22XvToOCP9VlPha98IDNcYUushAeLsWmhmouwYku2utHUjpAvGXK
vp57HVjiytM3zLCgcm2BORhhuUCWqCqEbc1/WAkg/Mys11Ub4bf2xIVVXOtQ
JxViMG2KhQEhprZb/P3AmE7k7GJeDHq/i/I4+fqBpPv6pphZMVyN/eEReaIf
Brlpf8SbFIOKgpTGUa5+yLWv+3HEuBh6xx1X8H6Jn5f1fuL9EjFMIn1yJ+bh
Gq967o51Kxhw37Dq5J6H83KbGS8RbpOo8l3nmgf6e2qJaQ9a4YIqVyozwrPv
YqKa7reCyMH3vnjfhFVlVx/eN3G5QlxMyX4KKD5bRM1nNAC5wccBi4dToBD+
1vwAwm6vnokk2k5Bof72inxaA6zU3PAYt54CIRYTD+fUBkh8jzOzIZqE4IXd
HCKpesjckWa1P5YAdzNT2rrHUqGzOV2b2zMJur3jc9xnUuGoeVCdl14GJIuF
u2Bxqf95fx64jviB0G4qnI9TTMH7KpQS+5/gfRXL43tuLL86IJWnkHhNKRtm
66dVVug6QZj9rNyycTYM+xeMRV3ohNN0otmzAdnQqeH95JJ1J0iR8nNPFGfD
eW8fks1VHKjl93vh/RZlkVLvzrqvQ7meE/vJKix06/OTzT9dh1crlM7jlVgg
xpB/N0T4FM9Z3isIn1afKm96sg5ofWNscQUWGl9WWeF9F5aVr5zxvovn/h6B
Tu9mQXfj+uOK4WYopWbSWH07C7gMBisphPloVbtmwmdBjM0fhIeaIeToY+1v
YUj/G/ibmWOwGfIPJD1uvjIH9UHR6FdvW+AIq+HRXY05oI687lQZ3gKdqOpL
99Tn4F3WQDkHcty+vCPo1qM2Bydq9Uk9Q1pgaU6+E+/H4Es0i8L7MYoVbjNm
haxCK5o/Fe/DSLkusvYmvRZcXFDS1IMfweRM/OtR20noy0qLupxSD6I3Be2/
IJxm8JRFBOGWYmlpL4S1MriZFpPrweqI2p4iwlu7jQxZCAvwmebjfRgDQl5H
ODWLgKjE3Ij8wzLcln4X14T0FzpWL3YZr7UCuvnRmF9wOjzSlXqvaDwKXKyH
R2+lV0F/T2yjl+Eo+HblakWmVsH7XjseIp1RuCb7U/LAxyow6kN5elwZhePJ
wOb4oQo6Au1WK3LXoCMrvfLQHAaKrKSqUQiH5OcrBM0i10WOelPVnDVgqahV
C5nBwPW+a3ua2WuQ0NB9LXwaA5OrC6/xfgyju9lxeD+GwpcI/jflq7AVUGf7
agUNuoXK7ygRVoipFydB2C0hwfl36Sr01jVWLi0h19XXf3SXS1bh3nyzyuAi
GtT8M03kYBVycJ39XPZowKYTE7teWIXnbW85jR4i1/GVDJW186sg4s6T02KL
hj8S+kH35VYBeyrvgqwNGnYijX47HK+BvJfPn0deSoAhkvysuzQ10Cfw6YG+
dAK867lFr4urBm6zNTk0WwJopxx1gb5qsH4ve+jyXjzMvGIcfK+/BmKNZ0ev
vMaAo30t7ROET7x08RJD+JUsqzt+EZ7yGHEz7hUGkkkd5kkRXo4NpClCWOmp
wa/Pu83AcCda4ZR8Ghz183J6udkMVYffOS2LI9fPD713ZSeagd745ZAtdxpc
yxcymUM3g9UnT9QKTRo4eFcquBL3QVNMg6WbTBF0xH00fU/UB3TH2j5ESBeB
95B19MBuL3wk7aPuliqCc8fO9TP/6QUJ3HfXQ5JF4F59NTDt3QAICfsRp/SW
gMb8/di18AEwZqU/NYstgXn5xEavkAGo+iXgZ40ugYCQ4eUjQQPAPSz/fbW9
BO53fyJHUQ1AC/MDj/BLJSAWf0/yU9M0bNYLepmO1IONBrbPvmYaHtTdR7/7
Wg9zSSks5IXToF+9bEv+vR4s/zy+8f7TNChUOhx2Gq8HymrWtXN3F0BHRIeq
Jxo534zfldy2WgDaXoxryIcWMA57rVJ1awH8BbsMDsW2wGzriz4m0wUgftou
5RfXAhGRqsb0GgsQsyxK8mC6GZ6f82CKV1iAJ73KewuzzWDsOPemXHQBpBJd
tpeWmkE8R5eyn2kB1h6mra/gmsHDWkHndO80GH2zQfFw1MP4jxQLHew0jC37
BlkhnH7Q4UVk1zSsHi4VXkLYgUs++VvbNOyJLjQ6ctbDZzWylUIlHIw5J5N6
nsHA2yCRbcmLOLjJvSo+j7Ahs9thKRQODF19wqvOYoDtYzO3lDwOtLtYN0OR
fr/f5Zy19uAKRKPLSh8ldwA/yx+p9vYVEDtyhONeagfUG/u3f/m8AgahldNX
MjrAPPborc9RK/Ds6P2cszkd4Gv6LMKveRqyQsNf9Uk3ALVuyG/D2mm4naSd
UyrTAEflC+t+I/sf09S58UyuAdj4B1/Gpk3DhYXLspfkG0DHGpfxDKmPrR5G
dP6SL4SFUm02baQ+Ho72JRNDFQKb3tFfP3aGoAF3+2EEUo+rr2IHooeGIEft
4sieaiFsShhqtw/UAcvaDSGgToKxmJ/c0F0HerXf3BJIk6DUTqaGt6EOOm5P
sVv9SIRgBTcz6pI6UKZYqBaeS0T6Adn8ZdICYFHz538vGQ7CN5h8nP98hqNy
X0dpxMIhnbxb+PnqZ6Dg9NPc4UZez3mFDp74DL9JzlTPMYRDaJHiCcaHTbAl
zHf8Q3oqvFg/gj3o2gQUrKG7yfmpcNcmV2LLrwkEfj6oK61NBY2Zq+/GI5pA
pU/Zv6M7FYJTaSN8eIZAz10jin8pG/R2mkzfVQ7C5IZBn/OfbBhNll4rXhsA
4knbqwwcOXBHM817kG8AOLA+jfkXcuA2a9z4tRsbsOZj2nhLqwdyFFQXhIw3
wCeh1tENYeXgqJ1Bww3IHPDvTEN4eHTxcKD+BmBpNYUGEV4ecYvbXF2H+sUZ
5icC3aD3rUfMcnkdjud7YHMQ5p3011KaW4eO85mqLILdsD573p5nch08mgYr
vRH2Fmc/maFHeL68+g2xxHFo2vddX3J+dtnHCw2ejcOdHnne8G6b3Nv76gbo
iEmg4qSxwGwQmsFyeQNy2tqDthFmo8/YVVHYgDtkuyLZMljg6Ki79kh6Azjl
xVpuyWJh6KTkLxXk+xdp2r2Me4vM5ynltDQ6yPZePxM7hvCT982yURobwPuw
sM7yHRbOB/Za8l/agDofmbk5hGkD9W5aGWwAmsgi4/grLKjRm3bm625Ab/Bp
w6mXWKg84ZCtrrkBo7lteY8QPivuEzSptgEzaKuDlAgP65XmZV7fAJqodrIp
HBYi7NiDHZD9++I5YwL9OvL3E6c0dpHv81A/fdAGYf++LKpA5O/dELpA34Lw
7rjzoh4KeT1cx+k5thuChH+nDMhsQBPLWKwRwjml5MdHRTfAmXdnnRxhtCpd
yATnBvCIManld3fDt1a0ZDQy3mpWLIdkJXvgzidJZgdkvP0G6EuqEH5ksTfC
abQBJJXkdM1SPeDD0ZLQgewv94SfdzDneuDQpxJKD+T/awnovw9y6IbaJVMJ
T/kN6Mv1oL2BsNjj7Jcvzm3ARF3Xnx/23fCB7M+4/5kNWO3jWgxBuOvabHIu
8nrkaZ/zB32R40Eg6aou8rqOaYwbpV83yPawxUfzIttT0fGDJKAbvDzD1saZ
NkCOddfhT2A30FO+W6woJ/hrdy6MeYplEBid6ZLGEEjgJFaa7u27BP4+ykHM
gWw/7T0SP9jBwkbENc6XyHicV2MRS0GYjbrY2BIZ3wgqlWeOf7Cg+PT4e6Q4
A1yLUxdqFwuXhhjTRpHjVY3N85eqJRaCf2p/OYWMVx2mE2WLMK9k20LnlQ0o
vmDdUYJwif1FZjtV5PzJoLpOaoWM96WVdtatNrCuGf9Nc9ANxu70VZ+k7YDr
HePbvfc9Qd2IZ3adqg10Wmy/Bypaw4uRMx3m+O/rqLLYdafnf/ZT/K/rQ/+9
3hISVPknjzUBAouVmo4pewGVmOoRKa4EcBBz9E+97AVeXd28ZXwJcD01QU1G
y+v/e39vom9sMk1VB7K+J0y0ZBLh/CRLlee5z8DKQJp08k7Y/3w9/u/fr/2v
4/0xw/d277/8xqnytNZ4n/GND75ReJ+xxk07fuZ7wxDS/pn8EV056J/8OSkk
iYMjR80y8L7iOd5dOcFf60Dz7oTtpno3SAbn873cWoeHWJ2dLYQ7rdi2tXDr
IHG1iPkPwnfkX7QzLuF9wMySJBrdQK5kH09MOg9hiw8N8P7i3Yz69kfC8+Aa
VV6H9xVreyesU3l3gcioSeLLinS4WfKhdn5nFs73aOI415pA2HYnyDkbB/cC
Wc3w/uHMajoixfMbsHjZV9k0sxuKzTGKPZIbkM2s6NuKcEeidvNPwQ0QLyil
SMjqhokp7BX24xtQonkmwCW7G8xRfYpT//IR/3p4bV3hX/7h3FN2k3w+DWB8
Nvp83qVEOPDzdebp8mKodnsivekYBr3KVSqebLUQOnomXD3UC97KEKVKk9UC
mTdmKCDcCxi4X1im9dbAXDPd3QORXhBOScnLmloDHbRam94xXrC8/NHftisJ
rPrvqZCS+oD509qez6z5sJFza6bPxh22SbaDfQ7mg0T26xGZp+4gRGdFJbSU
BwXpE417793BiB3r1d6ZBxKpMnkNhe6w4H64YQMzBt7geO9NWTm0W7WmR+aO
AalKpu3XwnIIklZ6NGEzBmT6HI8fZ5aDFnWVwkmRMSC3DH5+LLEc6H6IxV0y
XwCyn4/9WXmbQVnoUnywCVI/Fln4WyGc+8bhI6/+AtwVP+e/izDbj9ikkqsL
oJdL7f+OrxnkZ27H1CCf/5K945M60gwZxE1HDt5cgCt8GEsuhH2lXv6KQT5/
n96Njx7pR2890Jw4rbkAfnsnpkiRfvRa4/xF9dMLsBhbYBFl2wRz2Q9y7DiQ
v5/gSB+MnEcsnja+Jb/mgT9lze6JUxOoXbMzIe6fh610uw4r5Ly6bZnVfAE5
3lWC2Q2OXO+Bio5ddiZkPgo780aIBGGb2ITmJmR+FHZy6Jww6gEOOyVHV2Q+
qi3Sc2xAmPtGd3vvxYX9+733J8QuGJwj8M7aPNd1DgIHEaWQmpAT2I/vZz7e
nyx642M83p/82Z1cf/P0DAStct64bd8EUZLJuIGzM5A9OH9Rw7EJ7iZPisRL
zcBKhufR4y5NIHWM1/KO7Ayc9mCcWHBrgviXKzR0xZNgHGJzg0sN6ReouBUy
xoZgdcbl/ixfOWie4kA//pcvWSkr6OGdf/mReyxnqXvIp4H4SyUL3ofM3zLJ
2sQ/CLr3jP1Qw0UggTI/3MI9DuUpPme5VCuh+rbg/CfGcTjWc356Rq0SBrKK
Eu8SjUPrqYwrLlcqYfWHkong/Bg8C2ApPHCtEl4kBvc9PzEO79Lru12yq+Bh
08Blc6Zx+Obj1h6fUQWqtPfEjxOPg8P56S+45CrgMthm60c+T75+7bNSfBXS
j4zae4ngQCpbb9KjqhPYnQsKlLlwsNaFE1Wu7oQEVzPZTioc3Dl6sgpT1wl8
Tw9V66+vwVf9O5o3mzrBPO3+CDsPDpojict/2XXBmYE1RTJmHHxYulGxbNsF
3ZJEdBhiHCi9Pfqlz7oLHMNovsYsrMHi+daaqvtdoLVabJ5biPdZtloxLWNA
2e97/pmiNZBzOnwkGuHThWf0JIrXYHor1PLjCtKvTHhvSZUg/f4T5rK0VQzc
2dNd0BcZghbN8LmXtmVw2939utPpIQhNtx/LcSwD08ed0jmyQ8BrJorZci+D
686cTPNKQ1B8bKYa/Msg2Rc1/DgUA1426wIsWQkQGjmbq26MgUTNEQrfogSw
/iJ9YwXp82eOZjeZtiTApSl/ijAPNIgMPPOTGUmA6LbNbE/hagi8pPoKK5QI
Iao/afKFcIB7bTmzwIIGnkzz77rI/mQ9mJOVwYqGS44Sz6MOI9dPckVZkeNo
uCdLwf/91xqU7PU0prOjoVzlQZGnMA5enFI/dNcBDZ3ibLOc3DjIVhaKiLJH
w+bBJ+o3aXBA7DSVuWeHhuOTQ9kffq+B7seE2tu26P/n+XuuQPWTVVYN3P7R
yu75KQSI/5A7XveoATrNRZeA3BA4/dqTUpGrBu4QiYl+qA4BU/bfMUK11VCV
Z/89Gx0CqC8XO0nKMyCoR9ecwjAM/BQvV8xbr8DBpcYq2fYuIHGwCjWzXYHH
/BryhR1d8BMde0HDfgX0I/Vl09FdsCzWP3/OaQXED5mfi8N2wbR5xpfQlV6I
iYhYj0kpBM/ZyKPX9kZgW/nwz5DLldAg/T7tEPFXYPcxi7qmUQlag+G9z8m+
ghUpxTe0ZiUMPw4h/k3xFXK8sni0dCphMaVa6I0iDj6fPHzf370bnnySfZ+K
9MPXe7iOpnh0Q2TGKecBZRykcatajnki81f2CW0qVRz8tLUrPv68GwzrX/So
/ss3LY8a7hZWGoacl1eybPlLwFYc9fHLuWFw9T5OaiNQAmgWaoU5jmGgfliW
fU+oBM4S9XyjPzAMH0wMr1uJlMDxO4ae1crDcLGs9PuHgRJwVfr066rMMFzO
wUpg+0rgzJ7KygOuYTCIRDYCv/5ROjPhTzEMt5/L+WR2lQBRxraCdEAemJr4
SRUXhEDuGmfxbHAeiB5caHGrCAHaZ9YJavF50GSVv0rcHgK2h0pfpufkwa3a
x0wBwyFQYVnXcNomATCG23MT770gNVacxdQuAXRe6NcxRnuBOx1DdbljAhgO
0z16GucFOt4bVsyuCWBytk1wPNELEoy6yy1HV+BowfkvPB87IL+/EPUCuwLr
3eNvSRI7oBf79HJN5Qp8phO7WJ3UAT+7Lurupq6A49Unq+4pHXBHR06Z6Qbh
/kqOsqnri+sEBuHU0ZsGBO6iWVGS0SXwxUH1pTvNhPsxJG/KNk/VEtjzV8SB
ykIC71o6M6mnEZi6g3LR6NoCcL56EBr3vRlOJmJ1X6ouwFy73arOWDM4ad7j
X5ZaAPerrpkV481Q/nt3S4tvAY6gPe/yTTZDSY9N5PTdBUDdZ1p0RuqVA1tq
gfR3FuCxl7ybFsKJciTW9y2Q+TMn4B05wijPsiu1pgtgOtqbX6bTCqp6P91b
5HDg46xXZXAGA5wdEW5K53DA3PzjSgHCnRV9jrKncCBvI2NjcxYD7pnHbE7z
4eAW/dNX/OIY2NuUdo7PXQMyzvtfVmYxEN4t8YgX4eKAAGd7hL+8FHUVzVkD
dtPmMacZDMwr8bpJZK/BcwmKq67TGOBlOJz/6frw/v2fDyN1/JJ6BH5pc4bu
6mUCu+3GblsBgdu0jQqedUwCSrUqjvFAA9yp893k7ZwEAVvRC5UIxyWNudh1
TcJmKeUTasoG6Pc5/6sUPQl15FPFBlQNUMNPk/07oA2sUEVBNZ0ZwHNnTwTv
+7buLfilTt0JS2mO04d0ViCMX4oT7/c+XjZId9m2B3QOGpHsJRRAY6IvCc0j
pK6QPzGASSsAFpfvU299eiDoY8Fzl9ICsL4s28L+tgfMKFTE2FsKAKsiHmHE
uASX7MiF5aw7gPzsW2ch5iVQWmYhobPtgHsHjM62sC0BLCrml9p3QMfQ8ZW7
nEsgO//gloVzBxh8MtHA+8DnbDCqeB+46S0bRbz/+1yb/Q6dTDEkqttTO36Y
hTN3zBrbkOMnsOhOPXfMLMxcH6N4ihw/hq90pZ/FzoKu4kTgyEQz8N9SyByJ
m4UaoanD8lMIq9p2lDPMw86wQwTeB87m5us/s7MK5w+42VMg409D1JOB+rMK
h7Ee2uYIrx+0w+EQzrH2ZziGHA8DTNQyyburcO1AaF8zwhsDYYF4P7j5ryBN
vB+c7oZPEl30BJh9P0CH7qqDeaaNa5diJ4DRn9XxELYOYt5gMvOSJsC7X0HG
Z7AOtChyKNgzJmCZ/97en291ILDd9Vpfbw4Ui2LFJWJaQPYYbptGfw7SG96W
LiJsdkWq+5nBHDyefN5gEtcCvj5u6auGc6BK6oDpjG+BncHF40u9OGhJUNpN
PYqFk9+Jzjf34UB452CsLwMWghYF3TwGcLCYG/X1IiMWVn5qlUgMIfOdpRA7
MTMW5Ju12vC+cbaT7nfwvvFb6KLIm9Hr8MT8Q9naCBY+uQ0yDiHc+GJ+0fYr
FiS7OI5gPqzDjTbfa/dGsVDLc/tgc8w64I7wFFl8w8LhrmoFvH8880KkEd4/
7qb7MxUtMQfBMQ/18L5xXa+1h7E868DY2osauoUFu59e880It9Bf0SW3wELK
9Mev7LzrcCpJ7okxwiM9dWgHhMMkhT/mIczycLzxu/sqsDYuTmJy0cDAYqLk
5bEKed2eZY55aKB+oA21z1bh+8Iny+J8NBBVqsiSeq0CLTmG7vdnNHyw9zvs
9KYGOGwa6QXKE2D6mxDlg/c1wC3DFiZRnwBPKmz/GKbXgACmPkyzLwHoIj/j
LlXUwMl7tuH3ZxNAKkGQBu8vZ5hrzsb7y01/4IZFKPtA9yFrnCKqCDgzudkS
EQ5kUWIqRzie6VlkCFUfbD6IoMhUKALu5yNMntR9YFY1v/VBsQjAmUJ/0nIZ
vm/zSOwi9aBbWANn9J1leJB0UKesthOpExU+Uz5YBq/toeCzjZ1w91DVpUe2
yxBxNb0zvaUTNJrnVbTIRoHOv/LLDe0qcGRqOoj3l/MKTBrg/eVDohLuNQZr
UOXNEPkkCAMeMhuU3xCeyNt7G4lwo7YTF6vhGpyk63vTizCN9cY5PYRdbDMD
jrzBgKn7SCTeZ65C/4QR7zPPk+ls4xZYBachl0/uhmgQriwSw/vKpw3R5Hhf
uXwx1e/BqV7wdtp2w/vJKwIOVtw2GwC0VWZ5SmwJjBZ5fWuf7wfygZ3JDWQ+
P/Imu7BuoR9uH/bgWRIsgUdAeStvqR+4XT3jvwqXwPCKOU3cSj+Mjnuzd4oi
8z8VJ4Ur6TSkFdKV4v3kRJpxZHgfuRX9Wwe8jzzc9+YrUZJ5mGS1D3KXbAFy
Edb315H+4VjB1Vq8b/yyIPltw2M4cKhbcMX7xUldab90XFkGZU53atx4B3yx
OP9UA+knPvw8EYj3h18fHCsgReqTX0q43PDnhf/X9RCrcoFzeP83+UcRb7z/
O+ZMkfzvuHW4GLSnivd9z+lYRijUrUPle1LhuovI+Zi+YmG5ug5KlBh+vM9b
8ot3HiPJBngp9nvg/d12cnR9BgwbwNVw+S7e161x2SbZynEdHI+P6eH93Hc2
DOXxPm4/KcZBvI/7hzHDH7LxdbgQrleA928PnU/TxtiuA9Pz5By8b9ts93Eu
2b/82hfe0c3dF94A1Y+W9nifts/0XVkS4g1gFXPgwPuz//v3Qf+/1///D2fa
53Q=
    "], {{{}, {}, {}, {}, {}}, {{}, {}, {}, {}, {}}, {{}, 
      {GrayLevel[0], AbsoluteThickness[1], LineBox[{1540, 1328}], 
       LineBox[{1440, 1693, 1692}], LineBox[CompressedData["
1:eJwV0utL01EAx+Fj0KA5L6yhVDZTal1EgphLWRFTU9ZtXVihImajHFGbiF0Q
u1AiJhUhs8ReVRr0r/lP2PN78eH5nh3YOLCeR/W7taYQwk+diYUwcTCE87zP
fiYSIfw74I7nnOO8wBn3fWx2nrJzdoWfucV1PmSWs1zjD37lNle5yRU2+IEl
Tuq6/YJHNGSfYIx5PbGP6ZBd07r9jG84wS4u8zjf8x0fc8/7pph2HlfaPqmi
3cvbnGY3r6k7eqNu2BneYYUfOctXLEWfa8m+pXt2mTdZYJGXmeJRJe0ujjHB
U4wzx3YeZpKXOMeqMnZVV+yLHIx+Tyn7i+btDf22F7jvnWVmnU+rw36tEbuT
V/mAozyrVntRQ3Ybn+qtPjm/5AILHI7u7Ab/cIu7fK7v9g63+Zd1fbN/cTT6
fg5ozG7xf8nb/wH+EzjF
        "]]}, 
      TagBox[
       TooltipBox[
        {GrayLevel[0], AbsoluteThickness[0.8], Opacity[1], CapForm["Butt"], 
         LineBox[CompressedData["
1:eJwV0stLlFEYwOHTwpyR0daGidZCkCxyWUEgIrQK7KblLSFFM6d06a0W9Q8k
2kqzpHQiW3hJhSgXbsxW5qKl1aaLYNoFIqPnLB5+5z3fwDkz3xQ3p6s794QQ
Msxlh/A+EcKMTiZDmOCuOTcnhNM6q3f0sKa5wbq5VVso4bb5HWtMmQf0pN7S
51qhT7VXT2iPZvStPtYnHGeSCU7ZT2iSbMZ5xIr9vbrorlk65a659oashzlE
I888X45naIP5IDfND+M5Oha/hxZrHfUsmS/rJYq4bn7NK3adc03/aof+03L9
re36U9v0l77U77rNFj/Yocz+N93kM1/5wgv73Trvrked22X9x16Se+ZBCqll
zGem7Y9ojflAfAfmB/E96Khe1QK9wEUy5vN6jnyazePsj7+pdjqvSdN6Rct0
lSN84COfqKM2/kbaQD19DNvr1/taqm84w1nzhlZrlbOmrWdYcvdUfG8k2Mcx
z/O0Mv63PJ+3XiCfIlKpEP4D7qVcZg==
          "]]},
        "50"],
       Annotation[#, 50, "Tooltip"]& ], 
      TagBox[
       TooltipBox[
        {GrayLevel[0], AbsoluteThickness[0.8], Opacity[1], CapForm["Butt"], 
         LineBox[CompressedData["
1:eJwV0jtI1VEcwPEjBZU5RNAUJIGFZNGQphYq0WAQaA8haaqrTdX1QQUJiq/K
ymdlUHpTey41lWgaDeESOLno0MOl6KnY1GPIzxm+fM7v3MP936P/zYmaI8m0
EMKopleF0JgewtiaED7rk46vDuEjD/EDD3Oe5XzlbB+L2ctJzvMlu1nKHk7w
IMfZySJe5xjf8b2uaYe5TvWaMyd5VjlqN7/gAT5nG7fxtM5oxnyK1dqqJvMz
7udTNnILE6rStPkkTyhLDeZH3Md0XuR993ygh3prHuGwhnTBvNK5Yq7geQ7a
T+mepsy32a9bqjEH5wr535xkH2/opl6bn/j791h365/+Kt/+b/5Rp7r02LnJ
+Ht51dyhRS1op/1LzPSsyzynK/EOzo7GO/EXl1Rnv1bZ8X/Ikvgd1q3xjs6l
rFP8yR/aa3+PsqwLWKjv1glu9LxBZ6usB/jN/lflm3fH32Odyzx9sT7GFt11
tpJ32BzfO59tZ442WR9lhcriuxffQ71xdh3Xa602aJfn5ylXsz6fiPfKCGEZ
Jw1kAQ==
          "]]},
        "40"],
       Annotation[#, 40, "Tooltip"]& ], 
      TagBox[
       TooltipBox[
        {GrayLevel[0], AbsoluteThickness[0.8], Opacity[1], CapForm["Butt"], 
         LineBox[CompressedData["
1:eJwV0j1IlHEcwPG/Dnoeonem2WoiFISIREsEFWZ60CKp0FAggZqkYIZGTSmU
W0OBb1fqIEi+XWlLoUaDg4su5itO4qggjWqfGz58f//f83DP3cOVNHXUtWeE
EJK8iIZQzmeSXIiEsJQTwm+W2WPfrsi1ixRyiWIS9nEtIEa/8w/d0Xx9p3n6
Xv+ywTPnVq7y1vmbVuucprhCC9+zQ1izn9GfzPkes8TsdnSXTbbZYtW+R7u1
W3v0ld7SFYbNQ3SZM33+Tc3QTh2wH2TKc/44T+oH53PzGafkUmuXoIZ11vhl
f0+rqeK6832t0Jfaov9oTr8brnHXvlxP6DX38cU86rnHesQdu9tEzPNkM5J+
X6T4yhTTlDHLU/fPpO/VLCbNpRxwmQnGqHRtXEs4pN65gUbzc32ij/Sxtmkr
H81v9JO+1np9oA+1TheIs+idFtlFzYVUpv9TekM3/L5obgj/Aa3RUSc=
          "]]},
        "30"],
       Annotation[#, 30, "Tooltip"]& ], 
      TagBox[
       TooltipBox[
        {GrayLevel[0], AbsoluteThickness[0.8], Opacity[1], CapForm["Butt"], 
         LineBox[CompressedData["
1:eJwV0U1LVFEcwOFjZTaOKJa4tk1QLpxh0CZQEFRcOLVoEUWJab4hYgshxY9Q
uEl68QURZHTIggpKDYpoIW1KtFIIDEpBP0BoYYueu3j4nf+598zMvXO64/bl
gbwQwjiLBSG852UshBecPBHCVGEIl/SCXtQfZDhlHtUF9z3lGd/Nr5x/Yn2E
PAKl9ifcP0maZvO4ruljbdQmbdJVPjPv3BxZilyLU2L/rblY3+kmG3w0f9Gv
rFNv/qa1usJf/jDIHftpbdEMB9aPdMHv/mCd02X+WR8S54HvrtHzWq1vSNLg
XEr3uWt9j2XrWed/62T03CxRGD2X66t8otJ83+dNa535rNbqOb2p7bRx3byn
nbqrt/Q5Z+gxt2qX3tBuPa6/yLc+Ri5619H7JMsMCftJqqgwp/Q1R7lq7tdr
9NHLCMN+6xAJHrpnh22u8DN6luh9uDYW/UdapjHnynUrrtbxohD+A2LcUTA=

          "]], LineBox[{2037, 3085}]},
        "25"],
       Annotation[#, 25, "Tooltip"]& ], {}, {}}}], {}},
  AspectRatio->1,
  AxesLabel->{None, None},
  AxesOrigin->{0., 0.},
  AxesStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  BaseStyle->{FontSize -> 16, FontFamily -> "Times"},
  DisplayFunction->Identity,
  Frame->True,
  FrameLabel->{{
     FormBox[
      StyleBox[
      "\"sin\[Theta]\"", FontFamily -> "Times", FontSize -> 18, StripOnInput -> 
       False], TraditionalForm], None}, {
     FormBox[
      StyleBox[
      "\"\\!\\(\\*SubscriptBox[\\(m\\), \\(S\\)]\\) [GeV]\"", FontFamily -> 
       "Times", FontSize -> 18, StripOnInput -> False], TraditionalForm], 
     None}},
  FrameStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  FrameTicksStyle->{
    Directive[FontFamily -> "Times", FontSize -> 18]},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    AbsoluteThickness[0.5], 
    Opacity[0.5]],
  ImageSize->{440.65802001953125`, Automatic},
  ImageSizeRaw->{{180}, {180}},
  LabelStyle->{
    GrayLevel[0]},
  Method->{"DefaultBoundaryStyle" -> Directive[
      GrayLevel[0]], 
    "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "GridLinesInFront" -> 
    True},
  PlotRange->{{0, 16}, {0, 0.32}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic},
  TicksStyle->Directive[
    GrayLevel[0], FontSize -> 12]]], "Output",
 CellLabel->
  "Out[583]=",ExpressionUUID->"cc8fdab0-d168-4a7f-a460-dfdb68148667"]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"RegionPlot", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"strength", "[", 
     RowBox[{"mS", ",", "sin\[Theta]"}], "]"}], "\[LessEqual]", "1."}], ",", 
   RowBox[{"{", 
    RowBox[{"mS", ",", "0", ",", "16"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"sin\[Theta]", ",", "0", ",", "0.32"}], "}"}], ",", 
   RowBox[{"BaseStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"FontSize", "\[Rule]", "16"}], ",", 
      RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}]}], "}"}]}], ",", 
   RowBox[{"Frame", "\[Rule]", "True"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"Style", "[", 
       RowBox[{"\"\<\!\(\*SubscriptBox[\(m\), \(S\)]\) [GeV]\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}], ",", 
      RowBox[{"Style", "[", 
       RowBox[{"\"\<sin\[Theta]\>\"", ",", 
        RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
        RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}]}], "}"}]}], ",", 
   RowBox[{"LabelStyle", "\[Rule]", 
    RowBox[{"{", "Black", "}"}]}], ",", 
   RowBox[{"FrameTicksStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"Directive", "[", 
      RowBox[{
       RowBox[{"FontFamily", "\[Rule]", "\"\<Times\>\""}], ",", 
       RowBox[{"FontSize", "\[Rule]", "18"}]}], "]"}], "}"}]}], ",", 
   RowBox[{"BoundaryStyle", "\[Rule]", 
    RowBox[{
     RowBox[{"ColorData", "[", 
      RowBox[{"97", ",", "\"\<ColorList\>\""}], "]"}], "[", 
     RowBox[{"[", "1", "]"}], "]"}]}]}], "]"}]], "Input",
 CellLabel->
  "In[584]:=",ExpressionUUID->"3795d3d9-6473-4fd4-8a33-887f33c10a0a"],

Cell[BoxData[
 GraphicsBox[{GraphicsComplexBox[CompressedData["
1:eJxdmXtYzOnfx6MDoaSjptK5JtF6VuUU962sY9rZtApr9GDTbiGxpVgbMg5L
VopynmSL0igltdP5aKZSUwpL4Vexic1IUtTT77qe9z3XNf7gel3XzP19v9/f
z/25P/ew3Ljd+8fRKioqE0f++u+/Y9PStIUL0xfIPsulhwy8F5RLbb47Iu8i
YK8yo+2H5QOM1/9Rt7bDRY2C9cikIYF8HOOOvuu8uBYdxtRrhVWHiz7j89e6
37vEGTHu+3yiQiDnMP7O56uzLbwpjJu/D3SOa7FgHD1qzkWxwIrx1zc1NTpc
bBg/83u8TavDlnGMWlqLS5w9Y7dbeyjfw4Fx1zrP6wL5VMYJY8x0RcJpjBff
fhPZwnNi/I5f+J/h4a+ocl4px43+mEhOEeQFRl5g5AVGXmDkBUZeYOQFRl5g
5AVGXmDkBUZeYOQFRl5g5AVGXmDkBUZeyvks5QcNaxMhyweMfMDIB4x8wMgH
jHzAyAeMfMDIB4x8wMgHjHzAyAeMfMDIB4x8wMhHOQ+3PonrmdbrLA8w8gAj
DzDyACMPMPIAIw8w8gAjDzDyACMPMPIAIw8w8gAjDzDyUPY/eYlBrjbJZP7B
8A+GfzD8g+EfDP9g+AfDPxj+wfAPhn8w/IPhHwz/YPhX9qurV+Eti8lhfsHw
C4ZfMPyC4RcMv2D4BcMvGH7B8AuGXzD8guEXDL9g+J3QtuttfGse8wuGXzD8
guEXDL9g+AXDLxh+wfALhl8w/ILhFwy/YPhV9qeRZnNsjVMh8weGPzD8geEP
DH9g+APDHxj+wPAHhj8w/IHhDwx/yn6i5nXFapMS5gcMP2D4AcMPGH7A8AOG
HzD8gOEHDD9g+AHDj7L+Xy1L/XsCyph+MPSDoR8M/WDoB0M/GPrB0A+GfjD0
g6FfWW/kmHNOspgKphcMvWDoBUMvGHrB0AuGXjD0gqEXDL3K+sLf7Ph8+04V
0weGPjD0gaEPDH1g6ANDHxj6wNAHhr5djcsk8a33mD4w9IGhDwx9YOgDQx8Y
+sDQB4Y+ZT078iwTwjVqmB4w9IChBww9YOgBQw8YesDQo/z8bZc//bjGqY49
H4zng/F8MJ4PxvPBeD4Yz1d+3p3SV4lmjvXseWA8D4zngfE8MJ4HxvOU18/O
aYrSJg1sfTDWB2N9MNYHY33l9bJSi7cMe8vYemCsB8Z6YKwHxnq3zqd79QQ0
svXAWA+M9cBYT/n7opgEl+eRTez7YHwfjO+DRU8WGByTPyYB/2zP1CH7SYL+
A2mHSy+5YBNedaY1kYyeM/1oXIsKdZ9qcl4Wc43cfbl9mmvcGOpaGilb45RB
tp3Jut/C06YG099E9wTcJrFJkXligS7VrPAUh2uIydtPJr9rdRjQA0uH0swc
i8nS7wp+4HsY031OzRHD3qUkOZXvJBKa0r36GUueR5aToeH//jGnuwcOGZQl
VTIOe7b+P8mSanIzrT6QK7KkoZUumQK5hKzzFTbw+NY0JF3rt0BOLdFUDZ0b
oWVLt8Z2eC53v09yM9yvCsV21EXM1zLMqSc/rtWbIAni0pmi5b0a9xoY/0+S
698fn8iInkb7LjlnKp0Rb1X6qqeR+ce8jPsD/CMPzJOYp5EH8sH8Bf84L3De
wT/yQH/GeaKch3J/hl/4x/7AfoR/5IH6xv5RzgP1eqv9yzuxYDSd8yRJbbdG
OhnndrlMq0Odcnwc6s0cb5GA0zSe76FJDeeOujTsnUWKu54HiIQTqL7Fo6Dn
kdkkdNCHKwmaTIfsPd98fFJKzFNsT8k5HPqZM7vpVU8ZqfXu+8SRmNBBLZu/
HqlVkMihyo0eEWb0k4pOkmRyJeHeOCsN4prT/veDR/6aVkWmPJxvIRLa0O7f
1xW719wnvdP/9PGIsKVhRVcbjUPqSc1B7aMtPDvae9d4ie5APUl6FFYQxB3p
/1kn88dFN5DIr9reDQ/bU3maupOqtozV68RPFyyy7+Sy+RiM/cF+HxkKvXhi
pO7akyx6dckOoi1caKqiokN1NDLnlCXdYeuHXtsjHDwrY+8P+uEP7xP+lM9T
6IN+vB+sz/rP/+9X7F/8XgLWNntxs9Olnazdf//V2daj7PP4vWB+wcyBuJYe
kuRX5dgYE8/qHfsfvyeA3weFXXSN+0jO9NhtWut0iX0e9+vrnDz6kDdEvlnk
UdsTcJXli3lc7Kl7iCuaRMfJ/B8I5PmsP6BfYF4H/9T/4g2Pr0fH3p19NZBT
wD6P+dco+fbqCC0Dqn5JZ8dy9yK2v9BvMB+Dy7+NLhKKDeko913mhjmKfoT5
E/U67OyvMeZeCes36D+YR8BdU505kiAL+uFR4+qLPvfY5zEflEZd8OJ7WNL3
NUXzYwIV/Qr9C/MDOKFZ7aCcY0XlRWk2v+2Vss/jfA+ZtjVXILeiPVlnx4f8
UcP2O/of+gF48YEHrzkSa/rvtYNy/+Ra1t/Q73AegXmHltjFtXBpz6V3+r2W
jezzOG+4f4vWckUONCRhy5HutEZSkpn9lCNxZP0S5xF4aMbkk2LBVNpz6ulg
u0sTq0fsH9QT6g/vG/WB94X3i7z7Ol9uT6fVLC/kC7/IB36wH6EH+rEfsL+x
H6AP9at8vqGfo56x/7EfsP/A2G/YH/CLeka/AWO/oL6RB+oV+wOM/YD6RV6o
R/QbMPoN6hN5ot5Q32DUM+oPeaP+UL9g1Cv6Nd4H6gn9DIx+hvrC+0J9oT7B
qMdmQVQZj+/I3mfy/vRc6YpO1t/QD2Mt1U+YNcSRhzpBC98v72O/j6K/qVzt
LQp+kEzEH6+0CO4ofh/BeVAfbLtaVTufBJQ1XJDv02N5o1/Vr0x3Pre4iBTa
1dUKFhuyfoL+s2KqcIu6eylpDea6R4gU9xn0k5M/ZG/VMJWQ9y+XG0kmKe4T
6A9/7z19RSiuIVX6i+dH7LRm87b7MV27uA9clgf7/T3j6qfmq43EcF/ISbGP
Yj6M65reGcTRoX2GzTfDuu8S7k/SoaKAh6QgUdDq8uwAkf86NFd7pZxMmW78
uxfvAjlvuTwsOGqIiAThegsTUkhSW1d7i+domp/3Z18sL514/ho22yRgLO06
eYB//lkWkXS+LJHv06a9v76e/3ViLpkYc5mKVujSfrsYY21vMTnRO2N1hK8B
/WIl68s7VEwsH4eki18Z09m/mHqPn1ROlj+snS56b0ozB845L91WSQKXtArk
3eZ02WkHr52F1WRJkk6Hy5A55c5bpSeRVhO1iWfCtKItqbH+xuDbDlJSl3XG
QDLWmo5JvVHy/HAtKb4ys473rS3VtTm/R36sniz4SSDlNdrRRb8kbziQ30Di
vhQOta/m0pD9L4lvl4yI1z/dovXYgV5Yutl0B6eJ1Wdbd13hYU/FvOy9b33G
lIZo0px6JSb4bg87H9FPjtz/NHeDx2UiWHbJ1KFtkGB+rvDLftlyVzFPoj84
N7U9vq15g82XnBPO50MXicgOm450j/2qbN4M9VuucXi5Bps3MX9ypH3SkPlZ
bP5M7Ns5a9TjbLLruGiDRE+T9S9zv2c2rk/Gs/6F+dS1SVTdtjOXxK6fky/m
KOZ39DvXHNeUsHQx0f/wbY74hT47T9GvTKStpQJ5MXGz+dq047oRm+9ffX/L
WaSlmHfRr2LsNa07Tyrm3zLpKal5WzkJXKV/UGuHYh4uNOWPUWkyYfMw5uPQ
lbPke5wU8/APdx0ramqqyc5Z9LLYX3E/QL/L6bYe4xMsJfx9S28IS63Y+Yj+
NkpLde6/E+rY+Wl/MbTRbf19YnbK9TXHRjFfh2rZJ4qbbNh8jX4nPN7u9K2v
Yt6+7XNh8JfhepK5fdxG7iHF/P2u4c3UOBfF/I15PCJ8WX9JqmL+bog+5aE7
ICPTo0vWcq85sPM0trijLn7LK6JWlV9+oy2WvDtYaOWwpZ8krVo649+2JLLa
JOhYnLUKNd412PBVfyrZnXErJmKNOt3c1xUcGn2LWHi6vRCHjqchDhxeuVsO
yS++x4ur0KFhMrHjlQ95xGth4Iq4Bj0a6ZqoWpJRSGpTPvnznxrSghstQXrZ
JaSf5gi0Nk2mgxMH+CdflJIZP2eJxDqT6Zf3h59mJpWSWWvM37hsM6EfcteV
XScVZNn87Uflu6fQt5yP1poXq0hCmOymsMSCPp9dE3Vz4B7ZmrKGRCy2orx/
P/rm+dWQaYWNrzlSa1rgv95xXG4dufC6+5y40ZomlhV3XyuvI+Ei1acuJ0bm
m6Zlp8I77xMV9ac1vEl29EbN5a35Vg3k59tWHZx4e+o3WHE/eOQ8SZ3uNTNi
sgNV67B12nygkUitJ3UJk7XZeau+0q8zaLQOzVZ/l/V71V2Sui5QFBpYTZZu
4c2waQsl+P/b5wa+kXM2lSzQObRpWtTKSjYfxMRW7qld0czmbdwPruy+/OJ5
0S7yxyJLSdBnxf21cGPeooh55vRe4u3Hkc7VxCNs9QTJN4r7Z0Ksfgfnoj2t
c1W7e8tPRgZO3dwt/2BPyz+sT/n5kowUfTH6x+WOGY2J8H2dp17F7heuXLHj
b3uryD3fVcWCV2bsvmG602mf1izFfQP3Dyt/tQ1+TZVkh9vu0SrnprB5APeR
pT07+6/ZVJO6JZM2cJ3tafp2t1lRGjJ2P6lcHLDvf7+TkRsfgw63zFfcV+Ia
o+LFHXbsfMf9pSil3kTzSwPJjqkXtMQr7i8LvN95X/V8QR4Mtyc/4h0nbvwd
Vo+i3hLdgX9SPp89S6o8y9tFaz4Qwxk2s3jRV0iVd3SNZOcg4b5ps5zsdo2s
4zg1C45Moq5rS698sv6LZLxUrRRq69H8jNMxW6JG7gvnauoFcfp0duGl3Zyn
RcTJK/7vII4hfdLpk9lpWkJ8X/JGmVwxohzP8OSNK0bO503TeXx/C6qaP1Ey
b9M94nM0fHyHlSUtyD0yyzxWQgos5zXxOkbmpdgp1rYl0pFzcvcNYYoVTV/R
cHxlTw3ZOzsnSv6TNd2oFt33i3kdsXqyLFacy6WP9/xm9L1rI5k94VElL2Jk
fm4rUN9Y10io6PlYydyptDJz3MyogCYyu3+MzuFKFepiMVH9597r5Fn86U2S
Ng3KLbvTIcjKJP7Dam+E/RPojEzzCrOgO2R/1AR7ia0x9c1MPNsaUUaKEsJO
yGeYUp81g4P9zRXENPX89x5vbajdGAsrI1JP/vyy6xt+qB098w+nKn5bA5lQ
VZvHbVSl7qsWvPorPYPkbF6jFdytSdN/8GoRaWeTNONqz4g0Dk0wFB9bm1pO
ZHeCNnOn2dD4c6tai73uk4BE6WL+LVs6RA0f1DbXk/8DZSGFPQ==
    "], {{
      {RGBColor[0.368417, 0.506779, 0.709798], AbsolutePointSize[4], 
       AbsoluteThickness[1.6], Opacity[0.3], EdgeForm[None], 
       GraphicsGroupBox[{PolygonBox[CompressedData["
1:eJxNmAn4l1MWx9+LSlkiu1CWGm2iJElpUSKRFq2UlPZdsqVkmUbZiUEkS1I0
WUoKpSwzZF9GdtkZxmCMZcycj/t5n5/n6XTucu4533Pvue/9/v57DxnfY9xm
RVE8Gv9tHrpHKor9Q76I9m4hVUKqhvwh5fmdQwZHe4vQu4TUCdkqZOuQXYs8
zpq9Qmo4t7s+qoXsEbJlSPWQPdXY1Xaeubr62yakZ8RqEvJltPcJ2TakZsi+
6u1C9lNvH7K367Cr51gt8IfsELJjSH3H6DdOOY8GIQ3NiTw6hhwU0iykkWPs
RxPzqe1cHfE2NTdyOlDNHhyg7R7628s1LcUNxub6AHsL8yS/Q9TYNTY+sQ91
HXm0Mwbx22oHxta/y+lwNXn8w/HDQto4xppW+mOfjtAHftvrm5w6qMnj1ZAN
Ic+FHOlekEfnkIPN4xix4vsox8ipi5o96G488jjaMdYcbz7MfWU+XUO6mQ/Y
j1Njd6w2zJ2gP7D3Mx6+e5oDZ9tLDfbe6k4hPVyH3YmOkVNfc8BXH8fo99c3
uX4tjgEhJ6WMdSBtcYP3ZDX5DTEeuE5Vg2uw+0Iep6ixG+o8eIeLA1wj1OQ6
zHmwn6bGbpAx8TtSW7CfmjK+MSETjE3MseIG7zg1PsarsZuoLXnUTjnX0SGT
HCOnyWqwn64GIzVCvVPrUxwD7xlq8puqJr9zxQSWXhHr4JB/Rrt5ytjPDrks
ZEbI+fhI2f4c+5PEMk0f5DHdnCe6bqJ25znP3EzXgf0SceDj0JR9zwq5SNzg
vVh9ZsiF5sZc77BvFfJNtC/QH3N/EiN+Z+ub+HPUYLxUDcbLzQdcd4T8OeSG
kKvFQfx/ieuakCuNB5ar1Nhdqw3xrzf2HP3NMeZc58n7OjV2Vxj/AmNf6t7f
JD7mvzXG/JB5joHlZjVYblGD/VbtiXmj/vC1wDHi36YG7+1q8N4pDtZxX/kO
UV8LHQPXkyGrQx4JWSQmsNytBssyfbHuPvVdId+Jb3HIyJQxLQnZJ9o7hlQJ
WSomzuQvanzdK27m7tcf8VdpR8wTY33bkO+j3SZl38uLXEPUPvf14SLHYM0D
+gD7QyH3GGdMyu0V5rrMPFY6v1TbBWJ6xHlwParG93rjgXF8yjHWhjxuPPx1
SDnOmpB1jrHmqJRt/xryhD7A8pTxiNM3bLqE/BDtv/3O97P6Iv5Lzj1T5Hdm
vf6eU3OefcJHx5B/R/s/xnwh5Hnnifmi4/h6WX/EeUW9wbHHjX2X9TLPsWdc
u057+n/XLxjfKvJb+FrIG44RZ6OaOO+HvKntm46x5j3tGKsXeewWslXIpJR9
vRPyrn6x+8nYb4f0C5uuIT9H+/OQj0I+DpmS8poPHXvXOB+rwXJMyr43hfSP
dveQX6L92e/WHZ9y+9Mi80H4305F5n5wQLggXAYuACeA98GR4E9wPXgKfAX+
V0dbuB08CD5U8j/61RxjTQ1tmSu5HRwIPlhXX3A73g3ej5IL0q9iTHDA++q5
Ht9wmSbOw3fgPfBK+N/+4oYjwY3ghfCg+vpq5Fwtx+qJtYm+yKGpMUr+R/9H
6+R1cTdzv8DN96m5vuA7Lc2hhXMl/2thHHhQW/Np6Vxjx9qIq5W+Ss7XWuxt
7Ddw7DCxttNvycmOKio8r71YO9oH+5H2mzrfzhw6OddCHyUH7OwceI82BviO
sb9t7P2vRYXfwZ26irWbfWLAf7obs4d98obPHGcO3e2Dtad2xIb/9BZLH/sl
/+tZVPhfr6LCKbuZR1/XgBWu1K+ocDi40XaBv0iZ98H3Boi5m/MDxTTINV30
0Ve8g50red5gcQ+xz+8YaqVuUeGIQ8Q61D7Yh9nfI2UeNkrccKfTigpHpA+v
gVdME8tI54g9wXk44Wjxl/xvbFHhf+PMj7ExYp3oengZnOysosLhzigq/G+y
eKfYLznipKLC/6aY5yT9gu9MfcHrzjYu8eFN54hpmn1wTDfPkrvMNR6cZaYx
LrRPHDjVDLHMtA+Wi7QjBlxoljEusV9yPuzgd/DCPxYVTjldXLNdU/K5OcaG
C12h7WXOEfty+zOdv7yocLgriwpHnG3sq50DHzzqGvHNtY/d9e4FbxXv6PNF
hfMxBz/gDV9pbLjTDca+yX7JvW5z3R32S243r6hwu5vFPs/1YIF73FpU+NwC
8TE233h36pd5eMs9rudtXqi/RfbL3znD9IXtErEu1NdOcrTvigrnWyTupcaA
1y0Wxw78rSNlDlvyuaX6WmafOHChFfqCN8Gf4HFwuAfFxxgcC262XGz04UKr
jb1SX/hd5Rksc35VUeFnrIGPrTEufTjOk66Hr6zV3zr7JZddpi+4U8nv1mu3
Wh/MwcvgcE/rC+4Df4JnwfkeK3KN3uD51y3yG9BMf/ApeNKKsJ0rp8UHvOhF
41B31N/OMV8tZe7GHNzmdePRh8vtGvPVU35L4V4v64s4cKhX9PWa/ZLblXZv
6BeOBfd5z3n6G11T8rLd5V8/FRUOx3r42Nv62aiPkpeh4VJwMXjVB0XmYJsc
39L8XtAWbgWvgnPBsT6J9oGhm6b8txnm4XNwsb1ibM+U/xb3oWPwsUYx1jDl
v4kl3x/env8W2Qfv+/+K/K7ynvImfW0b7B8bp1msPcjf0e1Dt5NLcwe4K/ND
Nkv59yTfE+7pt7YPj/HWzm3hXcG+ZbQP8XduVe8c96pztDvJ+Tlf9oQag0+e
IqfkzH90r9jrn6yHbVLOjby29nw4j7op938WF3cWbMeF7ianrZVyDte4H7/q
h/3g7WZP6qS8j+xhzZT3jnF+25HbfLFQG+DZN+VvCvmdqF5sTtW8O+RdxbW7
pFwH5AXHr6Efcq3ueP2U6498+qnfFi/5l2cKZ/pK7DXFD8YdPANi7pgqZ7e9
Z8Te1HIfqoofjGDZ1ftFfruJDbw7W7tg3EWc7NN21tzWYmZvqKPa8gZ4SB33
c1DKfdo9rGf+3oymvuPfb1j29huN/cn6wvd+3knOur77MhB/Kf8mae25sj8N
HaOWGqT8u4XfLANtE3eI9dZY24ba9LTNWuqXOj5Ln420aew95fdCU3000GcT
16IP8G4e4Pj+3sum2g9LmSeNNhZ8CU7D37UOchzN/YTL1XZ8lGuHugZ/B+q/
hT7BjKZf/n3sEMd72ybWiNDDU+Yr/B7nbq4pMoc5LOX7yxx3HD7DWCvbI2zj
h9rCBt7B94O3j3fuiJTfwgedP9y1o8lDu9bGn+U49twj5tq5FlwdxNbH8eXG
aq/NuNBjtTsy5beKdwp/4Phem46u7agd37oO5v6YsTrb5l786PdnYrQnpPwu
8p3nfnI3j075vXxamy7i5Dc13x/eADTfI96F00NPTvmd6KbGhnX44tuI/2NT
/qaVduXarsadGvqMlN8P6rOn9XZCyu/KJ9p01z+1z334wvrsaa1ih4/yNz73
kzvCOnz9Yi31sib31AY/x4rnA2u1l/VGXfS2NrpYW2uKyrfye/e/j2fR1vHF
nmkfz4W8+4uhk35+KCrfx5/F2d8cuct8EzY3F+48Y3x7+KZwh377XZbye8g3
hvZA7yCYv/Gunez4SbYHaD/INn/rGGwsviW0+fvBAOOyh9QE9cJbtzr0tpvl
t6H8G9YGz4JvR2PPE97AO097qnkd71zJGaZ67uz96dbQFG02ubfUOnU91vby
otKHN/L7bqhnih5me37oOdYSv7OGe45j7Y+xHrBv7jzj/Da6OfRsawOeAyb4
z2RxUrfPhX4+5TeJnE9x30qOxZu/ZezTgym/aV/iL+X39hvqK+W/k30V8nXK
b0NzcxljrEnGZV8m24Ybjzd39ARrcpI2YGOMubWe2wTb1DB3uxyb6NrPQ3+W
8t+AT3JPB/rNHO6eTAv9TMrf7XNDX2L7/NAbUv4eTg/9bMrf5yqR+50pv/NV
o313yu/t+yHXpnx35oW+KeV7sSr09cl7EfaLos2DulW070iZO2wR7QUp86c3
Q29M+U68F/JWym/w+pAbU75T24T9AynzBWr2kZTrtma016bMBc4JeSrld4f7
y7ebb/YYa4x6G+U493qscyv8PvB+lPs0whrj9xQ2S6zJUz1T/I5yHD3S9lmh
n0z5XV4Z+rqUvxv8BsTm1tBLQ1+V8r2+J/SVKX+jRupzsfeXO8j3kP14xz2p
FvkuS5lDrTF/cn+Dc0v5u0e9v5tyzc9yjnHO5wnPiPM/13PnnM/zrBemfN6c
NbiWiO3l0Bem/AZxF29J+T7e7rlypuR9prnXiLHHU36vXgqZmfKbyLm96tlR
Iw9bJ9XD/tGU+d1FKccj1gWuZy17+VC5n6Gftn4uZV9T5jKc/9nWwF0p1x+1
90rIxSl//25Luf5YS73PsOY5k3s9l/vca/Z5rrVOnVOPm6zJT0J/mjJHvD/l
GqU+F6d8T7gjH4V8mDJPXJfyvrAnj6WcM/kuT/n7wrfl/+ALvkY=
          "]], 
         PolygonBox[CompressedData["
1:eJwtlNlL1VEUhc8xm6nMoYmobKIIJAIJRJAggsAE/RN67KFeM7G5rKABqoeC
JtLmuWwmG8F6aIAGKRposrLJSpvM+lZnPXx3nd+5557f3mvvffNnzqmYnRFC
mAGZUB9DuA05kAmfYTMMhjGwEH5BEdRw/jL6AAZBL9gB7TARPvD9BvSk78qF
C3AfBkAPWAldCKAcbeb8UfQm9IcMuASbYCD0hBfwCkbAb86fQNdDtu7RO6HV
sfaGa7qH+8vQDs5fR2uhk3VfdD78gEKoZO8A+h0mKXeeW9A7kAej4AbsgQBZ
jq3BuSjGs/AXpukejp1G/8AU+cLzWnS3zrDu59iPOddse33FXsrzbfANJsB7
fnMI/QmTlSfPr9F3MBK6wtuYPM/13hrYqRg42wc9It+gGL6wd8t7WsvzdbAf
IgyH53qvc+8Gc+ErjJcnIXlVDbNC8ky9c9y1yLG3B+1loWu9xbUody4L7L1y
+uQaqnajY+q1w861yLEviqn3iu3tEqgP9hiqYFxIPfhSucFY98xD56jchsIT
eBz/lzMM4+NNTJ7Ky3zXvtnrPOe+FVpC8qAJPjrWIdAG212rAteu0l6phvN8
psAxNqKrY5oV1aiT9Rn3ylR0aUx7WjeENDvnoXtMM3QP7kJjSHvLlQgeT/fe
rph6Xr2untvomqqWbSF50+TY5dG+mHo8yz3QAad8tgRd7L0Se67ZvxjTbOo/
QLO2Qu8IaeYiz1djmtVSx3YuptlQjJqVZc5NM7MKHrk26qm9Mc2MZkVFqvGd
pX5HXUw9rt7WTD9Fn8VUS9XoH2W8pMk=
          "]]}]}, {}, {}, {}, {}}, 
     {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], 
      LineBox[CompressedData["
1:eJwd0zdTFHEAxuE9wpEPFVsLKj8DX4HSUjs7GmgxYkAMYERATCgYwYCCggkj
oKKiSDAgioGkqBgwC89RPPfbd2d2/jd7c5mL8xbkhoIgaGWei+j1Cx8D9POK
lxSzn6/k08YSvlFFCQf4zlLaWcYkB9nMbaqJiwmCGt3CYVLsQ3qHI8TbR3Ur
x0m1j+ldagnbdbqNZzxnE8vpYAU/OMF2TvKTldyjgF+cYgen+c0q7rOaP9TT
yQPOkODcs1pKI2l2gz7kHIn2ed3JKGMMMcIwj3jPB15TxiDlNPGXNXSxln80
U8EF/lPIY9YxxUV2cSn6wzm/SLvp4QoR9y5rLy0k2Ve1kuuk29e0jxsk2zd1
N5+ZYJxPfOQJX3jLO96wh1uEPLten7KBGLtV97KPWHtj9F34qiOMMsYQg/TT
SyfttNBILVWUUkQ+OSwkmyzmM5cM5jCbWaQTIY1UUkgmiUQSCBNPHLEsCs28
ypn/wzRJCncJ
       "]]}}], {}},
  AspectRatio->1,
  Axes->{False, False},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  AxesStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  BaseStyle->{FontSize -> 16, FontFamily -> "Times"},
  DisplayFunction->Identity,
  Frame->{{True, True}, {True, True}},
  FrameLabel->{{
     FormBox[
      StyleBox[
      "\"sin\[Theta]\"", FontFamily -> "Times", FontSize -> 18, StripOnInput -> 
       False], TraditionalForm], None}, {
     FormBox[
      StyleBox[
      "\"\\!\\(\\*SubscriptBox[\\(m\\), \\(S\\)]\\) [GeV]\"", FontFamily -> 
       "Times", FontSize -> 18, StripOnInput -> False], TraditionalForm], 
     None}},
  FrameStyle->Directive[
    AbsoluteThickness[1], 
    GrayLevel[0], FontSize -> 14],
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  FrameTicksStyle->{
    Directive[FontFamily -> "Times", FontSize -> 18]},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    AbsoluteThickness[0.5], 
    Opacity[0.5]],
  ImagePadding->All,
  ImageSize->{416.228515625, Automatic},
  ImageSizeRaw->{{180}, {180}},
  LabelStyle->{
    GrayLevel[0]},
  Method->{
   "DefaultGraphicsInteraction" -> {
     "Version" -> 1.2, "TrackMousePosition" -> {True, False}, 
      "Effects" -> {
       "Highlight" -> {"ratio" -> 2}, "HighlightPoint" -> {"ratio" -> 2}, 
        "Droplines" -> {
         "freeformCursorMode" -> True, 
          "placement" -> {"x" -> "All", "y" -> "None"}}}}, "ScalingFunctions" -> 
    None, "TransparentPolygonMesh" -> True, "AxesInFront" -> True},
  PlotRange->{{0, 16}, {0, 0.32}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.02], 
     Scaled[0.02]}},
  Ticks->{Automatic, Automatic},
  TicksStyle->Directive[
    GrayLevel[0], FontSize -> 12]]], "Output",
 CellLabel->
  "Out[584]=",ExpressionUUID->"6fa16567-f449-488b-812b-f4f1c0b69542"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Output", "Section",ExpressionUUID->"2607223f-1f90-4c58-a23a-94eeeea51c7a"],

Cell[BoxData[
 RowBox[{
  RowBox[{"DumpSave", "[", 
   RowBox[{
    RowBox[{
     RowBox[{"NotebookDirectory", "[", "]"}], "<>", 
     "\"\</output/highT.m\>\""}], ",", 
    RowBox[{"{", 
     RowBox[{
     "strength", ",", "Tc", ",", "Positive\[ReversePrime]factor", ",", 
      "VhighT", ",", "VhighT2d"}], "}"}]}], "]"}], ";"}]], "Input",
 CellLabel->
  "In[585]:=",ExpressionUUID->"27ea4158-c135-40b7-b3de-1a768f0c3e7e"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{960, 1027},
WindowMargins->{{Automatic, 0}, {Automatic, -30}},
FrontEndVersion->"13.0 for Mac OS X x86 (64-bit) (February 4, 2022)",
StyleDefinitions->"Cambria Article.nb",
ExpressionUUID->"f368553e-9bcf-4380-a00f-61230b80b9d7"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 111, 0, 76, "Title",ExpressionUUID->"6cd35c81-39c0-48fd-830b-7d456a3fdc89"],
Cell[694, 24, 159, 3, 34, "Input",ExpressionUUID->"33b22c75-f5aa-47cf-a757-aba03fa7960e"],
Cell[CellGroupData[{
Cell[878, 31, 96, 0, 71, "Section",ExpressionUUID->"09784b56-5c1a-48b1-ae64-9874125f0427"],
Cell[CellGroupData[{
Cell[999, 35, 112, 0, 41, "Subsection",ExpressionUUID->"eeecd131-c144-47e2-b47a-2c75fcfa74c5"],
Cell[1114, 37, 146, 2, 26, "TextIndent",ExpressionUUID->"e6d25092-ca41-45b6-8ad4-525823e653cd"],
Cell[CellGroupData[{
Cell[1285, 43, 2240, 73, 143, "Input",ExpressionUUID->"69d3d85d-3aa9-4f5b-8164-a8603914bbae"],
Cell[3528, 118, 564, 17, 54, "Output",ExpressionUUID->"f615bf16-5f90-47aa-87b7-d7fadc145000"],
Cell[4095, 137, 359, 11, 54, "Output",ExpressionUUID->"a238d36f-4ce5-4cef-a3e9-3239616a2953"],
Cell[4457, 150, 553, 16, 54, "Output",ExpressionUUID->"1da25cff-277e-4fbd-8bc6-8a6dec03a7e9"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[5059, 172, 235, 5, 42, "Subsection",ExpressionUUID->"9c08ab25-f758-456f-83bb-d617fea16fa1"],
Cell[CellGroupData[{
Cell[5319, 181, 2539, 80, 164, "Input",ExpressionUUID->"64989a1f-b188-48e1-ae1b-76b6648eea81"],
Cell[7861, 263, 435, 14, 54, "Output",ExpressionUUID->"d2ce8b13-ffa6-4811-b9f0-3dea23d3a418"],
Cell[8299, 279, 685, 22, 57, "Output",ExpressionUUID->"395fa0e5-48ae-4841-a4cd-91bbecadee05"],
Cell[8987, 303, 179, 5, 31, "Output",ExpressionUUID->"38497e21-7299-4ba5-8c4b-6665c41f4c95"],
Cell[9169, 310, 295, 9, 52, "Output",ExpressionUUID->"7beb5017-34e2-405b-9181-8b6cb65c0cdc"],
Cell[9467, 321, 957, 30, 54, "Output",ExpressionUUID->"4ab905f3-58b1-46ed-8624-8243c274cfc7"]
}, Open  ]],
Cell[CellGroupData[{
Cell[10461, 356, 651, 20, 34, "Input",ExpressionUUID->"6bfa97e4-7db2-4078-9334-574123f43310"],
Cell[11115, 378, 455, 15, 54, "Output",ExpressionUUID->"dbbca74e-433b-4ee5-8080-dd2298de43fb"]
}, Open  ]],
Cell[11585, 396, 123, 0, 26, "TextIndent",ExpressionUUID->"39cec182-a215-4bf2-81e2-baa4eb199a6e"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11757, 402, 277, 7, 73, "Section",ExpressionUUID->"9661d2b9-84dc-431f-850a-75d73e467719"],
Cell[12037, 411, 9487, 308, 893, "Input",ExpressionUUID->"0a8b36df-21ed-4130-9d65-d0a68edb84ce"],
Cell[21527, 721, 466, 16, 58, "Input",ExpressionUUID->"5499e233-fd83-4c7d-9b85-367412aebf6a"],
Cell[CellGroupData[{
Cell[22018, 741, 479, 13, 38, "Input",ExpressionUUID->"3db30cd6-516d-49ae-b089-bcfff54becc1"],
Cell[22500, 756, 10701, 200, 220, "Output",ExpressionUUID->"0295452a-14a4-4bd4-8843-cc9ab896f769"]
}, Open  ]],
Cell[33216, 959, 228, 6, 31, "Input",ExpressionUUID->"ae3b750b-14d5-4e52-89d2-75e3af983480"],
Cell[33447, 967, 163, 2, 26, "TextIndent",ExpressionUUID->"46a938ad-3f76-4d1d-ab00-60484f047b0f"],
Cell[CellGroupData[{
Cell[33635, 973, 170, 4, 31, "Input",ExpressionUUID->"c6c6b4a4-95e3-4309-bdab-5f694857b383"],
Cell[33808, 979, 173, 4, 31, "Output",ExpressionUUID->"11eea865-aa48-4f79-a17e-dbf42ee47b37"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34018, 988, 219, 6, 31, "Input",ExpressionUUID->"68bd2c68-b2ae-4b89-88b2-c94ef87b250b"],
Cell[34240, 996, 128, 2, 31, "Output",ExpressionUUID->"d908dacc-ed12-4cca-a23a-f186f6b5918a"]
}, Open  ]],
Cell[CellGroupData[{
Cell[34405, 1003, 1134, 34, 73, "Input",ExpressionUUID->"58613e4b-6b4e-4947-bd69-0984f9e02323"],
Cell[35542, 1039, 2614, 71, 127, "Output",ExpressionUUID->"ddd7316e-dd4a-4a7c-a0d3-f5642bffca19"]
}, Open  ]],
Cell[CellGroupData[{
Cell[38193, 1115, 266, 7, 31, "Input",ExpressionUUID->"b4347670-6e0f-46f7-a474-bcb023deec4e"],
Cell[38462, 1124, 229, 3, 31, "Output",ExpressionUUID->"0d2e608e-67f9-4a9a-9d4d-0f37d19d8a32"]
}, Open  ]],
Cell[38706, 1130, 125, 0, 26, "TextIndent",ExpressionUUID->"04790823-db8b-4fcc-bfe9-f13cf98b72ad"],
Cell[CellGroupData[{
Cell[38856, 1134, 2207, 66, 146, "Input",ExpressionUUID->"1bec00ec-271c-4f3f-a8dd-06f7273ac3fd"],
Cell[41066, 1202, 255, 6, 31, "Output",ExpressionUUID->"0a556395-5d99-4195-9e07-f7df47f957ee"],
Cell[41324, 1210, 128, 2, 31, "Output",ExpressionUUID->"064e7518-94eb-4a76-a1a1-2705c993e0a3"],
Cell[41455, 1214, 269, 7, 31, "Output",ExpressionUUID->"3f20aae7-1e7a-4d65-9b6c-80cd327afa35"],
Cell[41727, 1223, 129, 2, 31, "Output",ExpressionUUID->"8280af7f-fe32-4453-b489-314afa9b6111"],
Cell[41859, 1227, 113, 2, 31, "Output",ExpressionUUID->"15e6467c-2604-47db-a6d0-5ef5c8d0956d"],
Cell[41975, 1231, 150, 3, 32, "Output",ExpressionUUID->"2647fa55-9030-4f81-aef4-1fd313543705"]
}, Open  ]],
Cell[42140, 1237, 85, 0, 26, "TextIndent",ExpressionUUID->"2f0eac60-0248-401f-b530-dfc6cedeb210"]
}, Open  ]],
Cell[CellGroupData[{
Cell[42262, 1242, 79, 0, 71, "Section",ExpressionUUID->"db68c9f9-089a-4f7f-bb60-1fc7f5ccaaeb"],
Cell[CellGroupData[{
Cell[42366, 1246, 594, 17, 34, "Input",ExpressionUUID->"1bd87dd9-c963-4e10-b94d-c3399d7e51b8"],
Cell[42963, 1265, 7747, 144, 378, "Output",ExpressionUUID->"34d94513-a59e-404b-81d4-70d7aa067f60"]
}, Open  ]],
Cell[CellGroupData[{
Cell[50747, 1414, 1684, 40, 127, "Input",ExpressionUUID->"8bd1f3f7-d846-4b3b-803c-80a05f3bab69"],
Cell[52434, 1456, 48106, 831, 445, "Output",ExpressionUUID->"cc8fdab0-d168-4a7f-a460-dfdb68148667"]
}, Open  ]],
Cell[CellGroupData[{
Cell[100577, 2292, 1685, 41, 108, "Input",ExpressionUUID->"3795d3d9-6473-4fd4-8a33-887f33c10a0a"],
Cell[102265, 2335, 13574, 252, 420, "Output",ExpressionUUID->"6fa16567-f449-488b-812b-f4f1c0b69542"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[115888, 2593, 80, 0, 71, "Section",ExpressionUUID->"2607223f-1f90-4c58-a23a-94eeeea51c7a"],
Cell[115971, 2595, 424, 12, 31, "Input",ExpressionUUID->"27ea4158-c135-40b7-b3de-1a768f0c3e7e"]
}, Open  ]]
}, Open  ]]
}
]
*)