Cross sections ~10^9 pb at 14TeV

Asked by Ciara Byers

I am using the model 2HDMtII_NLO ( found here http://feynrules.irmp.ucl.ac.be/wiki/2HDM) and running the process bq > qth_i (where q is defined as the 'light' quarks and h_i is defined as all the higgs in the 2HDM except the charged one). I am using version 2.7.2 of madgraph with python 2.7.5. , lhapdf6 and madanalysis5. Clearly the cross sections coming out as 10^9 pb are not correct, but I am unsure as to what could be causing this.

One thing I have noticed is a confusing error displayed at the start of the "run_01_tag_1_banner.txt" file:

#*********************************************************************
# MadGraph/MadEvent *
# http://madgraph.hep.uiuc.edu *
# *
# proc_card.dat *
#*********************************************************************
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#*********************************************************************
#*********************************************************************
# Process(es) requested : mg2 input *
#*********************************************************************
# Begin PROCESS # This is TAG. Do not modify this line
ql b > ql h_i t / z g #Process
# Be carefull the coupling are here in MG5 convention

I am not actually concerned about running madevent, and I try to turn off MadAnalysis using my proc_card. My proc_card looks like this:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.7.2 2020-03-17 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************

#set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
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+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model /scratch/cb27g11/MG5_aMC_v2_7_2/models/2HDMtII_NLO
define ql = u u~ d d~ s s~ c c~
define h_i = h1 h2 h3
generate ql b > ql h_i t / z g
output /scratch/cb27g11/bq_tqhi_-0.998586_144813.0
launch /scratch/cb27g11/bq_tqhi_-0.998586_144813.0
    3
    done
    set param_card frblock 1 -0.998586 #tan(beta) value
    set param_card frblock 2 144813.0 #sin(beta-alpha) value
    set param_card mass 35 125.0 #H mass value
    set param_card mass 36 705.649 #A mass value
    set param_card mass 37 692.03 #charged higgs mass value
    set ebeam1 7000
    set ebeam2 7000
    done

Everything else is left to the default value, the pdf used is the default one etc. The parameter point I am running is one that has been checked in HiggsBounds/HiggsSignals and 2HDMC. Do you have any thoughts about what might be causing this?

Kind regards,
Ciara

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Ciara Byers
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

I have run this with our LTS version (2.9.x)
and I got the following cross-section:

> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 0.028 +- 8.708e-05 pb
> Nb of events : 10000
>
> INFO: Running Systematics computation
> INFO: Idle: 2, Running: 2, Completed: 0 [ current time: 13h32 ]
> INFO: Idle: 0, Running: 3, Completed: 1 [ 30.2s ]
> INFO: # events generated with PDF: NNPDF23_lo_as_0130_qed (247000)
> INFO: #Will Compute 145 weights per event.
> INFO: #***************************************************************************
> #
> # original cross-section: 0.027996459901840486
> # scale variation: +3.71% -6.04%
> # central scheme variation: + 0% -4.38%
> # PDF variation: +1.78% -1.78%
> #
> # dynamical scheme # 1 : 0.0271424 +2.82% -4.68% # \sum ET
> # dynamical scheme # 2 : 0.0277143 +1.9% -3.4% # \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 3 : 0.026771 +3.52% -5.67% # 0.5 \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 4 : 0.0277708 +1.07% -2.18% # \sqrt{\hat s}
> #***************************************************************************

Since you use a quite old version, this is likely an old bug that I do not remember.

Cheers,

Olivier

PS: Madeven is the package that compute the LO cross-section within MG5aMC, so you do care about it.
But the warning is just stating that we do include a file which follow MG4 convention because it is needed by madevent but this file should not be used within MG4

> On 15 Jun 2021, at 12:45, Ciara Byers <email address hidden> wrote:
>
> set param_card frblock 1 -0.998586 #tan(beta) value
> set param_card frblock 2 144813.0 #sin(beta-alpha) value
> set param_card mass 35 125.0 #H mass value
> set param_card mass 36 705.649 #A mass value
> set param_card mass 37 692.03 #charged higgs mass value
> set ebeam1 7000
> set ebeam2 7000

Revision history for this message
Ciara Byers (ciarab93) said :
#2

Hi Olivier,

Thanks for getting back to me so quickly. I only use this older version because I inherited a toolbox from a predecessor that used it. Am I correct in thinking that the most recent version is written using python that is compatible with python2 and 3? I'm happy to upgrade, but I run madgraph on a supercomputer and it can be very finicky when it comes to different versions of python. I can only have one version activated at a time, so all my code has been written to work with python 2.7.5, I'm under a bit of pressure to get data at the moment so I don't want to convert all my code to python 3 right now if it can be avoided.

Best wishes,
Ciara

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

Yes we are still python2.7 compatible.
But you have to specify that you want to use it since we have set our default python executable to python3.

so you need to run mg5aMC like
python2.7 ./bin/mg5_aMc

Cheers,

Olivier

> On 15 Jun 2021, at 15:40, Ciara Byers <email address hidden> wrote:
>
> Question #697543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/697543
>
> Status: Answered => Open
>
> Ciara Byers is still having a problem:
> Hi Olivier,
>
> Thanks for getting back to me so quickly. I only use this older version
> because I inherited a toolbox from a predecessor that used it. Am I
> correct in thinking that the most recent version is written using python
> that is compatible with python2 and 3? I'm happy to upgrade, but I run
> madgraph on a supercomputer and it can be very finicky when it comes to
> different versions of python. I can only have one version activated at a
> time, so all my code has been written to work with python 2.7.5, I'm
> under a bit of pressure to get data at the moment so I don't want to
> convert all my code to python 3 right now if it can be avoided.
>
> Best wishes,
> Ciara
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Ciara Byers (ciarab93) said :
#4

Hi,

Great! Thanks for your help.

Best wishes,
Ciara

Revision history for this message
Ciara Byers (ciarab93) said :
#5

Hi again,

I'm actually having difficulty updating my install using the "install update" command in madgraph itself. I get the following:

move Template/SubProcesses/addmothers.f to madgraph/iolibs/template_files/addmothers.f
Command "install update" interrupted with error:
OSError : [Errno 2] No such file or directory: '/scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

When I look at MG5_debug I have the following:

Traceback (most recent call last):
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/master_interface.py", line 470, in do_install
    self.cmd.do_install(self, *args, **opts)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6041, in do_install
    self.install_update(['update']+install_options['update_options'],wget=program)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6757, in install_update
    need_binary = apply_patch(filetext)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6495, in apply_patch
    files.cp(pjoin(MG5DIR, orig), pjoin(MG5DIR, new), error=True)
  File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/iolibs/files.py", line 144, in cp
    shutil.copytree(path1, path2)
  File "/local/software/python/2.7.5/lib/python2.7/shutil.py", line 171, in copytree
    names = os.listdir(src)
OSError: [Errno 2] No such file or directory: '/scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f'
Related File: /scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f
                          MadGraph5_aMC@NLO Options
                          ----------------
        complex_mass_scheme : False
    default_unset_couplings : 99
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0
               stdout_level : 20 (user set)

                         MadEvent Options
                          ----------------
     automatic_html_opening : True
nb_core : None
        notification_center : True
                   run_mode : 2

                      Configuration Options
                      ---------------------
                        OLP : MadLoop
                    amcfast : amcfast-config
                   applgrid : applgrid-config
                auto_update : 7
         cluster_local_path : None
           cluster_nb_retry : 1
              cluster_queue : None (user set)
         cluster_retry_wait : 300
               cluster_size : 100
      cluster_status_update : (600, 30)
          cluster_temp_path : None
               cluster_type : condor
                    collier : ./HEPTools/lib
               cpp_compiler : None
             crash_on_error : False
               delphes_path : ./Delphes
                 eps_viewer : None
        exrootanalysis_path : ./ExRootAnalysis
              f2py_compiler : None
                    fastjet : None (user set)
           fortran_compiler : None
                      golem : None (user set)
                 hepmc_path : None (user set)
                  hwpp_path : None (user set)
                     lhapdf : /scratch/cb27g11/MG5_aMC_v2_7_2/HEPTools/lhapdf6/bin/lhapdf-config (user set)
          madanalysis5_path : /scratch/cb27g11/MG5_aMC_v2_7_2/HEPTools/madanalysis5/madanalysis5 (user set)
           madanalysis_path : ./MadAnalysis
  mg5amc_py8_interface_path : None (user set)
                      ninja : ./HEPTools/lib
        output_dependencies : external
            pythia-pgs_path : ./pythia-pgs
               pythia8_path : None (user set)
                    samurai : None
               syscalc_path : ./SysCalc
                    td_path : ./td
                text_editor : None
                thepeg_path : None (user set)
                    timeout : 60
         web_browser : None

if the problem is just a missing file can I just download MG5 again and copy it from there, or am I best to just start with a clean install?

Best wishes,
Ciara

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#6

Starting clean is probably better here.
The move to python3 (and to version3) was creating a lot of issue for that function.

Sorry,

Olivier

> On 15 Jun 2021, at 16:55, Ciara Byers <email address hidden> wrote:
>
> Question #697543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/697543
>
> Status: Answered => Open
>
> Ciara Byers is still having a problem:
> Hi again,
>
> I'm actually having difficulty updating my install using the "install
> update" command in madgraph itself. I get the following:
>
> move Template/SubProcesses/addmothers.f to madgraph/iolibs/template_files/addmothers.f
> Command "install update" interrupted with error:
> OSError : [Errno 2] No such file or directory: '/scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> When I look at MG5_debug I have the following:
>
> Traceback (most recent call last):
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1515, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
> return func(arg, **opt)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/master_interface.py", line 470, in do_install
> self.cmd.do_install(self, *args, **opts)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6041, in do_install
> self.install_update(['update']+install_options['update_options'],wget=program)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6757, in install_update
> need_binary = apply_patch(filetext)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/interface/madgraph_interface.py", line 6495, in apply_patch
> files.cp(pjoin(MG5DIR, orig), pjoin(MG5DIR, new), error=True)
> File "/scratch/cb27g11/MG5_aMC_v2_7_2/madgraph/iolibs/files.py", line 144, in cp
> shutil.copytree(path1, path2)
> File "/local/software/python/2.7.5/lib/python2.7/shutil.py", line 171, in copytree
> names = os.listdir(src)
> OSError: [Errno 2] No such file or directory: '/scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f'
> Related File: /scratch/cb27g11/MG5_aMC_v2_7_2/Template/SubProcesses/addmothers.f
> MadGraph5_aMC@NLO Options
> ----------------
> complex_mass_scheme : False
> default_unset_couplings : 99
> gauge : unitary
> group_subprocesses : Auto
> ignore_six_quark_processes : False
> loop_color_flows : False
> loop_optimized_output : True
> low_mem_multicore_nlo_generation : False
> max_npoint_for_channel : 0
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : True
> nb_core : None
> notification_center : True
> run_mode : 2
>
> Configuration Options
> ---------------------
> OLP : MadLoop
> amcfast : amcfast-config
> applgrid : applgrid-config
> auto_update : 7
> cluster_local_path : None
> cluster_nb_retry : 1
> cluster_queue : None (user set)
> cluster_retry_wait : 300
> cluster_size : 100
> cluster_status_update : (600, 30)
> cluster_temp_path : None
> cluster_type : condor
> collier : ./HEPTools/lib
> cpp_compiler : None
> crash_on_error : False
> delphes_path : ./Delphes
> eps_viewer : None
> exrootanalysis_path : ./ExRootAnalysis
> f2py_compiler : None
> fastjet : None (user set)
> fortran_compiler : None
> golem : None (user set)
> hepmc_path : None (user set)
> hwpp_path : None (user set)
> lhapdf : /scratch/cb27g11/MG5_aMC_v2_7_2/HEPTools/lhapdf6/bin/lhapdf-config (user set)
> madanalysis5_path : /scratch/cb27g11/MG5_aMC_v2_7_2/HEPTools/madanalysis5/madanalysis5 (user set)
> madanalysis_path : ./MadAnalysis
> mg5amc_py8_interface_path : None (user set)
> ninja : ./HEPTools/lib
> output_dependencies : external
> pythia-pgs_path : ./pythia-pgs
> pythia8_path : None (user set)
> samurai : None
> syscalc_path : ./SysCalc
> td_path : ./td
> text_editor : None
> thepeg_path : None (user set)
> timeout : 60
> web_browser : None
>
> if the problem is just a missing file can I just download MG5 again and
> copy it from there, or am I best to just start with a clean install?
>
> Best wishes,
> Ciara
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Ciara Byers (ciarab93) said :
#7

Hi Olivier,

I've reinstalled MadGraph, I'm using v3.1.0 now but I'm having the same problem I was having when using v2.7.2

I am once again using the model 2HDMtII_NLO, and the same proc card (with the path altered to lead to v3.1.0 of course, in fact I have completely removed v2.7.2 in order to make sure it didn't accidentally get called).

Could this be a problem with the model itself? I had redownloaded it after reinstalling MadGraph, so I would've thought it would be fine?

Best wishes,
Ciara

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#8

Hi,

I'm really sorry but I actually reproduce your error (also with the same version as the one mentioned above).
I should have made a mistake in my test at that time. I'm so sorry for that.

Now the issue seems to be in the model itself since the coupling used for the computation are:

       GC_118 = -0.32420E+00 0.00000E+00
       GC_119 = 0.32420E+00 0.00000E+00
       GC_159 = 0.00000E+00 0.74953E+07
       GC_161 = 0.00000E+00 0.00000E+00
       GC_199 = 0.19062E-01 0.00000E+00
       GC_200 = -0.00000E+00 -0.27643E+04
       GC_201 = -0.00000E+00 0.27604E+04
       GC_202 = 0.69956E+00 0.00000E+00
       GC_206 = 0.00000E+00 -0.10116E+06
       GC_207 = -0.00000E+00 -0.10130E+06
        GC_21 = 0.00000E+00 0.98932E+00
        GC_22 = 0.00000E+00 0.26957E-01
        GC_77 = -0.00000E+00 -0.00000E+00
        GC_78 = 0.00000E+00 0.00000E+00
         GC_9 = 0.00000E+00 0.45849E+00
        GC_92 = -0.00000E+00 -0.46949E+05
        GC_93 = 0.00000E+00 0.46949E+05

A lot of time are very large and likely the reason why your cross-section is so large.

Cheers,

Olivier

Revision history for this message
Ciara Byers (ciarab93) said :
#9

Hi Olivier,

Ah I see. I guess I may have to consider using a different model then.

Before giving up on this one though, could I ask a few questions to help me decode how the expressions for couplings are written in couplings.py ?

e.g.

GC_159 is the coupling for h1-WW, so I'm expecting an expression along the lines of (gv)^2 * WW/4

What I see is:
'(complex(0,1)*gw**2*TH1x1*vev)/2.'

It's probably very obvious, but I'm not sure what complex(0,1) and TH1x1 are, or where the latter might be defined?

It seems like it would be as easy for me to try and check these coupling expressions as it would be to create a new model, or try to alter another 2HDM model file to allow me direct access to sin(beta - alpha) and tan(beta)? I did note there's a feynrules file on the website that would be suitable, but I would need to learn how to convert it to a normal madgraph model, do you think I'm right in thinking I should try to fix my current model first?

Best wishes,
Ciara

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#10

Hi,

complex(0,1) is the pure imaginary complex number.

For the rest you can do
"display parameters" to see what they refer to:
You will see for example that
mdl_TH1x1 = mdl_sinbma

and the default value for mdl_sinbma (that you can change in the param_card) is
mdl_sinbma = 0.1

Cheers,

Olivier

> On 17 Jun 2021, at 18:50, Ciara Byers <email address hidden> wrote:
>
> Question #697543 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/697543
>
> Status: Answered => Open
>
> Ciara Byers is still having a problem:
> Hi Olivier,
>
> Ah I see. I guess I may have to consider using a different model then.
>
> Before giving up on this one though, could I ask a few questions to help
> me decode how the expressions for couplings are written in couplings.py
> ?
>
> e.g.
>
> GC_159 is the coupling for h1-WW, so I'm expecting an expression along
> the lines of (gv)^2 * WW/4
>
> What I see is:
> '(complex(0,1)*gw**2*TH1x1*vev)/2.'
>
> It's probably very obvious, but I'm not sure what complex(0,1) and TH1x1
> are, or where the latter might be defined?
>
> It seems like it would be as easy for me to try and check these coupling
> expressions as it would be to create a new model, or try to alter
> another 2HDM model file to allow me direct access to sin(beta - alpha)
> and tan(beta)? I did note there's a feynrules file on the website that
> would be suitable, but I would need to learn how to convert it to a
> normal madgraph model, do you think I'm right in thinking I should try
> to fix my current model first?
>
> Best wishes,
> Ciara
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#11

Maybe a naive question here (since I do not know the model)
is why you do the following
set param_card frblock 2 144813.0 #sin(beta-alpha) value

I would expect that a sinus as input would be between -1 and 1.

Cheers,

Olivier

Revision history for this message
Ciara Byers (ciarab93) said :
#12

You may have solved the problem in spotting that!

My parameter values are generated in a MCMC and then run through HB/HS and 2HDMC. Looks like I forgot that a new column is added to the start of my csv files, so the values are all off. All the real sinbma values are between -1 and 1.

Thank-you for having keener eyes than me!

Best wishes,
Ciara