Getting a warning: "Missing information for NLO type of reweighting. Keeping the mode to LO"
I am trying to run a process with SMEFTatNLO with reweighting. I saw in some other threads that NLO reweighting is broken for MG version 3, but I am getting this same issue for both 3.5.1 and 2.9.17. I am not getting any crashing, and I am able to generate events with updated weights. However, the log file contains the line
```[1;34mWARNING: Missing information for NLO type of reweighting. Keeping the mode to LO. [0m
```
My reweight card is this:
```
change output default
change helicity False
change mode NLO
launch --rwgt_name=test_1
set DIM6 5 1.9820755275058226
launch
```
I think that LHAPDF has been properly installed, since I am able to run systematics.
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Radha Mastandrea
- Solved:
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Ok, I think I have solved this by looking at this resource : https:/
Specifically with the line "If you want to perform the NLO-accurate reweighting, you might want to have the parameter "store_rwgt_info" of the run_card on True. Note that this parameter will be automatically switched to True if you request at the time of the generation a NLO accurate re-weighting (or the loop-improved one)."
One thing that I found confusing is that the default version of the MG run card does not have this parameter in there. So I had to add a fully new line to the run card "True = store_rwgt_info", which felt a little hacky.
However,I am now running into another problem:
My reweight card is this
```
# Global setup
change output default
change helicity False
change mode NLO
launch --rwgt_
set DIM6 5 -0.7610114176229708
set DIM6 4 1.5521615637101362
launch --rwgt_
set DIM6 5 0.704670826665224
set DIM6 4 -0.9370202314941363
```
While the first reweight works, I get an error with the second reweight:
"launch --rwgt_
InvalidCmd : mg5_aMC: error: no such option: --rwgt_name
Revision history for this message
|
#2 |
This is the full bug report:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
launch
Traceback (most recent call last):
File "/global/
return self.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
self.
File "/global/
self.
File "/global/
original_
File "/global/
self.
File "/global/
stop = Cmd.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
reweight_
File "/global/
self.
File "/global/
stop = Cmd.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
out = f(self, *args, **opt)
File "/global/
cross[name] += weight[name]
KeyError: ''
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 2023-07-11 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set low_mem_
set complex_mass_scheme False
set include_
set gauge unitary
set loop_optimized_
set loop_color_flows False
set max_npoint_
set default_
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model sm
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~
import model SMEFTatNLO-NLO
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate g g > H H QED=2 QCD=2 NP=2 [QCD]
output test
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR DIM6
#######
BLOCK DIM6 #
1 1.000000e+03 # lambda
4 1.000000e-05 # set of param :1*cdp, 1*cp
#######
## INFORMATION FOR DIM62F
#######
BLOCK DIM62F #
19 1.000000e-05 # ctp
#######
## 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
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 9.118000e+01 # mueft
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
2 1.166370e-05 # gf
3 1.300000e-01 # as (note that parameter not used if you use a pdf set)
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
6 1.720000e+02 # ymt
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.354080e+00 # wt
DECAY 23 2.428230e+00 # wz
DECAY 24 2.028440e+00 # ww
DECAY 25 3.514960e-03 # wh
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.428230e+00 # ghz : wz
DECAY 9000003 2.028440e+00 # ghwp : ww
DECAY 9000004 2.028440e+00 # ghwm : ww
#######
## INFORMATION FOR QNUMBERS 9000001
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000002
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000003
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000004
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000005
#######
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/
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update to_full *
#******
#
#******
# 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 *
#******
10000 = nevents ! Number of unweighted events requested
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, *
# 2=elastic photon of proton/ion beam *
# +/-3=PDF of electron/positron beam *
# +/-4=PDF of muon/antimuon beam *
#******
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 alpha_s and its evol. *
# pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
# iww=Improved Weizsaecker-
# eva=Effective W/Z/A Approx. [2111.02442] *
# edff=EDFF in gamma-UPC [eq.(11) in 2207.03012] *
# chff=ChFF in gamma-UPC [eq.(13) in 2207.03012] *
# none=No PDF, same as lhapdf with lppx=0 *
#******
nn23lo1 = pdlabel ! PDF set
230000 = 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_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = 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.
2 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
! 2 is new strategy (using only the denominator)
# To see advanced option for Phase-Space optimization: type "update psoptim"
#******
# Customization (custom cuts/scale/
# list of files containing fortran function that overwrite default *
#******
= custom_fcts ! List of files containing user hook function
#******
# Parton level cuts definition *
#******
0.0 = dsqrt_shat ! minimal shat for full process
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
{} = 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) *
#******
{} = 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 *
#******
#******
# Minimum and maximum invariant mass for pairs *
#******
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#******
# Inclusive cuts *
#******
0.0 = ptheavy ! minimum pt for at least one heavy final state
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = 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
True = store_rwgt_info ! Store info for reweighting in LHE file
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
Revision history for this message
|
#3 |
Hi,
Sorry I did post an answer earlier today but looks like you preceed me by a couple of second, and my answer was discarded.
But yes, indeed you need to have that parameter setup to True if you want NLO accurate re-weighting.
>One thing that I found confusing is that the default version of the MG run card does not have this parameter in there. So I had to add a fully new line to the run card "True = store_rwgt_info", which felt a little hacky.
In my test today, the parameter was present by default in the run_card. So I'm confused. But in any case, we indeed allow a lot of parameter to be include in the run_card which are not present in the default card. The main reason is to try to keep the run_card as short as possible to avoid to hide the relevant information within very technical and advanced parameter.
For the error, I'm trying to reproduce the problem to see what is going on
Cheers,
Olivier
Revision history for this message
|
#4 |
If it helps, I am working with an interference process g g > h h. It seems like there have been issues with MadGraph and NLO interference reweighting, though I can't find a good solution in launchpad.
Revision history for this message
|
#5 |
Sorry for the information overload -- the bug report from #2 was generated using MG 3.5.1. if I switch to 2.9.17, I get this report
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.9.17 2023-10-19 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events run_01
Traceback (most recent call last):
File "/global/
return self.onecmåd_
File "/global/
return func(arg, **opt)
File "/global/
self.
File "/global/
original_
File "/global/
self.
File "/global/
stop = Cmd.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
self.
File "/global/
self.
File "/global/
return misc.compile(
File "/global/
raise MadGraph5Error(
madgraph.
The compilation fails with the following output message:
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
cd PDF; make
cd CERNLIB; make
cd BIAS/dummy; make
make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
make[1]: Entering directory '/pscratch/
make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
make[1]: Entering directory '/pscratch/
make[1]: Entering directory '/pscratch/
run_
355 | REWEIGHT_PDF = ' False '
| 1
Error: Symbol 'reweight_pdf' at (1) has no IMPLICIT type
run_
353 | REWEIGHT_SCALE = ' False '
| 1
Error: Symbol 'reweight_scale' at (1) has no IMPLICIT type
run_
359 | RW_FSCALE = ' 1.0 '
| 1
Error: Symbol 'rw_fscale' at (1) has no IMPLICIT type; did you mean 's_scale'?
run_
357 | RW_RSCALE = ' 1.0 '
| 1
Error: Symbol 'rw_rscale' at (1) has no IMPLICIT type; did you mean 's_scale'?
run_
351 | STORE_RWGT_INFO = ' True '
| 1
Error: Symbol 'store_rwgt_info' at (1) has no IMPLICIT type
gfortran -w -fPIC -O -ffixed-
make: *** [<builtin>: setrun.o] Error 1
make: *** Waiting for unfinished jobs....
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
ar cru ../../.
ranlib ../../.
make[1]: Leaving directory '/pscratch/
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
gfortran -w -fPIC -O -ffixed-
ar cru ../../lib/
gfortran -w -fPIC -O -ffixed-
ranlib ../../lib/
make[1]: Leaving directory '/pscratch/
gfortran -w -fPIC -O -ffixed-
ar cru ../../lib/libpdf.a Ctq6Pdf.o pdfwrap.o opendata.o pdf.o PhotonFlux.o pdg2pdf.o NNPDFDriver.o
ranlib ../../lib/libpdf.a
make[1]: Leaving directory '/pscratch/
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
automatic_
cluster_
Revision history for this message
|
#6 |
Hi,
I actually do not reproduce your error with the -rwgt_name=
This option is correctly handle for me (3.5.2) and produce results:
INFO: All event done (nb_event: 10000) 1m 22s
INFO: Event /Users/
INFO: new cross-section is : 0.175649 pb (indicative error: 0.2103 pb)
INFO: Original cross-section: 0.028482 +- 0.00011268 pb (cross-section from sum of weights: 0.0284819999999
INFO: Computed cross-section:
INFO: reweight_1 : 0.7969480101064925 +- 1.1660301774088946 pb
INFO: reweight_2 : 0.17564940234461995 +- 0.21029975624666383 pb
INFO: quit rwgt
decay_events -from_cards
I guess this means that you do not run the code in the same way as I do. So could you specify how you run the code (which sequence of executable / which command you type /...)
Concerning your 2.9.17 error, this seems an error of your setup, the error is typical of a missmatch of run_card (here it seems that you provide a NLO run_card while the code expect a LO run_card). So clearly you are not running the process that you expect in that version of the code...
Cheers,
Olivier
Revision history for this message
|
#7 |
Hi Olivier,
Thank you for this -- interesting that you are not getting the same error! Here is my entire pipeline (I just ran it so the errors are up-to-date.
Here are the commands I run (in 3.5.1)
import model SMEFTatNLO-NLO
generate g g > H H QED=2 QCD=2 NP=2 [QCD]
Then I set reweighting on, with the run cards and reweight card in the following message.
I get the same error message as before, with this log:
Perhaps I should switch to 3.5.2?
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 2023-07-11 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events run_01
Traceback (most recent call last):
File "/global/
return self.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
self.
File "/global/
original_
File "/global/
self.
File "/global/
stop = Cmd.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
reweight_
File "/global/
self.
File "/global/
stop = Cmd.onecmd_
File "/global/
return func(arg, **opt)
File "/global/
out = f(self, *args, **opt)
File "/global/
cross[name] += weight[name]
KeyError: ''
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 3.5.1 2023-07-11 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set low_mem_
set complex_mass_scheme False
set include_
set gauge unitary
set loop_optimized_
set loop_color_flows False
set max_npoint_
set default_
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model sm
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~
import model SMEFTatNLO-NLO
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate g g > H H QED=2 QCD=2 NP=2 [QCD]
output test
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR DIM6
#######
BLOCK DIM6 #
1 1.000000e+03 # lambda
4 1.000000e-05 # set of param :1*cdp, 1*cp
#######
## INFORMATION FOR DIM62F
#######
BLOCK DIM62F #
19 1.000000e-05 # ctp
#######
## 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
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 9.118000e+01 # mueft
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
2 1.166370e-05 # gf
3 1.300000e-01 # as (note that parameter not used if you use a pdf set)
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
6 1.720000e+02 # ymt
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.354080e+00 # wt
DECAY 23 2.428230e+00 # wz
DECAY 24 2.028440e+00 # ww
DECAY 25 3.514960e-03 # wh
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.428230e+00 # ghz : wz
DECAY 9000003 2.028440e+00 # ghwp : ww
DECAY 9000004 2.028440e+00 # ghwm : ww
#######
## INFORMATION FOR QNUMBERS 9000001
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000002
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000003
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000004
#######
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/
#######
## INFORMATION FOR QNUMBERS 9000005
#######
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/
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# 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 *
#******
#******
# Running parameters *
#******
#******
# 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 *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
5000 = nevents ! Number of unweighted events requested
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
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
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
90900 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# 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_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format *
#******
False = gridpack ! True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant lifetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
sum = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only) *
#******
0 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 (turn off for VBF and single top processes)
0.0 = xqcut ! minimum kt jet measure between partons
#******
# handling of the helicities: *
# 0: sum over all helicities *
# 1: importance sampling over helicities *
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
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. Define on/off-shell for "$" and decay *
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
10.0 = ptj ! minimum pt for the jets
35.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.0 = ptl ! minimum pt for the charged leptons
25.0 = misset ! minimum missing Et (sum of neutrino's momenta)
30.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-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})
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
{} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
{} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
2.5 = etaj ! max rap for the jets
2.5 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
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 = drjj ! min distance between jets
0.4 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.4 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.4 = drjl ! min distance between jet and lepton
0.4 = drab ! min distance between gamma and b
0.4 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
-1.0 = drjjmax ! max distance between jets
-1.0 = drbbmax ! max distance between b's
-1.0 = drllmax ! max distance between leptons
-1.0 = draamax ! max distance between gammas
-1.0 = drbjmax ! max distance between b and jet
-1.0 = drajmax ! max distance between gamma and jet
-1.0 = drjlmax ! max distance between jet and lepton
-1.0 = drabmax ! max distance between gamma and b
-1.0 = drblmax ! max distance between b and lepton
-1.0 = dralmax ! max distance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#******
0.0 = mmjj ! min invariant mass of a jet pair
80.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
160.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma 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_
#******
# Minimum and maximum invariant mass for all leptons *
#******
0.0 = mmnl ! min invariant mass for all leptons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all leptons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
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 = ptheavy ! minimum pt for at least one heavy final state
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth 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
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin ! inclusive Ht for all partons (including b)
-1.0 = ihtmax ! inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#******
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# WBF cuts *
#******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#******
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
False = use_syst ! Enable systematics studies
none = systematics_program ! none, systematics [python], SysCalc [deprecated, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
#******
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#******
True = store_rwgt_info ! Store info for reweighting in LHE file
False = reweight_scale
False = reweight_PDF
1.0 = rw_rscale
1.0 = rw_fscale
#******
Revision history for this message
|
#8 |
Run card:
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# 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 *
#******
#******
# Running parameters *
#******
#******
# 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 *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
5000 = nevents ! Number of unweighted events requested
0 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
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
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
nn23lo1 = pdlabel ! PDF set
90900 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# 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_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format *
#******
False = gridpack ! True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant lifetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
sum = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only) *
#******
0 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0 (turn off for VBF and single top processes)
0.0 = xqcut ! minimum kt jet measure between partons
#******
# handling of the helicities: *
# 0: sum over all helicities *
# 1: importance sampling over helicities *
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
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. Define on/off-shell for "$" and decay *
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
10.0 = ptj ! minimum pt for the jets
35.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.0 = ptl ! minimum pt for the charged leptons
25.0 = misset ! minimum missing Et (sum of neutrino's momenta)
30.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-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})
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
{} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
{} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
2.5 = etaj ! max rap for the jets
2.5 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
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 = drjj ! min distance between jets
0.4 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.4 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.4 = drjl ! min distance between jet and lepton
0.4 = drab ! min distance between gamma and b
0.4 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
-1.0 = drjjmax ! max distance between jets
-1.0 = drbbmax ! max distance between b's
-1.0 = drllmax ! max distance between leptons
-1.0 = draamax ! max distance between gammas
-1.0 = drbjmax ! max distance between b and jet
-1.0 = drajmax ! max distance between gamma and jet
-1.0 = drjlmax ! max distance between jet and lepton
-1.0 = drabmax ! max distance between gamma and b
-1.0 = drblmax ! max distance between b and lepton
-1.0 = dralmax ! max distance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#******
0.0 = mmjj ! min invariant mass of a jet pair
80.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
160.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma 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_
#******
# Minimum and maximum invariant mass for all leptons *
#******
0.0 = mmnl ! min invariant mass for all leptons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all leptons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
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 = ptheavy ! minimum pt for at least one heavy final state
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth 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
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin ! inclusive Ht for all partons (including b)
-1.0 = ihtmax ! inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#******
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# WBF cuts *
#******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#******
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
False = use_syst ! Enable systematics studies
none = systematics_program ! none, systematics [python], SysCalc [deprecated, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
#******
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#******
True = store_rwgt_info ! Store info for reweighting in LHE file
False = reweight_scale
False = reweight_PDF
1.0 = rw_rscale
1.0 = rw_fscale
#******
Revision history for this message
|
#9 |
Reweight card:
change output default
change helicity False
change mode NLO
launch --rwgt_name=r_1
set DIM6 5 -0.9803267412636116
set DIM6 4 1.8131848793805339
launch --rwgt_name=r_2
set DIM6 5 0.7150108206212078
set DIM6 4 -1.8407642184090864
Revision history for this message
|
#10 |
I am a little more convinced that something is wrong with my MG version -- to switch on reweighting, I hit "5" at this screen
/======
| 1. Choose the shower/
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = ExRoot | OFF |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\======
But this question (https:/
/================ Description =======
| 1. Type of perturbative computation | order = NLO | LO |
| 2. No MC@[N]LO matching / event generation | fixed_order = OFF | ON |
| 3. Shower the generated events | shower = PYTHIA8 | HERWIG6|
| 4. Decay onshell particles | madspin = OFF | ON|onshell |
| 5. Add weights to events for new hypp. | reweight = Not Avail | Please install module |
| 6. Run MadAnalysis5 on the events generated | madanalysis = OFF | ON |
\======
Revision history for this message
|
#11 |
Hi,
Ok, I think I understand your issue now (or at least one).
So the first question is for LO type of process (that includes loop-induced process) while the second is for NLO process.
So the point that I did not catch before is that your process is actually a loop-induced process and therefore LO accurate.
Note that the notion of "LO" is ambiguous in this context since this loop-diagram (squrared) can be seen as a part of the NNLO computation for "p p > H H" but is the LO contribution for "g g > H H". Putting this nomanclature aside they are only one order of alpha_s involve in this computation and therefore we do not have the typical complication of substracting infinities like we do have for NLO computation, so for MG5aMC we treat them within the LO framework (with the only difference is that the matrix-element is a one loop process).
This also explains that you do not have the flag "store_rwgt_info" in the run_card.dat since that flag makes only sense for NLO computation not for LO.
Now given that the computation is "LO" like, we can only do "LO" type re-weigthing (which again will be a one-loop squared matrix-element at the numerator/
Cheers,
Olivier
Revision history for this message
|
#12 |
Thank you! This clears things up.