The issue with attractive pairwise interactions is that the integral of the intensity of the process is infinite as each new point always increases the intensity for others in the neighbourhood.
If your process is mildly attractive on medium distances and strongly repulsive on very short ones (like physical particles), you can consider a pairwise interaction like Lennard-Jones or similar. See Van Lieshout (2000) Markov point processes and their applications or in a more applied fashion using R: Baddeley (2010) Analysing spatial point patterns in 'R'.
Else, the best way out is to use higher-order interactions such as area-interaction process (e.g. Baddeley and van Lieshout, 1995) or one of the processes in Geyer (1999) Likelihood inference for spatial point processes. Higher order interactions allow for both attractive and repulsive behaviour, e.g. points could be pairwise attractive, but triplets of points could be repulsive.
Depending on your application, the inhomogeneity of space could play a very important role. However, whether the intensity depends on the properties of underlying space or on the presence of other points should be dictated by theory. Even if these effects potentially could be hard to tell apart empirically.
