interrupted with error: KeyError : 'scale_upp'

Hi,

The version of the code that you are using is about two years old.
I guess that such bug was solved at least 20 month ago (probably in 2.3.0)
And that part of the code has been deeply restructured 5 month ago.
So the only advise, that I can give to you is to update your code to the latest version.

Cheers,

Olivier

> On Dec 22, 2016, at 17:47, Sylvie Braibant <email address hidden> wrote:
>
> New question #408179 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/408179
>
> Dear Madgraph authors,
>
> I would like to generate HH events using the reweighting techniche for NLO.
> When I tried to simulate 100000 events, the generation quits with the error
> KeyError : 'scale_upp'.
> I tired to lower the number of events per job per still does not work.
> Thanks for any help.
> Sylvie
>
>
> Command "launch -p --name=run100k -f" interrupted with error:
> KeyError : 'scale_upp'
> Please report this bug on https://bugs.launchpad.net/madgraph5
> More information is found in '/gpfs_data/local/cms/braibant/FCC/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/run100k_tag_1_debug.log'.
> Please attach this file to your report.
> quit
>
>
> ************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.2.2.3 *
> #* *
> #* 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 -p --name=test -f
> Traceback (most recent call last):
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/extended_cmd.py", line 879, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/extended_cmd.py", line 872, in onecmd_orig
> return func(arg, **opt)
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/amcatnlo_run_interface.py", line 1192, in do_launch
> evt_file = self.run(mode, options)
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/amcatnlo_run_interface.py", line 1592, in run
> return self.reweight_and_collect_events(options, mode, nevents, event_norm)
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/amcatnlo_run_interface.py", line 2318, in reweight_and_collect_events
> self.print_summary(options, 2, mode, scale_pdf_info)
> File "/home/CMS-T3/braibant/ggHH_NLO/PROCNLO_SMEFT_UFO-no_b_mass_1/bin/internal/amcatnlo_run_interface.py", line 1796, in print_summary
> (scale_pdf_info['scale_upp'], scale_pdf_info['scale_low'])
> KeyError: 'scale_upp'
> Value of current Options:
> text_editor : None
> web_browser : None
> cluster_temp_path : None
> timeout : 60
> cluster_queue : madgraph
> madanalysis_path : None
> lhapdf : /cvmfs/cms.cern.ch/slc6_amd64_gcc530/external/lhapdf/6.1.6-giojec7/bin/lhapdf-config
> mg5_path : None
> cluster_memory : None
> cluster_status_update : (600, 30)
> cluster_time : None
> hepmc_path : None
> pythia8_path : None
> hwpp_path : None
> automatic_html_opening : False
> cluster_retry_wait : 300
> stdout_level : None
> pythia-pgs_path : None
> td_path : None
> delphes_path : None
> thepeg_path : None
> cluster_type : condor
> cluster_nb_retry : 1
> fortran_compiler : None
> nb_core : 8
> auto_update : 7
> exrootanalysis_path : None
> eps_viewer : None
> syscalc_path : None
> fastjet : /cvmfs/cms.cern.ch/slc6_amd64_gcc530/external/fastjet/3.1.0/bin/fastjet-config
> cpp_compiler : None
> run_mode : 2
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.2.3 2015-02-10 *
> #* *
> #* 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 loop_optimized_output True
> set complex_mass_scheme False
> 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~
> set complex_mass_scheme True
> import model SMEFT_UFO-no_b_mass
> define p = p b b~
> generate p p > H H [real= QCD]
> output
> ######################################################################
> ## 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 MASS
> ###################################
> Block mass
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.118760e+01 # MZ
> 24 7.982436e+01 # MW
> 25 1.250000e+02 # MH
> ## 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
> 5 0.000000 # b : 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
> 16 0.000000 # vt : 0.0
> 21 0.000000 # g : 0.0
> 22 0.000000 # a : 0.0
> 250 91.187600 # g0 : MZ
> 251 79.824360 # g+ : MW
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 1.279000e+02 # aEWM1
> 3 1.184000e-01 # aS
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.508336e+00 # WT
> DECAY 23 2.495200e+00 # WZ
> DECAY 24 2.085000e+00 # WW
> DECAY 25 5.753088e-03 # WH
> ## 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
> DECAY 250 2.495200 # g0 : WZ
> DECAY 251 2.085000 # g+ : WW
> #===========================================================
> # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
> #===========================================================
>
> Block QNUMBERS 250 # g0
> 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 251 # g+
> 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)
> #***********************************************************************
> # 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 *
> #***********************************************************************
> #
> #*******************
> # 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 *
> #***********************************************************************
> 100000 = nevents ! Number of unweighted events requested
> -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 100 = 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 ! average or sum
> #***********************************************************************
> # 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)
> 50000 = ebeam1 ! beam 1 energy in GeV
> 50000 = ebeam2 ! beam 2 energy in GeV
> #***********************************************************************
> # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> #***********************************************************************
> lhapdf = pdlabel ! PDF set
> 11000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> #***********************************************************************
> # 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!!!! *
> #***********************************************************************
> HERWIGPP = parton_shower
> #***********************************************************************
> # Renormalization and factorization scales *
> # (Default functional form for the non-fixed scales is the sum of *
> # the transverse masses of all final state particles and partons. This *
> # can be changed in SubProcesses/set_scales.f) *
> #***********************************************************************
> F = fixed_ren_scale ! if .true. use fixed ren scale
> F = fixed_fac_scale ! if .true. use fixed fac scale
> 91.188 = muR_ref_fixed ! fixed ren reference scale
> 91.188 = muF1_ref_fixed ! fixed fact reference scale for pdf1
> 91.188 = muF2_ref_fixed ! fixed fact reference scale for pdf2
> #***********************************************************************
> # Renormalization and factorization scales (advanced and NLO options) *
> #***********************************************************************
> F = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
> 91.188 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
> 1 = muR_over_ref ! ratio of current muR over reference muR
> 1 = muF1_over_ref ! ratio of current muF1 over reference muF1
> 1 = muF2_over_ref ! ratio of current muF2 over reference muF2
> 1 = QES_over_ref ! ratio of current QES over reference QES
> #***********************************************************************
> # Reweight flags to get scale dependence and PDF uncertainty *
> # For scale dependence: factor rw_scale_up/down around central scale *
> # For PDF uncertainty: use LHAPDF with supported set *
> #***********************************************************************
> .true. = reweight_scale ! reweight to get scale dependence
> 0.5 = rw_Rscale_down ! lower bound for ren scale variations
> 2.0 = rw_Rscale_up ! upper bound for ren scale variations
> 0.5 = rw_Fscale_down ! lower bound for fact scale variations
> 2.0 = rw_Fscale_up ! upper bound for fact scale variations
> .false. = reweight_PDF ! reweight to get PDF uncertainty
> 11001 = PDF_set_min ! First of the error PDF sets
> 11052 = PDF_set_max ! Last of the error PDF sets
> #***********************************************************************
> # 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 more details. *
> #***********************************************************************
> 0 = ickkw ! 0 no merging, 3 FxFx merging, 4 UNLOPS
> #***********************************************************************
> #
> #***********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) *
> #***********************************************************************
> 15 = 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 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 0.4 = jetradius ! The radius parameter for the jet algorithm
> 20 = ptj ! Min jet transverse momentum
> -1 = 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 gen cuts can be specified in SubProcesses/cuts.f) *
> #***********************************************************************
> 0 = ptl ! Min lepton transverse momentum
> -1 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 0 = drll ! Min distance between opposite sign lepton pairs
> 0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 30 = 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 *
> #***********************************************************************
> 20 = ptgmin ! Min photon transverse momentum
> -1 = 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)
> #***********************************************************************
> # Maximal PDG code for quark to be considered a jet when applying cuts.*
> # At least all massless quarks of the model should be included here. *
> #***********************************************************************
> 4 = maxjetflavor
> #***********************************************************************
> # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> #***********************************************************************
> 0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
> #***********************************************************************
>
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.