2HDM cross section and branching ratio for CMS

In this page we describe the production of cross sections and branching ratios for Higgs Bosons in the context of Two Higgs Doublet Model (2HDM).

News:

  • 2015-01-07 I have produced a new, extended set of trees: details below. What's new:
    • three new variables:
      • gammah width of 125GeV scalar boson h
      • ggh gluon fusion production cross section for h
      • bbh b-associated production cross section for h
    • More mass points: the full grid 200 to 1200 GeV has been produced for mH and mA.

2HDM models considered

We consider here only TypeI and TypeII 2HDM [cit.needed]

  • TypeI all quarks are coupled to the first Higgs doublet
  • TypeII up-type and down-type quarks are coupled to different doublet
  • TypeIII and TypeIV differs from TypeI and II only for the coupling of the leptons

There are 5 physical scalars: two CP-even (h and H), one CP-odd (A), two charged (H+-) We assume that the light CP-even (h) is the SM-like Higgs bosons with mass mh=125 GeV.

Parameters

  • pp collisions at $\sqrt{s}=8 TeV$

Masses

  • $m_h=125$ GeV
  • $m_H$ ranges from 200 to 1200 GeV
  • $m_A$ ranges from 200 to 1200 GeV
  • $m_H^{\pm} == m_A$
  • $m_{12}^2=m_A^2*\frac{\tan(\beta)}{1+\tan(\beta)^2}$ ($Z_2$ symmetry broken as in MSSM) as suggested by Higgs WG

Angles

  • $\tan(\beta)$ (ratio of doublet vev) range 1 to 10, with step 1
  • $\alpha$: mixing angles between the two doublets
  • $\cos(\beta-\alpha)\in[-1, 1]$ with step 0.1
  • angle convention $0<\beta-\alpha<\pi$

Other

  • $\lambda_6, \lambda_7=0$ (avoid tree-level CP violation)
  • PDF: MSTW2008lo90cl, nlo, nnlo
  • Renormalization/factorization scale:
    • $\mu_F=\mu_R=0.5 m_{\Phi(A/H)}$ for gluon-fusion process
    • $\mu_F=m_{\Phi}, \mu_R=m_{\Phi/4}$.

Tool

Cross section is calculated with =!SuShi+ ver. 1.2.0

Branching Ratio is computed with !2HDMC ver. 1.6.4

Tree structure

Branch Description
Quality
validity as returned by !SuShi
stability as returned by !SuShi
perturbativity as returned by !SuShi
unitarity as returned by !2HDMC
Masses
gammah Total width of the h [GeV] mh=125GeV new
mA Mass of the pseudo-scalar Higgs bosons [GeV]
gammaA Total width of the A [GeV]
mH Mass of the heavy scalar Higgs bosons [GeV]
gammaH Total width of the H [GeV]
mHC Mass of the charged Higgs bosons [GeV]
gammaHC Total width of charged H [GeV]
m12 $m_{12}$ parameter used [GeV]
Angles
tanbeta $\tan(\beta)$ (vev ratios)
alpha higgs doublets mixing angle
sinba $\sin(\beta-\alpha)$ ($0<\beta-\alpha<pi$)
cosba $\cos(\beta-\alpha)$ ($0<\beta-\alpha<pi$)
Cross Sections
ggh inclusive cross section for h production in the gluon-fusion channel [pb] new
bbh likewise for b-associated production [pb] new
ggA inclusive cross section for A production in the gluon-fusion channel [pb]
bbA likewise for b-associated production [pb]
ggH gluon-fusion cross-section for H [pb]
bbH b-assiociated production cross-section for H [pb]
Branching Ratio
BRh BR of light h ($m_h=125$ GeV)
BRh.ss  
BRh.cc  
BRh.bb  
BRh.ee  
BRh.mumu μμ
BRh.tautau ττ
BRh.gamgam γγ
BRh.ZZ  
BRh.WW  
BRh.Zgam
BRh.gluglu gluons
BRA BR of pseudoscalar Higgs
BRA.ss  
BRA.cc  
BRA.bb  
BRA.tt  
BRA.ee  
BRA.mumu μμ
BRA.tautau ττ
BRA.gamgam γγ
BRA.Zgam
BRA.gluglu  
BRA.Zh  
BRH BR of heavy higgs
BRH.ss  
BRH.cc  
BRH.bb  
BRH.tt  
BRH.ee  
BRH.mumu μμ
BRH.tautau ττ
BRH.gamgam γγ
BRH.ZZ  
BRH.WW  
BRH.Zgam
BRH.gluglu  
BRH.hh  
BRH.AA  
BRH.HCHC two charged Higgs
BRH.ZA
BRH.WHC W+charged Higgs
BRHC BR of charged Higgs
BRHC.us  
BRHC.ub  
BRHC.cd  
BRHC.cs  
BRHC.cb  
BRHC.td  
BRHC.ts  
BRHC.tb  
BRHC.enue e
BRHC.munumu μνμ
BRHC.taunutau τντ
BRHC.Wh  

Scenarios

In paper [arXiv:1312.5571v2] three scenarios for 2HDM are proposed, named A, B, and C. The parameters are the following:

Parameters Scenario A Scenario B Scenario C
Type I II II
Mh (GeV) 125 125 125
MH (GeV) 300 300 400
MA (GeV) 330 270 500
MH± (GeV) 230 335 550
$M^2_{12}$ ($GeV^2$) 25600 1798 15800
tanβ 1.5 50 10
sin(β − α) 0.901314 0.999001 0.999
λ6 0 0 0
λ7 0 0 0
α −0.14 0.0247 −0.0549436

The values for x-sections and BR found with the machinery described in this twiki are the following:

variable Scenarios
A B C
quality
validity 1 1 1
stability 1 1 1
perturbativity 1 1 1
unitarity 1 1 1
masses and width [GeV]
mA 330 270 500
gammaA 1.6347411 13.844805 1.3491742
mH 300 300 400
gammaH 1.8080502 15.117569 0.8999441
mHC 230 335 550
gammaHC 1.3918705 10.993239 4.5759229
m12 25600 1798 15800
angles [rad]
tanbeta 1.5 50 10
alpha -0.139999 0.0247053 -0.054943
sinba.sin 0.9013140 0.9990010 0.9990000
cosba.cos 0.4331663 0.0446878 0.0447101
x-sections [pb]
ggA 5.4635095 5.3575649 0.0359630
bbA 0.0021459 28.941694 0.0666051
ggH 0.0982994 2.7602992 0.0384601
bbH 0.0002069 18.394941 0.2003229
Branching fractions
BR h  
ss 0.0002488 0.0002388 0.0001311
cc 0.0305571 0.0202786 0.0520103
bb 0.6808501 0.6790268 0.3529983
ee 5.453e-09 5.524e-09 2.807e-09
μμ 0.0002331 0.0002361 0.0001200
ττ 0.0658571 0.0667158 0.0339068
γγ 0.0011304 0.0013916 0.0037891
ZZ 0.0154554 0.0178500 0.0454556
WW 0.1236669 0.1428274 0.3637140
0.0008748 0.0010108 0.0026303
gluons 0.0811258 0.0704234 0.1452440
BR A  
ss 5.399e-07 0.0003186 0.0002171
cc 7.218e-05 6.511e-09 2.773e-06
bb 0.0015685 0.8975482 0.6117505
tt 0.0140210 1.202e-09 0.1689434
ee 1.537e-11 8.355e-09 6.351e-09
μμ 6.573e-07 0.0003572 0.0002715
ττ 0.0001858 0.1010153 0.0767890
γγ 2.383e-05 1.125e-07 1.761e-06
5.725e-06 5.260e-08 4.851e-07
gluons 0.0088753 0.0005647 0.0008383
Zh 0.5745874 0.0001956 0.0442508
BR H  
ss 3.078e-08 0.0003179 0.0002724
cc 3.786e-06 9.829e-09 1.043e-06
bb 8.493e-05 0.8946859 0.7671372
tt 5.982e-06 1.552e-08 0.0141631
ee 0 8.500e-09 7.670e-09
μμ 3.419e-08 0.0003634 0.0003279
ττ 9.668e-06 0.1027553 0.0927279
γγ 1.642e-05 1.043e-07 1.187e-06
ZZ 5.300e-05 1.031e-07 1.764e-06
WW 0.2606424 0.0003317 0.0167151
0.5872392 0.0007475 0.0358604
gluons 8.383e-05 0.0004476 0.0005202
hh 0.1492024 0.0003493 0.0722715
AA 0 0 0
$H^\pm H^\pm$ 0 0 0
ZA 0 7.890e-07 0
W$H^\pm$ 0.0013291 0 0
BR $H^\pm$  
us 3.078e-08 0.0003179 0.0002724
ub 3.786e-06 9.829e-09 1.043e-06
cd 8.493e-05 0.8946859 0.7671372
cs 5.982e-06 1.552e-08 0.0141631
cb 0 8.500e-09 7.670e-09
td 3.419e-08 0.0003634 0.0003279
ts 9.668e-06 0.1027553 0.0927279
tb 1.642e-05 1.043e-07 1.187e-06
e 5.300e-05 1.031e-07 1.764e-06
μνμ 0.2606424 0.0003317 0.0167151
τντ 0.5872392 0.0007475 0.0358604
Wh 8.383e-05 0.0004476 0.0005202

How to use the Tree

The variables described above are directly available in the tree, so it's possible just draw them, eg as:

Tree2HDM.Draw("ggA+bbA:mA")

Or if you want to get a graph (and not a scatter plot), you can create a TGraph as follows:

TGraph *gr=new TGraph(Tree2HDM.GetSelectedRows(),
Tree2HDM.GetV2(),
Tree2HDM.GetV1())
gr->GetYaxis()->SetTitle("#sigma(gg#rightarrowA+bb#rightarrowbbA) [pb]")
gr->GetXaxis()->SetTitle("m_{A} [GeV]")

gr->Draw("al")

For 2D plots the procedure is similar:

Tree2HDM.Draw("(ggA+bbA)*BRA.Zh*BRh.bb:tanbeta:cosba")
TGraph2D *gr2=new TGraph2D(Tree2HDM.GetSelectedRows(),
           Tree2HDM.GetV2(), Tree2HDM.GetV3(), Tree2HDM.GetV1())
gr2->Draw("lego2")
gr2->GetXaxis()->SetTitle("#tan#beta")
gr2->GetYaxis()->SetTitle("#cos(#beta-#alpha)")
gr2->GetZaxis()->SetTitle("#sigma(pp#rightarrowA)*BR(A#rightarrowZh)*BR(h#rightarrowbb)")
gr2->Draw()

Macros:

  • plot.C: A simple macro to draw a 1D or 2D TGraph
You should edit the macro to load the tree you are interested in.

Usage within root:

.L plot.C
plot2d("(ggH)*BRH.gamgam:cosba:tanbeta", "abs(cosba)<0.6", "COLZ")
plot("(ggH)*BRH.gamgam:cosba","","")

  • plotChain.C: Example macro to load a chain of all trees and plot 1D and 2D

You need to create a list with the trees you want to put in the chain: the idea is to have all the TypeI and TypeII separately.

List with all Trees from afs:

Usage:

.L plotChain.C
plot2d("(bbA+ggA)*BRA.Zh*BRh.bb:cosba:tanbeta","mA==300 && mH==300","lego2")

Trees produced

Trees produced are available in afs at : ~slacapra/public/2HDM/

  • StdGrid contains the standard scan for a range of masses mH and mA

  • 13TeV contains a limited scan for 13 TeV

  • StepTB0.05 contains a detailed scan for mH=mA=300 with step in TanBeta 0.05, required for HIG-14-034
    • Summary_2HDM_tyY_mA300_mH300.root is the complete tree, the individual ones are available as well

  • StepTB0.05_CBA0.02 contains a detailed scan for mH=mA=300 with step in TanBeta 0.05 and step in cos(b-a) 0.02, also required for HIG-14-034
    • Summary_2HDM_tyY_mA300_mH300.root is the complete tree, the individual ones are available as well

Summary_2HDM_ty2_mAXXX_mHXXX.root 
mA==mH==XXX in 200 - 350 with 10 GeV step and 400 - 1200 range with step 100

Distribution of mass points generated with the previous trees.

  • mAVsMH.png:
    mAVsMH_new.png

Distribution of angles (tan(beta) and cos(beta-alpha) ) for each mass points generated.

  • tanBetaVsCosBA.png:
    tanBetaVsCosBA.png

Documentation

Bibliography

  • Harlander et al, Interim recommendations for the evaluation of Higgs production cross sections and branching ratios at the LHC in the Two-Higgs-Doublet Model, http://arxiv.org/abs/1312.5571
  • Branco et al, Theory and phenomenology of two-Higgs-doublet models, http://dx.doi.org/10.1016/j.physrep.2012.02.002
  • Robert V. Harlander, Stefan Liebler, Hendrik Mantler, SusHi: A program for the calculation of Higgs production in gluon fusion and bottom-quark annihilation in the Standard Model and the MSSM Comp. Phys. Commun. 184 (2013) 1605-1617 [arXiv:1212.3249] DOI: 10.1016/j.cpc.2013.02.006
  • D. Eriksson, J. Rathsman, O. Stål, 2HDMC - Two-Higgs-Doublet Model Calculator Comput.Phys.Commun.181:189-205 (2010); Comput.Phys.Commun.181:833-834 (2010) [arXiv:0902.0851]
  • N. Craig, J. Galloway, S. Thomas, "Searching for Signs of the Second Higgs Doublet", http://arxiv.org/abs/1305.2424

-- StefanoLacaprara - 2014-11-12

Topic attachments
I Attachment History Action Size Date Who Comment
Texttxt list_ty1_afs.txt r1 manage 46.5 K 2015-01-07 - 17:07 StefanoLacaprara  
Texttxt list_ty2_afs.txt r1 manage 46.5 K 2015-01-07 - 17:07 StefanoLacaprara  
PNGpng mAVsMH.png r1 manage 8.5 K 2014-11-18 - 11:31 StefanoLacaprara  
PNGpng mAVsMH_new.png r1 manage 27.2 K 2015-01-07 - 16:22 StefanoLacaprara  
C source code filec plot.C r1 manage 1.0 K 2014-11-18 - 11:36 StefanoLacaprara Example macros
C source code filec plotChain.C r2 r1 manage 1.1 K 2015-01-07 - 17:09 StefanoLacaprara Example macros
PNGpng tanBetaVsCosBA.png r1 manage 9.0 K 2014-11-18 - 11:31 StefanoLacaprara  
Edit | Attach | Watch | Print version | History: r12 < r11 < r10 < r9 < r8 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r12 - 2015-09-15 - StefanoLacaprara
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2019 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback