#************************************************************ #* MadGraph5_aMC@NLO/MadEvent * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.9.4 20xx-xx-xx * #* * #* The MadGraph5_aMC@NLO Development Team - Find us at * #* https://server06.fynu.ucl.ac.be/projects/madgraph * #* * #************************************************************ #* * #* Command File for MadEvent * #* * #* run as ./bin/madevent.py filename * #* * #************************************************************ launch Traceback (most recent call last): File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 1541, in onecmd return self.onecmd_orig(line, **opt) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 1490, in onecmd_orig return func(arg, **opt) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/madevent_interface.py", line 2677, in do_launch self.do_generate_events(line, *args, **opt) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/madevent_interface.py", line 2466, in do_generate_events switch_mode = self.ask_run_configuration(mode, args) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/madevent_interface.py", line 6266, in ask_run_configuration self.ask_edit_cards(cards, plot=False, first_cmd=first_cmd) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/common_run_interface.py", line 985, in ask_edit_cards banner=banner) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/common_run_interface.py", line 1057, in ask_edit_card_static cards=cards, mode=mode, **opt) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 1150, in ask fct=question_instance, fct_timeout=fct_timeout) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 1757, in timed_input result = fct(question) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 2120, in __call__ return self.cmdloop() File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 2288, in cmdloop super(SmartQuestion,self).cmdloop(intro) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/extended_cmd.py", line 179, in cmdloop stop = self.postcmd(stop, line) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/common_run_interface.py", line 6103, in postcmd self.do_update('dependent', timer=20) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/common_run_interface.py", line 6154, in do_update self.paths['param'], timer) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/common_run_interface.py", line 6253, in update_dependent modify = param_card.update_dependent(model, restrict_card, log_level) File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/check_param_card.py", line 520, in update_dependent if not misc.equal(abs(model_value), param_value, 4): File "/Users/souadsemlali/Desktop/MG5_aMC_v2_9_4/mg1/bin/internal/misc.py", line 1529, in equal return ( a==b or abs(int(a*10**power) - int(b*10**power)) < 10) ValueError: invalid literal for int() with base 10: 'autoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoautoaut Run Options ----------- stdout_level : None MadEvent Options ---------------- automatic_html_opening : True notification_center : True run_mode : 2 cluster_queue : None cluster_time : None cluster_size : 100 cluster_memory : None nb_core : 8 (user set) cluster_temp_path : None Configuration Options --------------------- pythia8_path : None (user set) hwpp_path : None (user set) thepeg_path : None (user set) hepmc_path : None (user set) madanalysis_path : None (user set) madanalysis5_path : None (user set) pythia-pgs_path : None (user set) td_path : None (user set) delphes_path : None (user set) exrootanalysis_path : None (user set) syscalc_path : None (user set) lhapdf : lhapdf-config lhapdf_py2 : None lhapdf_py3 : None 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 : 7 (user set) cluster_type : condor cluster_status_update : (600, 30) cluster_nb_retry : 1 cluster_local_path : None cluster_retry_wait : 300 #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.9.4 2021-05-30 * #* * #* 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 group_subprocesses Auto set ignore_six_quark_processes False set low_mem_multicore_nlo_generation False set complex_mass_scheme False set gauge unitary set loop_optimized_output True set loop_color_flows False set max_npoint_for_channel 0 set default_unset_couplings 99 set max_t_for_channel 99 set zerowidth_tchannel 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+ define l- = e- mu- define vl = ve vm vt define vl~ = ve~ vm~ vt~ import model THDM_type1_UFO__hgg_plugin generate g g > h2 z > h3 mu+ mu-, (h3 > z > h1 j j, h1 > b b~) output mg1 ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 #### ###################################################################### ######### ## INFORMATION FOR FRBLOCK ################################### BLOCK FRBLOCK # 1 6.300000e+00 # tb 2 -1.600000e-01 # sinbma ################################### ## INFORMATION FOR HIGGS ################################### BLOCK HIGGS # 1 2.640000e-01 # l2 2 5.920000e-01 # l3 3 0.000000e-00 # lr7 ############## ## INFORMATION FOR MASS ################################### BLOCK MASS # 1 5.040000e-03 # md 2 2.550000e-03 # mu 3 1.010000e-01 # ms 4 1.270000e+00 # mc 5 4.700000e+00 # mb 6 1.720000e+02 # mt 11 5.110000e-04 # me 13 1.056600e-01 # mmu 15 1.777000e+00 # mta 23 9.118760e+01 # mz 25 1.268000e+01 # mh1 35 1.250000e+02 # mh2 36 6.887000e+01 # mh3 37 1.350700e+02 # mhc 9000006 1.200000e+02 # mp 12 0.000000e+00 # ve : 0.0 14 0.000000e+00 # vm : 0.0 16 0.000000e+00 # vt : 0.0 21 0.000000e+00 # g : 0.0 22 0.000000e+00 # a : 0.0 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))) ################################### ## INFORMATION FOR SMINPUTS ################################### BLOCK SMINPUTS # 1 1.279000e+02 # aewm1 2 1.166390e-05 # gf 3 1.180000e-01 # as (note that parameter not used if you use a pdf set) ################################### ## INFORMATION FOR YUKAWA ################################### BLOCK YUKAWA # 1 5.040000e-03 # ymdo 2 2.550000e-03 # ymup 3 1.010000e-01 # yms 4 1.270000e+00 # ymc 5 4.700000e+00 # ymb 6 1.720000e+02 # ymt 11 5.110000e-04 # yme 13 1.056600e-01 # ymm 15 1.777000e+00 # ymtau ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.508336e+00 # wt DECAY 23 2.495200e+00 # wz DECAY 24 2.085000e+00 # ww DECAY 25 Auto # wh1 DECAY 35 4.520000e-03 # wh2 DECAY 36 1.130000e-04 # wh3 DECAY 37 1.000000e+00 # whc DECAY 9000006 5.753088e-03 # wh1__2 DECAY 1 0.000000e+00 # d : 0.0 DECAY 2 0.000000e+00 # u : 0.0 DECAY 3 0.000000e+00 # s : 0.0 DECAY 4 0.000000e+00 # c : 0.0 DECAY 5 0.000000e+00 # b : 0.0 DECAY 11 0.000000e+00 # e- : 0.0 DECAY 12 0.000000e+00 # ve : 0.0 DECAY 13 0.000000e+00 # mu- : 0.0 DECAY 14 0.000000e+00 # vm : 0.0 DECAY 15 0.000000e+00 # ta- : 0.0 DECAY 16 0.000000e+00 # vt : 0.0 DECAY 21 0.000000e+00 # g : 0.0 DECAY 22 0.000000e+00 # a : 0.0 ################################### ## INFORMATION FOR QNUMBERS 37 ################################### BLOCK QNUMBERS 37 # h+ 1 3 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 35 ################################### BLOCK QNUMBERS 35 # h2 1 0 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 36 ################################### BLOCK QNUMBERS 36 # h3 1 0 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 0 # particle/antiparticle distinction (0=own anti) ################################### ## INFORMATION FOR QNUMBERS 9000006 ################################### BLOCK QNUMBERS 9000006 # h1__1 1 0 # 3 times electric charge 2 1 # number of spin states (2s+1) 3 1 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # particle/antiparticle distinction (0=own anti) #********************************************************************* # MadGraph5_aMC@NLO * # * # run_card.dat MadEvent * # * # This file is used to set the parameters of the run. * # * # Some notation/conventions: * # * # Lines starting with a '# ' are info or comments * # * # mind the format: value = variable ! comment * # * # To display more options, you can type the command: * # update full_run_card * #********************************************************************* # #********************************************************************* # Tag name for the run (one word) * #********************************************************************* tag_1 = run_tag ! name of the run #********************************************************************* # Number of events and rnd seed * # Warning: Do not generate more than 1M events in a single run * #********************************************************************* 10000 = nevents ! Number of unweighted events requested 0 = iseed ! rnd seed (0=assigned automatically=default)) #********************************************************************* # Collider type and energy * # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, * # 3=photon from electron, 4=photon from muon * #********************************************************************* 1 = lpp1 ! beam 1 type 1 = lpp2 ! beam 2 type 6500.0 = ebeam1 ! beam 1 total energy in GeV 6500.0 = ebeam2 ! beam 2 total energy in GeV # To see polarised beam options: type "update beam_pol" #********************************************************************* # PDF CHOICE: this automatically fixes also alpha_s and its evol. * #********************************************************************* nn23lo1 = pdlabel ! PDF set 230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number # To see heavy ion options: type "update ion_pdf" #********************************************************************* # Renormalization and factorization scales * #********************************************************************* False = fixed_ren_scale ! if .true. use fixed ren scale False = fixed_fac_scale ! if .true. use fixed fac scale 91.188 = scale ! fixed ren scale 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2 -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices 1.0 = scalefact ! scale factor for event-by-event scales #********************************************************************* # Type and output format #********************************************************************* False = gridpack !True = setting up the grid pack -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. Normalization of the weight in the LHEF # To see MLM/CKKW merging options: type "update MLM" or "update CKKW" #********************************************************************* # #********************************************************************* # Phase-Space Optimization strategy (basic options) #********************************************************************* 0 = nhel ! using helicities importance sampling or not. ! 0: sum over helicity, 1: importance sampling 2 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx) ! 1 is old strategy (using amp square) ! 2 is new strategy (using only the denominator) # To see advanced option for Phase-Space optimization: type "update psoptim" #********************************************************************* # Generation bias, check the wiki page below for more information: * # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' * #********************************************************************* None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-] {} = bias_parameters ! Specifies the parameters of the module. # #******************************* # Parton level cuts definition * #******************************* # # #********************************************************************* # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay #********************************************************************* 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma) #********************************************************************* # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply) #********************************************************************* False = cut_decays ! Cut decay products #********************************************************************* # Standard Cuts * #********************************************************************* # Minimum and maximum pt's (for max, -1 means no cut) * #********************************************************************* 20.0 = ptj ! minimum pt for the jets 0.0 = ptb ! minimum pt for the b 10.0 = ptl ! minimum pt for the charged leptons -1.0 = ptjmax ! maximum pt for the jets -1.0 = ptbmax ! maximum pt for the b -1.0 = ptlmax ! maximum pt for the charged leptons {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) # # For display option for energy cut in the partonic center of mass frame type 'update ecut' # #********************************************************************* # Maximum and minimum absolute rapidity (for max, -1 means no cut) * #********************************************************************* 5.0 = etaj ! max rap for the jets -1.0 = etab ! max rap for the b 2.5 = etal ! max rap for the charged leptons 0.0 = etabmin ! min rap for the b 0.0 = etalmin ! main rap for the charged leptons {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5}) #********************************************************************* # Minimum and maximum DeltaR distance * #********************************************************************* 0.4 = drjj ! min distance between jets 0.0 = drbb ! min distance between b's 0.4 = drll ! min distance between leptons 0.0 = drbj ! min distance between b and jet 0.4 = drjl ! min distance between jet and lepton 0.0 = drbl ! min distance between b and lepton -1.0 = drjjmax ! max distance between jets -1.0 = drbbmax ! max distance between b's -1.0 = drllmax ! max distance between leptons -1.0 = drbjmax ! max distance between b and jet -1.0 = drjlmax ! max distance between jet and lepton -1.0 = drblmax ! max distance between b and lepton #********************************************************************* # Minimum and maximum invariant mass for pairs * #********************************************************************* 0.0 = mmjj ! min invariant mass of a jet pair 0.0 = mmbb ! min invariant mass of a b pair 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair -1.0 = mmjjmax ! max invariant mass of a jet pair -1.0 = mmbbmax ! max invariant mass of a b pair -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250}) {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only ! to pairs of particle/antiparticle and not to pairs of the same pdg codes. #********************************************************************* # Minimum and maximum invariant mass for all letpons * #********************************************************************* 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl) -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl) #********************************************************************* # Minimum and maximum pt for 4-momenta sum of leptons / neutrino * # for pair of lepton includes only same flavor, opposite charge #********************************************************************* 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl) -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl) #********************************************************************* # Inclusive cuts * #********************************************************************* 0.0 = xptj ! minimum pt for at least one jet 0.0 = xptb ! minimum pt for at least one b 0.0 = xptl ! minimum pt for at least one charged lepton #********************************************************************* # Control the pt's of the jets sorted by pt * #********************************************************************* 0.0 = ptj1min ! minimum pt for the leading jet in pt 0.0 = ptj2min ! minimum pt for the second jet in pt -1.0 = ptj1max ! maximum pt for the leading jet in pt -1.0 = ptj2max ! maximum pt for the second jet in pt 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts #********************************************************************* # Control the pt's of leptons sorted by pt * #********************************************************************* 0.0 = ptl1min ! minimum pt for the leading lepton in pt 0.0 = ptl2min ! minimum pt for the second lepton in pt -1.0 = ptl1max ! maximum pt for the leading lepton in pt -1.0 = ptl2max ! maximum pt for the second lepton in pt #********************************************************************* # Control the Ht(k)=Sum of k leading jets * #********************************************************************* 0.0 = htjmin ! minimum jet HT=Sum(jet pt) -1.0 = htjmax ! maximum jet HT=Sum(jet pt) 0.0 = ihtmin !inclusive Ht for all partons (including b) -1.0 = ihtmax !inclusive Ht for all partons (including b) 0.0 = ht2min ! minimum Ht for the two leading jets 0.0 = ht3min ! minimum Ht for the three leading jets 0.0 = ht4min ! minimum Ht for the four leading jets -1.0 = ht2max ! maximum Ht for the two leading jets -1.0 = ht3max ! maximum Ht for the three leading jets -1.0 = ht4max ! maximum Ht for the four leading jets #********************************************************************* # WBF cuts * #********************************************************************* 0.0 = xetamin ! minimum rapidity for two jets in the WBF case 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case #********************************************************************* # maximal pdg code for quark to be considered as a light jet * # (otherwise b cuts are applied) * #********************************************************************* 4 = maxjetflavor ! Maximum jet pdg code #********************************************************************* # #********************************************************************* # Store info for systematics studies * # WARNING: Do not use for interference type of computation * #********************************************************************* True = use_syst ! Enable systematics studies # systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++] ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule # Syscalc is deprecated but to see the associate options type'update syscalc'