The Airy$_β$ point process, $a_i \equiv N^{2/3} (λ_i-2)$, describes the eigenvalues $λ_i$ at the edge of the Gaussian $β$ ensembles of random matrices for large matrix size $N \to \infty$. We study the probability distribution function (PDF) of linear statistics ${\sf L}= \sum_i t \varphi(t^{-2/3} a_i)$ for large parameter $t$. We show the large deviation forms $\mathbb{E}_{{\rm Airy},β}[\exp(-{\sf L})] \sim \exp(- t^2 Σ[\varphi])$ and $P({\sf L}) \sim \exp(- t^2 G(L/t^2))$ for the cumulant generating function and the PDF. We obtain the exact rate function $Σ[\varphi]$ using four apparently different methods (i) the electrostatics of a Coulomb gas (ii) a random Schrödinger problem, i.e. the stochastic Airy operator (iii) a cumulant expansion (iv) a non-local non-linear differential Painlevé type equation. Each method was independently introduced to obtain the lower tail of the KPZ equation. Here we show their equivalence in a more general framework. Our results are obtained for a class of functions $\varphi$, the monotonous soft walls, containing the monomials $\varphi(x)=(u+x)_+^γ$ and the exponential $\varphi(x)=e^{u+x}$ and equivalently describe the response of a Coulomb gas pushed at its edge. The small $u$ behavior of the excess energy $Σ[\varphi]$ exhibits a change at $γ=3/2$ between a non-perturbative hard wall like regime for $γ<3/2$ (third order free-to-pushed transition) and a perturbative deformation of the edge for $γ>3/2$ (higher order transition). Applications are given, among them: (i) truncated linear statistics such as $\sum_{i=1}^{N_1} a_i$, leading to a formula for the PDF of the ground state energy of $N_1 \gg 1$ noninteracting fermions in a linear plus random potential (ii) $(β-2)/r^2$ interacting spinless fermions in a trap at the edge of a Fermi gas (iii) traces of large powers of random matrices.