visualization.nb

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

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

(* CreatedBy='Mathematica 11.0' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     24964,        789]
NotebookOptionsPosition[     22370,        706]
NotebookOutlinePosition[     22783,        724]
CellTagsIndexPosition[     22740,        721]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Sweep Visualizations", "Section"],

Cell["Set the working directory.", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"SetDirectory", "[", 
   RowBox[{"NotebookDirectory", "[", "]"}], "]"}], ";"}]], "Input"],

Cell["Import data.", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"Import", "[", "\"\<plots-math.csv\>\"", "]"}]}], ";"}]], "Input"],

Cell["\<\
Data comes in 17 columns:
1: phi0
2: phi0dot
3: H
4: rho
5:  deltarho2
6: phi2pt
7: psirms
8: efolds
9: kappa
10: type (0 -> Hartree Off, 1 -> Bunch Davies, 2-> Perturbed)
11: filename
12: runtime (in seconds)
13: \[Delta]\[Rho]2/\[Rho]
14: inflationstart
15: slowrollstart
16: kappacrossing
17: lastcrossing\
\>", "Text"],

Cell["Which index do you want to work with?", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"index", "=", "8"}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{
  RowBox[{"names", "=", 
   RowBox[{"{", 
    RowBox[{
    "\"\<\[Phi]0\>\"", ",", 
     "\"\<\!\(\*OverscriptBox[\(\[Phi]\), \(.\)]\)0\>\"", ",", "\"\<H\>\"", 
     ",", "\"\<\[Rho]\>\"", ",", 
     "\"\<\!\(\*SubscriptBox[\(\[Delta]\[Rho]\), \(2\)]\)\>\"", ",", 
     "\"\<<\!\(\*SuperscriptBox[\(\[Delta]\[Phi]\), \(2\)]\)>\>\"", ",", 
     "\"\<\!\(\*SubscriptBox[\(\[Psi]\), \(RMS\)]\)\>\"", ",", 
     "\"\<\!\(\*SubscriptBox[\(N\), \(e - folds\)]\)\>\"", ",", 
     "\"\<\[Kappa]\>\"", ",", "\"\<Type\>\"", ",", "\"\<Filename\>\"", ",", 
     "\"\<runtime\>\"", ",", "\"\<\[Delta]\[Rho]2/\[Rho]\>\"", ",", 
     "\"\<Inflation Start\>\"", ",", " ", "\"\<Slow roll start\>\"", ",", 
     " ", "\"\<Kappa crossing\>\"", ",", " ", "\"\<Last crossing\>\""}], 
    "}"}]}], ";"}]], "Input"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"criticalEfolds", "=", "65"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"critplane", "=", 
   RowBox[{"Plot3D", "[", "\[IndentingNewLine]", 
    RowBox[{"criticalEfolds", ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "50"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"y", ",", 
       RowBox[{"-", "1"}], ",", "1"}], "}"}], ",", "\[IndentingNewLine]", 
     RowBox[{"ColorFunction", "\[Rule]", 
      RowBox[{"Function", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
        RowBox[{"RGBColor", "[", 
         RowBox[{"0.5", ",", "0.5", ",", "0.5", ",", "0.5"}], "]"}]}], 
       "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ";"}]}], "Input"],

Cell["Select all \[Phi]0 and \[Phi]0dot values.", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"\[Phi]0vals", "=", 
   RowBox[{"Sort", "[", 
    RowBox[{"DeleteDuplicates", "[", 
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", "1"}], "]"}], "]"}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"\[Phi]0dotvals", "=", 
   RowBox[{"Sort", "[", 
    RowBox[{"DeleteDuplicates", "[", 
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{";;", ",", "2"}], "]"}], "]"}], "]"}], "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Length", "/@", 
  RowBox[{"{", 
   RowBox[{"\[Phi]0vals", ",", "\[Phi]0dotvals"}], "}"}]}]}], "Input"],

Cell[CellGroupData[{

Cell["Functions", "Subsection"],

Cell["Select specific kinds of data", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"(*", " ", 
   RowBox[{"Selects", " ", "Hartree", " ", "off", " ", "data", " ", "only"}], 
   " ", "*)"}], "\[IndentingNewLine]", 
  RowBox[{
   RowBox[{
    RowBox[{"SelectOff", "[", "data_", "]"}], ":=", 
    RowBox[{"Select", "[", 
     RowBox[{"data", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"#", "[", 
         RowBox[{"[", "10", "]"}], "]"}], "\[Equal]", "0"}], "&"}]}], "]"}]}],
    "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{
    "Selects", " ", "Bunch", " ", "Davies", " ", "data", " ", "only"}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SelectBD", "[", "data_", "]"}], ":=", 
    RowBox[{"Select", "[", 
     RowBox[{"data", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"#", "[", 
         RowBox[{"[", "10", "]"}], "]"}], "\[Equal]", "1"}], "&"}]}], "]"}]}],
    "\[IndentingNewLine]", 
   RowBox[{"(*", " ", 
    RowBox[{"Selects", " ", "perturbed", " ", "data", " ", "only"}], " ", 
    "*)"}], "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"SelectPerturbed", "[", "data_", "]"}], ":=", 
    RowBox[{"Select", "[", 
     RowBox[{"data", ",", 
      RowBox[{
       RowBox[{
        RowBox[{"#", "[", 
         RowBox[{"[", "10", "]"}], "]"}], "==", "2"}], "&"}]}], 
     "]"}]}]}]}]], "Input"],

Cell["\<\
Extract specific types of data (specify the column number that you want)\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"SelectQuantity", "[", 
   RowBox[{"data_", ",", "index_"}], "]"}], ":=", 
  RowBox[{"data", "[", 
   RowBox[{"[", 
    RowBox[{";;", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "2", ",", "index"}], "}"}]}], "]"}], "]"}]}]], "Input"],

Cell["Helper function to find index for a value in a list.", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"FindIndex", "[", 
   RowBox[{"value_", ",", "list_"}], "]"}], ":=", 
  RowBox[{
   RowBox[{"Position", "[", 
    RowBox[{"list", ",", "value"}], "]"}], "[", 
   RowBox[{"[", 
    RowBox[{"1", ",", "1"}], "]"}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Hartree Off", "Subsection"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"workingoff", "=", 
   RowBox[{"SelectOff", "[", "data", "]"}]}], ";"}], "\n", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Plotting \>\"", " ", "<>", " ", 
   RowBox[{"names", "[", 
    RowBox[{"[", "index", "]"}], "]"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"offplot", "=", 
   RowBox[{"ListPlot3D", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"SelectQuantity", "[", 
      RowBox[{"workingoff", ",", "index"}], "]"}], ",", "\[IndentingNewLine]", 
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"names", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"names", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"names", "[", 
         RowBox[{"[", "index", "]"}], "]"}]}], "}"}]}], ",", 
     RowBox[{"ColorFunction", "\[Rule]", 
      RowBox[{"Function", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
        RowBox[{
         RowBox[{"ColorData", "[", "\"\<NeonColors\>\"", "]"}], "[", "z", 
         "]"}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"offplot", ",", "critplane"}], "]"}]}], "Input"],

Cell["Same thing, but in ContourPlot form.", "Text"],

Cell[BoxData[
 RowBox[{"ListContourPlot", "[", 
  RowBox[{
   RowBox[{"SelectQuantity", "[", 
    RowBox[{"workingoff", ",", "index"}], "]"}], ",", 
   RowBox[{"Contours", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"45", ",", "50", ",", "55", ",", "60", ",", "65", ",", "70"}], 
     "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"names", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"names", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}]}], "]"}]], "Input"],

Cell["This is the critical line.", "Text"],

Cell[BoxData[
 RowBox[{"Offline", "=", 
  RowBox[{"ListContourPlot", "[", 
   RowBox[{
    RowBox[{"SelectQuantity", "[", 
     RowBox[{"workingoff", ",", "index"}], "]"}], ",", 
    RowBox[{"Contours", "\[Rule]", 
     RowBox[{"{", "criticalEfolds", "}"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"names", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"names", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ",", 
    RowBox[{"ColorFunction", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "1", ",", "1", ",", "0"}], "]"}]}], "]"}]}], ",", 
    RowBox[{"ContourStyle", "\[Rule]", "Blue"}]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Bunch-Davies", "Subsection"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"workingbd", "=", 
   RowBox[{"SelectBD", "[", "data", "]"}]}], ";"}], "\n", 
 RowBox[{"Print", "[", 
  RowBox[{"\"\<Plotting \>\"", " ", "<>", " ", 
   RowBox[{"names", "[", 
    RowBox[{"[", "index", "]"}], "]"}]}], "]"}], "\[IndentingNewLine]", 
 RowBox[{
  RowBox[{"bdplot", "=", 
   RowBox[{"ListPlot3D", "[", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"SelectQuantity", "[", 
      RowBox[{"workingbd", ",", "index"}], "]"}], ",", "\[IndentingNewLine]", 
     
     RowBox[{"AxesLabel", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"names", "[", 
         RowBox[{"[", "1", "]"}], "]"}], ",", 
        RowBox[{"names", "[", 
         RowBox[{"[", "2", "]"}], "]"}], ",", 
        RowBox[{"names", "[", 
         RowBox[{"[", "index", "]"}], "]"}]}], "}"}]}], ",", 
     RowBox[{"ColorFunction", "\[Rule]", 
      RowBox[{"Function", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
        RowBox[{
         RowBox[{"ColorData", "[", "\"\<NeonColors\>\"", "]"}], "[", "z", 
         "]"}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], 
  ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{"bdplot", ",", "critplane"}], "]"}]}], "Input"],

Cell["Same thing, but in ContourPlot form.", "Text"],

Cell[BoxData[
 RowBox[{"ListContourPlot", "[", 
  RowBox[{
   RowBox[{"SelectQuantity", "[", 
    RowBox[{"workingbd", ",", "index"}], "]"}], ",", 
   RowBox[{"Contours", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"45", ",", "50", ",", "55", ",", "60", ",", "65", ",", "70"}], 
     "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"names", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"names", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}]}], "]"}]], "Input"],

Cell["This is the critical line.", "Text"],

Cell[BoxData[
 RowBox[{"BDline", "=", 
  RowBox[{"ListContourPlot", "[", 
   RowBox[{
    RowBox[{"SelectQuantity", "[", 
     RowBox[{"workingbd", ",", "index"}], "]"}], ",", 
    RowBox[{"Contours", "\[Rule]", 
     RowBox[{"{", "criticalEfolds", "}"}]}], ",", 
    RowBox[{"FrameLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"names", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"names", "[", 
        RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ",", 
    RowBox[{"ColorFunction", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
       RowBox[{"RGBColor", "[", 
        RowBox[{"1", ",", "1", ",", "1", ",", "0"}], "]"}]}], "]"}]}], ",", 
    RowBox[{"ContourStyle", "\[Rule]", "Red"}]}], "]"}]}]], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Perturbed", "Subsection"],

Cell[CellGroupData[{

Cell["Obtaining statistics", "Subsubsection"],

Cell[BoxData[
 RowBox[{
  RowBox[{"workingperturbed", "=", 
   RowBox[{"SelectPerturbed", "[", "data", "]"}]}], ";"}]], "Input"],

Cell["\<\
Workhorse function that takes in the desired data, and outputs a list of {\
\[Phi]0, \[Phi]0dot, mean, \[Sigma]} values.\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"ExtractStats", "[", 
   RowBox[{"data_", ",", "index_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"grid", ",", "i", ",", "efolds", ",", "mean", ",", "sd"}], "}"}],
     ",", "\[IndentingNewLine]", 
    RowBox[{"(*", " ", 
     RowBox[{
     "Start", " ", "by", " ", "sorting", " ", "all", " ", "the", " ", "data", 
      " ", "by", " ", "\[Phi]0", " ", "and", " ", "\[Phi]0dot", " ", 
      "value"}], " ", "*)"}], "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"grid", "=", 
      RowBox[{"Table", "[", 
       RowBox[{
        RowBox[{"{", 
         RowBox[{"i", ",", "j", ",", 
          RowBox[{"{", "}"}]}], "}"}], " ", ",", 
        RowBox[{"{", 
         RowBox[{"i", ",", "\[Phi]0vals"}], "}"}], ",", 
        RowBox[{"{", 
         RowBox[{"j", ",", "\[Phi]0dotvals"}], "}"}]}], "]"}]}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"Do", "[", "\[IndentingNewLine]", 
      RowBox[{
       RowBox[{"AppendTo", "[", "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"grid", "[", 
          RowBox[{"[", 
           RowBox[{
            RowBox[{"FindIndex", "[", 
             RowBox[{
              RowBox[{"i", "[", 
               RowBox[{"[", "1", "]"}], "]"}], ",", "\[Phi]0vals"}], "]"}], 
            ",", 
            RowBox[{"FindIndex", "[", 
             RowBox[{
              RowBox[{"i", "[", 
               RowBox[{"[", "2", "]"}], "]"}], ",", "\[Phi]0dotvals"}], "]"}],
             ",", "3"}], "]"}], "]"}], ",", "\[IndentingNewLine]", "i"}], 
        "]"}], "\[IndentingNewLine]", ",", 
       RowBox[{"{", 
        RowBox[{"i", ",", "data"}], "}"}]}], "]"}], ";", 
     "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"Flatten", " ", "into", " ", "a", " ", "list"}], " ", "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"grid", "=", 
      RowBox[{"Flatten", "[", 
       RowBox[{"grid", ",", "1"}], "]"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"Compute", " ", "statistics"}], " ", "*)"}], 
     "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "\[LessEqual]", 
        RowBox[{"Length", "[", "grid", "]"}]}], ",", 
       RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{
        RowBox[{"efolds", "=", 
         RowBox[{
          RowBox[{"grid", "[", 
           RowBox[{"[", "i", "]"}], "]"}], "[", 
          RowBox[{"[", 
           RowBox[{"3", ",", ";;", ",", "index"}], "]"}], "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{"mean", "=", 
         RowBox[{"Mean", "[", "efolds", "]"}]}], ";", "\[IndentingNewLine]", 
        RowBox[{"sd", "=", 
         RowBox[{"StandardDeviation", "[", "efolds", "]"}]}], ";", 
        "\[IndentingNewLine]", 
        RowBox[{
         RowBox[{"grid", "[", 
          RowBox[{"[", "i", "]"}], "]"}], "=", 
         RowBox[{
          RowBox[{"grid", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", 
             RowBox[{"1", ";;", "2"}]}], "]"}], "]"}], "~", "Join", "~", 
          RowBox[{"{", 
           RowBox[{"mean", ",", "sd"}], "}"}]}]}]}]}], "\[IndentingNewLine]", 
      "]"}], ";", "\[IndentingNewLine]", 
     RowBox[{"(*", " ", 
      RowBox[{"Report", " ", "results"}], " ", "*)"}], "\[IndentingNewLine]", 
     "grid"}]}], "\[IndentingNewLine]", "]"}]}]], "Input"],

Cell["Obtain all the statistics.", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"stats", "=", 
   RowBox[{"ExtractStats", "[", 
    RowBox[{"workingperturbed", ",", "index"}], "]"}]}], ";"}]], "Input"],

Cell["Find any points with suspiciously large standard deviations.", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"FindPoints", "[", 
   RowBox[{"data_", ",", "entry_"}], "]"}], ":=", 
  RowBox[{"data", "[", 
   RowBox[{"[", 
    RowBox[{"Flatten", "[", 
     RowBox[{"Position", "[", 
      RowBox[{
       RowBox[{"data", "[", 
        RowBox[{"[", 
         RowBox[{";;", ",", 
          RowBox[{"{", 
           RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "]"}], ",", 
       RowBox[{"entry", "[", 
        RowBox[{"[", 
         RowBox[{"1", ";;", "2"}], "]"}], "]"}]}], "]"}], "]"}], "]"}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"badpoints", "=", 
  RowBox[{"Select", "[", 
   RowBox[{"stats", ",", 
    RowBox[{
     RowBox[{
      RowBox[{"#", "[", 
       RowBox[{"[", "4", "]"}], "]"}], ">", "1.0"}], "&"}]}], 
   "]"}]}], "\[IndentingNewLine]", 
 RowBox[{"Map", "[", 
  RowBox[{
   RowBox[{
    RowBox[{"FindPoints", "[", 
     RowBox[{"data", ",", "#"}], "]"}], "&"}], ",", "badpoints"}], 
  "]"}]}], "Input"]
}, Open  ]],

Cell[CellGroupData[{

Cell["Plotting", "Subsubsection"],

Cell["\<\
stats now contains a list of {phi0, phi0dot, mean, \[Sigma]} values.\
\>", "Text"],

Cell["\<\
Here is a function to add a specified number of standard deviations to the \
means, returning data ready for plotting.\
\>", "Text"],

Cell[BoxData[
 RowBox[{
  RowBox[{"Add\[Sigma]", "[", 
   RowBox[{"data_", ",", "\[Sigma]s_"}], "]"}], ":=", 
  RowBox[{"Module", "[", 
   RowBox[{
    RowBox[{"{", "result", "}"}], ",", "\[IndentingNewLine]", 
    RowBox[{
     RowBox[{"result", "=", 
      RowBox[{"{", "}"}]}], ";", "\[IndentingNewLine]", 
     RowBox[{"For", "[", 
      RowBox[{
       RowBox[{"i", "=", "1"}], ",", 
       RowBox[{"i", "\[LessEqual]", 
        RowBox[{"Length", "[", "data", "]"}]}], ",", 
       RowBox[{"i", "++"}], ",", "\[IndentingNewLine]", 
       RowBox[{"AppendTo", "[", 
        RowBox[{"result", ",", 
         RowBox[{
          RowBox[{"data", "[", 
           RowBox[{"[", 
            RowBox[{"i", ",", 
             RowBox[{"{", 
              RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "]"}], "~", "Join", 
          "~", 
          RowBox[{"{", 
           RowBox[{
            RowBox[{"data", "[", 
             RowBox[{"[", 
              RowBox[{"i", ",", "3"}], "]"}], "]"}], "+", 
            RowBox[{"\[Sigma]s", " ", 
             RowBox[{"data", "[", 
              RowBox[{"[", 
               RowBox[{"i", ",", "4"}], "]"}], "]"}]}]}], "}"}]}]}], "]"}]}], 
      "\[IndentingNewLine]", "]"}], ";", "\[IndentingNewLine]", "result"}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input"],

Cell["Here are the means.", "Text"],

Cell[BoxData[
 RowBox[{"means", "=", 
  RowBox[{"ListPlot3D", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"Add\[Sigma]", "[", 
     RowBox[{"stats", ",", "0"}], "]"}], ",", "\[IndentingNewLine]", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"names", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"names", "[", 
        RowBox[{"[", "2", "]"}], "]"}], ",", 
       RowBox[{"names", "[", 
        RowBox[{"[", "index", "]"}], "]"}]}], "}"}]}], ",", 
    RowBox[{"ColorFunction", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
       RowBox[{"Hue", "[", 
        RowBox[{".65", 
         RowBox[{"(", 
          RowBox[{"1", "-", "z"}], ")"}]}], "]"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input"],

Cell["Here are the standard deviations.", "Text"],

Cell[BoxData[
 RowBox[{"\[Sigma]s", "=", 
  RowBox[{"ListPlot3D", "[", "\[IndentingNewLine]", 
   RowBox[{
    RowBox[{"stats", "[", 
     RowBox[{"[", 
      RowBox[{";;", ",", 
       RowBox[{"{", 
        RowBox[{"1", ",", "2", ",", "4"}], "}"}]}], "]"}], "]"}], ",", 
    "\[IndentingNewLine]", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"names", "[", 
        RowBox[{"[", "1", "]"}], "]"}], ",", 
       RowBox[{"names", "[", 
        RowBox[{"[", "2", "]"}], "]"}], ",", "\"\<\[Sigma]\>\""}], "}"}]}], 
    ",", 
    RowBox[{"ColorFunction", "\[Rule]", 
     RowBox[{"Function", "[", 
      RowBox[{
       RowBox[{"{", 
        RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
       RowBox[{"Hue", "[", 
        RowBox[{".65", 
         RowBox[{"(", 
          RowBox[{"1", "-", "z"}], ")"}]}], "]"}]}], "]"}]}]}], 
   "\[IndentingNewLine]", "]"}]}]], "Input"],

Cell["\<\
We can compare the means to Bunch-Davies and the critical plane.\
\>", "Text"],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"means", ",", "bdplot", ",", "critplane"}], "]"}]], "Input"],

Cell["Here\[CloseCurlyQuote]s a contour plot of the means.", "Text"],

Cell[BoxData[
 RowBox[{"ListContourPlot", "[", 
  RowBox[{
   RowBox[{"Add\[Sigma]", "[", 
    RowBox[{"stats", ",", "0"}], "]"}], ",", 
   RowBox[{"Contours", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"45", ",", "50", ",", "55", ",", "60", ",", "65", ",", "70"}], 
     "}"}]}], ",", 
   RowBox[{"PlotLegends", "\[Rule]", "Automatic"}], ",", 
   RowBox[{"FrameLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"names", "[", 
       RowBox[{"[", "1", "]"}], "]"}], ",", 
      RowBox[{"names", "[", 
       RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}]}], "]"}]], "Input"],

Cell["\<\
Here\[CloseCurlyQuote]s the critical line for the means, means +- 2\[Sigma], \
the Bunch-Davies result in red, and the Hartree off result in blue.\
\>", "Text"],

Cell[BoxData[{
 RowBox[{
  RowBox[{"pm", "=", "2"}], ";"}], "\[IndentingNewLine]", 
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"Table", "[", 
    RowBox[{
     RowBox[{"ListContourPlot", "[", 
      RowBox[{
       RowBox[{"Add\[Sigma]", "[", 
        RowBox[{"stats", ",", "i"}], "]"}], ",", 
       RowBox[{"Contours", "\[Rule]", 
        RowBox[{"{", "criticalEfolds", "}"}]}], ",", 
       RowBox[{"FrameLabel", "\[Rule]", 
        RowBox[{"{", 
         RowBox[{
          RowBox[{"names", "[", 
           RowBox[{"[", "1", "]"}], "]"}], ",", 
          RowBox[{"names", "[", 
           RowBox[{"[", "2", "]"}], "]"}]}], "}"}]}], ",", 
       RowBox[{"ColorFunction", "\[Rule]", 
        RowBox[{"Function", "[", 
         RowBox[{
          RowBox[{"{", 
           RowBox[{"x", ",", "y", ",", "z"}], "}"}], ",", 
          RowBox[{"RGBColor", "[", 
           RowBox[{"1", ",", "1", ",", "1", ",", "0"}], "]"}]}], "]"}]}]}], 
      "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"i", ",", 
       RowBox[{"{", 
        RowBox[{"0", ",", "pm", ",", 
         RowBox[{"-", "pm"}]}], "}"}]}], "}"}]}], "]"}], ",", "BDline", ",", 
   "Offline"}], "]"}]}], "Input"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1440, 1030},
WindowMargins->{{0, Automatic}, {Automatic, 0}},
TrackCellChangeTimes->False,
Magnification:>1.5 Inherited,
FrontEndVersion->"11.0 for Mac OS X x86 (32-bit, 64-bit Kernel) (July 28, \
2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 39, 0, 95, "Section"],
Cell[622, 24, 42, 0, 46, "Text"],
Cell[667, 26, 123, 3, 48, "Input"],
Cell[793, 31, 28, 0, 46, "Text"],
Cell[824, 33, 125, 3, 48, "Input"],
Cell[952, 38, 332, 19, 548, "Text"],
Cell[1287, 59, 53, 0, 46, "Text"],
Cell[1343, 61, 71, 2, 48, "Input"],
Cell[1417, 65, 816, 16, 86, "Input"],
Cell[2236, 83, 757, 19, 173, "Input"],
Cell[2996, 104, 57, 0, 46, "Text"],
Cell[3056, 106, 626, 19, 111, "Input"],
Cell[CellGroupData[{
Cell[3707, 129, 31, 0, 65, "Subsection"],
Cell[3741, 131, 45, 0, 46, "Text"],
Cell[3789, 133, 1305, 39, 204, "Input"],
Cell[5097, 174, 96, 2, 46, "Text"],
Cell[5196, 178, 272, 8, 48, "Input"],
Cell[5471, 188, 68, 0, 46, "Text"],
Cell[5542, 190, 264, 8, 48, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[5843, 203, 33, 0, 65, "Subsection"],
Cell[5879, 205, 1267, 33, 266, "Input"],
Cell[7149, 240, 52, 0, 46, "Text"],
Cell[7204, 242, 597, 16, 80, "Input"],
Cell[7804, 260, 42, 0, 46, "Text"],
Cell[7849, 262, 844, 22, 111, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8730, 289, 34, 0, 65, "Subsection"],
Cell[8767, 291, 1268, 34, 266, "Input"],
Cell[10038, 327, 52, 0, 46, "Text"],
Cell[10093, 329, 596, 16, 80, "Input"],
Cell[10692, 347, 42, 0, 46, "Text"],
Cell[10737, 349, 841, 22, 111, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[11615, 376, 31, 0, 65, "Subsection"],
Cell[CellGroupData[{
Cell[11671, 380, 45, 0, 52, "Subsubsection"],
Cell[11719, 382, 128, 3, 48, "Input"],
Cell[11850, 387, 144, 3, 46, "Text"],
Cell[11997, 392, 3411, 88, 638, "Input"],
Cell[15411, 482, 42, 0, 46, "Text"],
Cell[15456, 484, 155, 4, 48, "Input"],
Cell[15614, 490, 76, 0, 46, "Text"],
Cell[15693, 492, 947, 31, 111, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16677, 528, 33, 0, 52, "Subsubsection"],
Cell[16713, 530, 92, 2, 46, "Text"],
Cell[16808, 534, 142, 3, 46, "Text"],
Cell[16953, 539, 1295, 35, 235, "Input"],
Cell[18251, 576, 35, 0, 46, "Text"],
Cell[18289, 578, 868, 24, 142, "Input"],
Cell[19160, 604, 49, 0, 46, "Text"],
Cell[19212, 606, 912, 27, 142, "Input"],
Cell[20127, 635, 88, 2, 46, "Text"],
Cell[20218, 639, 107, 2, 48, "Input"],
Cell[20328, 643, 68, 0, 46, "Text"],
Cell[20399, 645, 585, 16, 80, "Input"],
Cell[20987, 663, 170, 3, 46, "Text"],
Cell[21160, 668, 1170, 33, 111, "Input"]
}, Open  ]]
}, Open  ]]
}, Open  ]]
}
]
*)