TrajectoryMonteCarlo_IntensityScan_Heating.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[    155071,       3013]
NotebookOptionsPosition[    152104,       2959]
NotebookOutlinePosition[    152480,       2975]
CellTagsIndexPosition[    152437,       2972]
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_\
IntensityScan_Heating.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.8024604615931597`*^9, 3.802460462602696*^9}, {
   3.8024606886005917`*^9, 3.802460693720007*^9}, {3.8025241997037287`*^9, 
   3.8025242010632315`*^9}, {3.802601041649931*^9, 3.802601041835507*^9}, 
   3.8026022235362706`*^9, {3.8029665316506553`*^9, 3.8029665521514*^9}, 
   3.8029708884425173`*^9, {3.8029729099896936`*^9, 3.8029729114660296`*^9}, {
   3.8029799037006073`*^9, 3.8029799101612024`*^9}, {3.802981706395821*^9, 
   3.802981708037957*^9}, {3.802984167405116*^9, 3.8029841705401955`*^9}, {
   3.8029844858643*^9, 3.802984488527766*^9}, {3.8029852443625526`*^9, 
   3.8029852459628906`*^9}, 3.802985502581078*^9, {3.8029874407260776`*^9, 
   3.802987440814825*^9}, {3.8029875606705155`*^9, 3.8029876249280014`*^9}, 
   3.8029923468356586`*^9, {3.8030564073566775`*^9, 3.8030564105868587`*^9}, {
   3.803057712798812*^9, 3.8030577150973263`*^9}, 3.803058247813226*^9, {
   3.8030585026698112`*^9, 3.8030585035741014`*^9}, {3.803059444022292*^9, 
   3.8030594450832415`*^9}, 3.8030598687318463`*^9, {3.8030698525587435`*^9, 
   3.803069853425481*^9}, {3.8030829321865726`*^9, 3.8030829352561846`*^9}, {
   3.803133963241295*^9, 3.8031339815530157`*^9}, {3.8031393592448792`*^9, 
   3.803139359878765*^9}, {3.8033058399234657`*^9, 
   3.8033058432892623`*^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.8030571816931486`*^9, 
   3.8030571828000093`*^9}, 3.803057895782549*^9, {3.8030580168970795`*^9, 
   3.803058018284564*^9}, 3.803059449363638*^9, {3.8030598720455894`*^9, 
   3.8030598732495127`*^9}, {3.8030599339635344`*^9, 3.803059937215263*^9}, {
   3.8030599865906*^9, 3.8030599910080233`*^9}, {3.803060436261153*^9, 
   3.803060437505432*^9}, {3.803060898648136*^9, 3.803060899542336*^9}, 
   3.803063944186997*^9, {3.803091733879302*^9, 3.8030917618520513`*^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.8030609030161147`*^9, 3.8030609030671263`*^9}, {
   3.8031340021236196`*^9, 3.8031340022586703`*^9}},
 CellLabel->"In[23]:=",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.8024780585779605`*^9, 3.802478058827526*^9}, {
   3.8031339898644867`*^9, 3.803133990167058*^9}, {3.8031553623894634`*^9, 
   3.803155362583469*^9}},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.802478068766712*^9, 
   3.8024780737609825`*^9}, {3.802478137009843*^9, 3.802478162609681*^9}, {
   3.802987390658433*^9, 3.8029873910361977`*^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.8024607225509233`*^9, 
   3.802460722781674*^9}, {3.8024727213521166`*^9, 3.8024727221712637`*^9}, {
   3.802478078722689*^9, 3.8024780787928314`*^9}, 3.8026010670277543`*^9, {
   3.802966571249402*^9, 3.8029665712724714`*^9}, {3.8029672296118884`*^9, 
   3.8029672298993974`*^9}, {3.802973165738646*^9, 3.802973165798636*^9}, {
   3.8029734785536375`*^9, 3.8029734939981737`*^9}, 3.8030564246685567`*^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[{"ainterp32", "[", 
      RowBox[{"1.4", ",", 
       RowBox[{"QuantityMagnitude", "[", 
        RowBox[{"2", "*", 
         RowBox[{"6000", "/", "isat"}]}], "]"}], ",", 
       RowBox[{"{", 
        RowBox[{
         RowBox[{
          RowBox[{"-", "15"}], "/", "\[CapitalGamma]"}], ",", "0", ",", 
         RowBox[{"25", "/", "\[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.8024727296553793`*^9, 
   3.802472733400488*^9}, {3.802473373301007*^9, 3.802473373358121*^9}, {
   3.802967477569319*^9, 3.802967479113695*^9}, {3.8029679771608706`*^9, 
   3.802967981080007*^9}, {3.80296912357218*^9, 3.8029691278240566`*^9}, {
   3.802970319359497*^9, 3.8029703194635897`*^9}, {3.802973834661152*^9, 
   3.8029738350220966`*^9}, {3.802974084774898*^9, 3.8029740900961375`*^9}, {
   3.8029745901990256`*^9, 3.8029745905640163`*^9}, {3.802974675116275*^9, 
   3.802974675524158*^9}, {3.802976772246581*^9, 3.802976784546658*^9}, {
   3.802982413240527*^9, 3.8029824137940683`*^9}, {3.8029838328951683`*^9, 
   3.8029839015199986`*^9}, {3.8029841504226356`*^9, 
   3.8029842001373396`*^9}, {3.802984264054482*^9, 3.802984264504074*^9}, {
   3.802987423069471*^9, 3.802987435655458*^9}, {3.80298757580832*^9, 
   3.802987579306053*^9}, {3.8030596289801397`*^9, 3.803059629579274*^9}, {
   3.8030597447508287`*^9, 3.8030597449878817`*^9}, {3.8030598329479556`*^9, 
   3.803059833017971*^9}, {3.8030635789804564`*^9, 3.8030635790864806`*^9}, {
   3.803134063391732*^9, 3.8031340634533777`*^9}},
 CellLabel->"In[39]:=",ExpressionUUID->"e65b6191-864b-4f9a-bf41-1a27cb9ce5bb"],

Cell[BoxData[
 TemplateBox[{GraphicsBox[{{{{}, {}, 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.368417, 0.506779, 0.709798], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJxE13c01+/7OHB77733qtCiVOp1vW0NpCjKKiKhYVRCQkZJ9t4rK3uWvC4i
pQihQrbsHdm+fuf8zqfXP6/zOOd+Xs/7ea/rukWv3da3IiMhIWmmJCH5f/8n
epymdnbkkOT//35Upcxpb/yzS3TLUtjKPxcayGxKz/yzRGsPg173PzM0qMgn
l/5zXz7LXbj5z0+88lY822X/5wDuW58znu/7n4XqZdIcjff+z1sSNpyn9u/5
n5VOXWrZt/ef713S8OOV/OfxAIm1P3z/3D470JtL9c+ZlUbpvH0y/7PO2TOH
lwP/2d/q+PSQ3z+/99iT0er1z0pF1Nx5D/5ZjKt+/ZrlP6/0K9e1nvjnREf5
83mT0v/z0YA9+fOj/9yaKEGvOPjPpE18De++/7OlCPWxr/X/vLe5X2Qx8Z/f
D3W7H4n956urnd2uEf/8QvxLOHngP88/rKThfPjP5VIh80cv/LOucuA5N51/
Hj/vl4Pa/8zr7m6pDf/s9u3mD2O5f+acsFRMkvnn/G2z0GHxf+7fY3jWju+f
VTxV0J36n3siTgrWkf2zU66SK9W21P/MUHv4++m1f87okld4+eefT03vCfk2
989dpJKz3FP/fJtb5MzV3/9MI8eflTwohX4bajSUdQaYrMJFOdorhQ8O3elX
ZDfAY5dZr+35IYX+N+PKb1heRFtvaoHill1/X7D8SHUBm38uvXz/VgobSxJr
n5/RQwF6TrfUcikMe3ik77CALtopH7F5UiSFvc/3H9NLOYcMSQ+A8EoK9QtF
5/xzT+NZy635N6FS+HKN8so6UQ3jI4R+xbyQwgT6nTIudVWc/kD49MBfCsmE
VlkON/2Hz/d4pRzxkEJrn44X37NOYtMslX7xTSlc4EywetN2EPlEZE6FWErh
qpPQ0d8Vcmh7XnvvHTMp5KshS7v0RRrfeNly6RpLoTL1GNOnAVGkLX1OJm8g
he8SI5edJrnx8mjeLIOeFEYqhe24vB4mZHG1dE+dlkLWD6bO+edFQeshS0k2
SOHaMgfkE5UgOudgkv8JKfTqvz4YaAUw3qP/3PqIFJ5Se5eqcksN/E5FXJeU
lUIzpjuS+S5n4QM5H0cKtxRut9pG1GkYApfiCRJPNinU8eZ2SdC+DDduXJ02
ZZTCXMV6wwdnjYHyU2K9ALkUFsQK8sjrm4JakLhz9Kwk7j9k0/xF/jqEE9Us
7k9I4o3nNSePH7CEkXmrc4Yjkhg0wpH/6qAVeF/IkuTolkTBqNogLwVraPf5
xLrUIYku89w7c4o2IFY+udX2VRITtR1umxy9CffG6CcLmyTxQ2p9f5OSLdTx
yHW9bJDExQ0+PaXjt4DttE6dA0qiIcWHivILdnDt0e38c293+1MiMbu5bQ/F
ecGxsuWSyHLNW0I15zaQ9hX50hdJYjrN3df1OXcA0o4c9YuXxKc35PnOkzpB
jRxEezhKor5Y9tLNfFfot/rs4uwgiXcJjunSno+AJNHQwO6mJIZfPWkwet4N
VJjs2YzNJLEzsq3cfNkdPszGBCqekcTXlFvV+xI9oblgyXNaTBIlHjzPhhVv
mB3zMBsWlMTxDz/Jn9zwASYRulPdPJL4nlPGtLbLB/SCRTYamSTxXEk9m0rF
U+i4c845bV0CzWe3Hqk+8IPeg69uGrdL4FPLOzoa689g6+YhrfPNEni6pCbL
1/Y5CKW+k9L6KIHMZIzkjd3PwZytY1ixRgKjk7IrNN8GwpPTZnWyVRJ4dfYv
q/++F5DqNZksXiqBYic17D7GvYDhRRJT1hwJrO0eEtF2CwKKfYHKNBm779tz
8FHAdBBIXOfm30mSQIsHjzs/XX0J6nGpa8uxEjj79PTLtvcv4cY3uR/TERJ4
tuFnTCpFMPjRV5UPB0tgqiaT3VXpYMhSVYvofi6Bfz+pnOI6HQyTJcb6H59I
YHJL7lDAy2Cgnx49QHSTwBW9gVLV4mCQlbjLXH5fAs984/Db6ggGh3C/L2n2
Erj83X3fXb4QWHQq1XhsJIHag7pOlzNCYE2R6cSRgxK4vrREKf02FB7YZrzz
lpXArourps97QmE1URnapCUwt2yzYm4jFFaob6nbCUmg1X1K20rlMHBWJvtY
xSuBKj9o3wuYhMGfOzGnqTklUOgYk8AT9zBwzDjQfJFlN34Mm/NoQhgs/mzU
TaWXwB/rXC3aNWFwl8msfY5KAm98njVTbgmDeZWViyfJJNCY0zTpamsY3L7/
4vuzLXHUMWvuc2sPg9lcCeMfq+Kolq0slNARBvYDb3sl/4jjsaVck3ddYTDN
ccHMcU4c5U/yJ/z6EQZT7k8smX6LI3fbGr/QrzCwLeYZuzIojvT8N6+c6g+D
yd8FN7N7xZHE6kes6WAYjOv2Oah9E8eJ1XLepNEwGBE98dAVxbFaxIlbZDYM
rhu2bzS+Fcdi22FDmA+DoWc3PTgrxDGzVD/SfDEMBpaivAtfi2Ow1kHOlJUw
+NXw58XvWHG8fmeWTWwnDK6uP2dRiBTHy29M9FVIw6FXXjzsSYg4nqNoDrlG
Hg5Xrr/h/Boojio6ym1elOHQHXU+WsBfHI9E57KkUYeD0ZdxPltvcdw3xKf3
njYcfpB4JlR4iKOI7LOXw/ThcEmRW4TSVRw5XNZayJnCoetmfqq+szgm5Ycj
L004GCSqSybfEUcphlL9ut3+drT3vpq5JY55N7+N2O6OT/sJ+tf+18SxSoKN
tro4DL4y2lTU6Irjj9iX501351dHZecYwxlxNFvNH6bSCoMWl8hqI43d8TJo
cS44FAZf+utr/yiL418WxlgSmjD4WCT6Ze9ecWT1ezaUWhwKtQY9A5EU4qjt
+NRxjHp3fWVEL3rsiOFL6TXR9/MhIL9iQGG9LoZtPXatiT9DwFKDnUt3WQzZ
ggc8XOtCIDayVfrovBheULsoZ5gbAq2/XxwTnhLD8NXGnoPhIUB19MwZ6t9i
2JF34hmjewgo+9GYzA2IIZUBj97HKyFw73uDw/ceMTzFRMeVoRgCWdLensQu
MXzQuNHzhDkE+u5D6Ks2MSzwnEkxnQiGMzxvy+43iqHMUqscT0IwSBgVbnKX
iCHHjyj157v71ffqU54YCzE0+35Du0DtOVTnqex9cnW3/9qnO21Wn8Hi5s6J
m5fEUKVazkIs7xmYJriaHTsnhuRL+YlJ1M/gaN/tzJ9KYrj35PXuRyn+MG5m
rMDHIobX2z6fV+h8CkKF3Bpk9GI4PlxpnvvoKVwk6bw0SSmGDisZt8XEnkJt
su6jNxuiGMDv+YLltg/EDqrVGY+JYoGVwqdpGm9oPUjWoTIkiooPRH9cL/AC
Ki/i6N5folj1jGms28AL7okdp91oF0UsGKf4lPoEsu795RtuFkWtuk520H4C
fXWlsp8/iuLnjjqxijlP4GC/d6rkvSjqjRUclI/0hDPX9+vF1YhirMpDu0rw
hIGCVcruKlGskSHnqqTzBKfN2rc8ZaKYJ8j3+Hb2Y4iPuCgdmbubt8v/k3Ro
84DZxofbwdGiOPc7+O4tUXcIka3P870nilla++lu1D4EqYcvrjXYi+IUh6qJ
lcFDeNNgyENxUxT5Bg0LLScewJDphPcTM1F0euhheJ31AbjkFh9DY1FMVQ/N
vpZxH+hXH81tG4jiBGvmpsWx+6AQwnTF7awomuQ0p5pbuMDHX99Z3mqK4jOX
wWWzZWcw2ZvyYU1FFFtVlrXMApxh0cXWTemUKHIw08abCjqD7/vDh+4fE0W1
HoE5kyIn4GPZGitTEMW7rw6omKg7Qf7VDwl/9otioqNaxNWfjqCS/fLC4X2i
qMxgklaR4gjfly/T3pMSxYcLToJHexzBTkWMWCi6Ox9dgVHlHE5A8nLKaU5A
FLffprMe0XECGRmPATt2UQzx7aBUrHUCV+O0CmsSUYzmV5g99MoZBIkzVibd
IvhGfalG3vE+XPf3qtMKEsHIl1kUUX9d4WNPmW1SgAh6dRXfVZd5BPL7J9hW
fETQSvBd39LlR7DeqXs97ZEI7sltq9R78wiCxQUptm1E8MOHdTs6dzdYcdHL
u2gpgtmMlN2Vr93gSpP3xVwzEYy4yKxp3ecGUncn0y8biqDlkLhoPcEdqmsq
1EtURZBs+2ynO5kHiLFNzdASRHBE7ZKK7GEP8LMSijQ/LoKtzy0Kuq97gD7D
09+MB0Uwidcl4GiDB1SaVQZZyoqgt7nnyuiyBwiWTB15Ky2yO77ProdLPYZx
I30/GyERPKWYdHLB7zHovH66n8grgiJu2TlJlY+hhKTqOyenCNK+L+HWmXgM
PBenH9uxiOAWhZnKNoMneLwSlnlPL4L2MiMPDCQ9YWRdv5WXWgS14lYua+z3
BG0d3wd3yERQn+yn+W33x2Ci6vhf0IYwHmipZ9aP9wCZyZ8hJxaFkSq2sEbh
rTssBcPQ+LgwDlvF23N3u0HN0VeHIvuFsfagv8D66iO4+NTp2/wXYXzwycyV
eMQVPEWYOTNfCWP/KtOR08UucKbR2epCojBW16+PyLY5A5dDbxlJhDDGBP8O
Y553gry32YZXvITxnEzNQoecI/wwVItmMRHGOwXdxkT925C2lTP+7oIwHtIs
rvQJtgeHdNZjt07vxu8L4DrdcgsoFvt+NhwVxotMx9o7tG3ga7T6Xkd5YaTN
ZDkQ63cDYgl5riKSwlh3cvyFWYMlHHjxkP8RmzAetI/SngQL2Dg8cEuGVhhz
aft0zLZM4UO3RnXnjhAePvaCbvLBFQh58prBe0UIo22UPzj9uQRXZThMDswI
IUn01JOd2xdB+qvr61/DQmjuezNJKlEfFp0Ht551C6F725OK2czzUCOgpaPU
JoRdArFfywv0IOB9fuJooxCKtFp+3bqjC3UpFowf3wkhhWX1Zd/L52DzMYdb
TokQKqyyDzHBGXBQdjW6nSSEyiLv/wgza4J+mTrHmosQ7h9QUfAxPgnPw1a9
euyEcIpmUSNizzGov5u78O6aEAYfSDHKXFWAI/KsX710hLCZZGcrWG4P8L3q
DWCQEkJ/d+afX0pY8MLTl6uz/EKY7bBHs/6tCAZeV7FuYxVCRqHPoqmCe3Bb
OEstcksQzQl+U/nzCnh027jkwZIgmnD/+JEgegzv9DKKXZkQRJ85mQ+B+icx
5w0Gn+wXRLXh6p/vagGHox13hDsF0VXXSk+xWwUF7ks5kH0WRPu3jI15i2po
aPCzdwQFUcRU/3bwa2V8eTjwTGO5IEr3mDSd3ZKAj6yEN9l5gugn8CPeoVML
SOcXZAJTd9tbBJI3KBrC8Zb0KIdoQdwxMrx6k8Qc8p7ROR/2EcSRPBfW+59v
wqjNu2FOV0G81kU1vVXoAEKad/RXbwvi6OMz11U+34Vg8q79764I4uXsnpOf
Cl3A2S1lUuPwbrwt94dTJz3hlIOSufGQANbpJcqv0TyD1vM3zjwBATTY4Y7Y
5xwJH2QOkrEeFcCDfLd77bcj4e3ORmWynAAKKH4QL/SLgszXIbuXOwEctXUq
VoiNBjfad6Q7f/jRuav160liLNwb8KsMmuLH8wvS3J5acWBdoX9baIgfTzE8
Nq1riwN9q/Hek1/5kVJFbkZjJB5katkr3bL5Mfx1AJ0ebRIIRfc5MCTz492P
A+dDQ5OA/Xa2ZHwkP+YvxJ2dZEwGGg2n3n0v+FH8sead5/rJsC1ACHvrzY+h
DEthslHJ8GeJ9vQZV36kCvyYJfw9GSaaOkh67vCjOfd/mYNlydCXklRha82P
b1Kr0tLCk6Hjga3Dugk/cskfSrFyTIYmXUXJZxf50aUqJ1F6Nz5RiqSX9ww/
tqiJx08cSIacjgjtY0r8yFghI9/TmwRPDtRKPGDjx0XWGd6p9kRwoQ7soabl
x0yZPo7bwYlg12cYGrXDt7vevjIvnUuEy4HTO+XTfPjDrohy42MCyI9z9yx/
4MPyRucFOmI8SBCHQp6+48O7fTemg9zigTfytRZnKR9KLF8aYzseDxRqquUK
KXwYJnb8F19ZHHQnOYQ4PuJDR/ftj3tzY6HV5ZgW+T0+lAufe59vEwsN5yh2
Qm34cCBnoOaQVCy8kfhaJmbGh5G1bbsXlxgo3IixKzbgw/M/6kqPp8RAZrul
uMrZ3ec/jYoXlsdAXPb+7jYVPmyXHhWc+RIDvpcaNBf28yHH6DCr9VoMuMkH
b3tK8WGZ6jB9OnMs3KW8UsYiyIeXUocoByVj4WrJvNgBOj6MMRtcM9aPhUPM
AtsOI7wo1N03qLPb30jBw98Od/NiGzz7dfN9LKztO521+pUX/V4p/vTpiQXU
eqDv9ZYX55wCW9/Qx4GuV0dmRCgvEpmPE6Xs4qA4eOqRsT8vuriMvvnPJw64
ksjOC3vwovSv4PKr8XHQ9/bARtZNXnyZM/Y69EscqDRpttqb8aIaa1j265E4
yPhhmnHIgBf/3j+V8XEzDuyWA3WrgRct1CPit2XjoZU8XeLJEV5kz4NoXvV4
OMz2dk1dlhcZUO50lVU8RIq0t9CJ8SLv4nPHe6HxsC4/kfaVmxfFJSbj9+3O
p+lJkofhjLwoZ6j1YWQqHmrPcOsYkfPiCf/MuQSeBPC3UV8dmuVB3elrKiz3
EiCjZyuZNoEHM90uvvj5PRFoJjlcWkJ5sDS/uDyUIgnsVvedCfPnwXcDLANn
DiaBAqfxsoATD35Vaz5U8ywJosXvNg3e5MEfLvuu3q9Igo2D/kmZZjw4mhXw
9MBIEtTplGsfOMODHlTidm1CySBl0iy0DDzoJdxOIq+YDAG3RpaqjvBggNKT
iGdnkmH64cZHD1kefHn+wN4xi2TQ82dLVBXjwQjb/hrVB8lQGrnHkYaHB2O8
gy4kByUDdwZoNTPyYFL8yfHN9GRwLbkkGEq+295cwD40enf/1TosGq5xY+l9
yRGbJ8mg0vq0kX+OGzuC5K8QbJMhoy8+fmCEGxcyjrZzXkgGu40mDdtWbpTv
0MY6iWTYUJCOc3/FjRctepWOBScBT1a/aroBN8I3++xkwUQQqm/aij7DjYcj
uO9JzSeAxEBZxYv/uFHqUu3xvLoEOMATuO++HDfS93B8qbBOAE1/JfbTlNz4
bejtbEtRPLjYhAzNlXKhPI+LhohtHLj5uMWP5HDhTWYHB1K1OPBKtjb8mcyF
GVQ3IocE4yDox8mmukAuFFg2+J3eFguZWpNFEZZcSNuu6LvnWCx0yqg+OcHB
hQPPlxsO0cRAt5r8iQN0XMjjPTPDPhQNA+a8yxIkXKjvOsq5/DYapqPnbJim
OfGjTadVxe1oWCztFiMf4kQK8+YX0drRsNra0Pv3OycSLjWUPRSPhu3pwsjp
Zk5003n3y3grCiho4/UG33NiuXoZpfL3KKCT9KPrquJEOe1UCv3d51n+u1ff
VMCJx/WbJpeWooHLxMSDmMGJ6lcWWyMSY0DgoZZSaRwn6lnyVRzVjgWFiT/C
fo6c6On/Rpe4FA/HKQd+PrrJic9DhhQtEhMARD+H3THjxOhYOgFy7UQ4a5RC
bXyGEwvyjMc1EpPg+qdzc7LinPjAxh2sJVLAZlQpW5SXEyntQw5XPkwBB1KJ
61zMnPjiboYU7dcUcD223rWzzoFRrs0MuQ9TITQ7s6a9nQMLAgV/zrakwcWa
bHT6xIHnkGX8sVA6cH3Lq+NCDpxaIv/L4pAOP8YK6ivLOTBA+i9l6rt0iNss
/mD8mgOFr0xyHGbMAFPW8o+baRxYFvRLvP5qBohIVTUlxnKg4KscX/uKDBg6
Xv0FQjgwdIHh22/dTEjXJbYM+XEgjfJtYfPxTLhhWdfq48GB7r5tt356voI9
DxvapZw5cKj1cKU+bxbkp37usr3Gge+sVvXUT2dDxfL3vkVlDvR6+cRrb14u
uNL1DIQf5sCxn8MtaWp5oCzcN3RkLwdeldDgF/yVt3vejfx25eZAjSq6Umam
fGiKm58lW2RH94VP99aPFkJg4dJ8+jg7VnBt8Bz1KQSdhpVFjX52XDohS3Rs
LYRvs5srz76w79ZvQfQzNkXQ+x/tDtsrdlThv5DeH1MMiYYMpGUJu/HA57TA
72KwuMVMfimcHQutyuYvHyoBCU82yrVn7Dj87HdkhEcJ/A7npI57wo6Chdwn
25tKICubh/bkA3Y06tQaZuIuBdsafvp+B3a8kTP95+DhUpD9JsT4xIodG52O
BQVolsLcmCiz+FV2lCT4Sg9eKYXiTQnWBn129KL9hkp3SsGJVYbdWpsdh74J
Gwf7lMLqcXnu3CPsmGZTFUh4XQpkQceFD7Oyo3Do9ELCdinUp54U7aRmR8+r
x54vs5WBbwWI399mw34pX4lz0mVAP6gh/XaKDRPfCl/a1C0DDoUL+1U/sKHQ
qH61cUoZdGkZHhytZkOPgiSD4rIyiDExOuxXwoa/Hk7P0jaVgZCf2dHPyWwY
y+QrWrVYBtLdtwgXHrEhr1KVL7daOVyWemu29y4b2tj4lEsbl4P/PTpPUms2
rIzW/X30TjlU1hgl/7jKhjSf+Li0fMthnC4bC/TZ0HhtVP1yfDnwXFod8NVi
Q5Nrc+V982WgnaZJanqKDVfmKa7oHS6DnBO//2PYy4ZcTAfSDleWQLef4rVh
YTaMj1fXzFgvBvoOH683nGyosO/KFNfJYrCzE39vQ8qGF7V8D6/XFoJ8grla
ww9WlKas8nsSlQumEwWW8S2sqCA28tVSIAdeKu74ONaz4g0CM49WShbMNyc0
iBayYsZDqyym3Awo2u7W8PBjxc0e8dCj7kkwdHqvtYE7K5pzuYpt+8YDW9RD
P1lHVnTVay2uD46Be/t5PnabsmJevdu38+lhkPrIerzoIit+2G6/zpMfDN8a
y2kCTrPixXv64+9uPwcKDqo95sCKwmPbNGe8n4KCuYH20SO77y83vFAo9xgs
89JvMsmyoqZapqbpMWeIWF0KGBVlRWNeiaV5b2/4oKaaU83Nii9YnOsO3QqE
leDQpjBGVtSnGNY6nh0KUr8GJ23JWdHI8/mBHoEoMNxzkF5ljQWtNg/zuAXH
QUXt1zPzIyxI6Vq31HciFcYZhe0au1mwtSLSusA5HXiMHQITW1mwa9m253FB
JjxYZGg+W82CA3fZ34uI58Jx8dM62WEsWKdwSPsgZTHY3o5x8AxgwfoLAt68
YiUQ93Y86NJjFsy5R/WOlFAKmxf8vlLeYsHwwu4D7Q/KQTa5a67XnAUjvtbb
vomogKvTkiylhiwYPZufnlpcCe986vV2KwH0lPXmcZx+A7Ot7HePKbGgzRl7
/Su01SAkeD2ERZ4FbW0vBapKvQOdm8VFY+IsaBfw34d9qjXgUUbaXsPLgney
9pFwmBMhn/T8YgQzC1YN0wl6nELoP5fMZk/JgguRfKbt0QjCCv1W/ivMqDTA
kJd/CMHzi3vmkSFm3NTfXuKMIYKx4dn7u6Urjjh2ldl218DhAT4trWpmfB+e
f5/IXwO/lyomjCKZMaDLdMMmsRrO8S/Kup9mRlFjIll1WhVIZeDW06PMuOMa
Vc8yWgkk+1+2BEkw41jcbV8rqUooVZG9k7zNhBW/hOmYs8uBz/ZGyftiJtSw
8GS9ll+yW88r+HxJYkI5r8vfyueKodmd3KAzkAkZ0g5E0B8sBs+Q5JXfN5iw
fWSAp6ykEIz5bzfOXWBCJ/p7ZeyTBaCQcTJ6FZhQ7ZzYsXd2+TBW2X2Mjn+3
fasXsNzLBVTJpmOnYcIgtsP1VX+yIfbL/R7+ZUY0vTisef1+FjgZauRJDDHi
/siwzwzrmaAzwOEu95URSX6o6pa7ZYCM7fC5I9WMKNMk+72OKwNI/xQJEbIZ
MSw4d4/b3wzodvec04xkxG3DvW6KPzKhlFoX9bwZ8apgdsts5SsIChEMMbrD
iG+GpUWzYrLAmn/a4poJI0rkZDpauGYD3/4AcqejjDhzJJ2n40QuPDF8rxE1
xYCOeknXqmrywXgghDv5BwP+4hIqu5e0+7225uNZDQyY6NBCzgWFMO6+FfAm
iQG3FM4ZfhUpAp2Mo82/LjBgQ+3pjf+GikFmP1Xib2DAPxLP0hP1S4CsqsNh
To4BJf0/nduoK4GyL3dZSGkY0F9HK7k0rRSCDGGQdpke3xX5abNylIHNAFMx
2xA9TnM0Ltr7lIHAn9wLEtX0eLpHXV3aqhxW3F0l5LLp0fXU01nvznL4Sq29
rBhJjwUp9VED6hWQFcL94ZQ3PY5QUPx3srwCvPh/R2reocfG0QQJeLq7PzJK
rfVM6DH+5EBwVWs5HNnvrWR0mh7vRohtHRIohwkVkW5bCXrkUH/VKVlSCokn
/MZ+MdNjvzvbsv9cCVxUmP2jt0GHZeXuHNP7SgAlq5mPttOhhczFC8XpRRBN
c1md3JMOnbnJUt1lX4MOKVHf5RYdigb8GFDNyAWKdUnzCUM6bFnPF6YTyoHb
U0uuX+XokKPvSkIk8yvQaHlZGNdDi8XpZdEFi8mw+WHlHVMjLdbI8hRMdSVC
MdHk85NiWuRmaZB98DEOBIv2/bYOoMXULiGZkNwIaM8OXep2osWOqs8Zggmh
4J+6RqJjTovGxxc0Lgq+hFNx5ky1Z2iRgfSMij1JAPwJa+RXOLobvzH9pO+w
F+QEyu95JUaLy3pH/7x97gPmTyOO8DHRosbrp35tKwHA5bGp+mKNBsNoO/jG
LF7CZ5fr50lGaZDjMOsljxuh4Hm7ydSxlQavWlSK0ApEwBGbg3a/39Lg05dm
k2FtUZBqtOP7JZQGv029ds9WToDL+jfCCR40uMVroKmwmARMZ5pTim/SILuE
h3//hRR4qBxXHQ00eFZfUt+WNh30hI8uWs5SY16+4+hjp2yg4knc+f6TGonq
mqfVv+ZANQsl45kGauzs5Sug25sH0mTfpA/GUyMjfd2DiP582B5xMNnSpsaK
wzS+BpLFUPqr0/aOIjX+t3NXqbJ99/ztUn4wLEKNn5t6Jvk8S0Dka9pTQwZq
vBipnuAuVwqdjXRhn/5S4ZhFge5Adyk8w7vJysNUeNHv8hCLfRlA1Y/XBS1U
WL0yoG5xrQxWighvxd5QoeiNm9lFl8ogLyfzY0QGFfp1LjCQnSuDa2mMXTQh
VDij5npHX6UMeOKdhh+57cYrJetIPVoGPi9Uti0uUKFIGHucqlgZXLnTd1ie
lAqr9BT3PP1bCoY+M3SkM5R48XLSFeqZUtCL3hz89oMSp81ogwKGSkEN+YIf
FlKi6O1fiy+aS0Ge5dJ0vRklBrx4Wh25u/9kJG68jzpDieIRc7M8MaUgruQc
a3uUEt/FG4nGBZUCj3mYFgszJS7myvkmPSwFNqdU4eF1CgwsiaoUvV0KjP5F
K2W/KVDyLelUmmUp0MRjs387BWLdLUEp4916q/Br+pUaCjRu6tTN0t3NX+/7
HsnnUOCfNoLXXvVS+Pt9Rp80kgID8o/3NSqUwuLU5p6OJxQ48+uXUSxXKUzv
0JO+sqdAPcYnnXaru/UjO/+Ph0YUWKwscZ7QXQID0nsLzqpTIKdd4xfW6hLo
1NUyWRSgwJ4mxvryx7v15RxyBjwixyV16ZTTQcWQSt46fcWaHA2cmwQE7xRD
Anf/e/kL5FiVbh89d74YQglbdzv2kqMHeenLcM5iCLzAoP2KixyHDl2it/5b
BL7W/CKuZOSoem3d99jPInj08liLcDcZ0tWCR198EbikaWUsNpChw/zwWqFH
EdypuOTWUESGbcJ+zt7mRWD7+caF6AQyPKS7d8FApQgs+5333gogw0iPZjsZ
iSIwXfIhPeVMhquv74yvUxaBEXX4DxYLMrz6i92yeawQLvKnFQyfJcMIYnnt
TNduvb6/2LdciQxn/X6UKc4VgpZqrUmABBlq6q1nu1MXgeqlVoWrLGSYxCOQ
2CBcBEqPZ4dIx0nxfLaZr4FeEYh+3Us4F0WKZU2/DZ2LiuCPQ9rf4WUSfMS4
/LfUsRisnYZMWYdJkGf7Q/eVZ8XQ/VD0w6lWEiyfjX5HllIMtT5J4TE5JDj/
Vdlbr6UYFJ71bXyIIsEXyHTj70gxZL0UvP7HhwT3Fg1oJW4UQ1Bs3EFdMxK0
CvVhmpYpAdLk7mi3syRI5m24EEooAacMXpKcYySY7CjTccywBMZzLlt/lyJB
Zcv18gG7ErhSGNVCwUGCPy9+ifHzLoGWsi7FQ6Qk2Mt/em75QQl4u+w7Tuzd
Ie7r15DulCsB2yfZNHwVO8QHqapmpbv5SS9Q+rtTyA6xwQqiwqKK4UhURsbX
WztE1j0nv947WwyCqeJOezV2iObTx6j1SYthskKYdWB9m7hx7/B9lltF4DPM
feaM5TZxv6fUpRenCqHyOA1R9OgW8YgirzzT51xIVPd94ca6RWx5y5b4RSEX
nupRXP0+tUm8pcLA9DwxB/RvkKwFJm8Sk3V3ZqjvZYPSXY/GMddNolLXqsmH
niwQdtuMUDHYJHZdXWz2Uc+C6ZDVQ6u0m0R629HXpLyvoD3ehfTCyAYxY75P
EL0yoerVn6+vazaI/93/8cJjOgOSi+8l0sRsEPu22raUDTPA99283XXHDWJb
m3AA6+d0sPvocKLm3AZx8Putq6/t0+HCt2laXpkN4vKviv3aLOkgMjGe2dK7
TuSb1O3yNEgDqj83nPdUrBP3zcdlC6ymwsz2iKpPyDoRVsbcKmNT4Q3H4MAx
jXWiJZmn+Hx/ClyEn3wZBWvEfL7CO2a2yeAX9Snokdcq0TNShMttOxakPB90
5RitErMJ0cfXDsdCvY20UPeBVWLzOLOZy80YIDv+9LVS/1+i8ImdV3c6oiBZ
TGHZuvwvUW3E5csMbRQQ6IeVo178Jdq/mJ23JUSCWy98WT7xl/i2/9cxq5xw
4G+YY5dk/0sc9jcwHeoPgzevE69cnFwhMh5q9jLjDIPLkefSvGtXiLT7Rahe
7w+Gvx6bk8XRK0TjVZWGPwLPINI699DQ7RXiq1B1j5fd3qCgZ+zKqrlCjDAK
drh44BHcE62kvbO8TFTg6gpfV/wPWOmszyd9WSbqsqILV7UlFi5yxrSkLRM5
Oow/dQ654Mx7Rxk5/WVijGLWPcUVP7S5sV9zMv8P8cxvHc3Rx9FompvpY2mz
REw7lXrwXX0WboUZfAkjLBHZoi28DQ1yMN6NguM91xLRY0Gkc240F49blVxZ
nFkkLp0ekApweY0/z11LE21YJF5LT3ogRl2AD46wTunFLxLHPfaP/jpTiFzC
eMjTcZGYto85aC9HEZZR33YtOL1I9NLovjm9VIQX5wXr+kQXid+intI4qpTg
0o8vtExrC0TJiQOv1oJLMbT20fmTrQtEl+O96p79Zfg19OdAnMcCUfTXIe8X
bpV4WjFi53LPPPG22MDWq/J3OCGopulfPE+sdXweL09JRH+qpaCKgHliwLW4
hO49iA3f9YS4lOaJfSUayhYPatHzD49r5p454lPDhIcOnPVY0+BCS/ljlrh/
Z/FjUXU9bkR2RF/3nSX+fKXFs3y9AV2OBVeIDM8Q5daWyh8Vf8BbHtTLsfHT
xK6403/8dT9hIKdkw3W9aSJFp3JzydQnzMtTiZClmCYeYdqf2efXhM2q5lbL
FVNEK03Rx7QSn3G2212x5tYUcSHA8Zndr8/IfC+O0k94ilg+mWmiyNSMB2ir
OnW/TRIfnek+sH2qBe8d/eM8eHyS6O2bPKWp14phLawaObMTxNEjxw3yxlux
1Go/l2PqBFFz7FsNy5M2XAmzLaegmyAyalGH/ixqx4fzAyvSP8eIMbQOSrfY
OzHWb/vDfOAYcbOKOvVrbie+ERKIqoIxooVtCr2CWhdunr109EzWb6LU546+
Tafv6JndfN/h/iixKPDE06Cun5gCU5pH940S2U52zi46dGPddxoekv4RovOM
w+VL1D1IQaVWGaIxQjyuk7pPRKkXJRMs/I3Xh4kJ2ycifFp7UUPh8WXx/GHi
Tn7n9rjNL/S79ma1lHOY2MhM214U04dZa98/un8aIu7B1BNch/rxY/BytIb7
EPHFHeUM16Z+nJBmv8l8cIg4Fkx1uufCANIRDxz7MTJI3PZU2xSbH8B9hjq0
KdGDRPa7Xvm2gYN4dubWz5tnB4l7LNC8WGYIg/hfPdwoGSBe/E+54T+LYSwo
rteutx4g3jroej9gaxhbtYd4X/APED1FK/e0xYyg+Pxxp5NU/cTXpApB5u2j
uDdbetk7vZdIVbfv0mPVcbT7fi2Ck6GXyF98k7axfxzzKRMVXzn2EA+lvnrL
5DaBB69xuDSpdBNNvcRFEksnUYmf9C/L4HdipQr/RLXkDLpqK0elan0nNh8y
iqOom8Hq+/ePKhR2EYfFos6dNZ1FQsfM/UuPO4ks5OzFPZFz6Em2h2fidwcx
/KFMcQjvPNYdsKx01ekgJt7VfF5mPY8UZkmXGcq/EbNv3rD8WTaPGi+6VxME
vxFLLZ6e3CJfwLwIw99y/u3Exgu6qtF6C8ie2P7t3WIbcc1H9y1H3QK6ZurU
njNpI+o8TbJ+pbWIA/lN+b8aW4npT+fYj39dxJVk50/HR74Qo4xMDi4n/8H6
Ivt046rPxJLr7qofjyxjSJ3VY9egJuJX+wSDuC/LKDtioPhG6SNxVJOCj2t2
BS2kj6Qce/Ge6HXKL+zmxCrKK8m7GV2rI1LpzbMTjddwU0vq0sOjtcSP+7PE
Ln5dw09GQodiGJBY7aoyyKu2jrKdsnz33r8j5jf0JvVXrmPjikBrkf1bon/a
Y4E0+3W8xsPou8BdRaw2FT569s8aRl2ZWbhjV0qMkjEKE2Vaw0Puv14VchUT
Ic3t7NNXqxh9Mu7twJMC4pRgMuUErGKL5yNp864c4seXTqeUvFeQvP5KWP/e
TKJY8NGVnfAlVKJW3jF9nEJ8FLz+ujFzHnvClbqqX8UQTa17vs9RTaIG/9rr
t3rBxJAVKSONFwOIvr+5Jn8/IrJzN9101PmCskR5DsXJawT1y8xiVgdKYCX1
WcRLGz9CdsW0uh3xO7B8pljSoo8gpHzJSTEzGoVMNTX+09cSCW/ft334aDML
Qv21KRm9aYSngxm58WKLMB43qCjln0XYNzi01MGyDCWXSZsyD78mtA8IKzPu
/AXiM5aCsLlCAm3Rz8e9TqsgcVNk7YJ6CcH0SPmHBz2r8EzzgCpHXBmh7G0o
I6fKGsxJwouO+QpCDseslmXlGlyk0PservGG4GKkI//24Dq8GTITNYivJqgm
5rOz5a6DcO3tW5yLNQT7Z0tR6L4OdV+J3ZNpSHhEMPq5TrMOyzWNO00ttQS/
PzV8iuFrIJP/VSJ3rY7AVPOozER0DYwTvms/l6gnPM9jU6QpXIUXgf0Ot3Qb
CFRx2aXFp1YBH42FnXH9QKiZa1C5W/sXFm3nKvdlNBJGwq+3Lxz+C5LGf3/R
t34k0B4nvXYvcwWeKVHLfJH8TND8umaBnsvwTpr5XJ7eF8Lehp+6M5N/YJ6L
+17go2YCY51o3abxHxCjEo6yy2wheJpcctGaWQKDZanqs21fCQurgXvDHy+B
2r2oO1werYRf0Biaw7oEv27ye/uytxH84f5v7ueL4GKRFLGS1UY4BNLHn1Is
AouReNaNU+2EX4TvLxbdF+Djz+FOn4F2Am+9r6au+AI8MU4nT/P6RpCwT1ZS
aJ2HxSuSZv0fOggnid/4tvfMwzdTaa7Lut8Jgx+casTfzkCk5T4fbfNewutL
zK33AsZBZ3SqyJrsF+Fjtkqx8vAYUN3I63+a/oswsOEcTnVyDO5by52oG+8j
sCf3Xo6dH4UD4zM2A8/6CfsWmU9Ynh4FJ4VLZx5zDRDU1VQF5dNHQLaM+2FC
5QDh4Xj2YO2lYfBf/vsyj3aQEHb81/vnRUMwovgj863xIOF1IEumAf0QgEvl
u6bcQcKHPlV/YatBSCiP7vi5OUjoO3DfdqJmANZXHkyNnxsirHjlnC3hGQDD
o0Zkq4lDBK3h9mMCfv1QfP8YL/X8EEEh9Gw6zWgfMFXyHuD6b5gg8t8HpmWV
PrBdXdOQDB0mMM0THg4m/4IPSt0mCsPDhPXEquHm7V7wqIp9rv90hPBtM68y
800PzL3hHwv6NErwrgy7dvnQTzi7sbGVwPebYG/N2KwW/AOyTvRyvL71m3CF
y+/owdnvYF4dr/KZcYxw2MmVgS6nC1rfCSZS648Tfu+3KH8r2gGy21ulXGnj
hPa+bpGsx9/A/1TfZ8k/44SaFxefh/9qh/+IiWuqkROEiClNc/uYNkjc8WC5
MDZB8IytbTL62wrrBDPpa0qThFvaJxQ1DFqhGIUNPLsnCc90TDhOibUAEymJ
3ct9U4TreXoTXH1fwPa/Aa9EtymCj5vBUlLCZ/jwBGNeN08RCvdI+/j/bQKx
uuTCaqFpQm/XKsfd803gQfak8fPtaQK1T1OGUe4n6Fax6OvGacKhg/FHVCg/
Qdh7UYa1azME/+eEy+xVjcBwSmUoIX+GIJ/lZMYh1Ag+ldcqVdZnCJ312Tc4
vT+Ay+s0y8CQWYLwNrsz97kGuBwl+U4Y5whv7oyH8A/WwThNmvYQ5TzB/IVg
jIBGHUhqlJSW7psnkOfoJwvm1kIifntmdH+eoDNcnS/sjPCyjONoGtMCoTVX
LEqTrRq+LEqkOissELhHAzO0e98AzQFFRi3jBcINob8lZzKrwDPHYHg6Y4FA
HvK5Vfd4BVSPWekQPy8Qzjcp9p+nKIc1CZeqkIUFQjR58syFllI4cs1XwpJ7
kdCvTLdpEF0C95IiXx45uUiQc3Giu3ytGPJ7M9dpri8SgkXNQ1Yti2CKt8Kq
x3+RUPI6xvNqbiFYhn8/4dmxSFh9zlujdScfDAJ4xj/dWiJIzowetUnLgtAP
MhfiQ5YItyiNuA4ceAUt5MdqHCqWCBWCX/78rc4ADQ+jcHbyPwR9neIiv65U
UHKMBZPYP4RZ0y1DHfZYcCrMyd2PfwjPXg6z98lGQdHMGy6y338InPip1V4j
HPba9ExnHlwm3FaRjeBbCgR+E4GY+Y/LhGiB+wsW2zfgcqwsxfvZZcLlkIJc
9uB1Yvh35dsRHCuEI92SkGpjjQz6JurHzVcIYqV34iccn6DWS/tCBt8VwqR+
ulgT+OHTL+78/bkrhK6F71k5jC+wjjbIt6hthXD949ModfsQ3NZIXPD+u0II
bu1Y+2wYhsd98q8aCv4lsL9ZL3diC8f7tTWNMqp/CXFpIk6CLeE4f6I/oTno
LyHTxW7WQT0SBw+InOaWWiW8EaAYZlWIwRvBOyc4dVcJKjdPiznvxODkXJ8c
+4NVQn1ZsMWPplhczE9gZW5aJXzQERxItIhHMjn+n1QOa4RmzyO/9r1MQt/A
9SaKqDXCqxungjTlkpFu+mc1Ga4RxgMH5dZCk5E9Nzp5m3WdYPdT2v6qaQpK
7OGyWSlbJ9y7VzoltpyKWf7LRn/61gkl0Zeedxinoex4x5lF6g3CUs36Xl9M
w0LN0pPzBzYIiqMJTUel0lHhVdj+WaMNgjP9f7YTz9OxkspRdNprg3Cm6t3E
xlY6Kt/QZ5/M3SCwet/njEnJQHVJ1r+jWxuEDK4pu2/jmdjkMz8+LLVJuD2Q
Hn0n8BXqjHztHtTdJBzJMa1nPJCFhmlBNb9SNgkNJ9v5te5n4w1RBt9O9S3C
gGVVUzVlHk56Tt3/5rBFyJZ3XDHKzkOHgaabbVFbhDursmJ/z75Gl6SAc80T
W4Tt58kPD4bno58gNWfDi21C8E9Vjgshhcjg/pvqffk2wblGJEd0uBCDextW
sX+bYJy2RZhXKMLoOJ/e6oM7BGH7ilsvfhShwMb1ljfGOwQK/XDSq3uLMdlY
FSu9dwhjR+5G7XUrRok3YsXleTuEz/w6cmvNxZjFS5Ze2rlDKCDZ975RuARl
Hw5GFG/vEEJHqY0i75bgZbsB1R/sJOBXpqbwzrsE+QRalS/JkECAtP5ntlsl
2PuZqNilTALPYsyu2eiXYOKjAnmD87umt197d6wEzfclSXdYkUCgu2swu2gJ
jjzz4G0PIoGXFhE1NXPFmHncge18GgkEf0s14PhejDaTJvStFbtWL5y+WVOM
09ont5oHSCB8z2c+zhfFuEyzOfDpEClEL5Dcv7X7ffcSVFkmNEihVtX0ogJb
Mc4ffE6guUIKkxFvD26tFeGkMV+ihjcpKB93mQ76VIQ2cxbNVlGkYBX47dOl
oiIc9c7e9MklhZd9B16JxBThQJ7SlfffSGHAa8qi2KYITf/zfDY0Rgo0nVqE
R3pF2NPZWEW6SQqHpDMF1JSK8LIt84QICxlceUi+ziBShF3bhjwgQQben82/
d1Lv3n/DEjXNlMggT7CmNHGuENukf7t4nCWDztv8odbfC1GnWi4zwZwMZFau
cIuOFeIXPefOaqfdeGeTSo9RFuHp0WqKXn8yCEodOq8vXoSNDykUNuLJgLgq
OWf7XxHWpYaFHmsgA4mMvD3x7kVY+adIt5yTHGa2qhK3l4swLXr2U0E5OTwu
6PB//7sYReSOrLU0kUMhJY9UL8XufNW6y8z2kcPwlSvv/4iVYPQkg58sNQWo
0wxtSZqVINfji+Vn+Cnggalk3Cn3Egxjjx+13U8B2aU2SpfiSjBIeZ9a9iUK
YLaYu+f/vQQZ2u45frxFAf9VHGJJWS7BAKs3qWOPKcCJ0eV1FXspUq2TtlOF
U0DG9arT7QdL0SdIm1QqiwJoBWT7+NVKkUw85IB6NQX0b9hbh+qU4uOKH2aW
rRRQ2lMwT2NUiq4D1jWpqxRgGneYfMW+FFedCmZqGSjh8CPnQLsHpehM+1dg
UIQSaK9UcA57leKdQ76PhLUpoYTvhExrVCne8ElXio+hBOpYlXM5taV4Xmaw
OOI4FRREn48y5SxDw8+r543PUUGqvtOUiHAZGjuwLAiZU0E4QxRhWKYMLUsJ
+7OeUsEDz94xa+UytLl0qcUumgpsj++c2KtRhvbrDvYHc6ng6h/R4Gndst31
6cuw8o4KdPLVRvKNytAFEnPftFIBwcZa6e71MnQdLjv9eJgKDoo9CzxsX4aP
fZsnVFd241PoDVeM7Y7PnlF/GlpqGGQeVDLTLkX/L5vSzfzUoM9/L4gqtwQD
b3M0hshTQ70U+chrhhIMYZO9YfgfNcgdCj9m4FCMsZevpPffoAYa7fKRtINF
mLhxTzX9ITW4XtQ8fiasENMSnw3ZBFJDA7FQeHSiAPNGqkQWi6jh4cwfZ687
r5F4hyeBYosa2rXcxd65Z+GIf2fEnjAakNnh/2t5IAHH980ozGbQwIys+Tmy
ilicaaHoKK6kgRDj9LSkk9G4wqHAptxHAxtlsnrdp8NwveJMEckCDdQO38m8
3xaMO8bX9RrIaeGTmO3ansJApNh2nQvgogU2j8UMy3d+SJMSGqSzhxYSf7rq
JzV5IYNajhy7Mi0okfpIPzT0Qdax2i/fdWjhRa3m/bSpAOR69vNWvAUtDD2h
b2z2fIl8cgt0Fk60sJRsTb/5LBSFW2lyJP1oQfCJ9uEQgwiUcBTRnoyhBS2L
fVekRKJRrkrXz5FICy6iczl6ZQl46Kq1lFI7LbR62kwsEpLxyI5Hw+YILbTd
MXI/mZGCoJ5P7ktHB5t6JzLab6ejfhuDZ7QBHciw7vzZpMhBQycJYRMbOthL
Ou+vaZuLV7iVa0Qf0QEsDAiEfs1DS5NbmznJdGDYVqcmHVuALhOfHrybooM6
LdWmrKhiTBTU7uEmoYe6fjYLS6MSbDz/8eQ9Dnr44jL0V5i/FOefaiZ/kaGH
TsbioJ7eUuR584FM+uSu059IRiWWIcyqWz05v+sT56v1zcvxplhDY48VPfS3
i1xgEqvAUEO1vUdc6WHi5vzEp+EKfPPsfWBwED3s/52uI+dYgUM1KnOTqfRg
sh0P0evlSL9Ue169gh58uSIOkXuV49UryLXeTw89Gk+5fgaX4fa+d8nXDzJA
wCurPp/8EpQyVyavUWeACqJJ65xiCeqGv7XiMWaA/u8GdcbvdvPTZtXeZi8G
OEKjkXngSxGqNJeXHulggPq2s72V3wrQlvQId8g4A4iQncSgoXwMUyx7OLXJ
APcOyaVbLrzGkYSSU8mSjCAcxmjHwpyHTx0KP9LeZwSfiy0b1qezMD91/z7L
54zQ4VPTf9LoFXZ15b+oSWKEA2X579ltMnGbTn6ep5QRAkcTX024ZKAU4bW+
40dGYMiuWru3nI46jrJlzb2M0MUWYFcQlIn3X+VyyywwQqnb5f4p6Sz8yJzT
28vLBFf0/tZbGufivOoewlF5JpB48+FoylIe8j7ISglRYQIyicicX4H5aDuQ
eUPjFhN0PuLTyu3fze9FafMF1UzgrlslJ8FWigqjYhfo2pjgGoWodV1fKZrw
ppZZjjKBXpV/snluGeY/SXblZWIGCfHL7AlqFainn0DhbcYMX1ZWVjnvv8EH
vgLWvxyZoSjX9FCp6ltMeRP36ag/MySZf7ilz1KNi2KxQdOFzODUFNEXlP0O
+S7xLmg0MIPR4y1uOZcaVH0efSHlJzNoKVid/6xCRDsid/nGDDPEOd1oKKVG
bHp4hmEfKQsw6fINkTsgivE7LufQsYDmux9b97wRXatj+/ZysMDPx/YbjzeJ
2G5S15gjyAKST/ds8C7W4J6dicK90iwQ+Gx0vWTsHXoms8bmHGCBlZcp6+d+
VeOP/4557z3OAv9FmKyPtb9Ffx9/g73nduMlda7x11Sh5suGxD0xLGDUVPZ3
ynb3PDsw45edygL49e7fp+aFuNLGcXdPHgvw5FyZJeTlYyaHpeoeIgtoE8/M
vPLKQapYsnGZURbImdg3dX9/Kpqe2NuWNcsCqn/SJnP5krG89/wbmVUWKNs6
2kGkj0cmj4dpWaSsMMYm8016KAqthFMCZehZwU+Gp/1lZRhW40fnLA5W0HNj
SZ0TfIkc1+ZNZYRYoSRzPonY5It25DxaWdKsIJmackJ4zg3r0wkHZQ6yAq28
q80DzYcooGHNl3WcFYKPR5nkynmh01gQuYwaK/CYLvjGVvih+N6+TulLrKDV
3m5XbRWCXgVurtLurMDM60eMfRiP/72pei/1iRUcv90vztXLxugrg3mZ7axw
8GqLm3RJDs5v0kRK9bIC26iEZhpnHiYTLttIzbFC10pbT2x3Pq4OPD6fucoK
vQSOC8+FClHP69VxKTI2OOdDUPr0tRBJGlYYpDjZYD9jJJn6oWK8fENoJUOI
DVLO1455DRdjIbVGv6QMG7BFTn/B8BKkybb/mHGQDby6uYu31UvR7HREkeQJ
NhgUUo1S/luK5VPVsRlqbADXHdxcs8qQ6cWIt6QOGxS+irGoNCrHG/IM9hmX
2GA7o3ozNKQca74eNpS0YAPLrbQ3C27l6HTdz9jZlg1q4+pdK23Kcc/fbtMG
RzYQPT56fLdUwXAhL2srXzbodZF6yyBbjrb2rQ8yctnAYD7rrXRPGYqQSrgv
l7JBcdCnR7MfyrAz/P4T9Ro2YJWbPFFWvHteVgs9G21lg+ab+6pVnpUhF719
nOQKG2gMv642OVaGXxIxyXmH7f+auNNwKpyoAeCWSgkVbtbLve4SKRRKIkdK
lghlX8oaQmRN1hSShGvfKRIpWUoLTmUrkS1kJ5SEJG1Urw/v+/d+/D3PzJyZ
ec4zc+bLQF5oawCVVomhu3lv1q/nAWbSnNL05kr8bPrkjr0QD1RbyFT7fazA
vBmu+5VUHhBa0gtU6apA41CbyrVSPHA+1UN5HVbg88INNfnAA/LdZdWMlfrE
T8ny+aImD8R7dQaahVXgzrbShsPHeeAr9zdl8tkKHLNd05xowQP693n/fjCr
wJQfxm8m7Hmg5Jh8zV31CtS9Wtwpf5YH2GcNg7x3VyCr6L+ey3480LWOxaCG
twKrygwG3obygGFikNNCUzm6qReM0K7yQBtlKUR8pb6j9v0a907gAY0y3xTL
XeX4zlVnqj6TB5rh2734iTI8lPht3r6UB+otZ4aWdcpQcPwga8FbHniaPKzU
LnIfG0KHaXQSL+jxRc4vhJRgAEF2u48ELxySCyVZTd3BXbfDpRp288IRff9j
Tfp3MKNt516Hw7ygfNWlJEOsGA3sQvc/0OWFXYX2A+uuFiHbzy6VdSa8IFdv
tdFj4TaeEw3QKHDmBcl/ek6H6wpRvPzN0e+evEAU1kq5t6MQB9Up+uqBvEDa
p9YokHgL4/t8DJPCeYFqpPQ97E8Bari9Mp28zgt8nvK0WfuV859ZxGpPKi+w
x0qdMGnNx/JED5vwPF7gDI23knHKR0eJeofuYl4QcN64dqtgPopU85+hV/KC
yPFLxUuvbqKDAElftZYXEgXXM9G9b2KJ97a9Fi95gTB29YSe6Mr7UmoPa/wg
L3B4MP7c9LyBZtnH0/+sJcByTN4xXeE8zFsyC+XfTABPI/pNv4ZcnDK2OS0r
SIAFYtHPPPdc9NvsIeskRYDpO/dzf9TlYGJoTHOXEQEaE0M/qPZn4cBg4v25
UwRool6JTZLOQqpiZjL7GQI8KI/dNx2WiWVfi20hmAAP27OvJuzIwN+65Zrm
VwhQcOqW3FRIOqoWP5b2YRDg1tzdQeW3adhm+3Kp+BYBEjlrpD8EpiI/to02
3CdAeEZ97/6OFDwp3Ns4+oQAlyVbQmPpKVjgN1yyXE+APktBvv3DyTjbNcng
ayOAXsRhpcS1Kbhn1+z53X0r8y11t56TTMHAa4sndcYJsK8vPVzTIAU51Nfs
CPtJAMrOr21/slLwRN5G7iyWrZBqRPxuXJ+C6X+5f1ZxbAWOEA2hsukU3P6Q
XDdL2grfO7LsHRRSUUNcxcJMcys0+Wn/bm1PxTB2/+hdaVvBue4mJWdNOv5u
nfs8vJ8PRLPCahcLM/Hy7uEoncN8oMMQIxV0ZuLmpFbxJ7p84B/5LMTobybS
LUvskqz5oPPcP9WHx7PQYNp54GgEH0QcuVDn9zcLB3XN/B/H8kGFkiBNQiIb
Hcs0+cXT+GB416PL745nY9B58ROsJXygKPxDXfF2NhaxTTY/6uADX79H438k
clD+zFunbQN8wLWwnmC0Uk9iax1b4gQf3HQ1OXz3RA5q767IZ5lb6f/xlvc6
5xzsTryh5v6TD9ptfuRbBeeg9a/40UFmfjg9qN79ICEHP1tcDNbeyA9/jZPW
bSrKQV/0ID7i5YeKrnCdgcQcZKFaP6GL8APRoX0JQnPwWoSeacI2foj8IVSU
75KD/NMqP5h38cN8pIMJu0kOSpWJyA6q8UNd8VJFp1QOmsFAEcOJH3rGtAde
V2XjxI1mDWZPfihaO6qQdD0b3dmeTLoF8EOguE/iSYdsDG9NFdO6vvIEdsvR
/cqTjWUWxmlMlfzg8esb8p3NQmVUV3Cr4Qd14Sji6OEsbKLs6e5v5Ad+FVH/
IuEsHPrEy/2wjx9qLmnKKr/KRPbzHVdcmQSAfXNWvg0tExMHnm3r3yAAQ7tl
WXYsZyAJ7tdr8AjAfcMmq8WODCy+kWP7QFgALvlZPqm5nYF72GKZqXQBMEn/
yhcZkoHPnIOz46QFQLImwkvfOAOPtrop/1MQgJ+P+aQVVTKwZ5dVv8tBAdiU
Nh/EKpGBNok65/u0BYB+vrn1NXcGzvxU4tMwFID9JvkiScvp6Gexo7LSSgAM
9ga7nZxMxxjKxvlYDwEIWtzN9fVROj653+R4JFkAujzSdPis07F38MzH4FwB
6NFoSmfRTMdvGzY5VhULwIDo96kZmXSUsjE8LYECMNGiH17HlI553GP2G6cE
YE5iQ61HThrWHrg8obYgAAtMezktI9NwwFncPuCPAPzosTPXcE/DrS/c7Ga2
CMLSZfwhAmkoO7dlnCYsCEyWs4fZxdNQT6jS1oouCMxywgmLm9IwynPJ5o2i
IKwf85V5PZyKt7KzxtgOCwLXo/ygh42pWNesagPHVhzb+TrvXiqO/hgf9TMV
BJu1LvodT1LxLyXS+r6tIOR6mma0saahkJ7k6JSrIAyPqn9o1U5DhYDWU2J+
gkDUk9vdwkjDc128pxjRgpC6Y1PTS0o6XmeqGm5OEoTetGXupjPpeGeH+ck1
uYLAv+GTZUN5Ok5eyrXyrhQExkTdwnO1DDSTl7I0HhKEa1n+ok9tM9HXun3g
+kdBaOZwdH5cnIkJ17wsmr4KwsYLhpVVC5nYOvHYfB+bEEQYy2hXhmXh9Bar
fo8tQtBYT0wqb87C9QeYzYuEhGCt3MbR+yv5ejDpiJmQjBCEbJ70uXsjG08+
//TuuKIQPAvqfHZnOhsDZq+ZRh8SgmCfU+FvBXIwVXDXuzpdIcje2VPcpbxS
P6p3mfwxEYKacZ32Tusc7Dzn2ytvKwRD6XXfOy7n4JcsQRM3VyH4a7BfuON2
DnI0V/cU+AqBAHuZantLzko9c8p4OFQICuUZRuU3c1CdsqaHL1oI3p7ndrUx
y0GbY7eM9JKEYH1NXNiWzTmYcWvG8FmRENRGk6yHv2fhZhP/Ezc6hYB0V5Bg
cjMDf1Yx9E/ThGHsC1eg5r9kNPkscuWrlDCsnXGLD/ZKxkeiRRioIAz8U623
KqeS0D8cpRO0hUFnNKaD3JWIy8dnOJ95CK/cH1zivwsZyDSn/kqwRhgYT6eZ
v5+JwVNiHcwFjcLQfiRDU+vQNXxmaLlvV7swzHYcjcsSjsaLTz0Lj4wLg8ZU
CelIyxVcE5UT7sVOhD8Ed5UUqcu4nvb7YKsREd4MffzlNeaLjiaX/U1PEeHI
ouYVIsMbX17dXDbuRITQjcX8DWqeGPWVTl4KIIK3e5+blJYrfqKVmYSHE8Ey
M+NacLEjapsqx26JJcLYS6s7bRx2yIHH/227QYQA9YbOiSxDdFkY2lN+hwjs
qolim19oYAvd2e3AAyJsbvb1MSBLoJTZYv7LWiI0PyhrbbFTgJhrIYMnXhKB
f1+WTZH2YZjFjYSRDiK4lr8JvqSvBse+JR09M0AEx53MmVbGqrDZvORJ6BwR
Lqg13aU1KIJ7jMICxy8i0BQkp96174H2Zy+2p7CIwKhkDCVmcDfEi/el3iWI
gIBEoczMFA14v6/z698vAhmOIT9H3XlQIO6UnHykCGgoXTLUVjmIvQ/4ImLj
RKBx03f5y9KHMGmg9d10mgiIvT+9tVZUHbnFlYPzSkRgV5R2z24mLWzT+db+
54EIuFtWP3T5oo0xnsVUUxSBWRnplIIRHWSvFXi1qVMExs9c2XFtRA+bxtuE
zwyIQJbhc+GiXn0MZ4882zAhAi4qSxyNbQZ4SEblOXlOBKrqU1QLJQ2R1eg7
b+BPEXB63ld9stwEn18oOd3LLApytcKKfPstMCTX7rHsRlGIemL1oPX5SVz+
3HHyE1EURM08ekScbPExd1TZ4W2icNQkPvH3Dnv0U1BdmysjCqFG5ce7vzjg
Yti9ImM1UeA4vvjmmp8zzgpGL9Y5isLPfCW7vy/ccUTj8NXsclEYyUIL105f
zHJbHvz9VBS+E9v3m1f6oUVCuYxRgyhwZo4Kaiafx3fD5G6Od6JgkM7yjmp+
AVPWvhN3HFuJJ8hTxa0cgEaSsRdeTItCciolmUk0EDt8/pL8/4rCy+RDhgNj
QajR2vfoEysJxrYayr2qC0Zx6z22KetJIPracO2wfAiyfYvjUOckge6NJVUP
YghOhs9ULmwhQah/bhDr2hCsF9A8mbuVBNY3Ynwk3wXjzTs31x8TIsEL/b65
Y4XBaNNhbl5EJUHEXY+xFPVgVLV/uMZEggQTFtXm1YRgJP3kvrtWigTaGze8
HR0PwiHiKybrvSTgdMxp3B4WhGZOCgV8WiRgkN/ceVobiMeZCPOXzpIgMWGH
74hCAAb4am1M8SIBccOB2FHWALw5G0wrPk+C0sBjt8daL+C3gU+m7WEk6HY4
1z9udwETHuFzYgoJyPuqVKZi/bHrnEvCg5Vb12Nm0zOL2364PJV792U9CWrf
czAxufgh1bqnaeAVCbj7NqjckPJDr2MH/7C8JcHdBtbqqXJf5N3B73BsigRf
s39U+dT4oNINnRCbWRIcTPr2QyDUB+0Ew9K8F0iQED2/p1rNByvWz7am/yGB
3PnpCtaX3mg48WLvxy1kcDYYLo3t9MIgi1/6S1vJcE9jYE42yQsLOqVcuITJ
wKzyTqrHxAtbtezCyWQy6Mp3u/oLeeH3Z6k5cnQy5Ep23iEOeaLIvjePj0iS
YZ7cNo05nqheuuatmQwZknM6lN3GPNFtm+KcqzwZxvl56PcNvTAp6+yGUEUy
yMQd5/rW5IU1hHxKggoZQjckfN+j5I2T0X3Ktw6RoSW0a+j8PW/cc+HwuRZd
MtifM7z3L9EHu42F+gVsyZBw3Fhvo60f/m3VW9zhSIbR5mQF3W4/pKuHbwJX
Muw81EuK0zyPPnLzag6+ZGiQM53fKuOPW7c03im7SoYfBPN4sT8X8EDkckN9
LBkOxaT727sHoMO/XaO9iWSIWzdgU/g+AB/MZBD+ZZNB/IeFrNSrQByy65Dm
zieD59lMIQ/lIFw3wKZFKyLD0w+DrBWlQWj86lyQdgUZzvZade1LXsm3wxnf
2B6RQeK30hkH5hD8VFvvXFdNBnqk1mfpvSFYoTg3EvKcDGN/Wl9TQ0MwqJLf
WLmRDKa3bw/ad4WghszBll/NZGg/cWm2YFsochefUXvQRgZNppP/PviHYkFO
jbR0Hxn2G/OSnckX8ZsP+yg3kxj0lU3dOvfkEsZT8tT6tMVgv3NCZdjcFbTI
an6UpCcGZYSzDXUHo5AusCh93HCl/cPh5fHoKHzMpSH82mqlfUGxdY3gVRz9
NfOt2kMMpi+qSXrsiUaZtn0FucliYK7kWdPtGoO/tWyFrTLEIEr2nfuLihis
q4+OF8wVg8fbVSilSzFo8mQ4iFEkBvwCGyOvXLmOwQWXjS9Xi8GEf7t4lkEs
apDvtag+FwMfnzUXeX1ikTujV+1vgxhwntvbH5UaiwVxkjK+bWKg5JR5zXck
Ft05DAtk34pBl13b5AxrHCpGBAl/eScGLqdYwW5bHK5hKYy/MyQGayz2pPZp
xWFrQPt6p/dikG7s9FXPLQ5H+7vmoj6KQenr1rHzkXH4TbG3+86MGDzdm3Gx
OTsOBX4N3fzyQww6OPc+d2uNQzvvz+DHQQGmxy6qlZLx6Ns1ty1tCwU4aYqj
69TiMUp2gevpVgoIxLKFmpjFY+n8r4F/JArIOuThUmQ8LrmynY+Qp4ANd6/K
wcl45GphP3VbkQJugfnDjD/xSNrBdaRZhQI+H88FT/AyUP0TL2GTFgViazhr
I9UYaKrJvyxzjAKZEv1WfWYMdCkUem9wggJFCYV/Jc8xMN5B7H6SFQXqnNUO
vMlhYH49LaXKlgJtbzcPkaoYWEWVCO5zpMAgDAWee8PAV2E7HJZdKTBeXEys
m2Tg4Ji0jognBea3nq8m/GXgF1VZOfCjQPood0sbVwKy5O4RsgmkwExMNSmI
IwEJTIosly5SQEXJ0UuSPQHFTypP5UdQIG6Ku6mXLQH310BbY/TK+EnVQuFr
E1CXeOjhVBwF5A85npVlTUCvfq3LOzMo0JdVvTVmJV6tsQm5spQClBOO9mqL
DOx4YL6+p3JlP5h5Hn35ysBxwsm5n48p8PJuNUfWFwayd9lXK9VTwHUDT8XP
aQYa6Xma1r2jQEVNNcvdMQY63fOBySEKsLk4GpmPMDCAy3/b+nEKmAjwFK0f
YmDu65BF7VkK/PF01Ld9x8AKyUsDrgsUOEbmyd/cw8CGqIgX139SIK+1+md1
FwPfTUUV3f9DgcULjkfPdDDws0ZMXCcLFdQleHL42xj471ac3yIbFVK7qxfq
WxjIzZZ4ko+TCqfg9k7lCgbSHFLU93FTgTtpd2RvLAMV6tN3mvNRoW76yZin
KwO1qdm8gcJU8FE9rLxJi4FWYXlLWWQqSCS3JhfRGRimevvlmCQVEg6OHh0d
jse61IM1n2Wo0HfkxitydTyu+9pf9l2eCnw69ho2aSv5mMeVwQ5UYBhPqb0/
EY+xrF5uu09QIdllXnHiZRx2mHPaKplS4a1HeRXtVhzyVhQYq1tRgdfXe4/D
pThMtnsHZo4r7UN/7fpwIA6z6g/wXAxYiZ/ELDFVHosjxF626FAqtKe/KJCI
i0Wyj8dyYvhK/NzLVGe3WLRpZZ/PvkoF84IjeUXasXiTfnPiduzK/hRvIE2L
x+JkkHJfeSIVekubMyXXxaJ4T3drddrKel2E9vY2XEdnafcXjdlUeJ5rH5Bu
fR2LIzZUtd+kwmcjz1qN+Rj8PJx3p//2ynw5Qlm+h8SglIJS7sRdKtCfxRy+
sTkGyz66Rf2qosJRyaLmP9LXsEqzy4LcsjLev5GeMrer+DvXVV+yY2X9FbOC
p/5GodLvderyPVT46LRsyRkThbVF+6Q1R6lA6eIbP11yBRs2ZjN7LFIhrFD3
i/DnCOxscSpAIg0s9J+uj3C6hLP6/+ZOutFAw5KkUhh2ATtKooHHkwZhN5ne
j9T6Y+UGwbgGXxpMTY+ECyyfx4Bncrt3htIgxj+3NcrLD9l3OXv9ZtCgO1nM
ytXeG2ejf9SVpNBgyzALS0GuF7Z/vESwzqSBH/19/tCgJ6bkZD9oLKABsfLG
zDGjc3hheQfbheKVeMthcZHxHmhl8thYqpQG02p28s9a3ZG2+e3PhEc0EM7X
/G0X44brXWw0NWpo8EGW/WD3Olf83DiXuvScBpfDJq8VPjmDbyiBn+420uDr
FaHF2jfOWBbMvt/mNQ3MYvUset47YVJ/8lVCOw2uJl1+MfvDEc/vpQ00vaXB
csbj7es4HBHmVAKlh2mQeJv6S07OAS1ai5hfraEDvSrVfuS2JcJ2BYPADXSI
rGOV29RhjpTw+jwZLjqYtbuyHPhtip+UR9SS+ehQ+Ek1O13LCP1KCOF22+kA
il33S7p10GLDjW4+aTqQjxfqbtymhSr2MtuaZelwyyVg2tFXHdcRtZt2KdMh
8YCmlobMAfzk18s/oUqHY+7zYazte7Cly94pRZ0OL9RlLwS+l8JSmYVH2tp0
0CsJm12gkjAhOoT93zE6bOBUJRY4TNb6feQ0LztBB+0HClqsLVvB/FB6sb0p
HR5sPDd7spMEB3LEl/mt6KD4tHHvAaCD2HLl0dc2dLBxJYYKl0jCOhO1zODT
dPC8vs/+T7QMTJW3zex2oUPo6fErLmqyUHpmOibVmw7N/IqdWqUK0KxlxCOX
vhLfTWN5dNcR+BlV3XU0hw7We5ZGrYo0gNpMTbLPp4Pq35LGfjEtCNBe4E++
R4e9MTyMtzw6UHjVrO9eBR10jRr8DKJ14W3zs/SmR3RQv6507zSLHkgdjRX5
/YIOHsZzCeP++mAW/WOY+yUd3J40sVH2GkDEa6tcyVY6pIvm+VsvGEAFR4PN
oU46PAq7MJN97ziMHN1JteylQ8OHE6eGzpwAjmsJE96DdKjXluoUFjeEfS1L
BTFjdBi7x6ZuPm4IDpy2jrc+0OH//sON13klgZ9XXXNt13Tv/KqnW1LuzH9f
9SHd09IU5m3/eVF3jXw4YdUmegfg6IFVX4rNZ7ZXW/X9No4XgRqrZtfvP3zP
YNVP9f20uU+vWvR4mVHv9VVPGtJd7UZXXZEaGkWfXPXFwYFbHz/9v/72jDGX
b/9vfl5MJt5s4v/5VXzfwfCdqy55E8NfeH7VXpo/XkxvkfjPVZHqIT5K2/+z
Yh+Xp2WI5H+W4u21Sb+74z9rRfz+NlK+aoffwhH0R6vOHrEuvv9i1ZtLphca
e1c9f/jfpW8sO1fX60u/pWO86rZPGvsYFqv+bHmmudd61bRDpXO2LqtO2qyo
4B+6ar/bOi8LiledIOxu/rl01aXX42d2PVj1lHfPlmpctdlBa7POt6v2qQz7
zN+/6njxgiCrkVW/4prOm/q06slQTnnpL6tmWZRu9FpctYijgenj36v+3//L
//P/AOfAHKg=
          "]]}, Annotation[#, "Charting`Private`Tag$5572#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$5572#2"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.560181, 0.691569, 0.194885], 
          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$5572#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.8024014038220644`*^9, 3.8024014462289486`*^9, {3.8024015995619917`*^9, 
   3.8024016011977034`*^9}, 3.8024016434732265`*^9, {3.8024019932563252`*^9, 
   3.802402002397108*^9}, {3.8024039017556725`*^9, 3.802403953038184*^9}, 
   3.8024039997331743`*^9, 3.8024040353894863`*^9, 3.8024040711283407`*^9, {
   3.8024402724085836`*^9, 3.8024402836382747`*^9}, 3.8024607434387436`*^9, 
   3.8024619217099557`*^9, {3.802472726908188*^9, 3.8024727337617497`*^9}, {
   3.8024733705349693`*^9, 3.802473373799682*^9}, 3.8024737957305717`*^9, 
   3.802478092525135*^9, 3.8025242346594334`*^9, 3.802601077532339*^9, {
   3.802602251946018*^9, 3.80260227745997*^9}, 3.8029665879470177`*^9, 
   3.8029674795280247`*^9, {3.8029679757230563`*^9, 3.802967981600068*^9}, 
   3.802968296817315*^9, 3.8029687295761766`*^9, 3.8029687959981074`*^9, {
   3.802969121890779*^9, 3.802969145373878*^9}, 3.8029694168199434`*^9, 
   3.8029694734614067`*^9, 3.8029696107145567`*^9, 3.802969976269147*^9, 
   3.802970032875019*^9, 3.802970319834959*^9, 3.8029709262354665`*^9, 
   3.802970974585147*^9, 3.8029710063948755`*^9, 3.8029729670059357`*^9, 
   3.802973168461779*^9, {3.802973512807762*^9, 3.802973519461662*^9}, 
   3.802973725744069*^9, 3.8029738355568533`*^9, {3.8029739243018894`*^9, 
   3.8029739533144197`*^9}, {3.80297408530324*^9, 3.802974090797762*^9}, {
   3.8029745874507174`*^9, 3.8029745913411255`*^9}, {3.802974672041303*^9, 
   3.8029746760393763`*^9}, {3.802975137169429*^9, 3.802975141481453*^9}, 
   3.802975853867917*^9, 3.802976543068199*^9, {3.8029767734838405`*^9, 
   3.8029767853323045`*^9}, 3.8029784587071514`*^9, 3.802978520093154*^9, 
   3.8029788328760366`*^9, 3.8029799539280834`*^9, 3.802980044539624*^9, {
   3.802982408791937*^9, 3.802982414320069*^9}, 3.8029825665803432`*^9, 
   3.8029832876103344`*^9, {3.8029838341483326`*^9, 3.802983852301343*^9}, {
   3.8029838825646124`*^9, 3.802983902018034*^9}, {3.8029841509364934`*^9, 
   3.8029842204931808`*^9}, 3.8029842651180744`*^9, 3.80298449443804*^9, {
   3.8029874183166347`*^9, 3.8029874525192504`*^9}, {3.8029875739488344`*^9, 
   3.8029875799306946`*^9}, 3.8029881061848893`*^9, 3.8030564471159*^9, 
   3.803056552250534*^9, 3.8030568611698017`*^9, {3.803057362952325*^9, 
   3.803057375942271*^9}, 3.8030576502115116`*^9, 3.803057705325587*^9, 
   3.803057743903199*^9, {3.803057898781231*^9, 3.803057917895503*^9}, 
   3.8030580261283503`*^9, 3.8030581784057837`*^9, 3.803058255242637*^9, 
   3.8030586088595085`*^9, {3.803059465918867*^9, 3.8030594790683317`*^9}, 
   3.8030595433037653`*^9, {3.80305962649969*^9, 3.80305963005138*^9}, 
   3.803059666417262*^9, {3.8030596998144646`*^9, 3.80305971096183*^9}, 
   3.8030597454969964`*^9, 3.803059833492078*^9, 3.8030599163932285`*^9, 
   3.803060035405817*^9, 3.8030606723155203`*^9, 3.803060960180106*^9, 
   3.8030691381239977`*^9, 3.8031340059102*^9, 3.8031393795100307`*^9},
 CellLabel->"Out[40]=",ExpressionUUID->"5ca0a150-0496-45d5-8dff-c91e3146c076"]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Scan the intensity from 0 to 3000 ",
 Cell[BoxData[
  FormBox[
   RowBox[{"mW", "/", 
    SuperscriptBox["cm", "2"]}], TraditionalForm]],ExpressionUUID->
  "b5579f26-0243-4cd8-9fb1-4614368af08f"]
}], "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.803069143535514*^9, 
   3.8030691512885995`*^9}, {3.8030692098716726`*^9, 
   3.803069210746869*^9}},ExpressionUUID->"c44a9341-abbe-437f-a7b8-\
e777560bd4ec"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Is", "=", 
  RowBox[{"Table", "[", 
   RowBox[{"x", ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "6000", ",", "500"}], "}"}]}], "]"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "power", " ", "is", " ", "2", "x", " ", "of", " ", "that", " ", "listed", 
     " ", "in", " ", 
     RowBox[{"Fig", ".", " ", "S6"}]}], ",", " ", 
    RowBox[{
    "as", " ", "intensity", " ", "defined", " ", "here", " ", "is", " ", 
     StyleBox["peak",
      FontSlant->"Italic"], " ", "intensity"}], ",", " ", 
    RowBox[{
    "whereas", " ", "that", " ", "given", " ", "in", " ", "paper", " ", "is", 
     " ", 
     StyleBox["average",
      FontSlant->"Italic"], " ", "intensity"}]}], "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.802378015195921*^9, 3.8023780618768625`*^9}, {
   3.802399859883916*^9, 3.8023998603848224`*^9}, {3.8030691557826223`*^9, 
   3.8030692064249*^9}, 
   3.803305875206585*^9},ExpressionUUID->"765a6f84-5baa-41bd-a875-\
c58ef6d443ae"],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0", ",", "500", ",", "1000", ",", "1500", ",", "2000", ",", "2500", ",", 
   "3000", ",", "3500", ",", "4000", ",", "4500", ",", "5000", ",", "5500", 
   ",", "6000"}], "}"}]], "Output",
 CellChangeTimes->{{3.802378057419541*^9, 3.8023780681395884`*^9}, 
   3.80239993288593*^9, 3.8024021412926517`*^9, 3.80244042565911*^9, 
   3.8024607801825094`*^9, 3.802461927920744*^9, 3.802472836851163*^9, 
   3.8024733799319644`*^9, 3.802473811105195*^9, 3.802478099526806*^9, 
   3.8024781885495105`*^9, 3.802524297771593*^9, 3.8026010807655983`*^9, 
   3.802966602847678*^9, 3.8029680539569335`*^9, 3.802968301953377*^9, 
   3.802973468206977*^9, 3.802988958419587*^9, 3.803060802354598*^9, 
   3.8030610837815933`*^9, 3.8031340061760507`*^9, 3.8031393796920233`*^9},
 CellLabel->"Out[41]=",ExpressionUUID->"5dd4a8b1-dab2-4ce6-a7d8-e66962f5a6f5"]
}, 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.803069908131694*^9, 
   3.803069933888533*^9}, 3.8031286523487816`*^9, {3.803133998229158*^9, 
   3.8031340176633515`*^9}, {3.8031343034423475`*^9, 
   3.8031343034851246`*^9}, {3.803139372102103*^9, 3.803139372207983*^9}, 
   3.8033058796999826`*^9},ExpressionUUID->"5ceef669-cf9d-4c9f-996a-\
85e2b663751e"],

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", " ", "lower", " ", "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.8031340076664524`*^9, 3.803134023792755*^9}, {3.803305887397764*^9, 
  3.803305890149012*^9}},ExpressionUUID->"248a659e-539d-4868-8b3f-\
dda1f95b21ed"],

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[{"#", ",", "1.4", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", 
              RowBox[{"Intensity", "/", "isat"}]}], "]"}], ",", 
            RowBox[{
             RowBox[{"-", "15"}], "/", "\[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[{"Intensity", ",", "Is"}], "}"}]}], "]"}]}], "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.8024735600002484`*^9, 3.802473560058296*^9}, {
   3.8029669679940057`*^9, 3.8029669680549335`*^9}, {3.8029679917810287`*^9, 
   3.802967994034485*^9}, {3.8029681456571207`*^9, 3.8029681812982955`*^9}, {
   3.802973470774012*^9, 3.8029734722470617`*^9}, {3.802988264875614*^9, 
   3.8029882970454865`*^9}, 3.8029890051427517`*^9, 3.803063814117606*^9, {
   3.803069214020602*^9, 3.8030692476459937`*^9}, {3.8030698638483257`*^9, 
   3.8030698662651987`*^9}, {3.8031224981800427`*^9, 3.8031225027900076`*^9}, 
   3.8031343172942295`*^9, {3.803305914282776*^9, 
   3.8033059144282274`*^9}},ExpressionUUID->"ccc1bdba-34c7-4530-8f76-\
36d6baa364a0"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"widths", "=", 
    RowBox[{"{", "}"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"produce", " ", "width", " ", "from", " ", "beam", " ", "images"}],
    "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dist", "=", 
     RowBox[{"SmoothKernelDistribution", "[", 
      RowBox[{
       RowBox[{"image", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"InterpolationPoints", "\[Rule]", "4000"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"pdf", "=", 
     RowBox[{"PDF", "[", "dist", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"iF", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Cases", "[", 
       RowBox[{"pdf", ",", "_InterpolatingFunction", ",", "\[Infinity]", ",", 
        RowBox[{"Heads", "\[Rule]", "True"}]}], "]"}], "]"}]}], ";", "\n", 
    RowBox[{"Peaks", "=", 
     RowBox[{"Pick", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"PeakDetect", "[", 
        RowBox[{"iF", "[", "\"\<ValuesOnGrid\>\"", "]"}], "]"}], ",", "1"}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"PeakValue", "=", 
     RowBox[{"Max", "[", 
      RowBox[{"pdf", "[", "Peaks", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"GreaterThreshold", "=", 
     RowBox[{"Position", "[", 
      RowBox[{
       RowBox[{"pdf", "[", 
        RowBox[{
         RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", 
       RowBox[{"_", "?", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"#", ">", 
           RowBox[{"PeakValue", "/", 
            RowBox[{"Sqrt", "[", "E", "]"}]}]}], "&"}], ")"}]}]}], "]"}]}], 
    ";", "\[IndentingNewLine]", 
    RowBox[{"totalDist", "=", 
     RowBox[{
      RowBox[{"Last", "[", 
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "-", 
      RowBox[{"First", "[", 
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"left", "=", 
     RowBox[{
      RowBox[{"First", "[", "GreaterThreshold", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"right", "=", 
     RowBox[{
      RowBox[{"Last", "[", "GreaterThreshold", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", 
    RowBox[{"diameter", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"right", "-", "left"}], ")"}], "*", 
      RowBox[{"totalDist", "/", "4000"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"diameter", "/", "2"}], "]"}], ";", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"widths", ",", 
      RowBox[{"diameter", "/", "2"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"image", ",", "images1"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{3.801347246666169*^9, 3.802139092022823*^9},
 CellLabel->"In[47]:=",ExpressionUUID->"bc6dac95-1b93-4ce8-9af6-755f797db9e8"],

Cell[BoxData[{
 RowBox[{
  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[{
   RowBox[{"rinterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"rinterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "0.625"}]}], ";"}]}], "Input",
 CellChangeTimes->{{3.8030604642359715`*^9, 3.803060478836073*^9}, {
   3.803060540495633*^9, 3.803060560226141*^9}, {3.803060612223385*^9, 
   3.803060612346784*^9}, {3.8030609173244953`*^9, 3.8030609174005117`*^9}, 
   3.8030636743995957`*^9, {3.803063835279997*^9, 3.8030638369861684`*^9}, {
   3.8031345851852107`*^9, 3.8031345906373816`*^9}},
 CellLabel->"In[49]:=",ExpressionUUID->"e1d13343-c95f-46a7-bd96-775b2205822d"],

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[{"#", ",", "1.4", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", 
              RowBox[{"Intensity", "/", "isat"}]}], "]"}], ",", 
            RowBox[{
             RowBox[{"-", "15"}], "/", "\[CapitalGamma]"}], ",", "50"}], 
           "]"}], "&"}], ",", "init"}], "]"}]}], "]"}]}], ";", 
    RowBox[{"(*", 
     RowBox[{
     "2", "I", " ", "since", " ", "there", " ", "are", " ", "two", " ", 
      "beams"}], "*)"}], "\[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[{"Intensity", ",", "Is"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{{3.80213921204575*^9, 3.8021392423139696`*^9}, {
   3.8023782691198444`*^9, 3.8023782718159704`*^9}, {3.8023783477688518`*^9, 
   3.8023783711872344`*^9}, 3.802400156327544*^9, {3.8024003452528543`*^9, 
   3.802400346292104*^9}, {3.8024735655923023`*^9, 3.8024735656749964`*^9}, {
   3.8029885726028247`*^9, 3.8029885728517056`*^9}, {3.802992534566959*^9, 
   3.802992535371936*^9}, 3.80306062750261*^9, {3.803069224491373*^9, 
   3.803069253031706*^9}, {3.8030698713841105`*^9, 3.8030698737564335`*^9}, {
   3.803082956251447*^9, 3.8030829582363*^9}, {3.8031225060097294`*^9, 
   3.803122508124203*^9}, {3.8031355534325666`*^9, 3.803135553479582*^9}},
 CellLabel->"In[51]:=",ExpressionUUID->"ab82b3ea-5c09-40e3-a005-6f650494affb"],

Cell[BoxData[{
 RowBox[{
  RowBox[{
   RowBox[{"widths", "=", 
    RowBox[{"{", "}"}]}], ";"}], 
  RowBox[{"(*", 
   RowBox[{"produce", " ", "width", " ", "from", " ", "beam", " ", "images"}],
    "*)"}]}], "\[IndentingNewLine]", 
 RowBox[{"Do", "[", "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"dist", "=", 
     RowBox[{"SmoothKernelDistribution", "[", 
      RowBox[{
       RowBox[{"image", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", "1"}], "]"}], "]"}], ",", 
       RowBox[{"InterpolationPoints", "\[Rule]", "4000"}]}], "]"}]}], ";", 
    "\[IndentingNewLine]", 
    RowBox[{"pdf", "=", 
     RowBox[{"PDF", "[", "dist", "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"iF", "=", 
     RowBox[{"First", "[", 
      RowBox[{"Cases", "[", 
       RowBox[{"pdf", ",", "_InterpolatingFunction", ",", "\[Infinity]", ",", 
        RowBox[{"Heads", "\[Rule]", "True"}]}], "]"}], "]"}]}], ";", "\n", 
    RowBox[{"Peaks", "=", 
     RowBox[{"Pick", "[", 
      RowBox[{
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"PeakDetect", "[", 
        RowBox[{"iF", "[", "\"\<ValuesOnGrid\>\"", "]"}], "]"}], ",", "1"}], 
      "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"PeakValue", "=", 
     RowBox[{"Max", "[", 
      RowBox[{"pdf", "[", "Peaks", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"GreaterThreshold", "=", 
     RowBox[{"Position", "[", 
      RowBox[{
       RowBox[{"pdf", "[", 
        RowBox[{
         RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
         RowBox[{"[", "1", "]"}], "]"}], "]"}], ",", 
       RowBox[{"_", "?", 
        RowBox[{"(", 
         RowBox[{
          RowBox[{"#", ">", 
           RowBox[{"PeakValue", "/", 
            RowBox[{"Sqrt", "[", "E", "]"}]}]}], "&"}], ")"}]}]}], "]"}]}], 
    ";", "\[IndentingNewLine]", 
    RowBox[{"totalDist", "=", 
     RowBox[{
      RowBox[{"Last", "[", 
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}], "-", 
      RowBox[{"First", "[", 
       RowBox[{
        RowBox[{"iF", "[", "\"\<Coordinates\>\"", "]"}], "[", 
        RowBox[{"[", "1", "]"}], "]"}], "]"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"left", "=", 
     RowBox[{
      RowBox[{"First", "[", "GreaterThreshold", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"right", "=", 
     RowBox[{
      RowBox[{"Last", "[", "GreaterThreshold", "]"}], "[", 
      RowBox[{"[", "1", "]"}], "]"}]}], ";", 
    RowBox[{"diameter", "=", 
     RowBox[{
      RowBox[{"(", 
       RowBox[{"right", "-", "left"}], ")"}], "*", 
      RowBox[{"totalDist", "/", "4000"}]}]}], ";", "\[IndentingNewLine]", 
    RowBox[{"Print", "[", 
     RowBox[{"diameter", "/", "2"}], "]"}], ";", 
    RowBox[{"AppendTo", "[", 
     RowBox[{"widths", ",", 
      RowBox[{"diameter", "/", "2"}]}], "]"}]}], ",", "\[IndentingNewLine]", 
   RowBox[{"{", 
    RowBox[{"image", ",", "images2"}], "}"}]}], "]"}]}], "Input",
 CellChangeTimes->{
  3.802139208841636*^9, {3.802378377457013*^9, 3.802378377557906*^9}, 
   3.8024001591204295`*^9, {3.8030829691535006`*^9, 3.803082969223692*^9}},
 CellLabel->"In[53]:=",ExpressionUUID->"6ce184cd-c0f6-464e-919c-89cc33289447"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{2560, 1357},
WindowMargins->{{2552, 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, 20699, 394, 547, "Input",ExpressionUUID->"7824fddb-3943-46b5-acd8-799c34acddf8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22053, 437, 415, 6, 65, "Subsection",ExpressionUUID->"c541093d-3069-439e-a54a-d38289364e3b"],
Cell[22471, 445, 6290, 183, 372, "Input",ExpressionUUID->"a22cab01-84b6-4910-a8ee-e13c7c9bda3f"],
Cell[CellGroupData[{
Cell[28786, 632, 496, 7, 53, "Subsubsection",ExpressionUUID->"2086d6f6-0d0b-40d9-a62f-2851fbd41666"],
Cell[29285, 641, 4063, 95, 420, "Input",ExpressionUUID->"faefb1d2-c143-4cf4-bbad-8f8982efcabd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33385, 741, 524, 7, 53, "Subsubsection",ExpressionUUID->"20effcbb-203d-40b4-8029-107270944d9f"],
Cell[33912, 750, 9695, 192, 557, "Input",ExpressionUUID->"dccf9f9d-0271-4cee-a2a6-51b034d09fcc"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43644, 947, 578, 8, 53, "Subsubsection",ExpressionUUID->"6c519e3a-f10c-4538-afc2-24105e47eb71"],
Cell[44225, 957, 2756, 39, 35, "Input",ExpressionUUID->"9e84a53b-6a94-415f-8b47-20c1efb51727"],
Cell[46984, 998, 8852, 163, 372, "Input",ExpressionUUID->"c0c8e4fa-0227-4b9a-82f2-c524254ed188"],
Cell[55839, 1163, 2858, 43, 35, "Input",ExpressionUUID->"c5609752-0026-4bbc-8872-3b54364462a9"],
Cell[CellGroupData[{
Cell[58722, 1210, 16005, 234, 84, "Input",ExpressionUUID->"e65b6191-864b-4f9a-bf41-1a27cb9ce5bb"],
Cell[74730, 1446, 50968, 899, 450, "Output",ExpressionUUID->"5ca0a150-0496-45d5-8dff-c91e3146c076"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[125747, 2351, 1049, 19, 56, "Subsubsection",ExpressionUUID->"c44a9341-abbe-437f-a7b8-e777560bd4ec"],
Cell[CellGroupData[{
Cell[126821, 2374, 989, 25, 35, "Input",ExpressionUUID->"765a6f84-5baa-41bd-a875-c58ef6d443ae"],
Cell[127813, 2401, 882, 14, 39, "Output",ExpressionUUID->"5dd4a8b1-dab2-4ce6-a7d8-e66962f5a6f5"]
}, Open  ]],
Cell[128710, 2418, 1083, 23, 35, "Input",ExpressionUUID->"5ceef669-cf9d-4c9f-996a-85e2b663751e"],
Cell[129796, 2443, 1362, 34, 60, "Input",ExpressionUUID->"248a659e-539d-4868-8b3f-dda1f95b21ed"],
Cell[131161, 2479, 9166, 173, 252, "Input",ExpressionUUID->"ccc1bdba-34c7-4530-8f76-36d6baa364a0"],
Cell[140330, 2654, 3247, 83, 324, "Input",ExpressionUUID->"bc6dac95-1b93-4ce8-9af6-755f797db9e8"],
Cell[143580, 2739, 1414, 34, 60, "Input",ExpressionUUID->"e1d13343-c95f-46a7-bd96-775b2205822d"],
Cell[144997, 2775, 3726, 93, 228, "Input",ExpressionUUID->"ab82b3ea-5c09-40e3-a005-6f650494affb"],
Cell[148726, 2870, 3350, 85, 324, "Input",ExpressionUUID->"6ce184cd-c0f6-464e-919c-89cc33289447"]
}, Open  ]]
}, Open  ]]
}
]
*)