We present the first comparisons of experimental data with phenomenological results from 3+1d quasiparticle anisotropic hydrodynamics (aHydroQP). We compare charged-hadron multiplicity, identified-particle spectra, identified-particle average transverse momentum, charged-particle elliptic flow, and identified-particle elliptic flow produced in LHC 2.76 TeV Pb+Pb collisions. The dynamical equations used for the hydrodynamic stage utilize non-conformal aHydroQP. The resulting aHydroQP framework naturally includes both shear and bulk viscous effects in addition to higherorder non-linear transport coefficients. The 3+1d aHydroQP evolution obtained is self-consistently converted to hadrons using anisotropic Cooper-Frye freezeout performed on a fixed-energy-density hypersurface. The final production and decays of the primordial hadrons are modeled using a customized version of THERMINATOR 2. In this first study, we utilized smooth Glauber-type initial conditions and a single effective freeze-out temperature TFO = 130 MeV with all hadronic species in full chemical equilibrium. With this rather simple setup, we find a very good description of many heavy-ion observables.PACS numbers: 12.38. Mh, 24.10.Nz, 25.75.Ld, 47.75.+f Ultrarelativistic heavy-ion collision experiments at the Relativistic Heavy Ion Collider and Large Hadron Collider (LHC) were designed to create and study the quarkgluon plasma (QGP). Relativistic hydrodynamics has been quite successful in describing the collective behavior observed in high-energy heavy-ion collisions [1-3] and the current focus of the relativistic hydrodynamics community is on further improvements of the models to include e.g. bulk viscous effects and higher-order transport coefficients (see [29][30][31] for recent reviews). The goal of the relativistic viscous hydrodynamics program is to constrain key properties of the QGP such as its initial energy density, initial pressure anisotropies, shear viscosity, bulk viscosity, etc. and to also provide the soft-background evolution necessary to compute QGPmodification of hard probes such as jets and heavy quark bound states.One of the issues faced by practitioners of traditional second-order viscous hydrodynamics approaches is that, at early times after the nuclear impact, the QGP possesses a high degree of momentum-space anisotropy in the fluid local rest frame, P T /P L 1. The magnitude of the resulting momentum-space anisotropy is large at early times after the initial nuclear impact and also near the transverse/longitudinal "edges" of the QGP at all times. In these spacetime regions, traditional viscous hydrodynamics is being pushed to its limits, resulting in potentially negative total pressures and violations of positivity of the one-particle distribution function [32].As a way to address these problems, it was suggested that one should reorganize the expansion of the one-particle distribution function around a leadingorder form which possesses intrinsic momentum-space anisotropies but still guarantees positivity [33,34]. This met...