p1 = Plot3D[(3*E^(-(x^2))) + 0.05, {x, -4, 4}, {y, 0, 2},
   PlotRange -> {-1, 4}, RegionFunction -> Function[{x, z}, x < 0],
   Boxed -> False, Axes -> False, BoundaryStyle -> Yellow,
   PlotStyle -> RGBColor[1, 1, 0]];

p2 = Plot3D[4.05*E^(-(0.5*x^2)) - 1, {x, -4, 4}, {y, 0, 2},
   RegionFunction -> Function[{x, z}, x > 0], BoundaryStyle -> Yellow,
    PlotStyle -> RGBColor[1, 1, 0]];

p3 = Plot3D[E^(-(5*(x + 0.6)^2)), {x, -4, 4}, {y, 0, 2},
   RegionFunction -> Function[{x, z}, x < -0.4], Mesh -> None,
   Boxed -> False, BoundaryStyle -> Blue,
   PlotStyle -> RGBColor[0, 0, 1]];

p4 = Plot3D[0.5*E^(-(12.5*(x + 0.6)^2)) + 0.5, {x, -4, 4}, {y, 0, 2},
   RegionFunction -> Function[{x, z}, -0.4 < x < 0], Mesh -> None,
   BoundaryStyle -> Blue, PlotStyle -> RGBColor[0, 0, 1]];

p5 = Plot3D[0.5*E^(-(12.5*(x - 0.6)^2)) + 0.5, {x, -4, 4}, {y, 0, 2},
   RegionFunction -> Function[{x, z}, 0 <= x < 0.6], Mesh -> None,
   BoundaryStyle -> Blue, PlotStyle -> RGBColor[0, 0, 1]];

p6 = Plot3D[2*E^(-(2*(x - 0.6)^2)) - 1, {x, -4, 4}, {y, 0, 2},
   RegionFunction -> Function[{x, z}, x > 0.6], Mesh -> None,
   BoundaryStyle -> Blue, PlotStyle -> RGBColor[0, 0, 1]];

Show[p1, p2, p3, p4, p5, p6]

Graphics3D[hexTile[20, 20] /.
   Polygon[l_] :> {Directive[Orange, Opacity[0.8], Specularity[White, 30]],
     Polygon[l], Polygon[{Pi/5, 0} + {-1, 1} # & /@ l]} /.
  Polygon[l_List] :> Tube[top @@@ l],
 Boxed -> False, Axes -> False, PlotRange -> All,
 Lighting -> "Neutral", Background -> Black]

foo = {Cos[#], Sin[#] + Cos[#2], Sin[#2]} &;

Graphics3D[hexTile[20, 20] /.
  Polygon[l_List] :> {Directive[Orange, Opacity[0.8],
     Specularity[White, 30]], Tube[foo @@@ l]},
 Boxed -> False, Axes -> False, PlotRange -> All,
 Lighting -> "Neutral", Background -> Black]

pw1 = Piecewise[{{{3*E^(-(#^2)) + 0.05}, # <= 0},
     {{4.05*E^(-(0.5*#^2)) - 1}, # > 0}}] &;
pw2 = Piecewise[{{E^(-(5*(# + 0.6)^2)), # < -0.4},
     {0.5*E^(-(12.5*(# + 0.6)^2)) + 0.5, -0.4 <= # < 0},
     {0.5*E^(-(12.5*(# - 0.6)^2)) + 0.5, 0 <= # < 0.6},
     {2*E^(-(2*(# - 0.6)^2)) - 1, 0.6 <= #}}] &;

Plot3D[{pw1[x], pw2[x]}, {x, -4, 4}, {y, 0, 2},
 PlotRange -> All, Boxed -> False, Axes -> False,
 BoxRatios -> Automatic, PlotStyle -> {Yellow, Blue}, PlotPoints -> 80,
 Mesh -> None, Exclusions -> None, BoundaryStyle -> None]

top = Piecewise[{{{#, #2, 3*E^(-(#^2)) + 0.05}, # <= 0},
     {{#, #2, 4.05*E^(-(0.5*#^2)) - 1}, # > 0}}] &;
bottom = Piecewise[{{{#, #2, E^(-(5*(# + 0.6)^2))}, # < -0.4},
     {{#, #2, 0.5*E^(-(12.5*(# + 0.6)^2)) + 0.5}, -0.4 <= # < 0},
     {{#, #2, 0.5*E^(-(12.5*(# - 0.6)^2)) + 0.5}, 0 <= # < 0.6},
     {{#, #2, 2*E^(-(2*(# - 0.6)^2)) - 1}, 0.6 <= #}}] &;

ParametricPlot3D[{top[x, y], bottom[x, y]}, {x, -4, 4}, {y, 0, 2},
 PlotRange -> All, Boxed -> False, BoxRatios -> Automatic, Axes -> False,
 PlotStyle -> {Yellow, Blue}, PlotPoints -> 80, Mesh -> None, Exclusions -> None]

hexTile[n_, m_] := With[{hex = Polygon[Table[{Cos[2 Pi k/6] + #, Sin[2 Pi k/6] + #2},
  {k, 6}]] &},
        Table[hex[3 i + 3 ((-1)^j + 1)/4, Sqrt[3]/2 j], {i, n}, {j, m}] /.
            {x_?NumericQ, y_?NumericQ} :> 2 π {x/(3 m), 2 y/(n Sqrt[3])}]

gr1 = Graphics3D[hexTile[20, 20] /.
    Polygon[l_] :> {Yellow, Polygon[l], Polygon[{Pi/5, 0} + {-1, 1} # & /@ l]} /.
   Polygon[l_List] :> Polygon[top @@@ l],
  Boxed -> False, Axes -> False, PlotRange -> All, Lighting -> "Neutral"]

gr2 = Graphics3D[hexTile[20, 20] /.
    Polygon[l_] :> {Blue, Polygon[l], Polygon[{Pi/5, 0} + {-1, 1} # & /@ l]} /.
   Polygon[l_List] :> Polygon[bottom @@@ l],
  Boxed -> False, Axes -> False, PlotRange -> All, Lighting -> "Neutral"]

Show[gr1, gr2]