One plus two-body embedded Gaussian orthogonal ensemble of random matrices with parity [EGOE(1+2)-π] generated by a random two-body interaction (modeled by GOE in two particle spaces) in the presence of a mean-field, for spinless identical fermion systems, is defined, generalizing the two-body ensemble with parity analyzed by Papenbrock and Weidenmüller [Phys. Rev. C 78, 054305 (2008)], in terms of two mixing parameters and a gap between the positive (π = +) and negative (π = −) parity single particle (sp) states. Numerical calculations are used to demonstrate, using realistic values of the mixing parameters appropriate for some nuclei, that the EGOE(1+2)π ensemble generates Gaussian form (with corrections) for fixed parity eigenvalue densities (i.e. state densities). The random matrix model also generates many features in parity ratios of state densities that are similar to those predicted by a method based on the Fermi-gas model for nuclei.We have also obtained, by applying the formulation due to Chang et al [Ann. Phys. (N.Y.) 66, 137 (1971)], a simple formula for the spectral variances defined over fixed-(m 1 , m 2 ) spaces, where m 1 is the number of fermions in the +ve parity sp states and m 2 is the number of fermions in the −ve parity sp states. Similarly, using the binary correlation approximation, in the dilute limit, we have derived expressions for the lowest two shape parameters. The smoothed densities generated by the sum of fixed-(m 1 , m 2 ) Gaussians with lowest two shape corrections describe the numerical results in many situations. The model also generates preponderance of +ve parity ground states for small values of the mixing parameters and this is a feature seen in nuclear shell model results.