If the state of one qubit can be described by a ray in $\mathbb{C}^2$, then the combined state of an $n$-qubit system can be described by a ray in $(\mathbb{C}^2)^{\otimes n}=\mathbb{C}^{2 n}$.
However, if $G_1$ is the Pauli group of one qubit, with the 16 elements $$G_1=\{i,-1,-i,1\}\times\{I,X,Y,Z\}\,,$$ the Pauli group on $n$ qubits is defined by $$G_n=\{i,-1,-i,1\}\times\{I,X,Y,Z\}^{\otimes n}$$ which not the tensor product of $n$ Pauli groups $G_1$ (because $G_n$ contains $4\cdot 4^n$ elements, which does not equal $16^n$). My question thus is: what kind of tensor product do we use on the space of operators on a Hilbert space $\mathbb{C}^2$, to define $G_n$ using $G_1$?
(I do understand intuitively that we should disregard global phase (and that therefore the number of operators in the $n$-qubit Pauli group is not $16^n$), and that this can be done by introducing the projective Hilbert space, but how do tensor products work on the space of operators on a projective Hilbert space?)