brief_reportBv3.tex

\documentclass{article}[11pt]
\usepackage{longtable}
\usepackage{graphicx}
\usepackage{fix-cm}
\usepackage[margin=1in,paperwidth=8.5in,paperheight=11in]{geometry}
\usepackage[section]{placeins}
\usepackage{flafter}
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{etoolbox}
\usepackage{units}
\usepackage{multirow}
\usepackage{hyperref}
\usepackage{pdfcomment}
\usepackage{color}

\setcounter{topnumber}{3}
\setcounter{bottomnumber}{3}
\setcounter{totalnumber}{4}
\renewcommand{\topfraction}{0.9}
\renewcommand{\bottomfraction}{0.9}
\renewcommand{\textfraction}{0.1}
\renewcommand{\floatpagefraction}{0.7}

\newcommand{\rrangle}{\rangle\!\rangle} \newcommand{\llangle}{\langle\!\langle}
\newcommand{\ket}[1]{\ensuremath{\left|#1\right\rangle}}
\newcommand{\bra}[1]{\ensuremath{\left\langle#1\right|}}
\newcommand{\braket}[2]{\ensuremath{\left\langle#1|#2\right\rangle}}
\newcommand{\expec}[1]{\ensuremath{\left\langle#1\right\rangle}}
\newcommand{\ketbra}[2]{\ket{#1}\!\!\bra{#2}}
\newcommand{\braopket}[3]{\ensuremath{\bra{#1}#2\ket{#3}}}
\newcommand{\proj}[1]{\ketbra{#1}{#1}}
\newcommand{\sket}[1]{\ensuremath{\left|#1\right\rrangle}}
\newcommand{\sbra}[1]{\ensuremath{\left\llangle#1\right|}}
\newcommand{\sbraket}[2]{\ensuremath{\left\llangle#1|#2\right\rrangle}}
\newcommand{\sketbra}[2]{\sket{#1}\!\!\sbra{#2}}
\newcommand{\sbraopket}[3]{\ensuremath{\sbra{#1}#2\sket{#3}}}
\newcommand{\sproj}[1]{\sketbra{#1}{#1}}
\newcommand{\norm}[1]{\left\lVert#1\right\rVert}
\def\Id{1\!\mathrm{l}}
\newcommand{\Tr}[0]{\mathrm{Tr}}
\providecommand{\e}[1]{\ensuremath{\times 10^{#1}}}

%Command used for python automatic substitution
\newcommand{\putfield}[2]{#2}

\newtoggle{confidences}
\newtoggle{goodnessSection}
\toggletrue{confidences}
\toggletrue{goodnessSection}


\hypersetup{
  pdfinfo={ Author={pyGSTi},
Title={Brief GST report for MyDefaultReportName},
Keywords={GST},
pyGSTi_Version={0.9.3},
opt_long_tables={False},
opt_table_class={pygstiTbl},
opt_template_path={None},
opt_latex_cmd={pdflatex},
L_germ_tuple_base_string_dict={Dict[[0, Gx]: {}, [0, Gy]: {}, [0, Gi]: {}, [0, GxGy]: {}, [0, GxGyGi]: {}, [0, GxGiGy]: {}, [0, GxGiGi]: {}, [0, GyGiGi]: {}, [0, GxGxGiGy]: {}, [0, GxGyGyGi]: {}, [0, GxGxGyGxGyGy]: {}, [1, Gx]: [Gx], [1, Gy]: [Gy], [1, Gi]: [Gi], [1, GxGy]: {}, [1, GxGyGi]: {}, [1, GxGiGy]: {}, [1, GxGiGi]: {}, [1, GyGiGi]: {}, [1, GxGxGiGy]: {}, [1, GxGyGyGi]: {}, [1, GxGxGyGxGyGy]: {}, [2, Gx]: [Gx]2, [2, Gy]: [Gy]2, [2, Gi]: [Gi]2, [2, GxGy]: [GxGy], [2, GxGyGi]: {}, [2, GxGiGy]: {}, [2, GxGiGi]: {}, [2, GyGiGi]: {}, [2, GxGxGiGy]: {}, [2, GxGyGyGi]: {}, [2, GxGxGyGxGyGy]: {}, [4, Gx]: [Gx]4, [4, Gy]: [Gy]4, [4, Gi]: [Gi]4, [4, GxGy]: [GxGy]2, [4, GxGyGi]: [GxGyGi], [4, GxGiGy]: [GxGiGy], [4, GxGiGi]: [GxGiGi], [4, GyGiGi]: [GyGiGi], [4, GxGxGiGy]: [GxGxGiGy], [4, GxGyGyGi]: [GxGyGyGi], [4, GxGxGyGxGyGy]: {}, [8, Gx]: [Gx]8, [8, Gy]: [Gy]8, [8, Gi]: [Gi]8, [8, GxGy]: [GxGy]4, [8, GxGyGi]: [GxGyGi]2, [8, GxGiGy]: [GxGiGy]2, [8, GxGiGi]: [GxGiGi]2, [8, GyGiGi]: [GyGiGi]2, [8, GxGxGiGy]: [GxGxGiGy]2, [8, GxGyGyGi]: [GxGyGyGi]2, [8, GxGxGyGxGyGy]: [GxGxGyGxGyGy]]},
cptpPenaltyFactor={0},
defaultBasename={MyDefaultReportName},
defaultDirectory={temp_test_files},
distributeMethod={deriv},
fiducial_pairs={None},
gaugeOptParams={Dict[TPpenalty: 0, CPpenalty: 0, gatesMetric: frobenius, spamMetric: frobenius, itemWeights: Dict[spam: 0.001, gates: 1.0]]},
hessianProjection={optimal gate CIs},
linlogPercentile={5},
max_length_list={[0, 1, 2, 4, 8]},
memLimit={None},
minProbClip={1e-06},
minProbClipForWeighting={0.0001},
objective={chi2},
probClipInterval={[-1000000.0, 1000000.0]},
radius={0.0001},
weights={None}  }
}


\begin{document}

\title{Brief GST report for MyDefaultReportName}
\date{\vspace{-1cm}\today}
%\author{}

\begingroup
\let\center\flushleft
\let\endcenter\endflushleft
\maketitle
\endgroup

\section{Overview}
This report presents a brief gate-set tomography (GST) analysis of a dataset called ``MyDefaultReportName''.   It is intended for folks who are already familiar with GST output tables and plots.  If this isn't you, you should have GST generate a ``full report'' which give much more detail about what GST is and what the results in this ``brief report'' mean.

\section{Goodness-of-model}
\iftoggle{goodnessSection}{
Table \ref{progressTable} shows how well GST was able to fit you data through its iterations.
%, and Figure \ref{bestEstimateColorBoxPlot} shows in more detail the goodness of fit quantity for the final (best) gateset estimate.

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|c|c|c|}
\hline
L & $\chi^2$ & $k$ & $\chi^2-k$ & $\sqrt{2k}$ & $N_\sigma$ & $N_s$ & $N_p$ & Rating \\ \hline
0 & 53.25 & 52 & 1.245 & 10.2 & 0.12 & 92 & 40 & $\bigstar\bigstar\bigstar\bigstar\bigstar$ \\ \hline
1 & 53.25 & 52 & 1.245 & 10.2 & 0.12 & 92 & 40 & $\bigstar\bigstar\bigstar\bigstar\bigstar$ \\ \hline
2 & 106.4 & 128 & -21.64 & 16 & -1.35 & 168 & 40 & $\bigstar\bigstar\bigstar\bigstar\bigstar$ \\ \hline
4 & 353.9 & 401 & -47.15 & 28.32 & -1.66 & 441 & 40 & $\bigstar\bigstar\bigstar\bigstar\bigstar$ \\ \hline
8 & 690.8 & 777 & -86.18 & 39.42 & -2.19 & 817 & 40 & $\bigstar\bigstar\bigstar\bigstar\bigstar$ \\ \hline
\end{tabular}

\caption{\textbf{Comparison between the computed and expected $\chi^2$ for different values of $L$}.  $N_S$ and $N_p$ are the number of gate strings and parameters, respectively.  $\chi^2$ is expected to lie within $[k-\sqrt{2k},k+\sqrt{2k}]$ where $k = N_s-N_p$. $P$ is the p-value derived from a $\chi^2_k$ distribution (i.e.~when $P <$ some threshold like $0.05$ there is grounds for rejecting the GST fit). The rating from 1 to 5 stars gives a very crude indication of goodness of fit as explained in the text.\label{progressTable}}
\end{center}
\end{table}


%\begin{figure}
%\begin{center}
%\XXputfield{bestEstimateColorBoxPlot}{Color box plot for best gateset}
%\caption{\XXputfield{tt_bestEstimateColorBoxPlot}{}\textbf{\XXputfield{objective}{??} values for every individual experiment in the dataset}.  Each pixel represents a single experiment (gate sequence), and its color indicates whether GST was able to fit the corresponding frequency well.  Blue is typical; dark red squares indicating $\chi^2>10$ should appear only once per 638 experiments on average.  Square blocks of pixels correspond to base sequences (arranged vertically by germ and horizontally by length); each pixel within a block corresponds to a specific choice of pre- and post-fiducial sequences.  See text for further details.\label{bestEstimateColorBoxPlot}}
%\end{center}
%\end{figure}

}{ (You need to supply the GST Results object with more data in order to compute a goodness-of-fit section.) }

\FloatBarrier

\section{Gateset Estimates}
Below is a collection of GST's primary result tables, ordered in what we think is a natural way.  These tables all pertain to the best-estimate gateset, that is, the one corresponding to the final iteration of the GST algorithm.

\subsection{Derived Quantities}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Gate & \begin{tabular}{c}Process\\Infidelity\end{tabular} & \begin{tabular}{c}$\nicefrac{1}{2}$ Trace\\Distance\end{tabular} & $\nicefrac{1}{2}$ $\Diamond$-Norm \\ \hline
Gi & $ \begin{array}{c} 0.04 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.04 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.04 \\ \pm 0.002 \end{array} $ \\ \hline
Gx & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ \\ \hline
Gy & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ & $ \begin{array}{c} 0.038 \\ \pm 0.002 \end{array} $ \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|}
\hline
Gate & Error Generator \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
5\e{-7} & 3\e{-4} & 3\e{-4} & 8\e{-5} \\ 
-0.001 & -0.057 & -0.007 & -0.009 \\ 
0.001 & 0.001 & -0.054 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
-8\e{-6} & 7\e{-5} & 4\e{-5} & 5\e{-5} \\ 
4\e{-4} & -0.051 & -0.005 & -0.006 \\ 
3\e{-4} & 4\e{-4} & -0.052 & -0.002 \\ 
-2\e{-4} & -3\e{-4} & -4\e{-5} & -0.052
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
2\e{-5} & -5\e{-5} & -2\e{-5} & 3\e{-5} \\ 
2\e{-4} & -0.052 & 0.007 & 5\e{-4} \\ 
3\e{-4} & 0.001 & -0.05 & 0.001 \\ 
2\e{-4} & -7\e{-6} & -0.004 & -0.053
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{Comparison of GST estimated gates to target gates}.  This table presents, for each of the gates, three different measures of distance or discrepancy from the GST estimate to the ideal target operation.  See text for more detail.  The second table lists the ``Error Generator'' for each gate, which is the Lindbladian $\mathbb{L}$ that describes \emph{how} the gate is failing to match the target.  This error generator is defined by the equation $\hat{G} = G_{\mathrm{target}}e^{\mathbb{L}}$. \label{bestGatesetVsTargetTable}}
\end{center}
\end{table}

\begin{table}[h]
\small
\begin{center}
\begin{tabular}[l]{|c|c|c|c|c|c|}
\hline
Gate & Eigenvalues & Fixed pt & Rotn. axis & Diag. decay & Off-diag. decay \\ \hline
Gi & $ \begin{array}{c}
1 \\ 
0.949 \\ 
0.946e^{i0.00} \\ 
0.946e^{-i0.00}
 \end{array} $
 $\pm$ $ \begin{array}{c}
0.001 \\ 
0.013 \\ 
0.009e^{i0.68} \\ 
0.009e^{i0.68}
 \end{array} $
 & $ \begin{array}{c}
0.999 \\ 
-0.024 \\ 
0.011 \\ 
0.023
 \end{array} $
 & $ \begin{array}{c}
0.004 \\ 
-0.934 \\ 
0.149 \\ 
0.325
 \end{array} $
 & $ \begin{array}{c} 0.051 \\ \pm 0.013 \end{array} $ & $ \begin{array}{c} 0.054 \\ \pm 0.007 \end{array} $ \\ \hline
Gx & $ \begin{array}{c}
0.949e^{i1.57} \\ 
0.949e^{-i1.57} \\ 
1 \\ 
0.951
 \end{array} $
 $\pm$ $ \begin{array}{c}
0.003e^{i0.72} \\ 
0.003e^{i0.72} \\ 
0.001 \\ 
0.004
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
0.008 \\ 
-2\e{-5} \\ 
2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
0.001 \\ 
-1 \\ 
2\e{-5} \\ 
-3\e{-4}
 \end{array} $
 & $ \begin{array}{c} 0.049 \\ \pm 0.004 \end{array} $ & $ \begin{array}{c} 0.051 \\ \pm 0.002 \end{array} $ \\ \hline
Gy & $ \begin{array}{c}
1 \\ 
0.951 \\ 
0.949e^{i1.57} \\ 
0.949e^{-i1.57}
 \end{array} $
 $\pm$ $ \begin{array}{c}
0.001 \\ 
0.004 \\ 
0.003e^{i0.79} \\ 
0.003e^{i0.79}
 \end{array} $
 & $ \begin{array}{c}
1 \\ 
-8\e{-6} \\ 
0.007 \\ 
-2\e{-4}
 \end{array} $
 & $ \begin{array}{c}
-5\e{-4} \\ 
0.005 \\ 
-1 \\ 
0.001
 \end{array} $
 & $ \begin{array}{c} 0.049 \\ \pm 0.004 \end{array} $ & $ \begin{array}{c} 0.051 \\ \pm 0.002 \end{array} $ \\ \hline
\end{tabular}

\vspace{2em}
\begin{tabular}[l]{|c|c|c|c|c|}
\hline
\multirow{2}{*}{Gate} & \multirow{2}{*}{Angle} & \multicolumn{3}{c|}{Angle between Rotation Axes} \\ \cline{3-5}
 & & Gi & Gx & Gy \\ \hline
Gi & $ \begin{array}{c}(0.001 \\ \pm 0.002)\pi \end{array} $ &  & $ \begin{array}{c}(0.117 \\ \pm 0.336)\pi \end{array} $ & $ \begin{array}{c}(0.549 \\ \pm 0.426)\pi \end{array} $ \\ \hline
Gx & $ \begin{array}{c}(0.5 \\ \pm 0.001)\pi \end{array} $ & $ \begin{array}{c}(0.117 \\ \pm 0.336)\pi \end{array} $ &  & $ \begin{array}{c}(0.502 \\ \pm 0.001)\pi \end{array} $ \\ \hline
Gy & $ \begin{array}{c}(0.5 \\ \pm 0.001)\pi \end{array} $ & $ \begin{array}{c}(0.549 \\ \pm 0.426)\pi \end{array} $ & $ \begin{array}{c}(0.502 \\ \pm 0.001)\pi \end{array} $ &  \\ \hline
\end{tabular}

\caption{\textbf{Eigen-decomposition of estimated gates}.  Each estimated gate is described in terms of: (1) the eigenvalues of the superoperator; (2) the gate's fixed point (as a vector in $\mathcal{B}(\mathcal{H})$, in the Pauli basis); (3)  the axis around which it rotates, as a vector in $\mathcal{B}(\mathcal{H})$; (4) the angle of the rotation that it applies; (5) the decay rate along the axis of rotation (``diagonal decay''); (6) the decay rate perpendicular to the axis of rotation (``off-diagonal decay''); and (7) the angle between each gate's rotation axis and the rotation axes of the other gates.  ``X'' indicates that the decomposition failed or couldn't be interpreted. \label{bestGatesetDecompTable}}
\end{center}
\end{table}

\FloatBarrier

\subsection{Raw Estimates}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|c|}
\hline
Operator & Matrix & Hilbert-Schmidt vector (Pauli basis) & 95\% C.I. $\nicefrac{1}{2}$-width \\ \hline
$\rho_{0}$ & $ \left(\!\!\begin{array}{cc}
0.99 & 0.001e^{-i0.13} \\ 
0.001e^{i0.13} & 0.01
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.707 \\ 
0.002 \\ 
2\e{-4} \\ 
0.693
 \end{array} $
 & $ \begin{array}{c}
0.003 \\ 
0.004 \\ 
0.005 \\ 
0.004
 \end{array} $
 \\ \hline
$E_{0}$ & $ \left(\!\!\begin{array}{cc}
0.083 & 0.002e^{i1.84} \\ 
0.002e^{-i1.84} & 0.918
 \end{array}\!\!\right) $
 & $ \begin{array}{c}
0.708 \\ 
-0.001 \\ 
-0.002 \\ 
-0.591
 \end{array} $
 & $ \begin{array}{c}
0.001 \\ 
0.003 \\ 
0.004 \\ 
0.004
 \end{array} $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the SPAM operations}.  Compare to Table \ref{targetSpamTable}.\label{bestGatesetSpamTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
 & $E_{0}$ & $E_C$ \\ \hline
$\rho_{0}$ & $ \begin{array}{c} 0.091 \\ \pm 0.004 \end{array} $ & $ \begin{array}{c} 0.909 \\ \pm 0.006 \end{array} $ \\ \hline
\end{tabular}

\caption{\textbf{GST estimate of SPAM probabilities}.  Computed by taking the dot products of vectors in Table \ref{bestGatesetSpamTable}.\label{bestGatesetSpamParametersTable}}
\end{center}
\end{table}

\begin{table}[h]
\begin{center}
\begin{tabular}[l]{|c|c|c|}
\hline
Gate & Superoperator (Pauli basis) & 95\% C.I. $\nicefrac{1}{2}$-width \\ \hline
Gi & $ \left(\!\!\begin{array}{cccc}
1 & 3\e{-4} & 3\e{-4} & 7\e{-5} \\ 
-0.001 & 0.945 & -0.006 & -0.009 \\ 
0.001 & 0.001 & 0.947 & 0.004 \\ 
0.001 & 2\e{-5} & 0.001 & 0.949
 \end{array}\!\!\right) $
 & $ \left(\!\!\begin{array}{cccc}
0.001 & 0.002 & 0.002 & 0.001 \\ 
0.002 & 0.004 & 0.004 & 0.004 \\ 
0.002 & 0.004 & 0.004 & 0.004 \\ 
0.002 & 0.004 & 0.004 & 0.004
 \end{array}\!\!\right) $
 \\ \hline
Gx & $ \left(\!\!\begin{array}{cccc}
1 & 7\e{-5} & 4\e{-5} & 4\e{-5} \\ 
4\e{-4} & 0.951 & -0.005 & -0.005 \\ 
2\e{-4} & 3\e{-4} & 4\e{-5} & -0.949 \\ 
3\e{-4} & 3\e{-4} & 0.949 & -0.002
 \end{array}\!\!\right) $
 & $ \left(\!\!\begin{array}{cccc}
0.001 & 0.001 & 0.004 & 0.004 \\ 
0.002 & 0.004 & 0.001 & 0.002 \\ 
0.001 & 0.003 & 0.003 & 0.002 \\ 
0.001 & 0.003 & 0.002 & 0.003
 \end{array}\!\!\right) $
 \\ \hline
Gy & $ \left(\!\!\begin{array}{cccc}
1 & -5\e{-5} & -2\e{-5} & 3\e{-5} \\ 
2\e{-4} & -9\e{-6} & -0.004 & 0.949 \\ 
3\e{-4} & 0.001 & 0.951 & 0.001 \\ 
-2\e{-4} & -0.949 & -0.006 & -4\e{-4}
 \end{array}\!\!\right) $
 & $ \left(\!\!\begin{array}{cccc}
0.001 & 0.001 & 0.001 & 0.001 \\ 
0.001 & 0.003 & 0.001 & 0.002 \\ 
0.002 & 0.003 & 0.004 & 0.003 \\ 
0.001 & 0.002 & 0.002 & 0.002
 \end{array}\!\!\right) $
 \\ \hline
\end{tabular}

\caption{\textbf{The GST estimate of the logic gate operations}.  Compare to Table \ref{targetGatesTable}.\label{bestGatesetGatesTable}}
\end{center}
\end{table}

\end{document}