Test of CP invariance in vector-boson fusion production of the Higgs boson using the Optimal Observable method in the ditau decay channel with the ATLAS detector

A test of CP invariance in Higgs boson production via vector-boson fusion using the method of the Optimal Observable is presented. The analysis exploits the decay mode of the Higgs boson into a pair of τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau $$\end{document} leptons and is based on 20.3 fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{fb}^{-1}$$\end{document} of proton–proton collision data at s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s}$$\end{document} = 8 TeV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\,\mathrm{TeV}$$\end{document} collected by the ATLAS experiment at the LHC. Contributions from CP-violating interactions between the Higgs boson and electroweak gauge bosons are described in an effective field theory framework, in which the strength of CP violation is governed by a single parameter d~\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{d}$$\end{document}. The mean values and distributions of CP-odd observables agree with the expectation in the Standard Model and show no sign of CP violation. The CP-mixing parameter d~\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{d}$$\end{document} is constrained to the interval (-0.11,0.05)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(-0.11,0.05)$$\end{document} at 68% confidence level, consistent with the Standard Model expectation of d~=0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tilde{d}=0$$\end{document}.


Introduction
The discovery of a Higgs boson by the ATLAS and CMS experiments [1,2] at the LHC [3] offers a novel opportunity to search for new sources of CP violation in the interaction of e-mail: atlas.publications@cern.ch the Higgs boson with other Standard Model (SM) particles. C and CP violation is one of the three Sakharov conditions [4][5][6] needed to explain the observed baryon asymmetry of the universe. In the SM with massless neutrinos the only source of CP violation is the complex phase in the quark mixing (CKM) matrix [7,8]. The measured size of the complex phase and the derived magnitude of CP violation in the early universe is insufficient to explain the observed value of the baryon asymmetry [9] within the SM [10,11] and, most probably, new sources of CP violation beyond the SM need to be introduced. No observable effect of CP violation is expected in the production or decay of the SM Higgs boson. Hence any observation of CP violation involving the observed Higgs boson would be an unequivocal sign of physics beyond the SM.
The measured Higgs boson production cross sections, branching ratios and derived constraints on coupling-strength modifiers, assuming the tensor structure of the SM, agree with the SM predictions [12,13]. Investigations of spin and CP quantum numbers in bosonic decay modes and measurements of anomalous couplings including CP-violating ones in the decay into a pair of massive electroweak gauge bosons show no hints of deviations from the tensor structure of the SM Higgs boson [14,15]. Differential cross-section measurements in the decay H → γ γ have been used to set limits on couplings including CP-violating ones in vector-boson fusion production in an effective field theory [16]. However, the observables, including absolute event rates, used in that analysis were CP-even and hence not sensitive to the possible interference between the SM and CP-odd couplings and did not directly test CP invariance. The observables used in this analysis are CP-odd and therefore sensitive to this interference and the measurement is designed as a direct test of CP invariance.
In this paper, a first direct test of CP invariance in Higgs boson production via vector-boson fusion (VBF) is presented, based on proton-proton collision data corresponding to an integrated luminosity of 20.3 fb −1 collected with the ATLAS detector at √ s = 8 TeV in 2012. A CP-odd Optimal Observable [17][18][19] is employed. The Optimal Observable combines the information from the multi-dimensional phase space in a single quantity calculated from leadingorder matrix elements for VBF production. Hence it does not depend on the decay mode of the Higgs boson. A direct test of CP invariance is possible measuring the mean value of the CP-odd Optimal Observable. Moreover, as described in Sect. 2, an ansatz in the framework of an effective field theory is utilised, in which all CP-violating effects corresponding to operators with dimensions up to six in the couplings between a Higgs boson and an electroweak gauge boson can be described in terms of a single parameterd. Limits ond are derived by analysing the shape of spectra of the Optimal Observable measured in H → τ τ candidate events that also have two jets tagging VBF production. The event selection, estimation of background contributions and of systematic uncertainties follows the analysis used to establish 4.5σ evidence for the H → τ τ decay [20]. Only events selected in the VBF category are analysed, and only fully leptonic τ lep τ lep or semileptonic τ lep τ had decays of the τ -lepton pair are considered.
The theoretical framework in the context of effective field theories is discussed in Sect. 2 and the methodology of testing CP invariance and the concept of the Optimal Observable are introduced in Sect. 3. After a brief description of the ATLAS detector in Sect. 4, the simulated samples used are summarised in Sect. 5. The experimental analysis is presented in Sect. 6, followed by a description of the statistical method used to determine confidence intervals ford in Sect. 7. The results are discussed in Sect. 8, following which conclusions are given.

Effective Lagrangian framework
The effective Lagrangian considered is the SM Lagrangian augmented by CP-violating operators of mass dimension six, which can be constructed from the Higgs doublet and the U(1) Y and SU(2) I W ,L electroweak gauge fields B μ and W a,μ (a = 1, 2, 3), respectively. No CP-conserving dimension-six operators built from these fields are taken into account. All interactions between the Higgs boson and other SM particles (fermions and gluons) are assumed to be as predicted in the SM; i.e. the coupling structure in gluon fusion production and in the decay into a pair of τ -leptons is considered to be the same as in the SM.
The effective U(1) Y -and SU(2) I W ,L -invariant Lagrangian is then given by (following Refs. [21,22]): with the three dimension-six operators and three dimensionless Wilson coefficients fB B , fW W and fB; is the scale of new physics. Here Only two of the four couplingsg H V V (V = W ± , Z , γ ) are independent due to constraints imposed by U(1) Y and SU(2) I W ,L invariance. They can be expressed in terms of two dimensionless couplingsd andd B as: Hence in general W W , Z Z, Z γ and γ γ fusion contribute to VBF production. The relations betweend and fW W , andd B and fB B are given by: As the different contributions from the various electroweak gauge-boson fusion processes cannot be distinguished experimentally with the current available dataset, the arbitrary choiced =d B is adopted. This yields the following relation for theg H V V : The parameterd is related to the parameterκ W = κ W /κ SM  The effective Lagrangian yields the following Lorentz structure for each vertex in the Higgs bosons coupling to two identical or charge-conjugated electroweak gauge bosons H V ( p 1 )V ( p 2 ) (V = W ± , Z , γ ), with p 1,2 denoting the momenta of the gauge bosons: The first terms (∝ g μν ) are CP-even and describe the SM coupling structure, while the second terms (∝ ε μνρσ p 1ρ p 2σ ) are CP-odd and arise from the CP-odd dimension-six operators. The choiced =d B gives the same coefficients multiplying the CP-odd structure for H W + W − , H Z Z and H γ γ vertices and a vanishing coupling for the H Zγ vertex. The matrix element M for VBF production is the sum of a CP-even contribution M SM from the SM and a CPodd contribution M CP-odd from the dimension-six operators considered: The differential cross section or squared matrix element has three contributions: The first term |M SM | 2 and third termd 2 · |M CP-odd | 2 are both CP-even and hence do not yield a source of CP violation. The second termd · 2Re(M * SM M CP-odd ), stemming from the interference of the two contributions to the matrix element, is CP-odd and is a possible new source of CP violation in the Higgs sector. The interference term integrated over a CP-symmetric part of phase space vanishes and therefore does not contribute to the total cross section and observed event yield after applying CP-symmetric selection criteria. The third term increases the total cross section by an amount quadratic ind, but this is not exploited in the analysis presented here.

Test of CP invariance and Optimal Observable
Tests of CP invariance can be performed in a completely model-independent way by measuring the mean value of a CP-odd observable O CP . If CP invariance holds, the mean value has to vanish O CP = 0. An observation of a nonvanishing mean value would be a clear sign of CP violation. A simple CP-odd observable for Higgs boson production in VBF, the "signed" difference in the azimuthal angle between the two tagging jets φ j j , was suggested in Ref. [22] and is formally defined as: Here b μ + and b μ − denote the normalised four-momenta of the two proton beams, circulating clockwise and anticlockwise, and p μ + (φ + ) and p μ − (φ − ) denote the fourmomenta (azimuthal angles) of the two tagging jets, where p + ( p − ) points into the same detector hemisphere as b This ordering of the tagging jets by hemispheres removes the sign ambiguity in the standard definition of φ j j .
The final state consisting of the Higgs boson and the two tagging jets can be characterised by seven phase-space variables while assuming the mass of the Higgs boson, neglecting jet masses and exploiting momentum conservation in the plane transverse to the beam line. The concept of the Optimal Observable combines the information of the highdimensional phase space in a single observable, which can be shown to have the highest sensitivity for small values of the parameter of interest and neglects contributions proportional tod 2 in the matrix element. The method was first suggested for the estimation of a single parameter using the mean value only [17] and via a maximum-likelihood fit to the full distribution [18] using the so-called Optimal Observable of first order. The extension to several parameters and also exploiting the matrix-element contributions quadratic in the parameters by adding an Optimal Observable of second order was introduced in Refs. [19,27,28]. The technique has been applied in various experimental analyses, e.g. Refs. [15,[29][30][31][32][33][34][35][36][37][38][39].
The analysis presented here uses only the first-order Optimal Observable OO (called Optimal Observable below) for the measurement ofd via a maximum-likelihood fit to the full distribution. It is defined as the ratio of the interference term in the matrix element to the SM contribution: Figure 1 shows the distribution of the Optimal Observable, at parton level both for the SM case and for two non-zerod values, which introduce an asymmetry into the distribution and yield a non-vanishing mean value. The values of the leading-order matrix elements needed for the calculation of the Optimal Observable are extracted from HAWK [41][42][43]. The evaluation requires the fourmomenta of the Higgs boson and the two tagging jets. The momentum fraction x 1 (x 2 ) of the initial-state parton from the proton moving in the positive (negative) z-direction can be derived by exploiting energy-momentum conservation from the Higgs boson and tagging jet four-momenta as: where m Hjj (y Hjj ) is the invariant mass (rapidity) obtained from the vectorially summed four-momenta of the tagging jets and the Higgs boson. Since the flavour of the initialand final-state partons cannot be determined experimentally, the sum over all possible flavour configurations i j → kl H weighted by the CT10 leading-order parton distribution functions (PDFs) [44] is calculated separately for the matrix elements in the numerator and denominator:

The ATLAS detector
The ATLAS detector [45] is a multi-purpose detector with a cylindrical geometry. 1 It comprises an inner detector (ID) surrounded by a thin superconducting solenoid, a 1 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical coordinates (r, φ) are used in the transverse plane, φ being the azimuthal angle around the z-axis. The pseudorapidity is defined in terms of the polar angle θ as η = − ln tan(θ/2).
calorimeter system and an extensive muon spectrometer in a toroidal magnetic field. The ID tracking system consists of a silicon pixel detector, a silicon microstrip detector, and a transition radiation tracker. It provides precise position and momentum measurements for charged particles and allows efficient identification of jets containing b-hadrons (b-jets) in the pseudorapidity range |η| < 2.5. The ID is immersed in a 2 T axial magnetic field and is surrounded by high-granularity lead/liquid-argon sampling electromagnetic calorimeters which cover the pseudorapidity range |η| < 3.2. A steel/scintillator tile calorimeter provides hadronic energy measurements in the central pseudorapidity range (|η| < 1.7). In the forward regions (1.5 < |η| < 4.9), the system is complemented by two end-cap calorimeters using liquid argon as active material and copper or tungsten as absorbers. The muon spectrometer surrounds the calorimeters and consists of three large superconducting eight-coil toroids, a system of tracking chambers, and detectors for triggering. The deflection of muons is measured in the region |η| < 2.7 by three layers of precision drift tubes, and cathode strip chambers in the innermost layer for |η| > 2.0. The trigger chambers consist of resistive plate chambers in the barrel (|η| < 1.05) and thin-gap chambers in the end-cap regions (1.05 < |η| < 2.4). A three-level trigger system [46] is used to select events. A hardware-based Level-1 trigger uses a subset of detector information to reduce the event rate to 75 kHz or less. The rate of accepted events is then reduced to about 400 Hz by two software-based trigger levels, named Level-2 and the Event Filter.

Simulated samples
Background and signal events are simulated using various Monte Carlo (MC) event generators, as summarised in Table 1. The generators used for the simulation of the hardscattering process and the model used for the simulation of the parton shower, hadronisation and underlying-event activity are listed. In addition, the cross-section values to which the simulation is normalised and the perturbative order in QCD of the respective calculations are provided.
All the background samples used in this analysis are the same as those employed in Ref.
[20], except the ones used to simulate events with the Higgs boson produced via gluon fusion and decaying into the τ τ final state. The Higgs-plus-one-jet process is simulated at NLO accuracy in QCD with Powheg-Box [47][48][49]73], with the MINLO feature [74] applied to include Higgs-plus-zero-jet events at NLO accuracy. This sample is referred to as HJ MINLO. The Powheg-Box event generator is interfaced to Pythia8 [51], and the CT10 [44] parameterisation of the PDFs is used. Higgs boson events produced via gluon fusion and decay-  [51]. For these simulated events, the shape of the generated p T distribution is matched to a NNLO + NNLL calculation HRes2.1 [82,83] in the inclusive phase space. Simultaneously, for events with two or more jets, the Higgs boson p T spectrum is reweighted to match the MINLO HJJ predictions [84]. The overall normalisation of the gluon fusion process (ggF) is taken from a calculation at next-to-next-to-leading order (NNLO) [75][76][77][78][79][80] in QCD, including soft-gluon resummation up to nextto-next-to-leading logarithm terms (NNLL) [85]. Next-toleading-order (NLO) electroweak (EW) corrections are also included [86,87]. Higgs boson events produced via VBF, with SM couplings, are also simulated with Powheg interfaced with Pythia8 (see Table 1 and Ref. [20]). Production by VBF is normalised to a cross section calculated with full NLO QCD and EW corrections [41,42,52] with an approximate NNLO QCD correction applied [53]. The NLO EW corrections for VBF production depend on the p T of the Higgs boson, and vary from a few percent at low p T to ∼ 20% at p T = 300 GeV [88]. The p T spectrum of the VBF-produced Higgs boson is therefore reweighted, based on the difference between the Powheg-Box+Pythia calculation and the Hawk [41][42][43] calculation which includes these corrections.
In the case of VBF-produced Higgs boson events in the presence of anomalous couplings in the H V V vertex, the simulated samples are obtained by applying a matrix element (ME) reweighting method to the VBF SM signal sample. The weight is defined as the ratio of the squared ME value for the VBF process associated with a specific amount of CP mixing (measured in terms ofd) to the SM one. The inputs needed for the ME evaluation are the flavour of the incoming partons, the four-momenta and the flavour of the two or three final-state partons and the four-momentum of the Higgs boson. The Bjorken x values of the initial-state partons can be calculated from energy-momentum conservation. The leading-order ME from HAWK [41][42][43] is used for the 2 → 2 + H or 2 → 3 + H process separately. This reweighting procedure is validated against samples generated with MadGraph5_aMC@NLO [40]. As described in Ref. [89], MadGraph5_aMC@NLO can simulate VBF production with anomalous couplings at next-to-leading order. The reweighting procedure proves to be a good approximation to a full next-to-Leading description of the BSM process.  For all samples, a full simulation of the ATLAS detector response [90] using the Geant4 program [91] was performed. In addition, multiple simultaneous minimum-bias interactions are simulated using the AU2 [92] parameter tuning of Pythia8. They are overlaid on the simulated signal and background events according to the luminosity profile of the recorded data. The contributions from these pile-up interactions are simulated both within the same bunch crossing as the hard-scattering process and in neighbouring bunch crossings. Finally, the resulting simulated events are processed through the same reconstruction programs as the data.

Analysis
After data quality requirements, the integrated luminosity of the √ s = 8 TeV dataset used is 20.3 fb −1 . The triggers, event selection, estimation of background contributions and systematic uncertainties closely follow the analysis in Ref.
[20]. In the following a short description of the analysis strategy is given; more details are given in that reference.
Depending on the reconstructed decay modes of the two τ leptons (leptonic or hadronic), events are separated into the dileptonic (τ lep τ lep ) and semileptonic (τ lep τ had ) channels. Following a channel-specific preselection, a VBF region is selected by requiring at least two jets with p Inside the VBF region, boosted decision trees (BDT) 2 are utilised for separating Higgs boson events produced via VBF from the background (including other Higgs boson production modes). The final signal region in each channel is defined by the events with a BDT score value above a threshold of 0.68 for τ lep τ lep and 0.3 for τ lep τ had . The efficiency of this selection, with respect to the full VBF region, is 49% (51%) for the signal and 3.6% (2.1%) for the sum of background processes for the τ lep τ lep (τ lep τ had ) channel. A non-negligible number of events from VBF-produced H → W W events survive the τ lep τ lep selection: they amount to 17% of the overall VBF signal in the signal region. Their contribution is entirely negligible in the τ lep τ had selection. Inside each signal region, the Optimal Observable is then used as the variable with which to probe for CP violation. The BDT score does not affect the mean of the Optimal Observable, as can be seen in Fig. 2.
The modelling of the Optimal Observable distribution for various background processes is validated in dedicated control regions. The top-quark control regions are defined by the same cuts as the corresponding signal region, but inverting the veto on b-tagged jets and not applying the selection on the BDT score (in the τ lep τ had channel a requirement of the transverse mass 3 m T > 40 GeV is also applied). In the τ lep τ lep channel a Z → control region is obtained by requiring two same-flavour opposite-charge leptons, the invariant mass of the two leptons to be 80 < m < 100 GeV, and no BDT score  requirement, but otherwise applying the same requirements as for the signal region. These regions are also used to normalise the respective background estimates using a global fit described in the next section. Finally, an additional region is defined for each channel, called the low-BDT score control region, where a background-dominated region orthogonal to the signal region is selected by requiring the BDT score to be less than 0.05 for τ lep τ lep and less than 0.3 for τ lep τ had . The distribution of the Optimal Observable in these regions is shown in Figs. 3 and 4, demonstrating the good description of the data by the background estimates.
The effect of systematic uncertainties on the yields in signal region and on the shape of the Optimal Observable is eval-uated following the procedures and prescriptions described in Ref. [20]. An additional theoretical uncertainty in the shape of the Optimal Observable is included to account for the signal reweighting procedure described in Sect. 5. This is obtained from the small difference between the Optimal Observable distribution in reweighted samples, compared to samples with anomalous couplings directly generated with MadGraph5_aMC@NLO. While the analysis is statistically limited, the most important systematic uncertainties are found to arise from effects on the jet, hadronically decaying τ and electron energy scales; the most important theoretical uncertainty is due to the description of the underlying event and parton shower in the VBF signal sample.   The best estimate ofd is obtained using a maximumlikelihood fit performed on the Optimal Observable distribution in the signal region for each decay channel simultaneously, with information from different control regions included to constrain background normalisations and nuisance parameters. The normalisation of the VBF H → τ τ and H → W W signal sample is left free in the fit, i.e. this analysis only exploits the shape of the Optimal Observable and does not depend on any possibly model-dependent information about the cross section of CP-mixing scenarios. The relative proportion of the two Higgs boson decay modes is assumed to be as in the SM. All other Higgs boson production modes are treated as background in this study and normalised to their SM expectation, accounting for the corresponding theoretical uncertainties. A binned likelihood function L(x; μ, θ ) is employed, which is a function of the data x, the free-floating signal strength μ, defined as the ratio of the measured cross section times branching ratio to the Standard Model prediction, and further nuisance parameters θ . It relies on an underlying model of signal plus background, and it is defined as the product of Poisson probability terms for each bin in the distribution of the Optimal Observable. A set of signal templates corresponding to different values of the CP-mixing parameterd is created by reweighting the SM VBF H → τ τ and H → W W signal samples, as described in Sect. 5. The likelihood function is then evaluated for eachd hypothesis using the corresponding signal template, while keeping the same background model. The calculation profiles the nuisance parameters to the best-fit valuesθ, including information about systematic uncertainties and normalisation fac-tors, both of which affect the expected numbers of signal and background events.
After constructing the negative log-likelihood (NLL) curve by calculating the NLL value for eachd hypothesis, the approximate central confidence interval at 68% confidence level (CL) is determined from the best estimatord, at which the NLL curve has its minimum value, by reading off the points at which NLL=NLL−NLL min = 0.5. The expected sensitivity is determined using an Asimov dataset, i.e. a pseudo-data distribution equal to the signal-plus-background expectation for given values ofd and the parameters of the fit, in particular the signal strength μ, and not including statistical fluctuations [93].
In both channels, a region of low BDT score is obtained as described in the preceding section. The distribution of the BDT score itself is fitted in this region, which has a much larger number of background events than the signal region, allowing the nuisance parameters to be constrained by the data. This region provides the main constraint on the Z → τ τ normalisation, which is free to float in the fit. The event yields from the top-quark (in τ lep τ lep and τ lep τ had ) and Z → (in τ lep τ lep only) control regions defined in the previous section are also included in the fit, to constrain the respective background normalisations, which are also left free in the fit.
The distributions of the Optimal Observable in each channel are shown in Fig. 5, with the nuisance parameters, background and signal normalisation adjusted by the global fit performed for thed = 0 hypothesis. Table 2 provides the fitted yields of signal and background events, split into the various contributions, in each channel. The number of events observed in data is also provided.

Results
The mean value of the Optimal Observable for the signal is expected to be zero for a CP-even case, while there may be deviations in case of CP-violating effects. A mean value of zero is also expected for the background, as has been demonstrated. Hence, the mean value in data should also be consistent with zero if there are no CP-violating effects within the precision of this measurement. The observed values for the mean value in data inside the signal regions are 0.3 ± 0.5 for τ lep τ lep and −0.3 ± 0.4 for τ lep τ had , fully consistent with zero within statistical uncertainties and thus showing no hint of CP violation. As described in the previous section, the observed limit on CP-odd couplings is estimated using a global maximumlikelihood fit to the Optimal Observable distributions in data. The observed distribution of NLL as a function of the CP-mixing parameterd for the individual channels sepa- The 68% CL interval presented in this work is a factor 10 better than the one obtained in Ref. [15].
As a comparison, the same procedure for extracting the CP-mixing parameterd was applied using the φ sign j j observable, previously proposed for this measurement and defined in Eq. 11, rather than the Optimal Observable. The expected NLL curves for a SM Higgs boson signal from the combination of both channels for the two CP-odd observables are shown in Fig. 7, allowing a direct comparison, and clearly indicate the better sensitivity of the Optimal Observable. The observed NLL curve derived from the φ sign j j distribution is also consistent withd = 0, as shown in Fig. 8, along with the expectation for a signal withd = 0 scaled to the best-fit signal-strength value (μ = 2.02 +0. 87 −0.77 ).

Conclusions
A test of CP invariance in the Higgs boson coupling to vector bosons has been performed using the vector-boson fusion  Fig. 8 Observed (black) and expected (red) NLL for the combination of both channels as a function of thed values defining the underlying signal hypothesis when using the φ sign j j parameter as the final discriminating variable. An Asimov dataset with SM backgrounds plus pure CP-even VBF signal (d = 0), scaled to the best-fit value of the signal strength in the combined fit when using the φ sign j j parameter (μ = 2.02 +0. 87 −0.77 ) was used to calculate the expected values. The markers indicate the points where an evaluation was made -the lines are only meant to guide the eye production mode and the H → τ τ decay. The dataset corresponds to 20.3 fb −1 of √ s = 8 TeV proton-proton collisions recorded by the ATLAS detector at the LHC. Event selection, background estimation and evaluation of systematic uncertainties are all very similar to the ATLAS analysis that provided evidence of the H → τ τ decay. An Optimal Observable is constructed and utilised, and is shown to provide a substantially better sensitivity than the variable traditionally proposed for this kind of study, φ sign j j . No sign of CP violation is observed. Using only the dileptonic and semileptonic H → τ τ channels, and under the assumptioñ d =d B , values ofd less than −0.11 and greater than 0.05 are excluded at 68% CL. This 68% CL interval is a factor of 10 better than the one previously obtained by the ATLAS experiment from Higgs boson decays into vector bosons. In contrast, the present analysis has no sensitivity to constrain a 95% CL interval with the dataset currently available -however larger data samples in the future and consideration of additional Higgs boson decay channels should make this approach highly competitive.