Commit fd222a3a authored by Mohammad Akhlaghi's avatar Mohammad Akhlaghi

Plots for difference in RA and Dec created

The plots for the difference in RA and Dec as a function of magnitude
for the mosaic and the UDF10 have been created in the final PDF.
parent 21afff81
......@@ -75,20 +75,19 @@ The resulting plots for the UDF field (combined from all 9 subfields)
can be seen in Figures \ref{udff606w} to \ref{udff850lp} and those for
the deep UDF-10 field are shown in Figures \ref{udf10f606w} to
\ref{udf10f850lp}. The UDF-10 and mosiac magnitudes mostly agree with
each other. The F606W \new{and F814W} plot is mostly similar to a
similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\new{The degraded F814W image in the XDF survey has a strong bias in
its Sky value which dramatically affects the magnitudes of objects
below the 24th magnitude (Figure \ref{f814w-demo}). To correct this
residual, the Sky value generated by NoiseChisel during the creation
of segmentation maps has been subtracted from the image. \tonote{The
fact that the big problem was removed once we subtracted the Sky
value maybe one confirmation that the problem was in the XDF
images, not in our MUSE processing. If you can think of any other
test, I would be happy to apply it. This maybe an important issue
with the XDF survey and we should be very sure to confirm this
claim.}}
each other. The F606W and F814W plot is mostly similar to a similar
plot in the HDFS field, see Figure 10 in \citet{bacon15}.
The degraded F814W image in the XDF survey has a strong bias in its
Sky value which dramatically affects the magnitudes of objects below
the 24th magnitude (Figure \ref{f814w-demo}). To correct this
residual, the Sky value generated by NoiseChisel during the creation
of segmentation maps has been subtracted from the image. \tonote{The
fact that the big problem was removed once we subtracted the Sky
value maybe one confirmation that the problem was in the XDF images,
not in our MUSE processing. If you can think of any other test, I
would be happy to apply it. This maybe an important issue with the
XDF survey and we should be very sure to confirm this claim.}
\begin{figure}
\centering
......@@ -113,41 +112,41 @@ similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f606w}
\input{./tex/udf-f606w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure0.pdf}
\fi
\caption{\label{udff606w} UDF F606W comparison}
\caption{\label{udff606w} UDF F606W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f775w}
\input{./tex/udf-f775w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure1.pdf}
\fi
\caption{\label{udff775} UDF F775W comparison}
\caption{\label{udff775} UDF F775W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f814w}
\input{./tex/udf-f814w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure2.pdf}
\fi
\caption{\label{udff814w} UDF F814W comparison}
\caption{\label{udff814w} UDF F814W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f850lp}
\input{./tex/udf-f850lp-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure3.pdf}
\fi
\caption{\label{udff850lp} UDF F850LP comparison}
\caption{\label{udff850lp} UDF F850LP magnitude comparison}
\end{figure}
......@@ -156,41 +155,41 @@ similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f606w}
\input{./tex/udf10-f606w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure4.pdf}
\fi
\caption{\label{udf10f606w} UDF10 F606W comparison}
\caption{\label{udf10f606w} UDF10 F606W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f775w}
\input{./tex/udf10-f775w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure5.pdf}
\fi
\caption{\label{udf10f775} UDF10 F775W comparison}
\caption{\label{udf10f775} UDF10 F775W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f814w}
\input{./tex/udf10-f814w-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure6.pdf}
\fi
\caption{\label{udf10f814w} UDF10 F814W comparison}
\caption{\label{udf10f814w} UDF10 F814W magnitude comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f850lp}
\input{./tex/udf10-f850lp-mag}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure7.pdf}
\fi
\caption{\label{udf10f850lp} UDF10 F850LP comparison}
\caption{\label{udf10f850lp} UDF10 F850LP magnitude comparison}
\end{figure}
......@@ -205,41 +204,41 @@ similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f606w-dist}
\input{./tex/udf-f606w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure8.pdf}
\fi
\caption{\label{udff606w} UDF F606W comparison}
\caption{\label{udff606w} UDF F606W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f775w-dist}
\input{./tex/udf-f775w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure9.pdf}
\fi
\caption{\label{udff775} UDF F775W comparison}
\caption{\label{udff775} UDF F775W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f814w-dist}
\input{./tex/udf-f814w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure10.pdf}
\fi
\caption{\label{udff814w} UDF F814W comparison}
\caption{\label{udff814w} UDF F814W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f850lp-dist}
\input{./tex/udf-f850lp-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure11.pdf}
\fi
\caption{\label{udff850lp} UDF F850LP comparison}
\caption{\label{udff850lp} UDF F850LP RA comparison}
\end{figure}
......@@ -248,41 +247,133 @@ similar plot in the HDFS field, see Figure 10 in \citet{bacon15}.
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f606w-dist}
\input{./tex/udf10-f606w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure12.pdf}
\fi
\caption{\label{udf10f606w} UDF10 F606W comparison}
\caption{\label{udf10f606w} UDF10 F606W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f775w-dist}
\input{./tex/udf10-f775w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure13.pdf}
\fi
\caption{\label{udf10f775} UDF10 F775W comparison}
\caption{\label{udf10f775} UDF10 F775W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f814w-dist}
\input{./tex/udf10-f814w-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure14.pdf}
\fi
\caption{\label{udf10f814w} UDF10 F814W comparison}
\caption{\label{udf10f814w} UDF10 F814W RA comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f850lp-dist}
\input{./tex/udf10-f850lp-ra}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure15.pdf}
\fi
\caption{\label{udf10f850lp} UDF10 F850LP comparison}
\caption{\label{udf10f850lp} UDF10 F850LP RA comparison}
\end{figure}
\newpage
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f606w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure16.pdf}
\fi
\caption{\label{udff606w} UDF F606W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f775w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure17.pdf}
\fi
\caption{\label{udff775} UDF F775W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f814w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure18.pdf}
\fi
\caption{\label{udff814w} UDF F814W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf-f850lp-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure19.pdf}
\fi
\caption{\label{udff850lp} UDF F850LP Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f606w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure20.pdf}
\fi
\caption{\label{udf10f606w} UDF10 F606W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f775w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure21.pdf}
\fi
\caption{\label{udf10f775} UDF10 F775W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f814w-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure22.pdf}
\fi
\caption{\label{udf10f814w} UDF10 F814W Dec comparison}
\end{figure}
\begin{figure}
\centering
\ifdefined\makepdf
\input{./tex/udf10-f850lp-dec}
\else
\includegraphics[width=\linewidth]{./tikz/paper-figure23.pdf}
\fi
\caption{\label{udf10f850lp} UDF10 F850LP Dec comparison}
\end{figure}
\end{document}
......@@ -59,11 +59,14 @@ tex/pipeline.tex: $(foreach t, versions statistics, $(mtexdir)/$(t).tex)
# BibLaTeX references
# -------------------
#
# Before building the pdf, we want to build all the data-products, so
# fill in the targts into the data products variable.
#
# To build the BibLaTeX references, we need to run LaTeX once and then
# biber, the rule for building the final PDF will build the final PDF.
description.bbl: tex/ref.tex tex/pipeline.tex | tikz
data-products = $(radec2d) $(mag2d) $(deghst-demo)
echo; echo "Wants to make bib"; echo; exit 1
description.bbl: $(data-products) tex/ref.tex tex/pipeline.tex | tikz
if pdflatex -shell-escape -halt-on-error description.tex; then \
echo "LaTeX (in preparation for BibLaTeX) is done."; \
......@@ -81,20 +84,13 @@ description.bbl: tex/ref.tex tex/pipeline.tex | tikz
# Build final PDF
# ---------------
#
# Before building the pdf, we want to build all the data-products, so
# fill in the targts into the data products variable.
#
# Note that eventhough `tex/pipeline.tex' is in the `tex/' directory
# `tex/*' will not detect it necessarily because it might not be
# present (it is a target to be built by Make).
data-products = $(fullmosaic)
description.pdf: $(data-products) description.tex tex/* tex/pipeline.tex \
description.bbl | tikz
echo; echo "Wants to make PDF"; echo; exit 1
# Delete some outputs for TeX to rebuild (if needed)
# rm tikz/description-figure0*
......
......@@ -33,43 +33,45 @@ twodxnbins=40
twodynbins=42
twodmagmin=-1.8125
twodmagmax=1.8125
twoddistmin=0
twoddistmax=0.5
twoddistmin=-0.4
twoddistmax=0.4
twoddir = $(BDIR)/tex/2D-histograms
two-d-hists = $(foreach field, udf udf10, \
$(foreach filter, $(filters), \
$(foreach type, mag dist, \
$(twoddir)/$(field)-$(filter)-$(type).txt) ) )
$(two-d-hists): $(twoddir)/%.txt: \
$(mcatdir)/$$(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; \
scripts2d = reproduce/scripts/two-dim-hist.awk reproduce/scripts/checks.awk
mag2d = $(foreach field, udf udf10, \
$(foreach filter, $(filters), $(twoddir)/$(field)-$(filter)-mag.txt) )
$(mag2d): $(twoddir)/%-mag.txt: $(mcatdir)/%.txt $(scripts2d) | $(twoddir)
awk -vxcol=5 -vxmin=$(twodxmin) -vxmax=$(twodxmax) \
-vycol=8 -vymin=$(twodmagmin) -vymax=$(twodmagmax) \
-vxnumbins=$(twodxnbins) -vynumbins=$(twodynbins) \
-vshowemptylast=1 -vextraout=0 -ireproduce/scripts/checks.awk \
-freproduce/scripts/two-dim-hist.awk $< > $@;
radec2d = $(foreach field, udf udf10, \
$(foreach filter, $(filters), \
$(foreach d, ra dec, $(twoddir)/$(field)-$(filter)-$(d).txt) ) )
$(radec2d): $(twoddir)/%.txt: $(acatdir)/$$(word 1,$$(subst -, , $$*))-$$(word 2,$$(subst -, , $$*)).txt $(scripts2d) | $(twoddir)
# Pull out the magnitude and distance in arcseconds.
if [ $(word 3,$(subst -, , $*)) = ra ]; then \
awk '!/^#/{print $$5, 3600*$$6}' $< > $(@D)/$*-tmp.txt; \
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 > $@
# Clean up
rm $(twoddir)/$*_fh.txt
awk '!/^#/{print $$5, 3600*$$7}' $< > $(@D)/$*-tmp.txt; \
fi
# Make the 2D histogram.
awk -vxcol=1 -vxmin=$(twodxmin) -vxmax=$(twodxmax) \
-vycol=2 -vymin=$(twoddistmin) -vymax=$(twoddistmax) \
-vxnumbins=$(twodxnbins) -vynumbins=$(twodynbins) \
-vshowemptylast=1 -vextraout=0 -ireproduce/scripts/checks.awk \
-freproduce/scripts/two-dim-hist.awk $(@D)/$*-tmp.txt > $@
# Clean up.
rm $(@D)/$*-tmp.txt
......
......@@ -31,13 +31,14 @@
# Variables:
# ==========
#
# xcol: Column for the X axis.
# xmin: Minimum value for the histogram.
# xmax: Maximum value for the histogram.
# xnumbins: Number of bins in the histogram.
# ycol: Column for the Y axis.
# ymin: Minimum value for the histogram.
# ymax: Maximum value for the histogram.
# ynumbins: Number of bins in the histogram.
# col: The column that should be used.
# showemptylast: If ==1, then show the last bins if
# they empty (with a zero value).
# extraout: Add an extra empty bin on each side.
......@@ -56,8 +57,10 @@
# Do the preparations:
BEGIN {
# Check if all the variables are set:
checkifset(xmax, "xcol")
checkifset(xmax, "xmax")
checkifset(xmin, "xmin")
checkifset(xmin, "ycol")
checkifset(xmax, "ymax")
checkifset(xmin, "ymin")
checkifint(xnumbins, "xnumbins")
......@@ -84,19 +87,17 @@ $0 !~ /^#/ {
# If the value is smaller or larger than the desired region, then
# go to the next record. This is not of interest.
if($1 < xmin || $1 > xmax || $2 < ymin || $2 > ymax) next
if($xcol < xmin || $xcol > xmax || $ycol < ymin || $ycol > ymax) next
# Find the middle point of the bin point for row.
xmidbin = xmin + int( ($1-xmin)/xwidth )*xwidth + xwidth/2
ymidbin = ymin + int( ($2-ymin)/ywidth )*ywidth + ywidth/2
# Find the middle point of the bin that this object falls in.
xmidbin = xmin + int( ($xcol-xmin)/xwidth )*xwidth + xwidth/2
ymidbin = ymin + int( ($ycol-ymin)/ywidth )*ywidth + ywidth/2
# If $1==xmax, or $2==ymax then the midbins will be larger than
# If $xcol==xmax, or $ycol==ymax then the midbins will be larger than
# the maximum value. In that case, we want it to be added to the
# last bin.
if(xmidbin>xmax) {
if($1==xmax) xmidbin-=xwidth
else next
}
if(xmidbin>xmax) { if($xcol==xmax) xmidbin-=xwidth; else next; }
if(ymidbin>ymax) { if($ycol==ymax) ymidbin-=ywidth; else next; }
# Convert the midbin value to a string, so floating point errors
# don't bother with the counting.
......
\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\delta$ (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-dec.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -16,7 +16,7 @@
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel=MUSE $-$ HST $\Delta\alpha$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f606w-dist.txt};
] file {\buildtexdir/2D-histograms/udf-f606w-ra.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\delta$ (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-dec.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
......@@ -16,7 +16,7 @@
view={0}{90},
colorbar horizontal,
xlabel=HST Magnitude,
ylabel=MUSE $-$ HST $\Delta{r}$ (arc seconds),
ylabel=MUSE $-$ HST $\Delta\alpha$ (arc seconds),
ylabel shift=-0.1cm,
colorbar style={xlabel={Number}, at={(0,1.1)}, anchor=north west,
xticklabel pos=upper},
......@@ -27,7 +27,7 @@
mesh/ordering=rowwise,
mesh/rows=\twodxnbins,
mesh/cols=\twodynbins,
] file {\buildtexdir/2D-histograms/udf-f775w-dist.txt};
] file {\buildtexdir/2D-histograms/udf-f775w-ra.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\delta$ (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-dec.txt};
\addplot[mark=none, thin, black] coordinates {
(\twodxmin,0) (\twodxmax,0)};
\end{axis}
\end{tikzpicture}
\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}
\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,