Issue to generate charged Higgs going to cb

Asked by Adrián

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:

**************************************Process**************************************
 - For masses below 140GeV, the following process is simulated in aMC@NLO:
    set group_subprocesses Auto
    set ignore_six_quark_processes False
    set loop_optimized_output True
    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_six_quark_processes False
    set loop_optimized_output True
    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_six_quark_processes False
    set loop_optimized_output True
    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

***********************************Error Message***********************************
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-a/adberroc/Analysis_Hplus/plus_cb_the_revenge/MCGeneration/run/Samples/Run2__r23.6.10/Hminus/20__Recommended5FSPDF/NLO_ComplexMassScheme/tt/PROCNLO_2HDMtypeII-nobmass_0/run_01_tag_1_debug.log
10:41:37 Py:MadGraphUtils ERROR If you need help with this issue, please, contact us on https://answers.launchpad.net/mg5amcnlo
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_glob_infos: n_orderstags= 2
10:41:37 Py:MadGraphUtils ERROR get_orderstags_glob_infos: orderstags_glob 4 6
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.

********************************run_01_tag_1_debug********************************
I copied the contain of run_01_tag_1_debug.log file for the failed 20GeV process, the other failed masses have similar output messages except for the numbers between “step # 0” and “Value of current Options:”:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.3.1 2021-12-04 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch --name=run_01
Traceback (most recent call last):
  File "/cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/LCG_102b_ATLAS_12/MCGenerators/madgraph5amc/3.3.1.atlas4/x86_64-centos7-gcc11-opt/madgraph/interface/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/LCG_102b_ATLAS_12/MCGenerators/madgraph5amc/3.3.1.atlas4/x86_64-centos7-gcc11-opt/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/eos/home-a/adberroc/Analysis_Hplus/plus_cb_the_revenge/MCGeneration/run/Samples/Run2__r23.6.10/Hminus/20__Recommended5FSPDF/NLO_ComplexMassScheme/tt/PROCNLO_2HDMtypeII-nobmass_0/bin/internal/amcatnlo_run_interface.py", line 1784, in do_launch
    evt_file = self.run(mode, options)
  File "/eos/home-a/adberroc/Analysis_Hplus/plus_cb_the_revenge/MCGeneration/run/Samples/Run2__r23.6.10/Hminus/20__Recommended5FSPDF/NLO_ComplexMassScheme/tt/PROCNLO_2HDMtypeII-nobmass_0/bin/internal/amcatnlo_run_interface.py", line 2029, in run
    return self.reweight_and_collect_events(options, mode, nevents, event_norm)
  File "/eos/home-a/adberroc/Analysis_Hplus/plus_cb_the_revenge/MCGeneration/run/Samples/Run2__r23.6.10/Hminus/20__Recommended5FSPDF/NLO_ComplexMassScheme/tt/PROCNLO_2HDMtypeII-nobmass_0/bin/internal/amcatnlo_run_interface.py", line 3785, in reweight_and_collect_events
    raise aMCatNLOError('An error occurred during event generation. ' + \
amcatnlo_run_interface.aMCatNLOError: An error occurred during event generation. The event file has not been created:
  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_glob_infos: n_orderstags= 2
 get_orderstags_glob_infos: orderstags_glob 4 6
 found 5 files, bunch number is 1

Value of current Options:
             pythia8_path : None
                hwpp_path : None
              thepeg_path : None
               hepmc_path : None
         madanalysis_path : None
        madanalysis5_path : None
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
      exrootanalysis_path : None
             syscalc_path : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/MCGenerators/syscalc/1.1.7-6c332/x86_64-centos7-gcc11-opt
                   lhapdf : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/MCGenerators/lhapdf/6.5.3-d2254/x86_64-centos7-gcc11-opt/bin/lhapdf-config
               lhapdf_py2 : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/MCGenerators/lhapdf/6.5.3-d2254/x86_64-centos7-gcc11-opt/bin/lhapdf-config
               lhapdf_py3 : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/MCGenerators/lhapdf/6.5.3-d2254/x86_64-centos7-gcc11-opt/bin/lhapdf-config
                  timeout : 60
            f2py_compiler : None
        f2py_compiler_py2 : None
        f2py_compiler_py3 : None
              web_browser : None
               eps_viewer : None
              text_editor : None
         fortran_compiler : None
             cpp_compiler : None
              auto_update : 0
             cluster_type : condor
    cluster_status_update : (600, 30)
         cluster_nb_retry : 1
       cluster_local_path : None
       cluster_retry_wait : 300
             stdout_level : None
   automatic_html_opening : False
      notification_center : True
                 run_mode : 0
            cluster_queue : None
             cluster_time : None
             cluster_size : 100
           cluster_memory : None
                  nb_core : 1
        cluster_temp_path : None
                  fastjet : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/LCG_102b_ATLAS_12/fastjet/3.4.0p1/x86_64-centos7-gcc11-opt/bin/fastjet-config
                    ninja : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/LCG_102b_ATLAS_12/MCGenerators//gosam_contrib/2.0-779ba/x86_64-centos7-gcc11-opt/lib
                  collier : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/LCG_102b_ATLAS_12/MCGenerators//collier/1.2.5-e3f73/x86_64-centos7-gcc11-opt/lib
      output_dependencies : internal
                 mg5_path : /cvmfs/atlas.cern.ch/repo/sw/software/23.6/sw/lcg/releases/MCGenerators/madgraph5amc/3.3.1.atlas4-1a556/x86_64-centos7-gcc11-opt

*************************************Run Card*************************************
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/conventions: *
# *
# 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=default))
#***********************************************************************
# 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/set_scales.f or via *
# dynamical_scale_choice option) *
#***********************************************************************
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_scale_choice ! Choose one (or more) of the predefined
! 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_scale_choice to
! 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://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 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*Gamma). Determines which resonances are *
# 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://zenodo.org/record/3992765#.X2EWy5MzbVo) *
#***********************************************************************
False = pineappl ! PineAPPL switch
#***********************************************************************
3.0 = lhe_version
1 = muf1_over_ref
1 = muf2_over_ref
systematics = systematics_program
['--weight_info=MUR%(mur).1f_MUF%(muf).1f_PDF%(pdf)i','--remove_wgts=".*MUR0.5_MUF2.0.*|.*MUR2.0_MUF0.5.*"','--pdf=NNPDF30_nlo_as_0118_mc,PDF4LHC21_40_pdfas,NNPDF30_nlo_as_0119@0,NNPDF30_nlo_as_0117@0,NNPDF30_nnlo_as_0118_hessian@0,MSHT20nnlo_as118@0,MSHT20nlo_as118@0,CT18NNLO@0,CT18NLO@0,NNPDF31_nnlo_as_0118_hessian@0,NNPDF31_nlo_as_0118_hessian@0,NNPDF40_nnlo_as_01180_hessian@0,NNPDF40_nlo_as_01180@0,CT18ANNLO@0,CT18XNNLO@0,CT18ZNNLO@0,CT14nlo@0,MMHT2014nlo68clas118@0,MMHT2015qed_nlo@0,NNPDF31_nlo_as_0118_mc@0','--muf=0.5,1.0,2.0','--mur=0.5,1.0,2.0','--dyn=-1'] = systematics_arguments

*******************************Additional Information*******************************

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
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Do you have an ATLAS group that can support you on this?

They should be able to tell you if the lines are expected or not:
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

If not, then they should be the real cause of your issue (and likely related to the way you run MG5aMC in the atlas stack)
If they always occur then we can ignore those.

Did you also try to run such point outside of the atlas wrapper?

Cheers,

Olivier

PS: Sorry I have no idea what the atlas wrapper is doing (and this might be the issue here). If not do not hesitate to tell me

Revision history for this message
Sihyun Jeon (shjeon) said (last edit ):
#2

@Olivier stty: standard input: Inappropriate ioctl for device - this line also shows up in CMSSW environments but usually this is just ignored and mostly works fine for CMS. But I agree that the questioner should run this MadGraph standalone, as it is without relying on ATLAS wrapper and then see if the same error happens

Revision history for this message
Adrián (berrocal) said :
#3

Dear Olivier,

The message seems common, because it also appears in the cases where everything runs.

However, I'll try to contact with people in ATLAS because I understand from your message that it will be complicated to reproduce the issue since we use Madgrapgh integrated in ATLAS wrapper.

Thanks!

Kind regards,

Adrià

Can you help with this problem?

Provide an answer of your own, or ask Adrián for more information if necessary.

To post a message you must log in.