Skip to content

Commit

Permalink
[smeftfr] Add missing blocks to card
Browse files Browse the repository at this point in the history
  • Loading branch information
DavidMStraub committed Jul 8, 2019
1 parent 1265a28 commit 4bb1031
Show file tree
Hide file tree
Showing 2 changed files with 13 additions and 4 deletions.
2 changes: 1 addition & 1 deletion wcxf/cli.py 100755 → 100644
Expand Up @@ -239,5 +239,5 @@ def wcxf2smeftfr():
info_dict = smeftfr.analyze_smeftfr_param_card(card)
wc = wcxf.WC.load(args.FILE)
wc.validate()
smeftfr.wcxf2smeftfr(wc, info_dict, stream=args.output)
smeftfr.wcxf2smeftfr(wc, info_dict, card, stream=args.output)
return 0
15 changes: 12 additions & 3 deletions wcxf/converters/smeftfr.py
Expand Up @@ -62,7 +62,7 @@ def wcxf2smeftfr_dict(wc, info_dict):
return card


def wcxf2smeftfr(wc, info_dict, stream):
def wcxf2smeftfr(wc, info_dict, card_dict, stream):
"""Take a WC instance and dump it into an LHA file (or return as string)
in a format suitable for use as `param_card.dat` file for MadGraph
with SmeftFR.
Expand All @@ -83,5 +83,14 @@ def wcxf2smeftfr(wc, info_dict, stream):
wc_m = wc.translate('Warsaw mass')
except ValueError:
raise ValueError("wcxf2smeftfr requires a basis that can be translated to the Warsaw mass basis")
return pylha.dump({'BLOCK': wcxf2smeftfr_dict(wc_m, info_dict)},
fmt='lha', stream=stream)
card = {}
card['DECAY'] = card_dict['DECAY'].copy()
card['BLOCK'] = wcxf2smeftfr_dict(wc_m, info_dict)
for name, block in card_dict['BLOCK'].items():
if name not in card['BLOCK']:
card['BLOCK'][name] = card_dict['BLOCK'][name].copy()
# else:
# for k, v in block.items():
# if k not in card['BLOCK'][name]:
# card['BLOCK'][name][k] = v
return pylha.dump(card, fmt='lha', stream=stream)

0 comments on commit 4bb1031

Please sign in to comment.