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.