Skip to content
This repository was archived by the owner on Jan 30, 2026. It is now read-only.

Commit 1428543

Browse files
bbilinbbilin
authored
Merge pull request #3691 from pmastrap/scriptsAndTemplate
Scripts and template to produce H->ZA->lltautau
2 parents 9eedbf5 + ed63e38 commit 1428543

17 files changed

+1927
-0
lines changed

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/README.md‎

Lines changed: 40 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,40 @@
1+
# Full RunII H/A->ZA/H->lltautau
2+
## GridPacks production for full run2 ULegacy
3+
This README is mutuated from Khawla's HToZATo2L2B repo. Please refer to that for further explanation.
4+
5+
## Setup Your Enviroment and Prepare Template cards:
6+
This needs 2HDMC, which is a general-purpose calculator for the two-Higgs doublet model.
7+
It allows parametrization of the Higgs potential in many different ways, convenient specification of generic Yukawa sectors, the evaluation of decay widths (including higher-order QCD corrections), theoretical constraints and much more.
8+
9+
You can install Calculators42HDM in a CMSSW release (recommended)
10+
or a conda environment (which requires a few changes to the script), see the [installation instructions](https://github.com/kjaffel/Calculators42HDM/blob/master/README.md).
11+
## GridPacks Preparation:
12+
### Example of Template Cards:
13+
- For each process we have a dir template, HToZA or AToZH!
14+
```bash
15+
# gg Fusion; LO Loop Induced 4F-scheme
16+
# H -> ZA (mH > mA)
17+
run2Template_cards/template_HToZATo2L2Tau_X_X_X_ggH_TuneCP5_13TeV_pythia8
18+
# A -> ZH (mA > mH)
19+
run2Template_cards/template_AToZHTo2L2Tau_X_X_X_ggH_TuneCP5_13TeV_pythia8
20+
```
21+
``X_X_X`` : MH_MA_tanbeta parameters
22+
### How to Run:
23+
```python
24+
./prepare_MG5_cards.py --process ggH --test --templates run2Template_cards/ --mode H
25+
./prepare_MG5_cards.py --process ggH --test --templates run2Template_cards/ --mode A
26+
```
27+
```bash
28+
cd ./scripts/
29+
# compute the decay BR and width
30+
{}/MG5_aMC_vX_X_X/bin/mg5_aMC run_madwidths.sh
31+
./write_paramcards.sh
32+
33+
```
34+
## GridPacks Generation:
35+
Inside the ``scripts`` directory a simple shell script is generated to produce all the gridpacks for each process.
36+
It should be run from the genproductions/bin/MadGraph5_aMCatNLO directory in a clean environment.
37+
```bash
38+
./prepare_example_lo_htoza_gridpacks.sh
39+
./prepare_example_lo_atozh_gridpacks.sh
40+
```

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/prepare_MG5_cards.py‎

Lines changed: 671 additions & 0 deletions
Large diffs are not rendered by default.

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/run2Template_cards/template_AToZHTo2L2Tau_200_50_1_ggH_TuneCP5_13TeV_pythia8/AToZHTo2L2B_200_50_1_ggH_TuneCP5_13TeV_pythia8_extramodels.dat‎

Lines changed: 2 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,2 @@
1+
# http://feynrules.irmp.ucl.ac.be/attachment/wiki/2HDM/2HDMtII_NLO.tar.gz
2+
2HDMtII_NLO.tar.gz

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/run2Template_cards/template_AToZHTo2L2Tau_200_50_1_ggH_TuneCP5_13TeV_pythia8/AToZHTo2L2B_200_50_1_ggH_TuneCP5_13TeV_pythia8_madspin_card.dat‎

Lines changed: 6 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,6 @@
1+
set ms_dir ./madspingrid
2+
set spinmode none # Use one of the madspin special mode
3+
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
4+
# specify the decay for the final state particles
5+
decay h2 > ta+ ta-
6+
launch

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/run2Template_cards/template_AToZHTo2L2Tau_200_50_1_ggH_TuneCP5_13TeV_pythia8/AToZHTo2L2B_200_50_1_ggH_TuneCP5_13TeV_pythia8_param_card.dat‎

Lines changed: 152 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,152 @@
1+
######################################################################
2+
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
3+
######################################################################
4+
## ##
5+
## Width set on Auto will be computed following the information ##
6+
## present in the decay.py files of the model. ##
7+
## See arXiv:1402.1178 for more details. ##
8+
## ##
9+
######################################################################
10+
11+
###################################
12+
## INFORMATION FOR FRBLOCK
13+
###################################
14+
Block frblock
15+
1 1.537480e+00 # tanbeta
16+
2 1.000000e-01 # sinbma
17+
18+
###################################
19+
## INFORMATION FOR HIGGS
20+
###################################
21+
Block higgs
22+
1 5.000000e-01 # l2
23+
2 1.000000e+00 # l3
24+
3 1.000000e-01 # lR7
25+
26+
###################################
27+
## INFORMATION FOR LOOP
28+
###################################
29+
Block loop
30+
1 9.118800e+01 # MU_R
31+
32+
###################################
33+
## INFORMATION FOR MASS
34+
###################################
35+
Block mass
36+
5 4.700000e+00 # MB
37+
6 1.720000e+02 # MT
38+
23 9.118760e+01 # MZ
39+
25 1.200000e+02 # mh1
40+
35 1.300000e+02 # mh2
41+
36 1.400000e+02 # mh3
42+
37 1.500000e+02 # mhc
43+
## Dependent parameters, given by model restrictions.
44+
## Those values should be edited following the
45+
## analytical expression. MG5 ignores those values
46+
## but they are important for interfacing the output of MG5
47+
## to external program such as Pythia.
48+
1 0.000000e+00 # d : 0.0
49+
2 0.000000e+00 # u : 0.0
50+
3 0.000000e+00 # s : 0.0
51+
4 0.000000e+00 # c : 0.0
52+
11 0.000000e+00 # e- : 0.0
53+
12 0.000000e+00 # ve : 0.0
54+
13 0.000000e+00 # mu- : 0.0
55+
14 0.000000e+00 # vm : 0.0
56+
15 0.000000e+00 # ta- : 0.0
57+
16 0.000000e+00 # vt : 0.0
58+
21 0.000000e+00 # g : 0.0
59+
22 0.000000e+00 # a : 0.0
60+
24 7.982466e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
61+
9000002 9.118760e+01 # ghz : MZ
62+
9000003 7.982466e+01 # ghwp : MW
63+
9000004 7.982466e+01 # ghwm : MW
64+
65+
###################################
66+
## INFORMATION FOR SMINPUTS
67+
###################################
68+
Block sminputs
69+
1 1.279000e+02 # aEWM1
70+
2 1.166390e-05 # Gf
71+
3 1.180000e-01 # aS
72+
73+
###################################
74+
## INFORMATION FOR YUKAWA
75+
###################################
76+
Block yukawa
77+
5 4.700000e+00 # ymb
78+
6 1.720000e+02 # ymt
79+
15 1.777000e+00 # ymtau
80+
81+
###################################
82+
## INFORMATION FOR DECAY
83+
###################################
84+
DECAY 6 1.508336e+00 # WT
85+
DECAY 23 2.495200e+00 # WZ
86+
DECAY 24 2.085000e+00 # WW
87+
DECAY 25 1.000000e+00 # Wh1
88+
DECAY 35 1.000000e+00 # Wh2
89+
DECAY 36 1.000000e+00 # Wh3
90+
DECAY 37 1.000000e+00 # whc
91+
## Dependent parameters, given by model restrictions.
92+
## Those values should be edited following the
93+
## analytical expression. MG5 ignores those values
94+
## but they are important for interfacing the output of MG5
95+
## to external program such as Pythia.
96+
DECAY 1 0.000000e+00 # d : 0.0
97+
DECAY 2 0.000000e+00 # u : 0.0
98+
DECAY 3 0.000000e+00 # s : 0.0
99+
DECAY 4 0.000000e+00 # c : 0.0
100+
DECAY 5 0.000000e+00 # b : 0.0
101+
DECAY 11 0.000000e+00 # e- : 0.0
102+
DECAY 12 0.000000e+00 # ve : 0.0
103+
DECAY 13 0.000000e+00 # mu- : 0.0
104+
DECAY 14 0.000000e+00 # vm : 0.0
105+
DECAY 15 0.000000e+00 # ta- : 0.0
106+
DECAY 16 0.000000e+00 # vt : 0.0
107+
DECAY 21 0.000000e+00 # g : 0.0
108+
DECAY 22 0.000000e+00 # a : 0.0
109+
#===========================================================
110+
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
111+
#===========================================================
112+
113+
Block QNUMBERS 9000001 # gha
114+
1 0 # 3 times electric charge
115+
2 1 # number of spin states (2S+1)
116+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
117+
4 1 # Particle/Antiparticle distinction (0=own anti)
118+
Block QNUMBERS 9000002 # ghz
119+
1 0 # 3 times electric charge
120+
2 1 # number of spin states (2S+1)
121+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
122+
4 1 # Particle/Antiparticle distinction (0=own anti)
123+
Block QNUMBERS 9000003 # ghwp
124+
1 3 # 3 times electric charge
125+
2 1 # number of spin states (2S+1)
126+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
127+
4 1 # Particle/Antiparticle distinction (0=own anti)
128+
Block QNUMBERS 9000004 # ghwm
129+
1 -3 # 3 times electric charge
130+
2 1 # number of spin states (2S+1)
131+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
132+
4 1 # Particle/Antiparticle distinction (0=own anti)
133+
Block QNUMBERS 82 # ghg
134+
1 0 # 3 times electric charge
135+
2 1 # number of spin states (2S+1)
136+
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
137+
4 1 # Particle/Antiparticle distinction (0=own anti)
138+
Block QNUMBERS 37 # h+
139+
1 3 # 3 times electric charge
140+
2 1 # number of spin states (2S+1)
141+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
142+
4 1 # Particle/Antiparticle distinction (0=own anti)
143+
Block QNUMBERS 35 # h2
144+
1 0 # 3 times electric charge
145+
2 1 # number of spin states (2S+1)
146+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
147+
4 0 # Particle/Antiparticle distinction (0=own anti)
148+
Block QNUMBERS 36 # h3
149+
1 0 # 3 times electric charge
150+
2 1 # number of spin states (2S+1)
151+
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
152+
4 0 # Particle/Antiparticle distinction (0=own anti)

‎bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/HToZATo2L2Tau/run2Template_cards/template_AToZHTo2L2Tau_200_50_1_ggH_TuneCP5_13TeV_pythia8/AToZHTo2L2B_200_50_1_ggH_TuneCP5_13TeV_pythia8_proc_card.dat‎

Lines changed: 4 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,4 @@
1+
import model 2HDMtII_NLO
2+
#gg fusion -loop induced
3+
generate p p > h3 > h2 l+ l- $$ a [QCD] @0
4+
output AToZHTo2L2B_200_50_1_ggH -nojpeg

0 commit comments

Comments
 (0)