Polarisations of the $Z$ and $W$ bosons in the processes $pp \to ZH$ and $pp \to W^{\pm}_{}H$

The $Z$ boson in the process $pp \to ZH$ and the $W^+_{}$ and $W^-_{}$ in the process $pp \to W^{\pm}_{}H$ can be in polarised states. The polarisation density matrix of the $Z$ ($W$) boson contains the complete information about a state of polarisation of the $Z$ ($W$) boson, and $HZZ$, $HZ\gamma$ and $HWW$ interactions may be studied in detail from a careful analysis of these matrices. In this paper, a systematic approach to analyse these polarisation density matrices is presented. With the aim of making maximum use of the polarisation information, all of the elements of the polarisation density matrices are related with observables, which are measurable at the environment of $pp$ collisions. Consequences of non-standard $HZZ$, $HZ\gamma$ and $HWW$ interactions for these observables are discussed.

An arbitrary polarised state of a massive particle with spin s is entirely described by a (2s + 1) × (2s + 1) polarisation density matrix; see e.g. refs. [57][58][59][60]. When a particle with spin 1 is in an either completely or partially polarised state, it has 9 independent angular distributions of its decay products at its rest frame and the coefficients of these distributions are written in terms of all the elements of the polarisation density matrix; see e.g. refs. [59,60]. An experimental determination of the complete decay angular distributions of a particle's decay products is, therefore, identical to an attempt to determine the polarisation density matrix, from which detailed information in reactions which produce that particle can be extracted; see e.g. ref. [61] for a theoretical study and e.g. refs. [62][63][64] for experimental studies 1 .
As the Z boson in the process e + e − → ZH can be in a polarised state [6,[72][73][74], so the Z boson in the process pp → ZH can be. The W + and W − in the process pp → W ± H can be also in polarised states. It would be, therefore, possible to study HZZ, HZγ and HW W interactions in detail from a careful analysis of these states of polarisation. As we mentioned in the previous paragraph, the polarisation density matrix of the Z (W ) boson contains the complete information about a state of polarisation of the Z (W ) boson. Hence all of the elements of the polarisation density matrix should be made use of in such a careful analysis. In this paper, we present a systematic approach to analyse the polarisation density matrices of the Z boson and W ± in the processes pp → ZH and pp → W ± H. With the aim of making maximum use of the information about states of polarisation, we relate all of the elements of the polarisation density matrices with observables which we can measure at the environment of proton-proton (pp) collisions. We discuss consequences of non-standard HZZ, HZγ and HW W interactions for these observables. A detailed analysis of the polarisation of the Z boson in the process e + e − → ZH is present in ref. [74]. Our approach in the process pp → ZH is an extension of that work to pp collisions.
The paper is organised as follows. In Section 2, we introduce all the ingredients needed for our analyses in the following sections. At first, we give non-standard HZZ, HZγ and HW W couplings. We present the complete helicity amplitudes for the sub-processes qq → ZH, ud → W + H and dū → W − H. These amplitudes are given for the non-standard couplings. We derive the relations between the helicity amplitudes imposed by the CP and CP T symmetries. In Section 3, we analyse the polarisation density matrices of the Z boson in the process pp → ZH in detail. First of all, we define the polarisation density matrices consisting of the helicity amplitudes of the previous section. Then we derive the 4 different differential cross sections for the process pp → ZH followed by Z → ff with respect to the Z decay angles. Among the 36 coefficients (= 9 × 4) of these 4 different differential angular distributions, only 15 coefficients can be non-zero. These coefficients are written in terms of the elements of the polarisation density matrices. The restrictions on these coefficients imposed by the CP and CP T symmetries are clarified. We focus on the coefficients which are strictly zero in the SM due to CP invariance or small in the SM due to the smallness of re-scattering effects, and study the influences of the non-standard HZZ and HZγ couplings. In Section 4, we analyse the polarisation density matrices of the W + and W − in the process pp → W ± H following the same procedure as the previous section. Section 5 gives a summary.

Constituents of the polarisation density matrices
As the base vectors for the density matrices 2 , we choose the eigenfunctions of the helicity operator. The density matrices are, therefore, constructed of helicity amplitudes. In this section we derive the complete helicity amplitudes for the sub-processes qq → ZH, ud → W + H and dū → W − H.

Effective Lagrangian with CP-odd operators
We derive non-standard HZZ, HZγ and HW W couplings from the following effective Lagrangian [74,75]: sin θ w , g Z = g/ cos θ w with θ w being the electroweak mixing angle and e being the proton charge, and 0123 = +1 in our convention. All of the eleven coefficients a Z , a W , b V , c V and b V (V = Z, γ, W ) are real so that L eff = L † eff and zero at the tree level in the SM. The operators whose coefficients are a Z , a W , b V and c V (V = Z, γ, W ) are CP-even and the operators whose coefficients are b V (V = Z, γ, W ) are CP-odd. The Lagrangian with only the CP-even operators indicates that H is a CP-even scalar field, while the Lagrangian with only the CP-odd operators indicates that H is a CP-odd scalar field. In both cases, the theory is CP invariant. If both the CP-even operator(s) and the CP-odd operator(s) exist, H is no longer a CP eigenstate and the theory is not CP invariant. We call a Z , a W , b V and c V (V = Z, γ, W ) CP-even coefficients and b V (V = Z, γ, W ) CP-odd coefficients.
If we assign the momenta as shown in Figure 1 3 , Feynman rules for the HZV (V = Z, γ) vertex and the HW W vertex can be expressed as Terms proportional to k ν vanish for the on-shell Z boson and the on-shell W boson, and terms proportional to p µ also vanish in our processes where the intermediate off-shell vector boson (Z, γ or W ) couples to the fourvector consisting of the two massless quarks. These terms are, therefore, not included in the above formulae. All the form factors are constant and expressed in terms of the coefficients in the effective Lagrangian in eq. (2.1): The form factors h V i (V = Z, γ, W ) (i = 1, 2, 3) consist of the CP-even coefficients and h V 4 (V = Z, γ, W ) consist of the CP-odd coefficients. We call h V i (V = Z, γ, W ) (i = 1, 2, 3) CP-even form factors and h V 4 (V = Z, γ, W ) CP-odd form factors.

Helicity amplitudes
Feynman diagrams for the sub-processes qq → ZH, ud → W + H and dū → W − H are shown in Figure 2. The circles denote the non-standard couplings derived in the previous section. We evaluate production helicity amplitudes in the following two frames. Let us assume that a direction of the z-axis is fixed along pp collisions. In one frame, which we call the qq ( ) centre-of-mass (c.m.) frame and is shown in the left picture of Figure 3, the quark (q) moves along the positive direction of the z-axis. In another frame, which we call the q ( ) q c.m. frame and is shown in the right picture of Figure 3, the antiquark (q ( ) ) moves along the positive direction of the z-axis. In Figure 3, our notation for helicities and the polar angle θ are also shown. We neglect the masses of q and q ( ) , thus the helicity of q ( ) is always opposite to that of q in our sub-processes. Helicity amplitudes are, therefore, given for the quark helicity σ(= ±1) and the weak boson helicity λ(= ±1, 0); the antiquark helicity is automatically fixed to −σ 4 . We denote the helicity amplitudes evaluated in q, σ The production amplitudes M λ σ (qq ( ) ) are evaluated in this frame. Right: the q ( ) q c.m. frame, where the antiquark (q ( ) ) moves along the positive direction of the z-axis. The production amplitudes M λ σ (q ( ) q) are evaluated in this frame. The quark helicity σ(= ±1), the antiquark helicity −σ and the weak boson helicity λ(= ±1, 0) are shown, and the polar angle θ of the weak boson V from the z-axis is denoted in both the frames.
the qq ( ) c.m. frame by M λ σ (qq ( ) ) and those evaluated in the q ( ) q c.m. frame by M λ σ (q ( ) q). It should be emphasised that σ always denotes the helicity of the quark and θ is the polar angle of the weak boson from the z-axis (not the quark momentum direction). The helicity amplitudes M λ σ (qq) for the sub-process qq → ZH are Here √ŝ is the qq c.m. energy, w is the energy of the Z boson: w = (ŝ + m 2 Z − m 2 H )/(2 √ŝ ), k is the momentum of the Z boson: k = w 2 − m 2 Z , Q q is the electric charge of the quark in units of e, g Z+ = g Z (−Q q sin 2 θ w ), and g Z− = g Z (T 3 q −Q q sin 2 θ w ) where T 3 u = 1/2 and T 3 d = −1/2.
Although the actual calculation of the helicity amplitudes M λ σ (qq) is easy, it is also possible to estimate them from M λ σ (qq) as follows. If we denote the helicity ofq by σ and the helicity of q by −σ (this is opposite to our helicity notation), the angular part is the same as M λ σ (qq) and the amplitudes are given by where the helicity of each particle is explicitly shown in parenthesis. Because the helicity of a quark is equal to its chirality and the helicity of an antiquark is opposite to its chirality in the massless limit, q with helicity −σ has the chirality −σ. Hence the coupling g Zσ must be replaced by g Z−σ in eq. (2.5), and this replacement is expressed byM λ −σ in eq. (2.6). By the simple replacement σ → −σ in eq. (2.6), we obtain the helicity amplitudes M λ σ (qq) in our helicity notation: (2.7b) We note that the relative sign among the amplitudes in eq. (2.4) and that among the amplitudes in eq. (2.7) are important, because they appear in the off-diagonal elements of the density matrices.
Similarly, the helicity amplitudes for the sub-process ud → W + H are given by and those for the sub-process dū → W − H are given by Here w is the energy of the W boson: , k is the momentum of the W boson: k = w 2 − m 2 W , and V ud is the element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix.
By looking at the helicity amplitudes, we can already discuss a difference between the Z boson and the W boson with regard to states of polarisation. Since the W boson couples to only a fermion with chirality −1, the initial quark u or d always has the helicity σ = −1. Therefore, while the Z boson is in a partially polarised state (often called a mixed state; see e.g. refs. [57][58][59][60]), the W boson (both the W + and W − ) is in a completely polarised state (often called a pure state).

Requests from symmetries
Conditions imposed by symmetries lead to certain relations between the helicity amplitudes. The upper figures in Figure 4 show the original states in the process qq → ZH (left), the states after CP transformation (middle) and the states after the rotation around the z-axis by π (right). From these figures, we find that the invariance under CP leads to the relation where M λ σ (qq)(θ) = M λ σ (qq) given in eq. (2.4), and M −λ σ (qq)(π−θ) is obtained by setting θ → π−θ in M −λ σ (qq). It is easy to see that a non-zero value of the CP-odd form factors h Z 4 and/or h γ 4 violates this relation. Similarly, the lower figures in Figure 4 show the original states in the process ud → W + H (left), the states after CP transformation (middle) and the states after the rotation around the z-axis by π (right). From these figures, we find that CP invariance requests the relation where M λ − (ud)(θ) = M λ − (ud) given in eq. (2.8), and M −λ − (dū)(π − θ) is obtained by setting θ → π − θ in M −λ − (dū) given in eq. (2.9). This relation is violated by the imaginary part of the element V ud of the CKM matrix, even when the CP-odd form factor h W 4 is zero. However, this CP violation phase in V ud [77] is always regarded as an overall common phase among the amplitudes. For instance, the phase in V ud is an overall common phase among M λ=+ . Therefore, the effect of the CP violation phase in the CKM matrix never appears even in the off-diagonal elements of the density matrices. Only a non-zero value of the CP-odd form factor h W 4 violates the relation in eq. (2.12) and affects a state of polarisation of the W boson. We can write the relations in eqs. (2.11) and (2.12) in the following convenient forms, respectively: CPT invariance and the unitarity condition conclude [61] that the following relations hold at the tree level approximation and violation of the relations indicates the existence of re-scattering effects: It is easy to see that the relations hold if all of the form factors are real. Following ref. [61], we call the invariance which leads to these relations CP T invariance. We can write the relations in eq. (2.14) in the following convenient forms, respectively: 3 Polarisation of the Z boson

Polarisation density matrices of the Z boson
In this section we define the density matrices of the Z boson by using our notation of the helicity amplitudes of the previous section. We consider the sub-process where the helicity of each particle is shown in parenthesis. We neglect the masses of the final fermion f and the final antifermionf , thus the helicity off is always opposite to that of f . We express the full helicity amplitude as where the production amplitude M λ σ (qq) is given in eq. (2.4), D τ λ is the decay helicity amplitude, and denotes the propagator factor of the Z boson. We evaluate the decay amplitude in the following four-momentum frame: The decay amplitude is The explicit forms of the coupling g τ Zff are, g + Zll = g Z sin 2 θ w and g − Zll = g Z (−1/2 + sin 2 θ w ) for a charged lepton pair, g + represents the elements of the density matrix in the helicity basis of the Z boson in the qq c.m. frame, and at the last equality the following 3 × 3 matrix form is employed: The density matrix is a 3 × 3 Hermitian matrix: σ ρ σ (qq) = ( σ ρ σ (qq)) † 5 .
The full helicity amplitude T τ σ (qq), in which the production amplitude is evaluated in theqq c.m. frame shown in the right picture of Figure 3, can be treated in the same manner: represents the elements of the density matrix in the helicity basis of the Z boson in theqq c.m. frame. The production amplitude M λ σ (qq) is given in eq. (2.7).

Label
Integration over y Integration over cos θ Table 1: The 4 different approaches labelled A, B, C and D for performing integration over the rapidity y and the polar angle cos θ in the differential cross section of eq. (3.12). A value of y cut is ln s/ŝ; see eq. (3.14). A smaller value can be also chosen for y cut . A value of (0 ≤ < 1) will be determined according to experimental conditions.

Decay angular distributions of the polarised Z boson
In terms of the density matrices defined in eqs. (3.8) and (3.11), the complete differential cross section for the process pp → ZH followed by Z → ff in the narrow width approximation can be expressed as follows: where s is the c.m. energy squared of the pp collisions,ŝ is the c.m. energy squared of the qq collisions, y is the rapidity of the qq c.m. frame and theqq c.m. frame (see Figure 3) on the pp c.m. frame (i.e. the experimental frame), C f is an effective colour factor: C f = 1 for Z decays into a charged lepton pair and C f = 3 for Z decays into a quark pair, q(x i ) andq(x i ) are the quark and antiquark parton distribution functions (PDFs) with energy fraction x i . Summations are performed over the final fermion flavor f , the helicity τ of the final fermion, and the initial quark and antiquark flavor q. Averages are performed over the helicity and colour of the initial quark, and those of the initial antiquark. Recall that k is the momentum of the Z boson and θ is the polar angle of the Z boson in the qq c.m. frame and theqq c.m. frame. The decay angles θ and φ are defined in eq. (3.4). Here x 1 (x 2 ) denotes the energy fraction of the initial quark or antiquark in the proton that moves along the positive (negative) direction of the z-axis, once a direction of the z-axis is fixed along the pp collisions. The integration variablesŝ and y are related to x 1 and The allowed region is (3.14) Now we derive the differential cross sections with respect toŝ, cos θ and φ by integrating over y and cos θ in eq. (3.12). We consider the 4 different integration approaches summarised in Table 1. In the complete differential cross section of eq. (3.12), the polar angle θ dependence appears only in ρ σ (qq) and ρ σ (qq). The cos θ integration in the approaches A and B results in where c i (i = 1, 2, 3) are constant values depending on and summarised in Table 2, andM λ=±,0 σ are defined in eq. (2.5). Notice the difference between eq. (3.15a) and eq. (3.15b) that there is a minus sign in front of the elements that have the overall σ in eq. (3.15b). This result is actually obvious from the comparison between the amplitudes in eq. (2.4) and those in eq. (2.7). After integration over y and cos θ in the approach A, the differential cross section can be expressed as is introduced to simplify our writing. The PDFs labelled A in the second equality (i.e. q(x 1 )q(x 2 ) and q(x 2 )q(x 1 )) give the same numerical contribution after integration over y, therefore they are combined in the third (i.e. last) equality. The same is done for the PDFs labelled B. In the last equality, ρ σ (qq) and ρ σ (qq) are added, and the elements that have the overall σ vanish due to the sign difference in ρ σ (qq) and ρ σ (qq) . The vanished elements in the approach A in eq. (3.16) revive by performing integration in the approach B: where, in the last equality, ρ σ (qq) and ρ σ (qq) are subtracted. As a result, in contrast to the approach A, only the elements that have the overall σ survive.
The cos θ integration in the approaches C and D results in where the constant values c i (i = 4, 5) are summarised in Table 2. Notice again that the sign in front of the elements that have the overall σ is different between ρ σ (qq) and ρ σ (qq) . The integration over y and cos θ in the approaches C and D proceeds in the same manners as eq. (3.16) and eq. (3.18), respectively: (3.21) Only the elements that do not have the overall σ survive in eq. (3.20), while only the elements that have the overall σ survive in eq. (3.21).
The differential cross sections with respect toŝ and the decay angles θ and φ have 9 independent angular distributions and can be expressed as where the coefficients F ia (i = A, B, C, D) (a = 1, 2, · · · , 9) are functions ofŝ and written in terms of the non-vanishing elements in eqs. (3.16), (3.18), (3.20) and (3.21). Note that there are in total 36 (= 4 × 9) coefficients. It is straightforward to obtain the explicit form of the coefficients F ia (i = A, B, C, D) (a = 1, 2, · · · , 9): Among the 36 coefficients, only the 15 coefficients can be non-zero. It is easy to notice that there are the 10 combinations of the elements of the density matrix in total. However, only 9 of them are independent 6 . Some of them are strictly zero, if the amplitudes satisfy the restriction in eq. (2.13a) from CP invariance and/or the restriction in eq. (2.15a) from CP T invariance. These symmetry properties can be explicitly checked in the following way. By applying the restrictions in eqs. (2.13a) and (2.15a) to the combination in f A9 , we find where the former means that CP invariance requires it to be zero, while the latter means that CP T invariance does not require it to be zero. This indicates that observation of a nonzero value in F A9 signals CP violation. Similarly, by applying the restrictions in eqs. (2.13a) and (2.15a) to the combination in f B8 and f C7 , we find where CP T invariance requires it to be zero, while CP invariance does not. This indicates that observation of a nonzero value in F B8 or F C7 signals CP T violation. Finally, by applying the restrictions in eqs. (2.13a) and (2.15a) to the combination in f B5 and f C4 , we find where both CP invariance and CP T invariance require it to be zero. This indicates that observation of a nonzero value in F B5 or F C4 signals both CP violation and CP T violation. These symmetry properties of the 10 combinations are summarised in Table 3. The symbol + means that the symmetry does not require the combination to be zero, while − means that the symmetry requires the combination to be zero. We also show the coefficients F ia by measuring which the combinations of the density matrix elements can be determined. The symbol • in the column "f charge" means that observation of the coefficient F ia requires distinguishing the fermion f from the antifermionf . In other words, it requires identification of the charge (flavor) of the final fermion f . This can be confirmed by performing a translation θ → π − θ and φ → φ + π in eq. (3.22) and observing the change of the sign. Among the 15 coefficients, only 9 of them do not require the charge identification of the final fermion f . It should be emphasised that these 9 coefficients are necessary and sufficient to determine all of the 9 independent combinations of the density matrix elements.
The coefficients F ia (i = A, B, C, D) (a = 1, 2, · · · , 9) will be experimentally determined by measuring the decay angles θ and φ in the rest frame of the Z boson (see eq. (3.4)), since we completely know the differential angular distributions as eq. (3.22). We just do not know the coefficients F ia which uniquely depend on a state of polarisation of the Z boson. With appropriate integration over cos θ and φ, it is possible to isolate the angular distributions:  Table 3: Symmetry properties of the 10 combinations of the density matrix elements. The coefficients F ia of the differential angular distributions, by measuring which the combinations of the density matrix elements can be determined, are also shown. The symbol − means that the symmetry (CP or CP T) requires the combination to be zero, while the symbol + means that the symmetry does not; observation of a nonzero value in the combination with the symbol − under CP, for instance, signals CP violation. The symbol • in the column "f charge" means that observation of the coefficient F ia requires the charge (or flavor) identification of the final fermion f . The symbol • in the column "HZγ" indicates that the coefficient F ia has a good sensitivity to the HZγ coupling; see a discussion at the paragraph of eq. (3.30). The symbol • in the last column indicates that the coefficient F ia is weakened according to the degree of polarisation of the Z boson; see a discussion at the last paragraph of Section 4.3. By the simple replacementsM λ σ →N λ and F ia → F W ia , the table corresponding to the process pp → W ± H is immediately obtained; see Section 4.   3,4,5,6,7,8,9). For example, eqs. (3.28b) and (3.29b) gives F 5 . Only F 2 is not determined in this method. By a fitting procedure in the differential cross section with respect toŝ and cos θ in eq. (3.29a), F 2 may be determined. and −F C4 (right panel) divided by F A1 are shown in the same manner as Figure 5. These coefficients are doubly constrained to be identically zero by CP invariance and CP T invariance. In each panel, 2 points give predictions for h Z 4 = 0 + 0.1i and × points give those for h γ 4 = 0 + 0.1i.

Influences of non-standard HZZ and HZγ interactions
In the previous section, we have clarified the restrictions on the coefficients of the differential angular distributions imposed by the CP and CP T symmetries; see Table 3. Some of the coefficients are strictly zero in the SM due to CP invariance and some of them are small in the SM due to the smallness of re-scattering effects. These coefficients are in particular interesting as observables at the LHC, since observation of a non-zero or large value in these coefficients immediately signals the existence of physics beyond the SM. In this Section, we focus on these coefficients and study the influences of the non-standard HZZ and HZγ couplings.
Our numerical results are produced for the 14 TeV LHC. We set m H = 125.5 GeV, and = 0 in the cos θ integration; see eqs. (3.15) and (3.19), and Table 2. As the final fermion flavor f summed in eq. (3.17), only the electron and the muon are considered when calculating the coefficients F ia (i = A, B, C, D) (a = 3, 4, 7) (Recall that observation of these coefficients requires the charge identification of f ), and all the quark flavors but the top quark are additionally considered when calculating the other coefficients. MSTW PDFs [78] are used. The phase space integration is performed with the program BASES [79].
The coefficients F C8 and F B7 measure the same combination of the density matrix elements (see Table 3), thus show the same dependence onŝ 1/2 . However, F C8 's sensitivity to h Z 4 is stronger than its sensitivity to h γ 4 , while F B7 's sensitivity to h Z 4 is comparable to its sensitivity to h γ 4 . This difference originates from the overall σ in f B7 and is explained as follows. Eq. (2.5) shows thatM λ + andM λ − can be written in the following simplified form as the sum of the Z boson contribution and the γ contribution: By using this expression, we derive The couplings g Z+ and g Z− have opposite signs to each other, both for the up-type quarks and for the down-type quarks. The interference term between the Z boson contribution and the γ contribution in eq. (3.31a) and that in eq. (3.31b), therefore, have opposite signs to each other, too. Considering f λ γ = 0 at the tree level approximation in the SM, the interference terms may give a larger contribution than the (f λ γ ) * f λ γ term. The coefficient F C8 has a structure that eq. (3.31a) and eq. (3.31b) are added and the interference terms tend to cancel to each other. On the other hand, F B7 has a structure that eq. (3.31b) and eq. (3.31a) are subtracted due to the overall σ in f B7 , and the interference terms contribute in the same direction. Therefore, as long as f λ γ (i.e the HZγ coupling) is small, the effects of it appear larger in F B7 than in F C8 . Notice that this discussion is independent on explicit choices for λ and λ . The above finding is true not only in F B7 but also in the other coefficients which have the overall σ in f ia . These coefficients are indicated by the symbol • in the column "HZγ" in Table 3.
In Figure 6, the coefficients F B8 (left panel) and F C7 (right panel) divided by F A1 are shown in the same manner as Figure 5. CP T invariance requires these 2 coefficients to be identically zero, hence observation of a non-zero value indicates the existence of re-scattering effects. Re-scattering effects can be approximately included by allowing imaginary parts in the form factors [61]. In each panel, 2 points give predictions for h Z 3 = 0 + 0.1i, and × points give those for h γ 3 = 0 + 0.1i.
In Figure 7, the coefficients F A3 (left panel), F D1 (middle-left panel), −F B5 (middle-right panel) and −F C4 (right panel) divided by F A1 are shown in the same manner as Figure 5. These 4 coefficients are doubly constrained to be identically zero by CP invariance and CP T invariance: they are strictly zero if CP is conserved, even if CP T is violated, for instance. Observation of a non-zero value indicates CP violation and the existence of re-scattering effects. In each panel, 2 points give predictions for the CP-odd from factor h Z 4 = 0 + 0.1i, and × points give those for the CP-odd form factor h γ 4 = 0 + 0.1i.

Polarisation of the W boson
In this section, we analyse the density matrices of the W + and W − in the process pp → W ± H following the same procedure as Section 3. Because (1) a determination of the density matrix of the W boson is difficult in pp collisions when the W boson decays into a charged lepton and a neutrino and (2) we cannot distinguish W + from W − in view of the difficulty of flavor identification of both W + and W − decay products when the W boson decays into two quarks, we consider the process pp → W ± H (W ± → jj) as the sum of the process pp → W + H (W + → jj) and the process pp → W − H (W − → jj). Just as the process pp → ZH, only the 15 coefficients among the 36 coefficients of the 4 different differential angular distributions of the dijets can be non-zero. However only 9 coefficients among these 15 coefficients are actually measurable.

Polarisation density matrices of the W + and W −
First of all, we define the density matrices of the W + and W − . We consider the sub-processes and where the helicity of each particle is shown in parenthesis. We neglect the masses of the quark and the antiquark from the W decay, hence the helicity of the antiquark is always opposite to that of the quark. We prepare the 4 full helicity amplitudes T (ud), T (dū), T (du) and T (ūd) in which the production amplitude is evaluated in the ud c.m. frame, the dū c.m. frame, thedu c.m. frame and theūd c.m. frame, respectively. Recall that these c.m. frames are shown in Figure 3 and the production amplitudes are given in eqs. (2.8) and (2.9). These 4 full helicity amplitudes are given by denotes the propagator factor of the W boson, D λ is the helicity amplitude for the decay process W + → ud and D λ is the helicity amplitude for the decay process W − → dū: where d − λ is given in eq. (3.6). The decay amplitudes are evaluated in the following four-momentum frame: Then we obtain represent the elements of the density matrix in the helicity basis of the W + in the ud c.m. frame and those in thedu c.m. frame, respectively, and The relations indicate that the W + and W − are always in the same state of polarisation. This fact makes our analysis in the next section simpler, because we do not have to distinguish the W + from the W − with regard to states of polarisation.
4.2 Decay angular distributions of the polarised W + and W − The complete differential cross section for the process pp → W + H followed by W + → jj in the narrow width approximation can be expressed in terms of the density matrices in eq. (4.7) as and that for the process pp → W − H followed by W − → jj can be expressed in terms of the density matrices in eq. (4.8) as where the 3×3 matrix form in eq. (3.9) is used. For the definition of variables, see below eq. (3.12).
As the quark flavors u and d summed in the above equations, we consider all the quark flavors but the top quark. As a result, the unitarity of the CKM matrix gives u,d |V ud | 2 = 2. By using the relations in eq. (4.9), we obtain the complete differential cross section for the process pp → W ± H followed by W ± → jj as the sum of the above two cross sections in the following compact form: (4.12) As in Section 3.2, we derive the differential cross section with respect toŝ, cos θ and φ by integrating over y and cos θ in eq. (4.12). It is straightforward to confirm that the integration over y and cos θ in the 4 different approaches summarised in Table 1 can be performed in the same manner as eqs. (3.16), (3.18), (3.20) and (3.21). The resulting differential cross sections are   where the constant values c i (i = 1, 2, 3, 4, 5) are summarised in Table 2, andN λ=±,0 are defined in eq. (2.10). By comparing the above 4 differential cross sections with eq. (3.22), we obtain the 36 coefficients F W ia (i = A, B, C, D) (a = 1, 2, · · · , 9) in total: where f W ia (i = A, B, C, D) (a = 1, 2, · · · , 9) are obtained by the following replacements in those coefficients of the process pp → ZH in eq. (3.24): (4.16) Table 3 corresponding to the process pp → W ± H is obtained by the following simple replacements in Table 3: Just as the process pp → ZH, there are in total 15 coefficients which can be non-zero. Observation of the 6 coefficients F W has difficulty, since it requires the charge (or flavor) identification of the parent quark of the jet. Only the remaining 9 coefficients are, therefore, actually measurable. We emphasise that these 9 coefficients are necessary and sufficient to determine all of the 9 independent combinations of the density matrix elements.

Influences of non-standard HW W interaction
In this Section, with the same motivation as Section 3.3, we focus on the measurable coefficients which are strictly zero or small in the SM, and study the influences of the non-standard HW W coupling. In the left panel of Figure 8, the 2 coefficients −F W A9 (2) and F W C8 (×) divided by F W

A1
are shown for the CP-odd form factor h W 4 = 0.1. CP invariance requires these 2 coefficients to be identically zero; observation of a non-zero value signals CP violation. In the middle panel of Figure 8, the coefficient F W B8 divided by F W A1 is shown for the CP-even form factor h W 3 = 0 + 0.1i.
CP T invariance requires this coefficient to be identically zero; observation of a non-zero value indicates the existence of re-scattering effects. In the right panel of Figure 8, the 2 coefficients F W
We discuss the differences between the coefficients F ia in the process pp → ZH and those F W ia in the process pp → W ± H. From the comparisons between the results for the non-standard HZZ coupling in Figures 5, 6 and 7 and the results for the non-standard HW W coupling in Figure 8, we notice that F W A9 and F W C8 are comparable with F A9 and F C8 , respectively, while F W B8 , F W D1 and F W B5 are consistently larger than F B8 , F D1 and F B5 , respectively. As we have discussed at the last paragraph of Section 2.2, the Z boson is in a partially polarised state, while the W + and W − are in a completely polarised state. The degree of polarisation affects the magnitudes of the coefficients, hence it is expected that the coefficients F W ia are always equal or larger than the coefficients F ia . The existence of the overall σ in F B8 , F D1 and F B5 (see eq. (3.24)) indicates that these coefficients are weakened according to the degree of polarisation of the Z boson. These coefficients are denoted by the symbol • in the last column in Table 3.

Summary
The measurements of the Higgs boson couplings to the SM particles are essential tests of the SM. The Z boson in the process pp → ZH and the W + and W − in the process pp → W ± H can be in polarised states, and it would be possible to study the Higgs boson couplings to the weak bosons (HZZ, HZγ and HW W ) in detail from a careful analysis of these states of polarisation. A density matrix contains the complete information about a state of polarisation, and all of the elements of the density matrix should be made use of in such a careful analysis. In this paper, such an analysis approach has been presented.
Determination of the density matrix of the Z (W ) boson requires measurements of the angular distributions of the Z (W ) decay products. In pp collisions, this is difficult when the W boson decays into a charged lepton and a neutrino. When the W boson decays into two quarks, we cannot distinguish W + from W − in view of the difficulty of flavor identification of both W + and W − decay products. Therefore, we have considered the process pp → W ± H (W ± → jj) as the sum of the process pp → W + H (W + → jj) and the process pp → W − H (W − → jj). We have found that the W + and W − are always in the same state of polarisation (i.e. the density matrix of the W + is always the same as that of W − ). By using this fact, we have developed an analysis approach which can be applied both to pp → ZH (Z → ff ) and to pp → W ± H (W ± → jj) in the same manner.
We have derived the 4 different differential cross sections with respect toŝ (the invariant mass squared of the Z (W ) boson and the Higgs boson) and cos θ and φ (the Z (W ) decay angles), by integrating the complete differential cross section over the other phase space variables in the 4 different approaches. Among the 36 coefficients (= 9 × 4) of these 4 differential angular distributions, only 15 coefficients can be non-zero (these coefficients are summarised in Table 3). These 15 coefficients are written in terms of the elements of the density matrices, and there exist 9 independent combinations of the elements of the density matrices. Observation of the 9 coefficients among the 15 coefficients does not require the charge (or flavor) identification of the Z (W ) decay products. In the analysis of the W + and W − , only these 9 coefficients are measurable. We have found that these 9 coefficients are necessary and sufficient to determine all of the 9 independent combinations of the density matrix elements (i.e. one coefficient corresponds to one combination). We have clarified the restrictions on the 15 coefficients imposed by the CP and CP T symmetries. Some of the coefficients are required to be identically zero by CP invariance; observation of a non-zero value in these coefficients signals CP violation. Similarly, some of the coefficients are required to be identically zero by CP T invariance; observation of a non-zero value in these coefficients indicates the existence of re-scattering effects. These coefficients are in particular interesting as observables at the LHC, since observation of a non-zero or large value in these coefficients immediately signals the existence of non-standard HZZ, HZγ and/or HW W interactions.