Even though Planck data released in 2013 (P13) is not compatible with Background Imaging of Cosmic Extragalactic Polarization (B2) and some local cosmological observations, including Supernova Legacy Survey (SNLS) samples and H0 prior from Hubble Space Telescope (HST) etc, Wilkinson Microwaves Anisotropy Probe 9-year data (W9) is consistent with all of them in the base six-parameter ΛCDM+tensor cosmology quite well. In this letter, we adopt the combinations of B2+W9 and B2+W9+SNLS+BAO+HST to constrain the cosmological parameters in the base six-parameter ΛCDM+tensor model with nt = −r/8, where r and nt are the tensor-to-scalar ratio and the tilt of relic gravitational wave spectrum, and BAO denotes Baryon Acoustic Oscillation. We find that the Harrison-Zel'dovich (HZ) scale invariant scalar power spectrum is consistent with both data combinations, chaotic inflation is marginally disfavored by the data at around 2σ level, but the power-law inflation model and the inflation model with inverse power-law potential can fit the data nicely.PACS numbers: 98.70. Vc,98.80.Cq, More than thirty years ago inflation [1][2][3] was proposed to solve the puzzles, such as the flatness problem, horizon problem, monopole problem and so on, in the hot big bang model. In fact, the spatial flatness can be taken as a prediction of inflation model which has been confirmed by Wilkinson Microwaves Anisotropy Probe 9-year data (W9) [4] and Planck data released in 2013 (P13) [5]. On the other hand, the quantum fluctuations generated during inflation [6][7][8][9][10][11] can seed the anisotropies in the cosmic microwaves background (CMB) radiation and the formation of large-scale structure. Since the Hubble parameter during inflation is roughly a constant, the spectrum of scalar perturbations is nearly scale-invariant. An adiabatic, Gaussian and nearly scale-independent scalar power spectrum has also been confirmed by W9 and P13. In addition, the quantization of the gravitational field during inflation produces a primordial background of stochastic gravitational waves [12][13][14][15][16]. In the last decades, many group tried their best to hunt for the signal of relic gravitational waves.Recently discovery of relic gravitational waves was reported by Background Imaging of Cosmic Extragalactic Polarization (B2) [17], and the tensor-to-scalar ratio is given by r = 0.20with r = 0 disfavored at 7.0σ. It is certainly a breakthrough of basic science in these years. Before B2 claimed its discovery, some hints of relic gravitational waves around r ∼ 0.2 were illustrated in [18,19] 2 ), angular scale of the sound horizon at last-scattering (θ MC ), optical depth (τ ), scalar spectrum power-law index (n s ) and log power of the primordial curvature perturbations (ln(10 10 A s )). Including the perturbations of primordial gravitational waves, the six-parameter ΛCDM model is extended to be ΛCDM+tensor model. However, combining with WMAP Polarization data [4], ACT [21] and SPT [22], P13 [5] imply a much smaller tensor-to-scalar ratio, compared to tha...