We study the shape of eight dynamically relaxed galaxy clusters observed with the Hubble Space Telescope and Chandra X-Ray Observatory. Using strong and weak gravitational lensing, the shape of the Brightest Cluster Galaxy and the X-ray isophote we study the ellipticity of the cluster halo at four different radii. We find that the proxies probing the inner regions of the cluster are strongly coupled with the BCG shape correlated with both the shape predicted by strong gravitational lensing and the X-ray isophote. Conversely we find no such correlation between the shape as predicted by the weak lensing and the other three probes suggesting any coupling between the inner and outer regions is weak. We also present in this paper the public release of the HST weak lensing shape measurement code pyRRG, directly available from PyPi (https://pypi.org/project/pyRRG/). This python3.7 code, based on Rhodes et al. (2000) adopts an automated star-galaxy classifier based on a Random Forest and outputs scientifically useful products such as weak lensing catalogues suitable for the mass mapping algorithm Lenstool.