I would set up a new coordinate system in the plane determined by those three points.
For example, a vector from the point (2, 0, 0) to (0, 2, 0) is <2, -2, 0> so the line from (2, 0, 0) to (0, 2, 0) can be written as x= 2- 2u, y= 2u, z= 0 and we can take u to be one of the parameters for this new coordianate system. We also get that 2x- 2y= C for any plane perpendicular to that line and, in particular, 2x- 2y= 0 is a plane perpendicular to that line containing the point (0, 0, 1). The line x= 2- 2u, y= 2u, z= 0 intersects that plane where 2(2- 2u)- 2(2u)= 4- 8u= 0 or u= 1/2. That is the point (1, 1, 0) and so the line through (1, 1, 0) and (0, 0, 1) is a line in that plane perpendicular to the line. A vector from (1, 1, 0) to (0, 0, 1) is <1, 1, -1> and so the line is given by x= v+ 1, y= v+ 1, z= -v and we can take v as the other parameter.
That is, the line from (2, 0, 0) to (0, 2, 0) is the u-axis, v= 0, with u going from 0 to 1. The point (2, 0, 0) is given in uv coordines as (0, 0, (0, 2, 0) by (1, 0), and (0, 0, 1) by (1/2, 1). Integrate the density function (changed to u,v coordinates, of course) over the triangle with those coordinates.