(*start*)
f[t_] = D[RiemannSiegelTheta[t], t];
nnn = 60
cc = 10^-10;
g1 = Plot[(f[t] + cc + EulerGamma), {t, 0, nnn},
PlotStyle -> {Thickness[0.004], Red}, PlotRange -> {-2, cc + 5}];
c = 1 + 1/cc;
g2 = Plot[
Re[Zeta[1/2 + I*t]*Zeta[c]/Zeta[1/2 + I*t + c - 1]], {t, 0, nnn},
PlotRange -> {-2, cc + 5}, PlotStyle -> Thickness[0.004]];
Show[g2, g1, ImageSize -> Large]
(*end*)