We study the massive scalar field Sorkin-Johnston (SJ) Wightman function WSJ restricted to a flat 2D causal diamond D of linear dimension L. Our approach is two-pronged. In the first, we solve the central SJ eigenvalue problem explicitly in the small mass regime, up to order (mL) 4 . This allows us to formally construct WSJ up to this order. Using a combination of analytical and numerical methods, we obtain expressions for WSJ both in the center and the corner of D, to leading order. We find that in the center, WSJ is more like the massless Minkowski Wightman function W mink 0 than the massive one W mink m , while in the corner it corresponds to that of the massive mirror W mirror m . In the second part, in order to explore larger masses, we perform numerical simulations using a causal set approximated by a flat 2D causal diamond. We find that in the center of the diamond the causal set SJ Wightman function W c SJ resembles W mink 0 for small masses, as in the continuum, but beyond a critical value mc it resembles W mink m , as expected. Our calculations suggest that unlike W mink m , WSJ has a well-defined massless limit, which mimics the behavior of the Pauli Jordan function underlying the SJ construction. In the corner of the diamond, moreover, W c SJ agrees with W mirror m for all masses, and not, as might be expected, with the Rindler vacuum. * abhishekmathur@rri.res.in arXiv:1906.07952v2 [hep-th]