MadGraph5_aMC@NLO

Generation of cpg and other EFT operators in one single reweight card

Asked by Zef Wolffs

Dear Madgraph experts,

I am trying to generate reweighting of a ggF sample with some EFT operators enabled through the SMEFTatNLO model, however I believe this question is not model-specific. I am relying on the reweighting module of madgraph for this. Now, cpg is LO in SMEFT and the other operators (and SM) in ggF are thus all NLO due to the quark loop. This means that it is not possible to directly generate cpg x SM or cpg x <other operators> directly, but it can be done through reweighting as described here:

https://feynrules.irmp.ucl.ac.be/attachment/wiki/SMEFTatNLO/tree-loop_instructions.txt

So it can be done by reweighting any sample to (for SM x cpg):

```change process p p > h NP=2 QCD=0 QED=1 QCD^2==2 NP^2==2 [virt=QCD]
     launch --rwgt_name=SMxcpG
         set cpG 1
```

Now for the other samples, I reweight them as well, but to the loop-induced process, so by (for example for ctp):

```change process p p > h NP=2 QCD=0 QED=1 QCD^2==2 NP^2==2 [noborn=QCD]
     launch --rwgt_name=SMxcpt
         set cpt 1
```

What as I was wondering was if it is possible to do both of these in a single reweight card. Essentially this would require either being able to specify two "change process" commands in a reweight card, or if the "noborn=QCD" and "virt=QCD" parts are not strictly necessary I guess I could do without those parts in the process definition, and then I would not have to do "change process ...", but in that case I don't know if I can reweight to LO and NLO in one reweight card.

Any help would be much appreciated, thanks!

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

You can have multiple "change process" in the reweight_card as far as I remember.
Did you try?

Cheers,

Olivier

Revision history for this message
Zef Wolffs (zeff020) said : 		#2

I did try and it does not seem to work. It runs the reweighting successfully for every "launch" command after the first invocation of "change process" until the second invocation of "change process". It doesn't run reweighting for any of the launch commands after the second invocation of "change process".

The relevant part of the logs are below, note that some of the logs are ATLAS-specific, as I use MadGraph through an ATLAS software layer:

```
08:16:00 # content of reweight_card.dat ============
08:16:00
08:16:00 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @0 $$ z
08:16:00 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==2 [noborn=QCD] @1 $$ z --add
08:16:00 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @2 $$ z --add
08:16:00 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==2 [noborn=QCD] @3 $$ z --add
08:16:00 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @4 $$ z --add
08:16:00 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==2 [noborn=QCD] @5 $$ z --add
08:16:00
08:16:00 launch --rwgt_name=ctp --rwgt_info=ctp
08:16:00 set dim62f 19 1
08:16:00
08:16:00 launch --rwgt_name=ctG --rwgt_info=ctG
08:16:00 set dim62f 24 1
08:16:00
08:16:00 launch --rwgt_name=cpW --rwgt_info=cpW
08:16:00 set dim6 9 1
08:16:00
08:16:00 launch --rwgt_name=cpDc --rwgt_info=cpDc
08:16:00 set dim6 2 1
08:16:00
08:16:00 launch --rwgt_name=cdp --rwgt_info=cdp
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=c3pl1 --rwgt_info=c3pl1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=c3pl2 --rwgt_info=c3pl2
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=cll1221 --rwgt_info=cll1221
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=cp --rwgt_info=cp
08:16:00 set dim6 5 1
08:16:00
08:16:00
08:16:00
08:16:00
08:16:00 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @0 $$ z
08:16:00 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==4 [noborn=QCD] @1 $$ z --add
08:16:00 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @2 $$ z --add
08:16:00 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==4 [noborn=QCD] @3 $$ z --add
08:16:00 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @4 $$ z --add
08:16:00 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==4 [noborn=QCD] @5 $$ z --add
08:16:00
08:16:00
08:16:00
08:16:00
08:16:00 launch --rwgt_name=ctp_squared --rwgt_info=ctp_squared
08:16:00 set dim62f 19 1
08:16:00
08:16:00 launch --rwgt_name=ctG_squared --rwgt_info=ctG_squared
08:16:00 set dim62f 24 1
08:16:00
08:16:00 launch --rwgt_name=cpW_squared --rwgt_info=cpW_squared
08:16:00 set dim6 9 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_squared --rwgt_info=cpDc_squared
08:16:00 set dim6 2 1
08:16:00
08:16:00 launch --rwgt_name=cdp_squared --rwgt_info=cdp_squared
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=c3pl1_squared --rwgt_info=c3pl1_squared
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=c3pl2_squared --rwgt_info=c3pl2_squared
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=cll1221_squared --rwgt_info=cll1221_squared
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=cp_squared --rwgt_info=cp_squared
08:16:00 set dim6 5 1
08:16:00
08:16:00
08:16:00
08:16:00
08:16:00 launch --rwgt_name=ctp_ctG --rwgt_info=ctp_ctG
08:16:00 set dim62f 19 1
08:16:00 set dim62f 24 1
08:16:00
08:16:00 launch --rwgt_name=ctp_cpW --rwgt_info=ctp_cpW
08:16:00 set dim62f 19 1
08:16:00 set dim6 9 1
08:16:00
08:16:00 launch --rwgt_name=ctp_cpDc --rwgt_info=ctp_cpDc
08:16:00 set dim62f 19 1
08:16:00 set dim6 2 1
08:16:00
08:16:00 launch --rwgt_name=ctp_cdp --rwgt_info=ctp_cdp
08:16:00 set dim62f 19 1
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=ctp_c3pl1 --rwgt_info=ctp_c3pl1
08:16:00 set dim62f 19 1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=ctp_c3pl2 --rwgt_info=ctp_c3pl2
08:16:00 set dim62f 19 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=ctp_cll1221 --rwgt_info=ctp_cll1221
08:16:00 set dim62f 19 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=ctp_cp --rwgt_info=ctp_cp
08:16:00 set dim62f 19 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=ctG_cpW --rwgt_info=ctG_cpW
08:16:00 set dim62f 24 1
08:16:00 set dim6 9 1
08:16:00
08:16:00 launch --rwgt_name=ctG_cpDc --rwgt_info=ctG_cpDc
08:16:00 set dim62f 24 1
08:16:00 set dim6 2 1
08:16:00
08:16:00 launch --rwgt_name=ctG_cdp --rwgt_info=ctG_cdp
08:16:00 set dim62f 24 1
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=ctG_c3pl1 --rwgt_info=ctG_c3pl1
08:16:00 set dim62f 24 1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=ctG_c3pl2 --rwgt_info=ctG_c3pl2
08:16:00 set dim62f 24 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=ctG_cll1221 --rwgt_info=ctG_cll1221
08:16:00 set dim62f 24 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=ctG_cp --rwgt_info=ctG_cp
08:16:00 set dim62f 24 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=cpW_cpDc --rwgt_info=cpW_cpDc
08:16:00 set dim6 9 1
08:16:00 set dim6 2 1
08:16:00
08:16:00 launch --rwgt_name=cpW_cdp --rwgt_info=cpW_cdp
08:16:00 set dim6 9 1
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=cpW_c3pl1 --rwgt_info=cpW_c3pl1
08:16:00 set dim6 9 1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=cpW_c3pl2 --rwgt_info=cpW_c3pl2
08:16:00 set dim6 9 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=cpW_cll1221 --rwgt_info=cpW_cll1221
08:16:00 set dim6 9 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=cpW_cp --rwgt_info=cpW_cp
08:16:00 set dim6 9 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_cdp --rwgt_info=cpDc_cdp
08:16:00 set dim6 2 1
08:16:00 set dim6 4 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_c3pl1 --rwgt_info=cpDc_c3pl1
08:16:00 set dim6 2 1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_c3pl2 --rwgt_info=cpDc_c3pl2
08:16:00 set dim6 2 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_cll1221 --rwgt_info=cpDc_cll1221
08:16:00 set dim6 2 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=cpDc_cp --rwgt_info=cpDc_cp
08:16:00 set dim6 2 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=cdp_c3pl1 --rwgt_info=cdp_c3pl1
08:16:00 set dim6 4 1
08:16:00 set dim62f 4 1
08:16:00
08:16:00 launch --rwgt_name=cdp_c3pl2 --rwgt_info=cdp_c3pl2
08:16:00 set dim6 4 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=cdp_cll1221 --rwgt_info=cdp_cll1221
08:16:00 set dim6 4 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=cdp_cp --rwgt_info=cdp_cp
08:16:00 set dim6 4 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=c3pl1_c3pl2 --rwgt_info=c3pl1_c3pl2
08:16:00 set dim62f 4 1
08:16:00 set dim62f 5 1
08:16:00
08:16:00 launch --rwgt_name=c3pl1_cll1221 --rwgt_info=c3pl1_cll1221
08:16:00 set dim62f 4 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=c3pl1_cp --rwgt_info=c3pl1_cp
08:16:00 set dim62f 4 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=c3pl2_cll1221 --rwgt_info=c3pl2_cll1221
08:16:00 set dim62f 5 1
08:16:00 set dim64f4l 7 1
08:16:00
08:16:00 launch --rwgt_name=c3pl2_cp --rwgt_info=c3pl2_cp
08:16:00 set dim62f 5 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 launch --rwgt_name=cll1221_cp --rwgt_info=cll1221_cp
08:16:00 set dim64f4l 7 1
08:16:00 set dim6 5 1
08:16:00
08:16:00 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==2 [virt=QCD] @0 $$ z
08:16:00 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==2 [virt=QCD] @1 $$ z --add
08:16:00 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==2 [virt=QCD] @2 $$ z --add
08:16:00 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==2 [virt=QCD] @3 $$ z --add
08:16:00 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==2 [virt=QCD] @4 $$ z --add
08:16:00 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==2 [virt=QCD] @5 $$ z --add
08:16:00
08:16:00 launch --rwgt_name=cpG --rwgt_info=cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_renorscaled.dat
08:16:00
08:16:00 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==4 [virt=QCD] @0 $$ z
08:16:00 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==4 [virt=QCD] @1 $$ z --add
08:16:00 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==4 [virt=QCD] @2 $$ z --add
08:16:00 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==4 [virt=QCD] @3 $$ z --add
08:16:00 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==4 [virt=QCD] @4 $$ z --add
08:16:00 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==4 [virt=QCD] @5 $$ z --add
08:16:00
08:16:00 launch --rwgt_name=cpG_squared --rwgt_info=cpG_squared
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=ctp_cpG --rwgt_info=ctp_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_ctp_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=ctG_cpG --rwgt_info=ctG_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_ctG_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=cpW_cpG --rwgt_info=cpW_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_cpW_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=cpDc_cpG --rwgt_info=cpDc_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_cpDc_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=cdp_cpG --rwgt_info=cdp_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_cdp_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=c3pl1_cpG --rwgt_info=c3pl1_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_c3pl1_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=c3pl2_cpG --rwgt_info=c3pl2_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_c3pl2_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=cll1221_cpG --rwgt_info=cll1221_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_cll1221_renorscaled.dat
08:16:00
08:16:00 launch --rwgt_name=cp_cpG --rwgt_info=cp_cpG
08:16:00 /project/atlas/users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom/param_cards/param_card_cpG_cp_renorscaled.dat
08:16:00
08:16:00 # ====================================
08:16:00 Py:MadGraphUtils INFO Found madevent/bin/gridrun, starting generation.
08:16:00 Py:MadGraphUtils INFO Now generating 60 events with random seed 32289 and granularity 1
08:16:03 No handlers could be found for logger "madevent.cards"
08:16:03 P0_gg_llvlvl
08:16:04 P1_gg_llvlvl
08:16:05 P2_gg_llvlvl
08:16:06 P3_gg_llvlvl
08:16:07 P4_gg_llvlvl
08:16:07 P5_gg_llvlvl
08:16:44 WRITE GRIDCARD /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent
08:16:44 DONE
08:16:44 Py:MadGraphUtils INFO Running reweighting module on existing events
08:16:45 ************************************************************
08:16:45 * *
08:16:45 * W E L C O M E to *
08:16:45 * M A D G R A P H 5 _ a M C @ N L O *
08:16:45 * M A D E V E N T *
08:16:45 * *
08:16:45 * * * *
08:16:45 * * * * * *
08:16:45 * * * * * 5 * * * * *
08:16:45 * * * * * *
08:16:45 * * * *
08:16:45 * *
08:16:45 * VERSION 2.9.3 20xx-xx-xx *
08:16:45 * *
08:16:45 * The MadGraph5_aMC@NLO Development Team - Find us at *
08:16:45 * https://server06.fynu.ucl.ac.be/projects/madgraph *
08:16:45 * *
08:16:45 * Type 'help' for in-line help. *
08:16:45 * *
08:16:45 ************************************************************
08:16:45 INFO: load configuration from /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Cards/me5_configuration.txt
08:16:45 INFO: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/input/mg5_configuration.txt
08:16:45 INFO: load configuration from /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Cards/me5_configuration.txt
08:16:45 INFO: set f2py compiler to f2py2
08:16:45 INFO: set f2py compiler to f2py2
08:16:45 Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
08:16:45 No valid eps viewer found. Please set in ./input/mg5_configuration.txt
08:16:45 No valid web browser found. Please set in ./input/mg5_configuration.txt
08:16:45 Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
08:16:45 In general this optimization speed up the computation be a factor of two.
08:16:45 reweight GridRun_32289 -f
08:16:45 Do you want to edit a card (press enter to bypass editing)?
08:16:45 /------------------------------------------------------------\
08:16:45 | 1. reweight : reweight_card.dat |
08:16:45 \------------------------------------------------------------/
08:16:45 you can also
08:16:45 - enter the path to a valid card.
08:16:45 [0, done, 1, reweight, enter path][90s to answer]
08:16:45 >Helicity recycling optimization requires Python3. This optimzation is therefore deactivated automatically.
08:16:45 In general this optimization speed up the computation be a factor of two.
08:16:46 INFO: Extracting the banner ...
08:16:46 INFO: process: g g > e+ mu- ve vm~ QCD=2 QED=4 NP=0
08:16:46 INFO: options: $$ z
08:16:46 INFO: Running Reweighting
08:16:46 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @0 $$ z
08:16:46 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==2 [noborn=QCD] @1 $$ z --add
08:16:46 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @2 $$ z --add
08:16:46 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==2 [noborn=QCD] @3 $$ z --add
08:16:46 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==2 [noborn=QCD] @4 $$ z --add
08:16:46 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==2 [noborn=QCD] @5 $$ z --add
08:16:46 launch --rwgt_name=ctp --rwgt_info=ctp
08:16:46 INFO: detected model: /project/atlas/Users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom-NLO_renorscaled. Loading...
08:16:58 INFO: generating the square matrix element for reweighting
08:16:58 INFO: generate g g > e+ mu- ve vm~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @0 $$ z;add process g g > e- mu+ ve~ vm QCD=2 QED=4 NP=0 [sqrvirt=QCD] @1 $$ z;add process g g > e+ ta- ve vt~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @2 $$ z;add process g g > e- ta+ ve~ vt QCD=2 QED=4 NP=0 [sqrvirt=QCD] @3 $$ z;add process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @4 $$ z;add process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP=0 [sqrvirt=QCD] @5 $$ z;
08:17:19 INFO: Done 21.73
08:17:19 INFO: generating the square matrix element for reweighting (second model and/or processes)
08:17:19 INFO: generate g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @0 $$ z;add process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @1 $$ z;add process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @2 $$ z;add process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @3 $$ z;add process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @4 $$ z;add process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==2 [sqrvirt=QCD] @5 $$ z;
08:18:56 INFO: Done 96.32
08:19:39 Do you want to edit a card (press enter to bypass editing)?
08:19:39 /------------------------------------------------------------\
08:19:39 | 1. param : param_card.dat |
08:19:39 | 9. plot : plot_card.dat |
08:19:39 \------------------------------------------------------------/
08:19:39 you can also
08:19:39 - enter the path to a valid card or banner.
08:19:39 - use the 'set' command to modify a parameter directly.
08:19:39 The set option works only for param_card and run_card.
08:19:39 Type 'help set' for more information on this command.
08:19:39 - call an external program (ASperGE/MadWidth/...).
08:19:39 Type 'help' for the list of available command
08:19:39 [0, done, 1, param, 9, plot, enter path]
08:19:39 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:39 INFO: ctp
08:19:39 INFO: Event nb 0 current time: 08h19
08:19:40 INFO: Event nb 10 0.61s
08:19:40 INFO: Event nb 20 0.89s
08:19:40 INFO: Event nb 30 1.2s
08:19:41 INFO: Event nb 40 1.4s
08:19:41 INFO: Event nb 50 1.7s
08:19:41 INFO: All event done (nb_event: 60) 2s
08:19:41 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:41 INFO: new cross-section is : 0.000581752 pb (indicative error: 6.55964e-05 pb)
08:19:41 launch --rwgt_name=ctG --rwgt_info=ctG
08:19:41 Do you want to edit a card (press enter to bypass editing)?
08:19:41 /------------------------------------------------------------\
08:19:41 | 1. param : param_card.dat |
08:19:41 | 9. plot : plot_card.dat |
08:19:41 \------------------------------------------------------------/
08:19:41 you can also
08:19:41 - enter the path to a valid card or banner.
08:19:41 - use the 'set' command to modify a parameter directly.
08:19:41 The set option works only for param_card and run_card.
08:19:41 Type 'help set' for more information on this command.
08:19:41 - call an external program (ASperGE/MadWidth/...).
08:19:41 Type 'help' for the list of available command
08:19:41 [0, done, 1, param, 9, plot, enter path]
08:19:41 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:41 INFO: ctG
08:19:41 INFO: Event nb 0 current time: 08h19
08:19:42 INFO: Event nb 10 0.29s
08:19:42 INFO: Event nb 20 0.58s
08:19:42 INFO: Event nb 30 0.85s
08:19:43 INFO: Event nb 40 1.1s
08:19:43 INFO: Event nb 50 1.4s
08:19:43 INFO: All event done (nb_event: 60) 1.7s
08:19:43 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:43 INFO: new cross-section is : -0.0028739 pb (indicative error: 0.000246268 pb)
08:19:43 launch --rwgt_name=cpW --rwgt_info=cpW
08:19:43 Do you want to edit a card (press enter to bypass editing)?
08:19:43 /------------------------------------------------------------\
08:19:43 | 1. param : param_card.dat |
08:19:43 | 9. plot : plot_card.dat |
08:19:43 \------------------------------------------------------------/
08:19:43 you can also
08:19:43 - enter the path to a valid card or banner.
08:19:43 - use the 'set' command to modify a parameter directly.
08:19:43 The set option works only for param_card and run_card.
08:19:43 Type 'help set' for more information on this command.
08:19:43 - call an external program (ASperGE/MadWidth/...).
08:19:43 Type 'help' for the list of available command
08:19:43 [0, done, 1, param, 9, plot, enter path]
08:19:43 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:43 INFO: cpW
08:19:43 INFO: Event nb 0 current time: 08h19
08:19:44 INFO: Event nb 10 0.28s
08:19:44 INFO: Event nb 20 0.56s
08:19:44 INFO: Event nb 30 0.84s
08:19:44 INFO: Event nb 40 1.1s
08:19:45 INFO: Event nb 50 1.4s
08:19:45 INFO: All event done (nb_event: 60) 1.7s
08:19:45 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:45 INFO: new cross-section is : 0.00255995 pb (indicative error: 0.000222189 pb)
08:19:45 launch --rwgt_name=cpDc --rwgt_info=cpDc
08:19:45 Do you want to edit a card (press enter to bypass editing)?
08:19:45 /------------------------------------------------------------\
08:19:45 | 1. param : param_card.dat |
08:19:45 | 9. plot : plot_card.dat |
08:19:45 \------------------------------------------------------------/
08:19:45 you can also
08:19:45 - enter the path to a valid card or banner.
08:19:45 - use the 'set' command to modify a parameter directly.
08:19:45 The set option works only for param_card and run_card.
08:19:45 Type 'help set' for more information on this command.
08:19:45 - call an external program (ASperGE/MadWidth/...).
08:19:45 Type 'help' for the list of available command
08:19:45 [0, done, 1, param, 9, plot, enter path]
08:19:45 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:45 INFO: cpDc
08:19:45 INFO: Event nb 0 current time: 08h19
08:19:45 INFO: Event nb 10 0.28s
08:19:46 INFO: Event nb 20 0.56s
08:19:46 INFO: Event nb 30 0.83s
08:19:46 INFO: Event nb 40 1.1s
08:19:47 INFO: Event nb 50 1.4s
08:19:47 INFO: All event done (nb_event: 60) 1.7s
08:19:47 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:47 INFO: new cross-section is : 0.000286262 pb (indicative error: 3.23214e-05 pb)
08:19:47 launch --rwgt_name=cdp --rwgt_info=cdp
08:19:47 Do you want to edit a card (press enter to bypass editing)?
08:19:47 /------------------------------------------------------------\
08:19:47 | 1. param : param_card.dat |
08:19:47 | 9. plot : plot_card.dat |
08:19:47 \------------------------------------------------------------/
08:19:47 you can also
08:19:47 - enter the path to a valid card or banner.
08:19:47 - use the 'set' command to modify a parameter directly.
08:19:47 The set option works only for param_card and run_card.
08:19:47 Type 'help set' for more information on this command.
08:19:47 - call an external program (ASperGE/MadWidth/...).
08:19:47 Type 'help' for the list of available command
08:19:47 [0, done, 1, param, 9, plot, enter path]
08:19:47 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:47 INFO: cdp
08:19:47 INFO: Event nb 0 current time: 08h19
08:19:48 INFO: Event nb 10 0.28s
08:19:48 INFO: Event nb 20 0.56s
08:19:48 INFO: Event nb 30 0.84s
08:19:49 INFO: Event nb 40 1.1s
08:19:49 INFO: Event nb 50 1.4s
08:19:49 INFO: All event done (nb_event: 60) 1.7s
08:19:49 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:49 INFO: new cross-section is : -0.00114275 pb (indicative error: 9.39133e-05 pb)
08:19:49 launch --rwgt_name=c3pl1 --rwgt_info=c3pl1
08:19:49 Do you want to edit a card (press enter to bypass editing)?
08:19:49 /------------------------------------------------------------\
08:19:49 | 1. param : param_card.dat |
08:19:49 | 9. plot : plot_card.dat |
08:19:49 \------------------------------------------------------------/
08:19:49 you can also
08:19:49 - enter the path to a valid card or banner.
08:19:49 - use the 'set' command to modify a parameter directly.
08:19:49 The set option works only for param_card and run_card.
08:19:49 Type 'help set' for more information on this command.
08:19:49 - call an external program (ASperGE/MadWidth/...).
08:19:49 Type 'help' for the list of available command
08:19:49 [0, done, 1, param, 9, plot, enter path]
08:19:49 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:49 INFO: c3pl1
08:19:49 INFO: Event nb 0 current time: 08h19
08:19:50 INFO: Event nb 10 0.28s
08:19:50 INFO: Event nb 20 0.56s
08:19:50 INFO: Event nb 30 0.84s
08:19:50 INFO: Event nb 40 1.1s
08:19:51 INFO: Event nb 50 1.4s
08:19:51 INFO: All event done (nb_event: 60) 1.7s
08:19:51 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:51 INFO: new cross-section is : -0.0294462 pb (indicative error: 0.000193708 pb)
08:19:51 launch --rwgt_name=c3pl2 --rwgt_info=c3pl2
08:19:51 Do you want to edit a card (press enter to bypass editing)?
08:19:51 /------------------------------------------------------------\
08:19:51 | 1. param : param_card.dat |
08:19:51 | 9. plot : plot_card.dat |
08:19:51 \------------------------------------------------------------/
08:19:51 you can also
08:19:51 - enter the path to a valid card or banner.
08:19:51 - use the 'set' command to modify a parameter directly.
08:19:51 The set option works only for param_card and run_card.
08:19:51 Type 'help set' for more information on this command.
08:19:51 - call an external program (ASperGE/MadWidth/...).
08:19:51 Type 'help' for the list of available command
08:19:51 [0, done, 1, param, 9, plot, enter path]
08:19:51 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:51 INFO: c3pl2
08:19:51 INFO: Event nb 0 current time: 08h19
08:19:52 INFO: Event nb 10 0.28s
08:19:52 INFO: Event nb 20 0.56s
08:19:52 INFO: Event nb 30 0.84s
08:19:52 INFO: Event nb 40 1.1s
08:19:53 INFO: Event nb 50 1.4s
08:19:53 INFO: All event done (nb_event: 60) 1.7s
08:19:53 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:53 INFO: new cross-section is : -0.0301945 pb (indicative error: 0.000198414 pb)
08:19:53 launch --rwgt_name=cll1221 --rwgt_info=cll1221
08:19:53 Do you want to edit a card (press enter to bypass editing)?
08:19:53 /------------------------------------------------------------\
08:19:53 | 1. param : param_card.dat |
08:19:53 | 9. plot : plot_card.dat |
08:19:53 \------------------------------------------------------------/
08:19:53 you can also
08:19:53 - enter the path to a valid card or banner.
08:19:53 - use the 'set' command to modify a parameter directly.
08:19:53 The set option works only for param_card and run_card.
08:19:53 Type 'help set' for more information on this command.
08:19:53 - call an external program (ASperGE/MadWidth/...).
08:19:53 Type 'help' for the list of available command
08:19:53 [0, done, 1, param, 9, plot, enter path]
08:19:53 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:53 INFO: cll1221
08:19:53 INFO: Event nb 0 current time: 08h19
08:19:53 INFO: Event nb 10 0.29s
08:19:54 INFO: Event nb 20 0.56s
08:19:54 INFO: Event nb 30 0.84s
08:19:54 INFO: Event nb 40 1.1s
08:19:55 INFO: Event nb 50 1.4s
08:19:55 INFO: All event done (nb_event: 60) 1.7s
08:19:55 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:55 INFO: new cross-section is : 0.0447074 pb (indicative error: 0.000691317 pb)
08:19:55 launch --rwgt_name=cp --rwgt_info=cp
08:19:55 Do you want to edit a card (press enter to bypass editing)?
08:19:55 /------------------------------------------------------------\
08:19:55 | 1. param : param_card.dat |
08:19:55 | 9. plot : plot_card.dat |
08:19:55 \------------------------------------------------------------/
08:19:55 you can also
08:19:55 - enter the path to a valid card or banner.
08:19:55 - use the 'set' command to modify a parameter directly.
08:19:55 The set option works only for param_card and run_card.
08:19:55 Type 'help set' for more information on this command.
08:19:55 - call an external program (ASperGE/MadWidth/...).
08:19:55 Type 'help' for the list of available command
08:19:55 [0, done, 1, param, 9, plot, enter path]
08:19:55 >INFO: starts to compute weight for events with the following modification to the param_card:
08:19:55 INFO: cp
08:19:55 INFO: Event nb 0 current time: 08h19
08:19:55 INFO: Event nb 10 0.28s
08:19:56 INFO: Event nb 20 0.56s
08:19:56 INFO: Event nb 30 0.83s
08:19:56 INFO: Event nb 40 1.1s
08:19:56 INFO: Event nb 50 1.4s
08:19:57 INFO: All event done (nb_event: 60) 1.7s
08:19:57 INFO: Event /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/Events/GridRun_32289/unweighted_events.lhe have now the additional weight
08:19:57 INFO: new cross-section is : 4.60099e-07 pb (indicative error: 9.97455e-09 pb)
08:19:57 change process g g > e+ mu- ve vm~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @0 $$ z
08:19:57 change process g g > e- mu+ ve~ vm QCD=2 QED=4 NP^2==4 [noborn=QCD] @1 $$ z --add
08:19:57 change process g g > e+ ta- ve vt~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @2 $$ z --add
08:19:57 change process g g > e- ta+ ve~ vt QCD=2 QED=4 NP^2==4 [noborn=QCD] @3 $$ z --add
08:19:57 change process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP^2==4 [noborn=QCD] @4 $$ z --add
08:19:57 change process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP^2==4 [noborn=QCD] @5 $$ z --add
08:19:57 launch --rwgt_name=ctp_squared --rwgt_info=ctp_squared
08:19:57 INFO: detected model: /project/atlas/Users/zwolffs/madgraph/MG5_aMC_v3_5_1/models/SMEFTatNLOcustom-NLO_renorscaled. Loading...
08:20:08 INFO: generating the square matrix element for reweighting
08:20:08 INFO: generate g g > e+ mu- ve vm~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @0 $$ z;add process g g > e- mu+ ve~ vm QCD=2 QED=4 NP=0 [sqrvirt=QCD] @1 $$ z;add process g g > e+ ta- ve vt~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @2 $$ z;add process g g > e- ta+ ve~ vt QCD=2 QED=4 NP=0 [sqrvirt=QCD] @3 $$ z;add process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP=0 [sqrvirt=QCD] @4 $$ z;add process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP=0 [sqrvirt=QCD] @5 $$ z;
08:20:20 Command "reweight GridRun_32289 -f" interrupted with error:
08:20:20 InvalidCmd : Could not remove directory /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/rw_me.
08:20:20 Please report this bug on https://bugs.launchpad.net/mg5amcnlo
08:20:20 More information is found in '/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/GridRun_32289_tag_1_debug.log'.
08:20:20 Please attach this file to your report.
08:20:20 Do you want to continue? [y, n][60s to answer]
08:20:20 >INFO: storing files of previous run
08:20:20 gzipping output file: unweighted_events.lhe
08:20:20 INFO: Done
08:20:20 INFO:
08:20:20 INFO: Original cross-section: 0.032057 +- 0.0028508 pb (cross-section from sum of weights: 0.1843617)
08:20:20 INFO: Computed cross-section:
08:20:20 INFO: c3pl1 : -0.02944615117 +- 0.000193707873032 pb
08:20:20 INFO: c3pl2 : -0.0301944716947 +- 0.000198414312412 pb
08:20:20 INFO: cdp : -0.0011427476328 +- 9.39132642185e-05 pb
08:20:20 INFO: cll1221 : 0.0447074071999 +- 0.0006913171431 pb
08:20:20 INFO: cp : 4.60099392099e-07 +- 9.97454738333e-09 pb
08:20:20 INFO: cpDc : 0.000286262032439 +- 3.23214212361e-05 pb
08:20:20 INFO: cpW : 0.00255995459664 +- 0.00022218935543 pb
08:20:20 INFO: ctG : -0.00287389714602 +- 0.000246267912282 pb
08:20:20 INFO: ctp : 0.000581751765446 +- 6.55963839343e-05 pb
08:20:20 quit
08:20:20 INFO:
08:20:21 launch in debug mode
08:20:21 Py:MadGraphUtils INFO Finished reweighting
08:20:21 Py:MadGraphUtils INFO Events were already in place
08:20:21 Py:MadGraphUtils INFO For your information, ls of /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1:
08:20:21 Py:MadGraphUtils INFO ['100005.1_4091599_4091599.condorlog', '100005.1_4091599_4091599.log', '100005.1_4091599_4091599.stderr', 'Cards_bkup', 'ME5_debug', 'MGC_LHAPDF', 'events.lhe.gz', 'generate_1.sh', 'lhapdf-config', 'log.generate', 'madevent', 'mem.full.generate', 'mem.summary.generate.json_snapshot', 'py.py', 'run.sh', 'runargs.generate.py', 'runwrapper.generate.sh']
08:20:21 Py:MadGraphUtils INFO Moving generated events to be in correct format for arrange_output().
08:20:21 Py:MadGraphUtils INFO Unzipping generated events.
08:20:21 Py:MadGraphUtils INFO Moving file over to madevent/Events/GridRun_32289/unweighted_events.lhe
08:20:21 Py:MadGraphUtils INFO Re-zipping into dataset name madevent/Events/GridRun_32289/unweighted_events.lhe.gz
08:20:21 Py:MadGraphUtils INFO Finished at Thu Oct 3 08:20:21 2024
08:20:21 Py:MadGraphUtils INFO Restoring original LHAPDF env variables:
08:20:21 Py:MadGraphUtils INFO LHAPATH=/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current:/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/LCG_88b/MCGenerators/lhapdf/6.2.3/x86_64-centos7-gcc62-opt/share/LHAPDF
08:20:21 Py:MadGraphUtils INFO LHAPDF_DATA_PATH=/cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current:/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/LCG_88b/MCGenerators/lhapdf/6.2.3/x86_64-centos7-gcc62-opt/share/LHAPDF
08:20:21 Py:MadGraphUtils INFO Unzipping generated events.
08:20:21 Py:MadGraphUtils INFO Putting a copy in place for the transform.
08:20:21 Py:MadGraphUtils INFO Removed 18 empty lines from LHEF
08:20:21 Py:MadGraphUtils INFO The following 155 weights have been written to the LHE file: MUR0.5_MUF0.5_PDF260000,MUR0.5_MUF1.0_PDF260000,MUR0.5_MUF2.0_PDF260000,MUR1.0_MUF0.5_PDF260000,MUR1.0_MUF2.0_PDF260000,MUR2.0_MUF0.5_PDF260000,MUR2.0_MUF1.0_PDF260000,MUR2.0_MUF2.0_PDF260000,MUR1.0_MUF1.0_PDF260000,MUR1.0_MUF1.0_PDF260001,MUR1.0_MUF1.0_PDF260002,MUR1.0_MUF1.0_PDF260003,MUR1.0_MUF1.0_PDF260004,MUR1.0_MUF1.0_PDF260005,MUR1.0_MUF1.0_PDF260006,MUR1.0_MUF1.0_PDF260007,MUR1.0_MUF1.0_PDF260008,MUR1.0_MUF1.0_PDF260009,MUR1.0_MUF1.0_PDF260010,MUR1.0_MUF1.0_PDF260011,MUR1.0_MUF1.0_PDF260012,MUR1.0_MUF1.0_PDF260013,MUR1.0_MUF1.0_PDF260014,MUR1.0_MUF1.0_PDF260015,MUR1.0_MUF1.0_PDF260016,MUR1.0_MUF1.0_PDF260017,MUR1.0_MUF1.0_PDF260018,MUR1.0_MUF1.0_PDF260019,MUR1.0_MUF1.0_PDF260020,MUR1.0_MUF1.0_PDF260021,MUR1.0_MUF1.0_PDF260022,MUR1.0_MUF1.0_PDF260023,MUR1.0_MUF1.0_PDF260024,MUR1.0_MUF1.0_PDF260025,MUR1.0_MUF1.0_PDF260026,MUR1.0_MUF1.0_PDF260027,MUR1.0_MUF1.0_PDF260028,MUR1.0_MUF1.0_PDF260029,MUR1.0_MUF1.0_PDF260030,MUR1.0_MUF1.0_PDF260031,MUR1.0_MUF1.0_PDF260032,MUR1.0_MUF1.0_PDF260033,MUR1.0_MUF1.0_PDF260034,MUR1.0_MUF1.0_PDF260035,MUR1.0_MUF1.0_PDF260036,MUR1.0_MUF1.0_PDF260037,MUR1.0_MUF1.0_PDF260038,MUR1.0_MUF1.0_PDF260039,MUR1.0_MUF1.0_PDF260040,MUR1.0_MUF1.0_PDF260041,MUR1.0_MUF1.0_PDF260042,MUR1.0_MUF1.0_PDF260043,MUR1.0_MUF1.0_PDF260044,MUR1.0_MUF1.0_PDF260045,MUR1.0_MUF1.0_PDF260046,MUR1.0_MUF1.0_PDF260047,MUR1.0_MUF1.0_PDF260048,MUR1.0_MUF1.0_PDF260049,MUR1.0_MUF1.0_PDF260050,MUR1.0_MUF1.0_PDF260051,MUR1.0_MUF1.0_PDF260052,MUR1.0_MUF1.0_PDF260053,MUR1.0_MUF1.0_PDF260054,MUR1.0_MUF1.0_PDF260055,MUR1.0_MUF1.0_PDF260056,MUR1.0_MUF1.0_PDF260057,MUR1.0_MUF1.0_PDF260058,MUR1.0_MUF1.0_PDF260059,MUR1.0_MUF1.0_PDF260060,MUR1.0_MUF1.0_PDF260061,MUR1.0_MUF1.0_PDF260062,MUR1.0_MUF1.0_PDF260063,MUR1.0_MUF1.0_PDF260064,MUR1.0_MUF1.0_PDF260065,MUR1.0_MUF1.0_PDF260066,MUR1.0_MUF1.0_PDF260067,MUR1.0_MUF1.0_PDF260068,MUR1.0_MUF1.0_PDF260069,MUR1.0_MUF1.0_PDF260070,MUR1.0_MUF1.0_PDF260071,MUR1.0_MUF1.0_PDF260072,MUR1.0_MUF1.0_PDF260073,MUR1.0_MUF1.0_PDF260074,MUR1.0_MUF1.0_PDF260075,MUR1.0_MUF1.0_PDF260076,MUR1.0_MUF1.0_PDF260077,MUR1.0_MUF1.0_PDF260078,MUR1.0_MUF1.0_PDF260079,MUR1.0_MUF1.0_PDF260080,MUR1.0_MUF1.0_PDF260081,MUR1.0_MUF1.0_PDF260082,MUR1.0_MUF1.0_PDF260083,MUR1.0_MUF1.0_PDF260084,MUR1.0_MUF1.0_PDF260085,MUR1.0_MUF1.0_PDF260086,MUR1.0_MUF1.0_PDF260087,MUR1.0_MUF1.0_PDF260088,MUR1.0_MUF1.0_PDF260089,MUR1.0_MUF1.0_PDF260090,MUR1.0_MUF1.0_PDF260091,MUR1.0_MUF1.0_PDF260092,MUR1.0_MUF1.0_PDF260093,MUR1.0_MUF1.0_PDF260094,MUR1.0_MUF1.0_PDF260095,MUR1.0_MUF1.0_PDF260096,MUR1.0_MUF1.0_PDF260097,MUR1.0_MUF1.0_PDF260098,MUR1.0_MUF1.0_PDF260099,MUR1.0_MUF1.0_PDF260100,MUR1.0_MUF1.0_PDF90400,MUR1.0_MUF1.0_PDF90401,MUR1.0_MUF1.0_PDF90402,MUR1.0_MUF1.0_PDF90403,MUR1.0_MUF1.0_PDF90404,MUR1.0_MUF1.0_PDF90405,MUR1.0_MUF1.0_PDF90406,MUR1.0_MUF1.0_PDF90407,MUR1.0_MUF1.0_PDF90408,MUR1.0_MUF1.0_PDF90409,MUR1.0_MUF1.0_PDF90410,MUR1.0_MUF1.0_PDF90411,MUR1.0_MUF1.0_PDF90412,MUR1.0_MUF1.0_PDF90413,MUR1.0_MUF1.0_PDF90414,MUR1.0_MUF1.0_PDF90415,MUR1.0_MUF1.0_PDF90416,MUR1.0_MUF1.0_PDF90417,MUR1.0_MUF1.0_PDF90418,MUR1.0_MUF1.0_PDF90419,MUR1.0_MUF1.0_PDF90420,MUR1.0_MUF1.0_PDF90421,MUR1.0_MUF1.0_PDF90422,MUR1.0_MUF1.0_PDF90423,MUR1.0_MUF1.0_PDF90424,MUR1.0_MUF1.0_PDF90425,MUR1.0_MUF1.0_PDF90426,MUR1.0_MUF1.0_PDF90427,MUR1.0_MUF1.0_PDF90428,MUR1.0_MUF1.0_PDF90429,MUR1.0_MUF1.0_PDF90430,MUR1.0_MUF1.0_PDF90431,MUR1.0_MUF1.0_PDF90432,MUR1.0_MUF1.0_PDF266000,MUR1.0_MUF1.0_PDF265000,MUR1.0_MUF1.0_PDF13100,MUR1.0_MUF1.0_PDF25200,ctp,ctG,cpW,cpDc,cdp,c3pl1,c3pl2,cll1221,cp
08:20:21 Py:MadGraphUtils INFO Checking whether the following expected weights are in LHE file: ctp,ctG,cpW,cpDc,cdp,c3pl1,c3pl2,cll1221,cp,ctp_squared,ctG_squared,cpW_squared,cpDc_squared,cdp_squared,c3pl1_squared,c3pl2_squared,cll1221_squared,cp_squared,ctp_ctG,ctp_cpW,ctp_cpDc,ctp_cdp,ctp_c3pl1,ctp_c3pl2,ctp_cll1221,ctp_cp,ctG_cpW,ctG_cpDc,ctG_cdp,ctG_c3pl1,ctG_c3pl2,ctG_cll1221,ctG_cp,cpW_cpDc,cpW_cdp,cpW_c3pl1,cpW_c3pl2,cpW_cll1221,cpW_cp,cpDc_cdp,cpDc_c3pl1,cpDc_c3pl2,cpDc_cll1221,cpDc_cp,cdp_c3pl1,cdp_c3pl2,cdp_cll1221,cdp_cp,c3pl1_c3pl2,c3pl1_cll1221,c3pl1_cp,c3pl2_cll1221,c3pl2_cp,cll1221_cp,cpG,cpG_squared,ctp_cpG,ctG_cpG,cpW_cpG,cpDc_cpG,cdp_cpG,c3pl1_cpG,c3pl2_cpG,cll1221_cpG,cp_cpG,MUR1.0_MUF1.0_PDF260000,MUR1.0_MUF1.0_PDF260001,MUR1.0_MUF1.0_PDF90401,MUR1.0_MUF1.0_PDF266000,MUR1.0_MUF1.0_PDF265000,MUR1.0_MUF1.0_PDF13100,MUR1.0_MUF1.0_PDF25200,MUR0.5_MUF0.5_PDF260000,MUR0.5_MUF1.0_PDF260000,MUR0.5_MUF2.0_PDF260000,MUR1.0_MUF0.5_PDF260000,MUR1.0_MUF1.0_PDF260000,MUR1.0_MUF2.0_PDF260000,MUR2.0_MUF0.5_PDF260000,MUR2.0_MUF1.0_PDF260000,MUR2.0_MUF2.0_PDF260000
08:20:21 Shortened traceback (most recent user call last):
08:20:21 File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/AthGeneration/21.6.77/InstallArea/x86_64-centos7-gcc62-opt/jobOptions/EvgenJobTransforms/skel.GENtoEVGEN.py", line 271, in <module>
08:20:21 include(jofile)
08:20:21 File "/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/100005/mc.100005_MadgraphPythia8EvtGen_gg2l2v_DF.py", line 489, in <module>
08:20:21 arrange_output(process_dir=process_dir,runArgs=runArgs,lhe_version=3,saveProcDir=saveRundetails)
08:20:21 File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/AthGeneration/21.6.77/InstallArea/x86_64-centos7-gcc62-opt/python/MadGraphControl/MadGraphUtils.py", line 1506, in arrange_output
08:20:21 raise RuntimeError("Did not find expected weight "+w+" in lhe file. Did the reweight or systematics module crash?")
08:20:21 RuntimeError: Did not find expected weight ctp_squared in lhe file. Did the reweight or systematics module crash?
08:20:21 Py:Athena INFO leaving with code 8: "an unknown exception occurred"
```

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

So the interesting line in the log is this one:

08:20:20 InvalidCmd : Could not remove directory /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/rw_me.
08:20:20 Please report this bug on https://bugs.launchpad.net/mg5amcnlo
08:20:20 More information is found in '/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/GridRun_32289_tag_1_debug.log'.

Do you have the content of that file? (might be usefull if I do not reproduce the issue which is the case so far)

I'm trying to reproduce the issue with this example:

generate p p > e+ e-
output
launch
reweight =ON
~/rwt_card.dat
use_syst=False

with ~/rwt_card.dat:

change mode LO # Define type of Reweighting. For LO sample this command
                   # has no effect since only "LO" mode is allowed.

change process p p > e+ e- [noborn=QCD]
launch
change process p p > e+ e- [virt=QCD]
launch

Those lines do work with 2.9.21, will test now with 3.5.6 and 3.6.0 (will test with python3.10)
(maybe you can give it test too, such that we can compare)

Cheers,

Olivier

Revision history for this message
Zef Wolffs (zeff020) said : 		#4

I see, I use v2.9.3. Unfortunately due to being behind this ATLAS infrastructure it is not easy for me to switch to a newer version. For the "change mode LO", is that the right mode to use for both? Since the former is at NLO in the SMEFT model that I am using. So I would think that I would have to "change mode NLO" for the "[noborn=QCD]" (NLO) and "change mode LO" for the "[virt=QCD]" (LO) process.

The contents of that file that you requested is the following, I find it hard to tell why it couldn't remove that directory though from these logs:

```
#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.3 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
reweight GridRun_32289 -f
Traceback (most recent call last):
  File "/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/bin/internal/extended_cmd.py", line 1541, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/bin/internal/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/bin/internal/common_run_interface.py", line 2175, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/extended_cmd.py", line 1689, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/reweight_interface.py", line 503, in do_launch
    self.create_standalone_directory()
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/reweight_interface.py", line 1693, in create_standalone_directory
    has_nlo = self.create_standalone_tree_directory(data, second)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/various/misc.py", line 109, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/reweight_interface.py", line 1427, in create_standalone_tree_directory
    mgcmd.exec_cmd(commandline, precmd=True)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/master_interface.py", line 299, in do_output
    self.cmd.do_output(self, line, *args, **opts)
  File "/cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/madgraph5amc/2.9.3.atlas2-f2a0f/x86_64-centos7-gcc62-opt/madgraph/interface/loop_interface.py", line 461, in do_output
    %str(self._export_dir))
InvalidCmd: Could not remove directory /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/madevent/rw_me.
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : None
             cluster_memory : None
               cluster_size : 100
              cluster_queue : None
                    nb_core : 1 (user set)
               cluster_time : None
                   run_mode : 0 (user set)

                      Configuration Options
                      ---------------------
                text_editor : None
          f2py_compiler_py2 : f2py2 (user set)
          f2py_compiler_py3 : None
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : /cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/syscalc/1.1.4p3-61111/x86_64-centos7-gcc62-opt (user set)
           madanalysis_path : None (user set)
                     lhapdf : /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/lhapdf-config (user set)
         cluster_local_path : None
              f2py_compiler : f2py2 (user set)
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
                 lhapdf_py3 : /cvmfs/atlas.cern.ch/repo/sw/software/21.6/sw/lcg/releases/MCGenerators/lhapdf/6.2.3-0ff65/x86_64-centos7-gcc62-opt/bin/lhapdf-config (user set)
                 lhapdf_py2 : /data/atlas/users/zwolffs/HWWOffshel-EFT/processes/ggF-0j-forSubmission/run-100005/results/100005.1/lhapdf-config (user set)
           fortran_compiler : None
                auto_update : 0 (user set)
        exrootanalysis_path : None (user set)
                    timeout : 60
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.3 2021-03-25 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set group_subprocesses Auto
set ignore_six_quark_processes False
set max_t_for_channel 99
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set default_unset_couplings 99
set zerowidth_tchannel True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme True
set max_npoint_for_channel 4
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 vl = ve vm vt
define vl~ = ve~ vm~ vt~
define p = g u c d s b u~ c~ d~ s~ b~
define j = g u c d s b u~ c~ d~ s~ b~
import model /project/atlas/Users/zwolffs/madgraph/MG5_aMC_v3_5_1/mode\
ls/SMEFTatNLOcustom-NLO_renorscaled
generate g g > e+ mu- ve vm~ QCD=2 QED=4 NP=0 [noborn=QCD] @0 $$ z
add process g g > e- mu+ ve~ vm QCD=2 QED=4 NP=0 [noborn=QCD] @1 $$ z
add process g g > e+ ta- ve vt~ QCD=2 QED=4 NP=0 [noborn=QCD] @2 $$ z
add process g g > e- ta+ ve~ vt QCD=2 QED=4 NP=0 [noborn=QCD] @3 $$ z
add process g g > mu+ ta- vm vt~ QCD=2 QED=4 NP=0 [noborn=QCD] @4 $$ \
z
add process g g > mu- ta+ vm~ vt QCD=2 QED=4 NP=0 [noborn=QCD] @5 $$ \
z
output -f
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR DIM6
###################################
Block dim6
   1 1000 # Lambda
   2 1e-05 # cpDC
   3 1e-05 # cpWB
   4 1e-05 # cdp
   5 1e-05 # cp
   6 1e-05 # cWWW
   9 1e-05 # cpW
   10 1e-05 # cpBB

###################################
## INFORMATION FOR DIM62F
###################################
Block dim62f
   1 1e-05 # cpl1
   2 1e-05 # cpl2
   3 1e-05 # cpl3
   4 1e-05 # c3pl1
   5 1e-05 # c3pl2
   6 1e-05 # c3pl3
   7 1e-05 # cpe
   8 1e-05 # cpmu
   9 1e-05 # cpta
   10 1e-05 # cpqMi
   11 1e-05 # cpq3i
   12 1e-05 # cpQ3
   13 1e-05 # cpQM
   14 1e-05 # cpu
   15 1e-05 # cpt
   16 1e-05 # cpd
   19 1e-05 # ctp
   22 1e-05 # ctZ
   23 1e-05 # ctW
   24 1e-05 # ctG

###################################
## INFORMATION FOR DIM64F
###################################
Block dim64f
   1 1e-05 # cQq83
   2 1e-05 # cQq81
   3 1e-05 # cQu8
   4 1e-05 # ctq8
   6 1e-05 # cQd8
   7 1e-05 # ctu8
   8 1e-05 # ctd8
   10 1e-05 # cQq13
   11 1e-05 # cQq11
   12 1e-05 # cQu1
   13 1e-05 # ctq1
   14 1e-05 # cQd1
   16 1e-05 # ctu1
   17 1e-05 # ctd1
   19 1e-05 # cQQ8
   20 1e-05 # cQQ1
   21 1e-05 # cQt1
   23 1e-05 # ctt1
   25 1e-05 # cQt8

###################################
## INFORMATION FOR DIM64F2L
###################################
Block dim64f2l
   1 1e-05 # cQlM1
   2 1e-05 # cQlM2
   3 1e-05 # cQl31
   4 1e-05 # cQl32
   5 1e-05 # cQe1
   6 1e-05 # cQe2
   7 1e-05 # ctl1
   8 1e-05 # ctl2
   9 1e-05 # cte1
   10 1e-05 # cte2
   13 1e-05 # cQlM3
   14 1e-05 # cQl33
   15 1e-05 # cQe3
   16 1e-05 # ctl3
   17 1e-05 # cte3

###################################
## INFORMATION FOR DIM64F4L
###################################
Block dim64f4l
   1 1e-05 # cll1111
   2 1e-05 # cll2222
   3 1e-05 # cll3333
   4 1e-05 # cll1122
   5 1e-05 # cll1133
   6 1e-05 # cll2233
   7 1e-05 # cll1221
   8 1e-05 # cll1331
   9 1e-05 # cll2332

###################################
## INFORMATION FOR MASS
###################################
Block mass
    6 1.720000e+02 # MT
   23 9.118760e+01 # MZ
   24 7.982440e+01 # MW
   25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000e+00 # d : 0.0
  2 0.000000e+00 # u : 0.0
  3 0.000000e+00 # s : 0.0
  4 0.000000e+00 # c : 0.0
  5 0.000000e+00 # b : 0.0
  11 0.000000e+00 # e- : 0.0
  12 0.000000e+00 # ve : 0.0
  13 0.000000e+00 # mu- : 0.0
  14 0.000000e+00 # vm : 0.0
  15 0.000000e+00 # ta- : 0.0
  16 0.000000e+00 # vt : 0.0
  21 0.000000e+00 # g : 0.0
  22 0.000000e+00 # a : 0.0
  9000002 9.118760e+01 # ghz : MZ
  9000003 7.982440e+01 # ghwp : MW
  9000004 7.982440e+01 # ghwm : MW

###################################
## INFORMATION FOR RENOR
###################################
Block renor
    1 1.750000e+02 # mueft

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    2 1.166370e-05 # Gf
    3 1.184000e-01 # aS (Note that Parameter not used if you use a PDF set)

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.470800e+00 # WT
DECAY 23 2.416023e+00 # WZ
DECAY 24 2.002950e+00 # WW
DECAY 25 4.07e-03
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 9000002 2.416023e+00 # ghz : WZ
DECAY 9000003 2.002950e+00 # ghwp : WW
DECAY 9000004 2.002950e+00 # ghwm : WW
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 9000001 # gha
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000002 # ghz
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000003 # ghwp
        1 3 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000004 # ghwm
        1 -3 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
Block QNUMBERS 9000005 # ghg
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#*********************************************************************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#*********************************************************************
60 = nevents ! Number of unweighted events requested
32289 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron, 4=photon from muon *
#*********************************************************************
1 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
6500.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
lhapdf = pdlabel ! PDF set
260000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
True = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# Phase-Space Optimization strategy (basic options)
#*********************************************************************
1 = nhel ! using helicities importance sampling or not.
! 0: sum over helicity, 1: importance sampling
1 = sde_strategy ! default integration strategy (hep-ph/2021.xxxxx)
! 1 is old strategy (using amp square)
! 2 is new strategy (using only the denominator)

#*********************************************************************
# Phase-Space Optim (advanced)
#*********************************************************************
0 = hard_survey ! force to have better estimate of the integral at survey for difficult mode like interference
-1.0 = tmin_for_channel ! limit the non-singular reach of --some-- channel of integration related to T-channel diagram (value between -1 and 0), -1 is no impact
-1 = survey_splitting ! for loop-induced control how many core are used at survey for the computation of a single iteration.
2 = survey_nchannel_per_job ! control how many Channel are integrated inside a single job on cluster/multicore
-1 = refine_evt_by_job ! control the maximal number of events for the first iteration of the refine (larger means less jobs)
#*********************************************************************
# Compilation flag. No automatic re-compilation (need manual "make clean" in Source)
#*********************************************************************
-O = global_flag ! fortran optimization flag use for the all code.
= aloha_flag ! fortran optimization flag for aloha function. Suggestions: '-ffast-math'
= matrix_flag ! fortran optimization flag for matrix.f function. Suggestions: '-O3'

#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
10.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
-1.0 = ptlmax ! maximum pt for the charged leptons
-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
#
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
5 = etal ! max rap for the charged leptons
0.0 = etalmin ! main rap for the charged leptons
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
0.4 = drll ! min distance between leptons
-1.0 = drllmax ! max distance between leptons
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
140 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons / neutrino *
# for pair of lepton includes only same flavor, opposite charge
#*********************************************************************
0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
0.0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
20.0 = ptl1min ! minimum pt for the leading lepton in pt
10.0 = ptl2min ! minimum pt for the second lepton in pt
-1.0 = ptl1max ! maximum pt for the leading lepton in pt
-1.0 = ptl2max ! maximum pt for the second lepton in pt
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
5 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
True = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1.0,2.0','--dyn=-1','--muf=0.5,1.0,2.0','--weight_info=MUR%(mur).1f_MUF%(muf).1f_PDF%(pdf)i','--pdf=NNPDF30_nlo_as_0118,PDF4LHC15_nlo_30_pdfas,NNPDF30_nlo_as_0119@0,NNPDF30_nlo_as_0117@0,CT14nlo@0,MMHT2014nlo68clas118@0','--remove_wgts=".*MUR0.5_MUF2.0.*|.*MUR2.0_MUF0.5.*"'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'
5.0 = ptj
5.0 = ptb
10.0 = pta
2 = job_strategy
32289 = python_seed
5 = etaj
```

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#5

Just to update
- 3.5.6 does work as well
- 3.6.0 does not (but this is an unrelated issue that I'm going to fix for 3.6.1)

Tell me if you have the file, and if that simple test case does work for you or not in 3.5.1 (and if not maybe try 3.5.6)

Cheers,

Olivier

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#6

HI,

[noborn=QCD] is identical to [sqvirt=QCD]
which just say to take the loop square, while
[virt=QCD]
takes the interference between the tree amplitude and the loop

So both have only one type of matrix element and this is not enough for doing NLO re-weighting
So I would say that doing "NLO" does not make sense here.

For the version, if you are stuck with that one, then the only solution is to do multiple reweighting or something like that I guess...
Will not be able to fix that version (and that version is not even runing anymore on my machines)

Cheers,

Olivier

Revision history for this message
Zef Wolffs (zeff020) said : 		#7

Hello,

I see, for posterity I was able to try one more version: v2.9.9, for which I received the same error.

Thanks for all the help though, I'll try to see if I can find a way to have access to a newer version, and otherwise will indeed run it with separate reweigh cards.

All the best,

Zef

Can you help with this problem?

Provide an answer of your own, or ask Zef Wolffs for more information if necessary.

To post a message you must log in.