TrajectoryMonteCarlo_IntensityScan_Cooling.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[    155754,       3024]
NotebookOptionsPosition[    152786,       2970]
NotebookOutlinePosition[    153160,       2986]
CellTagsIndexPosition[    153117,       2983]
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_Cooling.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.8030827195662594`*^9, 
   3.8030827225163555`*^9}, {3.803131210995144*^9, 3.803131212978491*^9}, {
   3.803135791241442*^9, 3.803135793074301*^9}, {3.8031379630141897`*^9, 
   3.803137963715014*^9}, {3.803138407541874*^9, 3.8031384086419964`*^9}, 
   3.8031388795895243`*^9, {3.8031390157900543`*^9, 3.8031390162667294`*^9}, {
   3.8032144803938303`*^9, 3.8032144816131024`*^9}, 
   3.803305692006796*^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.8030917390104523`*^9, 3.8030917391608152`*^9}, {
   3.803091786308816*^9, 3.8030917863848333`*^9}},
 CellLabel->"In[15]:=",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[28]:=",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.8031312735525904`*^9, 3.8031313015696964`*^9}, {3.8031313388965616`*^9, 
   3.803131346898845*^9}, {3.803135904376921*^9, 3.803135907203696*^9}},
 CellLabel->"In[29]:=",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.80313593795884*^9, 3.803135938224117*^9}, {3.8031363882907095`*^9, 
   3.8031363885769367`*^9}, {3.8031372390561175`*^9, 3.803137239440917*^9}, {
   3.8033057007848897`*^9, 
   3.803305701007286*^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[32]:=",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[32]:=",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[{"2.4", ",", 
       RowBox[{"QuantityMagnitude", "[", 
        RowBox[{"2", "*", 
         RowBox[{"800", "/", "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.803121444037858*^9, 3.803121447285568*^9}, {3.803135862127427*^9, 
   3.8031358625615396`*^9}},
 CellLabel->"In[33]:=",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: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$5005#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$5005#2"]& ], 
       TagBox[{
         Directive[
          Opacity[1.], 
          RGBColor[0.560181, 0.691569, 0.194885], 
          AbsoluteThickness[1.6]], 
         LineBox[CompressedData["
1:eJxE13c41e//OHAre++9d0bJLDp3KsmOSla8Q4Q0iBBRiJKojIREJJWVbM7T
3hx7771HxjlGfr7X9bs+nX/O9bhe537e133fz3O/nk+hW/eM7UmIiIgKjxER
/d/3mUH3xcNDWSD6/5++os+rl/f+2SO29c+77X/Ovia5L7H8z6K4QVqjgX+m
rdGUS8r755FMxgfozj8HPPux7d8h8z+Hcjg3pb46/j/zV0umuJlL/88Hoo5s
Z+Wl/mfVs6atx6X/+aGp1gsusX+eCxUlbHL/c8fK2NB38n9OKzT7wjUi+T8b
6Ome2gr75xD700sTL/65yk8qFffsn1VzKDh+PP5nYfbq3Vt2/7w9ql6JO/PP
iW5yV34sSPzPKqFSmWvT/4xLFKVRGv9n4kbumrLef7YTpFBrq/5n6ZZRwY3E
f66aGPBVjvtnS3z3gHfUP78WaX5PGvbPa16FlGxe/5wvHrmmYvLPhuph+k8M
/nnuyosMuPzPXL6+dpfRPz/pvNNnLvvPbPN2Sp8k/znzr/XbSZF/HpW6rufC
/c+a/prgS/HPg1EafJUk/+z+XdWb/K/4/0xbcapXh/DPqT1yim82//nsklRk
5+o/9xCLrXAs/vM9DkFdy5l/ppTlSU8aF4eyvTpXvul8SNJkPzY9JA7Mnkp/
/PfyQe0G0y2pPnH4byPZc4qpAJyeU/DmtorDsXlf/+9nC6Cl/8+bqhJxuNRz
6q1abAHw0rA9Sc4XhxJ3n/Un2QXgoq7sGJAjDlIsVUZQXwC0nx4jzFdxIDcy
ZrhEKAA9u4O14rfiAA33w66aF0J8FP/wh9fiIOtYuBjzsBCWajENj0PEIZ6c
SHfwZSG8knr2WdlPHHzPv6G6VVIIjSvkxrl3xEG97GfwPd4i4BaUPBtpJw4Z
FlvTuYpF4HTlsvR9a3Hg2lW/uK1XBMXPnNgNzcXhRWzgFzX7IqDKe0Uid00c
dpWbSX19i+DG9I8VWiNxqLrhY/QyvgjS2VsHFnWOxqfmPdAvLwJtL8Zf35A4
VGtI5HUSF0MtKTfrZw5x8E9aL1tiKAF2pTNE/szi0LskPZp1sgRu37Zcukl3
tB41OyI3kxI41pBYzUsqDqPtPZqE6BK4EC7yKHZFDM6SltYd4y+FSk7Znjc1
YrBnH9zOL1MG5bIo1s9NDPJ0CzRIb2Nh1L7J45GrGKimz2ewhWKBKPH6NZc7
YlBCxssp+QMLmvR3mc2txaCqzH9DbwMLtSsfwpR0xWCeI9VEFwMwKy7mLKMl
BtJaZ192GQFQWmdfFjknBi7uvWB1C0Aq5owkt7oYZCXf35lxA9BpqyVnUhGD
NRyV3P0gAGcK42kKBTFQOEy2I0QDhGGGq/7KiIHO6uv5j6kAPzwdk7ckxCBu
J9qNPhOgJeuP/5KwGMwdfjrwzwdYmfWznuQTA2XKby/+lAPQC1KfHeAUg2DG
XKbbdQBGEYJ7dfRiICxULabbB9B1X/9Ryq4oWGjO6bPtAgyd/HrHvEMUWgJP
JH7TqICDOwraV1pEgee1miSvVgXwJ5eJa9eLwp0ozdxwgwqwYe6aVCoXhWNp
V2vdrCsgQMe6UqZIFK5mWhnNOFRA8rOFJJE8UUjOvz1w434FTG4Q3WTKEAVU
93hFI6ACyI6HqVOmikJ4W8Dj7NAKELXl4Dn8JAojvS9JRN5WwMWPyYStOFHA
5ai7FMRWwO1O2b6lKFHQG71W7XfkFzRF+ZMRolBHe49X68jp5y9EDbwShfOn
Q9zpjrzwy9y4PkAUzkQVi8bHVADN0vQJ7BNRyK/sfGJ7ZBnRBwz5nqKgsLbU
JX1k1/cvmlPuioK0rkBQUXQFbLjnaT01EwVO4sDpnqgKICjRn1E+KQrbLvpp
TO8r4LFTatlzGVGYdT/RMvSuAvCJ6qhdQhR6n7Bsfj3yNoXzRRd+USgM6z93
9siP1Enqi7hE4dv7UkeqI2/e/6BDwSYKH+I/vek6Wq9b6omWq4yi8PLLs/xP
R97orzNMphEFnx/2w05HfkBv3bFKfrRfpHOVbgkVsKa5fVWDRBRoyd87rOZU
wD3P170vD0TgFSWida6tgJXvouZ9eBGgolnKnhmogEXfADv6GREgZos/3y1T
CU65nLMW4yIQwKE9Z3SuEhZmsu58GxKBPa7NsOZrlTBnOOJ6oVMECPx6PVV+
lTAldMbLG0RgU3LfMRdXCcM1m69n4kRgTt0y/I17FVjuvmJUjBYBBwzlKZrQ
KhiSE3kXECkCE+fyeoMTqmAg5kosb4gIjGnRCvnVVoFZ8xy303MRuHm5sIYw
UAV9RP4JBX4iMKxr5/RotQpMlTgEj3mLgJUBI8M6aTX03MlMNn4kAgU3L6+3
C1XDtcSLYkn3RUDzAE0UXamGro6hr8vOItD8UbXzc0A1dJyh+RlySwRG+yR+
PxivhjY6x4Jyw6P1GZN7MX2uAQPNQzVaXREIWj+4Q8DVQKtHdKmZlgiQRWyZ
jxPVQvNodcWmughwtUyp59jUQn2OULO0tAhgLlURXRGsg4prg2PRZCLw7Ix/
SMSnesCnxm74HQrDC5ErtM8760Fu+xqZw64wRNMIRzyiaAA7LRZ2wy1h+PJn
g8XxTAPEReMkVNaEIW+wKsb8XgPgZl6rCSwKQ0XVex79lAYgV9HVpZgRhvbv
9p8wvQ2g/oLSanVMGMT2Hj9WEWiEh701rr2DwrB5uRC13m+EdInn/tieo/Gx
O5T2lY0w4onefm0XhvBZ5fY9libQ5Sz57VknDNxBv/+TKmgCUbPsfY5fwqAJ
uc8CD5rB/NtdeqJMYWBk2NDhMWyBCIK04Fy6MAxZnWTJTWqB/Q+p5wsTheHx
XlbK6PlW6B6Ie3njpTAkKv+sPv2yDWiP34hHQcJgEbQU1j7YBpo+bJmS/sLA
thAlnceOg0yeiHb8I2GYVO3e0HHEQbBlEOeH/4TBr8fk+WOadij9oSkdYCkM
075us+cV22Fj//DMHVNh0BN7p8tg1Q6S+mX6V4yF4VdzbtZAUDvcTPC2VtMX
BoVRsgv+Ne3wflnlgZC2MLQctNz+TNEBjRpbz6jOC4MDb0xo5eUOUBm5l9av
erQfZlJtZK0dMGdtrsjNKAxGBg+n6Ua7gD+bQ4uERhgWXM5Qygt1w1WibtOF
Y8IQ9IrsuJFtN1QkGfoU7wnBz/ro+29nuyFu/EKl+awQUJ4v2WNf6wHcSZIu
zQkhSPkvkF9VoRfIn2GnpYeFQN1f/5zZ0TvgofBpqr0OIXhQNhr8Ed8L6Q93
uCdbhIBmKP1b6ek+GKnMk2mqF4LU3QfNw0/6gJXl4dlfVUKA4Tqz+re8D3Rt
5Y0+lh/Nn1PzxedvH4xl4Y8NFAlBBq3OKqVUP7jvV5Rw/haCPIdWtWiT/qP3
/VWJ6O9CUM3X25qd3g8rdV5/I2KFoCN0Dj9tPACRMtU/gh8KAeE/Gn2jr4Mg
7vX6Vs1dIcCXvokZbh+E4prrnGR3hICIk3XCaX8QJm7OPw+wFgKmVl7PoCtD
4PE9Vw3MhYBL6nMF85MhoMH7rP69JgScgWK0SWlDoBhJb/FETwgETst9Lt4b
gvrhXsaSS0IgHpW7eEl8GKykP9cSNIVAbE1ZudtoGDY8nJ6onj16rlvif8tn
GIKrTil4qgmBbBqmaTV1GLgZD2Z/KwqBEnE1my9uGDItaxM25Y/2x1Lbhnpv
GDS/vTE5dVwIzlbHnT5ONwK9WzeoHooLAcmwMm3QwjC4aApjs4WEoG6rY3ik
bhiI3iy6r/IKQSD9vSzVo/iSkn5jLixC8Bfz1XjJehi8zVMKHIiOxqdwfDDp
HgI+7LK91YAgHLuDFxOGQcilKeCJ7xaEpoD3Oz4Jg6B1w799ACcIL+JONHR7
D4LrOrPGjTpBIG12dAlVGgSSs0Mb0ZWCUDFFdnaSeRBiXqamd5cJQuBBEoPG
2gCAiCqrSZ4gUMj3565lDMDV+8SNkVmCUHPpUaBOyADMlzY+xWUIQqgN0/Uv
9gPgR/VeiSFNEAy8fkr81RwA5utWi/qfBYHm7WWCqeAApCWLfw6LF4TGjOnG
nIN+OL26er0p5uj3r3aUhtr6ofVMES31O0EYJopACdH9YBvyrFI7XBC8PCR1
b1r1Q/3gb6dPoYLAcb/nM0asH+Tk55m3AwWh8k4gXnC5D3a7DW1TfASB3Wos
dfIoXyNE+Mj+OgrCt/OxRHfie6G0vODir/OC0MpIme822A3CzIvLVJij+NT5
tNdSuuGFPX+0zWlB4CSzs1V27gZj2qAZupOCYLODZSTsdsGcmfELR35B+DPs
6ezH3QUGP4PksVyCkNArVmk90Qm/iIp62dgE4Ux7J+e5jE7gvLr01IVREAYb
A+4JP+wEv68CklU0ghBcLV9LeroTpnaNcVwUgqD+adySnawTLhsEP75PIgjd
7js+bWMdYHXe7Vz4ngA81KH7GFLWAZIL/ZFnNgSO6iiR4nNxHfAnAk3MzQlA
xpZq/65HB5SrfFWIHhUArSYD/C+TDrga5N651iwATz28lcXpOsBfkIEt7asA
rDW351L7t0Pf9QuxjFYC0OHl/73THgcpBxlzZSYC4GoY3RimiQPXL0xqzjoC
QCP2Y/6iAA7INkb6a1QEIDgpy/B+Thu0xV6UdpMTgHLdwdcTTm0Qh/nhLSgm
ANvb5M1XRdvgxGsvHh9mAXDUv6mtFtMKe6fGnCWpBCAJHxr83agVage0SrsP
+aEv5Xc1H3UrRAb8pH2+zQ+MhuMkEVUtYCnJanVimR90dmnPkfi2gESb98/h
SX4gvUqbf+NKC2w8Gj94OcAP/ZG96sqnWqCcV9tAtZ0fctuSq1jYWiC0KjNx
uo4fQulcdda3m6Hy83909WX8kMyqQYPpaYb9p6xPMn7xQzkPbXPY72ZwVfc2
u/eJH1alMvQl3ZvB+PdFVoIHP1zRmjGyXm+CV+/wzwZd+MFJ/zfTT1wTVD/4
vl52ix+eXw3s2M1qAmU5prZnBvyQd0v4arRrE3B/HQqlFecHjqc3TVsWG8Ek
6A1+hYcf5INlOXmaGiHMVtOhnYkftF/v9zlmNMJfgfQL0Qd84PUxzpzsTiOo
/DX/9fgPH7xNvsNjrN0I94fohC3m+SDjm+rQJ4lGyCiGCI1RPqjLpkhYJm+E
yVi3Q4FuPhgt6LE6M9MAvJ7iriRNfIAvT+UPrWmA69f6h6aAD9hjs7+m5zXA
m1NhunX5fKDgNr9kVdIA9UyY4m8/+MDAQFiBpbIBiNfWJcOS+cBZysKzvr4B
Trd+iXGN5YMXZO9Lfdsa4MdL6kenAvkAio9dmhtqgGnHskk2bz4YjDobljDZ
APyX7hvj7/HBzn3PduOFBogg7ZEvs+ADeYl5i7KdBnj05POC1ik++PSu+f5b
5kbINL9qJiXFB6Wux/IvcTXCnCpFPY0AH/RdPru3f1RfmG+5pOKo+YCBKDvI
QbYRzrqq2phP8MJTl3dx6kf746m32KrexwsJl5pH1w0aIVs6UUOglRcKhY+J
fb3WCPOUV34QV/NC14GGk6VlIwjPkvJMFfHCap9HFpNtI1jU5IfWZvECXV7W
Zu3R/r9PuYNPT+UFyTdzak+O6pnWAF6HVx95oUqu5sf+q0agsGnrvhvJC0S9
lLmpVY2Azj67YPSCFzBP9QoM9hrBm1fpl4IvL3hKRJTuKDTBUl9cxI4jL2x6
ctRdTm4C3JXbugGIFxY5xYeW/ZuhVvIkCZMKLwwmaTsM5zVDyeFeYZIsL7RJ
Om80zzVD2s9IceDmhV+q2ZQ/jFrgCVUZ8eEmD/jeOK3kJNgKD8deFIYv8oDz
mCWYXW0FhwLje/wTPGDj+FT3ckgrGNvPDWm08YDW4+r/JNdaQbKCpfDJNx5g
iDUIn8G2AX/siCttEg8cCDzg6vnTBiz3vonFR/MAofNyxwQLDii13IeOv+aB
PxyDMQOyOPjLi3lX8pwHFi1crDou4WDzD5WOrjcPHE5knAm0xcF8YxfR4H0e
uJfp7ceSgoORz58KnBx4YNRLpyJ5Agddj51cd614QOciN5mCcDs0GiqJvbzK
A3mMC1oV/7UDVpxoiEuXB0SHikKNPrdDRlfUZTVVHiByM2O4J9gBAScqRB8z
88CPpFjuH/yd4EERNkhBxQP8Lo5WZ252gsvI9bcxh9wQoqKa1JjQCTfClg7z
l7jBorVXbJ63C+TmOAa3armBaZ/9pDhPNwx8co108+GGJ6ZRl5I4ekGBgfev
6xQXuH1qEwlYGIBovlOdpwa4wPaurs/W8UEgHNdJx7dxwY0zdR137g4CaD82
flbCBWd7ywOM1wbB8FlXWtRbLiBnyBwV2R6C3IhFH/MQLtgfklKJVRkG9k8k
VwT8uGApIzWc1msYRkpO7KXf4YIWrQSNrf+rHxov4e5ac0EpK1fUHY0RSO27
mapwjQvyJt4vDfuNAOXsI+8dHS7IyGa8YAwj4LIVZliKuCDOL+xjLfEo4Ei/
iAYoc8FLPco/p8+PwinmEsJFGS5Ye+V/3ll1FKIFO1qphbngoqp386WKEdiV
m09p4+CCT1Nu10R1RuCmBpHXezouOIi4O0LUOQwVuhwGZqRcYKTh4DBsMQwh
jhfxEyucsB9l7hV9dwiWPCxbvk5xgpnmVRK37aP1B7kluwxwwvcV/VeGTweB
Pfmz3nYNJ5hfOpdAGTkAqYMHSVQJnFCNl6jwzusDygVWj9a3nMCaKqRjqtEH
Lvjjuu9COOHBFZ7OU7W9oMhmvsXrzglSGfTTS709UGmQf/mELie8tNiitN7r
AnGrFv4txAkzFKuRZwK7INR56k+RMieo5M1xc9J1gVEIc+J54aPntEMy7Uf5
khct5UbJyQmaRd2/f37tAI5UpN1Cxwnv7dvOvjzRAd6/TPneknKC9Dn+gZc9
7TBS4bpxncABh4IpSWFH9b4mLqiOZ5UD2okkHcKP+oPUkfj4sSkOSBj7IRsx
iQOXvUYtJxwHnEnKL353Dgd7ihIffb8ejT/7/pvSXgtwpo+e/3KNA0jVJ5iM
F+uBv7rxIFaXA1hYiuZ09epBdOx3wetzHCCw8AZ78WcdnOAMO+4pywFnYzXu
qt2rBWVFz2lnEQ7Qv8d64dRRv6hudOuTDRcHWGgtcsuerIFLIaosOsc4wHEz
tl7wTxXofRFpObvHDreb7n3ivloFxkD/4tQ6OzxK1vJg/V0JN4YISHKWHQK9
+PTp2SvhJn5ql3eYHV4VcywEs1WCHSsuj6mTHXji84j0RI/69RMlruQN7PDT
9woH06lK8HCMnFjNY4dW9OrCR6NKeBL4JH4qgx1uCkta2FhXwrMkh+v9Seyw
QlrzQMy1EsL7NBorw9iBqe5vYtarSkjTXsiJsmMHE/2HhKraSuiWPB9whpUd
hszMMpfUj/r1C3JnTlCzg+uZ7eoc3SoYs+HaEiViP7qH3g16mFfBUuyqI/0S
GwiMNVMSP66CjbwBYdIJNsiuuCNQG1wFeFzN0E4vG2imkCu/iqqCv0vZ0Ust
bNAZmKJn9KUKyKjijcar2MD2NrJl+1UF1GIvqHuK2KD9tm/hLH01MJ57WN2Y
xQb5Izkq5i+qgddLWzXvIxs8b+VWxjyuAcX5TYEXbmzAnbmkQGVeB6ePjfX7
3GEDInGhXJ+OOkBCTe/uW7PBTMK1kys69aBn9pnCXJcNcl+Xy3ce9d+2Dfqr
MiJscP1uhEwCXxO8/ZZW3tHBCm9lFMVPjrfC1fJv4N7AChZqyxlBV9qAvfNH
JTuwgohWmvxARRv0zWZVF+azwqNv8uYcXDj4uJ9ba/6TFbrkw+LPXcHBTab8
+v0UVjiVPzfiHIoDQfGixsQ4VhB1KJvYL8fBxOnSZhTJCg/drxayzODgiyG2
deIFK0DA4mtpuna4bVeJC/RjBfo3z2zPHeW7lFdNh/gjVjCN51K7YdEOmclN
PU63WGE1/9JU0Ld2KNjqHdlQZwXBpZTT4xeP/k/Ug2PvT7GCE+E0I96lA9QF
RiaUpVmhkLxjmv59x9H9OTXjzcEKJkJEkeoTR/3yx7UVkg0WeHbdcjbGrxPC
sv+sfZljgTbbP6WZqZ1gULO9oTXKAnwPXr6tae6EzpX97ZfNLPD7VaHGn6P6
fOgc1SHzVxaYAdb3BgVdkHidlvh3AgsotX6/Yz/SBf85M5CavmcBv0FNzJNj
3SDqz3yM8JIF6uf6Wd/JdMPMezaKjwEswLZ9f+GbSTekf+Ok0njMAraklADe
3eBUzkMz6soCdON6f5oMukGmk58uwJ4F3hgPrfzH1g2rs0IMIpYswF7tvLAz
2AW5+6JMNcYskKC0N/06uQvcmSRZHC6zAP/Xl+Mid7oAf1qO47syCwi9/NZn
uN0JJOGnBU4xscBJo9nar4ydUJ2sIdRNwQK/KjwqNXo7ILgAiXj+ZQbFU+Tl
nQkdQDOuJVGyyAwa7GK/D6U7gFXRRP58LTNoD91KMb3YDj3a109OlzJDg/5G
4jJNO3ywMjv14hczGGED4p534ID/hbVKUxIzmCUnRWZZ40BiwBlj4sMMdbjS
Q/rNVrghXmIt/YAZ6L9LLT73bIWQh9T+xA7McCkouge/2wKF5WZJfZbM8Mya
tPKuXwvMUX+DLGNmKFa7/3OC+Oi+M8WPBWszg1Dhkh5fVTNcTrlEfPMsM+Ca
C434XjZDxpmZc7TSzCC5bXSDj7MZBl4o3ZoUYIYOaj5LvtEmoOkKfFbMdhRf
YN6aL60JXFxEqhyJmWFAO8CBT7EJ5BJsLtT0MUHYxxwPPqNGuDmfZRffygRK
2b7efByN8EbpMNCtmgmmqy/78Y00wFpLQo1QNhNorIwH8bk0QM7fAS2/F0ww
j5ij+ELqYUJH2uGaLxPEXhuJ5TOsB+YYrxcybkyg6ZQRz8deDw/lOesHbjJB
1DvNL3xf6iDZx2Eu5yoTYNLp0/mc66CzLp8yVIcJVksHvvMp1AEZK7mUDWKC
hPa0LD5CLSjaXLusoswEl2Ye/uKDWrD78eUOvQwTiMf8Z60fVAtR+D+h00JM
4JfguDUhXwvbEW8b39EdPc/wEKIPqoGCijbdtSlGcKl662UnXw1zdAIudQOM
UNPwgWG3vwo4zV3DEnGMwIpLSn0TWAWPN2hb9EoZ4ddQZntxfyWcFtEx+PaO
EUi3G6WYAivA6d4HV/9QRjDeb8emyVXAx5K5cNOnjPCdpP+aej/AvsmLtmPO
R88ZZgMc5ABkknpWh2wY4cLxIfndHCxYLokx5l1nhEmVBQsxCiyUBVYb3dJk
BIUr5Hk+OWWwgmN5oKbKCBNWrGNp5GXAz2cbySjHCAFOwrQdFqVgcCc3Z1aE
EYQ8T6geZJeA32/ijnIuRoDnZ+0kyUsgk/jKRhQDI9hE6EWYWBTDqH4S891j
jED7/ZXR2GERCCiO2odsM8DJwki+lw3F4N/sm6Y8wQBXa2IWTr0rAfPrep6Y
NgbALLPLuRqWwqkxbm3tUoajuiL6QTpNGcz8KZg3i2YAVvv3eN7ActDn2ZDx
1WGA28+0SZ8XA3A73f5VlUsPIsIXmVy+Vx3Vu4qBzZ/ogTOf7FQoVTW0+JJe
6w6jBxad6qtpDtXgH5m0PXObHkgfno8dE6kBc557dasm9LB7jLT4IODofZGq
EYtH9LD2oXKQe6wGZgsH1Kh56GGs4pzAtfhaAM1v1CyU9NB/jfjcw91aiGv2
HOTZogPcPNx6c6MO3K9r/RCdoIMaX//AH/l1YDDG6ivbRgclTCitgbUeJJ0m
9ZVL6WDD5hSP9st6IN7M4cd8o4MP+Vc2eHkaII/CEIye08HIf+FJteqNEB7J
F2l2nw78Cn54fmxpBAeepf9uWdGBKF2Twf2bR/2sfCipuwodOBVSHHAd9SsB
16u0YhZpYY0+wOxuVSuYj0VyJPXRQqzdpxOaV9tA0clmLr2GFpSLyyg4pttg
zvcgtPgTLRQrHJjbPMeBQapKy7AJLRCPqFFRp7aDpDx54gyihXbKBnf35XYg
KepyXZWlhaRTN8ZGlDrgd/MDRmJKWrgQ6lHwq7YDwq+jcaotGuDIOybCz9AJ
jmP0ucwTNDAz8j48xLQTeDe/m4iW0kCw4i97q7lO2Pb1FpX9RgPXrDXb6050
QRvF5S2laBqQfNmuruDVBemRHLVnn9PAQZ5NenxFFzzjmYm+dJ8GRHZaL/gH
dIFlap6DkRUNrClKCXN0dIKy/HNVMx0aKH34/PCncCfMawoOOInSgN6KSvFg
VTsknnkxO8xAA1tu00UNx9rhquLKptEeNXwlvC0qOOqXQKyUQaWDGliphx5I
xLRALOWNi6T+1PBFRfteTNrReRJjjT2cqYHifIqrm3QtkO2K2cxfpwZng793
DTOr4d7iH+82WWoQss9zocivAK3WN9kfB6ngjNi5rpWXBbBfu11GX0d1VF/u
dKS9zYNcrFVTQC4VNBr/bL8ZlwN8OcdnHEKpwJHoVGKCyVfo+Pb2z4A7FWTJ
fvpIZ/sZQpIJRAY2VOCfz66wNxUNZz/a0FfoUoGKyIhZPkUwbL6r41FUoYJC
R2aKj+m6KCNMTuqr8NF81Wf+aof5I5ugKGVueip48fOr3BnqV4jdb//8awIl
+DSFZrFXv0VNHrZXiKYpwQuwdQyksSjZ7DC4+S0ldDbosno/TUM3jG+/x/hR
wsnq20nb/N8QvW7L59w7lPCmPEDGrfw78lL/WBqLKOF3v4Vwvl82MhJQ2bBb
oYCA1sjTWWW/ETln4mFvPwVk637yv0tXgEoZj9Hp1lDASP2P2uNWhUiCpFPi
ZDwF6FXVGafvF6O/U65WB5cpgHbI7zeOD1DecLfTfSUK+BIX1JdVA8ipR/3x
pCAFqJuF7YXfrUCCbSlB12kpoIPjHb8rWyXqrqN+17BDDg97PpzTL6tEL+FB
kvokOWTRtJrKWVYhVNT3M6uVHNjMpkxO1FWh7RxMiXAxOfik7RoqKFSjHxlp
9VGp5DD+h1FPMaEa3Uqh66GMJAftcxLaypQ1iDPefdLnCTlkhmtcUHWrQYGv
Nf/+Z3I0XspJXf1yLbK4P3JKjpgcuK2rZC5s1aHrgcvUxMvHgDfnma3JzXpk
FLs/3tl3DIRJNeNu1dWjC8Ad4ZV9DE5+BYqADw1IjtF0qdr6GGiulk2UaTQh
SdHbVTG6x+DyOV+ulrQmJKL6KM5J5RjovlM3GmJoRpw277QZGY6BiUpJ2e54
M2J2TxaY3CUD81DvTSrdFkQXkrP9e4YMLAfVjnPltSCS7LYvFuVkYPe0MFYl
uBXtV434yGWQgXO7Z5vWaiva6V02Jo4mg9PTXQFVZ9vQxuK+VFcAGYSURfo7
GrehpUMa4q93yaA3yuAp3e02NMPC0+dlRgbirjR+uV5taExCOkvvIhl4aNU/
MX3dhroNta02eMmAeeec9+e8NoSzNVWsoSSD/9r+Ptaqb0ONnrdpYjdJIetr
iefiYBsq/xRYpNFMChI/4orjiHEofRXYQn1IocaStXBKEoeSSXFLFg6ksFj/
UKVAAYcSOEar5ExIgUWpPT9UHYfeYg4edEmTwi268N9yRjgUZkJ7+Ss7KYR6
LSkSm+NQsAOPoDcJKeRM6+R12uKQzxu1VoEBEjgso/j12BOHPFK0UzdqSEBc
+raCbgAO3S8wfVKTQwL60dU5fK9wyKnptklsAgk8IhE5ufYeh+xGH0k7h5JA
vGtAdmUiDt38E0h89hEJVA+Mykel45AZxfs+xv9IYEnrbJZDLg5d5UnJmtQj
gbcKX1SuTuKQgXxucL4qCRxwS+sqvW5H2ucrrEJFScCRNPsmu0oHOm+KU7Rk
JIGGRaWHO2MdSPXpygTxHDFEl577UKTYhYTapDH6McSgetNsVmOoB226puxM
bhHBw88hz7WODyEH94mbTJNE4OlLl1zRMIQGvIRqz+KI4InZW1B3HEYVgZ/e
f8ggAj+m+AOF1BEUHvfxpKE10VEl1Fvy5s0YMgqT6HWPPMTWC6zT6cxNIOWY
1NQ250OsOTEXU4zsJOJLFnGX1jrELkwg1qmHk2ihQIBpbPcvliYtgtvvYBIF
TnLo6tr9xerLnJDIZZ5GhacpsUIqB9h21XuIW2MWJV4Mfv2E6QBryx1z3uHZ
LAoyIrPsXdzHbuyVa+XVzSLj20SEsKR9LCuWXt/AeA6pPvCrm/Xex6YlKRt9
jJlDAk/2ozSv7WOVnt00mRuaQ0uReAU81T72+sVMs2cO86gj3oPYZGoPOyfe
Y9n2Yx4Vfd1s+1m+h/Wk/GvNuzGPknIfJlJ+2MOSL4jZ3lFZQMFlay62bntY
LcoI0snoBeRS73qmXH8Pe1fPUfXBrUVk0rlExSW5h30dge4eyi4hwfm5tNah
XWwH53oPT80yIt+8/UiqYBeLt2ygyYhYQct/p84HRu5iOT9/RqqWq6iYdXxM
TWsXm8kqdS3Daw1dRf3cqVkE7GZ8f3kt4wZ6EdMQ7vMMj3328/Tbw7NbSNz/
cU+GGR5rYetRcOPnFqp2lOAfOIHHKnPlDuXwbCOS00E/VUd3sNOBUhJ2+G2U
JKy45ZC/g4XT9nplt3cQhmZSPeb1DjZ2LekBe/cOejKEmrfO7GB1LTlL63Pw
iKdmlUWMZQcrwnx1XEiQgIp/JlpcXdjG7tW9IfcJJ6Ab0fopzyu2sY00yrmE
GALa8dtfyI3dxuYJ4k+bRhNQtMN3hYl729hEpeKqvPcEpGhk7s10aRsbovNE
j/kdAT0UKqS6v7WFtXQnutnyhoCYqB2ufGrewmqFVs5IH8XP3mD70JqyhT2R
GHgvJIyAlqvcJGWNt7Bk9ZQB50MJyPG2/KWFzE3sdx6mL0XPCIjCcCScK3gT
G32iU4YjgIDSVF73aFttYp9ejPrt/pSAJikX7b/SbGKv3uOqk39CQDe/pwXa
Of7B7lcKLaQ+IqCDd9ea32H+YKd7Jx+SuhNQ/BMy1ir2P9i2pdQ9m4cEdNr+
l8XG8ga2iNgxsPw+AfXr30oRqtnAprBL0/HeI6DHykyLRvEb2NfHl6K87hIQ
uwAo+LttYD1RJn+vMwH9prjnnaWzgV3nj35jeYeArq7xVY4IbWCjVBTipm8R
0J++Zip6wjpWzaj1i6slAb2t8LmigVvHDjk6Ze1cI6C2t/1jH/3WsSJxydVU
lwlIRynq8MbgGjZnn3lVVp6A5vkuXArJXcNeZcsi5EsSUAj5n/CC0DUsXlaX
DAkTUE2vET+76hr2rPVzLmM2AvLf5PROk1rFNlf8Oe+xh0flNR5Ux/pWsPcH
3hgcbuLRXnRXrG3wCpb1z3GzkBU88lCLKBCcXMZaiNq5xo3jkbMfxVZc/BJ2
PrgrtrwOj8LYxGpsjZaw0k27zr8r8OjHD80oGbIlrAuDEOZHCR61nLex3ypY
xGaZXGJO+Y1HKwO+SuXOi9i7yY8IC0f5xvDw47EXAovYqWCnqxkFeHSCqqjb
sHMBa+5snXWnDI8eqmw+Gj+9gNVWvGw/V49H71qZtDJW5rHlnGfhayse5dnL
s7slz2OVDxR4HLrwaPudUz4Z9TxWpJa3fXoUj7zWxrYl+mexf03X1Ce38Sju
xd/atbBZrLv6dEzyPh4V8/PGFKFZ7KLgwMZ/JAS0r2eqops+g+2br0ofozva
r28tnq6e09hfPjGsI6IE9BktXlI5Po09bhN2L0GagCp7KTmJRqewSRcCGi1P
EBAZ+YXCSK0p7Gs6Z//BMwQklvBfiPnuJPbYhvVg3DkC0lJ8ekMkcxLr03NV
2fwSAb24VYzPY5vE3vl0dqnPmIDSCb31vg0T2LHnpy7F3iCg+oitWC3fCayp
o2Sy6c2j85ZgucNwcgK7vcVj2MVPQNTYE2p9U+NYpoAZq3vP8ej4dQOqz7Hj
WBnaHBfquR2kt+zcf0dvHHspxscnVW8HhfN89dr7NYZ9+pPxwyDbNsrKrb5c
7TCG/ag6+NXDewvhLk9wveYZwxZUpeYzjW4ikbXT7hrko9il/tNdl9L/IOlv
ElvPvwxhzcjtGfJOryOX3ltRbLRD2EeR8vyGn9ZQ5rFEpa9ug9gI3l2ZBdI1
dPIWq0ej5gD2lT1/FIpbRqo8xDuM471Y+/KHNFzYOeR9WT0mWbsXm9BFIWUY
MotKPT1VFLN7sJ0L8VpBV2YQpmvZ0/RpNxbDURewPjmJ/EmkOOdnurDuspaf
JH5OoMoTdoXeBl3YH+fXS608xpHW6wF8Al8nVmFMee0d0wj6EXV9RjakAztw
oL2tmD2IWBI7Oss22rE2hl7Csz39yDvNoELfqh0rOnzF2l63D41lNmYO1+Gw
k07S8ZPYHnSpQCv+rgIOK7tGWZwT2I0sngax7Km1YXcrsrR84rqQZITvCjdD
K7bu3fXOC9mdaDvpUcPpqWbsT46TD06QdqLqnLtfzIuasGrS5ykJlR0ostL+
qXd4I7ZG/dqnimcd6GanlXmcbQOWEMRaQGLagWSmrikVq9ZjU+0nlXXdOhBh
U59xgK4Oe+Vibv67Nx2o9pjWImGiBhv0KsXd9XsHes9+tparsBqbsZQYGnYU
/z8J5c9qr6uwOP24xIz+DiSnKvfE7FYldjsz6lf9Wgfa1xY39VKpwJYeJIhb
HP2+wYxf4QMtYLffRioNHsWX6ZbhflhVhpWXDLpgYdWB6rZ5cTl3S7CrmLr/
ak51oFucdMHrHEVY9vMX7pvQHsVTOzhzsjIfq6FV4Tc+1Y5iLJbX77vkYfXO
3f0s+qodKfgOf81mz8U6DJiuOXG2o1iNjyVjAVnYZ26amJxUHGr195Gw6cnA
ttYIXWAJbUOk1RbvRqXTsKf3rtsWkbQgVQr1w5tPP2NTT4Q9s37SiAbfq/aU
fv2ALVEonjPvqUFaPISfJUYRWOEq+v7rtSUIgmfYF2Z8sF/OGVesYr8iRjOR
9NtnOzA3z0QaiH0agPr+ye7AsQ6M2TcxuXucwxBg/oU05Vknhu9R2MGzn6Ow
YSFmPVrbhTEilzXodp6AzpsS7DcMezHtCbtL+1GzEDY6e95jvRfzizSzK3hs
Di7YpD94/64P8/bOf6WMxxcg7z+pFlxvP8ZEqf6VKCxBtN3xwMs2QxjnGR02
8a01MJhezHEgGcZYfBClco5bB/LbP0aDvgxjdPX+7mdhNsDTQfZM5dwIRiY3
Z0rt5R84MbfsOPZyFMNn/6rPT34T3BVNdZ+yj2HoOO2bq7o2QeY3h1dC4Rhm
yZczT19gG0K2dt78oBrHDJ/Y+Pq2ehumlPrSSszHMc2TTR977+wA8igsa/w+
jimLTn3Dy4CHhPzYrv79cUzm5afP/8vDw+7248U5/QlM0v4NzzQzAlxXMSPB
J05gFjup5VssCJDrqcZFsTaBcS8NM8m+RgD6Qq4T7OcmMftf6B+/MySAE56g
JfZ2EhP0OiLe4zIBalUHrBQnJzH0HswVZucJ4FcU98o4aAojcImdWlCFAKvF
PLPhDdOYc/N88T/4CaC3t3eQwD2DaWj/BG84CZB+Zoj1p/MM5kqx0PRDZgLY
lMZrNtHNYm69EpNTIycArowvkcJ4DhMkKws1y3iQ+XuQx54yh6Fnz5pKn8VD
yNmRJrHNOUzU3xNUYeN4OIdNJJyPnsektSkaG3fjIfHQj9Fkdh4jX5jvodR2
tH6MtcQt1QVMQZLqR84GPOSCwDX/gQVMwwP1qZFSPNATE7m8Ob6IMTYvp6zM
x4PTubFniU8WMRJpItRr2XioDYAPP1sWMdN67o4WBXgQrkzKLuVfwqT8qa6t
LcODH0lAXdO9Jcx/cWxiCtV4GND8b2QAljCC524/T2jEw7sqIVrCrWVMfDgF
cu/FA+1ZzYmEzGVMKVVr94NhPAQW3irU3F3GDAW+d743iQePnyl2YZErGF4f
oRinVTzciBErE4BVzEeHMys2FASYo0y5PHFsDVMyQRx4k44AYlq/8vKOr2EG
req5LFkIkAidL8081zDcV69dNBUgwJvfrCop9OuYj+hevP7R+TRviCY/UlzH
FJcoKehqEIDyhBKdtvk6ZkB5v0776Dz9M65NLqWuY7hkQ/+cPzr/0ll7A2zT
OkYt3TD03FF+EEQ9iiLX1zFmIuwCmKP8Ub4VLGrHsYHxShzKU/+PAA8/Rb9R
1tjAfOBK0TntQIDMobRdStsNTFvdA305ZwIschXYD4ZsYBCWljXXhQCSpnW4
n5kbmNz8rwNKrgSwe997xr9rAyOaqfm56B4Bktpn04x3NzDRqcMOGg8IMESP
ZxIT/IOhSHgsV/GQANdCOecanP9gFl9lltQ/IsDbWkmT+Mg/GKvnl5/peRKg
lVSt3LXgD6bNe0ob95gAWn5m71lINzG/7nD39PoQQNUtDlnFbWK8LwWvzQQQ
wD0747s8bGIWzwoVOD0nQM5yMTvJzCbGSrnUdyWQANKOg0tpJ7cw58Q2qLdf
EIDHivfDWv0WhpL0ptixcALciJMhq1rZwnjt4hdD3hDgfa/6vSjWbczS+rtc
2sij/GbVH3A4vY2xmpfzinhLAFpjq4unbbYxbWMNGNb3BNB+czebNngbg/rs
yGOjCBDU7Msz+n0b86vtsJknhgCVVOHBOe3bGNG6uHefYgnwVytx/fnONqa7
r7apy4kApwMzLa/z7WCi1HRc3pERwLOivE7y/A7GJK6F1jgRD2tnRhNawncw
OIsu/bb2HRg/IajDIY7HtI6MtzWobsPtiMMzbIZ4TOX3KPvthi1YWB2RZXmM
xxQ/vrwnbLEFG5kJTAyNeEwGc474E99NIJHl6Sd3JWBeafn7ylVuQHDYbiNZ
DAHzjEWRxcJkA6iX+ktJgIDxHptNfzG1DizfY5P+Mu1iHLwNu8bI10FUit1x
+/cuJtRS47p4zAqkh2yZbY7sYry3vwQ8LVoGmbku3Q2KPYxjBO3PvsElyL6U
p7F2Yg9jJu3ed/LvIih+fSe/YraH0akeJH0luAiF5G5CS8/2MGwFIkV8Ngug
ftuYZeH7HsY4IS74K+sCQM3JY3Nde5iI50xXTzbMw0Uxpp3pgz1M050QoRLf
eWgMXJubFN/HkBsdrlxQmAeDqbaBccN9jKayR2nrzBxcTwkvH/68jyknsTOd
MJyD20K0wd0XDzA5HzGvqYRnYcF/0bPT9QCzGpBv/q5nBlzHGu+0xxxgjjvK
SvK9mgGPT6H6LfMHmDRFnqoTm9Pwgo+Creb1X8yHlm28ae0UpHORfMnrPsQ8
Pfxpl+Y0AVuU+2MNCsToDVm3Qoj0ADxMOM84r0WMht8INfIG9MPayVcYSgti
JMrj+l9Obx8smHMnaj0nRpUnKSIGAnvBcfW/FvsYYkRZZiJxb7AHpp9/2w/8
ToxMtZPKSRV6YOyHqkVVJzFavqm2LDPaBTfP+b+cmCVGpxaCAiuUumCwu66I
eJ8YPXvUwXM9rBNuODHMCzKSoHYigV8LEx3Q8/c6JxIlQRxhzjpP1Trg6rvE
S9aqJMiKo3CcJaId2iVmPPz0SNDPZDKv9BkcGJTKpiXYkKDPGsVVn3Jx0Gz0
qLvUnQTJylIyd6XjQGe6lGwohAQV8163ofyEgzovMsW9eBKkRfslUz0KB5XJ
796q1ZAg60WMbmoADgo3cwzz2UjRm2+egR7mOEiJXWnIyidFG2Kz25y0OBCU
VSa0NpKip2xKF/VJcJBY4Su5MkKK2KnJJUjX2yB2gfaFDAUZulSvYdPV0Abs
T6/m6/KQIeuv/fvq+W3wjiV+2kmeDD0OfvQhNbkNwtWPX/hmSoa+XfjZ4eHT
BrTtD93qnclQpcjle6MObRBqX5w8+5QMDZJM02hfbQPyXeIO8vdkaHPcPz0b
tUFg+GVi8XQyxMrhujNEaAUSkcgTF0vJECItCmKPbIWnBX3Wdjgy9HCVlNVI
shW8xxzKk/FkqL3uw4mq6y2Ad89arqA9hg5+TZXvLzfDI6od3nHBY0g+SV5f
OagZ7isE+whcPobCH9fc+farCW4HflGN/3AM6UqtJ79hbIQrkuO5UafJUeNL
bTbLxlq43oS/Yq5PjlLTvJSSHWvB3JVxnd+GHAVXZlybI68FuzyMfHoQOdLY
pY12P18DjqamrS6x5IiP/Wx+8Xg13N11vXvyOznaOXmvh8i/+ig/g2m3y8hR
h37SthZ/NXigxO/FOHL06047++vSKvCe/K3zdJIcvQ8iUek0r4KnwS3z57fJ
kcV8sHZpQyUESk2HUFJRIKnP5WUuJyohpHlfooWHAu3f2Fbgi6mAuBsWX0Zv
UyD/+xkM5bblkLj38PwXLwokPca1y3+yDFISX044hlGgKsPQqad/S+DHVJHg
Rg4FIpd3LDr3oQiw9zkTyA4okM+ymG1Nax5MhXRHSb2jRKtdmnc/uKXB3PFl
xZVUSnQhT0jTgv4LLLeSdeUWUiKPd0QcfN8+wzarIrP6CCXCpzwJOFf2AXYL
dHOI1inRp1unj9dsvIdDc1ujGlIqRCqE79KWjACyv96roexUqLfc+SC2JBgo
P78NN5CiQnqPWz8jfW+gvZAhy6JOhWx9aASK1awR02xFc68BFbpbGgszH0IR
t+w69X/uVOh6d4BrO+Yzki0yfOGGpUJUkfGPkNcvpGDpIK7aQYW+e5+yj5bP
R8qHfjX7U1RIwq7x6tJ0IUIXM0mDqamRowr+VKxxOTJup/WPvUaN2lekCDtR
1ei6u6iAlSM1emO7tSa8XoMsONTLhXyokXUfzOnr1SE7K+f9jCRqtFJh2pdC
0og85hsely1SI5qooALDe20oke/yIAcRDTpUv6Y044ZDdVfqNR6y0qCArNlw
0RPtaC3oUlKzJA0iFfaeu7XUjjiLa0kkNGiQ73tazc/pHQitXLQPuEKD8OSf
Po7adaI7wjV1g/Y0yNvr5BafUBd6e/2CtLI3DdpfrDKwHO5CxS+rwiLCaZCU
qZqkXXQXmijXXF1IpkFdD3ie/NTuQjR/Kq5cLKBBPmEHbdu7ncjSAth3R2kQ
rqLCM9S6E/09XpZke5IWecjo1GVKdiBxG3XS8ou0SOiSDA9+oB0Zvi+x5zSn
RQ3/0d8797odJe0XSbc8o0W8MR3sXes4pNmSn6fcRYuCSfObNx61ISdiZY7I
OVo0rOrtHjTQit4p/fZa3KdFMq5neTkxrWgq4dfZJDE61NlX66xB2YKCXLPr
qTzpkGRmL3VIXCPKTJY/bveKDj2e/JjLfdiAenoyX5d/okNtnDbmP20b0F9q
uTXOPDokYXB0ydbXI3HMT2O3ejp0T5Lhy3HlemTgJvO7ZYgOjQ8ojR9vr0Oe
X79zSK7TIf3XlvwyLnWoniFjaIiLHslsfIuVSa5Fa+elMCpy9CjxC65bRqMW
cT1O/xypSY+oTXeYZftqkNNY2m0tZ3o0W3LhtSx9DaLNSVnLKqVHX4JGA+S8
q5DitLAJdTs94lAlx8qxVSErruTfdtP06PmCzL5cdiXKDEjy5qJnQP8ZenvI
z1QgI+MEsufWDIifi835hDEc3aO8DsNuDEggxtPovwEs+lz8sUElhAFt71U3
/seIRRvCceFL2Qwor8Yae+tJGeI25VrXqmFAMdI/1WxzS9H5V7Emn/sZ0IM3
u79s50qQC5Yjf2+ZAV3evCRnx1+CGr10aY8TMyJZs6h0u6vFSJjHbSuDmhFJ
UtzmTqYsRt6lcSPSrIzI/Hf+qFV6MeqwqqzL4GNE4bYUqVzaJUjqcD5bWoIR
ZTHdcOqeLUH+SUxxGScYURs2XT7yRSnqO6f2XPo0IyJ1JWzqSZShkMCQa9L6
jMipMe5p9e1yNCaWfTbjOiOKf7x4wZ8ci1TreiWkbY7Gi6tTq6dh0RyVxK6U
GyOKq3Bt7H8I6NKbmkSpD4wIz8d47UpuBUo8sfziWzIjonxUqCJPUom221kf
SP1gRMLNNtx0xpUojdXuvBSWEV3zzhmtX69E5HEkc5LTjKi4x9gJnaxGN89I
t6evMKJO2T09/oBqlD90pVgSz4iWA1Pk93HVyF7gc5gkDRPiOLW5WXC/BpVC
/aN0ViYk9TK+NwpqEOuttZuS/ExIc/xCsRtjLXIh5dROl2BClqrL8VdsalH1
F8xJyZNMaF+PLqskpRbxajlwp59mQgoTugUU/bXIfTacVPICE3L0fIk1oa9D
ItIj3RKmTKj1MwVu8XEdepb1xFvClwklbBOv8xjUo3PFRVXiDUfzJW0rPME3
oFiL8R9pHUzogZLSmXrZRrS2TxktPsSEvja6nWe1bURJmBuO4qtMiHlrzeRH
SyPCjz29koZnQlov5S3xJE3I6NnX0+IkzMhXwNXugmoTIqrZphVnY0azlxfd
h5Kb0I3b/Nup/MyId1TKV7KvCWVTaI2KSTKjK+6OQe50zYjy29361JPMKIjq
62vQbEbWOlE5YmeYUUHidBTt42aUv1gal3qBGc2cEk288bMZ0b+eei5mcBSv
4Vbal4lmdFuO9m6qKTPKcYrnyjrbgsrbTl0X+48ZkasfLAv/akHuti/MHzkx
ozD50lRdsVYktTNws8aNGZGJ+Fi5x7Si9/zPHOyDmREJNaG5xufofrmLe5z6
nRkRuexeU0zGIUFiUd+tPGZUOd2foz2IQ93vPQMuljMjP+siOivWdoRK+V9O
45jRgfHjmqDgdsROc/ej2DYz2lPbVux16kDNifDp0SEzKv/VHbGY0oECFFi/
1FCyIB/Z30uHQx1oyazkhz0PC8ILun+RNOxEycv0Ob9FWVBZnAmRRkgnMg24
9fuYHAvyYj1leaWiE1WmU5WnIha0Q7HB4n2qCz1Wt6rcusyCigPa74W7dCFZ
XHbtRRMW5LmX3ZScenRf25I1RVmyoFOPIiQKRrpQ7I5p27Q9C9peufe8iaMb
Gbz63ql0jwUVORqOjhp1I1KBw96gxywow7V5McKrGxXmGg91B7CgMPU2hzzj
buSqlTYm9ooF3aHumOw93o1EBwhTj96zIO2+Lus9sm7Uf1d/viaBBUmn9Q7y
H813IWpz3T6bBS2eG+60j+hC3FOapGndLMjFa67qz2Ynqg0YFRMXZEUuc39T
GMo70BO2U9IeUqxIJ59Y8FRMBzr5LViuVoEVyQSSxV+/34HicbIqty+yomUB
6vcJwh3oocAT7TQnVuRqyhZ4/EU7kvzVprftxop0xTiJDG3a0bCWyBUtX1Yk
+4f7yUO1dvR2wONadDArYqzgw0cxtyNt10azmTesaCVc0L1oEYf+EvPfVP7A
iposRdaGqnHoV9SDW8HJrKhshbHixFsccpSqud3znRWJ8yrrmJvjEH8Zp7P4
b1YUetmi87kwDt3mErxyDsuKZo/71Pc3taGfjyRULBtYUU01ZuG7cRtSl1Mm
fTvMijy262WNbFqR+SeTjwfH2NCg2WBerE8zSt4zD+BkZENZG596nEia0bzp
LYdT3Gwo/JUdXj20CT1mfHDqjhwbQmXLZ8ajG1FUQHhT13U2lCNIVCWVW4+G
hqNyVm3YUERR9dSeaj0SPZ0QQ+3MhmyNQ8lbsXUod+O77VFxjoQCmXUettSi
XYNfly1C2dAf3v/XxJ2GQ/XHbQCPFmsl2Q1mM2mzFEWSO0qKUhLZKkSytSmS
UNZEWQpFRFJZS3aVMWObsWTLvlYkpI1W8cyL5/nP8/Jzzfmd8zv3da753ufN
6XTdYVaD7VllKhdixdFWmBgp0VuNZgfW36xH4gj5QG0pG6uCFL15uOaZOCwD
Pn6LOFWFo4Su2uFycWhK54od/cFEhvdgzmy1OARieXVn85iYah+NlWwWx/Uw
d2ViKhOb1KYubugRh5Bfl9yOGCYuR84c3fuec3/n9Jc6BzIhbLBoXeAvcUQf
kZ7Mc2TCLE1INJlXAqJmQb1t5kwkzon+KhGWwK3dn9k/dzGxpphUNUWUQLx6
9RPdNUwYKunaWO2WwD3BMydeTzEQFWSgd/6gBBQW9Jl/H2Sgc2ivUpStBFJn
DAwkWxhwumMzXX1GAulDBMWjzxkIFPSJULsrgZyiurdTXgywna6cNU6XgHL2
xlbRkwyIMsMOn8iVQF5qcuUmKwYsFaK2XS2RgFq8wDMrIwbuX4qn3mNIoCDC
877fVgbGOpMFSxokoHF18GbaegZU1DO+tHZIoNhrj3+NPAMXonI6Pg1JYIt7
ocf4cgZeTRa84J+QgKCg9+9TBznzZveLNMqMBA7WBAe5f2dg70Nm2LZ5CURe
jV3uGstEn22rmedKSYj8zVV0aqvCn6bPk4PakhCfHDlgc7AWwRsGw/fulITa
8RzLtLFaiMQ1KZXvk8T+Pk+7scucfmObczzOThLXGhae8XzMgumES59xqCT+
ZpOiIv7Vo3+flU9ZlCRkFT/Gt8Y0wDl/t5TSXUlo3HuaIqXUCL+LSmYLcyRx
KlI3L920CZl8o/WlrZL4XqR8YuBgMzRc35xc1ScJnnOirleSm0FvquK7PSKJ
5So/PKgfm2G0oeAh72dJKEx0n63d2IKO2w/0T/+SxLpHLy+4+LXA7nfMcD+P
FLQdUn2WslowaXPV30hICrsVgv2ermyFF/2MXKmYFMzLPV+Qr7SCl2pXTpPn
/P6InDI51YrI0P2Wt1Zx1sc2XymyaYPUhO5PHjUprPP3Ox7AboNyvvzGfn0p
iFj0rBbLaIcV+jJjT0qhY0lMgbZtB0Ye1BvynJMC65tu/OL6DpzmKx/18JVC
+cCni681OxHSdIe856YUkot26zqIdSHfxuLugkIpODnxsMPru6FDN9D0eCUF
C9O8bDOtHtRRNnX01krBaJvtTflHPRgYFxMt7pGCskSpWf7VXghebL3mvkAa
09WnB7u1+nG7r3JVr4A0Rp/JMx486gcRz6oNV0qj+15Durv4ALIe3HcoIkij
/ppP6ObAAWzii+Kh0qTx4rySC8/XAVS6+KdEq0gj267DuP7IIIybPHTmNaWh
t3tRtuHGQXSqHel105MGb32s9c8fA7C/vfdij5E0GMYUoYyyAXz6tVXS8JA0
rjbll5n5DcDbZl1h4RFpbN+v57JQbwA3KEJfo85wjj9oxzrG6kf5szrnXfGc
61ulhMp29aGr33XMP1UagT4e+iVxfZgWWO5ckiUN07s682aH+qBsf+jEaro0
3nX3XrjR2os00beOQh+lschSyomnsQfh5/7av94igz0W0Ts+VHbhUUryW76d
MhD1OrYgKKALVfXb7WEig644lRdE3S7MUcLsnjnIwHImo+BecSfOtosdi42Q
wSXvthcm6R2w0lC2tRiQwQu/dY3PL7fDy66l7+aYDIJSNpsMabbjVqSnTd03
GRjT9VqEOfOpaaTMWotPFn08lm8cXdqgF7fLSlZVFrNBwf2SFq04yhjvPrhF
FtUPo47uEOXMr6lIy4gdsoisSRw+3diCIoP2w/8Oy0KeP3+Epd+CtrNeXRoO
shhVennix1wzviTLHPZwl0XO7rqP5LJmCNe/7MzwkoWvS5uryflmTl85ZjF4
RRZf5/MSvhg3w4CyqFMyQhYD3wlzLarNsDd5ZL4/ThYNH645PBdrRtKjT4cq
M2WxK6XWZrLkNUQO+5g9aJNF+jaVsJ+dTfhVEnvghCIBxy7PDSyVaMThSflr
35QJCAn+TQkfbkCpQib9siYBGZHTzks4fcsnhK5yy4iAyXvj3+b0G9BTanRC
4RABKzNGNl9a3gDtTx3JmUcI2JA75Puzpx6zBz8trTxDwMVXHUu+nKmHTaj3
TuNLBKTUtBi56dTjZdnCy51BBFQ0NUSN8ddDfupGgf0NAoY6at8cb2fDjyQz
+SmegCWDDJnhFDYGzB5SLqYSYKibw6t/jg3dMFXrRVkEfFDrYU+5cfppeXnM
zQLO/ql8sXed2Fjw2YAt84oARUl1a4NjbBwjt/Jk1BJQLWBH+WbJRuUhWy21
FgIcZyMn7h1k4+qLc493vSfgwfCHSzMGbCwKvx/iKSgHYnZiUjqnD/Mr/tFr
MpcDYfvCAd6/LDgfDvaxPCaH8o2qGbnTLLCui+S/PykHK5qth9UUC+HfaKS/
vnJIECqeezbMwrhi/uGQEDlozr2rse1lwchSJ2pFlBw6vojcFHjDgjD94Pyq
B3IQ73BRsKtjwe37wKbn2XIorIv/IMxgoZHm4rGtSA5m5VV5JeUsKFvNPGRV
yOF7zlev44Us3IgM6DdjyaH6Q+ADAp2FKbqQ+FCrHAJvH+iYGWDBZDrO2LVP
Djv1FQRe/2NBxDqn/MpnOdQml3lc0Wbj9A3N78K/5RBsHJZqZcVGSyVzTQKv
PAz+HGrfeJGNGKWeO7ni8qg99FVrtIgNsR9LvHu15VEnrLTYaEM9pKOPqWuE
yaPR53ai2OoGdBVJhkZFy6PPsnredFcD4vqauifuyuOD5oxDtGMDRJV0/NNy
5DH7w2zdsrQGNO+dbvlXJA++jqAoY3oDpwdkUS3p8hAtLJgOH2iAYIU0e3mb
PFafE3vJJ9uIuvfNBNc+eaib7iAZaDUiRDDsVM2IPLapeQYHWTRih6oug/RZ
Ht1MMaei641YaP5D7PIveRhP8wQHtzaCcSnnRBePAiqonx+YSTchIPV42UYh
Baw51MegcPrY7GTr0XE5BSwoKlrA+NSEMtHw/J2rFOA6+kAhWv01vDW3L05V
VUCnRPS2Y5z3hZnAvEwLfQX0lmts81/UjCmZiJkqZwUUEkR85Th9cchw5/WU
5wpwGaj6aZnTgmSP2f4/LxSQckE7tLK5BTa3nqua1yigfVm+xOrpFnQPkjqE
uxWwQzdF49eWViQs7lZyfqsAv07x6qNHWmG+NuoSc0IBz09FmNVy5l3rhTmi
z5wCFO5fPBdX1wrDpp7S8YVEmGl+4Z2baIWS3SaHBH4izvrGTbxd2Qa+6Whh
g6VEmPuWxWgYt2E05FPh9xVEaPkOaIUFtaFaevfRVAki5Hx5h3tetCE9O53f
RJYIHl9a2PqZNti3WltnUomou+TR0erYju2OxYsOryYi+1LMZcXkdhB/ieYu
ViYi+lIR1bujHQNy7AV2m4mwuDR3Vm7XG1id1MyQ3EPEiM+N5c5KHTi4QPxr
0CkiRHxyDf+Nd8LXa49QgicRG09/URZX6kL6lL9i1kUi9jhtEF/v2IXpvnHL
lkAiLpkWDdsMdoFgRvR8H0ZEjOGvWk9CN3bUH7rxM5KIzG1bciMsu+Gmd/2x
YCwRleq+t9LjunGrlM6QSyCia80rnxdt3Xih+qNP9R4Rktpvw47w9uD9o7U/
9dOIMI2MGDi8qgfqcXFrXbKJYKsNX9t7pgftZ91uFdGJKPMJH9q0uBezH1Nz
WdVE/GxQ36S2phdUu866PjYRGxQGr6816YWnid4/3jdEZDE3bibe6YXYOikn
k49E3BXujxRY34etD/YG2E8R0XEk9N1C0z4clwm8e/47ESueqWnNXehDAf9U
U+I/IsIPhbz/Ru/DoRHm5rEVJFxJVtHuM+uHn83vA38lSKj+0h3VcbEfGW3K
bssIJCzQDxptTu5H057jISQSCVtvK2+tZ/bjR+Wd++o0Erw+dEVXj/VDXut1
2a61JBRoBX6oWDoAg6eL3lipkvClOT/T62c/PFZt+eyuwVlf64Fl7H7EJZ8S
uLKFhJCXazrSk/rxSvwh5ZYuCa+fj7pqn+rHaESPzqMdJEhnpvG0bu/Hpks7
zzbuIyEvTmb9/GgfOixke6UdSOg4JRhkqdaHuab9M+ucSSA61Uh/WdgHmkHI
criT4GJzNS+koxcX1L/qO3mRMG/4p+f5pV5IrKjNzr9OAo08rraspgfbwmZr
qqNIOCuVUZue0AOnebXhrtsklC+zt9V27UHRpyTx+RQS9v/tDnMW6YEF+6yf
UQEJ7u2sQYZVN4g7k6b5SkkIXneefUy1G+MV1S5VL0mIDSIWzi3uRsGWz0MB
DBIUp5qDw5q64FcoZaFTy9m/9fuvgtFdMFTVa/xdz8mD9dM28mAXRLNc9Yua
SdDcLMReLtGFjPuvVFR6SLgjuiFNLLET0xcEh0UXkLE1oJFXKe8NXn3daPF6
ERmfPg15ZJ55gzA328brAmQkW0/3rFN/A4L909JFK8mY3yT7XK20Hfp7zWN+
KJLR8snZQZvZBuE6f/7na8i4ae37+mVwGzr0nvidUiFDl3VTG4ZtcNGcdRnT
JCMtvWjlzsZWqD9XHE7XIeOgKNu/9mYr5tabWNjpkbEwoH9it2krYihp+j1G
ZDhaL6ra18n5f0quL43bT4YES1K15U4LaNIzKgcPcc6/xiGqyr0FX2LlM5Zb
kbG9IJr5ZlcLypYZEhqOkKG/jf5jlNSCoGtnYsIcOK6bWv3rbzP2LUrk3+lM
xk5TOVuBjmYM//40/fIMGQecfJhrw5uh2qyVkRpPRpVbf6lLYRP+7HEgHEki
ozZgCzHyZBOqqiNiZFLJYN+KD8mTa8Lh8kG/2EwyXr84cHA6uBH+GcEWwS/J
+ChcM3nZvAGGpLzG7Qwy/hEpB+8LcuZPUpf+XA0nT/WAUsaremREr1X1aiZj
ie2WEL5V9TgtfChj4xsy+M/ET67uZWNLqB/hSzcZgsHTpsY3Of2A93FM9gAZ
wncOlHros9Hk28J/8h0ZIjm5CtE/WRjubf8cPkbGPevpbYkvWZje0tWR/YmM
Q9k0fv1rLEj/Hkj/8pPMmXvX78YQWTh+fhLewhS8iDKrUDOpw193vouhGhSE
S3zgP51dg2WNgseebKFA/4R0q6R3DYjrlu2q16VgcYlRYoV+DQzGxcSX76Hg
guXT9cs57+uWu6VmVU0oUM58+6Mooxpuj2XfmZpR8OGPGP3I2WrEOJGfxR2h
4HDSRdNcgWo8rFZMKHGgQORTlqz5myqUUFf79zhTwNIZeP/vfhXYgeucZt0p
uH5DJPehWxX636rslT9HwY5BPa+9mlX4sn2jOrwp+OLBknwsVAXe1E2y9pcp
eKKsm2www4T4gi28QVcpsJ0qoI4MMKF0VOfjw1AKVuSuyQqsY0L7FZprIyio
dr+vRs5nYp/cjuKP0RRcXC9RQk9kwrN3T/D6JAqGs3mq/3kwUWFxmFT4lIKi
rK4BL3EmWous+TsLKTjpus9RYp6B9+JHP/8qo4C8tmqiYIwBwXbHl1urOXlm
5v36Ws6A+f5zllXdFFg9CRF1s2fgZN4FjA5QsOzkbIKgMQO+y3xW8b+ngKl0
VuGJBgOpDQEzRlOcPB/brhvlZ6BgbVCf+3cK3p1oyw/6Voma8FDmzV8UJK3a
rUXpq0T3x/DMZ/8oMPnwqqKyuhKThjei23ipWPhI3eBYXiXmH0V7z/BRUeyU
2TCXUAlRvttHJZdSkfUjV+9STCUUnRIMtESpkJXbw1PIsWZ14nprSSpu6I9U
THFsRE0Ru0ygYu5kgJ9SbCWOBKb9TSZRcSpKVsee48DtT1hv11JxoP9A2RuO
q+7ovZpUpaJoa4DrsluVWPKtN/+HBhWEpFzCLo7D05YlCYKKMSsh/xKOoxZ6
emwwo8K/q0o/8XYlWq2XOmy1pGJk8/fpNo7FCjIsDI5QYRRPyhCOq0T88W5Y
OVMhYe7H78dxcvW2lVd9qcht29x4LL4SQ3JdfBFXqFi50cnvDsekC2dmb4dQ
cTHmlkorx/ZNgl9TrlMx+JUxJMjJI52WPvIkioqdB75G63M86qfT8/w2J5+n
Cvq+HCt1djS9vEuFiMi+6QKOXVROM2tTqIhdwtftz8knK1SgpCWditnB+kmh
cE7+g2nZvU+ocCqNWpBwpRLKmltTR3KpaIw5JEa9WIn8MY/w3yVUJO0c1N7q
XImS3e02pEYqtmV9C3czqMSfVPcDa1upeBhUnPxLpxJb/ywx0OikYukR3/wg
jUpUZGqp7B6mYkBkSU8SlfM8CKXwnJmh4vIF6dWNvJXgP675w+cPFe9NBrZa
/qFj94uW8aB5KvatfrB/5CsdDW6L2hMEFCHft8773xAdbY0nM+hyini1HbXr
K+gQoy1MZJMVoSS7uLe0iA5zv6Sb7asUETXNmtqZS0dCh0bQwDpF/Gy8wdv6
kI4e5dfeY2qKCFzq3+oUTwch1Nn92ybO+oiUFz3BdNgO8tjPaivitSA9Y58n
HUM31Y1EDBQhx89zadMBOshjjboyRoqoDCE5ZurScRwn1Kn7FeG0WM9EXpmO
D1/uyGlaKaKAN5CyRIiOqQPzn496KCL2j0BC1rUKtOZEYOU5ReQU/O6ecq5A
oYBMdI2XIpo8PspuMKyAb6X6hvVXFMH/ri6lZEkFBNVcPP/EKuJSfeij6sBX
UBR58+tWqSJckxYXD15+CZumTB72IhpUtG72f75eDqzRNL0sQIO26Uvtebdy
UEKq01SX0bDXdeLOsn3lGNcZ0o+XpMH13i7z9SLl8M4RDzm+hoZHPAteu9wq
g43Agw5JFRryZZTX+5wvg66j6qr6jTRUb7S5fs28DEvkjOrUdGh451iy67FU
Gca9u6RGttMw5Tf6sPh3KRrbHU8mGNDAkyC2qLanFE9Vv5caGdGw8pmefUd5
KW5FBAjOm9CgwD5NH0kqhffYUut8MxoIpk+/VfuXwnpHYpajJQ1e2eoN106W
gjxbaNxgT8MqOx3/FTqc87lO3LhznoYw9r4R4vsS1O8xX6meSIPxvTPPdDeU
4Ff4y3bj+zQEtss+pciUgFpPjXN8SEOFUE0uH28JfI2+S8Xn0bDJRyb7dUsx
Hl+36skroOHss6rM56XFeFNfmVhXSsOTMY8n8anFUDaOkv/D5OzPgplx9Ewx
rCJ+DoqyaLC44f5Q37IYoQ1HUtc20XCjWjJ91fZiFAjX2O9o4+Q1W5kmtLoY
Q8brqbZdnDzU3VI/ixRDOPLWyPl+Gra6Stxv+1UErca/GTfe0nAujZ5cPFQE
p6UOzo8+0PB/32+N2cteTZ/k+lWk2kTXV64nGhOyv/7gese+EyoUnlX/eWbf
Io0Qca4P798G421cB0U95HHU5/pZszDzsiHXggd6d+aZcv3igLeR6AmuFQ7m
m3fd5Hr0EM39+DDXBXeuhNNGub7a3/dobPz/rXeMfes2/f/257ng8Hk+pf/M
junRC1nPdc7rG1KPL3Ltufsnc2LF6v9cEmYQcGHrmv+8pWfZOduAtf9ZWazL
PjF33X/eE/pneug5105/CKG0Uq5ThuyynjG5FsmZ+F7bxfXXnfNB07zruffr
RXu014Lr5nFDrVgbridtXeu77LhW3PH0s4Mb13EiWzR9rnDt/WQvKyOL61uE
09aTT7l+ejPmk1oR1x/Pd654SefaSs/Oqu0N1xcKAyelermOUcrwOzLENXvZ
RNrHca5HryzVUPnCNe+MSq3nDNfyzqaWZX+4/t/vbf/n/wFAfkvY
          "]]}, Annotation[#, "Charting`Private`Tag$5005#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.803121441134557*^9, 3.803121447596089*^9}, 
   3.8031312707191772`*^9, 3.803135773370302*^9, {3.8031358595635757`*^9, 
   3.803135863350645*^9}, 3.8031359531605597`*^9, 3.803136265909319*^9, 
   3.8031367026680746`*^9, 3.8031372539049225`*^9, {3.8032144940417876`*^9, 
   3.803214503540983*^9}},
 CellLabel->"Out[34]=",ExpressionUUID->"3cc2cc9e-eb36-4f36-90d0-6d782f8c2698"]
}, 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.8033057721610527`*^9, 
  3.8033057824915733`*^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.803121441310596*^9, 3.8031312708761587`*^9, 
   3.803135773668197*^9, 3.803135868749423*^9, 3.803135953324995*^9, 
   3.8031362660928574`*^9, 3.803137254090999*^9},
 CellLabel->"Out[48]=",ExpressionUUID->"b6b23034-bf59-496b-adbd-e014a8caf25c"]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"FScale", "=", "0.5"}], ";"}], 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{
    "simple", " ", "rate", " ", "analysis", " ", "suggests", " ", "a", " ", 
     "correction", " ", "factor", " ", "2", 
     RowBox[{"ne", "/", 
      RowBox[{"(", 
       RowBox[{"ne", "+", "ng"}], ")"}]}]}], "=", 
    RowBox[{
     RowBox[{
      RowBox[{"1", "/", "2"}], " ", "due", " ", "to", " ", "multi"}], "-", 
     RowBox[{"level", " ", "structure"}]}]}], "*)"}], " "}]], "Input",
 CellChangeTimes->{{3.803057618969963*^9, 3.803057619498081*^9}, {
   3.8030585409060755`*^9, 3.803058582983102*^9}, {3.8030603241077337`*^9, 
   3.8030604045619926`*^9}, {3.8030636999628677`*^9, 3.803063700069892*^9}, {
   3.8030637484673996`*^9, 3.8030638021806755`*^9}, {3.8030699528252354`*^9, 
   3.803069954377393*^9}, 3.8031286269305553`*^9, {3.80313130803823*^9, 
   3.803131308216285*^9}, {3.8031357540975137`*^9, 3.8031357541570935`*^9}, {
   3.803136307945779*^9, 3.803136308032483*^9}, {3.8031382052096834`*^9, 
   3.8031382052755938`*^9}, {3.803138415622952*^9, 3.803138415696938*^9}, {
   3.8031385548110213`*^9, 3.8031385549059258`*^9}, {3.80313888593921*^9, 
   3.803138886015087*^9}, 
   3.8033057223622704`*^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", "*", 
     "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.803131357627454*^9, 3.8031313744846964`*^9}, {
   3.803131553939128*^9, 3.8031315569707103`*^9}, {3.803131996909254*^9, 
   3.803131999145404*^9}, {3.803136032992083*^9, 3.803136034945009*^9}, {
   3.8031364786933947`*^9, 3.803136480824342*^9}, {3.8031373476633673`*^9, 
   3.8031373480658627`*^9}, {3.8031374520130043`*^9, 
   3.8031374541177683`*^9}, {3.803137539411834*^9, 3.803137541560197*^9}, {
   3.8031377527699604`*^9, 3.8031377555909896`*^9}, 
   3.8033057299747815`*^9},ExpressionUUID->"132a4f70-8a91-4466-aa4f-\
60aab504aaef"],

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[{"#", ",", "2.4", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", 
              RowBox[{"Intensity", "/", "isat"}]}], "]"}], ",", 
            RowBox[{"25", "/", "\[CapitalGamma]"}], ",", "50"}], "]"}], "&"}],
          ",", "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.8031214647259483`*^9, 
   3.80312147123936*^9}, {3.8031323040087533`*^9, 3.8031323067583685`*^9}, {
   3.8031379590169334`*^9, 3.8031379597901497`*^9}, {3.8031380041508675`*^9, 
   3.8031380042030983`*^9}, {3.803305715998209*^9, 
   3.803305716356146*^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[315]:=",ExpressionUUID->"bc6dac95-1b93-4ce8-9af6-755f797db9e8"],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"ainterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"ainterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "2.25"}]}], ";", 
  RowBox[{"(*", 
   RowBox[{
    RowBox[{"compute", " ", "the", " ", "upper", " ", "bound"}], ",", " ", 
    RowBox[{
    "corresponding", " ", "to", " ", "a", " ", "2.25", " ", "larger", " ", 
     "force", " ", 
     RowBox[{"(", 
      RowBox[{
       RowBox[{"~", "180"}], " ", "rather", " ", 
       RowBox[{"than", " ", "~", "80"}], " ", "photons"}], ")"}]}]}], "*)"}], 
  "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{"rinterp", "[", 
    RowBox[{"a_", ",", "b_", ",", "c_", ",", "d_"}], "]"}], ":=", 
   RowBox[{
    RowBox[{"rinterp32", "[", 
     RowBox[{"a", ",", "b", ",", "c", ",", "d"}], "]"}], "*", "FScale", "*", 
    "2.25"}]}], ";"}]], "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.8031332953746634`*^9, 3.8031333008714914`*^9}, {3.803305753996015*^9, 
   3.8033057545802984`*^9}},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[{"#", ",", "2.4", ",", 
            RowBox[{"QuantityMagnitude", "[", 
             RowBox[{"2", 
              RowBox[{"Intensity", "/", "isat"}]}], "]"}], ",", 
            RowBox[{"25", "/", "\[CapitalGamma]"}], ",", "180"}], "]"}], 
          "&"}], ",", "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.803121482514716*^9, 3.8031214859365683`*^9}, {
   3.803133305816197*^9, 3.8031333085654263`*^9}, {3.803133472968872*^9, 
   3.8031334730178823`*^9}, {3.8031339410236087`*^9, 
   3.8031339410981894`*^9}, {3.8031379836665435`*^9, 3.8031379913277683`*^9}},
 CellLabel->
  "In[319]:=",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.8031334753778987`*^9, 3.803133475449915*^9}},
 CellLabel->
  "In[321]:=",ExpressionUUID->"6ce184cd-c0f6-464e-919c-89cc33289447"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{2560, 1357},
WindowMargins->{{-8, Automatic}, {Automatic, -8}},
Magnification:>1.2 Inherited,
FrontEndVersion->"12.0 for Microsoft Windows (64-bit) (April 8, 2019)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[558, 20, 238, 4, 35, "Input",ExpressionUUID->"acf3f024-06f7-4b67-aa47-b28ae5a429d7"],
Cell[CellGroupData[{
Cell[821, 28, 493, 8, 53, "Subsubsection",ExpressionUUID->"5aaf2c87-6956-45b5-ba7b-95459d133cb4"],
Cell[1317, 38, 20847, 396, 547, "Input",ExpressionUUID->"7824fddb-3943-46b5-acd8-799c34acddf8"]
}, Open  ]],
Cell[CellGroupData[{
Cell[22201, 439, 415, 6, 65, "Subsection",ExpressionUUID->"c541093d-3069-439e-a54a-d38289364e3b"],
Cell[22619, 447, 6345, 184, 372, "Input",ExpressionUUID->"a22cab01-84b6-4910-a8ee-e13c7c9bda3f"],
Cell[CellGroupData[{
Cell[28989, 635, 496, 7, 53, "Subsubsection",ExpressionUUID->"2086d6f6-0d0b-40d9-a62f-2851fbd41666"],
Cell[29488, 644, 4063, 95, 420, "Input",ExpressionUUID->"faefb1d2-c143-4cf4-bbad-8f8982efcabd"]
}, Open  ]],
Cell[CellGroupData[{
Cell[33588, 744, 524, 7, 53, "Subsubsection",ExpressionUUID->"20effcbb-203d-40b4-8029-107270944d9f"],
Cell[34115, 753, 9793, 193, 557, "Input",ExpressionUUID->"dccf9f9d-0271-4cee-a2a6-51b034d09fcc"]
}, Open  ]],
Cell[CellGroupData[{
Cell[43945, 951, 578, 8, 53, "Subsubsection",ExpressionUUID->"6c519e3a-f10c-4538-afc2-24105e47eb71"],
Cell[44526, 961, 2859, 41, 35, "Input",ExpressionUUID->"9e84a53b-6a94-415f-8b47-20c1efb51727"],
Cell[47388, 1004, 8852, 163, 372, "Input",ExpressionUUID->"c0c8e4fa-0227-4b9a-82f2-c524254ed188"],
Cell[56243, 1169, 2858, 43, 35, "Input",ExpressionUUID->"c5609752-0026-4bbc-8872-3b54364462a9"],
Cell[CellGroupData[{
Cell[59126, 1216, 16054, 235, 84, "Input",ExpressionUUID->"e65b6191-864b-4f9a-bf41-1a27cb9ce5bb"],
Cell[75183, 1453, 50591, 892, 450, "Output",ExpressionUUID->"3cc2cc9e-eb36-4f36-90d0-6d782f8c2698"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[125823, 2351, 1049, 19, 56, "Subsubsection",ExpressionUUID->"c44a9341-abbe-437f-a7b8-e777560bd4ec"],
Cell[CellGroupData[{
Cell[126897, 2374, 1014, 25, 35, "Input",ExpressionUUID->"765a6f84-5baa-41bd-a875-c58ef6d443ae"],
Cell[127914, 2401, 1000, 16, 39, "Output",ExpressionUUID->"b6b23034-bf59-496b-adbd-e014a8caf25c"]
}, Open  ]],
Cell[128929, 2420, 1298, 27, 35, "Input",ExpressionUUID->"5ceef669-cf9d-4c9f-996a-85e2b663751e"],
Cell[130230, 2449, 1548, 37, 60, "Input",ExpressionUUID->"132a4f70-8a91-4466-aa4f-60aab504aaef"],
Cell[131781, 2488, 9209, 173, 252, "Input",ExpressionUUID->"ccc1bdba-34c7-4530-8f76-36d6baa364a0"],
Cell[140993, 2663, 3251, 84, 324, "Input",ExpressionUUID->"bc6dac95-1b93-4ce8-9af6-755f797db9e8"],
Cell[144247, 2749, 1392, 33, 60, "Input",ExpressionUUID->"e1d13343-c95f-46a7-bd96-775b2205822d"],
Cell[145642, 2784, 3759, 94, 228, "Input",ExpressionUUID->"ab82b3ea-5c09-40e3-a005-6f650494affb"],
Cell[149404, 2880, 3354, 86, 324, "Input",ExpressionUUID->"6ce184cd-c0f6-464e-919c-89cc33289447"]
}, Open  ]]
}, Open  ]]
}
]
*)