YtojjzMatrix(l) takes as input the angular momentum $l$ and returns a rotation matrix between the $l_z,s_z$ and $j,j_z$ basis. The dimension of the matrix returned is $2*(2*l+1)$. The format of the matrix can be found in the documentation of the different one particle basis states.