TrajectoryMonteCarlo_BFieldScan.nb

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

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

(* CreatedBy='Mathematica 12.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[    166676,       3266]
NotebookOptionsPosition[    163415,       3209]
NotebookOutlinePosition[    163789,       3225]
CellTagsIndexPosition[    163746,       3222]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{
Cell[BoxData[
 RowBox[{"Off", "[", 
  RowBox[{"General", "::", "munfl"}], "]"}]], "Input",
 CellChangeTimes->{{3.7782469813079433`*^9, 3.7782469827193623`*^9}},
 CellLabel->"In[1]:=",ExpressionUUID->"acf3f024-06f7-4b67-aa47-b28ae5a429d7"],

Cell[CellGroupData[{

Cell["Load force profiles, and set parameters and conversions", \
"Subsubsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
   3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
   3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, 
   3.776165628366666*^9, {3.801264123339542*^9, 3.8012641333422537`*^9}, {
   3.801275305033741*^9, 
   3.8012753227224874`*^9}},ExpressionUUID->"5aaf2c87-6956-45b5-ba7b-\
95459d133cb4"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
   "fname32", "=", 
    "\"\<C:/Users/Christian/Desktop/CaOCH3Data/uploaded_files/ptsJ32Jp12_\
BFieldScan.csv\>\""}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
   "file", " ", "path", " ", "for", " ", "force", " ", "profiles", " ", 
    "must", " ", "be", " ", "changed", " ", "accordingly"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"isat", "=", 
    RowBox[{
     RowBox[{"Quantity", "[", 
      RowBox[{"2.4", ",", 
       RowBox[{"\"\<Milliwatts\>\"", "/", 
        SuperscriptBox["\"\<Centimeters\>\"", "2"]}]}], "]"}], "*", "1.2"}]}],
    ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
     RowBox[{
     "simple", " ", "rate", " ", "analysis", " ", "suggests", " ", "a", " ", 
      "correction", " ", "factor", " ", "2", 
      RowBox[{"ng", "/", 
       RowBox[{"(", 
        RowBox[{"ne", "+", "ng"}], ")"}]}]}], "=", 
     RowBox[{
      RowBox[{
       RowBox[{"3", "/", "2"}], " ", "due", " ", "to", " ", "multi"}], "-", 
      RowBox[{"level", " ", "structure"}]}]}], ";", " ", 
    RowBox[{
    "we", " ", "here", " ", "find", " ", "1.2", " ", "appropriate"}]}], 
   "*)"}], " "}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"mCaOCH3", "=", 
    RowBox[{"(", 
     RowBox[{
      RowBox[{
       TemplateBox[{"\"calcium\"",RowBox[{"Entity", "[", 
           RowBox[{"\"Element\"", ",", "\"Calcium\""}], "]"}],
         "\"Entity[\\\"Element\\\", \\\"Calcium\\\"]\"","\"element\""},
        "Entity"], "[", 
       TemplateBox[{"\"atomic mass\"",RowBox[{"EntityProperty", "[", 
           RowBox[{"\"Element\"", ",", "\"AtomicMass\""}], "]"}],
         "\"EntityProperty[\\\"Element\\\", \\\"AtomicMass\\\"]\""},
        "EntityProperty"], "]"}], " ", "+", " ", 
      RowBox[{
       TemplateBox[{"\"oxygen\"",RowBox[{"Entity", "[", 
           RowBox[{"\"Element\"", ",", "\"Oxygen\""}], "]"}],
         "\"Entity[\\\"Element\\\", \\\"Oxygen\\\"]\"","\"element\""},
        "Entity"], "[", 
       TemplateBox[{"\"atomic mass\"",RowBox[{"EntityProperty", "[", 
           RowBox[{"\"Element\"", ",", "\"AtomicMass\""}], "]"}],
         "\"EntityProperty[\\\"Element\\\", \\\"AtomicMass\\\"]\""},
        "EntityProperty"], "]"}], " ", "+", " ", 
      RowBox[{
       TemplateBox[{"\"carbon\"",RowBox[{"Entity", "[", 
           RowBox[{"\"Element\"", ",", "\"Carbon\""}], "]"}],
         "\"Entity[\\\"Element\\\", \\\"Carbon\\\"]\"","\"element\""},
        "Entity"], "[", 
       TemplateBox[{"\"atomic mass\"",RowBox[{"EntityProperty", "[", 
           RowBox[{"\"Element\"", ",", "\"AtomicMass\""}], "]"}],
         "\"EntityProperty[\\\"Element\\\", \\\"AtomicMass\\\"]\""},
        "EntityProperty"], "]"}], "+", 
      RowBox[{"3", 
       RowBox[{
        TemplateBox[{"\"hydrogen\"",RowBox[{"Entity", "[", 
            RowBox[{"\"Element\"", ",", "\"Hydrogen\""}], "]"}],
          "\"Entity[\\\"Element\\\", \\\"Hydrogen\\\"]\"","\"element\""},
         "Entity"], "[", 
        TemplateBox[{"\"atomic mass\"",RowBox[{"EntityProperty", "[", 
            RowBox[{"\"Element\"", ",", "\"AtomicMass\""}], "]"}],
          "\"EntityProperty[\\\"Element\\\", \\\"AtomicMass\\\"]\""},
         "EntityProperty"], "]"}]}]}], ")"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"mass", " ", "of", " ", "molecule"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"tdecay", "=", 
   TemplateBox[{"35","\"ns\"","nanoseconds","\"Nanoseconds\""},
    "Quantity"]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vrecoil", "=", 
   RowBox[{"QuantityMagnitude", "@", 
    RowBox[{"UnitConvert", "[", 
     RowBox[{
      FractionBox[
       TemplateBox[{
        InterpretationBox["\[InvisibleSpace]", 1],StyleBox[
         "\"h\"", Italic, StripOnInput -> False],"Planck constants",
         "\"PlanckConstant\""},
        "Quantity"], "\[Lambda]"], "/", "mCaOCH3"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"ftoa", "=", 
    RowBox[{"QuantityMagnitude", "@", 
     RowBox[{"UnitConvert", "@", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          
          TemplateBox[{
           InterpretationBox["\[InvisibleSpace]", 1],"\"\[HBar]\"",
            "reduced Planck constants","\"ReducedPlanckConstant\""},
           "Quantity"], 
          RowBox[{"(", 
           FractionBox[
            RowBox[{"2", "\[Pi]"}], "\[Lambda]"], ")"}], 
          RowBox[{"(", 
           FractionBox["1", "tdecay"], ")"}]}], ")"}], "/", "mCaOCH3"}], 
       ")"}]}]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"conversion", " ", "between", " ", "force", " ", "and", " ", 
    RowBox[{"acceleration", ":", " ", 
     RowBox[{"m", "/", 
      RowBox[{"s", "^", "2"}]}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"rtor", "=", 
    RowBox[{"QuantityMagnitude", "@", 
     RowBox[{"UnitConvert", "@", 
      RowBox[{"(", 
       FractionBox["1", "tdecay"], ")"}]}]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"factor", " ", "for", " ", "scattering", " ", 
    RowBox[{"rate", ":", " ", "MHz"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"vtov", "=", 
    RowBox[{"QuantityMagnitude", "@", 
     RowBox[{"UnitConvert", "@", 
      RowBox[{"(", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{"(", 
          FractionBox["1", "tdecay"], ")"}], ")"}], "/", 
        RowBox[{"(", 
         FractionBox[
          RowBox[{"2", "\[Pi]"}], "\[Lambda]"], ")"}]}], ")"}]}]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"factor", " ", "for", " ", 
    RowBox[{"velocity", ":", " ", 
     RowBox[{"m", "/", "s"}]}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Beta]toB", "=", 
    RowBox[{"QuantityMagnitude", "@", 
     RowBox[{"(", 
      RowBox[{"UnitConvert", "[", 
       RowBox[{
        RowBox[{"(", 
         RowBox[{
          
          TemplateBox[{
           InterpretationBox["\[InvisibleSpace]", 1],"\"\[HBar]\"",
            "reduced Planck constants","\"ReducedPlanckConstant\""},
           "Quantity"], 
          RowBox[{
           RowBox[{"(", 
            FractionBox["1", "tdecay"], ")"}], "/", 
           
           TemplateBox[{
            InterpretationBox["\[InvisibleSpace]", 1],SubscriptBox[
             "\"\[Mu]\"", "\"B\""],"Bohr magnetons","\"BohrMagneton\""},
            "Quantity"]}]}], ")"}], ",", "\"\<Gauss\>\""}], "]"}], ")"}]}]}], 
   ";"}], 
  RowBox[{"(*", 
   RowBox[{"factor", " ", "for", " ", "B", " ", 
    RowBox[{"field", ":", " ", 
     RowBox[{"m", "/", "s"}]}]}], "*)"}], "\[IndentingNewLine]", 
  RowBox[{"(*", 
   RowBox[{"cooling", " ", "beam", " ", "parameters"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"\[Lambda]", "=", 
    RowBox[{"Quantity", "[", 
     RowBox[{"629.483", ",", "\"\<Nanometers\>\""}], "]"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"wavelength", " ", "of", " ", "cooling", " ", "light"}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"dFWHM", "=", 
    RowBox[{"0.001", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"3.14", "+", "3.63"}], ")"}], "/", "2"}]}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
   "FWHM", " ", "diameter", " ", "of", " ", "cooling", " ", "beam", " ", 
    RowBox[{"(", "mm", ")"}]}], "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"waist", "=", 
   RowBox[{"dFWHM", "/", 
    SqrtBox[
     RowBox[{"2", 
      RowBox[{"Log", "[", "2", "]"}]}]]}]}], 
  RowBox[{"(*", 
   RowBox[{"waist", " ", "of", " ", "cooling", " ", "beam"}], "*)"}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Clear", "@", "intensity"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"intensity", "[", 
   RowBox[{"x_", ",", "y_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"Exp", "[", 
    RowBox[{
     RowBox[{"-", "2"}], 
     FractionBox[
      RowBox[{
       SuperscriptBox["x", "2"], "+", 
       SuperscriptBox["y", "2"]}], 
      SuperscriptBox["waist", "2"]]}], "]"}], "/", 
   RowBox[{"QuantityMagnitude", "[", "isat", "]"}]}]}]}], "Input",
 CellChangeTimes->{{3.7761641907194366`*^9, 3.776164336741789*^9}, {
   3.7761653506587815`*^9, 3.7761653509707994`*^9}, {3.776165733958705*^9, 
   3.7761657855236545`*^9}, {3.7761660864288654`*^9, 3.776166096684452*^9}, {
   3.776166176176998*^9, 3.77616619823426*^9}, {3.776174190264378*^9, 
   3.776174356548889*^9}, {3.7761749521009526`*^9, 3.7761750385358963`*^9}, 
   3.7761751157303114`*^9, {3.7761754140093727`*^9, 3.7761755970208397`*^9}, {
   3.776175670767058*^9, 3.7761756742042546`*^9}, {3.776175798734377*^9, 
   3.7761758140042505`*^9}, {3.776175876496825*^9, 3.7761758797370105`*^9}, {
   3.776176225567791*^9, 3.776176247109023*^9}, {3.7761768649443607`*^9, 
   3.7761768726027994`*^9}, {3.7761787953527737`*^9, 3.7761787957998*^9}, 
   3.776179610088374*^9, 3.7761806957084684`*^9, {3.776181246312961*^9, 
   3.7761812481040635`*^9}, 3.77618714732848*^9, {3.7764500386062393`*^9, 
   3.776450043356511*^9}, {3.776453015047804*^9, 3.776453022894278*^9}, {
   3.776454290742094*^9, 3.776454291145746*^9}, {3.776512178631446*^9, 
   3.776512185775854*^9}, 3.7766090396807213`*^9, {3.7766108388396273`*^9, 
   3.7766108499922657`*^9}, 3.776611073420045*^9, 3.7766114070521274`*^9, {
   3.7766836914685593`*^9, 3.7766836917675767`*^9}, 3.7767639944156256`*^9, {
   3.7767641554678373`*^9, 3.776764156185878*^9}, 3.776771637345776*^9, {
   3.776782201463009*^9, 3.7767822197990575`*^9}, {3.776877846739606*^9, 
   3.776877846969619*^9}, {3.777021676241183*^9, 3.777021700023543*^9}, {
   3.7770217330344315`*^9, 3.777021762022089*^9}, {3.777022084594539*^9, 
   3.7770220859286156`*^9}, {3.7770222056024604`*^9, 
   3.7770222220113993`*^9}, {3.7770223350208626`*^9, 3.777022343628355*^9}, {
   3.777040229651377*^9, 3.777040270580718*^9}, {3.7771120748355474`*^9, 
   3.7771120752125688`*^9}, {3.777114031708474*^9, 3.7771140323655114`*^9}, {
   3.777124490246668*^9, 3.777124501578316*^9}, {3.77712482072057*^9, 
   3.7771248231207075`*^9}, 3.7772255723639936`*^9, {3.777225606641954*^9, 
   3.777225607937028*^9}, {3.7772264347433186`*^9, 3.777226485770237*^9}, {
   3.77737137364135*^9, 3.777371374136378*^9}, {3.777371871087802*^9, 
   3.7773718712168093`*^9}, {3.777372961831189*^9, 3.777372972679809*^9}, {
   3.7773735919552298`*^9, 3.7773736257271614`*^9}, {3.7773748914325557`*^9, 
   3.7773748935846786`*^9}, {3.7773765409949055`*^9, 3.7773765411869164`*^9}, 
   3.7773771379900513`*^9, {3.777377524087135*^9, 3.77737752418114*^9}, {
   3.77737832213278*^9, 3.777378323507859*^9}, 3.777378430727992*^9, {
   3.7773891257443*^9, 3.777389126039317*^9}, {3.7774610138511925`*^9, 
   3.777461020172554*^9}, {3.777677252838169*^9, 3.7776772531451864`*^9}, {
   3.778234441374858*^9, 3.778234480074332*^9}, {3.778234713548769*^9, 
   3.778234713812348*^9}, {3.778236292510146*^9, 3.778236324139155*^9}, {
   3.778248551681801*^9, 3.7782485597932653`*^9}, {3.778322789810979*^9, 
   3.778322790842038*^9}, 3.778322934753269*^9, {3.7784319254458*^9, 
   3.778431926013832*^9}, {3.7785050069928246`*^9, 3.778505095327877*^9}, {
   3.778505510204607*^9, 3.778505514157833*^9}, {3.7785056264672565`*^9, 
   3.778505627617323*^9}, {3.77850634806553*^9, 3.7785063761651373`*^9}, {
   3.7785068683292875`*^9, 3.7785068685483*^9}, {3.778508125528632*^9, 
   3.7785081344707017`*^9}, 3.7785089685196924`*^9, {3.778512730933878*^9, 
   3.778512732057865*^9}, 3.7785128451581945`*^9, {3.7785133717459326`*^9, 
   3.778513380488564*^9}, {3.7785135038623714`*^9, 3.7785135039511585`*^9}, {
   3.778513882983816*^9, 3.778513883052623*^9}, {3.7785139881792965`*^9, 
   3.778513989040989*^9}, {3.778580089267271*^9, 3.778580090294496*^9}, {
   3.778580193523275*^9, 3.778580198619627*^9}, {3.778580694010886*^9, 
   3.7785807067567725`*^9}, {3.778580977257877*^9, 3.7785809776647873`*^9}, {
   3.7785816253504715`*^9, 3.778581634876973*^9}, {3.7785818348308673`*^9, 
   3.778581834964508*^9}, {3.778582084632372*^9, 3.778582085116069*^9}, {
   3.7785839806908946`*^9, 3.7785839807437525`*^9}, {3.778584629331064*^9, 
   3.7785846330919995`*^9}, {3.77858620480285*^9, 3.778586204882641*^9}, {
   3.7785865438305635`*^9, 3.778586546577207*^9}, {3.778587063774083*^9, 
   3.778587079307542*^9}, {3.7943917315833445`*^9, 3.7943917462031775`*^9}, {
   3.7943993787477093`*^9, 3.79439940508922*^9}, {3.79439974047785*^9, 
   3.7943999102117624`*^9}, {3.7944000427045436`*^9, 3.794400069928785*^9}, {
   3.794400181990406*^9, 3.7944002039072456`*^9}, {3.7944003153392944`*^9, 
   3.7944003177707515`*^9}, {3.7944014839194455`*^9, 3.794401484215206*^9}, {
   3.7944030500942745`*^9, 3.7944030610929365`*^9}, {3.7944034255038767`*^9, 
   3.794403427639406*^9}, {3.7944036365259123`*^9, 3.7944036665855007`*^9}, {
   3.794413638482148*^9, 3.7944136399484124`*^9}, {3.794413711472114*^9, 
   3.7944137115628767`*^9}, {3.7944240465321627`*^9, 
   3.7944240531391163`*^9}, {3.7944289146814613`*^9, 
   3.7944289152098975`*^9}, {3.7944303389501343`*^9, 3.794430367759261*^9}, {
   3.7944355785451837`*^9, 3.794435604134654*^9}, {3.794440519009368*^9, 
   3.794440519526306*^9}, {3.7944405792578974`*^9, 3.7944405806060176`*^9}, {
   3.794440613511942*^9, 3.7944406136195626`*^9}, {3.7944415129220657`*^9, 
   3.7944415133168836`*^9}, {3.7944415450847225`*^9, 3.794441545232836*^9}, {
   3.794441630642764*^9, 3.7944416814348936`*^9}, {3.7944542963291955`*^9, 
   3.79445429639754*^9}, {3.794454591663227*^9, 3.7944545933616095`*^9}, {
   3.79445465437789*^9, 3.794454660437259*^9}, {3.7944556738825035`*^9, 
   3.7944556743542356`*^9}, {3.7944557208102355`*^9, 
   3.7944557211461773`*^9}, {3.7944557522253675`*^9, 
   3.7944557524654446`*^9}, {3.794455788580146*^9, 3.7944557891071796`*^9}, {
   3.794484261132514*^9, 3.794484261651701*^9}, {3.7944842976101856`*^9, 
   3.7944842977658577`*^9}, {3.7946054360981073`*^9, 
   3.7946054371171503`*^9}, {3.7946140583265257`*^9, 
   3.7946140644404287`*^9}, {3.7946181393658533`*^9, 3.794618141151183*^9}, 
   3.794655578350626*^9, {3.7946777806380777`*^9, 3.794677782018899*^9}, {
   3.7946780498786383`*^9, 3.794678050098687*^9}, {3.7947574879710608`*^9, 
   3.794757490437614*^9}, {3.794759610323228*^9, 3.7947596608794155`*^9}, {
   3.794775820754927*^9, 3.7947758294865265`*^9}, {3.7947759351731215`*^9, 
   3.7947759354891944`*^9}, {3.7947936440156813`*^9, 
   3.7947936500202117`*^9}, {3.7947975354674206`*^9, 3.794797536631342*^9}, {
   3.7947994308070927`*^9, 3.7947994325905094`*^9}, {3.7947994634169645`*^9, 
   3.794799467079652*^9}, {3.7948262468383207`*^9, 3.7948262484526834`*^9}, {
   3.7948264564614005`*^9, 3.794826460757368*^9}, {3.7948265945110607`*^9, 
   3.794826596095416*^9}, {3.7948358694855704`*^9, 3.794835872614273*^9}, {
   3.8012611162574234`*^9, 3.801261122825426*^9}, {3.801266892436345*^9, 
   3.8012669125845547`*^9}, {3.8012669477628183`*^9, 
   3.8012669493186655`*^9}, {3.8012671428375564`*^9, 
   3.8012671511164303`*^9}, {3.8012734036885853`*^9, 3.801273423059538*^9}, {
   3.801275191535077*^9, 3.801275225064637*^9}, {3.8012753113665752`*^9, 
   3.801275312173647*^9}, {3.80127746279447*^9, 3.8012774677496223`*^9}, {
   3.8012775148166676`*^9, 3.8012775157468753`*^9}, {3.801354980092315*^9, 
   3.8013549845919433`*^9}, 3.8013554408241215`*^9, {3.8014817866955547`*^9, 
   3.801481786951611*^9}, {3.80155691072915*^9, 3.8015569158790207`*^9}, {
   3.80155703681548*^9, 3.801557053096414*^9}, {3.801557230575369*^9, 
   3.801557258170372*^9}, {3.8018351988553405`*^9, 3.801835199228504*^9}, {
   3.801835251831071*^9, 3.801835253448409*^9}, {3.8018352993836184`*^9, 
   3.801835304281185*^9}, {3.8018435263363857`*^9, 3.801843527311673*^9}, {
   3.801843562268626*^9, 3.801843567225281*^9}, {3.8018604761108923`*^9, 
   3.8018605009090157`*^9}, {3.8018636503801975`*^9, 
   3.8018636924791365`*^9}, {3.801881478763352*^9, 3.8018814914698844`*^9}, {
   3.8018815346716337`*^9, 3.8018815372902107`*^9}, {3.8018816193803368`*^9, 
   3.801881631221342*^9}, {3.801881663856637*^9, 3.8018816826883764`*^9}, {
   3.801884502031747*^9, 3.8018845085764575`*^9}, {3.801885527116788*^9, 
   3.801885549893566*^9}, {3.801885591498884*^9, 3.801885609367434*^9}, {
   3.801886811254795*^9, 3.8018868266043053`*^9}, {3.8018871939628544`*^9, 
   3.80188720622418*^9}, {3.8018874912602634`*^9, 3.801887513696555*^9}, 
   3.8019126044408255`*^9, 3.8019126401020937`*^9, {3.801913703780259*^9, 
   3.8019137131538625`*^9}, {3.801914533535617*^9, 3.8019145465748844`*^9}, {
   3.801917715734974*^9, 3.8019177240653973`*^9}, {3.801918918158475*^9, 
   3.801918928524885*^9}, {3.8019189640064726`*^9, 3.80191896722381*^9}, {
   3.801919447247379*^9, 3.80191945816234*^9}, {3.80191967823343*^9, 
   3.8019196784274735`*^9}, {3.801919988170961*^9, 3.8019199882729845`*^9}, {
   3.801920080560355*^9, 3.801920080632371*^9}, {3.801920161054385*^9, 
   3.8019201630568333`*^9}, {3.80192058888855*^9, 3.801920588971072*^9}, {
   3.8019206919477973`*^9, 3.801920692909019*^9}, {3.801920855965064*^9, 
   3.8019208563286486`*^9}, {3.8019210370980864`*^9, 
   3.8019210372391176`*^9}, {3.8019214132849493`*^9, 3.801921419886386*^9}, 
   3.801921754357629*^9, 3.801921789960856*^9, {3.8019230402960997`*^9, 
   3.8019230460870066`*^9}, {3.801923687451543*^9, 3.8019237274919777`*^9}, {
   3.80192500696103*^9, 3.8019250193164663`*^9}, {3.8019251638999324`*^9, 
   3.8019251799082904`*^9}, {3.8019255208239965`*^9, 3.801925521139083*^9}, {
   3.8019256044624834`*^9, 3.8019256047814846`*^9}, 3.8019266520627556`*^9, {
   3.8019270021560082`*^9, 3.8019270088705397`*^9}, {3.80192810113472*^9, 
   3.801928105640748*^9}, {3.8019287570098066`*^9, 3.801928757096821*^9}, {
   3.8019287875099907`*^9, 3.80192879357106*^9}, {3.801929429615091*^9, 
   3.801929448310198*^9}, {3.801931321865648*^9, 3.8019313219472027`*^9}, {
   3.8019313699247885`*^9, 3.801931375989208*^9}, {3.8019314696501055`*^9, 
   3.8019314782739744`*^9}, {3.8019318193200593`*^9, 
   3.8019318490900807`*^9}, {3.801932472348931*^9, 3.8019324881582108`*^9}, {
   3.8019333927410126`*^9, 3.8019334579708796`*^9}, {3.8019341414483085`*^9, 
   3.8019341490697937`*^9}, {3.801934971747822*^9, 3.801934977863535*^9}, {
   3.8019350587902174`*^9, 3.8019350635305753`*^9}, {3.8019357693104725`*^9, 
   3.801935775399001*^9}, {3.8019361835221834`*^9, 3.8019362099822445`*^9}, {
   3.801938240069064*^9, 3.801938256336527*^9}, {3.801939291495164*^9, 
   3.8019393098834767`*^9}, {3.8019441450109024`*^9, 3.801944169672052*^9}, {
   3.8019446316809015`*^9, 3.8019446317390504`*^9}, {3.8019446855227337`*^9, 
   3.801944708947871*^9}, {3.801944898392271*^9, 3.801944957787821*^9}, {
   3.80194556260907*^9, 3.801945580814294*^9}, {3.8020237728458967`*^9, 
   3.8020238267919226`*^9}, {3.8020250819067717`*^9, 3.802025138044566*^9}, {
   3.8020251828445807`*^9, 3.8020252422776127`*^9}, {3.8020338864454813`*^9, 
   3.802033899292673*^9}, {3.8020340446779575`*^9, 3.8020340450990767`*^9}, {
   3.8020376917131567`*^9, 3.8020376986523886`*^9}, {3.802045550041193*^9, 
   3.8020455799706388`*^9}, {3.8020472326531725`*^9, 
   3.8020472365484123`*^9}, {3.80205583765518*^9, 3.8020558654649405`*^9}, {
   3.802192966157908*^9, 3.802193015278722*^9}, {3.802284200016081*^9, 
   3.802284203163167*^9}, {3.802286297496386*^9, 3.802286302325402*^9}, {
   3.8023075309357986`*^9, 3.8023075371552486`*^9}, {3.802439986721073*^9, 
   3.8024399924906583`*^9}, {3.8024424006494956`*^9, 3.802442405202862*^9}, {
   3.8024613586382003`*^9, 3.8024613641140327`*^9}, {3.802477604052747*^9, 
   3.8024776111179743`*^9}, {3.80247773728697*^9, 3.8024777379797955`*^9}, {
   3.802523413922249*^9, 3.80252341516829*^9}, {3.802597881809531*^9, 
   3.802597888394988*^9}, {3.8027983679197707`*^9, 3.8027983694781218`*^9}, {
   3.802798531496144*^9, 3.802798531574216*^9}, {3.802798608993863*^9, 
   3.802798612616717*^9}, 3.8029707021761255`*^9, {3.8029707416511803`*^9, 
   3.802970744210448*^9}, 3.8029707790708575`*^9, {3.802970809215195*^9, 
   3.8029708385121126`*^9}, {3.8029708701735277`*^9, 3.8029708826790156`*^9},
    {3.8029801451768913`*^9, 3.8029801497714043`*^9}, {3.8029812339240465`*^9,
    3.802981236644396*^9}, {3.8029849225272865`*^9, 3.802984925040943*^9}, {
   3.802985022764078*^9, 3.8029850239433947`*^9}, 3.8029855119921083`*^9, {
   3.802985959388447*^9, 3.8029859595366764`*^9}, {3.8029860418830748`*^9, 
   3.802986041979969*^9}, {3.8029862077133923`*^9, 3.8029862078491993`*^9}, {
   3.802987953765476*^9, 3.80298795749236*^9}, {3.8030082810208764`*^9, 
   3.8030083201661863`*^9}, {3.803082982544694*^9, 3.8030829905236225`*^9}, {
   3.803122351885906*^9, 3.8031223542794704`*^9}, {3.8031311992752495`*^9, 
   3.803131202742576*^9}, {3.803145242067628*^9, 3.803145242503354*^9}, {
   3.8033054136347218`*^9, 3.803305474081072*^9}, 
   3.803305612956951*^9},ExpressionUUID->"7824fddb-3943-46b5-acd8-\
799c34acddf8"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Load force data", "Subsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
  3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
  3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, {
  3.776165362372452*^9, 3.776165368372795*^9}, {3.8012614284506035`*^9, 
  3.80126143387082*^9}},ExpressionUUID->"c541093d-3069-439e-a54a-\
d38289364e3b"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"rawdata", "=", 
   RowBox[{"Import", "[", 
    RowBox[{"fname32", ",", "\"\<CSV\>\"", ",", 
     RowBox[{"HeaderLines", "\[Rule]", "1"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"indices", " ", "=", " ", 
   RowBox[{"DeleteCases", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"If", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"rawdata", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", 
             RowBox[{"1", ";;", "4"}]}], "]"}], "]"}], "!=", 
          RowBox[{"rawdata", "[", 
           RowBox[{"[", 
            RowBox[{
             RowBox[{"i", "+", "1"}], ",", 
             RowBox[{"1", ";;", "4"}]}], "]"}], "]"}]}], ",", "i", ",", 
         "None"}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", 
         RowBox[{
          RowBox[{"Length", "@", "rawdata"}], "-", "1"}]}], "}"}]}], "]"}], 
     ",", "\[IndentingNewLine]", "None"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"Append", "[", 
    RowBox[{
     RowBox[{"rawdata", "[", 
      RowBox[{"[", "indices", "]"}], "]"}], ",", 
     RowBox[{"rawdata", "[", 
      RowBox[{"[", 
       RowBox[{"-", "1"}], "]"}], "]"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"data", " ", "=", " ", 
   RowBox[{"Select", "[", 
    RowBox[{"data", ",", 
     RowBox[{
      RowBox[{
       RowBox[{"#", "[", 
        RowBox[{"[", "4", "]"}], "]"}], "\[NotEqual]", "0"}], "&"}]}], 
    "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"Join", "[", 
    RowBox[{"data", ",", 
     RowBox[{"(", 
      RowBox[{"data", "/.", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
         "a_", ",", "b_", ",", "c_", ",", "d_", ",", "e_", ",", "f_"}], "}"}],
         "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"a", ",", "b", ",", "c", ",", 
          RowBox[{"-", "d"}], ",", 
          RowBox[{"-", "e"}], ",", "f"}], "}"}]}]}], ")"}]}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"Bs", "=", 
   RowBox[{"\[Beta]toB", " ", 
    RowBox[{"data", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", "1"}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"s0s", "=", 
   RowBox[{"data", "[", 
    RowBox[{"[", 
     RowBox[{";;", ",", "2"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[CapitalDelta]s", "=", 
   RowBox[{"data", "[", 
    RowBox[{"[", 
     RowBox[{";;", ",", "3"}], "]"}], "]"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vs", "=", 
   RowBox[{"vtov", " ", 
    RowBox[{"data", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", "4"}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"as", "=", 
   RowBox[{"ftoa", " ", 
    RowBox[{"data", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", "5"}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"rs", "=", 
   RowBox[{"rtor", " ", 
    RowBox[{"data", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", "6"}], "]"}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"ainterp32", " ", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Interpolation", "[", 
      RowBox[{"#", ",", 
       RowBox[{"InterpolationOrder", "\[Rule]", "1"}], ",", 
       RowBox[{"\"\<ExtrapolationHandler\>\"", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.0", "&"}], ")"}], ",", 
          RowBox[{"\"\<WarningMessage\>\"", "\[Rule]", "False"}]}], "}"}]}]}],
       "]"}], "&"}], "@", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Bs", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"s0s", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"\[CapitalDelta]s", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"vs", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"as", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"Length", "@", "data"}]}], "}"}]}], "]"}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"rinterp32", "=", " ", 
   RowBox[{
    RowBox[{
     RowBox[{"Interpolation", "[", 
      RowBox[{"#", ",", 
       RowBox[{"InterpolationOrder", "\[Rule]", "1"}], ",", 
       RowBox[{"\"\<ExtrapolationHandler\>\"", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"(", 
           RowBox[{"0.0", "&"}], ")"}], ",", 
          RowBox[{"\"\<WarningMessage\>\"", "\[Rule]", "False"}]}], "}"}]}]}],
       "]"}], "&"}], "@", 
    RowBox[{"Table", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{
          RowBox[{"Bs", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"s0s", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"\[CapitalDelta]s", "[", 
           RowBox[{"[", "i", "]"}], "]"}], ",", 
          RowBox[{"vs", "[", 
           RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
        RowBox[{"rs", "[", 
         RowBox[{"[", "i", "]"}], "]"}]}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"i", ",", 
        RowBox[{"Length", "@", "data"}]}], "}"}]}], "]"}]}]}], 
  ";"}]}], "Input",
 CellChangeTimes->{{3.778505649547577*^9, 3.778505653284791*^9}, {
   3.7785063666755943`*^9, 3.7785063691387353`*^9}, 3.7947758423822217`*^9, 
   3.7948263883006372`*^9, {3.8019197966919465`*^9, 3.801919800500352*^9}, {
   3.8020252504946136`*^9, 3.8020252515925426`*^9}, {3.8029867876275077`*^9, 
   3.8029867882749166`*^9}, {3.8030917431460185`*^9, 
   3.8030917435971193`*^9}, {3.8030917940765934`*^9, 3.803091794136606*^9}},
 CellLabel->"In[9]:=",ExpressionUUID->"a22cab01-84b6-4910-a8ee-e13c7c9bda3f"],

Cell[CellGroupData[{

Cell["Define function to initialize particles", "Subsubsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
  3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
  3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, {
  3.776165362372452*^9, 3.776165381438542*^9}, {3.7761656187571163`*^9, 
  3.7761656250464754`*^9}, {3.8012752430432863`*^9, 
  3.801275254375457*^9}},ExpressionUUID->"2086d6f6-0d0b-40d9-a62f-\
2851fbd41666"],

Cell[BoxData[
 RowBox[{
  RowBox[{"InitializeParticles", "[", 
   RowBox[{
   "nparticles_", ",", "daperture_", ",", "dcell_", ",", "zaperture_", ",", 
    "vzmin_", ",", "vzavg_", ",", "vzstd_", ",", "vxstd_", ",", " ", "zcool_",
     ",", "xstd_"}], "]"}], ":=", 
  RowBox[{"ParallelTable", "[", 
   RowBox[{
    RowBox[{"Module", "[", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"vzinit", ",", "vxinit", ",", "vyinit", ",", 
        RowBox[{"xinit", "=", 
         RowBox[{"daperture", "+", "dcell", "+", "1"}]}], ",", 
        RowBox[{"yinit", "=", 
         RowBox[{"daperture", "+", "dcell", "+", "1"}]}]}], "}"}], ",", 
      "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"While", "[", 
        RowBox[{
         RowBox[{
          RowBox[{
           RowBox[{"Abs", "[", "xinit", "]"}], " ", ">", " ", 
           RowBox[{"daperture", "/", "2"}]}], " ", "||", " ", 
          RowBox[{
           RowBox[{"Abs", "[", "yinit", "]"}], " ", ">", " ", 
           RowBox[{"daperture", "/", "2"}]}]}], ",", "\[IndentingNewLine]", 
         RowBox[{"(*", 
          RowBox[{"Create", " ", "particles", " ", "at", " ", "cell"}], 
          "*)"}], "\[IndentingNewLine]", 
         RowBox[{
          RowBox[{"vzinit", "=", 
           RowBox[{"Map", "[", 
            RowBox[{
             RowBox[{
              RowBox[{"Max", "[", 
               RowBox[{"vzmin", ",", "#"}], "]"}], "&"}], ",", 
             RowBox[{"RandomVariate", "[", 
              RowBox[{"NormalDistribution", "[", 
               RowBox[{"vzavg", ",", " ", "vzstd"}], "]"}], "]"}]}], "]"}]}], 
          ";", "\[IndentingNewLine]", 
          RowBox[{"vxinit", "=", 
           RowBox[{"RandomVariate", "[", 
            RowBox[{"NormalDistribution", "[", 
             RowBox[{"0", ",", " ", "vxstd"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"vyinit", "=", 
           RowBox[{"RandomVariate", "[", 
            RowBox[{"NormalDistribution", "[", 
             RowBox[{"0", ",", " ", "vxstd"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"xinit", "=", 
           RowBox[{"RandomVariate", "[", 
            RowBox[{"NormalDistribution", "[", 
             RowBox[{"0", ",", " ", "xstd"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"yinit", "=", 
           RowBox[{"RandomVariate", "[", 
            RowBox[{"NormalDistribution", "[", 
             RowBox[{"0", ",", " ", "xstd"}], "]"}], "]"}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"(*", 
           RowBox[{"Propagate", " ", "to", " ", "aperture"}], "*)"}], 
          "\[IndentingNewLine]", 
          RowBox[{"xinit", "+=", 
           RowBox[{"vxinit", "*", 
            RowBox[{"(", 
             RowBox[{"zaperture", "/", "vzinit"}], ")"}]}]}], ";", 
          "\[IndentingNewLine]", 
          RowBox[{"yinit", "+=", 
           RowBox[{"vyinit", "*", 
            RowBox[{"(", 
             RowBox[{"zaperture", "/", "vzinit"}], ")"}]}]}], ";"}]}], 
        "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", 
       RowBox[{"(*", 
        RowBox[{
        "propagate", " ", "to", " ", "start", " ", "of", " ", "cooling", " ", 
         "region"}], "*)"}], "\[IndentingNewLine]", 
       RowBox[{"xinit", "+=", 
        RowBox[{"vxinit", " ", 
         RowBox[{"zcool", "/", "vzinit"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"yinit", "+=", 
        RowBox[{"vyinit", " ", 
         RowBox[{"zcool", "/", "vzinit"}]}]}], ";", "\[IndentingNewLine]", 
       RowBox[{"{", 
        RowBox[{
        "xinit", ",", "vxinit", ",", "vzinit", ",", "yinit", ",", "vyinit"}], 
        "}"}]}]}], "]"}], ",", "\[IndentingNewLine]", 
    RowBox[{"{", 
     RowBox[{"i", ",", "nparticles"}], "}"}], ",", 
    RowBox[{"Method", "\[Rule]", "\"\<CoarsestGrained\>\""}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.8020534368028145`*^9, 3.8020534564815226`*^9}, {
  3.80205354327876*^9, 3.802053544925171*^9}},
 CellLabel->"In[22]:=",ExpressionUUID->"faefb1d2-c143-4cf4-bbad-8f8982efcabd"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Define function to simulate trajectories", "Subsubsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
   3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
   3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, {
   3.776165362372452*^9, 3.776165395656355*^9}, 3.776165630526789*^9, {
   3.794530301368927*^9, 3.7945303021651297`*^9}, {3.8012752498082466`*^9, 
   3.8012752521148233`*^9}},ExpressionUUID->"20effcbb-203d-40b4-8029-\
107270944d9f"],

Cell[BoxData[
 RowBox[{
  RowBox[{"SimulateTrajectoryFast", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"x0_", ",", "vx0_", ",", "vz_", ",", "y0_", ",", "vy_"}], "}"}], 
    ",", "\[Beta]_", ",", "p_", ",", "\[CapitalDelta]_", ",", 
    "nscatteravg_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"eqns", ",", "tmax", ",", "tend", ",", "sol", ",", "nscatter"}], 
     "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"nscatter", "=", 
      RowBox[{"RandomVariate", "[", 
       RowBox[{"ExponentialDistribution", "[", 
        FractionBox["1", "nscatteravg"], "]"}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"tmax", "=", 
      RowBox[{"lcool", "/", "vz"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"decayed", "=", "0"}], ";", "\[IndentingNewLine]", 
     RowBox[{"eqns", "=", 
      RowBox[{"{", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{
         RowBox[{"vx", "[", "0", "]"}], "\[Equal]", "vx0"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"x", "[", "0", "]"}], "\[Equal]", "x0"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"n", "[", "0", "]"}], "\[Equal]", "0"}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"vx", "'"}], "[", "t", "]"}], "\[Equal]", 
         RowBox[{"ainterp", "[", 
          RowBox[{"\[Beta]", ",", 
           RowBox[{"p", " ", 
            RowBox[{"intensity", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"vz", " ", "t"}], "-", 
               RowBox[{"lcool", "/", "2"}]}], ",", " ", 
              RowBox[{"y0", "+", 
               RowBox[{"vy", " ", "t"}]}]}], "]"}]}], ",", "\[CapitalDelta]", 
           ",", 
           RowBox[{"vx", "[", "t", "]"}]}], "]"}]}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"x", "'"}], "[", "t", "]"}], "\[Equal]", 
         RowBox[{"vx", "[", "t", "]"}]}], ",", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{
          RowBox[{"n", "'"}], "[", "t", "]"}], "\[Equal]", 
         RowBox[{"rinterp", "[", 
          RowBox[{"\[Beta]", ",", 
           RowBox[{"p", " ", 
            RowBox[{"intensity", "[", 
             RowBox[{
              RowBox[{
               RowBox[{"vz", " ", "t"}], "-", 
               RowBox[{"lcool", "/", "2"}]}], ",", " ", 
              RowBox[{"y0", "+", 
               RowBox[{"vy", " ", "t"}]}]}], "]"}]}], ",", "\[CapitalDelta]", 
           ",", 
           RowBox[{"vx", "[", "t", "]"}]}], "]"}]}], ",", 
        "\[IndentingNewLine]", 
        RowBox[{"WhenEvent", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"n", "[", "t", "]"}], ">", "nscatter"}], ",", 
          "\"\<StopIntegration\>\""}], "]"}], ",", "\[IndentingNewLine]", 
        RowBox[{"WhenEvent", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"Mod", "[", 
            RowBox[{
             RowBox[{"n", "[", "t", "]"}], ",", "1"}], "]"}], ">", "0.5"}], 
          ",", 
          RowBox[{
           RowBox[{"vx", "[", "t", "]"}], "\[Rule]", 
           RowBox[{
            RowBox[{"vx", "[", "t", "]"}], "+", 
            RowBox[{"RandomReal", "[", 
             RowBox[{"{", 
              RowBox[{
               RowBox[{"-", "vrecoil"}], ",", "vrecoil"}], "}"}], "]"}]}]}]}],
          "]"}], ",", "\[IndentingNewLine]", 
        RowBox[{"WhenEvent", "[", 
         RowBox[{
          RowBox[{
           RowBox[{"n", "[", "t", "]"}], ">", "nscatter"}], ",", " ", 
          RowBox[{"decayed", "=", "1"}]}], "]"}]}], "\[IndentingNewLine]", 
       "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"sol", "=", 
      RowBox[{"First", "@", 
       RowBox[{"NDSolve", "[", 
        RowBox[{"eqns", ",", 
         RowBox[{"{", 
          RowBox[{"x", ",", "vx", ",", "n"}], "}"}], ",", 
         RowBox[{"{", 
          RowBox[{"t", ",", "0", ",", "tmax"}], "}"}], ",", 
         "\[IndentingNewLine]", 
         RowBox[{"AccuracyGoal", "\[Rule]", "5"}], ",", "\[IndentingNewLine]", 
         RowBox[{"PrecisionGoal", "\[Rule]", "3"}]}], "]"}]}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"tend", "=", 
      RowBox[{
       RowBox[{"(", 
        RowBox[{
         RowBox[{"x", "[", "\"\<Domain\>\"", "]"}], "/.", "sol"}], ")"}], "[", 
       RowBox[{"[", 
        RowBox[{"1", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{
        RowBox[{
         RowBox[{"x", "[", "tend", "]"}], "+", 
         RowBox[{
          RowBox[{"vx", "[", "tend", "]"}], "*", 
          RowBox[{"(", 
           RowBox[{"tmax", " ", "-", " ", "tend", " ", "+", " ", 
            RowBox[{"zdetector", "/", "vz"}]}], ")"}]}]}], ",", 
        RowBox[{"vx", "[", "tend", "]"}], ",", 
        RowBox[{"tmax", "-", "tend"}], ",", 
        RowBox[{"n", "[", "tend", "]"}], ",", "decayed", ",", 
        RowBox[{"x", "[", "0", "]"}], ",", 
        RowBox[{"vx", "[", "0", "]"}], ",", "y0", ",", "vy", ",", "vz"}], 
       "}"}], "/.", "sol"}]}]}], "\[IndentingNewLine]", "]"}]}]], "Input",
 CellChangeTimes->{{3.776174916418912*^9, 3.7761749431704416`*^9}, {
   3.7761753432633257`*^9, 3.776175343719352*^9}, {3.7761756976875973`*^9, 
   3.776175788865813*^9}, {3.7761758305431967`*^9, 3.776175863875103*^9}, {
   3.77617589861709*^9, 3.7761760192399893`*^9}, {3.776176176268971*^9, 
   3.7761762161692533`*^9}, {3.7761762609278135`*^9, 
   3.7761763252724934`*^9}, {3.776176423677122*^9, 3.776176429166436*^9}, {
   3.77617648990191*^9, 3.776176550190358*^9}, {3.77617658728448*^9, 
   3.7761766392574525`*^9}, {3.7761767284945564`*^9, 3.776176728813575*^9}, {
   3.7761767726930847`*^9, 3.776176835595682*^9}, {3.776176887459649*^9, 
   3.7761770346530676`*^9}, {3.7761771013998857`*^9, 3.776177125273251*^9}, {
   3.776177160488265*^9, 3.776177291839778*^9}, {3.776177416344899*^9, 
   3.776177424632373*^9}, {3.776178850516929*^9, 3.7761788921623116`*^9}, {
   3.776179057552771*^9, 3.7761790611119747`*^9}, {3.776179320394805*^9, 
   3.7761793664154367`*^9}, {3.776180333439748*^9, 3.7761803850937023`*^9}, {
   3.7761817269974546`*^9, 3.7761817663397045`*^9}, {3.7761818016497245`*^9, 
   3.7761818803122234`*^9}, {3.776182115720688*^9, 3.776182115839695*^9}, 
   3.776182195547254*^9, {3.7764541730239077`*^9, 3.776454251541815*^9}, {
   3.776454348026761*^9, 3.776454361679059*^9}, {3.776513553854104*^9, 
   3.7765135795985765`*^9}, {3.776683726294552*^9, 3.776683786049969*^9}, {
   3.776683962329052*^9, 3.776683969601468*^9}, {3.7766841501127925`*^9, 
   3.7766841599063525`*^9}, {3.776684714749088*^9, 3.776684741153598*^9}, {
   3.7771128820607176`*^9, 3.7771128872600155`*^9}, {3.777112920601922*^9, 
   3.7771129450683217`*^9}, {3.777123850644085*^9, 3.7771238832469497`*^9}, {
   3.777123924196292*^9, 3.7771239434173913`*^9}, {3.7771240840254335`*^9, 
   3.777124089256733*^9}, {3.777124204060299*^9, 3.7771242047783403`*^9}, {
   3.7771243989394455`*^9, 3.7771244157594075`*^9}, {3.7771246277135305`*^9, 
   3.7771246366570425`*^9}, {3.777124728838315*^9, 3.77712473172548*^9}, {
   3.7771248359814425`*^9, 3.7771248379875574`*^9}, {3.7772257683142014`*^9, 
   3.777225769659278*^9}, {3.7772258175640182`*^9, 3.7772258178970375`*^9}, {
   3.7772258898701534`*^9, 3.7772259413791*^9}, {3.777225991471965*^9, 
   3.777225993119059*^9}, {3.7772261294918594`*^9, 3.7772261556693563`*^9}, {
   3.777226351352549*^9, 3.7772264001163383`*^9}, {3.7773713668619614`*^9, 
   3.777371380557745*^9}, {3.777378502509097*^9, 3.777378530791715*^9}, {
   3.7773885607639847`*^9, 3.7773885869104805`*^9}, {3.7774602552878046`*^9, 
   3.777460344434904*^9}, {3.777460960595146*^9, 3.77746100123147*^9}, {
   3.7774610357414446`*^9, 3.777461074529663*^9}, {3.777553158625386*^9, 
   3.777553183258795*^9}, {3.778235013344737*^9, 3.77823501410914*^9}, {
   3.7782350693207827`*^9, 3.778235246047667*^9}, {3.7782356047904577`*^9, 
   3.7782356857981052`*^9}, {3.778235729419179*^9, 3.778235740906472*^9}, {
   3.7782357968991423`*^9, 3.7782358751388407`*^9}, {3.7782359642721443`*^9, 
   3.778236103984399*^9}, {3.778236216751554*^9, 3.7782362746022253`*^9}, {
   3.7785080714224567`*^9, 3.778508071629874*^9}, {3.7785812997398295`*^9, 
   3.778581304077222*^9}, {3.7785852223640146`*^9, 3.7785852225355225`*^9}, {
   3.7785856318281918`*^9, 3.7785856830401773`*^9}, {3.778585729782076*^9, 
   3.778585750770913*^9}, {3.7785858968280396`*^9, 3.778585917373057*^9}, {
   3.794398701583085*^9, 3.794398712997385*^9}, {3.794398817369732*^9, 
   3.794398840328929*^9}, {3.7943999492219834`*^9, 3.7943999538130655`*^9}, {
   3.7943999858247557`*^9, 3.794399986643861*^9}, {3.79440178016883*^9, 
   3.794401839307207*^9}, {3.7944022037112727`*^9, 3.7944022194774265`*^9}, {
   3.794402383940595*^9, 3.794402384172909*^9}, {3.7944865777073936`*^9, 
   3.7944865841324525`*^9}, {3.794487023673592*^9, 3.7944870248833704`*^9}, 
   3.7946779192618337`*^9, {3.7946927033924246`*^9, 3.7946927118438377`*^9}, {
   3.7947936150104666`*^9, 3.794793616164729*^9}, {3.8012754530356846`*^9, 
   3.801275453869988*^9}, {3.801299777839184*^9, 3.8012998056420946`*^9}, {
   3.801557108083917*^9, 3.8015571090401306`*^9}, 3.801557525411193*^9, {
   3.8015575822587523`*^9, 3.801557582541818*^9}, {3.8018854929418335`*^9, 
   3.8018854988534136`*^9}, {3.801885644081657*^9, 3.80188564416286*^9}, {
   3.802045571001005*^9, 3.8020455714890833`*^9}, {3.8020867031762295`*^9, 
   3.802086716510292*^9}, {3.803091796352104*^9, 3.803091796419119*^9}, {
   3.8031452938197436`*^9, 3.8031452939806395`*^9}},
 CellLabel->"In[44]:=",ExpressionUUID->"dccf9f9d-0271-4cee-a2a6-51b034d09fcc"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Initialize beam with experimental parameters", "Subsubsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
   3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
   3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, {
   3.776165362372452*^9, 3.776165395656355*^9}, 3.776165630526789*^9, {
   3.794530301368927*^9, 3.7945303021651297`*^9}, {3.8012752498082466`*^9, 
   3.801275295145632*^9}, {3.801275345493616*^9, 
   3.8012753500168557`*^9}},ExpressionUUID->"6c519e3a-f10c-4538-afc2-\
24105e47eb71"],

Cell[BoxData[
 RowBox[{
  RowBox[{"nparticles", "=", "10000"}], ";"}]], "Input",
 CellChangeTimes->{{3.7944169993043137`*^9, 3.794417003826792*^9}, {
   3.7946539606021285`*^9, 3.7946539610622315`*^9}, {3.794657794485712*^9, 
   3.7946578025453887`*^9}, {3.79465906243653*^9, 3.7946590626755843`*^9}, 
   3.794660235224152*^9, {3.7946779111817827`*^9, 3.7946779164801683`*^9}, {
   3.7946854581590242`*^9, 3.794685469409958*^9}, 3.794688410701245*^9, {
   3.794694316618695*^9, 3.7946943181510377`*^9}, {3.7946970506091723`*^9, 
   3.794697050681189*^9}, {3.79477585523313*^9, 3.7947758587915897`*^9}, {
   3.794798994054731*^9, 3.7947989945407934`*^9}, {3.794837864660307*^9, 
   3.7948378661526413`*^9}, {3.794938980254684*^9, 3.794938980310358*^9}, {
   3.794939356136031*^9, 3.7949393564931107`*^9}, {3.794945721548031*^9, 
   3.7949457218564386`*^9}, {3.7949558349650292`*^9, 
   3.7949558350520487`*^9}, {3.7949577288327723`*^9, 3.794957734302003*^9}, {
   3.794962856825825*^9, 3.7949628571042256`*^9}, 3.8012641961977863`*^9, {
   3.801271432973773*^9, 3.801271433071795*^9}, {3.801273458376217*^9, 
   3.8012734584772396`*^9}, {3.801275141422189*^9, 3.8012751414972415`*^9}, {
   3.801275265309698*^9, 3.801275265618518*^9}, {3.801276120940916*^9, 
   3.8012761212788353`*^9}, {3.8013461928857665`*^9, 
   3.8013461929787874`*^9}, {3.801355521808771*^9, 3.8013555218807735`*^9}, {
   3.8014252159236016`*^9, 3.8014252161426506`*^9}, {3.801560314350685*^9, 
   3.801560314454709*^9}, {3.801600195044063*^9, 3.801600195104892*^9}, {
   3.801600500977229*^9, 3.801600501028126*^9}, {3.801600757940977*^9, 
   3.801600758018566*^9}, {3.801601077941332*^9, 3.80160107802839*^9}, {
   3.8016011647230873`*^9, 3.80160116478224*^9}, {3.801822691374462*^9, 
   3.8018226914154596`*^9}, {3.801833157849103*^9, 3.8018331579576316`*^9}, {
   3.8018384581285086`*^9, 3.801838458645624*^9}, {3.801838861356563*^9, 
   3.80183886143058*^9}, {3.801912403662136*^9, 3.8019124073968315`*^9}, {
   3.8019143277007885`*^9, 3.8019143278248158`*^9}, {3.801940350953227*^9, 
   3.8019403509978724`*^9}, {3.8019408059768066`*^9, 3.801940806027597*^9}, {
   3.8019412168526597`*^9, 3.8019412169567757`*^9}, {3.8019419731039495`*^9, 
   3.801941973230673*^9}, {3.8020534868364706`*^9, 3.8020534869669247`*^9}, {
   3.8020538917116027`*^9, 3.802053891788604*^9}, {3.802054279700254*^9, 
   3.802054279874982*^9}, {3.802054344828485*^9, 3.802054344921379*^9}, {
   3.8020886898661337`*^9, 3.802088689919146*^9}, {3.8020890297853622`*^9, 
   3.8020890302554665`*^9}, {3.8024756629481645`*^9, 3.802475663135049*^9}, {
   3.803026581539429*^9, 3.8030265816434526`*^9}, {3.8031452595289226`*^9, 
   3.8031452595907764`*^9}},
 CellLabel->"In[24]:=",ExpressionUUID->"9e84a53b-6a94-415f-8b47-20c1efb51727"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"constants", " ", "parameterizing", " ", "the", " ", 
    RowBox[{"experiment", ".", "\[IndentingNewLine]", "The"}], " ", "beam", 
    " ", "propagates", " ", "along", " ", "the", " ", "z", " ", 
    RowBox[{"axis", ".", "The"}], " ", "x", " ", "axis", " ", "is", " ", 
    "the", " ", "transverse", " ", 
    RowBox[{"direction", "."}]}], " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dcell", "=", "0.0115"}], ";"}], 
   RowBox[{"(*", " ", 
    RowBox[{
     RowBox[{"Width", " ", "of", " ", "cell", " ", "exit"}], ";", "m"}], 
    "*)"}], "\n", 
   RowBox[{
    RowBox[{"vzavg", "=", "150.0"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Average", " ", "forward", " ", "velocity"}], ";", 
     RowBox[{"m", "/", "s"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"vzstd", "=", "30.0"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Standard", " ", "deviation", " ", "of", " ", "forward", " ", 
      "velocity"}], ";", 
     RowBox[{"m", "/", "s"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"vzmin", "=", "10.0"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Minimum", " ", "forward", " ", "velocity"}], ";", 
     RowBox[{"m", "/", "s"}]}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"vxstd", "=", "20"}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Standard", " ", "deviation", " ", "of", " ", "transverse", " ", 
      "velocity"}], ";", 
     RowBox[{"m", "/", "s"}]}], "*)"}], "\n", 
   RowBox[{
    RowBox[{"zaperture", "=", "0.355"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Distance", " ", "from", " ", "cell", " ", "to", " ", "aperture"}], ";", 
     "m"}], "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"daperture", "=", "0.003"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Diameter", " ", "of", " ", "aperture"}], ";", "m"}], "*)"}], 
   "\n", 
   RowBox[{
    RowBox[{"zcool", "=", "0.02"}], ";"}], " ", 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{"Length", " ", "of", " ", "cooling", " ", "region"}], ";", "m"}],
     "*)"}], "\n", 
   RowBox[{
    RowBox[{"lcool", "=", "0.01"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"zdetector", "=", "0.515"}], ";"}], 
   RowBox[{"(*", 
    RowBox[{
     RowBox[{
     "Distance", " ", "between", " ", "end", " ", "of", " ", "cooling", " ", 
      "region", " ", "and", " ", "detector"}], ";", " ", "m"}], "*)"}], 
   "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"rdetector", "=", "0.015"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"xstd", "=", "0.0038"}], ";"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"init", "=", 
     RowBox[{"InitializeParticles", "[", 
      RowBox[{
      "nparticles", ",", "daperture", ",", "dcell", ",", "zaperture", ",", 
       "vzmin", ",", "vzavg", ",", "vzstd", ",", "vxstd", ",", " ", "zcool", 
       ",", "xstd"}], "]"}]}], ";"}]}]}]], "Input",
 CellChangeTimes->{{3.7944038032008038`*^9, 3.794403827910453*^9}, {
   3.794403920340012*^9, 3.794403946682702*^9}, 3.7944039801799545`*^9, {
   3.794404078790533*^9, 3.794404085731656*^9}, {3.7944041379631567`*^9, 
   3.7944041509624677`*^9}, {3.794404224780773*^9, 3.7944042252380433`*^9}, {
   3.7944042806072035`*^9, 3.7944043037506313`*^9}, {3.794404355432145*^9, 
   3.7944043694529595`*^9}, {3.794404408351403*^9, 3.7944044111938033`*^9}, {
   3.7944044939003863`*^9, 3.7944044940756693`*^9}, {3.7944047408531065`*^9, 
   3.7944047422674427`*^9}, {3.794404777212698*^9, 3.7944047957653213`*^9}, {
   3.794404832928384*^9, 3.794404833193509*^9}, {3.794404869277589*^9, 
   3.7944048709071293`*^9}, {3.7944049223486595`*^9, 3.79440495658934*^9}, {
   3.794405183883483*^9, 3.794405183971407*^9}, {3.7944052288872166`*^9, 
   3.794405232583681*^9}, {3.794405440700643*^9, 3.794405440996295*^9}, {
   3.7944055240793505`*^9, 3.7944055249344473`*^9}, {3.794405576345329*^9, 
   3.794405576569295*^9}, {3.794405842103505*^9, 3.7944058541753097`*^9}, {
   3.79440589429871*^9, 3.794405895752221*^9}, {3.7944550459902945`*^9, 
   3.794455046060461*^9}, {3.794455096943061*^9, 3.794455097074314*^9}, {
   3.794455153145345*^9, 3.794455153320022*^9}, {3.794456061762029*^9, 
   3.7944560658981037`*^9}, {3.794491127841591*^9, 3.794491128653432*^9}, {
   3.794491633521778*^9, 3.794491633672577*^9}, {3.794653969111059*^9, 
   3.7946539691740737`*^9}, {3.7946673860088587`*^9, 
   3.7946673867930346`*^9}, {3.7946773687020636`*^9, 3.794677369117157*^9}, {
   3.794678291750182*^9, 3.794678291805195*^9}, 3.794678471016124*^9, {
   3.7946854384500413`*^9, 3.7946854527367783`*^9}, {3.794685908536018*^9, 
   3.794685909056135*^9}, {3.794686442002347*^9, 3.794686442520463*^9}, {
   3.794686849278145*^9, 3.7946868494861913`*^9}, {3.794687923601124*^9, 
   3.7946879236811414`*^9}, {3.7946884043512955`*^9, 3.794688406141698*^9}, {
   3.794690097667096*^9, 3.7946900986733217`*^9}, {3.7946913415571527`*^9, 
   3.794691347645033*^9}, {3.7946928054556594`*^9, 3.7946928064653893`*^9}, {
   3.7946931126780252`*^9, 3.7946931127500424`*^9}, {3.794693294290248*^9, 
   3.794693298174621*^9}, {3.7946943290981417`*^9, 3.7946943309806695`*^9}, {
   3.794694613098367*^9, 3.7946946159820213`*^9}, {3.7946946554204845`*^9, 
   3.7946946640011063`*^9}, {3.7947438292197638`*^9, 
   3.7947438292767773`*^9}, {3.7947990036830077`*^9, 3.794799003766839*^9}, {
   3.79482641928953*^9, 3.7948264197946434`*^9}, {3.7948365964156995`*^9, 
   3.794836596770779*^9}, {3.7948369100423627`*^9, 3.7948369101593885`*^9}, {
   3.7948548981224847`*^9, 3.7948548988553286`*^9}, {3.7949553878862762`*^9, 
   3.794955388329376*^9}, {3.7949638986837835`*^9, 3.7949639108986216`*^9}, {
   3.79496403222542*^9, 3.794964032595701*^9}, {3.7949670181691747`*^9, 
   3.7949670186502824`*^9}, {3.7949679922297387`*^9, 3.794967992892903*^9}, {
   3.7950160730758753`*^9, 3.79501607329093*^9}, {3.7950326637216578`*^9, 
   3.795032664504648*^9}, {3.8012606790717616`*^9, 3.8012606801910114`*^9}, 
   3.8012680957025213`*^9, {3.8012725782624454`*^9, 3.8012725783411307`*^9}, {
   3.80127346308727*^9, 3.8012734633513284`*^9}, {3.801274317945307*^9, 
   3.801274318569776*^9}, {3.8012761097027426`*^9, 3.8012761146139603`*^9}, {
   3.8012767042551403`*^9, 3.80127670438795*^9}, {3.801277395147415*^9, 
   3.801277395271443*^9}, {3.8012997747034316`*^9, 3.8012997748884726`*^9}, {
   3.8015591841791115`*^9, 3.8015591845361905`*^9}, {3.801838575083886*^9, 
   3.8018385764932175`*^9}, {3.8018391223894973`*^9, 
   3.8018391240368752`*^9}, {3.8018394101556597`*^9, 3.801839411915064*^9}, {
   3.801839892245731*^9, 3.801839893365906*^9}, {3.8018848100385623`*^9, 
   3.801884810168342*^9}, {3.801886136258749*^9, 3.801886136574029*^9}, {
   3.8019148627084208`*^9, 3.801914862864291*^9}, {3.8019304473360925`*^9, 
   3.801930447407585*^9}, {3.801939887274111*^9, 3.80193988787298*^9}, 
   3.8019400191077967`*^9, {3.8019402462791214`*^9, 3.801940252062274*^9}, {
   3.8019403003070107`*^9, 3.801940301300228*^9}, {3.80194046219112*^9, 
   3.8019404731366053`*^9}, {3.8019405275227976`*^9, 
   3.8019405279588947`*^9}, {3.801940797522706*^9, 3.8019408236412416`*^9}, {
   3.801941133025793*^9, 3.801941134554677*^9}, {3.801941292147977*^9, 
   3.801941292615161*^9}, {3.801941385485051*^9, 3.801941386295328*^9}, {
   3.801941960145096*^9, 3.8019419622646847`*^9}, 3.801942518635738*^9, {
   3.8019429725367956`*^9, 3.8019429726307974`*^9}, {3.8020261864534645`*^9, 
   3.802026186535722*^9}, {3.802026674478385*^9, 3.8020266747767267`*^9}, {
   3.802045898361003*^9, 3.802045898575239*^9}, {3.8020472616061616`*^9, 
   3.8020472616671653`*^9}, {3.8020493285148687`*^9, 3.802049328867948*^9}, 
   3.8020519498583937`*^9, 3.802052398002857*^9, {3.802053038466491*^9, 
   3.8020530412232695`*^9}, {3.802053290489812*^9, 3.8020532913135977`*^9}, {
   3.80205346624981*^9, 3.802053478849326*^9}, {3.8020536334883366`*^9, 
   3.8020536426363773`*^9}, {3.8020537285797825`*^9, 
   3.8020537292195463`*^9}, {3.8020538322847013`*^9, 3.802053832622427*^9}, {
   3.8020542492447386`*^9, 3.8020542834232807`*^9}, {3.802087341774027*^9, 
   3.802087342087096*^9}, {3.8020875816186676`*^9, 3.8020875818797264`*^9}, {
   3.802088672929885*^9, 3.8020886752659187`*^9}, {3.802088897723119*^9, 
   3.8020889008148165`*^9}, {3.8020892766623883`*^9, 
   3.8020892769714575`*^9}, {3.80212753404093*^9, 3.8021275343347106`*^9}, {
   3.802284650549883*^9, 3.8022846506366453`*^9}, {3.802285832044445*^9, 
   3.8022858357017956`*^9}, {3.8023095570312862`*^9, 3.802309583520257*^9}, {
   3.8024020835984106`*^9, 3.802402084013175*^9}, {3.802461431754057*^9, 
   3.802461431844738*^9}, {3.8024743796061964`*^9, 3.802474379933423*^9}, {
   3.8024747308862267`*^9, 3.8024747312400513`*^9}, {3.802475030374666*^9, 
   3.8024750304231825`*^9}, {3.8024753580340123`*^9, 3.802475358879683*^9}, {
   3.802475583598369*^9, 3.802475583955789*^9}, {3.8024756417380943`*^9, 
   3.8024756418228407`*^9}, {3.8024759253671017`*^9, 3.802475943224651*^9}, {
   3.8024773886867657`*^9, 3.802477390786579*^9}, {3.8024774311303773`*^9, 
   3.802477431221923*^9}, {3.802479544893146*^9, 3.802479547578958*^9}, {
   3.802604102479721*^9, 3.8026041028142214`*^9}},
 CellLabel->"In[25]:=",ExpressionUUID->"c0c8e4fa-0227-4b9a-82f2-c524254ed188"],

Cell[BoxData[
 RowBox[{
  RowBox[{"\[CapitalGamma]", "=", 
   RowBox[{
    RowBox[{"N", "[", 
     RowBox[{"QuantityMagnitude", "[", 
      RowBox[{"1", "/", 
       RowBox[{"(", 
        RowBox[{"2", "\[Pi]", "*", "tdecay"}], ")"}]}], "]"}], "]"}], "*", 
    SuperscriptBox["10", "3"]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.801822763882579*^9, 3.801822876113387*^9}, {
   3.801845561644772*^9, 3.801845562344228*^9}, {3.80184571561304*^9, 
   3.801845716465281*^9}, {3.801873180440552*^9, 3.8018732154812937`*^9}, {
   3.8019154827791853`*^9, 3.801915505830362*^9}, 3.8019158582058563`*^9, {
   3.8019159284595547`*^9, 3.801915929195719*^9}, 3.801918980582033*^9, 
   3.8019190160605307`*^9, 3.801920372204812*^9, {3.801922642943591*^9, 
   3.801922642986601*^9}, 3.8019229636286483`*^9, {3.80193092265907*^9, 
   3.80193093436341*^9}, {3.801931330950056*^9, 3.8019313322610683`*^9}, {
   3.8019335017452545`*^9, 3.801933532921264*^9}, {3.8019349606493587`*^9, 
   3.801934962904294*^9}, {3.8019350367082405`*^9, 3.801935055018607*^9}, {
   3.80193543362868*^9, 3.8019354340935073`*^9}, {3.8019357597629538`*^9, 
   3.801935760448981*^9}, {3.8019370461516795`*^9, 3.801937046274697*^9}, {
   3.801937101702786*^9, 3.8019371025757465`*^9}, {3.801944887627026*^9, 
   3.8019448877695875`*^9}, {3.802036519635071*^9, 3.8020365197130876`*^9}, {
   3.802036583931672*^9, 3.802036595046625*^9}, {3.802037735457374*^9, 
   3.802037735692346*^9}, {3.802047269724803*^9, 3.8020472889897685`*^9}, {
   3.8020493051448393`*^9, 3.8020493058260465`*^9}, {3.802049382699054*^9, 
   3.8020493828726735`*^9}, 3.802055829723775*^9, {3.8020596030370884`*^9, 
   3.802059603317359*^9}, {3.802060350816439*^9, 3.8020603508654504`*^9}, {
   3.802086732244137*^9, 3.802086732334158*^9}, {3.802088372698949*^9, 
   3.802088374003744*^9}, {3.8020900987251835`*^9, 3.8020901193156157`*^9}, 
   3.802189248385887*^9, {3.802284598193411*^9, 3.802284604689617*^9}, 
   3.8022848864490767`*^9, {3.8022851532065525`*^9, 3.802285168992221*^9}, {
   3.80228531136736*^9, 3.802285321090247*^9}, {3.8022856183095093`*^9, 
   3.8022856183863735`*^9}, {3.802301699454156*^9, 3.802301699592894*^9}, {
   3.80230451109772*^9, 3.802304511205656*^9}, 3.8023047924645724`*^9, {
   3.8023082652772846`*^9, 3.8023082653668337`*^9}, 3.8023083401535454`*^9, {
   3.8023778109036283`*^9, 3.802377811057498*^9}, {3.8024435815348577`*^9, 
   3.8024435896235*^9}, {3.8024741396470747`*^9, 3.8024741398179846`*^9}, {
   3.802477948523782*^9, 3.80247794903541*^9}, {3.8025983643086443`*^9, 
   3.8025983751662407`*^9}, {3.802599101998064*^9, 3.802599102197171*^9}, {
   3.80259929174094*^9, 3.8025992918879814`*^9}, {3.802599404180417*^9, 
   3.8025994074463634`*^9}, {3.802971263178914*^9, 3.802971263248232*^9}, {
   3.8029724736979327`*^9, 3.8029724737323055`*^9}, {3.802980169041366*^9, 
   3.8029801697562375`*^9}, {3.8029803825391197`*^9, 3.8029803830082784`*^9}, 
   3.8030829965576315`*^9},
 CellLabel->"In[38]:=",ExpressionUUID->"c5609752-0026-4bbc-8872-3b54364462a9"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", 
   RowBox[{"visual", " ", "check", " ", "of", " ", "force", " ", "profiles"}],
    "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"curves", " ", "=", 
     RowBox[{"ainterp", "[", 
      RowBox[{
       RowBox[{"-", "5"}], ",", 
       RowBox[{"QuantityMagnitude", "[", 
        RowBox[{"2", "*", 
         RowBox[{"800", "/", "isat"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "15"}], "/", "\[CapitalGamma]"}], ",", "0", ",", 
         RowBox[{"20", "/", "\[CapitalGamma]"}]}], "}"}], ",", "v"}], "]"}]}],
     ";"}], "\[IndentingNewLine]", 
   RowBox[{"Plot", "[", 
    RowBox[{"curves", ",", 
     RowBox[{"{", 
      RowBox[{"v", ",", 
       RowBox[{"-", "15"}], ",", "15"}], "}"}], ",", 
     RowBox[{"ImageSize", "\[Rule]", "Large"}], ",", 
     RowBox[{"PlotLegends", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
       "\"\<Heating\>\"", ",", "\"\<Resonant\>\"", ",", "\"\<Cooling\>\""}], 
       "}"}]}], ",", 
     RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]}]}]], "Input",
 CellChangeTimes->{{3.7944292168794184`*^9, 3.79442924734289*^9}, {
   3.7944294433179064`*^9, 3.794429444262847*^9}, {3.794429667520056*^9, 
   3.794429672372655*^9}, {3.7944297318387547`*^9, 3.794429789989504*^9}, {
   3.7944299590211124`*^9, 3.794429970832421*^9}, {3.794430010003907*^9, 
   3.794430033282007*^9}, {3.7944300885541506`*^9, 3.7944300895081315`*^9}, {
   3.7944301211262608`*^9, 3.794430122018141*^9}, {3.7944305708874083`*^9, 
   3.794430571692419*^9}, {3.7944312263760743`*^9, 3.79443122648833*^9}, {
   3.794431642360224*^9, 3.794431666895884*^9}, {3.794431755988937*^9, 
   3.7944317722270985`*^9}, {3.7944319768596716`*^9, 3.794431990160483*^9}, {
   3.7944320697057595`*^9, 3.7944321446200323`*^9}, {3.7944322115537214`*^9, 
   3.7944322133682013`*^9}, {3.794432243919973*^9, 3.794432246465009*^9}, {
   3.7944325638588467`*^9, 3.794432563947193*^9}, {3.7944327523300037`*^9, 
   3.7944328056614704`*^9}, {3.79443286552392*^9, 3.794432866813369*^9}, {
   3.794432980778113*^9, 3.7944329876315737`*^9}, {3.7944331662964187`*^9, 
   3.7944332261699977`*^9}, {3.794433562147122*^9, 3.7944335960883217`*^9}, {
   3.7944337507398605`*^9, 3.79443377420704*^9}, {3.7944339388827686`*^9, 
   3.7944339437564554`*^9}, {3.7944340044224443`*^9, 
   3.7944340208546658`*^9}, {3.7944343013561416`*^9, 
   3.7944343974877515`*^9}, {3.7944345273616743`*^9, 
   3.7944345345377007`*^9}, {3.794434721147224*^9, 3.794434786608362*^9}, {
   3.7944352213466597`*^9, 3.794435228832205*^9}, {3.794435368796913*^9, 
   3.7944353826397123`*^9}, {3.7944354134107738`*^9, 
   3.7944354181767097`*^9}, {3.7944370374387865`*^9, 
   3.7944370413617535`*^9}, {3.7944399400360465`*^9, 
   3.7944399401140633`*^9}, {3.794440127841544*^9, 3.794440168312812*^9}, {
   3.794440364473666*^9, 3.7944404524727025`*^9}, {3.7944405919719195`*^9, 
   3.7944406509609947`*^9}, {3.7944406864157286`*^9, 3.794440706679124*^9}, {
   3.794440741561021*^9, 3.7944407529959326`*^9}, {3.7944416101788206`*^9, 
   3.794441623091748*^9}, {3.794441672040619*^9, 3.7944417758268566`*^9}, {
   3.794442137051558*^9, 3.794442146003973*^9}, {3.794442361502695*^9, 
   3.79444239518848*^9}, {3.794453457585288*^9, 3.794453479257453*^9}, {
   3.794453608205964*^9, 3.7944536378703256`*^9}, {3.7944536920351048`*^9, 
   3.7944537225836773`*^9}, {3.7944538362883997`*^9, 3.794453862989058*^9}, {
   3.794454892932234*^9, 3.7944549580825305`*^9}, {3.794455178423169*^9, 
   3.7944553402693443`*^9}, {3.7944556296930947`*^9, 
   3.7944556633802495`*^9}, {3.794455738379207*^9, 3.794455747091197*^9}, {
   3.7944557788506117`*^9, 3.7944558127932634`*^9}, {3.794456032698623*^9, 
   3.7944560338538847`*^9}, {3.7944568895869493`*^9, 3.794456890241599*^9}, {
   3.794486377853352*^9, 3.7944863788808217`*^9}, {3.7944865019673615`*^9, 
   3.794486531116542*^9}, {3.7944866995368886`*^9, 3.7944867013328905`*^9}, {
   3.7944867782921295`*^9, 3.7944867783881493`*^9}, {3.7944868913053265`*^9, 
   3.7944869455424533`*^9}, {3.7944870070963454`*^9, 3.794487017172444*^9}, {
   3.7944912802458057`*^9, 3.7944912932224846`*^9}, {3.7944914657757635`*^9, 
   3.794491516238504*^9}, {3.7944915720995617`*^9, 3.7944915938110094`*^9}, {
   3.794529573326285*^9, 3.7945297797781143`*^9}, {3.7946028248126116`*^9, 
   3.794602824956644*^9}, {3.794614464068571*^9, 3.794614466014864*^9}, {
   3.794617688628434*^9, 3.794617695907029*^9}, {3.7946177466043406`*^9, 
   3.794617974775462*^9}, {3.7946181150902805`*^9, 3.794618162571054*^9}, {
   3.7946182155582323`*^9, 3.794618226840307*^9}, {3.794618260945631*^9, 
   3.7946182690562925`*^9}, {3.79466725486517*^9, 3.7946672937813396`*^9}, {
   3.794671990693123*^9, 3.794671991133561*^9}, {3.7946722243071985`*^9, 
   3.794672225115382*^9}, {3.794672256618478*^9, 3.7946722675164423`*^9}, {
   3.7946723250961933`*^9, 3.7946723276677823`*^9}, {3.7946732022922955`*^9, 
   3.7946732078655615`*^9}, {3.7946861840423374`*^9, 3.794686184082347*^9}, {
   3.7946876156381288`*^9, 3.794687628780711*^9}, {3.7946916885440454`*^9, 
   3.7946916887760973`*^9}, {3.7946926448311615`*^9, 3.79469269044518*^9}, {
   3.7946927390827312`*^9, 3.7946927431096363`*^9}, {3.7946937059893456`*^9, 
   3.7946937230237436`*^9}, {3.7946937619163194`*^9, 
   3.7946937706422896`*^9}, {3.794694551107609*^9, 3.7946945734466333`*^9}, {
   3.794757474363106*^9, 3.794757476680472*^9}, {3.7947575081418715`*^9, 
   3.7947575761072874`*^9}, {3.7947591154680986`*^9, 
   3.7947591266904917`*^9}, {3.794759917210444*^9, 3.794759950465684*^9}, {
   3.7947655559923697`*^9, 3.7947655700555425`*^9}, {3.794775867190753*^9, 
   3.7947759163637996`*^9}, {3.7947976326025057`*^9, 
   3.7947976355783076`*^9}, {3.7947990115471687`*^9, 
   3.7947990118771524`*^9}, {3.7947990801793427`*^9, 
   3.7947991162030973`*^9}, {3.7947992606671543`*^9, 
   3.7947994111055126`*^9}, {3.7947994528533535`*^9, 3.794799455253042*^9}, {
   3.794799499906472*^9, 3.79479959913795*^9}, {3.7948264436205034`*^9, 
   3.794826585605047*^9}, {3.7948266241804204`*^9, 3.7948266242234297`*^9}, {
   3.794835924335146*^9, 3.7948359549437733`*^9}, {3.7948361021597586`*^9, 
   3.794836162046491*^9}, {3.7948379612177644`*^9, 3.794837984788723*^9}, {
   3.794838031032382*^9, 3.7948380313884616`*^9}, {3.794850442120225*^9, 
   3.794850477877018*^9}, {3.7948505876622915`*^9, 3.794850612016514*^9}, {
   3.7948507731077075`*^9, 3.794850953895576*^9}, {3.7948547126531954`*^9, 
   3.7948548754870768`*^9}, {3.7948549212874117`*^9, 
   3.7948549214405203`*^9}, {3.794955679844499*^9, 3.7949557101345425`*^9}, {
   3.795091907249049*^9, 3.79509190738208*^9}, {3.795091956721109*^9, 
   3.795091977255028*^9}, {3.801264211592141*^9, 3.8012642292823505`*^9}, {
   3.801264337115242*^9, 3.801264339238716*^9}, {3.801268611655065*^9, 
   3.801268616837158*^9}, {3.801271129931566*^9, 3.8012712047630405`*^9}, {
   3.8012754867957754`*^9, 3.8012754913396645`*^9}, {3.801277217403945*^9, 
   3.8012772246591*^9}, {3.8013472095252814`*^9, 3.8013472134017982`*^9}, {
   3.8013564836776514`*^9, 3.801356495014716*^9}, {3.801826524753736*^9, 
   3.8018265272794857`*^9}, {3.8018449482138414`*^9, 
   3.8018449488054485`*^9}, {3.8018449947931657`*^9, 3.80184500886556*^9}, {
   3.801846532454213*^9, 3.80184666370492*^9}, {3.801850020632947*^9, 
   3.80185008790644*^9}, {3.8018507986055636`*^9, 3.801850841046332*^9}, {
   3.8018509572182293`*^9, 3.801850963990383*^9}, {3.8018627105724707`*^9, 
   3.8018627489466*^9}, {3.8018630497445116`*^9, 3.8018631021665745`*^9}, {
   3.8018633631561003`*^9, 3.8018633659897327`*^9}, {3.8018680433420076`*^9, 
   3.801868097348604*^9}, {3.801868132868515*^9, 3.80186816246635*^9}, {
   3.8018692404707727`*^9, 3.8018692529038377`*^9}, {3.8018697367686157`*^9, 
   3.80186978971426*^9}, {3.801870208238078*^9, 3.801870208329282*^9}, {
   3.8018703116168957`*^9, 3.8018703171090283`*^9}, {3.801870379586094*^9, 
   3.8018704064243746`*^9}, {3.801873165870328*^9, 3.8018731896516275`*^9}, {
   3.8018732226422205`*^9, 3.8018732446887646`*^9}, {3.801873306803525*^9, 
   3.801873322653352*^9}, {3.8018743559758654`*^9, 3.8018743845621805`*^9}, {
   3.801879527761982*^9, 3.801879540174575*^9}, {3.801879689071541*^9, 
   3.8018796896467066`*^9}, {3.801882290031912*^9, 3.801882291135378*^9}, {
   3.8018844627892456`*^9, 3.801884463844269*^9}, {3.8018845311715565`*^9, 
   3.801884542121196*^9}, {3.801885094340872*^9, 3.8018852029777007`*^9}, {
   3.801886748180935*^9, 3.8018867786069927`*^9}, {3.8018868498074713`*^9, 
   3.801886874707864*^9}, {3.8018875455564613`*^9, 3.8018875456410847`*^9}, {
   3.8018886234447436`*^9, 3.801888658278296*^9}, {3.8019118249308815`*^9, 
   3.8019118579138837`*^9}, {3.8019132989873466`*^9, 3.801913306565405*^9}, {
   3.8019136538105125`*^9, 3.8019136708474927`*^9}, {3.8019137473487635`*^9, 
   3.8019137518391223`*^9}, {3.801913901392661*^9, 3.8019139018518405`*^9}, {
   3.801915576751954*^9, 3.8019155836182966`*^9}, {3.8019158659190474`*^9, 
   3.801915865977061*^9}, {3.8019176712165775`*^9, 3.801917671407187*^9}, {
   3.8019190347731533`*^9, 3.8019191398190126`*^9}, {3.801919215231352*^9, 
   3.8019192736210833`*^9}, {3.8019196996029325`*^9, 
   3.8019197394377375`*^9}, {3.8019198317089796`*^9, 
   3.8019199094363146`*^9}, {3.801920187860766*^9, 3.801920335078472*^9}, {
   3.801920366695856*^9, 3.801920397443574*^9}, {3.8019206263992577`*^9, 
   3.801920684455345*^9}, {3.801920801011139*^9, 3.8019208435083637`*^9}, {
   3.8019209304070663`*^9, 3.801920981192596*^9}, {3.8019210978648834`*^9, 
   3.8019211091649113`*^9}, {3.8019213724767*^9, 3.8019213854674354`*^9}, {
   3.801921464644923*^9, 3.801921494930662*^9}, {3.8019218337275457`*^9, 
   3.80192184167323*^9}, {3.8019223608377495`*^9, 3.801922361454032*^9}, {
   3.801923075829405*^9, 3.801923096787888*^9}, {3.8019237139814157`*^9, 
   3.801923714447526*^9}, {3.8019240861754227`*^9, 3.8019241055302534`*^9}, {
   3.8019248915426216`*^9, 3.8019248978278913`*^9}, {3.801925292184738*^9, 
   3.8019253011038637`*^9}, {3.80192668857377*^9, 3.801926688978444*^9}, {
   3.801928158706963*^9, 3.8019281975411673`*^9}, {3.8019287228321943`*^9, 
   3.8019287324259934`*^9}, {3.801928913976655*^9, 3.801928943255873*^9}, {
   3.8019294134662976`*^9, 3.8019294248403587`*^9}, {3.801929485627909*^9, 
   3.801929489327449*^9}, {3.8019297902878513`*^9, 3.8019298608411803`*^9}, {
   3.8019299116529346`*^9, 3.801929911882105*^9}, {3.801930942606366*^9, 
   3.8019309501579776`*^9}, {3.801931432551155*^9, 3.801931441630267*^9}, {
   3.8019315243508773`*^9, 3.8019315244249344`*^9}, {3.8019319543734837`*^9, 
   3.801931977569975*^9}, {3.8019322809856367`*^9, 3.801932298626684*^9}, {
   3.8019323554153767`*^9, 3.8019323866821723`*^9}, {3.8019324246971073`*^9, 
   3.801932448966215*^9}, {3.801932534283118*^9, 3.801932562451703*^9}, {
   3.801932839499998*^9, 3.801932855219019*^9}, {3.8019351198504405`*^9, 
   3.8019351453014283`*^9}, {3.801935842742654*^9, 3.801935843290395*^9}, {
   3.8019362417483215`*^9, 3.8019362634885907`*^9}, {3.8019363126943665`*^9, 
   3.801936338732231*^9}, {3.801936630249535*^9, 3.8019367291985207`*^9}, {
   3.801937021565485*^9, 3.8019370271781917`*^9}, {3.801937364531278*^9, 
   3.801937368390484*^9}, {3.8019376362393203`*^9, 3.8019376637496834`*^9}, {
   3.8019383065593557`*^9, 3.8019383301921854`*^9}, {3.801939266031307*^9, 
   3.801939285079913*^9}, {3.8019394157012873`*^9, 3.801939468390563*^9}, {
   3.8019435168192854`*^9, 3.8019435269713173`*^9}, {3.801944250974132*^9, 
   3.8019442878920665`*^9}, {3.8019455280829277`*^9, 
   3.8019455446538677`*^9}, {3.8020172252160397`*^9, 
   3.8020172685863495`*^9}, {3.802025352710949*^9, 3.802025404046271*^9}, {
   3.8020255806367674`*^9, 3.8020255820874176`*^9}, {3.802028396150092*^9, 
   3.8020284054391403`*^9}, {3.802031607018956*^9, 3.8020316129618063`*^9}, {
   3.802033724908978*^9, 3.8020337416616087`*^9}, {3.8020338331127877`*^9, 
   3.8020338549836864`*^9}, {3.8020339234396467`*^9, 3.802033969118572*^9}, {
   3.8020340227333565`*^9, 3.8020340228957095`*^9}, {3.8020340529029665`*^9, 
   3.802034093795556*^9}, {3.802035748883009*^9, 3.802035750024272*^9}, {
   3.8020365898825765`*^9, 3.8020366075995464`*^9}, {3.8020377124634686`*^9, 
   3.8020377132083445`*^9}, {3.802038339890069*^9, 3.802038340384682*^9}, {
   3.802039203790666*^9, 3.802039207186472*^9}, {3.802039645032303*^9, 
   3.8020396486624994`*^9}, {3.802043898849061*^9, 3.8020439135977106`*^9}, {
   3.8020444498368835`*^9, 3.8020444761206636`*^9}, {3.8020445152969522`*^9, 
   3.8020445474986267`*^9}, {3.80204550439827*^9, 3.802045524292509*^9}, {
   3.802051727211583*^9, 3.8020517277927136`*^9}, {3.8020526939111147`*^9, 
   3.8020527023167715`*^9}, 3.802056336924117*^9, {3.8020593487851925`*^9, 
   3.802059349283702*^9}, {3.802059623032366*^9, 3.802059675132579*^9}, {
   3.802060005573262*^9, 3.8020600057064247`*^9}, {3.802090073109006*^9, 
   3.80209012716261*^9}, {3.802125032564892*^9, 3.802125039352764*^9}, {
   3.8021252886220427`*^9, 3.8021253014967356`*^9}, {3.8021311534461193`*^9, 
   3.802131176925544*^9}, {3.80213271111032*^9, 3.8021327233200655`*^9}, {
   3.80213704651925*^9, 3.8021370497176895`*^9}, {3.802172730167743*^9, 
   3.8021727518857164`*^9}, {3.8021870419251313`*^9, 
   3.8021870710444613`*^9}, {3.802188649673108*^9, 3.802188670491143*^9}, {
   3.8021923236072445`*^9, 3.802192330354847*^9}, {3.8021924460398345`*^9, 
   3.802192451772171*^9}, {3.8021924854528704`*^9, 3.8021925074673557`*^9}, {
   3.8021925827288017`*^9, 3.8021925828775835`*^9}, {3.8021927938333163`*^9, 
   3.802192928523094*^9}, {3.802193058338052*^9, 3.802193067236809*^9}, {
   3.802284300686278*^9, 3.8022843098492355`*^9}, {3.80228442144822*^9, 
   3.80228451614888*^9}, {3.802284547353369*^9, 3.802284611889417*^9}, {
   3.80228470895325*^9, 3.802284727719279*^9}, {3.8022854606013136`*^9, 
   3.8022855072631674`*^9}, {3.8022863124795194`*^9, 3.802286317301363*^9}, {
   3.802299530157411*^9, 3.8022995716277094`*^9}, {3.802301311263026*^9, 
   3.802301328301696*^9}, {3.8023021084399505`*^9, 3.8023021372453165`*^9}, {
   3.8023021970016317`*^9, 3.8023022034928637`*^9}, {3.8023041948615847`*^9, 
   3.8023042012494497`*^9}, {3.8023048130560827`*^9, 
   3.8023048374628477`*^9}, {3.8023078066403637`*^9, 
   3.8023078202470446`*^9}, {3.8023082201221232`*^9, 3.802308257117585*^9}, {
   3.8023083996502686`*^9, 3.8023084083431053`*^9}, {3.8023084685062017`*^9, 
   3.80230847744431*^9}, {3.802396758363549*^9, 3.802396758480236*^9}, {
   3.802399942549632*^9, 3.802399949772826*^9}, {3.80240130558005*^9, 
   3.8024013059147725`*^9}, {3.8024013950181885`*^9, 
   3.8024014034218955`*^9}, {3.8024016002481833`*^9, 3.802401600487198*^9}, {
   3.80244027949909*^9, 3.8024402831857224`*^9}, {3.8024752557102523`*^9, 
   3.802475272677368*^9}, {3.802604060265811*^9, 3.802604086095415*^9}, {
   3.8027483326073933`*^9, 3.8027483327943926`*^9}, {3.802798454868368*^9, 
   3.8027985278169622`*^9}, {3.8027985709111624`*^9, 
   3.8027985803945384`*^9}, {3.8029801956399345`*^9, 
   3.8029802067490005`*^9}, {3.80298181402057*^9, 3.802981825187064*^9}, {
   3.8029826304118185`*^9, 3.802982634652507*^9}, {3.8029831259435415`*^9, 
   3.8029831403931017`*^9}, {3.8029858778681602`*^9, 
   3.8029859182260327`*^9}, {3.8029868263242855`*^9, 3.802986833508671*^9}, {
   3.803008412005069*^9, 3.8030084490777874`*^9}, {3.8030088261384134`*^9, 
   3.803008878125034*^9}},
 CellLabel->"In[45]:=",ExpressionUUID->"9c21872c-cc44-4cd6-9620-54ca4e5713c1"],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.368417, 0.506779, 0.709798], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJxF1H0wE3AYB/BxnVwtRty4W7ste2NYXq/oWCph5864cUg1SrRCTUup3YY7
ui4uK1xKL2px5KboyFlqdRolzXWt7Y6h3M2m5VhYatUdz++5+973Pv99/3rI
2YXJxxwxGEziv/zvSH2xyW4PGMCsnbb3niX+F7KoYWRB9hNZwWOs0ueQKaN6
bJIOGfsmJvBuF/J4B+40Ox9ZWtb+U6LxB1/GC4YfXmGCia8ZzcIMP/BvSp5n
FMsX3CQM5LbP0sHvvyzUqPpoYGUAu0EspIJXwlwiw4Mo4Jc8vaFugw+Yk8NK
alSSwTlVZa/iqklon6qj6dsgETzKzeVI2QRwsCvhT8FXb7BXy8TeBzw8uLZV
rtRoPMB0nSA6pdQdLHl3SR4+5Qo2xpB0JyibwZlF4yGBDk7gxYLmpWkrBlxZ
r64uLVt+sW6BeKP15i0zOPaqbvn2tjEwLt2nJTdKE71umYqMXcmeA0/uICXg
actgq/OqQR3swF43lzH55EaEE1hkVJf0mzaBh85zsEwHHFj5ISSVyncHE/u9
BLRuD7BtxGKeiMSD+zrf5h2o9wYv98i4x6kEsPe1I6FhVURwCsZzvqKQBE4b
OiPmdJHBIxc/OudPbweX72lVTzEp4O9cu+VwARU8nJC6NbSRBp7h0U4dnaSD
i+OXVCY3X3BPVaxEtNsPHKFzEWZJmOBAD212Y4c/OKHStmh4ipxrI1TSepHv
GPhtnSpk3GPTwqAWeX6/vWLRMQDcdY72KDENeXQ2bpfsILI5SzCs5SNT9yks
OSeR63AROy9IkUtaE9XyNuTrhKJMswJZUVM7F/QM2Xj2s1v/AHJGDD9j7BOy
qLvc7KVHrmXIxYcMyEMupvvGWeQZ6ZYw1g9kRytrsNiKTMxLTn9uQ177X+C/
Vtkonw==
          "]]}, Annotation[#, "Charting`Private`Tag$6084#1"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.880722, 0.611041, 0.142051], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJxF1H0wE3AYB/BxnVwtRty4W7ste2NYXq/oWCph5864cUg1SrRCTUup3YY7
ui4uK1xKL2px5KboyFlqdRolzXWt7Y6h3M2m5VhYatUdz++5+973Pv99/3rI
2YXJxxwxGEziv/zvSH2xyW4PGMCsnbb3niX+F7KoYWRB9hNZwWOs0ueQKaN6
bJIOGfsmJvBuF/J4B+40Ox9ZWtb+U6LxB1/GC4YfXmGCia8ZzcIMP/BvSp5n
FMsX3CQM5LbP0sHvvyzUqPpoYGUAu0EspIJXwlwiw4Mo4Jc8vaFugw+Yk8NK
alSSwTlVZa/iqklon6qj6dsgETzKzeVI2QRwsCvhT8FXb7BXy8TeBzw8uLZV
rtRoPMB0nSA6pdQdLHl3SR4+5Qo2xpB0JyibwZlF4yGBDk7gxYLmpWkrBlxZ
r64uLVt+sW6BeKP15i0zOPaqbvn2tjEwLt2nJTdKE71umYqMXcmeA0/uICXg
actgq/OqQR3swF43lzH55EaEE1hkVJf0mzaBh85zsEwHHFj5ISSVyncHE/u9
BLRuD7BtxGKeiMSD+zrf5h2o9wYv98i4x6kEsPe1I6FhVURwCsZzvqKQBE4b
OiPmdJHBIxc/OudPbweX72lVTzEp4O9cu+VwARU8nJC6NbSRBp7h0U4dnaSD
i+OXVCY3X3BPVaxEtNsPHKFzEWZJmOBAD212Y4c/OKHStmh4ipxrI1TSepHv
GPhtnSpk3GPTwqAWeX6/vWLRMQDcdY72KDENeXQ2bpfsILI5SzCs5SNT9yks
OSeR63AROy9IkUtaE9XyNuTrhKJMswJZUVM7F/QM2Xj2s1v/AHJGDD9j7BOy
qLvc7KVHrmXIxYcMyEMupvvGWeQZ6ZYw1g9kRytrsNiKTMxLTn9uQ177X+C/
Vtkonw==
          "]]}, Annotation[#, "Charting`Private`Tag$6084#2"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.560181, 0.691569, 0.194885], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJxE13c8Fd7/OHB77733TipZZdwTRWkQJSklWUnqjRQhSkZJKUkpRFZEZJRw
X2bi2nvfa13XvFf2/N7f4/F7fLr/3Mfz8Tr3rHvG68g63rZypqOhoelgpKH5
f9/6Az4zu7t7geb/f3p/flo4sfnPvvHNf1+v/PO38ypbynP/rNA6wGHZ/88c
tcYayYX/PJzL8x+68c8hj3JWgtvV/+dI4ZuNac/2/M9SNSqp3nZq//O2gpug
0T7V/1nP6ELzHrV/9rpgGi6q+M9TkQrrS2L/3D6PH8xm+uf0Hxc/iw6r/M9n
Tp08uBz1zxHOh2dHw/+5Okg1rfXRP+vlMwvn3P9nOaGaDUenf14ZMahq1f/n
RG+NsznTyv+zbqRqLnnin1sTFdi1Cf9M2yBWW97zz04yzIdaav5ZrWlEZjHx
n6tH+wN13v/z5bWufv83//xcHhdLH/XPZL8fLIJ+/1ysFEPWtf5nC4Oo0wFn
/nnqbPgXOPHPooGBTifQPwd03Oi12/vPgiQn7SSVf87dufpqTP6fR1RtTnmI
/bNxsDEEMv/zwBtDySq6f/bJ1vNn2lH6nzkqD/aYr/9zWreG1oulfzaaVY3p
WPjnblrFeeGZf74tLHPy8uQ/s+wVz0wmKMHKbIUtmYyDZGMhxolBJfC7OPtL
YxoHh2x5HVV7lWC3VlT61hgO3B8zSxQ0KwFLks8EqQsHTX1/X1T/UgKJs2pe
Y6U4kGAXDEgpVoKqhXqmmO848DDQcQvJVwLXaNcEoxwccCTdR5gMJShq/Fzz
7iMOTjltk0tfKYGFqbSwZTAOPryRGnr3nNqf8fLs7fs4mK3D/LkfoQQJjy+j
7P9w8Ez10SedICWYgnc3mK7joGGeyarghhKE6AuUlR3DgZiMilGMkxKo9BdY
uhtR+3v2hNqdq0rQdP/shLAuDkofuQtZ2CmBjzDZr3YfDlgLn9FpnFcC0eJo
Lm8VHNhO5MxzWFLry/vjmqeAg0yh5v4ZcyV4lN146YIcDtbMFuoajilBaEaT
xa40Do778XzPQlSntphkSOIg/suBpAh9JXiS1KZrIY6DqQGrZ646ShCW0LFn
VQQH4UZvriuqK0H46x5+MwEc1NGLCXwSVoKoxyNjk2w4ENLWpwnmo/ohoSea
BQcuLpdnr3AqwfMHY406TDhg/JNYI0GvBC+8id/DaHFwNFr+bvy8Iry6vvBY
aa0RYrFHr90jUX2Vcq95uRHGyc6nbcYVIfbS35u+fxvhsXWmokC/IsRZr1rX
zTdCe+gf3r+dVFusm3nONoJc8fR2W4sivD25qS803QheRPbpbw2K8M5se18F
sRGqRPZ2v6hVhPcmu/IuE43AZ36myhMUgeIcGx853giOD27nnv6lCLefMo5k
UOsryHn5Xr1YEWZzfRXrlhqBdjg/jD1fETw6iDfHtxrBkrvDazpbEWZWbQvo
GKnrDy1d+ZOuCG4SDWsynDhAqTq64R8UwdE554k9dX4r9qL4IG9FsIrU/mlx
BQcjzo2+dz0VoZD/l8hxFxzQJNqc97ihCMKJ6D7yxIEx1y0+u6uKMFBwUvdA
EA6um66RLe0UwcCwLU41DAehQY9bzM4rQuJvmxXZaOr/Mf8uSvsktb3Ba0V8
iTggKineVDdVhBoXogB7Og5Yrn47IX9EEZQpHj70uThQfauvImagCJEPFjs2
i3Bg3lLHxKtLHS/j/YNL5Ti4yWw1waypCKdfbr+arcVBFGaoekddEX5aavuy
dOIg555byrKyIqxrVO7oL1D3U97f4Fk5RdDiPB1+m60J5olBV8ckqfM508ud
qtgEXDJsRv0iivD1j1N8N2oCy5cym7+5FEEuLCDL8F4TdN45fTd1QwE4d7Mb
2cebYPBAxg27dgUYJrPxZJY3w/YNzeNnmxRAtCUufqC3GaRSypWO1yvAua9y
stxLzeDA1zmmXaEAf24c1vRVa4EQ86tV6j8VgO54XWnWsRZIeTSdLF+oAIZK
ViZDDi0wtkhzhfeLApSM3jh39G0LMOyJMmBJUwAyLA/eK2gBhevC4rtJCqCZ
FOKc3dQCxxJS1pffK8ClfoaEkYUWcOnY2zv7RgGOBp8d/iPYCuHsP4vHXlLj
r2UWUg+2QqbJ0Tf9zxTAO31hN/BsK0x/t7OqD1GAZNxzWc3nrcA+O7EfG6AA
P0Yua3Jkt4K6wn/cxfcUoGVxj8lkfSt4xobjUm8pwI5Ig9N7hjZY9Ck0fXhR
AeycmTLpAtpgXZtLX+cAdT52gjQ63Nvhvnta+WN1BWDYOp6SHdYOa4kGqE1Z
AUjrfEKhKe2wwnzzmIeUAlQsp29r9bfDXQO6+p+iCpD9985/nCvtsHTnnTmz
oALEUg5PTPB2gHfa/qZzPAoQvMBwsWJvByz2/bZIYVeAm3PNuLgTHfAf19X2
BSYF0JLrpkv5rwPIxivnDOkUYPauGXEutQPmsxXsetfk4Yqk2rcnLJ0wExji
xDUpD9oeC0eK2jrBvUCEeIkgDwtYByVahi6Ynsy7kTUoDx/429lO63TBlMWw
59EOeeD6Vdgx/r4LxmX1/fxBHiis/s4CTt0wVLv0fPK9PORkMER6r/TA5Y1n
PFpx8uC0efcWVqUXBjXkX4fEyAO/BfEs+6Ve6H97Nl4iQh5CV/6IpWJ74SJu
Ssz9sTwYmuvvzlN6oZcm+GNJkDysfMwZO6zQBxe0hWUY/eUhlyJZH2bTB903
clOs7sqDUtb5A70v++B84jHF5Dvy0Jfp2/i8vw862wcz5m7KQ3TmW2cThX5o
12f/GuEoD8sZfe9yS/qhhdOtpMJCHgrT7emeDA9A5fkBfByDPEy36lgtGOBh
LS1+MWhXDvi0T7FMOuFBY+U8g+uGHBi/u1YxGIWH93GtyrpkOUi5FqXWMIgH
r55az54BOUCfFUfs7hJA4eK3LeHvchAVg1048moU7LJucdHkyoFuS3vbcNEo
vFxXk5nKlINRjsnvD/pGYetdmsmPRDnQi+C4XyQzBl3975/aPpWDoaCLNKp5
Y8Cxx/YDeiIHIWUeo7XtY2D8QDBXJVgO1Dce1jiujEGu+Mu2tbty8PhuesQH
w3EYv3l6DH9HDvZ+/+l+6No4iJWxLdfflIN+Mu5Ud+g4hF1+IvLumhzs8/jL
w4Mbh7IcY7WQy9R4FtPfnIVxWNza1b9xQQ6eEEW7TvBPgMrp8tNnreRgv+Le
kkmdCbjy0f/qodNyQCNhdmVdewJi53T/kz0uB/eumzbwH5iABsPlR6wmcrDw
5ZiOhvoE6A7fTu/To47n0FGu63IT4Kmx90flQTk4H2LiHyg5AZ+Dpv9kacgB
rt548q3IBPBIO8/5KcjBzwtHKnBcEzB11U5LjEcOYu8aeepsjoPUN2FTOnY5
YC837LdcGYdzNF0Xphmp42cwNL1JGYfKZIsHpZuy8N9rfalk4ji8JxytsiPK
wul8vWbWznFoPUDXaTwqCzVruocVWsaB6RF2Qm1IFvSRbrpRwzgYtAes8vXK
QkG4Dt/F2nHwkjvMutkuC8ot2kHeMA6ZXqtiY02ykCSkPf381zgMVxWqN9bL
gtAVLZvM4nEQ4Pcy+l4tC1FpB6uq8sfh5PV9lgkVsuC6V92fIWsc8HlrjP0/
ZeFP/geBJxHj4LNV+UukSBYUdDjzGN3GqfnLOeW4bFnoN5ofY1Ieh/nffjsv
42XBxvLbada0MYhRr8kJ85IFU28tWa5Po6Dk99yx9pYsZK5+/hUdPAqltTYi
DDdkgSFA0IbbYRRGr5Aeh1yVhdLQlac80qOgFcN1KeCULOjE/Vji+0iA+qEe
nl9mshAjpvridQAB7NU+1a0by8LfxHeqApcJsOjrHqBnRG1fnq0mVp8AYdUH
Ne8dkoX8DP8rguIEEOPZJhZpyQKn+szamw085F6u+7i0Txbcvl16LdSPB+Os
F9YH98hC2pPjiV8q8dCzbMvqpSQL6R3fDQpK8eBhLIf9JisLGbLSAz+/44Hm
xYzPggQ1fvupX2UOHlRUgvAe/FSzXytuTcSDv11qiSsNtbwJ1/6FR3jgyvC4
lb4pA+kxfs0rD/CQ+ldbfmJFBjJHxj12fPCAi6qPdpyj+sGvLE5XPEhi55zt
+6n+7iavfhIPBewl4h+6ZCCLtrPy4FE8mNoGt/W3UuMWGAd9Qzx4UvgMbX9T
PSP40XwfHuiMBhfjqmQg41CIvpUKHt4+TcvsKqfGw2f7LsriAeT1BKwLqfXJ
Vwvd4MfDuTu0DTF51Ph/GkV3OPBAKmt42PqFGse+s77PiIcg1lht7nQZOCfO
Nr67NQJ8NvYzpz/JgHUJXoZ7YQTSU5Q+RX2QASvrEnspwggcXliwaXwrA5xX
vs9bzI1As/5PDrbXMtC9ufVZkA0P1yMeVR2PloGU48lWjyXxUD9Q5J4UKQO1
+w1VXuynnmf7SHwroTIwJtK//d4EDxtdFtdTH8iALIk/q+AGHl7KSzLsuMlA
amT4LqEID2UVJce+m1D73+CZg84QQGW6L0Z/URpefgx+nR85Cn9fotGpKWkY
nT6ctS95FCp0MzTjRqRBS2+54mvxKJx74tNBxklDX/uN6ayxUQiW4RZMz5AG
GRbrIymGY9BrczSex14aCrwUF15RxiB1+8tUubU0MMAIAy/LOHh+5j1001wa
bDjfi72QGgeGxeG+Wl1p2MjgMn12chz2P/cTf8AnDSZDqx8ep41DZHVu4sRv
KTjy1113/r8JqPp0jbO+XAoqzj9JiI2ZgK2HAgFfvkuB3o+k3cP51PPIwP/i
7SQpOBDY+TucPAFWRccE1n2lQI7ZyFb29iSIZQxGcihJAZM4r985DyJYP3mx
Ni8uBaGBe4Y2nhEh6rqxaxuvFNDgjx35lE2EHenMo3HbkrD+2Z91fpoIY/He
u9JdkjClMf4u/MYU5Dxlu3swVBIi+XlWVhAJJtzKxwT9JaFidEa8xI4EUmZ3
rNZuSwIl//eRez4keEnfva/8kiTYWj58vppBgrsBn6ZND0qC4vN5hTWuaci1
O3dRVVUS7C43mP9QmYYpPeZ6dmlJeLon/c5942mwW/ZIa2WThMU/9mVrd6fB
yFPPwW5UAoCpyXp9cBrunZppNuiVgKWuTL+fK9PwTS3RULpZAtTSQpP8eGaA
xHI2h7ZGAi75ONQeUpsBOSK9+PhPCXhpYjCzbjIDl2qLI+vyJKCaT4S31H4G
YlNvrGWmScAq4a+O/70ZaA6RcH2WIAHeXwUZ3EJmgNmhpetWjARIfj0/MXV3
BpDRo6OW4RJQl/Om7ob7DPhLaH/XDJSA/3K6MqavzMBs7/uXq27UePZ591mz
GWg963IyBEnArS/nFxZEZ6BO5QAdr64EPAtQqNTmnoFfu5s/kvdKQMaZv68e
MMxA+tcYJRCToOaZL3WYF6YhgLWcdndJHC7q4wKlaqbBCx/+I3pGHO5zJlg4
/ZwG1xKr21Kj4vBm5Ibsl9xpsHKeGjRsEYe2x8y12u+mQaWS/0dAljgcxxlz
nPacBqn4YU+OZHFwSeQdfnV9GvhvZyl+iBOHR3fweb2208Bi6jO457k49V2c
FyJ1Zhp2JDCvfz0WB6xAkLWTyTQs/WU1P+kvDguzten+etNAauikGbgjDu/f
RyUVq05D5313zw17cfi7LBqjwTENDRbaik/PUetLxUe6b5MAq0QzKHpSHE6d
zXiUPk+CL51vThzSE4eMHC0fqTYShOyvVLjPJw521y0ucMWRwJc5aoCZVRwY
eYUszcNJ4DFs8+rtrhh8rxg8HnafBLZRs7vFs2LAJeZ+eIe6HjWmhAeW68QA
2xoqOS9NAgXsaMyTcjHwCDoppMZLAtG4r8cFC8VAVJ2P24WOBAxHTYq1PomB
V3gSzfD4FPQnecZ4PxADBaPSsaasKWj1PXSc3ksM2maCB1kTpqD2NMPuKzcx
CHxn1nUsagpKFVqK5K6KgaoZV3NI4BR823znUXBeDLqXOuvKPacgvd1J3viU
GKQ1UeI/uUxBQta+/jZjMXBY5X8Se2UKwi7UmlH2Ucub29r7n5mCAI2XO8FK
YhDj43/ilukU/Md4qYhHUgxOJn7QdjCagsvfyXL72cSgmoLnNNWYAk1uiR3P
cVHIi79Rwcs1BXGSBzsO9ouCiMdOwgXGKVjfY5651iIKDzGv/T5uEQGO37d6
9EsULCfLtFRniGDxqDP9zStRIB/kyUb1RCh4OfPALkIUbJnTIsKwRBBKojsr
HSQK0H/IBVdMhOFf+zczb4hCTMh1Wds0Ihg3mLXeuioK6+fWdj5+IEJa75U0
zfOi4KjyfHDsNRFYiHf9V81FoX5TtlSVev54LEdZlCFR0Gwpfnv7ERFa6T8r
hOiIwvuUk3eL/IhwkO/X+jF1UcBJnT39zJcIcTLtzWxyovCm9C+YPSDChgYp
tUVYFK7axGkxBBPhiiGNXyynKKgu6mVCKBEqTwqfuUgvCn+fD4gHRhIhwu3Y
2ui8CITVytCvvCLCrO/lpoxxEbC8Vu1b8JY6/ifeKR79IiC27TztSe2/UMqn
Uyu1IvBVK7uNSB1f2sB2MutHEehM1U7u+0Edz7SAb/MrEUjE9PLHlVPHs7bn
5OsIEXAb8A+3qiKClqDdsoSPCGzygWdjIxHi5f9rINwQgdpcx9HwViJsHohI
Sr8qAtHmjDZHu4hQdab4xP6TIiD3yNywfJgISvZNUstIBGYl5775jREh8ub4
3586IlD884WCzhS1/36b9UHqIvDwvGb84iwRLCP4Ek3kROAEpZM9j0KEwjhV
bxYREeB/fu/hzRUiCKeh402cIjCkIvZXeZMI/t8vSL6iFwGvG4mfuHinYLjS
c9FmXRj2RFvImahPgXHrk9/iC8IwWUCTes9sCtKGP3zAjwtDSk++fI7jFHhs
Npi6twqDiKyAouC7KdjUUk4IzBAG+sJeZSIjCUQyR0w+nxeGgb5re5+NUO+H
mobt+JPCsCAxYY02qPsRX1Ty/Igw0Dm4+S8LTMN+kag99/YKg8qk5++r5tOg
o3Vv4qa8MBioLs4JOE+DgaVjkoOoMJzxuCvQ8HAazCL0+M0ZheHe4oNrWkXT
cOqzfJPRphA8194NJ7VQzz/gCj9IEYLk+49yE6enwXZwHakQheD7L8Yua8YZ
uLI2viExJAQz78IYH+nNgJNAayFvB7V8gA+Fw3MG3Pf/8mT6IwTnrjoOvU2d
AV+3mNGFQiEoVzAq+so1CwGhAR/GvwiBF7P6J72js/Ao2dWmL1kIFKdFn1f7
zUJ0r2FDVZQQROctO/WNz0L68en8N05C4HD4qwDjrznoUjEJ0RcQAlZLCZ/L
rgvQf1RDfz+bEFRrsjkQPywA3kF0WYFGCO4Jrp30al+A2fgFN65ZQbgzR+Kt
kSLDYmG/HP2oIDQM48MjD5FhrbV2cLVHEBRae7bPnCPDzuy3uNkmQQiqbPYW
uE0GBtYPloRqQegtqCX1RZKBTTGcrfunIDVPUI21e0EGniNeNQ15gnAsX+zO
UhQZhOztg7Bp1PZ+sZ+KfkoGCb/jeoUJgpBQu6WsEkEGuTcHFzNjBKGuZY6+
6gkZVPKlcj6GCwKlb3jk0mMyaJGWpMO9BcF0Ht6+CCLDYUZ834MbgvDfWr63
agAZkGzj6ztXqfXRpVpU+5Hh1MVPzHYnqb8XesK86kOG639OL6jLC8IHzIl7
V26SwW1CL0tWVBB+nzhsveZGBk9ahetC3NTy1nv2vXIhg/+hje7dDQE47sZJ
rL1GhldZ6RXt7QLUPK/NZq8tGc5VZIHPHwHwP8cSRz5PHV9HTpUQCACdCKbr
uzUZeol5NT+KBeDZwF2Be2fJkLBVUGf3VQD4k3KsD1uQ4Qpvcf1WqgB8cBx7
tX2KDDJKPxsS3wuAkYrLa2NTMoweLsOhGAEYOfT72kMlMny2wDaPhgtAyEmV
/WVMZHBxqmoNDRKAo6V825SmBVD1q21XuisAz4c6392MXoDclMZud0cBkFC4
2HGFewFKlnuGFw0EQCP/mukx3nnwZxvAxx4UgHud8mMVbXNgID08qqMmAKWr
Ew/1Xs1Rz/PxSX9hAWpe4v5zD/8cNCSQ5+kW+cGu0UudT3AWor79JX+e4ofE
ea0/z7pm4EztyqLpCD/geVedGeNmoGN+a+Upjh9u2T5IWhOagcEjrLt8GfwQ
PPGYb0RkGhJtOGiLPvJDNYtpnm0fCa7d5Ka/EMtPfT+xnGp/R93PwXyM60/5
4YJFw9RJ6v07GSvInBDCD++8op7UipEgM0uE1fA+Pwy+OSOHGZgC9wpx9hFP
fkgwxqW7MJBAvUOKM8SZHyJ/nhsLNSLBAlGWW/4yP3juH5JOpd7vBVsKvLVW
1PoznC9XFpDAh1eF3/UEP5hIzcePzJJg7bCGcLYOP0hx7PJKXJsGuujD0gd5
+aGOKMPw5sQM1KQYynYx88P3K1no++MZCCtB8vd2+CCp60BgW/kMsBNMlX/N
8IF3tfEKp+YsCGhZ7zOp4wO5JKepMLE56D5uc2CijA+4heYU0s7NwTv7iwfD
v/PBStTda9XRcyAVflW3MZkP/viH9e/SzYNy/02M9QM+uG+Tibs/Mw+2Sr+u
qv3HB75xV5tD1BYgwostmNaVD/y6hVqf3liAHxUXk3sv88EDoea215kLMMWW
BXlWfPDQ5knHB+ICiFxYw4cd5wM0GfVdlIUMJ1LNaK8Y8UGl3PUXtuJk+KI/
eYRDjQ9qE7jNuo+QoT9c23FMmg/MeifkBKnrnb0z9FGpIB/UC5TtWFP3j4eH
fLUbLR80Pncrbosmg8ZHh6O1vbzQ8bBSsWmAut5JeU4fmnnhfNlbWo4FMrzQ
3g31ruGFrrVbQ+Z0FCA3fayV/cYLff+JvalXoUD+Tr9pUDgvjF73Zqi5S4FR
czXX84G8cD35BJ4ukgJ8b/3C1b15YWJQuuzIBwp47ROp77/CC9Pncd7YKgqk
PHCdyj/HC7depVjsdFGg43cxS6Q5L8w3399jSKIAgwCTqgPiBU92C+aALQpo
OZw/oavDCxQzxbFS7kVwyvl8g0udFwJGV9msZBfhzdrfyAlZXjDSdVRIFVuE
uqMmX8qEeYE2Cme4xL8IKy9fNbzm5IUavM6FY5yLoDREmHan54Un2p/uxDEt
go3qAXbjdR44/pT9KXGXAiWVLSfJ4zyAO4gvi1ykwBSntMfvfh54EWHe3T9D
ARE7z6jEVh44O1S4sGeCAvcXOZpOlfFAV1ikXHMPBQ7Lm5/Jes0D+D7NCL5y
CrjffucZHMkDKRofP10vpkDCr6noCw95wPkx86/CPApsWYe3MN7kAZL64JxN
CgXUk7sXBh14ICfElDkzgQKXZxV5Cm14wLP7m8x6LAXKQ2ssHY154O/DJ9YJ
4RSYb+X/75AeDxR1LnjMBlNASvJ6DI8GD9xTtQsz9KfAmRsF+UR5HjgUVJMU
7U2BoCLa9gpRHthq1/g54kGBXNqzi2+4eaBc+V37fhcKjJxO5rvFyAMJ3iWE
FhsKSGuNOEescEMyx4Fz1soUCMYFpuuMcsPntC913atksLM5dQ/Twg25R5wU
QxvJcBAvdvx4GTcUDUmGHkgiw+TfEtLFOG6oE4w5EkU9P0+LL6oHmnNDN/vC
/cfU/SHm7vK9uoALnHERThedZ6nvDa1QXBIX3GPYki3kmIWmQPrzXVFcEGZw
e4SrcAaCY5JXJl24ICvnvF0N3QzYid/+vWDNBWXj9cLSX6ZBK80wfg1xQauE
QZff2Wkg/ug/xCbOBYtRcpb7kkkAxlls/CxcQFP7hvOpGQne4+4NiC9zgsA2
S+P4/BT42JjmKIxygrJ2QAQmbgrO4AUC97ZwwuFbC8feG06BivvYaZ0yTvAJ
XHL+TTcFtEv5UpgsTnjWvCzQQX1P9wcGL5jFcUKK9Gr1cDs1P2O2AMvHnFB6
Z81rupQI0TGSMRfvcEJ75brsSgoRXMVnrznacwKJb7OVlpovi+2LpPfR5QSR
om0N0UtECLGpNn07wwH31+lmzDYnwQ4fI5zcywEvzBneW49Ngpa7w1RmLQdk
JDCeuNo4CVOB25GlSRzQY8iS4ZswCWfSdJuGrDlA5xHn9TT9SVDZx5Q4iTjA
ooOLL19+Euh+dnou7OUAFwWeyjL2SSjC/cdDy8IBcXV80p0DExBtgwisy+yQ
IyzQPFI9AW54rgK+UXaodRMMnMmeAImlbGuFMnZYZhMZoAuYgJVAf4W9WezA
dVn0KZfTBLQwn1jWjmMH+a9ih8ROTUBmjHCd0WN2MNgRn1LUmoBH4pNxZnfY
YbhIYDTVZAIupxW6WtqzQ8JEpL659wTo7Husd9GcHWwFd2PJKRNAMpbpd1dg
h04fkpkh3SQk6ocTh7jZgbiDFxzRnoRzWvNLlptssBrROxZ8YxJAsYxbt50N
RD7+fljTOgnxLLbH6IPZwLY2reRkEhFMm198SxhghUFBR2U7AxJs1a2Uc/1m
hdmki8ubt6n3Eda+MaSAFbZUz1Z/TCWBZP6eSddIVpA2QlcJbNT38evf4lq6
rHDNReqtW980pFzcDcO9YoHKH+Hzl91nwdbKJRYTxALfnULpb36YBa6TTZ8K
brDAZ54QEb/mWfAzSCiLRywQ6eZv/EZzDiyldRed5pnBUuRWXNPGHDCJJO72
9DHDkZob2QPq81DGw8h5spYZDtxxAdIV6v1D16F84AMzCNZfnWasnoedcU/7
7RPMMHTPysjo2QIUDnW539FmhmZ5C+tT5Qvg3m1wf0yGGaDlpJvdwgLItKQ+
seFghl7Dz0aFPGTo+s32+s8qE7S/ljvxS5UMT+G/ZIMxJjjUZmw2tI+az/3s
/ZrXzASByd+/0L1agJV8zC+5UiaISCTUmdfMQ86X9Po3aUzw6gP32KuVOXBM
5exmiWGCD+8Nd/tV5kDkg8/YgwAmyIi/KS5/aRZCnxvvXLNmgrLY39bfsdNw
6c7wQQ1aJuiPeljTEUEEm9A5Nto5RhD9MFzRXzEJlvFbhI5eRjiTbfiTsDQB
R0Hspd83Rij5s/l1wWEcNHguzNZcZYSPTPffsh8igIqCS/Xbk4zQLtgTw3cb
D/J6d9+76zKCnLmHyEHzERBxeH2ch5sRWMKnaC3kBoHPJ0V6bIMB7qSlr1yO
7gfOiPyVokkGmKt2mnHf6AW6by2fL1UwQAcNoTOsvQu2qocfaHxhACPppD+x
Rp2w2jNnRRvHAEcLd1a9xtthcWZLtTOEAfIDhd8n8LbD7C47bcYtBhAw229Q
Y9QGk/zivX4XGcCX58Tw7M1WwCur5Z06xgCF36tp9mi2QJfFcftFCQY416cU
2TTZCK3XL2jVsjDAXg/bJqkjDdBwz4U9fokelnYiee8k1ENFUuhPQxw9BCvM
veOzrAVqiiEY+YDqEZfzeV7lkELfOnvJlR64GRuGziqXwkfhkWoNa3qoUtvr
sjRQDK8w2/91qtEDne9f30PHvkGUNceJDCF6kND9ulPyNwvCXMVl/OnoQfOR
N/f25VR48OJQs3Q/HWw40tRIVESAb+rxtMVaOqBZc4lewJ1Ad0ouBNTm0wGL
SVNR6eoz5N7oYh3/kQ7s2O0OWasmIKeRu2o3I+kgvnTymdOez+jK31Bao7t0
4H9b5oF01Rd0kTm2l+caHXDx7dW+lfMNnRNPzRs7RQf11tVLXoYF6My+grBi
PTpoFn+nd3o2Bx03qbSPVKD250V2eKF9GjK50Kp1mYcOFFgDvuomJiC9h/Oj
tFO0wKPWlm9LcALZFjXM6be0cP9Ix5n7EXmw5Jm6OrZMA1+JL+5a7NaBq8/o
Fd4xGsg4zKnyqKoe+v1k64xaaSD5+dP+wicNUBmaFPvuCw0EHnyCEeNoBq2n
w5t1b2nAN4xu8VRLC2S+kLy+FEoDObOYdsyfVoh+n3DA4iq1Pkx5qFdUO9Am
98cHnKIBBdLmELtnB/ikidJ8OUQDn14d1k2z6ISpL7auPUo0wGfg99LoQBdc
+va2mUGABqInSkg9fN3QXNStrUlLA3dsj/7CneuGx757DmMHd7EWlUm5oQ5d
4B6SxSJWsovdr7b5ycCjEyyjlHt8Ynax7K9t3izd6wCdt2lpLTd3seOb+RE5
j9tBMkXeR810F1vnxBng9KINpkukefEbO9hdj4jazNMtEDomfPKk0w52+vt1
B7fs3/DjMAtWVncbG/RYS72l6RskHgt7HsC7jcWdPvrXfj0LnlgyXO6Z2cK2
x0pHS15IASsXmvWo5C2siGC+1GTsfZiNWdNcY93C5leWtTr5ZiMZ0lR68+AG
luFeUpBYfzViWnK5q1qygRUT31PprVuL5nbGTUJjNrCq2GL6ptg6VCpAwB8y
3cCeYW6OeGjxB51DfWJpeevY92+3Xo3XNKPwt3+iHzxaw0oM5lj39HQipeD7
3V8urmH1Lbv/O/SmC9W4KUv171/DXqnZfZFg1Y3oDj/5qjeyio3PscJdbepB
yXJay67Fq9hfMgHTlU97EYZ9zODt81Vsb2wai8LxPhQwiHDL+qtY0cC1o1NV
/Ui8doFfkX8Vq02RvW4ePIBKvyZeOje9gj3nfDIkx3AQ2cadTn1cuYKl588i
HH85iFaDtqYL4lewP/54TUn0DKI412zN0dsrWPdggwWy5BDSsrTz5zVbwcro
Mq3UOA0hL9kfrHeWl7ERn9/R31ocQrxsrmeTcMtYw0vX2Y4cGkbfFgXfNacu
Y8m8e3kFg4fRXLW3yl6rZazdQ5Aq5xxBbi77zKZzl7BCX/9+CBPAI2aL4WjR
sCXsnQjbZLIuHqXrPu8+br+Ebbhenmp3CY/GWGacM9iXsEFi4V80UvDoSnZ6
qJPbXywxXOxH9z4C2n59Hvca8xdrfP3hL2RNQB8CGASqhf5iPxiNV3zxJaDD
zt8vLc4tYldET1QJvCegvtOOqbK1i1jL5a+1QeUEdF+Hd8bywyL2Syvfnyk8
AQlJg2aw9yKWMecezophFBUx3/bPM1/EcrxRb7sgOIrOkSWrhmUXsbeZUhqD
Wwjoby+OlWudgu24J1yXFUlAryofnDVspWB1SFHQbkJALa/68AlBFOxm470i
hR94ZK79Ztd2gIzlkE8wijkygkiSR80iCsjYdsawg2NFwyiC6W90SSQZ+27q
joq22jCq7bGUEtIjY+VyTfn6+IZQ8JKIf7rqAvbwocUJmfF+VFHry8rYO4/d
ER/q97LtR5txnfHXw+axNTu/W2pwfcj30MsSmbE57Jmajz/dinrRzSDm5fcf
ZrFOFieefwvrRlGCirXXLWexa8yM0aYM3Sgnx/iNOsMs9ilA9GBIF2oycXBe
LpnBSvkFvPDa7UTz/YHaFTdnsJoKl78EPO9E3F4JjOHSM9gKWmaHqiudaD/r
zy6Ljmms+Ui+IMv+TuSlu3SXcHgae+09c8jr9g70upnX9Ms8CTt3r0CnL7UD
FTrvE/JOIWH9ztvPSt3tQCuv3YsZ2EjYVzzfL3wR7kB+ZPyKch8RC2FXNCov
t6P34Tt15Cgi9pQT6ziTRjsqlZJ4+xMRsT1HCt+d2m1DW6cu6J7MnMTOb1Gn
5VMbCs5quud5bwIr5V00MD/Zij6hGTPdPRPYLEuHGK0fraiqh0WEZmQcq63B
buYf2YoYmI7+iDEdx54mORQwqrcixY/XIuw2xrD5z4Vb4zZakKnWQ1v53DFs
/fFb+nIfW1C4Y+laoeAYdrVChO/QaDPKXO+pD/wziuXy9wysCW1G9S+X400D
R7FK2jVTFsrNiKTMf4P7wCj2V8nzNmXXZsSG3X+od5yALX2hIZtzrwXtsTnD
+imegB387qoncrQVnZq72XfjFNXSeg52uFYULZ7ht/kdjx1Y7c0bHmxDeQU1
J2pc8dj+61k9Ms7tqPXEqOhzcapb/HYd59qRPPmwjyHTCLY/Q8yCuN2B1LKU
lx9/HsT2Xrg8vyDVjTx6HN8Icgxie6rVhTQzulEuY6J2hvcAtmfftqHPvh50
wFHAt8G4H9vNnPh8DdOL9MRpV3kIPdiukhF12mv9yPR5/9pHyQ5sp/D1m3yp
wyjnjc3k3oh2rMlU+usMlRHEn9jeUb7Yhv08/XZTMHUE4XMbcod+t2KP5XPu
fXcMj8xKTD/c0mzF6utkSJ8PxKNLD5/wbx5qwdaWMDNXFuCRysvAeTHuZqzl
Ibd59Sk8Wkm+++fwOA57WkU2TFuUgGryb322+9mIPbD+ofaWMQHFVDk/9I9u
wAo2iDCm3ySgKx32du+v/8ESjRdnI14SkPr4ee1SvXqs1ns3X/V0AlpfOs3T
z/kbG0IZ3m35RUB1jKYz66O12GOTL5plqfs7VsioTvRHDfa48khlRz0BXVPW
+XToeTX2pJtG0ZMqAtLQ0wi46FiFPZMVmKlbRkBbx5Uu+OlWYknfHlyqfkFA
fy5Kab7jAGyKRDmG/xq1/S51Ma/qcuzliB3565oE9HtFojX/1i9srYOX6Osu
PHIU4QyjCP/Eth5oH8q7hUdbh7b1D1QVYwfoNFNwjHj09tIc5Y5HIZY0X7wQ
3DKMNAOHMr4JFWBFkq/sf2s/hOINE37hQ/KwZmeZ/vs6M4Cagx8oO3R/wS46
lDE/J/Uh+ppLr0fU0rGjwlvX3F17kB6zwe6Vh5+wHc36ZWYTnWggVq+7LOMd
lmtkpIB+qRWZiq9//WX5Env9TlY368taBGGTQtOTD7C2x5ReiNwpQOpYDQHt
aUfMVVtRg6e78bCS8vTNC7dwTChF8+vQQgXwNDL8Pc7+BuNXl+c829II6UeP
ips7JmImT/CKHcnsAKmRyk9pg6mYlVdt2Otq3TCVQNBWisjEMA7GOIfl9MJ3
W9qG9INfMYKKVuxZGgOAfcqT93rhG4YdW2u/kj4ICjdk1q2PfccMcEpJ4kWG
4anZfhOBhCJM9mXfwT9PR2BBET3vJJdgnkNKmvNhPJxjsOyJNS3FPFowJhxP
wUPp6FXZ8x/KMPelxiTU2QggXXn7puBiBUbLMtmDyZwAVS3Y/ulUwPiPX4ko
vk+A5Yrfuw3NlZjK+5KfXTKo+X1ui0L2ehWGfd+7HjcsAew+9px4plCDYZk/
MV5cS4DnUSOeNy1qMQxfN8gMOALAA+Lrk/51mCKF+MmgOgIsui/82JP2G+Pa
JGs/8I0AinarQ+yt9Rgx3+wO3QQC2J7YpZvd+IOhvVz8c/chAZ7qMavgFBsx
ks8O+/+2IUC5MvfpHEsc5lBpxeEXewlAFhL2inrQhCnJK2uKoiOAHJP0W4/0
Zoy6xuVJmSY8nF9WKjvV1oJJydncKXyLB56L8pkuRu0YBzfZMWu3EajvG+sK
xbdjzlz8xfmtbxhC7D7Tpz7qwPjOBWnWHh+GxUuKV0fqOjF1gkwBC0pD0HFF
WcjWogfj5iLAfXSsH6JGiCa+lB7Mi42egxet+uGoQ+Z/sa97MT+iE2w9K/ug
8JpqU2tPH4atRC4lPqkX4pz2hJ5wGMR8YT6gPXepG85MzOS70g1h2hOWLtI1
dgGTS87Ik89DmI19P4KED3fBPde9+lVTw5hTtpj6IyKdsH9qzg3/dATjPUs3
ZxNOzV+1Lpx8KITHvA+u4/VYaQf1ImG/jz/wmNnM05fiOtsgYnn1RQ4rASNo
yBucbdIG49q96b/sCBjDts7PUNAKyPdHeUM2AePiHP+nS7YVPhbHd/ZtETCi
f0sGu2RbYGPl/szU6VGMqzhGiYamGWx0L9KtJY5ijLhuXU22bIaCe4dEmcmj
GG6VyrOMhc3A9UN0v9CRMcz4EcGj7sIt4L62bqr4agxTeOmGTot/C9Tp9dtr
jY1hIu6Wq2gNt0DQz/fPrJ6MY+I6C5geWrXCQqk4MfrPBMaA5dLMcdE2OLW5
uf1RbBKjJPjEadO8DTL1BwW+3pzE8MjlDX0NaAOHsg/GjZxEzKg+fSsfvg1a
yyUTma2mMGG3swsH06nzsbNdKJQ6hbkT0LX3RW87RBgNNyouTWEuRu6mH2Hr
gCPYxHWTOBJGPdX6XfqtDkjcDeKxJpIwQt8C+S4mdcAG5qqyo940Zrcs4xl7
WwcUgPT54P5pTHv3ZuAdrU7goqXxeLFnBlM6prgq59IJ7kfwjxIDZjC/uk9J
k951Ql0IvPvaNINhebhY/kGhC+Sqkr+VSc1irJXjL1vmdUEQXcjvxtuzmKQW
w036w93Qb3xtuB9mMbP3xt6V1HTD62pZjnXHOcyTeo0eqf4e4DAyHv2YO4dh
SKB3j77UC6E/HH8Yb8xhQm71bm0P9oLv11SnqJh5TBDfI7mhkT6wfatYLg0L
GJ/LXbcSJgZgiiX1xCgjGbOo8YWG3XUQFE2/FxbuIWPu0D587T81CInQ8fTi
PTLGM13l58WZIXhRJKCbykXBXCc/YBChjABuUSHlrhYF43Uts9HkEB5Y9mtz
HrejYPIb+N5s+uEh+Mv5sdk0CkbjI1H55iYeyojOZ7CNFIwHkxVFzoAA6wq+
P2MoFEzO7bLS/gAC6DiGKTgJL2LmepVCX5UTwCsp7oWO4SJmj3HMafMdAuQO
pm+wXF/E7Nyd0/rANwozoiXOAxGLGO7/soQaKdTz6MLv1q+5ixgZD+e19VYC
OMX26Ad3LmIOusr2q1DPl+Q2YrrVxiLmmOPQrwsvCDDItcarKPMXY2P/7mOY
JwHOR4pM/bn5F+NvzXttXJ0Ar+pUrD/E/MU8O9NkzM9BgGb6QxWeJX8xH09E
KhjP4ME06GIsP/0SphJDO5WUhQc97/fI/v0Shlmd4nlJCQ8+375k74MljJjy
V8unjHjInysVoptcwpw72u7k9WcE1NwGZtMPLGM+7q+a2Ws6AuL2Eu/I9csY
ddaU9c+YYbB9r85QPb+MOb0rWk63MwSxPQa33wisYO4sxwQ7lA0Bh5X9scMO
K5hSQgizxKEhOP7i1jeOsBXMUM96g9/qIDzBBYqPZK9gaJr/i+4pGoQq1uiw
/LYVjHwN6ay2zyDsmCZSHq+uYPLqYmsCaQbhcGjuZRvJVUzFhOPJXcsBuFdZ
8VvFZBXzm+FAe/CnfiDrj3xsil7F4I2bRh4b9wFhv4y5sNIaRhjLuBkZ2g0u
L3f1BS3WMKwGZAabui6YXhjey39/DUP3s59LjqULFnM/8nI3rGEWC/LkSp91
AN1e8T4mz3VMR/pF8+lXrRAWtdHA8HYdE2cWkemR1QJss31ldLCOCT9iX8aw
txn4s+OTd3g3MPbazONaWo2goCrktlK0gRGWvHTwjXEtZEYsX1wa3sAwCO83
21tXDepTnScXmTcxKzyMl2pPVME3s0JD8v5NzARbv6d9E4BWxut98xc3Mbms
VxhnPlXADyZv2dlHmxi7gtQLDz1+gYGLFf909ibmjf4bzZCjJXBMkXd1YnsT
Q0ogt7KSvkJDKHlqTGkLo5rbEcb5Lg3OjLf0Eyy2MEYht/zkXr0Hm9ToiqFP
Wxhch0jZ4RJ/5CLLEdZ1bBtz3m0T43GyAE0Hz9zr8NzGiIywOrt7FyNPfMON
trfbmNHzIk/dEkqRb1Lk6SbSNqb4k5h7wd1KFC7JLFj7fAdz+oVuYW7tH8QR
OMlUXbyDKW+oqPva14heDtauwcgORpbJtC9nrgnFJ4QOlh3YxXgbqT3sYmtD
EpvXm0vtdjEKg7dYugTaUbKdCfx4vIvp8suP6ZTqQAqlcgXFObuYcOFlsU6V
TpQpSve5sGsXo1ek97lDswup+xHeFOzsYqatAvZ2GHQjWw+8SS8/DVo8KV47
J9CNxCRaDS6o0KAvAgFSptJdaLARq91tQIMchwbvJap2osQHeRrnz9IgkXTD
tpWDHchhT5JypzMNavVMVLMwakfjT4NE26NpkCGNw9CuVStKP+zJdzaVBpFi
Qm59K21GbtP27K0lNMhW3kTmWgsOzZ4w3G7C0yAt07onVWv1aJllC/9HkxZ1
PcNZPjlZibw+mvCQTGmRUebzE5fbKxD5wDMMyyVa9EnvEeUr5y80bSeWaPqY
FtVfvHnkbOh35LZwrcn5LS0624ltP/YwD008ztoKzaZF4cuZcqTlDITP0btU
3UGLej1XQq/WxaArR4KfjhJpkcPjY0oHPjqjga7fP2m3aFG8WIOvUkUE2Lpz
k2R46NBiw6sE3H8J0L1jI4IU6FDeUfy2sXEanHudaHZVjw7Ny69+TWrMgTbl
Sd+gU3TImZZ5X5BaAZwp25v+0YEOaYfoXjfAfgOc5d2uMh86lEafqnJm9guY
T5QxDEbQIRJDz1N5u8/w249Ba/MDHdriEOqnvEiAqpTXrw7V0iFQvabh63MO
fizlWxQL0iMuSlktU1oWSo2f/5NXTI9sDzHPxT2vRjJ7ddabG+jRkp5/rKxm
LUqsDFSZH6ZHn/Xm9HO661D8NEe4OjMDGtDteFop3YCEHp4rPinOgJ7qmmqe
qmlEr/k/TLjvY0CKuj/7ut2aULTBnqNZFxiQXoKFwIptK+Jo8/Kuv8mAeMRZ
3fT3tKFI59IU4kMGNPy+6lfwdhti2qBtZ4plQNliAdx1Le0oNPoErVImAwqw
KPsPL9iJ6ORj9h8rY0AyiZX5pZld6GFJ71WnVgaUO1tHidXvQf5414qUNQZE
iGz77wT1fbbmkzdXycGInvR258svDaC7rKsSBBlGxK88SNkOG0J3NMMeSJ9g
RPKOf/Uu2uCRS+hnvQ/vGJHN7vYPzoAxdFaFUPDmMBPKP8T/0/c3Cdk0rp21
O82E/ErjPtk8nkZ2njwUKQcmpKsv9lQHM4OcCjH7Mp8woXoD2UsrxbPI7cKF
Zo94JhRe/tmk22sO3drwvHUgmwmdNlJRL9aYp67PMI6VciYkgs0WiJueR74o
Mbu0lQl1YfZt301fQP5jReYPx5iQuUPfotoBMnoY1kQyWWFCVupPu97sIaNQ
1YkIFlZmtLQbF14qQEYRuC3lJnFmtDfDzLticwFF3Rb4HaPBjHzOXOthbl5A
MXzqLjZHmFHxsr/+2eQF9N720ucRF2ZkcjSXfvzoAkrc9DL57MeMImd+u+4V
XkCpiU9H3aKYEe4VodGXNI9yxn/KLOYzIxuCQCxr9DzC3hH5yLDNjNyCHijs
65lD4xFdb1Rfs6DoH8eP+d+YRVN75rTm01hQ51XHrGr9WTTXzNBZ8IMFiTAH
cHJyzaIVAS0+g2EW9Pl8XldSwQxi+fQq+owqKypZFHSp3ZxGe39ahHtjWREf
l7yOyGMS0rzsqqTXzoq4V+YmyNYkpLMbVLs1TvXwjzd/FEgIHculD2NjQ+y5
Z5b966aQVRtHcPx5NsRyxr9wiHUK2fgoSNu7sSE2nWNOxf1EdEnYoEL2ATUu
xSPwIpuInOxvbn1JZkNM82ne6DQR+ZL+3C+fYUMM0W2aqS8nUaLkiQFhGnaU
Ly7mGWw/iX6frTf0EmBH9lmOWfZqk4j8xCwZp8KOWHSzxw+vTiCR0jo6ZUN2
VFzzV1qkZgKh+WPOIWfZkZOVwaXllxPohlzt7wFndsSDD41rt59Ar2yOqun4
s6Nft5ra8tQmUOnT6qiX0ezosvDnuf+MJ9BohfHCdAo7GtBcxStdnkDsfyvP
HithR1fOmHcO3J1Aly+B0MYIO7r2hFx6LGsC7ewpT75+gAO5/jWIEOWZREoO
BvQVxzjQFNfLB82qk8gi9peziB0HuqE25vnYZBIlb/1Ua3rEgTyuRZ6f851E
xk3FhTqdHMirpUuucmgSudPqCMdMcaDFaRUh39VJ9Fq7yG9miwP5MAWw7uEl
ovGP342SFTnRfUM5cuxRInri+a2e9R4nepDtWe5Gnd/clH17nJ5xop26ym+S
tUTU3Z37vCKJEwWNCnxuHyaiHTYNskghJ9rddn0bvkZESpivVt71nGiixrJD
RWoKnfFWL2oa5EQc7zROmwdOoXsZ2cIqFE50+BbH75uDU6ie+8vgoCgXei5Y
X5r7noTIJqoYXQ0uVExK02pdJyHR+5mfYoy50Gj541yK7TRyx6e7mN7kQtou
mBQt4RnEkZ9KzivjQkPFxZGlr2eR1oScNVsbF2J6Fks/uDiL7EVTipwmuJDm
Va/A7bNzKDck2V+UixuFMWv8h3jmkaXVR4bHV7mRhm2abV3UArofJuE65M2N
bNQftxNnFtCn0oQ/uhHc6Ejw/vRQRTJalHsfPfuNG908xyWTHkpGYhdEKaa1
3Oih0fagTiEZmTyLt/7Ux43eqMy++z1GRh5Y4eLNOW6UzTdgY8tPQQ1+Jzn2
0PIg7NYffpIxBcmJey9/YeNByqLOqUuOFORf9n5YTYAHkVqvOWndoKB2+6rf
XyR5UE7EFUWf2xSkukv6pqbMg26hS5Pf71JQcDLv+y/7edC+tQsZfx9QUO+R
Q4/VDvMgct45t4OPKCgiNOK82mke5CN9ZrogmoLwit+MvtjwIN0e8+zFWArS
+92jrObAg9ajzTw0EyhoilV5Q9WbBwXtoPn8DAoye1GbqPqOB4UNHVz8BhSU
uH8uPCuFB5m92f+dXEdBK20C/6nm8CDW03t99jdRULqAk4kqlgdFlymt5PVR
0E7RM/Wseh5k6SNfsjBCQTYXvguqtvMgPnWZ+/smKYjpPd2UygQPiksQ3chd
pKAr+mptmfM8yNZa6Nf8GgUVD54tVVnjQeLs/AEauxTEFeSXmknLi4aquA09
GReRs/SnKBV2XpToz7HzlX0RlUH93UwBXuSgyYqd411EAo7kKypSvEhumjF4
r8gi8qAXOZ6pzIvGPtEduSW1iGo+Yw6oHOBFikFimpaii0jC1FUs8zAvuuDy
0/Us8yLyIUbTqxzlRU/P2H60WqYgebXhLuULvIgsFcdyvo2CHuUF+CsH8qIy
LO/IxUgKOlL6s1rpDy86S8sadX2djOIvEXLS23lRKCmj0olIRuQtljilQV70
o8101bmLjJIxtm5KC7xIIiXU0S2fjGhqVziUBPnQpDHNoVs3yMjWRWolTYoP
iexJ8vS8QEbfmE1HFFX4kDm/0efbx8iIJetWfdoBPhSwOdh35yAZXTV/k6+o
z4fyxx5we8mSUfFM2fu0o9T6GsWOeXOTEdfz8ceKZ/jQPjfv5amNBeSiwXEr
7QIfsmGKP9BTv4AqWg7aKF7jQ7qaZY8mny8gn+vhdnfd+dCMZNqm1LkFpLra
f6XWmw9lsEb7XBBbQLFSj1ydw6j9I1x1qU+fR+63Wu+nZfOhjRf05l+wc0iG
ViFwuZAP5T2YrRp9Moe6Yu+FHKvgQzddu/TFT80hVCb1dKKVDw0YZeyN6p1F
Quy3EhRX+FD5nDmvJ2UG4RIh6e4uH/LqO/g0vWQGhWgKfK5l4UcatRL0I4Ez
aPbirxxncX6U9GF+yYJtBqXMceUXKfCjcxE9tyJap9GFEMciRg1+xOMDkxA3
jaoyWSvSED8KOfmq94D8NLpvYF+1fIIf6eg+OOtOIqG9rd/qjlnzI4qcU0NK
HgmNXmdofHOZH33lOm0ycJeE4lcvtEw48yPHDe0yfgMSOvMsu0P7Nj8SmZTS
PkVHQvTSuz1P7vOjWAGlCfWhKfSjwGqwK4QffTiPDzknQUKe/9fEnQZC/UV9
AFcpkiTGvs2uRZZCKPUNSaFUipAiW4UWe1GJEIooRKkkkhYh0aLNVsaSfc1g
MIxlhJJKz7x4nr/n5efF3HPu+d1zfnfe/IwzmLQocWQlJCsG2g6C2vaL5XNd
HAXNlq/v8+Zfq4f5YOltcbyRXmZd2ToIoxuT48454qhODrkmy5t/siyDBRmN
4viR6s5fZMlBWXAXjU4kwCBz09CU+QgCJdat8l1JgAV7Okzhygg0ssJUy9YS
YLkyj7KVMYJbtWvWu2wlwC1b+eD1HaPY4xS8oWAnAb7D3b9eXx6FwHTD5kXW
BASuuZXQW8G77ygFmmQcIyD2mWjt2m1jWJFXY/bDi4Bb3C/uNmFj6DSm7DYO
IiBN49Lii6VjiGvz3ZcQxotXpFRvOZ8LE88vB/pjCGhvDbzH5p3H2XmK9to3
Cbz6tZ4I3MJF3o1TjmFpBJAJ5vR0My7cVpa6NGUTYNIZF9u5kwvFt9LH6S8I
8MxomZHczYWLDHH3lncEZNzf1nZhLxdPfJTX230moPNuQdHQPi42qmoviOsk
wDzlun+xDRc2d/am/F0oAZVYix/OTlyk/bYJlhaVwJEr7xprXLgYtHJ0XScr
gZRI1Re6vH7zFz217qiqBJZcEvZe6snFjeCrlQ37JTAYUDGe78dFR+eN52OH
JUDyW/9V8QwXVL3biULHJXDAOyMnIpCL3O/ZR3BeAhWeoSfsgrmY2Zm33fay
BPjcJ3aWhXCxJfuVmm+8BHSOOqqqh3FRe+Tz7+xMCWQcwciCKC6k39d2lz2X
QNfhZwyPK7x+lm8p734tAUl7xcfNMVxk+Hc9+VMqgdvipqYbee+f0Yb+eKla
CXw+26tmzstfW2M0YG2bBCZ7zxLsefULujJ1yJwlAYoZ4Zcnbz4IG/OrhExL
4JdOZCiJd9+0TFsiljpfEhaLX5sb8p5fyqzYdKGwJB63ciSdY8aw6iWpZJQo
icNnzLMeUsZgsmKznc12SZS9Fq1WMxtFbKixgc9eSShEb0ncLTmKZqb5itiD
kgiwO33YizkCl5t2k6WnJLHmb/33Au8RhAididZIlkTSpkSpTbeH8cUl+LRZ
uiQmRSqYh1yHIfYpwtr1qSTMuqazgjWGcfdsIvX2R0ksvGCzsbSUA3ZzqlAh
QxK2FlELB2I4UNPM4NY1SeIl8U21oA0HvrFPmkaYkiCMDyeuonJQPJz/RpAj
iWmOTwk/7366aPubNMqUJEx8u5fTmUMwf/ApYtM/SaTwmR/aVjqEjoN1lt7i
UtgsQZ25HDOEmeqx4a4NUri8qU5tGXUIl9Z2RZpvlULrY8Etf4WGIJpQveL1
Timsktu8hzM+CPrBJ04JDlJg/Mz2Ln83iD2cYx1m4VIQe36p6Dyv3zt32px5
FSuFI0pvvngYDMItd7v0imQp5F353m67chDnAlZYLngihf3HDs2u/8nGI4H+
yqI6KSRTdA25cWxoHW88qtwhheFrJyy/BbDxvrpE4EafFDbyZTgzDrNhujb/
wfwxKVzx7PB9tY2Nphv3DU9OS6GrQyzioSobDr/iujvnSUPddPvNBAk2hu0u
njddIo2LRecfhf4ZgN/7UwpFBGlMNLDfTvwcwHyqw2u6ojT+Ptp5KmFiAFfC
LQ5cV5bGouB8qu7YAKQ5m3/O05DGUivZlvahAajmKq7rNJSG4oL+TaSeAdig
41H8UWnANkj4ce0A+u5XmszzkoaJRu+7nYwBnBR43e8ZKI3dAtu9xssHEFZ9
k7wjRhqOeYQ27XcDyLWzSuZ7IY1LSx5nfngyAP33xjqexdKI7l5u65Q1gAqK
dlN7uTRuvPQTEXgwgG9DBLGXbdK8fjL0Mbs1AKGAussefDKoeN1m0Bw5gBsd
H5TbF8ug9hp+BIQNgIjnpSbiMmh1zciSvziA7Pt3jxTIy4CpL3zwXdAAtAVi
51HpMhgUPy3qGDCAD8fO37mmJoPvg82f+H0GYFbtqf9PRwZxmgbvG88OoFnD
vt3dQAbHdMyuMEMG4HjDPKDNVAZGG/fbcKIGMDK9UcpknwwUcFj5R/wA/O1U
Xrywl8FPw2OTfLx8r1KWjMeekkGW6bmrko8H8Pp5hdu2RBmI2GZ83MKrV0vn
cfb5ezLQe6Gf+qthAJOLl7kVZsvAeVljwPNOXv0d97mufC+DN58WaJB4zydN
rMd5yaAMPFUc7vCJsfFu06U+wwkZJIdNnymUZaPj2ArnwL8yKGXG7D9JYUPy
k6fTyHJZKN4oXsrUZGPd2HIWTV4WO8b2DSbqs2Eh9+KIPV0WfttHSnYZsxHp
9duxRk8W1X/lAt9ZsZF5J7VHYKssflvlWfnxzmNJ5RZH7JIFPXfHOrWjbHT/
ZHX7H5BF5n2tvAfebMxSIhyeH5FFiXIWNy2UDTmL1d2DHrL4li2vevc6GzqB
1YfJ/rL4oxZ7/HY6G6cbCIfjo2WhpePXn1jCho2W6kGrb7KQOiBeF2w0CD+H
rx0xbFms6wxbdt5yENeveNtVfJeFhcOMWaDTIKr7XtnqCsjh8tHuct/QQRgk
bLORU5fDTMCzt8dLBnHo41DrXj05SM6Sf7s1DCJw9MqBaCM5rLuYoOPCGkSB
cYP1X2s5HI8MyjvMP4T6034tWkfkEC4yzj1IGAI3Vdba00MO9+OcVG1580K4
8m1zhp8ciiVbjltrDvHuT4etuoLlME0KFtu2ZQjGFP5mqWg5NAl2MVhmQ3Dc
lbnfIkEOeWMbwy/y3ve3Mkf2fXgkhxNvp3+/PcGbN9ZnLO/Xy6HP+mSfZuoQ
pgvjd7vS5NEQc6iQzJtv1sOKl7+rykM/XfXi74UcFCk9eh+kI4+HhX9NG0Q5
OBP2Xu26qTwCu1O6LtE5aCsydVXaJ4/+qWNZ9hocbBhpSn1kL49dQnpe6zdy
8GfvyNIPp+RBXdeyaHA3B3bh/lvNzsojeltm7Qc7Dt6+WhDUHCqPn7a+ycmu
HCiOXs13vCoPx5NbnbxOc3COJDs8kiiPylCCqhnvvvbN8gEl4J48+Cb3vCJf
5GBzhLotf7Y8RGaOUJqPcnD39eu4mHx5yPL5REfx4vGNGX+RLZaH6qKwqc26
HBwm183LKOftTzjRfpLIwYd9B3U1vsrDTOxh+UNBDi6+8Xq4jSWP44pfbi5v
GQJ/5N0wbyEFZGmLbOr1HYIgbcager8CaK4JYWH5vPlrfenMgcMKWO+ROaZ3
axCfo0RzWUcVYORVaD0WMojI73TS70AFOJ1vW2W9dxDC7/f+U76vgPuJijUr
v7MhavvkdfCYAqxi1RT/SLIhc+2wplaEIsyEQzw0dfvRUiAVHntNEW7+3xgj
cv1I6Khu5SQrIoKlq5I52wexFfrn054oovQNd0impA9C72S+LKtXhIGHvdu/
nX34M1x3aEhBCc7pQZYjl1l4JRaZu1VZCf7rt7jZeLHgr7Nl4T11JUR+4Q8s
t2NhKuTZIytDJTzmRt2/p8bCqGz0VImbEiY2Jn/f19gLpsnWqDt5SrjQWBBb
TOxFquefzpk3SohzO5O+WqgXdtfz1PeXKSH9t35h0kQPWrtITcKtSigjlXad
LutBne8s8cysEpZ61qsqu/fAsc7W9hGViNm0BRdEqrqxxfklv/VKIhjbNi0c
L+oGcVrs6UJVIpKH/SPrM7rxTeELn8N6IrS1R68nXeiGzVGdDKkdRHh8bski
a3ZjL5/EeOgJIlrGntRrpzAR6LdjSZI3EbH3S7vOhDGRPnqelh1AhLlVJ6f4
FBOTHUMHvoYQUVYszL/NhAl5S6I3K4KIkNNU0ah1TBhV7rv68woRm+kb5WsU
mXA3iHooFE/E79a9K8SFmLhe9P6jQhIRJ+Wa3Z4OdeGN+o8O9dtEMKVetv3r
6QIrc/VPwzQiBBf99GOOd0FYyWG5VSYR0mdSIg/zMaGZkLD62GMiVp/a0For
wITdUsbWoOe8eG7tyluWMREaOu9wbAERew+d9X0uyUTDaffrBe+JOGP+WjyO
xsSfwXtPP5cScdXI1nH+GiaoDs0VHV+ISNvwO+e0JhPeuwz+zm8k4vPKDTv3
GjJBUJF22TXIi7/odQTBgYmN980vOI4SofLXpjnUjQkn2ZBknwle/MkZ2tQJ
JvIFR6tT/hLh2qP3qfEcE/v6Pq1nLyfhZfGrvwm3mDhn92v3b0kSql7YmAmk
M5FRr+ouIk8C8/FMsl82E9U7nMJIJBImeX9OBnKZ+PHh5l1NOgmCKXo61q+Y
UNStebVtNQlycW1hFR+YMM7hb7RRJ+Hs63zlrx1MeCrrjXlokfC9USTn/D8m
ElJPLA7WI+Eo101HldyNYokHlOubefGEPr3vMOpGf3SbfqYRCdY0he1Rrt3Q
Prv1dNVOErbafD3AftyNJiu5dpkjJJBLdM+mre/BbLXFlIobCUnf4hfstukB
3ThsGTxIEP01EvUvsAe+muOGLn4k/F2TdsvuYw8kl5c/zo0ioS1RqFjSvBeb
Iv6UlcaSsDvXybj0RC9c/ml0t9wgoYJRXO0V14uCkVsS/+6QUDDPq6u2uRff
nOrUxB6QoCJf5Xp+pheLOgR20B6RkK6tzF2jwILVl9PnTPNJuHm8nS/KkQXi
1luTAkUkNHdUugxeZGHoXemxkre8eu18w9h2n4V8vTHmhY8kzNbs0b2Vw8K5
F9JW+uUkeGb7dq58wYKJukHVr0oSvoUlB78sYkEs+7hhQS0J5o7FtK3FLGTc
LVZTayNhjcwiz8PlLN45Zj/gfCMhdXKl2GglC3oJy+Uf9pKwtNa84GwtC1VR
ToKkERKGw27MJrSyMOkr1C3GR0bWJDWmisNC8fg6qxp+MmRqTdbZclmIcD9Y
FbWYjMvZ7s3sSRbkHXOK+MXJcHPMV+KfZcHQfH/cDxoZyrWGOXqifRCuOC+Y
t4qMhGxXywpCH5oMss6dUCNjUXjU9D6ZPt696M8xtg4Z/fr1OEXug2YerTtd
n4z9Mj9Zs/Q+zK7ZZeVgQEbppOzl6NV9qHjoX6WwjQyt2k1rZNX7EEdJM2wz
JeNBtuPXTM0+2KVWFiVYkCERHuajpdsHusyU2t59ZJjIZlV3kfrAjVfMWGZD
hvuwte9Hvj68EjGRZ9iTca1YUPFBFwuhl0/FRRwhoyC2sDScV8+d/CmCW93I
aHd08zh2m4XuXyOTb0+RsUqg4o2aLQuPvKWOn/Elw6rVz0lMjwXvMXRrnyXj
bLay8JQ0C4J98VXPQnn57wq3fdXUC/Va3Yx7iWRYJxnPGO3pRXXgV8GjvWTk
vDL0OYFudLc3jEWyyfhkcDzir3Q3JvVamh6PkPH0itmNQi4TMr++pXN/knGF
suyZ6l1ef/sMw1+YgoupK0q1Zrvw20MgIFyLAhNz/fj+gE6IVAkdztKjQL73
32O7lZ0gqohsq9xMAX/Ax7K6lg4YDxEklu2gIDd920zx+g7EuZCfJ9hTIPxn
t0PiVBselNKSCo9Q0HONcHbpgzYUUleeb3OjoEi5+XqIZRu+hKi4/PGgIPXt
zae/FrShs0fNXNGLAp+9dhUn8lrB3bJOE/4U7Ft9Kcz/aCvm39OWcwyioKTo
QlbI2lasOKQ/+CCcgpgmv7GbJS3YUIza8mgKJp28xB5cacFOBaOXg9coODDh
qZWzvwXe7TsurblFAUXU5WwZuxnvrKxJL3Io+LjdUmBCuBl1BbaCzS8oWNGy
a9VsYxNYEofGpl9RkOBiar74ThOEGpzfbiylwCHEIF5Jown7LbwOlLRSMPVG
XclsXyOOPvNF/zdePFMVAyvFRgSKnFEWZFFQ3Krs7DjQgHuMC1OmoxTE/1DM
9g9oQP7q0A6PCQqmQ2WrQwwaUBYZ/ilmmgI3ccnxq0saMGxy9Vr9fCrWqi9d
/+B2Pf5lXvOfEqAisVjQJselHmICNw5JLaWiOWjcYAL1oLkkGeuKUXGxrc86
jFwPndKUNbZSVKivb/OU4a+HKfUOIUieivb46tDHfXWwD0n7nUqiIoL7MXlz
eR1CtmR97llNRX9WdplzZB1KbhoUD6tTYdOsXhfmXodF39tzf2hRUcFf0Jm5
sw6RaSK3hEBFxqH3k4NidYhd4O251pKKYxKNZM+Ur6izXXpk4wEqWg1sVGOC
voKQn2FlbE+Fycku3ZxDX5Ho1AobNypWVQ5afKd8RWrpJvGLgVQMXpg95/u4
FkyFFoHoYCqsn4ZEJl6tBcn31J8bYVQw2gUTCk/WwrFaaPxOFBWbF1+917qn
Fun09L6sWCqeaos/mdGsRf85/ba8G1QoOCUVyknVYkVzU/XbZCokKBsiAr7X
4JjayU/ld6gwKrWwaqquQXb44sKv6VQMZGywSc2rwXBX2uP2LCqiI+gHXZJq
oKqz8V7fUypUji0/rBpUg1y2Z+SvQiqOrhlwLjauQeH2BjtSFRVjeXGnvrVW
Y+aex+7VdVRcuRHknVFcjY0zi4y1mqlY4+fm53m/Gu8e6apt76bipN6moFn3
apQtuTPv1BQV3PeDEYp81RB00vlxZoaKmLSGqH5WFba/+ToU+o8KzdB3V59+
rgLDnb8haTENvttuXN8UX4X6qqMZ7xVomK7CXXt6FQj0BSlfyDSkPFt9n76k
CvvP3YppUKZB55pkxugYA0lNWqHfVGhoOT0vq6CBgTbVGn+2Bg3rtLx2KZYx
IB/u5vFdmwaw396jFjNwsGue458NNJinCE6uKmCAGaNpKmpMg+u81KT1GQyQ
2VWbZU1pOJ3PHtJPZcAJrppUCxrOu67TN0pgYIB7U0HHhoakqopuizAGRnf/
GzvkSUPlzYlVx1wYqHsSDXEvGoLpb1N/2TPwYrHstTI/Xr55YWKXrRgI/KC5
dk0wDekMmZlMEwaENI55z8Tz1v+3+XPfKt560T9LniTRoBW9WN+HwsBXdqiE
w20aRqXrc/jlefu9e6egPIMGm7UuSZSlDJz9oyJwNpuG5e/UhPMWMmBv/cpK
NYeGCtNf5w1mK2GQv+1hdz5v/ZaPE19/VIIm2jh9vYgGbedoV4exSgi6O243
KaZhZHxfO3egEsPlYzd/f6RB81mMhw6rEjWUoKGn5TRsimdY2PdXIve80AZH
Bg3b/RZrhrIrkdCeGCXxlQZLW2OpR0OVCFhP66hopOHQ5pCZmuFKYGxzkFoX
Dd4Cf9/LjVeCalpV1dPLy4ejm75lohICmTaKCWwaLtf4hrtOVaLa3ufdn3Ea
7iRyzfN+VcKu+tG8L/x0fFZmte2YzwBW6ewJWkxHwxJS8Ul+BihhpWnqInR8
Gzt4L2ERA0P6TMNEKTomXja79ggx4P9EIsxpFR2KJpWTAeIM2C2+3ySlRscq
FcGWOxIMbHZWV65cR8c60a2vS6UYWKRgWqGhT8eOluLg5bx6D/m3SPdtoWPv
m99O6xUZqGpwPppkTMfBuzomB4kM5KhPFJma0uEW6rM6hMzA9egLQv920eHl
liuSReXFZy+1zbWko7DWIvKjMgO2RinZzgfo8F/St9VflYFNd1f8kbanY71x
wDxVLd75+/PCjOFIx/SFpW97N/DWP865etOHjl/H7bNf7Wegcsd+cc0Unrfw
RbhFMzAd+bbB7C4dSTsoYs/iGaBWUhOcH9Chvdf41lQy77yZTkgnPqPD2yk6
52IWAw+jbNqe5dOx3OOZ3uccBhorP6RUFNGR61NXsqyQAVWzWMWZT3SMhUm3
3uL1j030zy6xz3TExGw40lvFQDjD/t7qajpUkuxHVjYykC9c5mhUT0f13WC/
kx28fjJbQz3YQod7Vvq8l70MCF+53ufTScfi3PKov0MM6Fb9zrjaQ0fWqyEJ
o+8MuCw94pY5QMf/fW83zvzLyvfDcy6+osFpGZ8zpyrp8fiPORvtdFWjzFP+
z1M7+bXCJOZsbbEJZpvmHBr7YJ6z4Zyf1wp/CjKZs9Du9q3P9sz5zW5/UzHX
OSvtzd3fEjPn/n10D6fuOeffDI6k98/5YmdHJnvo//3eOb7HffL/5efNZ+0j
sOI/f4lrMwhbM+cnNVelHwbM2Xv7z0+c5Sv/c2GE8QXfjav+s16biNfBC6v/
syqhxTHlqcp/3hE+M8nMm7PLjHw4vWjOd5gO2c8/zVn0CWeivGXO41v/hU7O
XzO3Xz96prnVnGuHTHTj7eY8fPB4ZYvDnGlGOWNH3OecIKqncyZ4zv5Z5p8z
sud8Xf6k7XDOnHNi4kY0CuY86NO8/O37OdsYONjUN87Z90XIsHT7nONWZJyz
Z875iwgnbXBozv3BS7XUuHOeP6VW7j01Z0W3PQdezcz5f7+P/p//B2gprEI=

          "]]}, Annotation[#, "Charting`Private`Tag$6084#3"]& ]}}, {}}, {
    DisplayFunction -> Identity, Ticks -> {Automatic, Automatic}, 
     AxesOrigin -> {0, 0}, FrameTicks -> {{Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}, {Automatic, 
        Charting`ScaledFrameTicks[{Identity, Identity}]}}, 
     GridLines -> {None, None}, DisplayFunction -> Identity, 
     PlotRangePadding -> {{
        Scaled[0.02], 
        Scaled[0.02]}, {
        Scaled[0.05], 
        Scaled[0.05]}}, PlotRangeClipping -> True, ImagePadding -> All, 
     DisplayFunction -> Identity, AspectRatio -> 
     NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, 
     AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, DisplayFunction :> 
     Identity, Frame -> {{False, False}, {False, False}}, 
     FrameLabel -> {{None, None}, {None, None}}, 
     FrameTicks -> {{Automatic, Automatic}, {Automatic, Automatic}}, 
     GridLines -> {None, None}, GridLinesStyle -> Directive[
       GrayLevel[0.5, 0.4]], ImageSize -> Large, 
     Method -> {
      "DefaultBoundaryStyle" -> Automatic, 
       "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 -> {All, All}, PlotRangeClipping -> 
     True, PlotRangePadding -> {{Automatic, Automatic}, {
       Automatic, Automatic}}, Ticks -> {Automatic, Automatic}}],FormBox[
    FormBox[
     TemplateBox[{"\"Heating\"", "\"Resonant\"", "\"Cooling\""}, "LineLegend",
       DisplayFunction -> (FormBox[
        StyleBox[
         StyleBox[
          PaneBox[
           TagBox[
            GridBox[{{
               TagBox[
                GridBox[{{
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}, {
                   GraphicsBox[{{
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {
                    LineBox[{{0, 10}, {20, 10}}]}}, {
                    Directive[
                    EdgeForm[
                    Directive[
                    Opacity[0.3], 
                    GrayLevel[0]]], 
                    PointSize[0.5], 
                    Opacity[1.], 
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    AbsoluteThickness[1.6]], {}}}, AspectRatio -> Full, 
                    ImageSize -> {20, 10}, PlotRangePadding -> None, 
                    ImagePadding -> Automatic, 
                    BaselinePosition -> (Scaled[0.1] -> Baseline)], #3}}, 
                 GridBoxAlignment -> {
                  "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, 
                 AutoDelete -> False, 
                 GridBoxDividers -> {
                  "Columns" -> {{False}}, "Rows" -> {{False}}}, 
                 GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}},
                  GridBoxSpacings -> {
                  "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, 
             GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, 
             AutoDelete -> False, 
             GridBoxItemSize -> {
              "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, 
             GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], 
            "Grid"], Alignment -> Left, AppearanceElements -> None, 
           ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> 
           "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], {
         FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> 
         False], TraditionalForm]& ), 
      InterpretationFunction :> (RowBox[{"LineLegend", "[", 
         RowBox[{
           RowBox[{"{", 
             RowBox[{
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.368417, 0.506779, 0.709798], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.24561133333333335`, 0.3378526666666667, 
                    0.4731986666666667], FrameTicks -> None, PlotRangePadding -> 
                    None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.368417`", ",", "0.506779`", ",", "0.709798`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.368417, 0.506779, 0.709798]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.368417, 0.506779, 0.709798], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.880722, 0.611041, 0.142051], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.587148, 0.40736066666666665`, 0.09470066666666668], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.880722`", ",", "0.611041`", ",", "0.142051`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.880722, 0.611041, 0.142051]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.880722, 0.611041, 0.142051], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}], 
               ",", 
               RowBox[{"Directive", "[", 
                 RowBox[{
                   RowBox[{"Opacity", "[", "1.`", "]"}], ",", 
                   InterpretationBox[
                    ButtonBox[
                    TooltipBox[
                    GraphicsBox[{{
                    GrayLevel[0], 
                    RectangleBox[{0, 0}]}, {
                    GrayLevel[0], 
                    RectangleBox[{1, -1}]}, {
                    RGBColor[0.560181, 0.691569, 0.194885], 
                    RectangleBox[{0, -1}, {2, 1}]}}, DefaultBaseStyle -> 
                    "ColorSwatchGraphics", AspectRatio -> 1, Frame -> True, 
                    FrameStyle -> 
                    RGBColor[
                    0.37345400000000006`, 0.461046, 0.12992333333333334`], 
                    FrameTicks -> None, PlotRangePadding -> None, ImageSize -> 
                    Dynamic[{
                    Automatic, 1.35 CurrentValue["FontCapHeight"]/
                    AbsoluteCurrentValue[Magnification]}]], 
                    StyleBox[
                    RowBox[{"RGBColor", "[", 
                    RowBox[{"0.560181`", ",", "0.691569`", ",", "0.194885`"}],
                     "]"}], NumberMarks -> False]], Appearance -> None, 
                    BaseStyle -> {}, BaselinePosition -> Baseline, 
                    DefaultBaseStyle -> {}, ButtonFunction :> 
                    With[{Typeset`box$ = EvaluationBox[]}, 
                    If[
                    Not[
                    AbsoluteCurrentValue["Deployed"]], 
                    SelectionMove[Typeset`box$, All, Expression]; 
                    FrontEnd`Private`$ColorSelectorInitialAlpha = 1; 
                    FrontEnd`Private`$ColorSelectorInitialColor = 
                    RGBColor[0.560181, 0.691569, 0.194885]; 
                    FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; 
                    MathLink`CallFrontEnd[
                    FrontEnd`AttachCell[Typeset`box$, 
                    FrontEndResource["RGBColorValueSelector"], {
                    0, {Left, Bottom}}, {Left, Top}, 
                    "ClosingActions" -> {
                    "SelectionDeparture", "ParentChanged", 
                    "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> 
                    Automatic, Method -> "Preemptive"], 
                    RGBColor[0.560181, 0.691569, 0.194885], Editable -> False,
                     Selectable -> False], ",", 
                   RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}]}], "]"}]}],
              "}"}], ",", 
           RowBox[{"{", 
             RowBox[{#, ",", #2, ",", #3}], "}"}], ",", 
           RowBox[{"LegendMarkers", "\[Rule]", "None"}], ",", 
           RowBox[{"LabelStyle", "\[Rule]", 
             RowBox[{"{", "}"}]}], ",", 
           RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), 
      Editable -> True], TraditionalForm], TraditionalForm]},
  "Legended",
  DisplayFunction->(GridBox[{{
      TagBox[
       ItemBox[
        PaneBox[
         TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, 
         BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], 
       "SkipImageSizeLevel"], 
      ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, 
    GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, 
    AutoDelete -> False, GridBoxItemSize -> Automatic, 
    BaselinePosition -> {1, 1}]& ),
  Editable->True,
  InterpretationFunction->(RowBox[{"Legended", "[", 
     RowBox[{#, ",", 
       RowBox[{"Placed", "[", 
         RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output",
 CellChangeTimes->{
  3.802798638744398*^9, 3.802971285931092*^9, 3.80297171438122*^9, 
   3.8029719939477596`*^9, 3.8029720839587283`*^9, 3.802972476302081*^9, {
   3.8029801902653875`*^9, 3.8029802076211166`*^9}, 3.8029803860980277`*^9, {
   3.8029812842386913`*^9, 3.802981303038907*^9}, {3.802981814715806*^9, 
   3.8029818258045416`*^9}, {3.8029826160868316`*^9, 
   3.8029826573802333`*^9}, {3.8029831280418425`*^9, 3.8029831411708956`*^9}, 
   3.802984641611767*^9, 3.8029847142051153`*^9, 3.802985573990672*^9, {
   3.802985875510088*^9, 3.8029859187481594`*^9}, {3.8029868207472744`*^9, 
   3.802986834113019*^9}, 3.8029880033625126`*^9, 3.8029935159942794`*^9, 
   3.8030083562936287`*^9, {3.803008415432871*^9, 3.803008449591169*^9}, {
   3.803008823503785*^9, 3.8030088785521936`*^9}, 3.803026242998581*^9, 
   3.8030267322757654`*^9, {3.8031452861803007`*^9, 3.803145298113488*^9}},
 CellLabel->"Out[46]=",ExpressionUUID->"af967df7-494f-4b14-947c-f734a8a9784e"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Scan the B-field from -5 G to 11 G", "Subsubsection",
 CellChangeTimes->{{3.7701329614446096`*^9, 3.7701329832348557`*^9}, {
   3.774612466089605*^9, 3.7746124720577526`*^9}, {3.7746226304517884`*^9, 
   3.774622638368577*^9}, {3.7746226690289564`*^9, 3.7746226695757027`*^9}, {
   3.776165362372452*^9, 3.776165395656355*^9}, 3.776165630526789*^9, {
   3.794530301368927*^9, 3.7945303021651297`*^9}, {3.8012752498082466`*^9, 
   3.801275295145632*^9}, {3.801275345493616*^9, 3.8012753500168557`*^9}, {
   3.8013471581823025`*^9, 3.8013471949624243`*^9}, {3.8018354238237987`*^9, 
   3.8018354252907934`*^9}, {3.801840382378128*^9, 3.801840383535265*^9}, {
   3.8018704362120605`*^9, 3.801870437440422*^9}, {3.8025983798226385`*^9, 
   3.8025983811701255`*^9}, {3.802598556275345*^9, 3.8025985575050144`*^9}, {
   3.8030083368478146`*^9, 
   3.803008342339774*^9}},ExpressionUUID->"c44a9341-abbe-437f-a7b8-\
e777560bd4ec"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Bfields", "=", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{"x", ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", 
       RowBox[{"-", "5"}], ",", "11"}], "}"}]}], "]"}], "+", "0.4", 
   RowBox[{"(*", 
    RowBox[{
    "we", " ", "add", " ", "an", " ", "offset", " ", "of", " ", "0.4", " ", 
     "G", " ", "to", " ", "account", " ", "for", " ", 
     RowBox[{"Earth", "'"}], "s", " ", "field"}], "*)"}]}]}]], "Input",
 CellChangeTimes->{{3.802378015195921*^9, 3.8023780618768625`*^9}, {
  3.802399859883916*^9, 3.8023998603848224`*^9}, {3.803008375918177*^9, 
  3.803008387992896*^9}, {3.8033055427446055`*^9, 
  3.8033055540222807`*^9}},ExpressionUUID->"765a6f84-5baa-41bd-a875-\
c58ef6d443ae"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"-", "4.6`"}], ",", 
   RowBox[{"-", "3.6`"}], ",", 
   RowBox[{"-", "2.6`"}], ",", 
   RowBox[{"-", "1.6`"}], ",", 
   RowBox[{"-", "0.6`"}], ",", "0.4`", ",", "1.4`", ",", "2.4`", ",", "3.4`", 
   ",", "4.4`", ",", "5.4`", ",", "6.4`", ",", "7.4`", ",", "8.4`", ",", 
   "9.4`", ",", "10.4`", ",", "11.4`"}], "}"}]], "Output",
 CellChangeTimes->{{3.802378057419541*^9, 3.8023780681395884`*^9}, 
   3.80239993288593*^9, 3.8024021412926517`*^9, 3.80244042565911*^9, {
   3.8024433875116386`*^9, 3.802443402101817*^9}, 3.8024436002123404`*^9, 
   3.802461535592783*^9, 3.8024741833102903`*^9, 3.8024744457399497`*^9, 
   3.8024747440804434`*^9, 3.8024747931504097`*^9, 3.8024750723683147`*^9, 
   3.8024751041491537`*^9, 3.8024753697470455`*^9, {3.802475964453163*^9, 
   3.8024759877631755`*^9}, 3.802477955373023*^9, 3.8024797254347687`*^9, 
   3.8025234781520853`*^9, 3.8025979961721478`*^9, 3.802598383670566*^9, 
   3.802598566853506*^9, 3.8025993049737854`*^9, 3.802599420688442*^9, 
   3.8026035479324446`*^9, 3.802604136557271*^9, 3.8029938393815603`*^9, 
   3.8030083888303456`*^9, 3.803008428199629*^9, 3.8030089755035496`*^9, 
   3.8030267626991615`*^9, 3.803031810780321*^9, 3.8030335749703083`*^9, 
   3.803145286420413*^9},
 CellLabel->"Out[41]=",ExpressionUUID->"7b6ac239-2586-43c2-8be4-68c57c1b26bc"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"FScale", "=", "0.50"}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "simple", " ", "rate", " ", "analysis", " ", "suggests", " ", "a", " ", 
     "correction", " ", "factor", " ", "2", 
     RowBox[{"ne", "/", 
      RowBox[{"(", 
       RowBox[{"ne", "+", "ng"}], ")"}]}]}], "=", 
    RowBox[{
     RowBox[{
      RowBox[{"1", "/", "2"}], " ", "due", " ", "to", " ", "multi"}], "-", 
     RowBox[{"level", " ", "structure"}]}]}], "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.803057618969963*^9, 3.803057619498081*^9}, {
   3.8030585409060755`*^9, 3.803058582983102*^9}, {3.8030603241077337`*^9, 
   3.8030604045619926`*^9}, {3.8030636999628677`*^9, 3.803063700069892*^9}, {
   3.8030637484673996`*^9, 3.8030638021806755`*^9}, {3.8030699528252354`*^9, 
   3.803069954377393*^9}, 3.803128637700624*^9, {3.8031452675433464`*^9, 
   3.8031452690495033`*^9}, {3.8033054870920877`*^9, 
   3.8033055023702726`*^9}},ExpressionUUID->"40ad01e0-8fb6-412a-8f61-\
4b35f584590c"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{
    RowBox[{"ainterp", "[", 
     RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
    RowBox[{
     RowBox[{"ainterp32", "[", 
      RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
     "2.25"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"compute", " ", "the", " ", "upper", " ", "bound"}], ",", " ", 
    RowBox[{
    "corresponding", " ", "to", " ", "a", " ", "2.25", " ", "times", " ", 
     "smaller", " ", "force", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"~", "180"}], " ", "rather", " ", 
       RowBox[{"than", " ", "~", "80"}], " ", "photons"}], ")"}]}]}], 
   "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{
   RowBox[{"rinterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"rinterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "2.25"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.8030604917691407`*^9, 3.8030605325130987`*^9}, {
  3.8030609143101935`*^9, 3.8030609143802094`*^9}, {3.8030636672580643`*^9, 
  3.803063668315302*^9}, {3.803063830642311*^9, 3.803063831469017*^9}, {
  3.803145279848048*^9, 3.803145284130925*^9}, {3.8033057327898016`*^9, 
  3.8033057382181454`*^9}},ExpressionUUID->"804e7a3d-4a13-49ca-81cd-\
0ea0ddd020d6"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"images1", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"First", "@", 
     RowBox[{"AbsoluteTiming", "[", 
      RowBox[{"cooled", "=", 
       RowBox[{"ParallelMap", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"SimulateTrajectoryFast", "[", 
           RowBox[{"#", ",", "Bfield", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", "*", 
              RowBox[{"800", "/", "isat"}]}], "]"}], ",", 
            RowBox[{"20", "/", "\[CapitalGamma]"}], ",", "180"}], "]"}], 
          "&"}], ",", "init"}], "]"}]}], "]"}]}], ";", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{
      "arguments", " ", "for", " ", "simulating", " ", "trajectories", " ", 
       RowBox[{"are", ":", " ", 
        RowBox[{"B", "-", "field"}]}]}], ",", " ", 
      RowBox[{"intensity", " ", 
       RowBox[{"(", 
        RowBox[{
        "2", "x", " ", "for", " ", "two", " ", "laser", " ", "beams"}], 
        ")"}]}], ",", " ", "detuning", ",", " ", 
      RowBox[{
      "and", " ", "average", " ", "#", " ", "photons", " ", "scattered"}]}], 
     "*)"}], "\[IndentingNewLine]", 
    RowBox[{"Mean", "[", 
     RowBox[{"cooled", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", "4"}], "]"}], "]"}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"cooledf", "=", 
     RowBox[{"Select", "[", 
      RowBox[{"cooled", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"#", "[", 
          RowBox[{"[", "5", "]"}], "]"}], "\[Equal]", "0"}], "&"}]}], "]"}]}],
     ";", "\[IndentingNewLine]", 
    RowBox[{"cooledimaged", "=", 
     RowBox[{"Select", "[", 
      RowBox[{"cooledf", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", 
          RowBox[{"#", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "rdetector"}], "&"}]}],
       "]"}]}], ";", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"images1", ",", "cooledimaged"}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"c1", "=", 
     RowBox[{"Length", "[", "cooledimaged", "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"SmoothHistogram", "[", 
      RowBox[{
       RowBox[{"cooledimaged", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", "0.001", ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", 
       RowBox[{"ImageSize", "\[Rule]", "Medium"}], ",", 
       RowBox[{"PlotRange", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "150"}], "}"}]}], "}"}]}], ",", 
       RowBox[{"Axes", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"True", ",", "True"}], "}"}]}], ",", 
       RowBox[{"ScalingFunctions", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"Identity", ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"#", "*", "c1"}], "&"}], ",", 
            RowBox[{
             RowBox[{"#", "/", "c1"}], "&"}]}], "}"}]}], "}"}]}]}], "]"}], 
     "]"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"Bfield", ",", "Bfields"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.7944038032008038`*^9, 3.794403827910453*^9}, {
   3.794403920340012*^9, 3.794403946682702*^9}, 3.7944039801799545`*^9, {
   3.794404078790533*^9, 3.794404085731656*^9}, {3.7944041379631567`*^9, 
   3.7944041509624677`*^9}, {3.794404224780773*^9, 3.7944042252380433`*^9}, {
   3.7944042806072035`*^9, 3.7944043037506313`*^9}, {3.794404355432145*^9, 
   3.7944043694529595`*^9}, {3.794404408351403*^9, 3.7944044111938033`*^9}, {
   3.7944044939003863`*^9, 3.7944044940756693`*^9}, {3.7944047408531065`*^9, 
   3.7944047422674427`*^9}, {3.794404777212698*^9, 3.7944047957653213`*^9}, {
   3.794404832928384*^9, 3.794404833193509*^9}, {3.794404869277589*^9, 
   3.7944048709071293`*^9}, {3.7944049223486595`*^9, 
   3.7944049733425093`*^9}, {3.7944050219432793`*^9, 
   3.7944050222423925`*^9}, {3.794405241266675*^9, 3.794405241352714*^9}, {
   3.794405291699999*^9, 3.794405295635968*^9}, {3.794405341608671*^9, 
   3.7944053424633093`*^9}, {3.7944053931197357`*^9, 3.7944053946027*^9}, {
   3.794406265296444*^9, 3.7944062654965453`*^9}, {3.7944064314303102`*^9, 
   3.794406431802267*^9}, {3.7944064758643723`*^9, 3.794406477226362*^9}, {
   3.7944070081982527`*^9, 3.794407008826415*^9}, {3.794407685798216*^9, 
   3.7944076876636014`*^9}, {3.794408702756246*^9, 3.794408720756444*^9}, {
   3.794409131122303*^9, 3.7944091312012205`*^9}, {3.7944092083583727`*^9, 
   3.794409211058243*^9}, 3.794409305889973*^9, {3.7944093859009943`*^9, 
   3.7944093862704945`*^9}, {3.7944094298659945`*^9, 
   3.7944094299601617`*^9}, {3.7944095382212815`*^9, 3.794409538368109*^9}, {
   3.794409685541462*^9, 3.794409685606153*^9}, {3.794409745826335*^9, 
   3.7944097459021335`*^9}, {3.7944098146212215`*^9, 
   3.7944098146771593`*^9}, {3.794409872469082*^9, 3.794409872545332*^9}, {
   3.794409949721157*^9, 3.794409956477316*^9}, {3.7944099983125935`*^9, 
   3.794409999117816*^9}, {3.7944100886040297`*^9, 3.794410101182096*^9}, {
   3.794410179842183*^9, 3.7944101801612816`*^9}, {3.7944102534876857`*^9, 
   3.7944102536181283`*^9}, {3.794410319055793*^9, 3.794410323690439*^9}, {
   3.79441082823357*^9, 3.79441083223265*^9}, {3.794410888563336*^9, 
   3.794410888645339*^9}, {3.7944109558748083`*^9, 3.794410956362916*^9}, {
   3.7944110169154344`*^9, 3.7944110172292547`*^9}, {3.7944156314039984`*^9, 
   3.794415632247884*^9}, {3.7944157809904475`*^9, 3.794415794746232*^9}, {
   3.7944159060668526`*^9, 3.7944159338799467`*^9}, {3.7944160245462465`*^9, 
   3.794416040202263*^9}, 3.794416274391674*^9, {3.7944164004988785`*^9, 
   3.7944164016408143`*^9}, {3.794417172930576*^9, 3.7944171742476263`*^9}, {
   3.7944172080348597`*^9, 3.7944172089930725`*^9}, {3.794417342365755*^9, 
   3.7944173425530586`*^9}, {3.794417400639179*^9, 3.79441740097908*^9}, {
   3.7944174710516615`*^9, 3.7944174796117516`*^9}, {3.7944175455501404`*^9, 
   3.7944175699649453`*^9}, {3.7944179061150904`*^9, 
   3.7944179064519157`*^9}, {3.7944179692596583`*^9, 
   3.7944179694270325`*^9}, {3.794418238394996*^9, 3.794418238470233*^9}, {
   3.79441832093019*^9, 3.794418322134024*^9}, {3.794418388936208*^9, 
   3.7944183892228327`*^9}, {3.7944184501769104`*^9, 3.794418451374069*^9}, {
   3.794418521794922*^9, 3.7944185225768595`*^9}, {3.794418588412101*^9, 
   3.7944185885945954`*^9}, {3.7944310886748285`*^9, 3.794431089484665*^9}, {
   3.794432573077078*^9, 3.794432576089105*^9}, {3.794433280923937*^9, 
   3.7944332842150507`*^9}, {3.794433682269907*^9, 3.794433682683097*^9}, {
   3.7944338156841516`*^9, 3.7944338157656584`*^9}, {3.7944338736185346`*^9, 
   3.7944338738794284`*^9}, {3.7944339482291546`*^9, 3.794433948972297*^9}, {
   3.794434556421958*^9, 3.7944345579455395`*^9}, {3.7944346067666216`*^9, 
   3.794434606864291*^9}, {3.7944353102988615`*^9, 3.7944353134568806`*^9}, {
   3.7944357159894037`*^9, 3.7944357167493553`*^9}, {3.7944357475452747`*^9, 
   3.794435782488765*^9}, {3.7944404712179456`*^9, 3.7944404712817264`*^9}, {
   3.794440982775674*^9, 3.794440986142047*^9}, {3.7944410285178275`*^9, 
   3.7944410291146646`*^9}, {3.794442439585224*^9, 3.7944424894807434`*^9}, {
   3.7944425217059593`*^9, 3.7944425259053965`*^9}, {3.7944425908812375`*^9, 
   3.7944426059566693`*^9}, {3.794442649831872*^9, 3.794442710327902*^9}, 
   3.794442781807767*^9, {3.794442988237514*^9, 3.7944429919029064`*^9}, {
   3.794448529796775*^9, 3.7944485424287424`*^9}, {3.794448895200401*^9, 
   3.794448895698408*^9}, {3.794492001722169*^9, 3.794492006918971*^9}, {
   3.794492725548942*^9, 3.794492725863717*^9}, {3.794527439855525*^9, 
   3.794527442545347*^9}, {3.794605039173044*^9, 3.794605103448554*^9}, {
   3.794648422840992*^9, 3.7946484230240335`*^9}, {3.794648550455897*^9, 
   3.7946485507009516`*^9}, {3.8012644175670223`*^9, 
   3.8012644184142113`*^9}, {3.8012700756628265`*^9, 3.801270077570756*^9}, 
   3.801347119141498*^9, {3.801356703875865*^9, 3.801356706094185*^9}, {
   3.801877911061778*^9, 3.8018779677003527`*^9}, {3.801958005764391*^9, 
   3.801958006387244*^9}, {3.802138195039628*^9, 3.8021381955969114`*^9}, {
   3.8021390514641275`*^9, 3.8021390799250855`*^9}, {3.8023780002612133`*^9, 
   3.8023780004977064`*^9}, {3.802443398665583*^9, 3.8024433990856752`*^9}, {
   3.802993673876672*^9, 3.8029936745584655`*^9}, {3.8029940721438527`*^9, 
   3.8029940736182594`*^9}, {3.8030083945685115`*^9, 
   3.8030084081631565`*^9}, {3.8031452870380826`*^9, 
   3.8031452882751083`*^9}, {3.80330555701101*^9, 
   3.803305595496028*^9}},ExpressionUUID->"176e990b-1be8-4872-88cb-\
8c4320ac0ca9"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"DoubleGaussianWithDip", "=", 
   RowBox[{"y0", "+", 
    RowBox[{"A1", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"Exp", "[", 
        RowBox[{"-", 
         RowBox[{"(", 
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"x", "-", "\[Mu]"}], ")"}], "2"], 
           RowBox[{"2", " ", 
            SuperscriptBox["\[Sigma]1", "2"]}]], ")"}]}], "]"}], "-", 
       RowBox[{"Exp", "[", 
        RowBox[{"-", 
         RowBox[{"(", 
          FractionBox[
           SuperscriptBox[
            RowBox[{"(", 
             RowBox[{"x", "-", "\[Mu]"}], ")"}], "2"], 
           RowBox[{"2", " ", 
            SuperscriptBox["\[Sigma]2", "2"]}]], ")"}]}], "]"}]}], ")"}]}], 
    "+", 
    RowBox[{"A2", " ", 
     RowBox[{"Exp", "[", 
      RowBox[{"-", 
       RowBox[{"(", 
        FractionBox[
         SuperscriptBox[
          RowBox[{"(", 
           RowBox[{"x", "-", "\[Mu]"}], ")"}], "2"], 
         RowBox[{"2", " ", 
          SuperscriptBox["\[Sigma]2", "2"]}]], ")"}]}], "]"}]}]}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"vars", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{"A1", ",", "50"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"A2", ",", "100"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"\[Mu]", ",", "0"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"\[Sigma]1", ",", "0.006"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"\[Sigma]2", ",", "0.001"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"y0", ",", "10"}], "}"}]}], "}"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"cons", "=", 
   RowBox[{
    RowBox[{"And", "@@", 
     RowBox[{"Thread", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"A1", ",", "A2", ",", "\[Sigma]1", ",", "\[Sigma]2"}], "}"}], 
       ">", "0"}], "]"}]}], "&&", 
    RowBox[{"(", 
     RowBox[{"\[Sigma]1", ">", "0.004"}], ")"}], "&&", 
    RowBox[{"(", 
     RowBox[{"\[Sigma]1", "<", "0.01"}], ")"}], "&&", 
    RowBox[{"(", 
     RowBox[{"\[Sigma]2", "<", "0.004"}], ")"}], "&&", 
    RowBox[{"(", 
     RowBox[{"\[Mu]", "<", "0.002"}], ")"}], "&&", 
    RowBox[{"(", 
     RowBox[{"\[Mu]", ">", 
      RowBox[{"-", "0.002"}]}], ")"}]}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.802009337376983*^9, 3.8020093374419975`*^9}, {
   3.802010151841791*^9, 3.8020101862889786`*^9}, {3.8020103270628843`*^9, 
   3.8020103512070637`*^9}, {3.802010521472079*^9, 3.8020105287677937`*^9}, {
   3.802010653114375*^9, 3.8020106562237*^9}, {3.80201078208329*^9, 
   3.8020107848931293`*^9}, {3.8020109498579597`*^9, 3.802010949932521*^9}, {
   3.8020110392182508`*^9, 3.802011039379287*^9}, {3.80201117538125*^9, 
   3.8020111782596684`*^9}, {3.8021784667594166`*^9, 
   3.8021784671823792`*^9}, {3.802178633721464*^9, 3.802178642932085*^9}, {
   3.802178864985571*^9, 3.802178868039163*^9}, {3.8021789156484275`*^9, 
   3.8021789198875136`*^9}, {3.802182035548355*^9, 3.8021820397462583`*^9}, {
   3.8021851850724673`*^9, 3.8021851893122177`*^9}, {3.8021861111142063`*^9, 
   3.802186171815131*^9}, {3.8021863196951156`*^9, 3.8021863199553423`*^9}, {
   3.802186647204237*^9, 3.8021866618144255`*^9}, {3.802187389622501*^9, 
   3.802187398162064*^9}, {3.8021874326719627`*^9, 3.80218743712914*^9}, 
   3.8030259402073555`*^9, 3.80302602936683*^9, {3.803026126224403*^9, 
   3.803026213119158*^9}, {3.8030262614552298`*^9, 3.8030262638144264`*^9}, {
   3.8030525559534807`*^9, 3.803052557860908*^9}},
 CellLabel->"In[50]:=",ExpressionUUID->"2f2b7678-2a05-401e-9ba8-46346cf1686a"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"cooledfractions", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"dist", "=", 
     RowBox[{"SmoothKernelDistribution", "[", 
      RowBox[{
       RowBox[{"image", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"InterpolationPoints", "\[Rule]", "10000"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"pdf", "=", 
     RowBox[{"PDF", "[", "dist", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"data", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"pdf", "[", "x", "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"-", "0.050"}], ",", "0.050", ",", "0.00001"}], "}"}]}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"fit", "=", 
     RowBox[{"NonlinearModelFit", "[", 
      RowBox[{"data", ",", 
       RowBox[{"{", 
        RowBox[{"DoubleGaussianWithDip", ",", "cons"}], "}"}], ",", "vars", 
       ",", "x", ",", 
       RowBox[{"MaxIterations", "\[Rule]", 
        SuperscriptBox["10", "6"]}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"fit", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"A2fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "2"}], "]"}], "]"}]}], ";", 
    RowBox[{"\[Mu]fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "2"}], "]"}], "]"}]}], ";", 
    RowBox[{"\[Sigma]2fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"a", "=", 
      RowBox[{"Integrate", "[", 
       RowBox[{
        RowBox[{"fit", "[", "x", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", 
          RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}]}], "]"}]}], "]"}], 
    ";", 
    RowBox[{"(*", 
     RowBox[{
     "The", " ", "total", " ", "area", " ", "of", " ", "the", " ", "trace"}], 
     "*)"}], "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"b", "=", 
      RowBox[{"Integrate", "[", 
       RowBox[{
        RowBox[{"A2fit", " ", 
         RowBox[{"Exp", "[", 
          RowBox[{"-", 
           RowBox[{"(", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"x", "-", "\[Mu]fit"}], ")"}], "2"], 
             RowBox[{"2", " ", 
              SuperscriptBox["\[Sigma]2fit", "2"]}]], ")"}]}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", 
          RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}]}], "]"}]}], "]"}], 
    ";", 
    RowBox[{"(*", 
     RowBox[{
     "The", " ", "area", " ", "of", " ", "the", " ", "narrow", " ", 
      "Gaussian", " ", "representing", " ", "the", " ", "cooled", " ", 
      "fraction"}], "*)"}], " ", "\[IndentingNewLine]", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"cooledfractions", ",", 
      RowBox[{"b", "/", "a"}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", "\[IndentingNewLine]", 
     RowBox[{"Show", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"Plot", "[", 
        RowBox[{
         RowBox[{"fit", "[", "x", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"x", ",", 
           RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}], ",", 
         RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
         RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"ListPlot", "[", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"x", ",", 
            RowBox[{"pdf", "[", "x", "]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"x", ",", 
            RowBox[{"-", "0.015"}], ",", "0.015", ",", "0.00001"}], "}"}]}], 
         "]"}], "]"}]}], "]"}], "\[IndentingNewLine]", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"image", ",", "images1"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.8020093640257497`*^9, 3.8020093871496444`*^9}, {
   3.8020094378748198`*^9, 3.802009439301165*^9}, {3.8020094710679483`*^9, 
   3.802009471304001*^9}, 3.802009528939002*^9, {3.8020101890726213`*^9, 
   3.8020101903826156`*^9}, {3.8021852976902723`*^9, 3.80218531546032*^9}, {
   3.802185597048528*^9, 3.8021856941299663`*^9}, {3.8021857745025454`*^9, 
   3.8021858196795197`*^9}, {3.8021862023730745`*^9, 
   3.8021862130277348`*^9}, {3.8021872792380075`*^9, 3.802187309429247*^9}, {
   3.8022857667939787`*^9, 3.802285777497543*^9}, {3.802361350066171*^9, 
   3.8023614171056604`*^9}, {3.8023618350398703`*^9, 3.802361835370702*^9}, {
   3.8023620560876904`*^9, 3.8023620564436803`*^9}, {3.8024066269262953`*^9, 
   3.8024066270406895`*^9}, 3.8030258733750916`*^9, {3.803026463791092*^9, 
   3.8030264638601074`*^9}, {3.8030317268642936`*^9, 3.803031726919306*^9}, {
   3.8030334719662538`*^9, 3.803033472060275*^9}, 3.8030388236525216`*^9},
 CellLabel->"In[53]:=",ExpressionUUID->"c0602803-9f61-43c8-9b56-6982ac3d8524"],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Print", "[", 
    RowBox[{"cooledfractions", "[", 
     RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"i", ",", 
     RowBox[{"Table", "[", 
      RowBox[{"x", ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", "1", ",", "17"}], "}"}]}], "]"}]}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8030263370963697`*^9, 3.8030264095937185`*^9}, {
  3.80302684384328*^9, 3.8030268442273674`*^9}},
 CellLabel->"In[55]:=",ExpressionUUID->"93e65e24-b5de-4980-a1bf-6b959bad002d"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ainterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"ainterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "0.625"}]}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"compute", " ", "the", " ", "lower", " ", "bound"}], ",", " ", 
    RowBox[{
    "corresponding", " ", "to", " ", "a", " ", "0.625", " ", "times", " ", 
     "smaller", " ", "force", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"~", "50"}], " ", "rather", " ", 
       RowBox[{"than", " ", "~", "80"}], " ", "photons"}], ")"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"rinterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"rinterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "0.625"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.803146173054251*^9, 3.8031461781790743`*^9}},
 CellLabel->"In[56]:=",ExpressionUUID->"de57bc3e-ca7d-4974-8fd9-12835307d824"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"images2", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"First", "@", 
     RowBox[{"AbsoluteTiming", "[", 
      RowBox[{"cooled", "=", 
       RowBox[{"ParallelMap", "[", 
        RowBox[{
         RowBox[{
          RowBox[{"SimulateTrajectoryFast", "[", 
           RowBox[{"#", ",", "Bfield", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", "*", 
              RowBox[{"800", "/", "isat"}]}], "]"}], ",", 
            RowBox[{"20", "/", "\[CapitalGamma]"}], ",", "50"}], "]"}], "&"}],
          ",", "init"}], "]"}]}], "]"}]}], ";", 
    RowBox[{"(*", 
     RowBox[{
      RowBox[{
       RowBox[{"arguments", " ", "are", " ", "B"}], "=", 
       RowBox[{"1.0", " ", "G"}]}], ",", " ", 
      RowBox[{"s0", "=", 
       RowBox[{"92", " ", 
        RowBox[{"(", 
         RowBox[{"corresponding", " ", "to", " ", "400", " ", 
          RowBox[{"mW", "/", 
           RowBox[{"cm", "^", "2"}]}]}], ")"}]}]}]}], "*)"}], 
    "\[IndentingNewLine]", 
    RowBox[{"Mean", "[", 
     RowBox[{"cooled", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", "4"}], "]"}], "]"}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"cooledf", "=", 
     RowBox[{"Select", "[", 
      RowBox[{"cooled", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"#", "[", 
          RowBox[{"[", "5", "]"}], "]"}], "\[Equal]", "0"}], "&"}]}], "]"}]}],
     ";", "\[IndentingNewLine]", 
    RowBox[{"cooledimaged", "=", 
     RowBox[{"Select", "[", 
      RowBox[{"cooledf", ",", 
       RowBox[{
        RowBox[{
         RowBox[{"Abs", "[", 
          RowBox[{"#", "[", 
           RowBox[{"[", "1", "]"}], "]"}], "]"}], "<", "rdetector"}], "&"}]}],
       "]"}]}], ";", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"images2", ",", "cooledimaged"}], "]"}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"c1", "=", 
     RowBox[{"Length", "[", "cooledimaged", "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"SmoothHistogram", "[", 
      RowBox[{
       RowBox[{"cooledimaged", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", "0.001", ",", 
       RowBox[{"PlotStyle", "\[Rule]", "Blue"}], ",", 
       RowBox[{"ImageSize", "\[Rule]", "Medium"}], ",", 
       RowBox[{"PlotRange", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{
            RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"0", ",", "150"}], "}"}]}], "}"}]}], ",", 
       RowBox[{"Axes", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"True", ",", "True"}], "}"}]}], ",", 
       RowBox[{"ScalingFunctions", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{"Identity", ",", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{
             RowBox[{"#", "*", "c1"}], "&"}], ",", 
            RowBox[{
             RowBox[{"#", "/", "c1"}], "&"}]}], "}"}]}], "}"}]}]}], "]"}], 
     "]"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"Bfield", ",", "Bfields"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.7944038032008038`*^9, 3.794403827910453*^9}, {
   3.794403920340012*^9, 3.794403946682702*^9}, 3.7944039801799545`*^9, {
   3.794404078790533*^9, 3.794404085731656*^9}, {3.7944041379631567`*^9, 
   3.7944041509624677`*^9}, {3.794404224780773*^9, 3.7944042252380433`*^9}, {
   3.7944042806072035`*^9, 3.7944043037506313`*^9}, {3.794404355432145*^9, 
   3.7944043694529595`*^9}, {3.794404408351403*^9, 3.7944044111938033`*^9}, {
   3.7944044939003863`*^9, 3.7944044940756693`*^9}, {3.7944047408531065`*^9, 
   3.7944047422674427`*^9}, {3.794404777212698*^9, 3.7944047957653213`*^9}, {
   3.794404832928384*^9, 3.794404833193509*^9}, {3.794404869277589*^9, 
   3.7944048709071293`*^9}, {3.7944049223486595`*^9, 
   3.7944049733425093`*^9}, {3.7944050219432793`*^9, 
   3.7944050222423925`*^9}, {3.794405241266675*^9, 3.794405241352714*^9}, {
   3.794405291699999*^9, 3.794405295635968*^9}, {3.794405341608671*^9, 
   3.7944053424633093`*^9}, {3.7944053931197357`*^9, 3.7944053946027*^9}, {
   3.794406265296444*^9, 3.7944062654965453`*^9}, {3.7944064314303102`*^9, 
   3.794406431802267*^9}, {3.7944064758643723`*^9, 3.794406477226362*^9}, {
   3.7944070081982527`*^9, 3.794407008826415*^9}, {3.794407685798216*^9, 
   3.7944076876636014`*^9}, {3.794408702756246*^9, 3.794408720756444*^9}, {
   3.794409131122303*^9, 3.7944091312012205`*^9}, {3.7944092083583727`*^9, 
   3.794409211058243*^9}, 3.794409305889973*^9, {3.7944093859009943`*^9, 
   3.7944093862704945`*^9}, {3.7944094298659945`*^9, 
   3.7944094299601617`*^9}, {3.7944095382212815`*^9, 3.794409538368109*^9}, {
   3.794409685541462*^9, 3.794409685606153*^9}, {3.794409745826335*^9, 
   3.7944097459021335`*^9}, {3.7944098146212215`*^9, 
   3.7944098146771593`*^9}, {3.794409872469082*^9, 3.794409872545332*^9}, {
   3.794409949721157*^9, 3.794409956477316*^9}, {3.7944099983125935`*^9, 
   3.794409999117816*^9}, {3.7944100886040297`*^9, 3.794410101182096*^9}, {
   3.794410179842183*^9, 3.7944101801612816`*^9}, {3.7944102534876857`*^9, 
   3.7944102536181283`*^9}, {3.794410319055793*^9, 3.794410323690439*^9}, {
   3.79441082823357*^9, 3.79441083223265*^9}, {3.794410888563336*^9, 
   3.794410888645339*^9}, {3.7944109558748083`*^9, 3.794410956362916*^9}, {
   3.7944110169154344`*^9, 3.7944110172292547`*^9}, {3.7944156314039984`*^9, 
   3.794415632247884*^9}, {3.7944157809904475`*^9, 3.794415794746232*^9}, {
   3.7944159060668526`*^9, 3.7944159338799467`*^9}, {3.7944160245462465`*^9, 
   3.794416040202263*^9}, 3.794416274391674*^9, {3.7944164004988785`*^9, 
   3.7944164016408143`*^9}, {3.794417172930576*^9, 3.7944171742476263`*^9}, {
   3.7944172080348597`*^9, 3.7944172089930725`*^9}, {3.794417342365755*^9, 
   3.7944173425530586`*^9}, {3.794417400639179*^9, 3.79441740097908*^9}, {
   3.7944174710516615`*^9, 3.7944174796117516`*^9}, {3.7944175455501404`*^9, 
   3.7944175699649453`*^9}, {3.7944179061150904`*^9, 
   3.7944179064519157`*^9}, {3.7944179692596583`*^9, 
   3.7944179694270325`*^9}, {3.794418238394996*^9, 3.794418238470233*^9}, {
   3.79441832093019*^9, 3.794418322134024*^9}, {3.794418388936208*^9, 
   3.7944183892228327`*^9}, {3.7944184501769104`*^9, 3.794418451374069*^9}, {
   3.794418521794922*^9, 3.7944185225768595`*^9}, {3.794418588412101*^9, 
   3.7944185885945954`*^9}, {3.7944310886748285`*^9, 3.794431089484665*^9}, {
   3.794432573077078*^9, 3.794432576089105*^9}, {3.794433280923937*^9, 
   3.7944332842150507`*^9}, {3.794433682269907*^9, 3.794433682683097*^9}, {
   3.7944338156841516`*^9, 3.7944338157656584`*^9}, {3.7944338736185346`*^9, 
   3.7944338738794284`*^9}, {3.7944339482291546`*^9, 3.794433948972297*^9}, {
   3.794434556421958*^9, 3.7944345579455395`*^9}, {3.7944346067666216`*^9, 
   3.794434606864291*^9}, {3.7944353102988615`*^9, 3.7944353134568806`*^9}, {
   3.7944357159894037`*^9, 3.7944357167493553`*^9}, {3.7944357475452747`*^9, 
   3.794435782488765*^9}, {3.7944404712179456`*^9, 3.7944404712817264`*^9}, {
   3.794440982775674*^9, 3.794440986142047*^9}, {3.7944410285178275`*^9, 
   3.7944410291146646`*^9}, {3.794442439585224*^9, 3.7944424894807434`*^9}, {
   3.7944425217059593`*^9, 3.7944425259053965`*^9}, {3.7944425908812375`*^9, 
   3.7944426059566693`*^9}, {3.794442649831872*^9, 3.794442710327902*^9}, 
   3.794442781807767*^9, {3.794442988237514*^9, 3.7944429919029064`*^9}, {
   3.794448529796775*^9, 3.7944485424287424`*^9}, {3.794448895200401*^9, 
   3.794448895698408*^9}, {3.794492001722169*^9, 3.794492006918971*^9}, {
   3.794492725548942*^9, 3.794492725863717*^9}, {3.794527439855525*^9, 
   3.794527442545347*^9}, {3.794605039173044*^9, 3.794605103448554*^9}, {
   3.794648422840992*^9, 3.7946484230240335`*^9}, {3.794648550455897*^9, 
   3.7946485507009516`*^9}, {3.8012644175670223`*^9, 
   3.8012644184142113`*^9}, {3.8012700756628265`*^9, 3.801270077570756*^9}, 
   3.801347119141498*^9, {3.801356703875865*^9, 3.801356706094185*^9}, {
   3.801877911061778*^9, 3.8018779677003527`*^9}, {3.801958005764391*^9, 
   3.801958006387244*^9}, {3.802138195039628*^9, 3.8021381955969114`*^9}, {
   3.8021390514641275`*^9, 3.8021390799250855`*^9}, {3.8023780002612133`*^9, 
   3.8023780004977064`*^9}, {3.802443398665583*^9, 3.8024433990856752`*^9}, {
   3.802993673876672*^9, 3.8029936745584655`*^9}, {3.8029940721438527`*^9, 
   3.8029940736182594`*^9}, {3.8030083945685115`*^9, 
   3.8030084081631565`*^9}, {3.8030089707946835`*^9, 
   3.8030089725499973`*^9}, {3.803148281659013*^9, 3.803148281807537*^9}},
 CellLabel->"In[57]:=",ExpressionUUID->"dd10abc8-8e7a-4d4f-b878-7b142f132232"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"cooledfractions", "=", 
   RowBox[{"{", "}"}]}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"dist", "=", 
     RowBox[{"SmoothKernelDistribution", "[", 
      RowBox[{
       RowBox[{"image", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"InterpolationPoints", "\[Rule]", "10000"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"pdf", "=", 
     RowBox[{"PDF", "[", "dist", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"data", "=", 
     RowBox[{"Table", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"pdf", "[", "x", "]"}]}], "}"}], ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", 
         RowBox[{"-", "0.050"}], ",", "0.050", ",", "0.00001"}], "}"}]}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"fit", "=", 
     RowBox[{"NonlinearModelFit", "[", 
      RowBox[{"data", ",", 
       RowBox[{"{", 
        RowBox[{"DoubleGaussianWithDip", ",", "cons"}], "}"}], ",", "vars", 
       ",", "x", ",", 
       RowBox[{"MaxIterations", "\[Rule]", 
        SuperscriptBox["10", "6"]}]}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"fit", "[", 
      RowBox[{"[", "1", "]"}], "]"}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"A2fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"2", ",", "2"}], "]"}], "]"}]}], ";", 
    RowBox[{"\[Mu]fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"3", ",", "2"}], "]"}], "]"}]}], ";", 
    RowBox[{"\[Sigma]2fit", "=", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<BestFitParameters\>\"", "]"}], "[", 
      RowBox[{"[", 
       RowBox[{"5", ",", "2"}], "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"a", "=", 
      RowBox[{"Integrate", "[", 
       RowBox[{
        RowBox[{"fit", "[", "x", "]"}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", 
          RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}]}], "]"}]}], "]"}], 
    ";", 
    RowBox[{"(*", 
     RowBox[{
     "The", " ", "total", " ", "area", " ", "of", " ", "the", " ", "trace"}], 
     "*)"}], "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"b", "=", 
      RowBox[{"Integrate", "[", 
       RowBox[{
        RowBox[{"A2fit", " ", 
         RowBox[{"Exp", "[", 
          RowBox[{"-", 
           RowBox[{"(", 
            FractionBox[
             SuperscriptBox[
              RowBox[{"(", 
               RowBox[{"x", "-", "\[Mu]fit"}], ")"}], "2"], 
             RowBox[{"2", " ", 
              SuperscriptBox["\[Sigma]2fit", "2"]}]], ")"}]}], "]"}]}], ",", 
        RowBox[{"{", 
         RowBox[{"x", ",", 
          RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}]}], "]"}]}], "]"}], 
    ";", 
    RowBox[{"(*", 
     RowBox[{
     "The", " ", "area", " ", "of", " ", "the", " ", "narrow", " ", 
      "Gaussian", " ", "representing", " ", "the", " ", "cooled", " ", 
      "fraction"}], "*)"}], " ", "\[IndentingNewLine]", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"cooledfractions", ",", 
      RowBox[{"b", "/", "a"}]}], "]"}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", "\[IndentingNewLine]", 
     RowBox[{"Show", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"Plot", "[", 
        RowBox[{
         RowBox[{"fit", "[", "x", "]"}], ",", 
         RowBox[{"{", 
          RowBox[{"x", ",", 
           RowBox[{"-", "0.015"}], ",", "0.015"}], "}"}], ",", 
         RowBox[{"PlotRange", "\[Rule]", "All"}], ",", 
         RowBox[{"PlotStyle", "\[Rule]", "Red"}]}], "]"}], ",", 
       "\[IndentingNewLine]", 
       RowBox[{"ListPlot", "[", 
        RowBox[{"Table", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"x", ",", 
            RowBox[{"pdf", "[", "x", "]"}]}], "}"}], ",", 
          RowBox[{"{", 
           RowBox[{"x", ",", 
            RowBox[{"-", "0.015"}], ",", "0.015", ",", "0.00001"}], "}"}]}], 
         "]"}], "]"}]}], "]"}], "\[IndentingNewLine]", "]"}]}], ",", 
   RowBox[{"{", 
    RowBox[{"image", ",", "images2"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.8020093640257497`*^9, 3.8020093871496444`*^9}, {
   3.8020094378748198`*^9, 3.802009439301165*^9}, {3.8020094710679483`*^9, 
   3.802009471304001*^9}, 3.802009528939002*^9, {3.8020101890726213`*^9, 
   3.8020101903826156`*^9}, {3.8021852976902723`*^9, 3.80218531546032*^9}, {
   3.802185597048528*^9, 3.8021856941299663`*^9}, {3.8021857745025454`*^9, 
   3.8021858196795197`*^9}, {3.8021862023730745`*^9, 
   3.8021862130277348`*^9}, {3.8021872792380075`*^9, 3.802187309429247*^9}, {
   3.8022857667939787`*^9, 3.802285777497543*^9}, {3.802361350066171*^9, 
   3.8023614171056604`*^9}, {3.8023618350398703`*^9, 3.802361835370702*^9}, {
   3.8023620560876904`*^9, 3.8023620564436803`*^9}, {3.8024066269262953`*^9, 
   3.8024066270406895`*^9}, 3.8030258733750916`*^9, {3.803026463791092*^9, 
   3.8030264638601074`*^9}, {3.8030317268642936`*^9, 3.803031726919306*^9}, {
   3.8030334719662538`*^9, 3.803033472060275*^9}, 3.8030388236525216`*^9, 
   3.8030830868914633`*^9},
 CellLabel->"In[59]:=",ExpressionUUID->"8114984e-8954-455e-9fc9-850c08ff99c0"],

Cell[BoxData[
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"Print", "[", 
    RowBox[{"cooledfractions", "[", 
     RowBox[{"[", "i", "]"}], "]"}], "]"}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"i", ",", 
     RowBox[{"Table", "[", 
      RowBox[{"x", ",", 
       RowBox[{"{", 
        RowBox[{"x", ",", "1", ",", "17"}], "}"}]}], "]"}]}], "}"}]}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.8030263370963697`*^9, 3.8030264095937185`*^9}, {
  3.80302684384328*^9, 3.8030268442273674`*^9}},
 CellLabel->"In[61]:=",ExpressionUUID->"60189437-ff43-4754-8c0b-9895904a8ad1"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{2560, 1357},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>1.2 Inherited,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 238, 4, 35, "Input",ExpressionUUID->"acf3f024-06f7-4b67-aa47-b28ae5a429d7"],
Cell[CellGroupData[{
Cell[821, 28, 493, 8, 53, "Subsubsection",ExpressionUUID->"5aaf2c87-6956-45b5-ba7b-95459d133cb4"],
Cell[1317, 38, 20831, 396, 547, "Input",ExpressionUUID->"7824fddb-3943-46b5-acd8-799c34acddf8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22185, 439, 415, 6, 65, "Subsection",ExpressionUUID->"c541093d-3069-439e-a54a-d38289364e3b"],
Cell[22603, 447, 5979, 179, 372, "Input",ExpressionUUID->"a22cab01-84b6-4910-a8ee-e13c7c9bda3f"],
Cell[CellGroupData[{
Cell[28607, 630, 496, 7, 53, "Subsubsection",ExpressionUUID->"2086d6f6-0d0b-40d9-a62f-2851fbd41666"],
Cell[29106, 639, 4063, 95, 420, "Input",ExpressionUUID->"faefb1d2-c143-4cf4-bbad-8f8982efcabd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33206, 739, 524, 7, 53, "Subsubsection",ExpressionUUID->"20effcbb-203d-40b4-8029-107270944d9f"],
Cell[33733, 748, 9691, 192, 557, "Input",ExpressionUUID->"dccf9f9d-0271-4cee-a2a6-51b034d09fcc"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43461, 945, 578, 8, 53, "Subsubsection",ExpressionUUID->"6c519e3a-f10c-4538-afc2-24105e47eb71"],
Cell[44042, 955, 2780, 39, 35, "Input",ExpressionUUID->"9e84a53b-6a94-415f-8b47-20c1efb51727"],
Cell[46825, 996, 9308, 169, 372, "Input",ExpressionUUID->"c0c8e4fa-0227-4b9a-82f2-c524254ed188"],
Cell[56136, 1167, 3032, 46, 35, "Input",ExpressionUUID->"c5609752-0026-4bbc-8872-3b54364462a9"],
Cell[CellGroupData[{
Cell[59193, 1217, 15582, 230, 84, "Input",ExpressionUUID->"9c21872c-cc44-4cd6-9620-54ca4e5713c1"],
Cell[74778, 1449, 48868, 871, 450, "Output",ExpressionUUID->"af967df7-494f-4b14-947c-f734a8a9784e"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[123695, 2326, 928, 13, 53, "Subsubsection",ExpressionUUID->"c44a9341-abbe-437f-a7b8-e777560bd4ec"],
Cell[CellGroupData[{
Cell[124648, 2343, 725, 17, 35, "Input",ExpressionUUID->"765a6f84-5baa-41bd-a875-c58ef6d443ae"],
Cell[125376, 2362, 1366, 23, 39, "Output",ExpressionUUID->"7b6ac239-2586-43c2-8be4-68c57c1b26bc"]
}, Open  ]],
Cell[126757, 2388, 1011, 22, 35, "Input",ExpressionUUID->"40ad01e0-8fb6-412a-8f61-4b35f584590c"],
Cell[127771, 2412, 1364, 34, 60, "Input",ExpressionUUID->"804e7a3d-4a13-49ca-81cd-0ea0ddd020d6"],
Cell[129138, 2448, 8861, 169, 252, "Input",ExpressionUUID->"176e990b-1be8-4872-88cb-8c4320ac0ca9"],
Cell[138002, 2619, 3594, 91, 111, "Input",ExpressionUUID->"2f2b7678-2a05-401e-9ba8-46346cf1686a"],
Cell[141599, 2712, 5337, 129, 424, "Input",ExpressionUUID->"c0602803-9f61-43c8-9b56-6982ac3d8524"],
Cell[146939, 2843, 608, 15, 84, "Input",ExpressionUUID->"93e65e24-b5de-4980-a1bf-6b959bad002d"],
Cell[147550, 2860, 1099, 28, 60, "Input",ExpressionUUID->"de57bc3e-ca7d-4974-8fd9-12835307d824"],
Cell[148652, 2890, 8756, 166, 228, "Input",ExpressionUUID->"dd10abc8-8e7a-4d4f-b878-7b142f132232"],
Cell[157411, 3058, 5365, 130, 424, "Input",ExpressionUUID->"8114984e-8954-455e-9fc9-850c08ff99c0"],
Cell[162779, 3190, 608, 15, 84, "Input",ExpressionUUID->"60189437-ff43-4754-8c0b-9895904a8ad1"]
}, Open  ]]
}, Open  ]]
}
]
*)