Constraining new physics in b → cℓν transitions

B decays proceeding via b → cℓν transitions with ℓ = e or μ are tree-level processes in the Standard Model. They are used to measure the CKM element Vcb, as such forming an important ingredient in the determination of e.g. the unitarity triangle; hence the question to which extent they can be affected by new physics contributions is important, specifically given the long-standing tension between Vcb determinations from inclusive and exclusive decays and the significant hints for lepton flavour universality violation in b → cτ ν and b → sℓℓ decays. We perform a comprehensive model-independent analysis of new physics in b → cℓν, considering all combinations of scalar, vector and tensor interactions occuring in single-mediator scenarios. We include for the first time differential distributions of B → D∗ℓν angular observables for this purpose. We show that these are valuable in constraining non-standard interactions. Specifically, the zero-recoil endpoint of the B → Dℓν spectrum is extremely sensitive to scalar currents, while the maximum-recoil endpoint of the B → D∗ℓν spectrum with transversely polarized D∗ is extremely sensitive to tensor currents. We also quantify the room for e-μ universality violation in b → cℓν transitions, predicted by some models suggested to solve the b → cτ ν anomalies, from a global fit to B → Dℓν and B → D∗ℓν for the first time. Specific new physics models, corresponding to all possible tree-level mediators, are also discussed. As a side effect, we present Vcb determinations from exclusive B decays, both with frequentist and Bayesian statistics, leading to compatible results. The entire numerical analysis is based on open source code, allowing it to be easily adapted once new data or new form factors become available.