Flavour physics in the LHC era

These lectures give a topical review of heavy flavour physics, in particular \CP violation and rare decays, from an experimental point of view. They describe the ongoing motivation to study heavy flavour physics in the LHC era, the current status of the field emphasising key results from previous experiments, some selected topics in which new results are expected in the near future, and a brief look at future projects.


Introduction
The concept of "flavour physics" was introduced in the 1970s [1] The term flavor was first used in particle physics in the context of the quark model of hadrons. It was coined in 1971 by Murray Gell-Mann and his student at the time, Harald Fritzsch, at a Baskin-Robbins ice-cream store in Pasadena. Just as ice cream has both color and flavor so do quarks.
Leptons also come in different flavours, and flavour physics covers the properties of both sets of fermions. Counting the fundamental parameters of the Standard Model (SM), the 3 lepton masses, 6 quark masses and 4 quark mixing (CKM) matrix [2,3] parameters are related to flavour physics. In case neutrino masses are introduced, the new parameters (at least 3 more masses and 4 more mixing parameters) are also related to flavour physics. This large number of free parameters is behind several of the mysteries of the SM: • Why are there so many different fermions?
• What is responsible for their organisation into generations / families? • What breaks the flavour symmetries?
Unfortunately these mysteries will not be answered in these lectures -they are mentioned here simply because it is important to bear in mind their existence. Instead the focus will be on specific topics in the flavour-changing interactions of the charm and beauty quarks, 1 with occasional digressions on related topics.
While our main interest is in the properties of the charm and beauty quarks, due to the strong interaction, experimental studies must be performed using one or more of the many different charmed or beautiful hadrons. These can decay to an even larger multitude of different final states, making learning the names of all the hadrons a big challenge for flavour physicists. Moreover, hadronic effects can often obscure the underlying dynamics. Nevertheless, it is the hadronisation that results in the very rich phenomenology that will be discussed, so one should bear in mind that [4] The strong interaction can be seen either as the "unsung hero" or the "villain" in the story of quark flavour physics.

Motivation to study heavy flavour physics in the LHC era
There are two main motivations for ongoing experimental investigations into heavy flavour physics: (i) CP violation and its connection to the matter-antimatter asymmetry of the Universe; (ii) discovery potential far beyond the energy frontier via searches for rare or SM forbidden processes. These will be discussed in turn below.
First let us consider one of the mysteries listed above (What breaks the flavour symmetries?) to see how it is connected to these motivations. In the SM, the vacuum expectation value of the Higgs field breaks the electroweak symmetry. Fermion masses arise from the Yukawa couplings of the quarks and charged leptons to the Higgs field, and the CKM matrix arises from the relative misalignment of the Yukawa matrices for the up-and down-type quarks. Consequently, the only flavourchanging interactions are the charged current weak interactions. This means that there are no flavour-changing neutral currents (the GIM mechanism [5]), a feature of the SM which is not generically true in most extended theories. Flavour-changing processes provide sensitive tests of this prediction; as an example, many new physics (NP) models induce contributions to the µ → eγ transition at levels close to (or even above!) the current experimental limit, recently made more restrictive by the MEG experiment [6], B(µ + → e + γ) < 5.7 × 10 −13 at 90% confidence level (CL). Improved experimental reach in this and related charged lepton flavour violation Flavour Physics in the LHC Era 3 searches therefore provides interesting and unique NP discovery potential (for a review, see, e.g., Ref. [7]).

CP violation
As mentioned above, the CKM matrix arises from the relative misalignment of the Yukawa matrices for the up-and down-type quarks: where U u and U d diagonalise the up-and down-type quark mass matrices respectively. Hence, V CKM is a 3 × 3 complex unitary matrix. Such a matrix is in general described by 9 (real) parameters, but 5 can be absorbed as unobservable phase differences between the quark fields. This leaves 4 parameters, of which 3 can be expressed as Euler mixing angles, but the fourth makes the CKM matrix complexand hence the weak interaction couplings differ for quarks and antiquarks, i.e. CP violation arises. The expression "CP violation" refers to the violation of the symmetry of the combined C and P operators, which replace particle with antiparticle (charge conjugation) and invert all spatial co-ordinates (parity) respectively. Therefore CP violation provides absolute discrimination between particle and antiparticle: one cannot simply swap the definition of which is called "particle" with a simultaneous redefinition of left and right. 2 There is a third discrete symmetry, time reversal (T ), and it is important to note that there is a theorem that states that CPT must be conserved in any locally Lorentz invariant quantum field theory [11]. Therefore, under rather reasonable assumptions, an observation of CP violation corresponds to an observation of T violation, and vice versa. Nonetheless, it remains of interest to establish T violation without assumptions regarding other symmetries [12,13].
The four parameters of the CKM matrix can be expressed in many different ways, but two popular choices are the Chau-Keung (PDG) parametrisation -(θ 12 , θ 13 , θ 23 , δ ) [14] -and the Wolfenstein parametrisation -(λ , A, ρ, η) [15]. In both cases a single parameter (δ or η) is responsible for all CP violation. This encapsulates the predictivity that makes the CKM theory such a remarkable success: it describes a vast range of phenomena at many different energy scales, from nuclear beta transitions to single top quark production, all by only four parameters (plus hadronic effects).
Let us digress a little into history. In 1964, CP violation was discovered in the kaon system [16], but it was not until 1973 that Kobayashi and Maskawa proposed that the effect originated from the existence of three quark families [3]. On a shorter time-scale, in 1967 Sakharov noted that CP violation was one of three conditions necessary for the evolution of a matter-dominated universe, from a symmetric initial state [17]: 1. baryon number violation, 2. C and CP violation, 3. thermal inequilibrium.
This observation evokes the prescient concluding words of Dirac's 1933 Nobel lecture, discussing his successful prediction of the existence of antimatter, in the form of the positron [18]: If we accept the view of complete symmetry between positive and negative electric charge so far as concerns the fundamental laws of Nature, we must regard it rather as an accident that the Earth (and presumably the whole solar system), contains a preponderance of negative electrons and positive protons. It is quite possible that for some of the stars it is the other way about, these stars being built up mainly of positrons and negative protons. In fact, there may be half the stars of each kind. The two kinds of stars would both show exactly the same spectra, and there would be no way of distinguishing them by present astronomical methods.
Dirac was not aware of the existence of CP violation, that breaks the complete symmetry of the laws of Nature. Moreover, modern astronomical methods do allow to search for antimatter dominated regions of the Universe, and none have been observed (though searches, for example by the PAMELA and AMS experiments, are ongoing). Therefore, CP violation appears to play a crucial role in the early Universe.
We can illustrate this with a simple exercise. Suppose we start with equal amounts of matter (X) and antimatter (X). The matter X decays to final state A (with baryon number N A ) with probability p and to final state B (baryon number N B ) with probability (1 − p). The antimatter,X, decays to final stateĀ (with baryon number −N A ) with probabilityp and final stateB (baryon number −N B ) with probability (1 −p). The resulting baryon asymmetry is So clearly ∆ N tot = 0 requires both p =p and N A = N B , i.e. both CP violation and baryon number violation.
It is natural to next ask whether the magnitude of the baryon asymmetry of the Universe could be caused by the CP violation in the CKM matrix. The baryon asymmetry can be quantified relative to the number of photons in the Universe, This can be compared to a dimensionless and parametrisation invariant measure of the amount of CP violation in the SM, J × P u × P d /M 12 , where • J = cos(θ 12 ) cos(θ 23 ) cos 2 (θ 13 ) sin(θ 12 ) sin(θ 23 ) sin(θ 13 ) sin(δ ) , • and M is the relevant scale, which can be taken to be the electroweak scale, O(100 GeV).
The parameter J is known as the Jarlskog parameter [19], and is expressed above in terms of the Chau-Keung parameters. Putting all the numbers in, we find a value of ∼ 10 −17 , much below the observed 10 −10 . This is the origin of the widely accepted statement that the SM CP violation is insufficient to explain the observed baryon asymmetry of the Universe. Note that this occurs primarily not because J is small, but rather because the electroweak mass scale is far above the mass of most of the quarks. Therefore, to explain the baryon asymmetry of the Universe, there must be additional sources of CP violation that occur at high energy scales. There is, however, no guarantee that these are connected to the CP violation that we know about.
The new sources may show up in the quark sector via discrepancies with CKM predictions (as will be discussed below), but could equally appear in the lepton sector as CP violation in neutrino oscillations. Or, for that matter, new sources could be flavour-conserving and be found in measurements of electric dipole moments, or could be connected to the Higgs sector, or the gauge sector, or to extra dimensions, or to other NP. In any case, precision measurements of flavour observables are generically sensitive to additions to the SM, and hence are well-motivated. In this context, it is worth noting the enticing possibility of "leptogenesis", where the baryon asymmetry is created via a lepton asymmetry (see, e.g., Ref. [20] for a review). In the case that neutrinos are Majorana particles -i.e. they are their own antiparticles -the right-handed neutrinos may be very massive, which provides an immediate connection with the needed high energy scale. Experimental investigation of this concept requires the determination of the lepton mixing (PMNS) [21,22] matrix, and proof whether or not neutrinos are Majorana particles. The recent determination of the neutrino mixing angle θ 13 [23,24] provides an important step forward; the next challenges are to establish CP violation in neutrino oscillations and to observe (or limit) neutrinoless double beta decay processes.

Rare processes
We have already digressed into history, and we should avoid doing so too much, but it is striking how often NP has shown up at the precision frontier before "direct" discoveries at the energy frontier. Examples include: the GIM mechanism being established before the discovery of charm; CP violation being discovered and the CKM theory developed before the discovery of the bottom and top quarks; the observation of weak neutral currents before the discovery of the Z boson. In particular, loop processes are highly sensitive to potential NP contributions, since SM contributions are suppressed or absent.
As a specific example of this we can consider the loop processes involved in oscillations of neutral flavoured mesons. (Rare decay processes will be discussed in more detail below.) There are four such pseudoscalar particles in nature (K 0 , D 0 , B 0 and B 0 s ) which can oscillate into their antiparticles via both short-distance (disper-sive) and long-distance (absorptive) processes, as illustrated in Fig. 1. Representing such a meson generically by M 0 , the evolution of the particle-antiparticle system is given by the time-dependent Schrödinger equation, where the effective 3 Hamiltonian H = M − i 2 Γ is written in terms of 2×2 Hermitian matrices M and Γ . Note that the CPT theorem requires that M 11 = M 22 and Γ 11 = Γ 22 , i.e. that particle and antiparticle have identical masses and lifetimes.
where p and q are complex coefficients that satisfy |p| 2 +|q| 2 = 1. Here the subscript labels L and H distinguish the eigenstates by their nature of being lighter or heavier; in some systems the labels S and L are instead used for short-lived and long-lived respectively (the choice depends on the values of the mass and width differences; the labels 1 and 2 are also sometimes used, usually to denote the CP eigenstates). CP is conserved (in mixing) if the physical states correspond to the CP eigenstates, i.e. if |q/p| = 1. Solving the Schrödinger equation gives with eigenvalues given by , corresponding to mass and width differences ∆ m = m H − m L and ∆Γ = Γ H − Γ L given by ∆ m∆Γ = 4 Re(M 12 Γ * 12 ) . Flavour Physics in the LHC Era 7 Note that with this notation, which is the same as that of Ref.
[25], ∆ m is positive by definition while ∆Γ can have either sign. 4 Rather than going into the details of the formalism (which can be found in, e.g., Ref. [27]) let us instead take a simplistic picture.
• The value of ∆ m depends on the rate of the mixing diagram of Fig. 1(left). This depends on CKM matrix elements, together with various other factors that are either known or (in the case of decay constants and bag parameters) can be calculated using lattice QCD. Moreover for the B mesons, these other factors can be made to cancel in the ∆ m d /∆ m s ratio, such that the measured value of this quantity gives a theoretically clean determination of |V td /V ts | 2 . • The value of ∆Γ , on the other hand, depends on the widths of decays of the meson and antimeson into common final states (such as CP-eigenstates). Therefore, ∆Γ is large for the K 0 system, where the two pion decay dominates, small for This simplistic picture is sufficient to explain qualitatively the experimental values of the mixing parameters given in Table 1. It should be noted that ∆Γ (B 0 s ) has become well-measured only very recently (as discussed below), and that the experimental sensitivity for the CP violation parameters in all of the D 0 , B 0 and B 0 s systems is still far from that of the SM prediction, making improved measurements very well motivated. Table 1 Qualitative expectations and measured values for the neutral meson mixing parameters. Experimental results are taken from Refs. [28][29][30]. The definition of a sl is given in footnote 6. Thus, neutral meson oscillations are rare processes described by parameters that can be both predicted in the SM and measured experimentally. All measurements to date are consistent with the SM predictions (though see below). These results 8 Tim Gershon can then be used to put limits on non-SM contributions. This can be done within particular models, but the model-independent approach, described in, e.g., Ref. [31] is illustrative. The NP contribution is expressed as a perturbation to the SM Lagrangian, where the dimension d of higher than 4 has an associated scale Λ and couplings c i . 5 Given the observables in a given neutral meson system, NP contributions described effectively as four-quark operators (d = 6) can be constrained, either by putting bounds on Λ for a fixed value of c i (typically 1), or by putting bounds on c i for a fixed value of Λ (typically 1 TeV). In the former case bounds of O(100 TeV) are obtained; in the latter case the bounds can be O(10 −9 ) or below [31], with the strongest (weakest) bounds being in the K 0 (B 0 s ) sectors. A similar analysis, but with more up-to-date inputs has been performed in Ref. [32], with results illustrated in Fig. 2. The mixing amplitude, normalised to its SM value, is denoted by ∆ , and experimental constraints give (Re∆ , Im∆ ) consistent with (1, 0) (i.e. with the SM) for both B 0 and B 0 s systems. This is a very puzzling situation. Limits on the NP scale give values of at least 100 TeV for generic couplings. But, as discussed elsewhere, we expect NP to appear at the TeV scale to solve the hierarchy problem (and to provide a dark matter candidate, etc.) If NP is indeed at this scale, NP flavour-changing couplings must be small. But why? This is the so-called "new physics flavour problem".
A theoretically attractive solution to this problem, known as minimal flavour violation (MFV) [33], exploits the fact that the SM flavour-changing couplings are also small. Therefore, if there is a perfect alignment of the flavour violation in a NP model with that in the SM, the tension is reduced. The MFV paradigm is highly predictive, stating that there are no new sources of CP violation and also that the correlations between certain observables share their SM pattern (the ratio of branching fractions of B 0 → µ + µ − and B 0 s → µ + µ − being a good example). Therefore, once physics beyond the SM is discovered, it will be an important goal to establish whether or not it is minimally flavour violating. This further underlines that the flavour structure carries information about physics at very high scales.
Nonetheless, it must be reiterated that there are several important observables that are not yet well measured, and that could rule out MFV. For example, the bounds on NP scales obtained above (from Ref. [31]) do not include results on CP violation in mixing in the B 0 and B 0 s sectors. In fact, the D0 collaboration has reported a measurement of an anomalous effect [34] of the inclusive same-sign dimuon asymmetry, which is 3.9σ away from the SM prediction (of very close to zero [35]). This measurement is sensitive to an approximately equal combination of the parameters of CP violation in B 0 and B 0 s mixing, a d sl and a s sl , 6 however some sensitivity to the source of the asymmetry can be obtained by applying additional constraints on the impact parameter to obtain a sample enriched in either oscillated B 0 or B 0 s candidates. In addition, a d sl and a s sl can be measured individually. The latest world average, shown in Fig. 3, gives a d sl = −0.0003 ± 0.0021, a s sl = −0.0109 ± 0.0040 [30]. Improved measurements are needed to resolve the situation. The green ellipse comes from the D0 inclusive same-sign dimuon analysis [34]; the blue shaded bands give the world average constraints on a d sl and a s sl individually; the red ellipse is the world average including all constraints [30].

The CKM matrix and the Unitarity Triangle
Much of the experimental programme in heavy quark flavour physics is devoted to measurements of the parameters of the CKM matrix. As discussed above, the CKM matrix can be written in terms of the Wolfenstein parameters, which exploit the observed hierarchy in the mixing angles: where the expansion parameter λ is the sine of the Cabibbo angle (λ = sin θ C ≈ V us ). It should be noted that although the hierarchy is highly suggestive, there is no underlying reason known for this pattern; moreover, the pattern in the lepton sector is completely different. Note also that at O(λ 3 ) in the Wolfenstein parametrisation, the complex phase in the CKM matrix enters only in the V ub and V td (top right and bottom left) elements, but this is purely a matter of convention -only relative phases are observable. The unitarity of the CKM matrix, , puts a number of constraints on the magnitudes and relative phases of the elements. Among these relations, one which has been precisely tested is where the measurements of |V ud | 2 from, e.g., super-allowed β decays and |V us | 2 from leptonic and semileptonic kaon decays are indeed consistent with the prediction to within one part in 10 3 [36]. 7 The unitarity condition also results in six constraints, which correspond to three complex numbers summing to zero, and hence can be represented as triangles in the complex plane. The triangles have very different shapes, but all of them have the same area, which is given by half of the Jarlskog parameter [19]. The particular triangle relation given in Eq. (10) is particularly useful to visualise the constraints from various different measurements, as shown in the iconic images from the CKMfitter [37] and UTfit [38] collaborations, reproduced in Fig. 4. Conventionally, this "Unitary Triangle" (UT) is rescaled by V cd V * cb so that by definition the position of the apex is where (ρ, η) [39] are related to the Wolfenstein parameters by Two popular naming conventions for the UT angles exist in the literature: (13) The (α, β , γ) set is used in these lectures. The lengths of the sides R u and R t of the UT are given by A major achievement of the past decade or so has been to significantly improve the precision of the parameters of the UT. In particular, the primary purpose of the so-called "B factory" experiments, BaBar and Belle, was the determination of sin 2β using B 0 → J/ψ K 0 S (and related modes). This was carried out using completely new experimental techniques to probe CP violation in a very different way to previous experiments in the kaon system. In particular, if we denote the amplitude for a B 0 meson to decay to a particular final state f as A f , and that for the charge conjugate process asĀf , then using the parameters p and q from Eq. (3), we define the parameter λ f = q pĀf A f and the following categories of CP violation in hadronic systems: 8 1. CP violation in mixing (|q/p| = 1), 3. CP violation in interference between mixing and decay (Im λ f = 0).
Additionally, in the literature the concepts of indirect and direct CP violation are often encountered: the former is where the effect is consistent with originating from a single phase in the mixing amplitude, while the latter cannot be accounted for in such a way. Following this categorisation, CP violation in mixing is indirect, CP violation in decay (the only category available to baryons or charged mesons) is direct, while CP violation in interference can be indirect so long as only one measurement is considered -but if two such measurements give different values, this also establishes direct CP violation.

Determination of sin(2β )
The determination of sin(2β ) from B 0 → J/ψ K 0 S [40,41], exploits the fact that some measurements of CP violation in interference between mixing and decay can be cleanly interpreted theoretically, since hadronic factors do not contribute. The full derivation of the decay-time-dependent decay rate of B 0 mesons to a CP eigenstate f is a worthwhile exercise for the reader, and can be found in, e.g., Refs. [42,43]. The result, for mesons that are known to be either B 0 or B 0 at time t = 0, is where In these expressions ∆Γ has been assumed to be negligible, as appropriate for the B 0 system. Assuming further |q/p| = 1, then for decays dominated by a single amplitude, C f = 0 and S f = sin(arg(λ f )), and so for B 0 → J/ψ K 0 S , S = sin(2β ), to a very good approximation.
Flavour Physics in the LHC Era 13 The experimental challenge for the measurement of sin(2β ) then lies in the ability to measure the coefficient of the sinusoidal oscillation of the decay-timeasymmetry. Until recently, the most copious sources of cleanly reconstructed B mesons came from accelerators colliding electrons with positrons at the ϒ (4S) resonance (a bb bound state just above the threshold for decay into pairs of B mesons). For symmetric colliders, the B mesons are produced at rest, and therefore lifetime measurements are not possible. A boost is necessary, which can be advantageously achieved by making the e + e − collisions asymmetric. 9 One strong feature of this approach is that the quantum correlations of the B mesons produced in ϒ (4S) decay are retained, so that the decay of one into a final state that tags its flavour (B 0 or B 0 ) can be used to set the clock to t = 0 and specify the flavour of the other at that time.
The concept of the asymmetric B factory was such a good one that two were built: PEP-II at SLAC, colliding 9.0 GeV e − with 3.1 GeV e + , and KEKB at KEK (8.0 GeV e − on 3.5 GeV e + ). These have achieved world record luminosities, with peak instantaneous luminosities above 10 34 cm −2 s −1 , and a combined data sample of over 1 ab −1 , corresponding to over 10 9 BB pairs. The detectors (BaBar [44,45] and Belle [46] respectively) share many common features, such as silicon vertex detectors, gas based drift chambers, electromagnetic calorimeters based on Tl-doped CsI crystals, and 1.5 T solenoidal magnetic fields. The main difference is in the technology used to separate kaons from pions: a system based on the detection of internally reflected Cherenkov light for BaBar, and a combination of aerogel Cherenkov counters and a time-of-flight system for Belle.
Through the measurement of sin 2β , BaBar [47] and Belle [48] were able to make the first observations of CP violation outside the kaon sector, thus validating the Kobayashi-Maskawa mechanism. The latest (and, excluding upgrades, most likely final) results from BaBar [49] and Belle [50] shown in Fig. 5 give a clear visual confirmation of the large CP violation effect. The world average value, using determinations based on b → ccs transitions, is [30] sin 2β = 0.682 ± 0.019 which gives β = (21.5 +0.8 −0.7 ) • .

Determination of α
Additional measurements are needed to over-constrain the UT and thereby test the Standard Model. The angle α can, in principle, be determined in a similar way as β , but using a decay mediated by the b → uūd tree-diagram which carries the relative weak phase γ (since π − (β + γ) = α by definition). However, in any such decay a contribution from the b → d loop ("penguin") amplitude, which carries a different weak phase, is also possible. This complicates the interpretation of the observables, since S = sin(2α); on the other hand direct CP violation becomes observable, if the relative strong phase is non-zero. Constraints on α can still be obtained using a is dominated by constraints from the B 0 → ρ + ρ − decay [52,53], which is consistent with having negligible penguin contribution.

The sides of the Unitarity Triangle
The lengths of the sides of the UT have also been constrained by various observables. The value of R t depends on |V td |, and can be determined from b → d transitions such as the rate of B 0 oscillations, i.e. ∆ m d , or the branching fraction B → ργ. In both cases, theoretical uncertainties are reduced if the measurement is performed relative to that for the corresponding b → s transition. The most precise constraint to date comes from the ratio of ∆ m d [30,54,55] and ∆ m s [30,56,57] and gives V td V ts = 0.211 ± 0.001 ± 0.005, where the first uncertainty is experimental and the second theoretical (originating from lattice QCD calculations).
The value of R u depends on |V ub | and can be determined from b → u tree-level transitions. Semileptonic decays allow relatively clean theoretical interpretation, 10 but still require QCD calculations to go from the parton level transition to the observed (semi-hadronic) final state (for a recent review, see Ref. [58]). Two approaches have been pursued: exclusive decays, such as B 0 → π − e + ν, and inclusive decays, B → X u e + ν. The theory of inclusive decays is based on the operator product expansion (discussed in Sec. 5.7) and would be extremely clean, were it not for the fact that experimentally cuts are needed to remove the more prevalent b → c transition. Exclusive decays tend to have less background from b → c processes. The differential branching fractions can be translated in constraints on |V ub | using knowledge of form-factors at the kinematic limit obtained from lattice QCD calculations, together with phenomenological models that extrapolate over the whole phase space. The most precise results use B → π + ν decays ( = e, µ) [59][60][61], and give an "exclusive" determination of |V ub | that is, however, in tension with the "inclusive" value [30]: where the first uncertainties are experimental and the second theoretical. Since the origin of the discrepancy, which is also seen in determinations of |V cb | from b → c ν transitions, is not understood, the uncertainty is scaled to give The results on β , α, R t and R u are the most constraining inputs to the CKM fits shown in Fig. 4 [37,38]. While the results are all consistent with the Standard Model prediction of a single source of CP violation, there are some tensions that deserve further investigation. Moreover, there are still certain important observables where large NP contributions are possible, as will be discussed in more detail below.

Flavour physics at hadron colliders
Results from the B factory experiments provided enormous progress in the understanding of heavy flavour physics (only a very brief selection has been discussed above). Nonetheless, many results remain statistically limited, and the B 0 s sector is relatively unexplored. To progress further, it is necessary to have a copious source of production of all flavours of b hadron. As shown in Table 2, high energy hadron colliders satisfy these criteria, but present significant experimental challenges: to be able to identify the decays of interest from the high multiplicity environment, and to reject the even more copious rate of minimum bias events. 11 The LHCb detector [63], shown in Fig. 6, has been designed to meet these challenges. It is in essence a forward spectrometer (covering the acceptance region that optimises its flavour physics capability), with a dipole magnet, a precision silicon vertex detector and strong particle identification capability. Tracks can be identified as different hadron species using information from ring-imaging Cherenkov detectors, while calorimeters and muon detectors enable charged leptons to be distinguished and also provide trigger signals. The trigger system [64] uses these hard- 16 Tim Gershon BB pair alone Many additional particles Production vertex Not reconstructed Reconstructed from many tracks B 0 -B 0 pair production Coherent (from ϒ (4S) decay) Incoherent Flavour tagging power εD 2 ∼ 30% εD 2 ∼ 5% Fig. 6 The LHCb detector [63].
ware level signals to reduce the rate from the maximum LHC bunch-crossing rate of 40 MHz to the 1 MHz rate at which the detector can be read out. A software trigger then searches for inclusive signatures of b-hadron decays such as high-p T signals and displaced vertices, and also performs reconstruction of several exclusive b and c decay channels, in order to further reduce the rate to a level that can be written to offline data storage (3 kHz in 2011, 5 kHz in 2012). During the LHC run, the detector operated with data taking efficiency above 90%, with instantaneous luminosity around 3 (4) × 10 32 cm −2 s −1 recording data samples of 1 (2) fb −1 at √ s = 7 (8) TeV in 2011 (2012). 12 The luminosity is less than that delivered to ATLAS and CMS, since the experimental design requires low pileup, i.e. a low number of pp collisions per bunch-crossing. However, this allows the luminosity to be "levelled" and remain at a constant value throughout the LHC fill, providing very stable data taking-conditions. 13 In addition, the polarity of the dipole magnet is reversed regularly, which enables cancellation of detector asymmetries in various measurements.
In addition to LHCb, it must be noted that the "general purpose detectors" AT-LAS and CMS at the LHC, and CDF and D0 at the Tevatron, have capability to study flavour physics. For most of these experiments, their programme is, however, restricted to decay modes triggered by high p T muons, but CDF benefited from a two-track trigger [65] that enabled a broader range of measurements to be performed.

Heavy flavour production and spectroscopy
The capabilities of the different experiments can be demonstrated from the measurements of production cross-sections that have been performed by each. Most have studied J/ψ production (e.g. Refs. [66][67][68][69][70][71]) as well as b hadron production using decay modes involving muons or J/ψ mesons [72][73][74][75][76][77]. However, only CDF and LHCb have been able to study prompt charm production [78,79]. 14 The crosssections measured confirm the theoretical predictions, and enable the values for integrated luminosity to be translated into more easily comprehensible terms. For example, with 1 fb −1 recorded at √ s = 7 TeV, and the measured bb production crosssection [77,80], it is easily shown that over 10 11 bb quark pairs have been produced in the LHCb acceptance. This compares to the combined BaBar and Belle data sample of ∼ 10 9 BB meson pairs. Consequently, for any channel where the efficiency, including effects from reconstruction, trigger and offline selection, is not too small, LHCb has the world's largest data sample. This further emphasises the need for an excellent trigger in order to perform flavour physics at hadron colliders.
Production measurements such as those mentioned above test QCD models, and are important (and highly-cited) results. However, since they are not within the remit of flavour-changing interactions of the charm and beauty quarks, they will not be discussed further here. Nonetheless, a brief digression into studies of another aspect of QCD, that of spectroscopy, will be worthwhile. This covers the study of properties of hadronic states such as lifetimes, masses, decay channels and quantum numbers, and also the discoveries of new states. Indeed, some of the most highly-cited papers from recent flavour physics experiments relate to such topics, including the discovery of the X(3872) particle by Belle [81] and of the D sJ states by BaBar [82] and CLEO [83]. The first new particles discovered at the LHC, prior to the Higgs boson, were hadronic states [84][85][86]. More recently, significant progress has been made in understanding the nature of the X(3872) [87]. New results are eagerly anticipated in several related areas, for example to clarify the situation regarding the existence of charged charmonium-like states, claimed by Belle [88][89][90] but not confirmed by BaBar [91,92], which would be smoking gun signatures for an exotic hadronic nature. 15 Recent claims of charged bottomonium-like states by Belle [95,96] seem to strengthen the case that such exotics can exist in nature, but one should note that history teaches us that not all claimed states turn out to be real [97].
The topic of spectroscopy also provides a useful illustration of the importance of triggering for flavour physics experiments at hadron colliders. In 2008, the BaBar experiment discovered the η b meson using the process e + e − →ϒ (3S) → η b γ, where only the photon is reconstructed and the signal is inferred from a peak in the photon energy spectrum [98]. The η b meson is the pseudoscalar bb ground state. It is the lightest bottomonium state, so why did it take more than 30 years after the discovery of the ϒ (1S) meson [99] (the lightest vector bb state) to see it in experiments? In particular, since hadron collisions produce particles with all possible quantum numbers, why was it not observed at, e.g., the Tevatron? The answer lies in the fact that the vector state decays to dimuons, which have a distinctive trigger signature. The dominant decay channels of the η b are expected to be multibody hadronic final states, which make its observation in a hadronic environment extremely challenging.

Direct CP violation
As mentioned above, the condition for direct CP violation is Āf /A f = 1. In order for this to be realised we need the amplitude to be composed of at least two parts with different weak and strong phases. This is often realised by tree (T ) and penguin (P) amplitudes (example diagrams are shown in Fig. 7), so that where the strong (weak) phases δ T,P (φ T,P ) keep the same (change) sign under the CP transformation by definition. The CP asymmetry is defined from the rate difference between the particle involving the quark (D orB) and that containing the antiquark (D or B). Using the definition for B decays, we trivially find Therefore, for large direct CP violation effects to occur, we need |P/T |, sin(δ T − δ P ) and sin(φ T − φ P ) to all be O(1). Charmless B decays, i.e. decays of B mesons to final states that do not contain charm quarks, provide good possibilities for the observation of direct CP violation, since many decays have both tree and penguin contributions with similar levels of CKM suppression. These are of interest to search for NP, since the penguin loop diagrams are sensitive to potential contributions from new particles. An excellent example is B 0 → K + π − , which provided the first observation of direct CP violation outside the kaon sector, and has a world average value of A CP (B 0 → K + π − ) = −0.086 ± 0.007 [30,[100][101][102][103]. Curiously, the CP violation effect observed in B + → K + π 0 decays is rather different: A CP (B + → K + π 0 ) = 0.040 ± 0.021 [30,100,101], although naïvely changing the spectator quark in Fig. 7 suggests that similar values should be expected. This is referred to as the "Kπ puzzle", and could in principle be a hint for NP, though the more mundane explanation of larger than expected QCD corrections cannot be ruled out at present. Several methods are available to test the QCD explanations, which motivate improved measurements of other Kπ modes (in particular, of A CP (B 0 → K 0 S π 0 )), of similar decay modes with three-body final states (Kρ, K * π), and of charmless two-body B 0 s decays. On this last topic, following pioneering work by CDF [103,104], LHCb has recently reported both the first decay time-dependent analysis of B 0 s → K + K − [105] and the first observation of CP violation in B 0 s → K − π + decays [106], which demonstrate good prospects for progress in the coming years.
With regard to three-body decays, it is worth noting that despite hundreds of measurents by the B factories, the significance of the world average in any other charmless B + or B 0 decay mode does not exceed 5 σ , though channels such as B + → ηK + and B + → ρ 0 K + approach this level. However, very recently, LHCb has demonstrated that large CP violation effects occur in specific regions of the phase space of three-body charmless decays such as B + → K + π + π − [107-109]. Further study is necessary to quantify the effect and identify its origin.

The UT angle γ from B → DK decays
The angle γ of the CKM Unitarity Triangle is unique in that it is the only CPviolating parameter that can be measured using only tree-level decays. This makes it a benchmark Standard Model reference point. Improving the precision with which γ is known is one of the primary goals of contemporary flavour physics, and this will only become more important after NP is discovered, since it will be essential to disentangle SM and NP contributions to CP-violating observables.
The phase γ can be determined exploiting the fact that in decays of the type B → DK, the b → cūs and b → ucs amplitudes can interfere if the neutral charmed meson is reconstructed in a final state that is accessible to both D 0 and D 0 decays. There are many possible such final states, with various experimental advantages and disadvantages. These include CP eigenstates, doubly-or singly-Cabibbo-suppressed decays and multibody final states. Moreover, decays of different b hadrons can all be used to provide constraints on γ. Two particularly interesting approaches are to study decay time-dependent asymmetries of B 0 s → D ∓ s K ± decays [110] and to study the Dalitz plot (i.e. phase-space) dependent asymmetries in B 0 → DK + π − decays [111,112]. First results from LHCb show promising potential for these decays [113,114]. All such measurements will help to improve the overall precision in a combined fit. Fig. 8 Illustration of the concept behind the determination of γ using B ± → DK ± decays. For B − decays the amplitudes add with relative phase δ − γ, while for B + the relative phase is δ + γ. Here the simplest case with D decays to CP eigenstates (such as K + K − ) is shown, but the method can be extended to any final state accessible to both D 0 and D 0 decays.
The basic concept behind the method is illustrated in Fig. 8 for B − → D CP K − decays. It must be emphasised that due to the absence of loop contributions to the decay it is extremely clean theoretically [115]. This, and the abundance of different final states accessible, means that all parameters can be determined from data. The relevant parameters are the weak phase γ, an associated strong (CP conserving) phase difference between the b → cūs and b → ucs decay amplitudes, labelled δ B , and the ratio of their magnitudes, r B . The small value of r B (B − → DK − ) ∼ 10 % means that large event samples are necessary to obtain good constraints on γ, and only recently has the first 5σ observation of CP violation in B → DK decays been achieved [116]. Larger values of r B are expected in B 0 → DK * 0 and B 0 s → D ∓ s K ± decays, but until now the samples available in these channels have not been been sufficient to give meaningful constraints on γ. . Significant progress in this area is anticipated from LHCb in the coming years. 16 5.3 Mixing and CP violation in the B 0 s system A complete analysis of the time-dependent decay rates of neutral B mesons must also take into account the lifetime difference between the eigenstates of the effective Hamiltonian, denoted by ∆Γ . This is particularly important in the B 0 s system, since the value of ∆Γ s is non-negligible. Neglecting CP violation in mixing, the relevant replacements for Eq. (15) are [122] where N is a normalisation factor and Note that, by definition, Also A ∆Γ f is a CP-conserving parameter, unlike S f and C f (since it appears with the same sign in equations for both B 0 s and B 0 s states). Nonetheless, it provides sensitivity to arg(λ f ), which means that interesting results can be obtained from untagged timedependent analyses (a.k.a. effective lifetime measurements [123]).
The formalism of Eq. (20) is usually invoked for the determination of the CP violation phase in B 0 s oscillations, φ s = −2β s , using B 0 s → J/ψ φ decays. However, in that case things are complicated even further by the fact that the final state, containing two vector particles, is an admixture of CP-even and CP-odd which must be disentangled by angular analysis. 17 Moreover, there is a potential contribution from S-wave K + K − pairs within the φ mass window used in the analysis. However, all of these features can be turned to the benefit of the analysis, providing better sensitivity and allowing to resolve an ambiguity in the results [125]. A compilation of the latest results is shown in Fig. 9. 18 Although great progress has been made over the last few years, the experimental precision does not yet challenge the theoretical uncertainty, and so further updates are of great interest.  Fig. 9 Compilation of the latest results [30] in the φ s vs. ∆Γ s plane from LHCb [128], CDF [129], D0 [130] and ATLAS [131].

Mixing-induced CP violation in hadronic b → s penguin decay modes
As discussed in Sec. 5.1, decay modes mediated by penguin diagrams are potentially sensitive to NP effects, although it is a considerable challenge to disentangle QCD effects. One useful approach is to study mixing-induced CP violation effects in channels that are dominated by the penguin transition, so that little or no tree (or other) contribution is expected. Such channels include 19 For the B 0 decays, the formalism is the same as given in Eq. (15), and the parameters are expected in the SM to be given, to good approximations, by C f ≈ 0, S f ≈ sin(2β ) (up to a sign, given by the CP eigenvalue of the final state). These channels have been studied extensively by BaBar and Belle: early results provided hints for discrepancies with the SM predictions, but the significance of the deviation diminshed with improved results [132][133][134][135][136]. For the B 0 s decays, the formalism is as given in Eq. (20) (though with modifications due to the vector-vector final states), and the SM expectation is that CP violation effects vanish, to a good approximation, since the very small phase in the b → s decay cancels that in the B 0 s -B 0 s oscillations. First results have been reported by LHCb [137][138][139], and will reach a very interesting level of sensitivity as more data is accumulated.

Charm mixing and CP violation
In the charm system the mixing parameters x = ∆ m/Γ and y = ∆Γ /(2Γ ) are both small, x, y 1. Therefore, a Taylor expansion can be performed on the generic expression of Eq. (20) to give Hence an untagged analysis of D 0 → K + K − can measure A ∆Γ f y (also known as y CP ), while a tagged analysis can additionally probe S f x. Since the mixing parameters are small, the focus until now has been to establish definitively oscillation effects, but in the coming years the main objective will be to observe or limit CP violation in the charm system, which is expected to be very small in the SM. Note that in case the source of D 0 mesons is either from D * + decays or semileptonic bhadron decays, the flavour tagging is very effectively achieved from the charge of the associated pion or lepton, respectively. Many other final states can be used to gain additional sensitivity to charm mixing and CP violation parameters, a recent example being the observation of charm mixing at LHCb using D 0 → K + π − decays [140]. The result of this analysis, and the world average constraints on the x and y parameters in the D 0 system, 20 are shown in Fig. 10.
Direct CP violation in the charm system can also be used to test the SM. One interesting recent result has been the measurement of ∆ A CP , which is the difference between the direct CP violation parameters of D 0 → K + K − and D 0 → π + π − decays. By measuring the difference, a cancellation of production and detection asymmetries can be exploited, while the physical CP asymmetry may be enhanced. 21 This method was first used by LHCb [141] and then by CDF [142] and Belle [143], all indicating a larger than expected effect. This prompted a great deal of theoretical activity, summarised in Ref. [144], with the conclusion that a SM origin of the CP violation, although unlikely, was not ruled out. Many further studies were proposed to test both SM and NP hypotheses, and these remain of great interest and will Tim Gershon Fig. 10 (Left) Decay-time evolution of the ratio, R, of D 0 → K + π − to D 0 → K − π + yields (points) with the projection of the mixing allowed (solid line) and no-mixing (dashed line) fits overlaid, from Ref. [140]; (Right) World average constraints on the x and y parameters in the D 0 system [30]. be pursued. However, the most recent results by LHCb [145,146] suggest that the central value is smaller than previous thought, and therefore the SM explanation becomes harder to rule out.

Photon polarisation in radiative B decays
The b → sγ transition is an archetypal flavour-changing neutral-current (FCNC) transition, and has been considered a sensitive probe for NP since the first measurements of its rate [147,148]. The latest results for the inclusive branching fraction [30] are consistent with the SM prediction [149] B (B → X s γ) exp E γ >1.7 GeV = (3.43 ± 0.21 ± 0.07) × 10 −4 , However, additional observables, such as CP and isospin asymmetries provide complementary sensitivity and still have experimental uncertainties much larger than those of the theoretical predictions of their values in the SM. One particularly interesting observable is the polarisation of the emitted photon in b → sγ decays, since the V − A structure of the SM weak interaction results in a high degree of polarisation, that is not necessarily reproduced in extended models. Until now, the most promising approach to probe the polarisation has been from time-dependent asymmetry measurements of B 0 → K 0 S π 0 γ [150,151] but the most recent measurements [152,153] still have large uncertainties. LHCb can use several different methods to study photon polarisation in b → sγ transitions, such as measuring the effective lifetime in B 0 s → φ γ decays [154]. Although all such measurements are highly challenging, the observed yields in B 0 s → φ γ [155] and other related channels such as B 0 → K * 0 e + e − [156] suggest there are good prospects for significant progress in the coming years.

Angular observables in
The b → sl + l − FCNC transitions offer similar, but complementary, sensitivity to NP as b → sγ, but are experimentally more convenient to study, in particular when the lepton pair is muonic, i.e. l + l − = µ + µ − . The multi-body final state provides access to a wide range of kinematic observables, several of which have clean theoretical predictions (especially at low values of the dilepton invariant mass squared, q 2 ). This makes these decays a superb laboratory for NP tests. The theoretical framework for these (and other) processes is the operator product expansion. This is an effective theory, applicable for b physics, which describes the weak interactions of the SM by integrating out the heavier (W , Z, t) fields. As such it can be considered a modern version of the Fermi theory of beta decay. Conceptually, it can be expressed as where O n represent the local interaction terms, and C n are coupling constants that are referred to as Wilson coefficients. 22 The Wilson coefficients encode information on the weak scale, and are calculable and known in the SM (at least to leading order). Moreover, they are affected by NP, so comparing the measured values with their expectations provides tests of the SM. A more detailed description of the operator product expansion can be found in, e.g. Ref. [157]. For the purposes of discussing b → sl + l − decays, the Wilson coefficients of interest are C 7 (which also affects b → sγ), C 9 and C 10 . The differential decay distribution, for the inclusive process, can be written [158] where θ l is the angle between the momentum vectors of the positively charged lepton and the opposite of the decaying b hadron in the dilepton rest frame. 23 The coefficients are given by Note that the term involving H A depends linearly on cos θ l and hence gives rise to a q 2 -dependent forward backward asymmetry, A FB . The shape of A FB , in particular the value of q 2 at which it crosses zero, can be predicted with low uncertainty in the SM. The expressions for exclusive processes, such as B 0 → K * 0 µ + µ − , are conceptually similar to those of Eqs. (27) and (28), but are more complicated as they also involve hadronic form factors. On the other hand, exclusive channnels also provide additional observables that can be studied (such as the longitudinal polarisation of the K * 0 meson, F L ), some of which can be precisely predicted in the SM, and are sensitive to NP contributions. The decay rates and angular distributions of B 0 → K * 0 µ + µ − decays have been studied by many experiments, with the most precise results to date, from LHCb [159], shown in Fig. 11. This analysis provides the first measurement of the A FB zero crossing point, q 2 0 = 4.9 ± 0.9 GeV 2 /c 4 , consistent with the SM prediction. Significant progress, including improved measurements of other NP-sensitive angular observables, can be expected in the coming years.

The very rare decay
The "killer app." for flavour physics as a tool to probe for (and potentially discover) NP is the very rare decay B 0 s → µ + µ − . The branching fraction is highly suppressed in the SM due to a combination of three factors, none of which are necessarily re-Flavour Physics in the LHC Era 27 produced in extended models: (i) the absence of tree-level FCNC transitions; (ii) the V − A structure of the weak interaction that results in helicity suppression of purely leptonic decays of flavoured pseudoscalar mesons; (iii) the hierarchy of the CKM matrix elements. In particular, in the minimally supersymmetric extension of the SM, the presence of a pseudoscalar Higgs particle alleviates the helicity suppression and enhances the branching fraction by a factor proportional to tan 6 β /M 4 A 0 , where tan β is the ratio of Higgs' vacuum expectation values, and M A 0 is the pseudoscalar Higgs mass. Therefore, in the region of phase-space where tan β is not too small, and M A 0 is not too large, the decay rate can be significantly enhanced above its SM expectation [160], 24 Due to the very clean signature of this decay, it has been studied by essentially all high-energy hadron collider experiments. The crucial components to obtain good sensitivity are high luminosity, a large B production cross-section within the acceptance, and good vertex and mass resolution to reject the background. Although AT-LAS [162] and CMS [163] have collected more luminosity, at present the strengths of the LHCb detector have allowed it to obtain the most precise results for this mode. Following a series of increasingly restrictive upper limits [164][165][166], LHCb has recently obtained the first evidence, with 3.5σ significance, for the decay [167], as shown in Fig. 12. The branching fraction is measured to be in agreement with the SM prediction. Further updates of this measurement are keenly anticipated, and are likely to appear at regular intervals throughout the lifetime of the LHC. It is worth noting that even in case the B 0 s → µ + µ − branching fraction remains consistent with the SM, the decay provides an additional handle on NP through its effective lifetime [168]. Moreover, it will be important to study also the even more suppressed B 0 → µ + µ − decay, since the ratio of the B 0 and B 0 s branching fractions is a benchmark test of MFV.

Future flavour physics experiments
As stressed in the previous sections, the first results from the LHC have already provided dramatic advances in flavour physics, and significant further progress is anticipated in the coming years. However, the instantaneous luminosity of LHCb is limited due to the fact that its subdetectors are read out at 1 MHz. As shown in Fig. 13 (left), increasing the luminosity requires tightening of the hardware trigger thresholds in order not to exceed this limit. This then results in lower efficiencies, especially for decay channels triggered by the calorimeter (i.e., channels without muons in the final state), so that there is no net gain in yield. Therefore, after several years of operation at the optimal instantaneous luminosity at √ s = 13 or 14 TeV, 25 the time required to double the accumulated statistics becomes excessively long. As should be clear from the discussions above, it remains of great importance to pursue a wide range of flavour physics measurements and improve their precision to the level of the theoretical uncertainty, and therefore it is of clear interest to get past the 1 MHz readout limitation. The concept of the LHCb upgrade [169,170] is to read out the full detector at 40 MHz (which corresponds to the maximum bunch crossing rate, with 25 ns spacing) and implement the trigger fully in software. This will allow the experiment to run at higher luminosities, up to 1 or 2 × 10 33 cm −2 s −1 , and will also significantly improve the efficiency for modes currently triggered by calorimeter signals at the hardware level. The accumulated samples in most key modes will increase by around two orders of magnitude compared to what was collected in 2011. Moreover, with a flexible trigger scheme, the capability to search for other signatures of NP will be enhanced, so that the upgraded experiment can be considered a general purpose detector in the forward region. The LHCb upgrade is planned to occur during the second long shutdown of the LHC, in 2018. Since its target luminosity is still below that which can be delivered by the LHC, it does not depend (though it is consistent with) the HL-LHC machine upgrade.
There are several other flavour physics experiments that will be coming online on a similar same timescale. The KEKB accelerator and Belle experiment are being upgraded [171], in order to allow luminosities almost two orders of magnitude larger than has previously been achieved. Compared to the LHCb upgrade, the e + e − environment is advantageous for decay modes with missing energy and for inclusive measurements. Some of the key channels for Belle2 are lepton flavour violating decays of τ leptons, mixing-induced CP asymmetries in decays such as B 0 → φ K 0 S and B 0 → η K 0 S , and the leptonic decay B + → τ + ν (which can be considered a counterpart of B 0 s → µ + µ − , and is sensitive to the exchange of charged Higgs particles) [1,172].
In addition, the NA62 [173] and K0T0 [174] experiments will search for the K + → π + νν and K 0 L → π 0 νν decays, respectively. Long considered the "holy grail" of kaon physics these decays are highly suppressed in the SM and have clean theoretical predictions. The new generation of experiments should be able to observe these channels for the first time, if they occur at around the SM rate.

Conclusion
Flavour physics continues to present many mysteries, and these demand continued experimental and theoretical investigation. Heavy flavour physics is complementary to other sectors of the global particle physics programme such as the high-p T experiments at the LHC, and neutrino oscillation and low energy precision experiments. The prospects are good for significant progress in the coming few years and, with upgraded experiments planned to come online in the second half of this decade, beyond.