Trouble with reweighting while using coupling orders

Asked by Avik Roy on 2019-05-07

A little bit of context:

The process I am trying to simulate is single production of vector-like quarks and their subsequent decay via w/z/h boson. I am trying to use it witin ATLAS EDM and the athena cache that has a stable support for the reweighting module uses MG version 2.6.1. So far I have been using MG to produce the vlq and decay it to the boson and quark, as in the following process -

add process p p > j VLQ bb / tp tp~ p t t~ y y~ bp bp~ x x~ z h a, (VLQ > ww bb)

Reweighting has worked just fine for this process. With this job option, decays of z and t are taken care of by Pythia. Everything works just fine for this.

Now, the problem:

Instead of specifying the production of vlq in the process, I now only specify the final state I want and constrain the process using coupling orders, as in the following commands-

add process p p > j w+ b~ b / p t t~ tp tp~ bp bp~ x x~ z h a YWB^2==4
add process p p > j w- b b~ / p t t~ tp tp~ bp bp~ x x~ z h a YWB^2==4

Allowing such processes are interesting for a number of studies, including interference. The coupling orders should make sure that I have two vertices with Y in association with W and b, essentially regenerating the same set of Feynman diagrams as before. However, with this command, reweighting fails with the following set of errors-

generate 18:56:57 REWEIGHT: generate p p > j w+ b~ b / p t t~ tp tp~ bp bp~ x x~ z h a YWB^2==4 ;add process p p > j w- b b~ / p t t~ tp tp~ bp bp~ x x~ z h a YWB^2==4 ;
generate 18:56:57 REWEIGHT: Done 7.828
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M2_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M4_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M6_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M14_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M31_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M26_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M11_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M34_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M28_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M38_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M0_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M3_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M5_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M7_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M24_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M15_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M30_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M35_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M27_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M25_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M29_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M10_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M39_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]
generate 18:56:57 CRITICAL: pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper() ) = /afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/rw_me/SubProcesses/MadLoop5_resources/M1_HelConfigs.dat [reweight_interface.py at line 1731]
generate 18:56:57 CRITICAL: os.path.exists(pjoin(path_me,onedir,'SubProcesses','MadLoop5_resources', '%sHelConfigs.dat' % prefix.upper())) = False [reweight_interface.py at line 1732]

The reweight card is as follows-
change mode LO
launch --rwgt_name=M10K05BR040303
 set KY 0.5
 set MY 1000
 set WY auto
launch --rwgt_name=M10K15BR040303
 set KY 1.5
 set MY 1000
 set WY auto
launch --rwgt_name=M11K05BR040303
 set KY 0.5
 set MY 1100
 set WY auto
launch --rwgt_name=M11K15BR040303
 set KY 1.5
 set MY 1100
 set WY auto

It is the exact same reweight card that does work for the first process I mentioned.

From the log file, I get the following error:

reweight run_01 -from_cards
Traceback (most recent call last):
  File "/afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/bin/internal/extended_cmd.py", line 1467, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/bin/internal/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/a/avroy/MG_Reweight_testinATLAS2/run_test_Reweight__10_madspin/PROC_VLQ_UFO_for_interference_0/bin/internal/common_run_interface.py", line 2035, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1609, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1494, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/extended_cmd.py", line 1421, in onecmd_orig
    return func(arg, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/reweight_interface.py", line 480, in do_launch
    self.load_module()
  File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.1.atlas/x86_64-slc6-gcc47-opt/madgraph/interface/reweight_interface.py", line 1749, in load_module
    hel = hel_dict[prefix]
KeyError: 'm4_'

Please note that the VLQ UFO model we use is a LO model, so we are not looking for any loop induced diagrams.

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2019-05-07
Last reply:
2019-05-07

Can you try with 2.6.5?

If this does not work, can you give me the full banner for the generated sample and potentially your model if it is not part of the list available in the command "display modellist" (in short everything I need to reproduce the problem).

Note that I do not know Athena (not even sure that I'm allowed to have access to it)

Cheers,

Olivier

Avik Roy (avroy) said : #2

Hi Olivier

Trying with 2.6.5 worked. I guess it is rather a version dependent issue. We will probably need a version of athena that supports 2.6.5, but that's a question for a different forum I suppose.
Thanks for your suggestion.

Regards,
Avik

Can you help with this problem?

Provide an answer of your own, or ask Avik Roy for more information if necessary.

To post a message you must log in.