#************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.4 2018-11-09 * #* * #* 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 auto Traceback (most recent call last): File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1501, in onecmd return self.onecmd_orig(line, **opt) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/usr/mg5amcnlo/madgraph/interface/amcatnlo_run_interface.py", line 1686, in do_launch self.exec_cmd('decay_events -from_cards', postcmd=False) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd stop = Cmd.onecmd_orig(current_interface, line, **opt) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/usr/mg5amcnlo/madgraph/interface/common_run_interface.py", line 3660, in do_decay_events madspin_cmd.mg5cmd.exec_cmd( 'set %s %s' %(key,value), errorhandling=False, printcmd=False, precmd=False, postcmd=True) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1530, in exec_cmd stop = current_interface.postcmd(stop, line) File "/usr/mg5amcnlo/madgraph/interface/madgraph_interface.py", line 243, in postcmd stop = super(CmdExtended, self).postcmd(stop, line) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1015, in postcmd stop = getattr(self, 'post_%s' % cmd)(stop, subline) File "/usr/mg5amcnlo/madgraph/interface/madgraph_interface.py", line 7449, in post_set self.exec_cmd('save options %s' % args[0] , log=False) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd stop = Cmd.onecmd_orig(current_interface, line, **opt) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig return func(arg, **opt) File "/usr/mg5amcnlo/madgraph/interface/master_interface.py", line 502, in do_save return self.cmd.do_save(self, *args, **opts) File "/usr/mg5amcnlo/madgraph/interface/madgraph_interface.py", line 7142, in do_save self.write_configuration(filepath, basefile, basedir, to_define) File "/usr/mg5amcnlo/madgraph/interface/extended_cmd.py", line 1985, in write_configuration writer = open(filepath,'w') IOError: [Errno 13] Permission denied: '/usr/mg5amcnlo/input/mg5_configuration.txt' Related File: /usr/mg5amcnlo/input/mg5_configuration.txt Value of current Options: text_editor : None notification_center : True pjfry : None cluster_local_path : None default_unset_couplings : 99 group_subprocesses : Auto ignore_six_quark_processes : False loop_optimized_output : True cluster_status_update : (600, 30) fortran_compiler : None hepmc_path : None collier : /usr/mg5amcnlo/HEPTools/lib auto_update : 7 pythia8_path : /usr/mg5amcnlo/HEPTools/pythia8 hwpp_path : None low_mem_multicore_nlo_generation : False golem : None pythia-pgs_path : None td_path : None delphes_path : None thepeg_path : None cluster_type : condor madanalysis5_path : /usr/mg5amcnlo/HEPTools/madanalysis5/madanalysis5 exrootanalysis_path : None OLP : MadLoop applgrid : applgrid-config eps_viewer : None fastjet : None run_mode : 2 web_browser : None automatic_html_opening : False cluster_temp_path : None cluster_size : 100 cluster_queue : None syscalc_path : None madanalysis_path : None lhapdf : /usr/mg5amcnlo/HEPTools/lhapdf6/bin/lhapdf-config stdout_level : 10 nb_core : 64 f2py_compiler : None ninja : /usr/mg5amcnlo/HEPTools/lib amcfast : amcfast-config cluster_retry_wait : 300 output_dependencies : external crash_on_error : False mg5amc_py8_interface_path : /usr/mg5amcnlo/HEPTools/MG5aMC_PY8_interface loop_color_flows : False samurai : None cluster_nb_retry : 1 mg5_path : /usr/mg5amcnlo timeout : 60 gauge : unitary complex_mass_scheme : False cpp_compiler : None max_npoint_for_channel : 0 #************************************************************ #* MadGraph5_aMC@NLO * #* * #* * * * #* * * * * * #* * * * * 5 * * * * * #* * * * * * #* * * * #* * #* * #* VERSION 2.6.4 2018-11-09 * #* * #* 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 /home/hadh1g17/singularity/2HDM_NLO generate p p > h2 [QCD] output try_madspin ###################################################################### ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL #### ###################################################################### ## ## ## Width set on Auto will be computed following the information ## ## present in the decay.py files of the model. ## ## See arXiv:1402.1178 for more details. ## ## ## ###################################################################### ################################### ## INFORMATION FOR CKMBLOCK ################################### Block ckmblock 1 2.277360e-01 # cabi ################################### ## INFORMATION FOR HIGGS ################################### Block higgs 1 5.000000e-01 # l2 2 1.000000e+00 # l3 3 1.000000e-01 # lR7 4 2.000000e-01 # lI7 5 3.000000e-01 # mixh 6 1.000000e-01 # mixh2 7 2.000000e-01 # mixh3 ################################### ## INFORMATION FOR LOOP ################################### Block loop 1 9.118800e+01 # MU_R ################################### ## INFORMATION FOR MASS ################################### Block mass 5 4.700000e+00 # MB 6 1.720000e+02 # MT 23 9.118760e+01 # MZ 25 1.250000e+02 # mh1 35 7.000000e+02 # mh2 36 7.000000e+02 # mh3 37 7.000000e+02 # mhc ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. 1 0.000000 # d : 0.0 2 0.000000 # u : 0.0 3 0.000000 # s : 0.0 4 0.000000 # c : 0.0 11 0.000000 # e- : 0.0 12 0.000000 # ve : 0.0 13 0.000000 # mu- : 0.0 14 0.000000 # vm : 0.0 15 0.000000 # ta- : 0.0 16 0.000000 # vt : 0.0 21 0.000000 # g : 0.0 22 0.000000 # a : 0.0 24 79.824660 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2))) 9000002 91.187600 # ghz : MZ 9000003 79.824660 # ghwp : MW 9000004 79.824660 # ghwm : MW ################################### ## INFORMATION FOR SMINPUTS ################################### Block sminputs 1 1.279000e+02 # aEWM1 2 1.166390e-05 # Gf 3 1.180000e-01 # aS ################################### ## INFORMATION FOR YUKAWA ################################### Block yukawa 5 4.700000e+00 # ymb 6 1.720000e+02 # ymt ################################### ## INFORMATION FOR YUKAWAGDI ################################### Block yukawagdi 1 1 0.000000e+00 # GDI1x1 1 2 0.000000e+00 # GDI1x2 1 3 0.000000e+00 # GDI1x3 2 1 0.000000e+00 # GDI2x1 2 2 0.000000e+00 # GDI2x2 2 3 0.000000e+00 # GDI2x3 3 1 0.000000e+00 # GDI3x1 3 2 0.000000e+00 # GDI3x2 3 3 0.000000e+00 # GDI3x3 ################################### ## INFORMATION FOR YUKAWAGDR ################################### Block yukawagdr 1 1 1.000000e+00 # GDR1x1 1 2 1.000000e-02 # GDR1x2 1 3 1.000000e-02 # GDR1x3 2 1 1.000000e-02 # GDR2x1 2 2 1.000000e+00 # GDR2x2 2 3 1.000000e-02 # GDR2x3 3 1 1.000000e-02 # GDR3x1 3 2 1.000000e-02 # GDR3x2 3 3 1.000000e+00 # GDR3x3 ################################### ## INFORMATION FOR YUKAWAGLI ################################### Block yukawagli 1 1 0.000000e+00 # GLI1x1 1 2 0.000000e+00 # GLI1x2 1 3 0.000000e+00 # GLI1x3 2 1 0.000000e+00 # GLI2x1 2 2 0.000000e+00 # GLI2x2 2 3 0.000000e+00 # GLI2x3 3 1 0.000000e+00 # GLI3x1 3 2 0.000000e+00 # GLI3x2 3 3 0.000000e+00 # GLI3x3 ################################### ## INFORMATION FOR YUKAWAGLR ################################### Block yukawaglr 1 1 1.000000e+00 # GLR1x1 1 2 1.000000e-02 # GLR1x2 1 3 1.000000e-02 # GLR1x3 2 1 1.000000e-02 # GLR2x1 2 2 1.000000e+00 # GLR2x2 2 3 1.000000e-02 # GLR2x3 3 1 1.000000e-02 # GLR3x1 3 2 1.000000e-02 # GLR3x2 3 3 1.000000e+00 # GLR3x3 ################################### ## INFORMATION FOR YUKAWAGUI ################################### Block yukawagui 1 1 0.000000e+00 # GUI1x1 1 2 0.000000e+00 # GUI1x2 1 3 0.000000e+00 # GUI1x3 2 1 0.000000e+00 # GUI2x1 2 2 0.000000e+00 # GUI2x2 2 3 0.000000e+00 # GUI2x3 3 1 0.000000e+00 # GUI3x1 3 2 0.000000e+00 # GUI3x2 3 3 0.000000e+00 # GUI3x3 ################################### ## INFORMATION FOR YUKAWAGUR ################################### Block yukawagur 1 1 1.000000e+00 # GUR1x1 1 2 1.000000e-02 # GUR1x2 1 3 1.000000e-02 # GUR1x3 2 1 1.000000e-02 # GUR2x1 2 2 1.000000e+00 # GUR2x2 2 3 1.000000e-02 # GUR2x3 3 1 1.000000e-02 # GUR3x1 3 2 1.000000e-02 # GUR3x2 3 3 1.000000e+00 # GUR3x3 ################################### ## INFORMATION FOR DECAY ################################### DECAY 6 1.508336e+00 # WT DECAY 23 2.495200e+00 # WZ DECAY 24 2.085000e+00 # WW DECAY 25 1.000000e+00 # Wh1 DECAY 35 1.000000e+00 # Wh2 DECAY 36 1.000000e+00 # Wh3 DECAY 37 1.000000e+00 # whc ## Dependent parameters, given by model restrictions. ## Those values should be edited following the ## analytical expression. MG5 ignores those values ## but they are important for interfacing the output of MG5 ## to external program such as Pythia. DECAY 1 0.000000 # d : 0.0 DECAY 2 0.000000 # u : 0.0 DECAY 3 0.000000 # s : 0.0 DECAY 4 0.000000 # c : 0.0 DECAY 5 0.000000 # b : 0.0 DECAY 11 0.000000 # e- : 0.0 DECAY 12 0.000000 # ve : 0.0 DECAY 13 0.000000 # mu- : 0.0 DECAY 14 0.000000 # vm : 0.0 DECAY 15 0.000000 # ta- : 0.0 DECAY 16 0.000000 # vt : 0.0 DECAY 21 0.000000 # g : 0.0 DECAY 22 0.000000 # a : 0.0 #=========================================================== # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE) #=========================================================== Block QNUMBERS 9000001 # gha 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) Block QNUMBERS 9000002 # ghz 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) Block QNUMBERS 9000003 # ghwp 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) Block QNUMBERS 9000004 # ghwm 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) Block QNUMBERS 82 # ghg 1 0 # 3 times electric charge 2 1 # number of spin states (2S+1) 3 8 # colour rep (1: singlet, 3: triplet, 8: octet) 4 1 # Particle/Antiparticle distinction (0=own anti) 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) 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) 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) #*********************************************************************** # 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 * #*********************************************************************** 100 = 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). #*********************************************************************** # Normalize the weights of LHE events such that they sum or average to * # the total cross section * #*********************************************************************** average = event_norm ! valid settings: average, sum, bias #*********************************************************************** # 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 * #*********************************************************************** 0 = 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. * #*********************************************************************** nn23nlo = pdlabel ! PDF set 244600 = 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) * #*********************************************************************** False = fixed_ren_scale ! if .true. use fixed ren scale False = 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 -1 = 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 1.0 = 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 True = 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 * #*********************************************************************** False = 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. * # - 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) #*********************************************************************** # 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 #*********************************************************************** # Photon-isolation cuts, according to hep-ph/9801442. 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 #*********************************************************************** # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)* #*********************************************************************** 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids) #***********************************************************************