Reweight pilotrun AssertionError
Dear experts,
I am using MG5 v2.6.5 together with SMEFTsim_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme False
define p = g u c d s b u~ c~ d~ s~ b~
define j = p
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model SMEFTsim_
generate p p > z w- j j QCD=0 SMHLOOP=0 NP=1, z > l+ l- NP=1, w- > j j NP=1 @0 NP=1
output ZWm_cHbox_
and by using additional customizecards and reweight_card. During the reweighting step I get this error:
Command "reweight pilotrun" interrupted with error:
AssertionError :
Please report this bug on https:/
More information is found in '/home/
Please attach this file to your report.
Here is the debug file (I have included the error part only, let me know if you need more details):
reweight pilotrun
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
reweight_
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
out = f(self, *args, **opt)
File "/home/
self.
File "/home/
out = f(self, *args, **opt)
File "/home/
has_nlo = self.create_
File "/home/
out = f(self, *args, **opt)
File "/home/
mgcmd.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
self.
File "/home/
return self.cmd.
File "/home/
me, self._curr_
File "/home/
proc_
File "/home/
matrix_element)
File "/home/
matrix_
File "/home/
assert wfin.get('number') in first.values()
AssertionError
automatic_
cluster_
I made some additional attempts:
1) generate a simpler process
2) generate the same process described above, but decaying vector bosons in Madspin
Both 1) and 2) reweighting steps succeded, but I would be very grateful if I could manage to generate the VBS process without Madspin, in order to include possible "New Physics" vertices both in bosons production and decay.
Thanks for your support.
Kind regards,
Riccardo
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- 2021-01-07
- Last query:
- 2021-01-07
- Last reply:
- 2021-01-07
hi,
- What is customize_card? Are you using a plugin version of MG5aMC, we do not have such cards in MG5aMC.
If you are using a plugin, could you contact the author of such plugin or try without such plugin?
- Your model is unknow from MG5aMC, could you put a link to the associate tarball?
- Did you try with 2.8.2? 2.6.5 is very old and that your issue might have been fixed in the meantime.
Cheers,
Olivier
| Riccardo (ricbrs) said : | #2 |
Dear Olivier,
- The "customizecards" just set the default value for frblock parameters, so the ones used to compute the nominal weights.
- You can find the model " SMEFTsim_
(I renamed it "_v3_2" just because I downloaded it and added some restrict_cards)
- I am going to try with another version, but if there is a way to let 2.6.5 properly work, it would be better.
Cheers,
Riccardo
Ok thanks,
I was succesfull to reproduce the error within our development version.
So I should be able to provide a patch for the lastest version.
Providing a patch compatible with 2.6.5 is quite unlikely since it diverge that much from the current version and sounds a quite deep issue.
Cheers,
Olivier
> On 6 Jan 2021, at 16:15, Riccardo <email address hidden> wrote:
>
> Question #694807 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Riccardo is still having a problem:
> Dear Olivier,
>
> - The "customizecards" just set the default value for frblock
> parameters, so the ones used to compute the nominal weights.
>
> - You can find the model " SMEFTsim_
> (I renamed it "_v3_2" just because I downloaded it and added some restrict_cards)
>
> - I am going to try with another version, but if there is a way to let
> 2.6.5 properly work, it would be better.
>
> Cheers,
>
> Riccardo
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Here is the patch:
https:/
I guess that this should also work for 2.6.5 but I did not test.
Cheers,
Olivier
PS: Note that re-weighting is know to not work with such syntax (overall coupling order).
This is set to be fixed in the next version of the code.
| Riccardo (ricbrs) said : | #5 |
Dear Olivier,
I replaced the "old" helas_objects.py with the new one you provided, but now I have this error:
Command "import /home/public/
"output ZWm_cHbox_
PhysicsObjectError : polarization is not a valid property for this object: Leg
Valid property are ['number', 'onshell', 'state', 'from_group', 'loop_line', 'id']
Please report this bug on https:/
More information is found in 'MG5_debug'.
Please attach this file to your report.
Below you can find the MG5_debug.
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set loop_color_flows False
set gauge unitary
set max_npoint_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme False
define p = g u c d s b u~ c~ d~ s~ b~
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 p = g u c d s b u~ c~ d~ s~ b~
define j = p
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model SMEFTsim_
define p = 21 2 4 1 3 -2 -4 -1 -3 # pass to 4 flavors
define j = p
generate p p > z w- j j QCD=0 SMHLOOP=0 NP=1, z > l+ l- NP=1, w- > j j NP=1 @0 NP=1
output ZWm_cHbox_
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
self.
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
self.
File "/home/
return self.cmd.
File "/home/
ndiags, cpu_time = generate_
File "/home/
subproc_
File "/home/
self.
File "/home/
tmp_
File "/home/
self.
File "/home/
decay_ids)
File "/home/
**matrix_
File "/home/
self.
File "/home/
for leg in process.
File "/home/
if 99 in leg.get(
File "/home/
return self[name]
File "/home/
self.
File "/home/
Valid property are %s""" % (name,self.
PhysicsObjectError: polarization is not a valid property for this object: Leg
Valid property are ['number', 'onshell', 'state', 'from_group', 'loop_line', 'id']
default_
ignore_
loop_
low_mem_
max_
automatic_
cluster_
mg5amc_
Many thanks,
Riccardo
No applying the patch might work, but the file has diverged a lot since 2.6.5
Cheers,
Olivier
> On 7 Jan 2021, at 10:35, Riccardo <email address hidden> wrote:
>
> Question #694807 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Riccardo is still having a problem:
> Dear Olivier,
>
> I replaced the "old" helas_objects.py with the new one you provided, but
> now I have this error:
>
> Command "import /home/public/
> "output ZWm_cHbox_
> PhysicsObjectError : polarization is not a valid property for this object: Leg
>
> Valid property are ['number', 'onshell', 'state', 'from_group', 'loop_line', 'id']
> Please report this bug on https:/
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> Below you can find the MG5_debug.
>
> #******
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.5 2018-02-03 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https:/
> #* *
> #******
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #******
> set default_
> set loop_color_flows False
> set gauge unitary
> set max_npoint_
> set group_subprocesses Auto
> set ignore_
> set loop_optimized_
> set complex_mass_scheme False
> define p = g u c d s b u~ c~ d~ s~ b~
> 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 p = g u c d s b u~ c~ d~ s~ b~
> define j = p
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model SMEFTsim_
> define p = 21 2 4 1 3 -2 -4 -1 -3 # pass to 4 flavors
> define j = p
> generate p p > z w- j j QCD=0 SMHLOOP=0 NP=1, z > l+ l- NP=1, w- > j j NP=1 @0 NP=1
> output ZWm_cHbox_
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.cmd.
> File "/home/
> self.import_
> File "/home/
> self.exec_cmd(line, precmd=True)
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.cmd.
> File "/home/
> self.export(nojpeg, main_file_name, group_processes, args)
> File "/home/
> return self.cmd.
> File "/home/
> ndiags, cpu_time = generate_
> File "/home/
> subproc_
> File "/home/
> self.get(
> File "/home/
> tmp_matrix_
> File "/home/
> self.generate_
> File "/home/
> decay_ids)
> File "/home/
> **matrix_
> File "/home/
> self.generate_
> File "/home/
> for leg in process.
> File "/home/
> if 99 in leg.get(
> File "/home/
> return self[name]
> File "/home/
> self.is_
> File "/home/
> Valid property are %s""" % (name,self.
> PhysicsObjectError: polarization is not a valid property for this object: Leg
>
> Valid property are ['number', 'onshell', 'state', 'from_group', 'loop_line', 'id']
> MadGraph5_aMC@NLO Options
> ----------------
> complex_mass_scheme : False
> default_
> gauge : unitary
> group_subprocesses : Auto
> ignore_
> loop_color_flows : False
> loop_optimized_
> low_mem_
> max_npoint_
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_
> nb_core : None
> notification_center : True
> run_mode : 2
>
> Configuration Options
> -------
> OLP : MadLoop
> amcfast : amcfast-config
> applgrid : applgrid-config
> auto_update : 0 (user set)
> cluster_local_path : /cvmfs/
> cluster_nb_retry : 1
> cluster_queue : None (user set)
> cluster_retry_wait : 300
> cluster_size : 100
> cluster_
> cluster_temp_path : None
> cluster_type : condor
> collier : ./HEPTools/lib
> cpp_compiler : None
> crash_on_error : False
> delphes_path : ./Delphes
> eps_viewer : None
> exrootanalysis_path : ./ExRootAnalysis
> f2py_compiler : None
> fastjet : None (user set)
> fortran_compiler : None
> golem : None (user set)
> hepmc_path : None (user set)
> hwpp_path : None (user set)
> lhapdf : /cvmfs/
> madanalysis5_path : None (user set)
> madanalysis_path : ./MadAnalysis
> mg5amc_
> ninja : ./HEPTools/lib
> output_dependencies : external
> pjfry : None (user set)
> pythia-pgs_path : ./pythia-pgs
> pythia8_path : None (user set)
> samurai : None
> syscalc_path : ./SysCalc
> td_path : ./td
> text_editor : None
> thepeg_path : None (user set)
> timeout : 60
> web_browser : None
>
> Many thanks,
>
> Riccardo
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Riccardo (ricbrs) said : | #7 |
Dear Olivier,
By applying the patch, I have the following error:
INFO: Creating files in directory P0_qq_zwmqq_
Command "import /afs/cern.
"output ZWm_cHbox_
KeyError : 6
Please report this bug on https:/
More information is found in 'MG5_debug'.
Please attach this file to your report.
MG5_debug file:
output ZWm_cHbox_
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
self.
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
self.
File "/home/
return self.cmd.
File "/home/
self.
File "/home/
me_group, fortran_model, group_number)
File "/home/
subproc_
File "/home/
matrix_element)
File "/home/
matrix_
File "/home/
outdated.
KeyError: 6
default_
ignore_
loop_
low_mem_
max_
automatic_
cluster_
mg5amc_
Can you try with 2.8.2?
i'm worried that this might be linked to a python version ordering issue.
In that case using python3 will certainly help since the ordering has been ensure in python3 compare to python2. (This being said i have tested this with python2 and 2.8.2, will try to doulble check this with python3)
Olivier
> On 7 Jan 2021, at 14:01, Riccardo <email address hidden> wrote:
>
> Question #694807 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> Riccardo is still having a problem:
> Dear Olivier,
>
> By applying the patch, I have the following error:
>
> INFO: Creating files in directory P0_qq_zwmqq_
> Command "import /afs/cern.
> "output ZWm_cHbox_
> KeyError : 6
> Please report this bug on https:/
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> MG5_debug file:
>
> output ZWm_cHbox_
> Traceback (most recent call last):
> File "/home/
> return self.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.cmd.
> File "/home/
> self.import_
> File "/home/
> self.exec_cmd(line, precmd=True)
> File "/home/
> stop = Cmd.onecmd_
> File "/home/
> return func(arg, **opt)
> File "/home/
> self.cmd.
> File "/home/
> self.export(nojpeg, main_file_name, group_processes, args)
> File "/home/
> return self.cmd.
> File "/home/
> self._curr_
> File "/home/
> me_group, fortran_model, group_number)
> File "/home/
> subproc_
> File "/home/
> matrix_element)
> File "/home/
> matrix_
> File "/home/
> outdated.
> KeyError: 6
> MadGraph5_aMC@NLO Options
> ----------------
> complex_mass_scheme : False
> default_
> gauge : unitary
> group_subprocesses : Auto
> ignore_
> loop_color_flows : False
> loop_optimized_
> low_mem_
> max_npoint_
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_
> nb_core : None
> notification_center : True
> run_mode : 2
>
> Configuration Options
> -------
> OLP : MadLoop
> amcfast : amcfast-config
> applgrid : applgrid-config
> auto_update : 0 (user set)
> cluster_local_path : /cvmfs/
> cluster_nb_retry : 1
> cluster_queue : None (user set)
> cluster_retry_wait : 300
> cluster_size : 100
> cluster_
> cluster_temp_path : None
> cluster_type : condor
> collier : ./HEPTools/lib
> cpp_compiler : None
> crash_on_error : False
> delphes_path : ./Delphes
> eps_viewer : None
> exrootanalysis_path : ./ExRootAnalysis
> f2py_compiler : None
> fastjet : None (user set)
> fortran_compiler : None
> golem : None (user set)
> hepmc_path : None (user set)
> hwpp_path : None (user set)
> lhapdf : /cvmfs/
> madanalysis5_path : None (user set)
> madanalysis_path : ./MadAnalysis
> mg5amc_
> ninja : ./HEPTools/lib
> output_dependencies : external
> pjfry : None (user set)
> pythia-pgs_path : ./pythia-pgs
> pythia8_path : None (user set)
> samurai : None
> syscalc_path : ./SysCalc
> td_path : ./td
> text_editor : None
> thepeg_path : None (user set)
> timeout : 60
> web_browser : None
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
| Riccardo (ricbrs) said : | #9 |
Dear Olivier,
I made a mistake applying the patch, now it works.
Thank you for your support.
Cheers,
Riccardo
| Riccardo (ricbrs) said : | #10 |
Thanks Olivier Mattelaer, that solved my question.
