Commit a5d7f860 authored by Mohammad Akhlaghi's avatar Mohammad Akhlaghi

Included the difference in distance plots for each filter

The difference in distance plots are now also plotted in a similar way
to the differences in magnitude. To do this generically and easily,
the $(two-d-hists) rule in `statistics.mk' was slightly modified, so
all the two-d histograms are generated with separate invokations of a
single AWK command. This greatly helps later corrections/additions.
parent 506e81d5
......@@ -193,4 +193,96 @@ similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\caption{\label{udf10f850lp} UDF10 F850LP comparison}
\end{figure}
\newpage
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f606w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure8.pdf}
\fi
\caption{\label{udff606w} UDF F606W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f775w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure9.pdf}
\fi
\caption{\label{udff775} UDF F775W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f814w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure10.pdf}
\fi
\caption{\label{udff814w} UDF F814W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f850lp-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure11.pdf}
\fi
\caption{\label{udff850lp} UDF F850LP comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f606w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure12.pdf}
\fi
\caption{\label{udf10f606w} UDF10 F606W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f775w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure13.pdf}
\fi
\caption{\label{udf10f775} UDF10 F775W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f814w-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure14.pdf}
\fi
\caption{\label{udf10f814w} UDF10 F814W comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f850lp-dist}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure15.pdf}
\fi
\caption{\label{udf10f850lp} UDF10 F850LP comparison}
\end{figure}
\end{document}
......@@ -30,23 +30,40 @@
twodxmin=20
twodxmax=30
twodxnbins=40
twodymin=-1.8125
twodymax=1.8125
twodynbins=42
twodmagmin=-1.8125
twodmagmax=1.8125
twoddistmin=0
twoddistmax=0.5
twoddir = $(BDIR)/tex/2D-histograms
two-d-hists = $(foreach field, udf udf10, \
$(foreach filter, $(filters), \
$(twoddir)/$(field)-$(filter).txt) )
$(two-d-hists): $(twoddir)/%.txt: $(ccatdir)/%.txt \
reproduce/scripts/two-dim-hist.awk \
reproduce/scripts/checks.awk | $(twoddir)
# Pull out the two columns to generate the histogram
awk '!/^#/{print $$5, $$8}' $(ccatdir)/$*.txt > $(twoddir)/$*_fh.txt
# Generate the 2D histogram.
awk -vxmin=$(twodxmin) -vxmax=$(twodxmax) -vymin=$(twodymin) \
-vymax=$(twodymax) -vxnumbins=$(twodxnbins) \
two-d-hists = $(foreach field, udf udf10, \
$(foreach filter, $(filters), \
$(foreach type, mag dist, \
$(twoddir)/$(field)-$(filter)-$(type).txt) ) )
$(two-d-hists): $(twoddir)/%.txt: \
$(ccatdir)/$$(word 1, $$(subst -, , $$*))-$$(word 2, $$(subst -, , $$*)).txt reproduce/scripts/two-dim-hist.awk reproduce/scripts/checks.awk \
| $(twoddir)
# See if the magnitude is required or the distance, then put
# the columns necessary for a 2D histogram in a respective
# temporary file. For the distance, the columns are in units
# of degrees, so to make the numbers more understandable, we
# will convert them to arcseconds.
tmp=$(twoddir)/$*_fh.txt; \
mord=$(word 3, $(subst -, , $*)); \
if [ $$mord = mag ]; then \
ymin=$(twodmagmin); \
ymax=$(twodmagmax); \
awk '!/^#/{print $$5, $$8}' $< > $$tmp; \
else \
ymin=$(twoddistmin); \
ymax=$(twoddistmax); \
awk '!/^#/{print $$5, 3600*sqrt( ($$6*$$6)+($$7*$$7) )}' \
$< > $$tmp; \
fi; \
\
awk -vxmin=$(twodxmin) -vxmax=$(twodxmax) -vymin=$$ymin \
-vymax=$$ymax -vxnumbins=$(twodxnbins) \
-vynumbins=$(twodynbins) -vcol=8 -vshowemptylast=1 \
-vextraout=0 -i reproduce/scripts/checks.awk \
-f reproduce/scripts/two-dim-hist.awk $(twoddir)/$*_fh.txt > $@
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f606w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f606w.txt};
] file {\buildtexdir/2D-histograms/udf-f606w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f775w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f775w.txt};
] file {\buildtexdir/2D-histograms/udf-f775w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f814w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f814w.txt};
] file {\buildtexdir/2D-histograms/udf-f814w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f850lp-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f850lp.txt};
] file {\buildtexdir/2D-histograms/udf-f850lp-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f606w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f606w.txt};
] file {\buildtexdir/2D-histograms/udf10-f606w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f775w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f775w.txt};
] file {\buildtexdir/2D-histograms/udf10-f775w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f814w-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f814w.txt};
] file {\buildtexdir/2D-histograms/udf10-f814w-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
\pgfplotsset{
/pgfplots/colormap={coldredux}{
[1cm]
rgb255(0cm)=(255,255,255)
rgb255(2cm)=(0,192,255)
rgb255(4cm)=(0,0,255)
rgb255(6cm)=(0,0,0)
}
}
\small
\begin{tikzpicture}
\begin{axis}[
width=\linewidth,
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
]
\addplot3[
surf,
shader=flat corner,
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f850lp-dist.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf10-f850lp.txt};
] file {\buildtexdir/2D-histograms/udf10-f850lp-mag.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment