Problems with reweighting model parameters

The typical issue is that f2py executable is not compatible with your python version that you are using.

To have more details you can do
cd rw_me/Subprocesses/P1_uux_epem
make matrix2py.so

Cheers,

Olivier

> On 14 Jul 2021, at 16:41, Philip Keicher <email address hidden> wrote:
>
> New question #697979 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/697979
>
> Dear authors,
>
> I'm trying to generate events and would like to use the reweight method to check the impact of a set of model parameters. However, I ran into problems with the reweighting. I'm using MG5_aMC_v3_1_0.
>
> In order to understand this, I generated p p > e+ e- at LO and tried to do the reweighting described here: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Reweight#ppee-cross-section
>
> I'm (still) using python2.7, so I did
>
> python2.7 bin/mg5_aMC
> generate p p > e+ e-
> output Xcheck_rwgt
>
> You can find the output in [1].
>
> Afterwards, I created the reweight_card.dat in the Cards directory with
>
> launch --rwgt_name=case1
> set aewm1 100
> launch --rwgt_name=case2
> set aewm1 200
> launch --rwgt_name=case3
> set aewm1 300
>
>
> You can find the output that I see when I launch the event generation in [2]. As you can see, there is an import error indicating that it doesn't find a certain module.
>
> Do you have any idea what I did wrong?
>
> Kind regards,
> Philip
>
> [1]
> WARNING:root:python2 support will be removed in last quarter 2021. If you use python2 due to issue with Python3, please report them on https://bugs.launchpad.net/mg5amcnlo
> ************************************************************
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 3.1.0 2021-03-30 *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * and *
> * http://amcatnlo.web.cern.ch/amcatnlo/ *
> * *
> * Type 'help' for in-line help. *
> * Type 'tutorial' to learn how MG5 works *
> * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
> * Type 'tutorial MadLoop' to learn how MadLoop works *
> * *
> ************************************************************
> load MG5 configuration from input/mg5_configuration.txt
> fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
> Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
> MG5_aMC> set fastjet /PATH/TO/fastjet-config
>
> set lhapdf to lhapdf-config
> None does not seem to correspond to a valid lhapdf-config executable.
> Please set the 'lhapdf' variable to the (absolute) /PATH/TO/lhapdf-config (including lhapdf-config).
> Note that you can still compile and run aMC@NLO with the built-in PDFs
> MG5_aMC> set lhapdf /PATH/TO/lhapdf-config
>
> Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> Loading default model: sm
> INFO: Restrict model sm with file models/sm/restrict_default.dat .
> INFO: Run "set stdout_level DEBUG" before import for more information.
> INFO: Change particles name to pass to MG5 convention
> Defined multiparticle p = g u c d s u~ c~ d~ s~
> Defined multiparticle j = g u c d s u~ c~ d~ s~
> Defined multiparticle l+ = e+ mu+
> Defined multiparticle l- = e- mu-
> Defined multiparticle vl = ve vm vt
> Defined multiparticle vl~ = ve~ vm~ vt~
> Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
> MG5_aMC>generate p p > e+ e-
> INFO: Checking for minimal orders which gives processes.
> INFO: Please specify coupling orders to bypass this step.
> INFO: Trying process: g g > e+ e- WEIGHTED<=4 @1
> INFO: Trying process: u u~ > e+ e- WEIGHTED<=4 @1
> INFO: Process has 2 diagrams
> INFO: Trying process: u c~ > e+ e- WEIGHTED<=4 @1
> INFO: Trying process: c u~ > e+ e- WEIGHTED<=4 @1
> INFO: Trying process: c c~ > e+ e- WEIGHTED<=4 @1
> INFO: Process has 2 diagrams
> INFO: Trying process: d d~ > e+ e- WEIGHTED<=4 @1
> INFO: Process has 2 diagrams
> INFO: Trying process: d s~ > e+ e- WEIGHTED<=4 @1
> INFO: Trying process: s d~ > e+ e- WEIGHTED<=4 @1
> INFO: Trying process: s s~ > e+ e- WEIGHTED<=4 @1
> INFO: Process has 2 diagrams
> INFO: Process u~ u > e+ e- added to mirror process u u~ > e+ e-
> INFO: Process c~ c > e+ e- added to mirror process c c~ > e+ e-
> INFO: Process d~ d > e+ e- added to mirror process d d~ > e+ e-
> INFO: Process s~ s > e+ e- added to mirror process s s~ > e+ e-
> 4 processes with 8 diagrams generated in 0.046 s
> Total: 4 processes with 8 diagrams
> MG5_aMC>output Xcheck_rwgt
> INFO: initialize a new directory: Xcheck_rwgt
> INFO: remove old information in Xcheck_rwgt
> INFO: Organizing processes into subprocess groups
> INFO: Generating Helas calls for process: u u~ > e+ e- WEIGHTED<=4 @1
> INFO: Processing color information for process: u u~ > e+ e- @1
> INFO: Combined process c c~ > e+ e- WEIGHTED<=4 @1 with process u u~ > e+ e- WEIGHTED<=4 @1
> INFO: Generating Helas calls for process: d d~ > e+ e- WEIGHTED<=4 @1
> INFO: Reusing existing color information for process: d d~ > e+ e- @1
> INFO: Combined process s s~ > e+ e- WEIGHTED<=4 @1 with process d d~ > e+ e- WEIGHTED<=4 @1
> INFO: Creating files in directory P1_qq_ll
> INFO: Generating Feynman diagrams for Process: u u~ > e+ e- WEIGHTED<=4 @1
> INFO: Generating Feynman diagrams for Process: d d~ > e+ e- WEIGHTED<=4 @1
> INFO: Finding symmetric diagrams for subprocess group qq_ll
> Generated helas calls for 2 subprocesses (4 diagrams) in 0.013 s
> Wrote files for 16 helas calls in 0.283 s
> ALOHA: aloha starts to compute helicity amplitudes
> ALOHA: aloha creates 5 routines in 2.407 s
> save configuration file to /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt/Cards/me5_configuration.txt
> INFO: Use Fortran compiler gfortran
> INFO: Use c++ compiler g++
> INFO: Generate jpeg diagrams
> gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_19/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
> gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_19/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
> gs: /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_19/external/slc7_amd64_gcc820/lib/libtiff.so.5: no version information available (required by /lib64/libgs.so.9)
> INFO: Generate web pages
> Output to directory /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt done.
> Type "launch" to generate events from this process, or see
> /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt/README
> Run "open index.html" to see more information about this process.
>
> [2]
>
> python2.7 bin/madevent
> launch in debug mode
> ************************************************************
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * M A D E V E N T *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 3.1.0 20xx-xx-xx *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * *
> * Type 'help' for in-line help. *
> * *
> ************************************************************
> INFO: load configuration from /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt/Cards/me5_configuration.txt
> INFO: load configuration from /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/input/mg5_configuration.txt
> INFO: load configuration from /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt/Cards/me5_configuration.txt
> Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> Xcheck_rwgt>launch
> The following switches determine which programs are run:
> /=================== Description ===================|============= values ==============|======== other options ========\
> | 1. Choose the shower/hadronization program | shower = Not Avail. | Please install module |
> | 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
> | 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
> | 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
> | 5. Add weights to events for new hypp. | reweight = ON | OFF |
> \=======================================================================================================================/
> Either type the switch number (1 to 5) to change its setting,
> Set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type 'help' for the list of all valid option
> Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>>
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> Do you want to edit a card (press enter to bypass editing)?
> /------------------------------------------------------------\
> | 1. param : param_card.dat |
> | 2. run : run_card.dat |
> | 3. reweight : reweight_card.dat |
> \------------------------------------------------------------/
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/MadWidth/...).
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, 3, reweight, enter path][90s to answer]
>>
> INFO: Update the dependent parameter of the param_card.dat
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> Generating 100 events with run name run_01
> survey run_01
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> INFO: compile directory
> compile Source Directory
> Running without a11y support!
> Using random number seed offset = 21
> INFO: Running Survey
> Creating Jobs
> Working on SubProcesses
> INFO: P1_qq_ll
> INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 16h12 ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 1.2s ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 1.2s ]
> INFO: End survey
> refine 100
> Creating Jobs
> INFO: Refine results to 100
> INFO: Generating 100.0 unweighted events.
> sum of cpu time of last step: 1 seconds
> INFO: Effective Luminosity 0.158673489627 pb^-1
> INFO: need to improve 1 channels
> - Current estimate of cross-section: 756.27 +- 30.8942449333
> P1_qq_ll
> INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 16h12 ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 1.2s ]
> INFO: Combining runs
> sum of cpu time of last step: 1 seconds
> INFO: finish refine
> refine 100 --treshold=0.9
> No need for second refine due to stability of cross-section
> INFO: Combining Events
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 819.8 +- 16.64 pb
> Nb of events : 100
>
> store_events
> INFO: Storing parton level results
> INFO: End Parton
> reweight -from_cards
> INFO: Extracting the banner ...
> INFO: process: p p > e+ e-
> INFO: options:
> INFO: Running Reweighting
> change mode NLO # Define type of Reweighting. For LO sample this command
> Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
> In general this optimization speed up the computation be a factor of two.
> WARNING: Missing information for NLO type of reweighting. Keeping the mode to LO.
> launch --rwgt_name=case1
> INFO: detected model: sm. Loading...
> INFO: generating the square matrix element for reweighting
> INFO: generate p p > e+ e- ;
> INFO: Done 4.086
> Command "launch" interrupted with error:
> ImportError : No module named P1_uux_epem2py
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt/run_01_tag_1_debug.log'.
> Please attach this file to your report.
> INFO: Original cross-section: 819.82 +- 16.64272 pb
> INFO: Computed cross-section:
> Xcheck_rwgt>INFO: storing files of previous run
> gzipping output file: unweighted_events.lhe
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear Olivier,

thanks for your fast reply! Indeed, after I switched to python 3.8, everything worked fine.

Just one more question: I now see these values after calculating 10000 events

INFO: Original cross-section: 840.95 +- 2.79872 pb (cross-section from sum of weights: 840.9500000001293)
INFO: Computed cross-section:
INFO: case1 : 1008.4427639054032 +- 102.84109201449088 pb
INFO: case2 : 977.6378802503561 +- 78.46595620204488 pb
INFO: case3 : 1191.8647126398423 +- 132.07035995981005 pb

I compared them to the values on the Reweight page:

The associated cross sections are

    1135.25 pb
    1095.28 pb
    1329.52 pb

As you can see, my values are smaller than what is quoted on the Reweight documentation page. Is that a problem?

Thanks again!

Cheers,
Philip

Dear all,

I switched from LO to NLO and tried the reweighting again. Here [1] are the steps I did In MG. I modified the reweight card I used in the LO tests to match the parameter names at NLO [2].

When I launch aMC@NLO, it first looks fine (to me) [3]. But in the end, an error occurs in the reweighting. I appended the debug.log file in [4]. Do you have an idea where this is coming from?

Cheers,
Philip

[1]
```
generate p p > e+ ve [QCD]
output Xcheck_rwgt_NLO
```

[2]
```
change mode NLO # Define type of Reweighting. For LO sample this command
                   # has no effect since only "LO" mode is allowed.

launch --rwgt_name=case1
set aEWM1 100
launch --rwgt_name=case2
set aEWM1 200
launch --rwgt_name=case3
set aEWM1 300
```

[3]

launch aMC@NLO
INFO: will run in mode: aMC@NLO
INFO: modify parameter parton_shower of the run_card.dat to HERWIG6
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. reweight : reweight_card.dat |
| 4. shower : shower_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, reweight, 4, enter path, ... ][90s to answer]
>2
open /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Cards/run_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. reweight : reweight_card.dat |
| 4. shower : shower_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, reweight, 4, enter path, ... ]
>
INFO: NLO reweighting is on ON. Automatically set store_rwgt_info to True
INFO: modify parameter store_rwgt_info of the run_card.dat to True
INFO: modify parameter extralibs of the shower_card.dat to pythia8 boost_iostreams z dl stdc++
INFO: modify parameter extrapaths of the shower_card.dat to ../lib /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/pythia8//lib /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_19/external/slc7_amd64_gcc820/lib /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/zlib/lib
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 10 cores
INFO: Compiling P0_udx_veep...
INFO: Compiling P0_dxu_veep...
INFO: P0_udx_veep done.
INFO: P0_dxu_veep done.
INFO: Checking test output:
INFO: P0_udx_veep
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_dxu_veep
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 10 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grids
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 3s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 3.5s ]
sum of cpu time of last step: 0 second
INFO: Determining the number of unweighted events per channel

      Intermediate results:
      Random seed: 34
      Total cross section: 1.088e+04 +- 7.1e+01 pb
      Total abs(cross section): 1.253e+04 +- 7.3e+01 pb

INFO: Computing upper envelope
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 7.4s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 8.2s ]
sum of cpu time of last step: 0 second
INFO: Updating the number of unweighted events per channel

      Intermediate results:
      Random seed: 34
      Total cross section: 1.086e+04 +- 4.5e+01 pb
      Total abs(cross section): 1.267e+04 +- 4.3e+01 pb

INFO: Generating events
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 1.2s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 1.2s ]
sum of cpu time of last step: 0 second
INFO: Doing reweight
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.57s ]
INFO: Idle: 0, Running: 0, Completed: 2 [ 0.66s ]
INFO: Collecting events
INFO:
   --------------------------------------------------------------
      Summary:
      Process p p > e+ ve [QCD]
      Run at p-p collider (6500.0 + 6500.0 GeV)
      Number of events generated: 1000
      Total cross section: 1.086e+04 +- 4.5e+01 pb
   --------------------------------------------------------------

INFO: The /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz file has been generated.

INFO: Events generated
reweight -from_cards
INFO: split the event file in bunch of 2500 events
INFO: Idle: 1, Running: 0, Completed: 0 [ 0.025s ]
INFO: Idle: 0, Running: 1, Completed: 0 [ 5m 0s ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 6m 11s ]
INFO: Collect and combine the various output file.
Command "launch aMC@NLO" interrupted with error:
Exception : Some of the run failed: Please read /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz_0.lhe_tag_1_debug.log
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/run_02_tag_1_debug.log'.
Please attach this file to your report.

[4]

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.3.1.0 20xx-xx-xx *
#* *
#* 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 *
#* *
#************************************************************
reweight /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz_0.lhe -from_cards --multicore=create
Traceback (most recent call last):
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/extended_cmd.py", line 1544, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/common_run_interface.py", line 2184, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1692, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 564, in do_launch
    weight = self.calculate_weight(event)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 1003, in calculate_weight
    return self.calculate_nlo_weight(event)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 1046, in calculate_nlo_weight
    event.parse_nlo_weight(threshold=self.soft_threshold)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 1496, in parse_nlo_weight
    self.nloweight = NLO_PARTIALWEIGHT(text, self, real_type=real_type,
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 3024, in __init__
    self.parse(input)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 3072, in parse
    wgt = OneNLOWeight(line, real_type=self.real_type)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 2662, in __init__
    self.parse(input)
  File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 2767, in parse
    self.born_related = int(data[flag+7])
ValueError: invalid literal for int() with base 10: '0.13998036e+01'
Value of current Options:
             pythia8_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/pythia8
                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 : None
                   lhapdf : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/lhapdf6_py3/bin/lhapdf-config
               lhapdf_py2 : None
               lhapdf_py3 : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/lhapdf6_py3/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 : 7
             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 : True
      notification_center : True
                 run_mode : 2
            cluster_queue : None
             cluster_time : None
             cluster_size : 100
           cluster_memory : None
                  nb_core : 10
        cluster_temp_path : None
                 mg5_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0
mg5amc_py8_interface_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/MG5aMC_PY8_interface
                    golem : None
                    ninja : ./HEPTools/lib
                  collier : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/collier
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.1.0 2021-03-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 include_lepton_initiated_processes 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
set nlo_mixed_expansion 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 loop_sm
tutorial aMCatNLO
generate p p > e+ ve [QCD]
output Xcheck_rwgt_NLO
######################################################################
## 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
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   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.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 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 8.041900e+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.325070e+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
    5 4.700000e+00 # ymb
    6 1.730000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-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.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
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 82 # gh
        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)
#***********************************************************************
# 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 *
#***********************************************************************
  1000 = 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 *
#***********************************************************************
  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 *
# - 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
#***********************************************************************

NLO re-weighting is currently broken within 3.x branch.
You can use the LTS version (2.9.4) which supports such feature.

Cheers,

Olivier

> On 15 Jul 2021, at 10:31, Philip Keicher <email address hidden> wrote:
>
> Question #697979 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/697979
>
> Status: Solved => Open
>
> Philip Keicher is still having a problem:
> Dear all,
>
> I switched from LO to NLO and tried the reweighting again. Here [1] are
> the steps I did In MG. I modified the reweight card I used in the LO
> tests to match the parameter names at NLO [2].
>
> When I launch aMC@NLO, it first looks fine (to me) [3]. But in the end,
> an error occurs in the reweighting. I appended the debug.log file in
> [4]. Do you have an idea where this is coming from?
>
> Cheers,
> Philip
>
>
> [1]
> ```
> generate p p > e+ ve [QCD]
> output Xcheck_rwgt_NLO
> ```
>
> [2]
> ```
> change mode NLO # Define type of Reweighting. For LO sample this command
> # has no effect since only "LO" mode is allowed.
>
> launch --rwgt_name=case1
> set aEWM1 100
> launch --rwgt_name=case2
> set aEWM1 200
> launch --rwgt_name=case3
> set aEWM1 300
> ```
>
> [3]
>
> launch aMC@NLO
> INFO: will run in mode: aMC@NLO
> INFO: modify parameter parton_shower of the run_card.dat to HERWIG6
> Do you want to edit a card (press enter to bypass editing)?
> /------------------------------------------------------------\
> | 1. param : param_card.dat |
> | 2. run : run_card.dat |
> | 3. reweight : reweight_card.dat |
> | 4. shower : shower_card.dat |
> \------------------------------------------------------------/
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/MadWidth/...).
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, 3, reweight, 4, enter path, ... ][90s to answer]
>> 2
> open /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Cards/run_card.dat
> Do you want to edit a card (press enter to bypass editing)?
> /------------------------------------------------------------\
> | 1. param : param_card.dat |
> | 2. run : run_card.dat |
> | 3. reweight : reweight_card.dat |
> | 4. shower : shower_card.dat |
> \------------------------------------------------------------/
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/MadWidth/...).
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, 3, reweight, 4, enter path, ... ]
>>
> INFO: NLO reweighting is on ON. Automatically set store_rwgt_info to True
> INFO: modify parameter store_rwgt_info of the run_card.dat to True
> INFO: modify parameter extralibs of the shower_card.dat to pythia8 boost_iostreams z dl stdc++
> INFO: modify parameter extrapaths of the shower_card.dat to ../lib /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/pythia8//lib /cvmfs/cms.cern.ch/slc7_amd64_gcc820/cms/cmssw/CMSSW_10_6_19/external/slc7_amd64_gcc820/lib /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/zlib/lib
> INFO: Update the dependent parameter of the param_card.dat
> INFO: Starting run
> INFO: Compiling the code
> INFO: Using built-in libraries for PDFs
> INFO: Compiling source...
> INFO: ...done, continuing with P* directories
> INFO: Compiling directories...
> INFO: Compiling on 10 cores
> INFO: Compiling P0_udx_veep...
> INFO: Compiling P0_dxu_veep...
> INFO: P0_udx_veep done.
> INFO: P0_dxu_veep done.
> INFO: Checking test output:
> INFO: P0_udx_veep
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: P0_dxu_veep
> INFO: Result for test_ME:
> INFO: Passed.
> INFO: Result for test_MC:
> INFO: Passed.
> INFO: Result for check_poles:
> INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> INFO: Starting run
> INFO: Using 10 cores
> INFO: Cleaning previous results
> INFO: Doing NLO matched to parton shower
> INFO: Setting up grids
> INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
> INFO: Idle: 0, Running: 1, Completed: 1 [ 3s ]
> INFO: Idle: 0, Running: 0, Completed: 2 [ 3.5s ]
> sum of cpu time of last step: 0 second
> INFO: Determining the number of unweighted events per channel
>
> Intermediate results:
> Random seed: 34
> Total cross section: 1.088e+04 +- 7.1e+01 pb
> Total abs(cross section): 1.253e+04 +- 7.3e+01 pb
>
>
> INFO: Computing upper envelope
> INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
> INFO: Idle: 0, Running: 1, Completed: 1 [ 7.4s ]
> INFO: Idle: 0, Running: 0, Completed: 2 [ 8.2s ]
> sum of cpu time of last step: 0 second
> INFO: Updating the number of unweighted events per channel
>
> Intermediate results:
> Random seed: 34
> Total cross section: 1.086e+04 +- 4.5e+01 pb
> Total abs(cross section): 1.267e+04 +- 4.3e+01 pb
>
>
> INFO: Generating events
> INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
> INFO: Idle: 0, Running: 1, Completed: 1 [ 1.2s ]
> INFO: Idle: 0, Running: 0, Completed: 2 [ 1.2s ]
> sum of cpu time of last step: 0 second
> INFO: Doing reweight
> INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 10h15 ]
> INFO: Idle: 0, Running: 1, Completed: 1 [ 0.57s ]
> INFO: Idle: 0, Running: 0, Completed: 2 [ 0.66s ]
> INFO: Collecting events
> INFO:
> --------------------------------------------------------------
> Summary:
> Process p p > e+ ve [QCD]
> Run at p-p collider (6500.0 + 6500.0 GeV)
> Number of events generated: 1000
> Total cross section: 1.086e+04 +- 4.5e+01 pb
> --------------------------------------------------------------
>
> INFO: The /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz file has been generated.
>
> INFO: Events generated
> reweight -from_cards
> INFO: split the event file in bunch of 2500 events
> INFO: Idle: 1, Running: 0, Completed: 0 [ 0.025s ]
> INFO: Idle: 0, Running: 1, Completed: 0 [ 5m 0s ]
> INFO: Idle: 0, Running: 0, Completed: 1 [ 6m 11s ]
> INFO: Collect and combine the various output file.
> Command "launch aMC@NLO" interrupted with error:
> Exception : Some of the run failed: Please read /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz_0.lhe_tag_1_debug.log
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/run_02_tag_1_debug.log'.
> Please attach this file to your report.
>
>
> [4]
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.3.1.0 20xx-xx-xx *
> #* *
> #* 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 *
> #* *
> #************************************************************
> reweight /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/Events/run_02/events.lhe.gz_0.lhe -from_cards --multicore=create
> Traceback (most recent call last):
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/extended_cmd.py", line 1544, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/Xcheck_rwgt_NLO/bin/internal/common_run_interface.py", line 2184, in do_reweight
> reweight_cmd.import_command_file(path)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1692, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1573, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/extended_cmd.py", line 1493, in onecmd_orig
> return func(arg, **opt)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/misc.py", line 109, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 564, in do_launch
> weight = self.calculate_weight(event)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 1003, in calculate_weight
> return self.calculate_nlo_weight(event)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/interface/reweight_interface.py", line 1046, in calculate_nlo_weight
> event.parse_nlo_weight(threshold=self.soft_threshold)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 1496, in parse_nlo_weight
> self.nloweight = NLO_PARTIALWEIGHT(text, self, real_type=real_type,
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 3024, in __init__
> self.parse(input)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 3072, in parse
> wgt = OneNLOWeight(line, real_type=self.real_type)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 2662, in __init__
> self.parse(input)
> File "/afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/madgraph/various/lhe_parser.py", line 2767, in parse
> self.born_related = int(data[flag+7])
> ValueError: invalid literal for int() with base 10: '0.13998036e+01'
> Value of current Options:
> pythia8_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/pythia8
> 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 : None
> lhapdf : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/lhapdf6_py3/bin/lhapdf-config
> lhapdf_py2 : None
> lhapdf_py3 : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/lhapdf6_py3/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 : 7
> 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 : True
> notification_center : True
> run_mode : 2
> cluster_queue : None
> cluster_time : None
> cluster_size : 100
> cluster_memory : None
> nb_core : 10
> cluster_temp_path : None
> mg5_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0
> mg5amc_py8_interface_path : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/MG5aMC_PY8_interface
> golem : None
> ninja : ./HEPTools/lib
> collier : /afs/cern.ch/work/p/pkeicher/public/MG5_aMC_v3_1_0/HEPTools/collier
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 3.1.0 2021-03-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 include_lepton_initiated_processes 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
> set nlo_mixed_expansion 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 loop_sm
> tutorial aMCatNLO
> generate p p > e+ ve [QCD]
> output Xcheck_rwgt_NLO
> ######################################################################
> ## 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
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.118800e+01 # MZ
> 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.000000e+00 # d : 0.0
> 2 0.000000e+00 # u : 0.0
> 3 0.000000e+00 # s : 0.0
> 4 0.000000e+00 # c : 0.0
> 11 0.000000e+00 # e- : 0.0
> 12 0.000000e+00 # ve : 0.0
> 13 0.000000e+00 # mu- : 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 8.041900e+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.325070e+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
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-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.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
> #===========================================================
> # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
> #===========================================================
>
> Block QNUMBERS 82 # gh
> 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)
> #***********************************************************************
> # 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 *
> #***********************************************************************
> 1000 = 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 *
> #***********************************************************************
> 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 *
> # - 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
> #***********************************************************************
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear Olivier,

thanks for the message. I'll try something else then.

Cheers,
Philip

Given the size on your statistical error. This sounds reasonable

Cheers,

Olivier

> On 14 Jul 2021, at 17:31, Philip Keicher <email address hidden> wrote:
>
> Question #697979 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/697979
>
> Status: Answered => Solved
>
> Philip Keicher confirmed that the question is solved:
> Dear Olivier,
>
> thanks for your fast reply! Indeed, after I switched to python 3.8,
> everything worked fine.
>
> Just one more question: I now see these values after calculating 10000
> events
>
> INFO: Original cross-section: 840.95 +- 2.79872 pb (cross-section from sum of weights: 840.9500000001293)
> INFO: Computed cross-section:
> INFO: case1 : 1008.4427639054032 +- 102.84109201449088 pb
> INFO: case2 : 977.6378802503561 +- 78.46595620204488 pb
> INFO: case3 : 1191.8647126398423 +- 132.07035995981005 pb
>
>
> I compared them to the values on the Reweight page:
>
> The associated cross sections are
>
> 1135.25 pb
> 1095.28 pb
> 1329.52 pb
>
> As you can see, my values are smaller than what is quoted on the
> Reweight documentation page. Is that a problem.
>
> Thanks again!
>
> Cheers,
> Philip
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.