Bose-Einstein correlations of identical hadrons reveal information about hadron creation from the sQGP formed in ultrarelativistic heavy ion collisions. The measurement of three particle correlations may in particular shed light on hadron creation mechanisms beyond thermal/chaotic emission. In this talk we show the status of PHENIX measurements of three pion correlations as a function of momentum differences within the triplets. We analyze their shape through the assumption of Levy sources and a proper treatment of the Coulomb interaction within the triplets.We determine Levy parameters $R$, $\alpha$ and $\lambda_3$, which give source size, shape and correlation strength, respectively. The $\lambda_3$ encodes information about hadron creation mechanism, together with two particle correlation strength $\lambda_2$. From a consistent analysis of two- and three-particle correlation strength we may be able to establish an experimental measure of thermalization and coherence in the source.