  ~In this paper we study a class of Fourier integral operators, whose symbols
lie in the multi-parameter H\"ormander class $S^{\vec m}( \mathbb{R}^{\vec
n})$, where ~$\vec m=(m_1,m_2,\dots,m_d)$ is the order. We show that if in
addition the phase function $\Phi(x,\xi)$ can be written as
$\Phi(x,\xi)=\sum_{i=1}^d\Phi_i(x_i,\xi_i)$, and each $\Phi_i(x_i,\xi_i)$
satisfies the non-degeneracy condition, then we show that such Fourier integral
operators with order ~$\vec m=(-(n_1-1)/2, -(n_2-1)/2,\dots, -(n_d-1)/2)$ are
actually bounded from rectangular Hardy space $H_{rect}^1(\mathbb{R}^{\vec n})$
to $L^1( \mathbb{R}^n )$. Moreover, they do not decay away from the support of
rectangle atoms. Hence the usual argument via a Journ\'e-type covering lemma to
deduce bounds on product Hardy space $H^1(\mathbb{R}^{\vec n})$ is not valid.
