Issue to generate charged Higgs going to cb
Dear aMC@NLO team,
I have trouble generating a process with a charged Higgs which decays into a quark c and b. Mass between 20GeV and 1000GeV are considered and depending on the mass different processes are simulated to take into account the most relevant diagrams:
*******
- For masses below 140GeV, the following process is simulated in aMC@NLO:
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme True
import model sm
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
import model 2HDMtypeII-nobmass
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
generate p p > t t~ [QCD]
output -f
The tops are decayed using Madspin:
set BW_cut %i
set seed %i
define j = g u c d s b u~ c~ d~ s~ b~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
decay t~ > h- b~
decay t > b w+, w+ > wdec wdec
launch
- For masses between 140GeV and 200GeV, the following process is simulated in aMC@NLO:
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme True
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
import model 2HDMtypeII
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
generate p p > w+ b h- b~ / h1 h2 h3 [QCD]
output -f
The W is decayed using Madspin:
set BW_cut %i
set seed %i
define j = g u c d s b u~ c~ d~ s~ b~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
decay w+ > wdec wdec
- For masses over 200GeV, the following process is simulated in aMC@NLO:
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme True
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
import model 2HDMtypeII
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
generate p p > t h- b~ [QCD]
output -f
The top is decayed using Madspin:
set BW_cut %i
set seed %i
define j = g u c d s b u~ c~ d~ s~ b~
define wdec = e+ mu+ ta+ e- mu- ta- ve vm vt ve~ vm~ vt~
decay t > b w+, w+ > wdec wdec
launch
*******
The issue is that for certain masses, the following error is raised:
10:39:48 INFO: Collecting events
10:41:36 quit
10:41:36 INFO:
10:41:37 Py:MadGraphUtils INFO Some errors detected by MadGraphControl - checking for serious errors
10:41:37 Py:MadGraphUtils INFO stty: standard input: Inappropriate ioctl for device
10:41:37 Py:MadGraphUtils INFO stty: standard input: Inappropriate ioctl for device
10:41:37 Py:MadGraphUtils ERROR [1;31mError detected in "launch --name=run_01"
10:41:37 Py:MadGraphUtils ERROR write debug file /eos/home-
10:41:37 Py:MadGraphUtils ERROR If you need help with this issue, please, contact us on https:/
10:41:37 Py:MadGraphUtils ERROR str : An error occurred during event generation. The event file has not been created:
10:41:37 Py:MadGraphUtils ERROR Overwrite the event weights?
10:41:37 Py:MadGraphUtils ERROR give '0' to keep original weights;
10:41:37 Py:MadGraphUtils ERROR give '1' to overwrite the weights to sum to the Xsec;
10:41:37 Py:MadGraphUtils ERROR give '2' to overwrite the weights to average to the Xsec (=default)
10:41:37 Py:MadGraphUtils ERROR give '3' to overwrite the weights to either +/- 1.
10:41:37 Py:MadGraphUtils ERROR 2
10:41:37 Py:MadGraphUtils ERROR step # 0
10:41:37 Py:MadGraphUtils ERROR get_orderstags_
10:41:37 Py:MadGraphUtils ERROR get_orderstags_
10:41:37 Py:MadGraphUtils ERROR found 5 files, bunch number is 1
10:41:37 Py:MadGraphUtils ERROR [0m
For example, for mass 20GeV, the error appears if complex mass scheme (CMS) is set to True. On the other hand, for mass 130GeV the error does not appear neither with CMS neither without it. Similarly, for 400GeV the code works with and without CMS but for mass 1000GeV it works only without it. However, I am starting to suspect that maybe the error is not related to CMS, but it happens randomly, because for 600GeV the code does not work neither without CMS.
*******
I copied the contain of run_01_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.3.1 2021-12-04 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* and *
#* http://
#* *
#******
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#******
launch --name=run_01
Traceback (most recent call last):
File "/cvmfs/
return self.onecmd_
File "/cvmfs/
return func(arg, **opt)
File "/eos/home-
evt_file = self.run(mode, options)
File "/eos/home-
return self.reweight_
File "/eos/home-
raise aMCatNLOError('An error occurred during event generation. ' + \
amcatnlo_
Overwrite the event weights?
give '0' to keep original weights;
give '1' to overwrite the weights to sum to the Xsec;
give '2' to overwrite the weights to average to the Xsec (=default)
give '3' to overwrite the weights to either +/- 1.
2
step # 0
get_orderstags
get_orderstags
found 5 files, bunch number is 1
Value of current Options:
exrootana
cluster_
automatic_
notificat
output_
*******
I copied the run card:
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#******
11000 = nevents ! Number of unweighted events requested
-1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
-1 = nevt_job! Max number of events per job in event generation.
! (-1= no split).
#******
# Output format
#******
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum/bias. Normalization of the weight in the LHEF
#******
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#******
0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
5000 = npoints_FO_grid ! number of points to setup grids
4 = niters_FO_grid ! number of iter. to setup grids
10000 = npoints_FO ! number of points to compute Xsec
6 = niters_FO ! number of iter. to compute Xsec
#******
# Random number seed *
#******
12345 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
#******
1 = lpp1 ! beam 1 type (0 = no PDF)
1 = lpp2 ! beam 2 type (0 = no PDF)
6500.0 = ebeam1 ! beam 1 energy in GeV
6500.0 = ebeam2 ! beam 2 energy in GeV
#******
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#******
lhapdf = pdlabel ! PDF set
260800 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
! numbers for central PDF sets are allowed. Can be a list;
! PDF sets beyond the first are included via reweighting.
#******
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#******
PYTHIA8 = parton_shower
1.0 = shower_scale_factor ! multiply default shower starting
! scale by this factor
#******
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/
# dynamical_
#******
F = fixed_ren_scale ! if .true. use fixed ren scale
F = fixed_fac_scale ! if .true. use fixed fac scale
91.118 = muR_ref_fixed ! fixed ren reference scale
91.118 = muF_ref_fixed ! fixed fact reference scale
2 = dynamical_
! dynamical choices. Can be a list; scale choices beyond the
! first are included via reweighting
0.5 = muR_over_ref ! ratio of current muR over reference muR
1.0 = muF_over_ref ! ratio of current muF over reference muF
#******
# Reweight variables for scale dependence and PDF uncertainty *
#******
1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
False = reweight_scale ! Reweight to get scale variation using the
! rw_rscale and rw_fscale factors. Should be a list of
! booleans of equal length to dynamical_
! specify for which choice to include scale dependence.
False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
! list booleans of equal length to lhaid to specify for
! which PDF set to include the uncertainties.
#******
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#******
True = store_rwgt_info ! Store info for reweighting in LHE file
#******
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#******
0 = ickkw
#******
#
#******
# BW cutoff (M+/-bwcutoff*
# written in the LHE event file *
#******
15.0 = bwcutoff
#******
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - If gamma_is_j, photons are also clustered with jets. *
# Otherwise, they will be treated as tagged particles and photon *
# isolation will be applied. Note that photons in the real emission *
# will always be clustered with QCD partons. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#******
1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
0.7 = jetradius ! The radius parameter for the jet algorithm
10.0 = ptj ! Min jet transverse momentum
-1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
True = gamma_is_j! Wether to cluster photons as jets or not
#******
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#******
0.0 = ptl ! Min lepton transverse momentum
-1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
0.0 = drll ! Min distance between opposite sign lepton pairs
0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#******
# Fermion-photon recombination parameters *
# If Rphreco=0, no recombination is performed *
#******
0.1 = Rphreco ! Minimum fermion-photon distance for recombination
-1.0 = etaphreco ! Maximum abs(pseudo-rap) for photons to be recombined (a value .lt.0 means no cut)
True = lepphreco ! Recombine photons and leptons together
True = quarkphreco ! Recombine photons and quarks together
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# Not applied if gamma_is_j *
# When ptgmin=0, all the other parameters are ignored *
# More specific cuts can be specified in SubProcesses/cuts.f *
#******
20.0 = ptgmin ! Min photon transverse momentum
-1.0 = etagamma ! Max photon abs(pseudo-rap)
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#******
{} = pt_min_pdg ! Min pT for a massive particle
{} = pt_max_pdg ! Max pT for a massive particle
{} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#******
# Use PineAPPL to generate PDF-independent fast-interpolation grid *
# (https:/
#******
False = pineappl ! PineAPPL switch
#******
3.0 = lhe_version
1 = muf1_over_ref
1 = muf2_over_ref
systematics = systematics_program
['--weight_
*******
For masses between 100 and 140GeV, the branching ratio (BR) of top going to the charged Higgs was small and Madgraph was setting it to zero. For this reason, when this decay was asked in Madspin, an error was raised. This was fixed computing the BR of the top going to Wb and BR of top going to charged Higgs and b and passing these parameters by hand to Madgraph.
The decay of charged Higgs was done through the parton shower with Pythia8.
The version of Magraph that was used is 3.3.1.atlas4.
Let me know if you need any additional information to reproduce the error.
Kind regards,
Adrià
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
|
Revision history for this message
|
#1 |
|
Revision history for this message
|
#2 |
|
Revision history for this message
|
#3 |
Can you help with this problem?
Provide an answer of your own, or ask Adrián for more information if necessary.
