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