Probing heavy charged Higgs bosons through bottom flavored hadrons in the H+→b¯t→B+X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^+\rightarrow \bar{b}t\rightarrow B+X$$\end{document} channel in 2HDM

Observing light or heavy charged Higgs bosons H±\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^\pm $$\end{document}, lighter or heavier than the top quark, would be instant evidence of physics beyond the Standard Model. For this reason, in recent years searches for charged Higgs bosons have been in the center of attention of current colliders such as the CERN Large Hadron Collider (LHC). In spite of all efforts, no signal has been yet observed. Especially, the results of CMS and ATLAS experiments have excluded a large region in the MSSM mH+-tanβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H^+}-\tan \beta $$\end{document} parameter space for mH+=80-160\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H^+}=80{-}160$$\end{document} GeV corresponding to the entire range of tanβ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tan \beta $$\end{document} up to 60. Therefore, it seems that one should concentrate on probing heavy charged Higgs bosons (mH±>mt\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_{H^\pm }>m_t$$\end{document}) so in this context each new probing channel is welcomed. In this work, we intend to present our proposed channel to search for heavy charged Higgses through the study of scaled-energy distribution of bottom-flavored mesons (B) inclusively produced in charged Higgs decay, i.e., H+→tb¯→B+X\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$H^+\rightarrow t\bar{b}\rightarrow B+X$$\end{document}. Our study is carried out within the framework of the generic two Higgs doublet model (2HDM) using the massless scheme where the zero mass parton approximation is adopted for bottom quark.


Introduction
Despite of all successes of the standard model (SM), this model does not represent a theory of everything since there remain many unsolved open questions such as the origin of dark matter, matter-antimatter asymmetry in the universe, the hierarchy problem, etc. To solve these problems many theories have been proposed which are generally qualified as the theories beyond the SM (BSM). Among the most important ones are those based on the supersymmetry. These extended models often contain an extended Higgs sector. As an overview, the minimal extensions known as two-Higgsa e-mail: mmoosavi@yazd.ac.ir (corresponding author) doublet models (2HDMs) [1] include a second complex Higgs doublet which, after spontaneous symmetry breaking, leads to five physical Higgs boson states, i.e., two neutral scalars (h and H , with the assumption m h < m H ), two charged Higgs bosons (H ± ) and one neutral pseudoscalar (A) [2]. Furthermore, after imposing a discrete symmetry that gives natural flavor conservation the 2HDMs can be also classified into four categories; Type I, II, III and IV, according to the couplings of the doublets to the fermions. The minimal supersymmetric standard model (MSSM) [3] is one of the most popular and very well-studied BSM scenarios where one doublet couples to up quarks and the other to down quarks and charged leptons. It should be noted that, the Higgs sector of the MSSM is a Type-II 2HDM which provides elegant solutions to some of the short comings of the SM. It does also predict rich and various phenomenology to be testable in colliders.
Since there is no fundamental charged scalar boson in the SM, then the discovery of a charged scalar boson would clearly represent unambiguous evidence for the presence of new physics beyond the standard model. In this context, searching for the charged Higgs bosons signal is unique and in this work we propose a new channel to search for them at the current and future colliders.
In all classes of 2HDM scenario, the charged Higgs bosons H ± can appear lighter or heavier than the top quark, while the lightest CP-even Higgs boson h can align with the properties of the SM. Therefore, looking for charged Higgs bosons H ± in various decay channels over a wide range of masses is a top priority program in the current LHC experiments and future colliders.
Experimental searches for light charged Higgs bosons (m H ± < m t ) have already been started at the Tevatron. For example, the CMS [4] and the ATLAS [5] collaborations have reported their results of proton-proton collision data recorded at √ s = 8 TeV using the τ + jets channel with a hadronically decaying τ lepton in the final state, i.e., t → bH + (→ τ + ν τ ). Last results on searching for charged Higgs bosons in the H ± → τ ± ν τ decay channel in proton-proton collisions at √ s = 13 TeV is reported by the CMS experiment [6]. According to reported results, the large region in the MSSM m H + − tan β parameter space is excluded for m H + = 80−160 GeV corresponding to the entire range of tan β up to 60, except a hole around m H + ≈ 150−160 GeV for tan β ≈ 10. Here, tan β is the ratio of the vacuum expectation values of the neutral components of the two Higgs doublets. Therefore, it seems that there is no much chance to find the light charged Higgs bosons and colliders should concentrate on probing the heavy charged Higgs bosons (m H ± > m t ).
Heavy charged Higgs bosons are mainly produced directly in association with a top quark (and also a bottom quark) [7,8]. Moreover, charged Higgs bosons can be produced in supersymmetric (SUSY) cascade decays via heavier neutralino and chargino production in squark and gluino decays, see Refs. [9,10]. On the other hand, in many models a heavy charged Higgs boson is predicted to decay predominantly either to a tau and its associated neutrino, or to a top and a bottom quark (H + → tb). However, the channel H + → tb suffers from large multi-jet background, but it dominates in the heavy mass region, see Refs. [11][12][13][14]. Searches for the signature H + → tb have been interpreted by the ATLAS and CMS Collaborations in proton-proton collisions at center-of-mass energies of 8 [12] and 13 TeV [13][14][15] and a small excluded region in the MSSM m H + − tan β parameter space has been presented. For example, the corresponding searches carried out by ATLAS at √ s = 13 TeV and the integrated luminosity L = 13.2 fb −1 have excluded m H + ≈ 300−900 GeV for a very low tan β(≈ 0.5 − 1.7) region [11], where as for high values of tan β > 44(60), m H + ≈ 300(366) GeV have been excluded. Therefore, large regions in the parameter space are still allowed and corresponding searches are in progress.
In the present work, we study the dominant decay mode H + → tb followed by b → B + X , where B is the bottom-flavored hadron and X collectively denotes the unobserved final state particles. Therefore, our proposed channel to search for heavy charged Higgs bosons at colliders is to study the energy distributions of B-hadrons inclusively produced in the decay mode H + → B + X . To this aim, our primary purpose is the evaluation of the next-to-leading order (NLO) QCD corrections to the differential partial decay width d (H + → tb(+g))/dx b , where x b stands for the scaled-energy of bottom quark. This differential width, which is presented for the first time, is needed to obtain the energy spectrum of B-mesons through heavy charged Higgs decays. Also, the hadronization process b → B is described by the nonperturbative fragmentation functions (FFs) which will be introduced in Sect. 3. The differential decay width at the parton level (d /dx b ), the nonperturbative FFs and the factorization theorem, introduced in Sect. 3, allow us to com-pute the desired physical quantity; the energy spectrum of Bhadrons. Beforehand, in Ref. [16] we have studied the energy spectrum of B-mesons produced form direct decay of top quarks in the SM, i.e., t → BW + + X . It would be expected that a comparison between the energy spectrum of B-mesons from charged Higgs decays and those from top decays at SM indicates a signal for new physics beyond the SM. This paper is organized as follows. In Sect. 2, we express our analytical results of the O(α s ) QCD corrections to the Born level rate of H + → tb. We shall apply the massless scheme where the bottom quark mass is ignored but the arbitrary value of charged Higgs mass is retained. In Sect. 3, we give our numerical analysis of inclusive production of B-hadrons from heavy charged Higgs decay considering the factorization theorem and the DGLAP evaluation equations. Sect. 4 is devoted to our summary and conclusions.

Parton level results in the general 2HDM
Assuming m H + > m t , we first study the NLO radiative corrections to the partial decay width in the general 2HDM, where H 1 and H 2 are the doublets whose vacuum expectation values (VEV's), i.e., v 1 and v 2 , give masses to the down and up type quarks, respectively. The squared sum of VEV's is fixed by the Fermi constant However, the ratio of two VEV's is a free parameter and can be characterized by the angle β by introducing tan β = v 2 /v 1 . A linear combination of the charged components of doublets H 1 and H 2 does also give the observable charged Higgs H ± , i.e., H ± = H ± 2 cos β − H ± 1 sin β. In a general 2HDM, tree-level flavor-changing neutral currents (FCNC) can be avoided if one does not couple the same Higgs doublet to up-and down-type quarks simultaneously. Therefore, for our purpose we need the specific models which naturally stop these problems by restricting the Higgs coupling. In this context, there are two possibilities (which are also called two models) for the two Higgs doublets to couple to the quarks.
In the first possibility (or model I), the Higgs doublet H 1 couples to all bosons and another doublet H 2 couples to all quarks in the same manner as in the SM. In this model, the Yukawa couplings between the top-and the bottomquark and the charged Higgs are given by the following Lagrangian [17] where, g 2 W = 4 √ 2m 2 W G F and the CKM matrix element is labeled by V tb .
In the second possibility (model II), the doublet H 1 couples only to the right chiral down-type quarks while the H 2 couples only to the right chiral up-type quarks. In this model, the charged Higgs boson couplings to fermions are given by the following Lagrangian These two models are also known as Type-I and Type-II 2HDM scenarios and, as mentioned in the Introduction, the MSSM [18][19][20] is a special case of a Type-II 2HDM.
For the process (1), considering the interaction Lagrangians (2) and (3) the current density is expressed as J μ ∝ ψ b (a + bγ 5 )ψ t so that the coupling factors in two models are given by and In next section, we describe the technical detail of our calculation for the O(α s ) radiative corrections to the tree-level decay rate of H + → tb using dimensional regularization to regularize all divergences.

Born decay width of H + → tb
The decay process (1) is analyzed in the rest frame of the charged Higgs boson. It is straightforward to calculate the Born term contribution to the partial decay rate of the process (1) in the 2HDM. According to the given Lagrangian in Eqs. (2) and (3), the coupling of the charged-Higgs to the fermions (top and bottom quark in (1)) can either be expressed as a superposition of scalar and pseudoscalar coupling factors or as a combination of right-and left-chiral coupling factors [17]. Therefore, the lowest order decay amplitude is of the form where, g t = a + b and g b = a − b. Therefore, the tree-level decay width reads where λ(x, y, z) = (x − y − z) 2 − 4yz is the Källén function and N c = 3 is a color factor. Here, for simplicity, we have defined: This result is in complete agreement with the one presented in Ref. [21]. In the limit of vanishing bottom quark mass, the tree-level decay width is of the form where, in both models I and II one has Since m b << m t , finite-m b corrections are expected to be negligible in the case at hand. This expectation has been actually confirmed in Ref. [16], by a comparative analysis of the partial width of the decay t → bW + in the general-mass variable-flavor-number scheme (GM-VFNS), where bottom-quark mass is preserved, and the zero-mass variable-flavor-number scheme (ZM-VFNS), where bottom is included among the massless quark flavors. Then, throughout this work we apply the ZM-VFNS or massless scheme. In next section, we compute the O(α s ) QCD corrections to the Born-level decay rate of H + → tb and present, for the first time, the analytical parton-level expressions for d (H + → B + X )/dx B at NLO in the ZM-VFNS. To this aim, we calculate the quantity d b /dx b where, is the scaled-energy of b-quark. It ranges as 0 ≤ x b ≤ 1.

O(α s ) virtual corrections
The QCD virtual one-loop corrections to the process H + → tb contain both infrared (IR) and ultraviolet (UV) divergences where the UV-divergences appear when the integration region of the internal momentum of the virtual gluon goes to infinity and the IR-divergences arise from the softgluon singularities. In this work, we adopt the "on-shell" mass renormalization scheme and apply the dimensional regularization scheme to regularize all divergences. Through this scheme, all singularities are regularized in D = 4 − 2 dimensions to become single poles in . Considering the twobody phase space for the virtual corrections the contribution of virtual radiations into the differential decay width reads where, Here, M 0 is the Born term amplitude (6) and the renormalized amplitude of the virtual corrections is given by where ct represents the counterterm and l arises from the one-loop vertex correction [22]. Following Refs. [23,24], the counterterm of the vertex includes the wave-function renormalizations of quarks as well as the top quark mass renormalization Since, we are working in the ZM-VFN scheme where m b = 0 is assumed, then the b-quark mass counterterm is δm b = 0. The wave function and the mass renormalization constants are given by [25] where, γ E = 0.5772 · · · is the Euler constant, C F = (N 2 c − 1)/(2N c ) = 4/3 for N c = 3 quark colors, and μ F is the factorization scale which is arbitrarily set as μ F = m H in our work. Conventionally, I R and U V represent the infrared and the ultraviolet divergences, respectively.
The real part of the vertex correction is given by where, B 0 and C 0 are the Passarino-Veltman 2-point and 3point integrals [26]. By summing all virtual corrections up, the UV-singularities are canceled so that the virtual differential decay rate is ultraviolet finite. But, the IR-divergences are remaining which are now labeled by . Eventually, the virtual one-loop contributions read where, Li 2 (y) = − y 0 (dt/t) ln(1−t) is the Spence function and

Real gluon corrections (Bremsstrahlung)
To obtain the infrared-finite physical results for d b /dx b one must include the contributions of real gluons emission. Considering two Feynman graphs including the real gluon emissions from the top and bottom quark, the O(α s ) real gluon emission (tree-graph) amplitude reads where ( p g , r ) refers to the polarization vector of the emitted real gluon with the spin r . The first and second expressions in the curly brackets are related to the real gluon emissions from the top and bottom quarks, respectively. In order to regulate the IR-divergences which arise from the soft and collinear real-gluon emissions, as before, we apply dimensional regularization scheme. According to this scheme, the real differential decay rate for the process H + → tbg is given by where, μ F is an arbitrary reference mass and the phase space element d R 3 is defined as To evaluate the differential decay rate d real b /dx b , we fix the momentum of bottom quark in Eq. (18) and integrate over the gluon energy which ranges as Note that, when we integrate over the phase space of the real gluon radiation, terms of the form (1−x b ) −1−2 appear which are due to the radiation of soft gluon, i.e., E g → 0 ≡ x b → 1. Thus, we employ the following prescription introduced in Ref. [27] ( where the plus distributions are defined as

Analytical results for d /dx i at parton level
The O(α s ) corrections to the differential decay rate of H + → tb is obtained by summing the Born, the virtual and the real gluon contributions. It reads where, by defining S = (1 − y)/2 (with y = m 2 t /m 2 H ) one has Our result of differential decay rate, which is presented for the first time, after integration over x b (0 ≤ x b ≤ 1) is in complete agreement with the result presented in [21]. Note that, our main purpose is to evaluate the energy distribution of B-hadrons produced in heavy charged Higgs boson decay: H + → tb(+g) → B + X , where B-hadrons can be produced from the fragmentation of b-quark as well as the emitted real gluons. Therefore, in order to obtain the most accurate energy spectrum of produced B-hadrons we have to consider the contribution of gluon fragmentation as well. It should be noted that, the gluon splitting contribution is important at the low energy of the observed B-hadron so this contribution decreases the size of decay rate at the threshold, see Refs. [28,29]. With this explanation, we also need to compute the NLO differential decay rate d nlo g /dx g , where y)) is the scaled-energy of emitted real gluon, as in (10). Ignoring the details of calculation, this differential decay rate is given by where, In Eqs. (23) and (25), the terms T 1 and T 2 are free of all IR-divergences. In order to subtract the singularities remaining in the differential decay widths, we employ the modified minimal-subtraction (M S) scheme, where the singularities are absorbed into the bare fragmentation functions (FFs). This renormalizes the FFs, endowing them with μ F dependence, and creates in the differential decay widths the finite terms of the form (α s /π ) ln(m 2 H /μ 2 F ) which are rendered perturbatively small by choosing μ F = O(m H ). Following the M S scheme, in order to have the finite coefficient functions we have to subtract from Eqs. (23) and (25), the O(α s ) term multiplying the characteristic M S constant, i.e., −1/ + γ E − ln 4π [27].

Numerical results
In this work, using the ZM-VFNS we study the decay process followed byb/g → B + X . In this process, top quark dominantly decays as: t → bW + → bl + ν l . In the narrow-width approximation (NWA), where we set p 2 t = m 2 t and p 2 W + = m 2 W + and ignore small terms of order O( 2 i /m 2 i )(i = t, W + ), the total decay rate reads where, for the branching ratios one has B(t → bW + ) = 96.2% and B(W + → l + ν l ) = 10.86% [30]. More details about the NWA can be found in Ref. [31].
Having the differential decay widths for the process (27), i.e., Eqs. (23) and (25), we are now in a situation to make our phenomenological predictions for the scaled-energy (x B ) distribution of B-hadrons inclusively produced in the decay of heavy charged Higgs bosons. To present our results for the x B -distribution, we consider the differential distribution d nlo /dx B of the partial width of the decay H (1 − y)) is the scaled-energy of Bhadrons in the charged Higgs rest frame. The x B -variable is defined as x b in (10).
Our tool to compute the scaled energy distribution of Bhadrons is the factorization theorem of QCD-improved parton model [32]. According to this theorem [33], the energy distribution of B-hadrons can be expressed as the convolution of the parton-level spectrum d a /dx a (a = b, g) with the nonperturbative FFs of a → B, describing the hadronization process of a → B. The a → B FFs are labeled by D B a (z, μ F ), where μ F is the factorization scale and z = E B /E a is the fragmentation variable which indicates the energy fraction of parent parton carried by the pro-duced hadron. The factorization theorem is expressed as where, μ R and μ F are the renormalization and factorization scales, respectively. The scale μ R is related to the renormalization of the QCD coupling constant. In this paper, we use the convention μ R = μ F = m H , a choice often made. Several searches for the signature H + → tb in the context of 2HDMs have been done by the ATLAS and CMS Collaborations in proton-proton collisions at center-of-mass energies of 8 and 13 TeV [12][13][14]. For example, in Ref. Although, a definitive search over the m H + − tan β plane is a program that still has to be carried out and this belongs to the LHC experiments and future colliders.
In this work, for our numerical analysis we restrict ourselves to the allowed regions of the m H + − tan β parameter space evaluated by the CMS experiments, see Fig.7 in Ref. [13]. Moreover, from Ref. [34] we adopt other input parameters as G F = 1.16637 × 10 −5 GeV −2 and m t = 172.98 GeV. We will also evaluate the QCD coupling constant α s at NLO in the MS scheme through the following relation where, is the QCD scale parameter. Also, b 0 and b 1 are given by where, n f is the number of active quark flavors. In this work, we adopt , it is seen that the QCD corrections decrease the charged Higgs boson decay width and their amounts depend on the charged Higgs mass. Note that, for the total decay rate of process H + →bt (→ bW + (→ l + ν l )) the above results should be multiplied by B(t → bW + ) and B(W + → l + ν l ), see Eq. (28). Now, we go back to our main aim: the evaluation of energy distribution of B-hadrons in heavy charged Higgs decays. For this purpose, we use the factorization relation (29) where to describe the splitting (b, g) → B, from Ref. [35] we employ the realistic nonperturbative B-hadron FFs determined at NLO in the ZM-VFN scheme. These FFs have been determined through a global fit to electronpositron annihilation data presented by ALEPH [36] and OPAL [37] at CERN LEP1 and by SLD [38] at SLAC SLC. According to the approach used in [35], the power ansatz splitting where the free parameters have been determined at the initial scale μ ini F = 4.5 GeV. The fit yielded N = 2575.014, α = 15.424, and β = 2.394. The gluon FF is assumed to be zero at the initial scale μ ini F and generated via the DGLAP evolution equations [39,40].
In Fig. 1  corrections at the parton level, we evaluated the LO distribution using the same NLO b → B FF. Our results show that the NLO corrections lead to a significant enhancement of the partial decay width in the peak region and above, while these corrections decrease the size of partial decay rate in the lower-x B range. It should be noted that, the contribution of gluon splitting is appreciable only in the low-x B region. For higher values of x B , the contribution of b-quark fragmentation dominates, as expected [16]. In Fig. 2, the dependence of x B spectrum on tan β is studied, taking m H + = 200 GeV. As is seen, when the value of tan β increases the decay rate decreases, because the Born rate 0 (8) is proportional to cot 2 β.
In Fig. 3, by fixing tan β = 2 we have investigated the dependence of x B spectrum on the charged Higgs mass taking m H + = 200 (solid line), m H + = 400 GeV (dashed line) and m H + = 600 GeV (dot-dashed line). This figure shows that, if m H + increases the size of partial decay width increases as well. Nevertheless, the peak position of x B -distribution is approximately independent of the charged Higgs mass.

Conclusions
The SM of particle physics predicts one neutral Higgs boson, whereas the Minimal Supersymmetric requires five Higgs particles, three neutral bosons and two charged bosons. The discovery of charged Higgs bosons would be proof of new physics beyond the SM. For this reason, searches for the charged Higgs bosons are strongly motivated so that in recent years it has been a goal of many high energy colliders such as the CERN LHC. Searches for light charged Higgs bosons (particles lighter than the top quark) has been inconclusive and no evidence has been yet found. In this regard, the results reported by the CMS and ATLAS Collaborations show the large excluded region in the MSSM m H + − tan β parameter space. Therefore, it sounds that most efforts should be concentrated on probing heavy charged Higgs bosons (heavier than the top quark). These scalar bosons are predicted to decay predominantly either to a tau and its associated neutrino (τν τ ), or to a top and a bottom quark (tb). In spite of the fact that the decay channel H + → tb suffers from large multi-jet background, but it dominates in the heavy mass region.
In this work, we studied the dominant decay channel H + → tb(+g) followed by the hadronization process (b, g) → B. At colliders, the bottom-flavored hadrons could be identified by a displaced decay vertex associated whit charged lepton tracks. On other words, B-hadrons decay to the J/ψ followed by the J/ψ → μ + μ − decays, see Ref. [41]. Then a muon in jet is associated to the b-flavored hadron. Furthermore, one can also explore an other way to associated the J/ψ with the corresponding isolated leptonby measuring the jet charge of identified b and not requiring the tagging muon. Therefore, at the LHC and future colliders the decay channel H + → B + X is proposed to search for the heavy charged Higgs bosons and evaluating the distribution in the scaled-energy (x B ) of B-mesons would be of particular interest. This distribution is studied by evaluating the quantity d /dx B . To present our phenomenological prediction of the x B -distribution, we first calculated an analytic expression for the NLO radiative corrections to the differential decay width d (H + → tb)/dx a (a = b, g) and then employed the nonperturbative (b, g) → B FFs, relying on their universality and scaling violations. Our results have been presented in the ZM-VFN scheme where the b-quark mass is ignored from the beginning. In this scheme, results are the same in both the type-I and II 2HDM scenarios.
Our analysis is expected to make a contribution to the LHC searches for charged Higgs bosons. In fact, a comparison between the energy spectrum of B-mesons produced from charged Higgs decays at 2HDM and those from top decays at SM (t → B + X ) would indicate a signal for new physics beyond the SM.
Our analysis can be also extended to the production of hadron species other than the B-hadron, such as pions, kaons and protons, etc. This would be possible by using the nonperturbative (b, g) → π/K /P/D + FFs presented in Refs. [42][43][44][45].
Data Availability Statement This manuscript has no associated data or the data will not be deposited. [Authors' comment: All Data are reported by the CMS and ATLAS collaborations and can be found in the website of CERN LHC.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .