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

When	Who	Action	New State
2021-07-14	Olivier Mattelaer	Answer	Answered
2021-07-14	Philip Keicher	Confirm	Solved
2021-07-15	Philip Keicher	Reopen	Open
2021-07-15	Olivier Mattelaer	Answer	Answered
2021-07-15	Philip Keicher	Comment	Answered
2021-07-16	Olivier Mattelaer	Answer	Answered

MadGraph5_aMC@NLO

History of question #697979