In this article we develop new techniques for studying concentration of Laplace eigenfunctions φ λ as their frequency, λ, grows. The method consists of controlling φ λ (x) by decomposing φ λ into a superposition of geodesic beams that run through the point x. Each beam is localized in phase-space on a tube centered around a geodesic whose radius shrinks slightly slower than λ − 1 2 . We control φ λ (x) by the L 2 -mass of φ λ on each geodesic tube and derive a purely dynamical statement through which φ λ (x) can be studied. In particular, we obtain estimates on φ λ (x) by decomposing the set of geodesic tubes into those that are non self-looping for time T and those that are. This approach allows for quantitative improvements, in terms of T , on the available bounds for L ∞ norms, L p norms, pointwise Weyl laws, and averages over submanifolds.