In this paper, we study the following asymptotically critical Choquard equation 0.1 \(\begin{aligned} \left\{ \begin{array}{ll} -\Delta u=\displaystyle \left( \int _{\Omega }\frac{u^{2_\alpha ^*-\varepsilon }(y)}{|x-y|^\alpha }dy\right) u^{2_\alpha ^*-1-\varepsilon }, \ \ & \text{ in }\ \Omega ,\\ u>0,\ \ & \text{ in }\ \Omega ,\\ u=0, \ \ & \text{ on }\ \partial \Omega , \end{array} \right. \end{aligned}\) where \(\Omega \) is a smooth bounded domain in \(\mathbb {R}^N\) , \(N\ge 3\) , \(\alpha \in (0,N)\) , \(2_\alpha ^*=\frac{2N-\alpha }{N-2}\) is the upper critical exponent in the sense of the Hardy–Littlewood–Sobolev inequality, and \(\varepsilon >0\) is a small parameter. We prove that any family of solutions to (0.1), under a suitable constrain on their gradients, blows up exactly at a critical point of the Robin function. Conversely, by a reduction argument, we construct a family of solutions to (0.1) concentrating around a critical point of the Robin function as \(\varepsilon \rightarrow 0\) .