AbstractWe analyze a random lozenge tiling model of a large regular hexagon, whose underlying weight structure is periodic of period 2 in both the horizontal and vertical directions. This is a determinantal point process whose correlation kernel is expressed in terms of non‐Hermitian matrix valued orthogonal polynomials (OPs). This model belongs to a class of models for which the existing techniques for studying asymptotics cannot be applied. The novel part of our method consists of establishing a connection between matrix valued and scalar valued OPs. This allows to simplify the double contour formula for the kernel obtained by Duits and Kuijlaars by reducing the size of a Riemann–Hilbert problem. The proof relies on the fact that the matrix valued weight possesses eigenvalues that live on an underlying Riemann surface of genus 0. We consider this connection of independent interest; it is natural to expect that similar ideas can be used for other matrix valued OPs, as long as the corresponding Riemann surface is of genus 0. The rest of the method consists of two parts, and mainly follows the lines of a previous work of Charlier, Duits, Kuijlaars and Lenells. First, we perform a Deift–Zhou steepest descent analysis to obtain asymptotics for the scalar valued OPs. The main difficulty is the study of an equilibrium problem in the complex plane. Second, the asymptotics for the OPs are substituted in the double contour integral and the latter is analyzed using the saddle point method. Our main results are the limiting densities of the lozenges in the disordered flower‐shaped region. However, we stress that the method allows in principle to rigorously compute other meaningful probabilistic quantities in the model.