Case 2) Z = X+Y without any constraint

is easy, you're adding two independent gammas, use the MGF technique here

Case 1) Z = X-Y with the constraint that X is strictly higher than Y

Here you can draw the region.

I would get the joint density and just integrate the infinite triangle in the first quadrant

where X>Y