Reweighting to points with vanishing cross-sections

Asked by Simonas Drauksas on 2020-11-09

Hi,

I have quite a few benchmarks to which I am reweighting my samples. Somewhat by mistake, some parameter points have a certain coupling set to 0, which makes the cross-section vanish. This all good and expected.

However, I have multiple parameter points with this vanishing cross-section and I think it causes some problems.
Even though in the log files I see that the reweighting seems to be done correctly and I only get zeroes where they are supposed to be, the resulting lhe file has zeroes for all the parameter points in nearly all the events.

As far as I have checked, this really seems to be caused by having a few parameter points with vanishing (identical?) cross-sections.
In principle, this is not a big deal - I will be able to avoid all but one such point with a vanishing cross-section, the problem does not occur in that case.

I am using MadGraph_aMC-2.8.0 with python2.

Sincerely,

Simonas

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Simonas Drauksas
Solved:
2020-11-16
Last query:
2020-11-16
Last reply:
2020-11-14

This question was reopened

Can you give your setup where you do have that problem such that I can reproduce your problem?

thanks,

Olivier

Simonas Drauksas (tehmeh) said : #2

Is it okay if I just copy paste a few events from the lhe file along with the header?
(I have trimmed the number of benchmarks, so that it is shorter)

<LesHouchesEvents version="3.0">
<header>
<!--
#*********************************************************************
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://madgraph.hep.uiuc.edu *
# http://madgraph.phys.ucl.ac.be *
# http://amcatnlo.cern.ch *
# *
# The MadGraph5_aMC@NLO team *
# *
#....................................................................*
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#*********************************************************************
-->
<MGVersion>
2.8.0
</MGVersion>
<MG5ProcCard>
<![CDATA[
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.8.0 2020-08-21 *
#* *
#* 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 False
set max_npoint_for_channel 0
import model VPrime_NLO
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 = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
define wp = wp+ wp-
define mu = mu+ mu-
define vmu = vm vm~
generate p p > wp, wp > mu vmu @0
add process p p > wp j, wp > mu vmu @1
add process p p > wp 2j, wp > mu vmu @2
output ./mg_processes_NLO_new/signal_pythia
]]>
</MG5ProcCard>
<MGProcCard>
<![CDATA[
#*********************************************************************
# MadGraph/MadEvent *
# http://madgraph.hep.uiuc.edu *
# *
# proc_card.dat *
#*********************************************************************
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#*********************************************************************
#*********************************************************************
# Process(es) requested : mg2 input *
#*********************************************************************
# Begin PROCESS # This is TAG. Do not modify this line
p p > wp , wp > mu vmu @0 #Process
# Be carefull the coupling are here in MG5 convention

end_coup # End the couplings input

done # this tells MG there are no more procs
# End PROCESS # This is TAG. Do not modify this line
#*********************************************************************
# Model information *
#*********************************************************************
# Begin MODEL # This is TAG. Do not modify this line
VPrime_NLO
# End MODEL # This is TAG. Do not modify this line
#*********************************************************************
# Start multiparticle definitions *
#*********************************************************************
# Begin MULTIPARTICLES # This is TAG. Do not modify this line

# End MULTIPARTICLES # This is TAG. Do not modify this line
]]>
</MGProcCard>
<MGRunCard>
<![CDATA[
#*********************************************************************
# 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 *
#*********************************************************************
  100 = nevents ! Number of unweighted events requested
 21 = 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 *
#*********************************************************************
  1 = lpp1 ! beam 1 type
  1 = lpp2 ! beam 2 type
  4000.0 = ebeam1 ! beam 1 total energy in GeV
  4000.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
  261000 = 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
#*********************************************************************
  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
#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
  1 = ickkw ! 0 no matching, 1 MLM
  1.0 = alpsfact ! scale factor for QCD emission vx
  False = chcluster ! cluster only according to channel diag
  5 = 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)
  10.0 = xqcut ! minimum kt jet measure between partons
#***********************************************************************
# 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
# 21, 1, 2, 3, 4, 5, 6 = pdgs_for_merging_cut ! PDGs for two cuts above

#*********************************************************************
#
#*********************************************************************
# 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.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)
 #*********************************************************************
 # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
 # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
 #*********************************************************************
  False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
  20.0 = ptj ! minimum pt for the jets
  10.0 = ptl ! minimum pt for the charged leptons
  0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
  -1.0 = ptjmax ! maximum pt for the jets
  -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.0 = etaj ! max rap for the jets
  2.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 = drjj ! min distance between jets
  0.4 = drjl ! min distance between jet and lepton
  -1.0 = drjjmax ! max distance between jets
  -1.0 = drjlmax ! max distance between jet and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
  0.0 = mmjj ! min invariant mass of a jet pair
  -1.0 = mmjjmax ! max invariant mass of a jet 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 *
 #*********************************************************************
  0.0 = 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 *
 #*********************************************************************
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
  0.0 = xptj ! minimum pt for at least one jet
 #*********************************************************************
 # 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
  -1.0 = ptj1max ! maximum pt for the leading jet in pt
  -1.0 = ptj2max ! maximum pt for the second jet in pt
  0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
 #*********************************************************************
 # 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)
 #*********************************************************************
 # 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
#*********************************************************************
# 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 *
#*********************************************************************
# False = use_syst ! Enable systematics studies # Commented out by MadMiner
#
# systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++] # Commented out by MadMiner
# ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset', '--alps=0.5,1,2'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule # Commented out by MadMiner
# Syscalc is deprecated but to see the associate options type'update syscalc'#*********************************************************************
# Additional hidden parameters
#*********************************************************************
  0.4 = drll # hidden_parameter
  -1.0 = ptj4max # hidden_parameter
  -1.0 = ptl1max # hidden_parameter
  0.0 = ht3min # hidden_parameter
  -1.0 = ht4max # hidden_parameter
  -1.0 = drllmax # hidden_parameter
  0.0 = ptj4min # hidden_parameter
  -1.0 = ht2max # hidden_parameter
  0.0 = ht2min # hidden_parameter
  0.0 = ptl1min # hidden_parameter
  -1.0 = ptl2max # hidden_parameter
  0.0 = ht4min # hidden_parameter
  -1.0 = ptllmax # hidden_parameter
  0.0 = mmll # hidden_parameter
  0.0 = ptllmin # hidden_parameter
  0.0 = ptj3min # hidden_parameter
  -1.0 = mmllmax # hidden_parameter
  -1.0 = ptj3max # hidden_parameter
  0.0 = xptl # hidden_parameter
  -1.0 = ht3max # hidden_parameter
  0.0 = ptl2min # hidden_parameter

#*********************************************************************
# MadMiner systematics setup *
#*********************************************************************
False = use_syst
]]>
</MGRunCard>
<slha>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.733000e+02 # mt
      23 9.118760e+01 # mz
      25 1.257000e+02 # mh
      32 3.000000e+03 # mzp
      34 2.000000e+02 # madminer
      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
      24 7.995123e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      9000002 9.118760e+01 # ghz : mz
      9000003 7.995123e+01 # ghwp : mw
      9000004 7.995123e+01 # ghwm : mw
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.279400e+02 # aewm1
      2 1.174560e-05 # gf
      3 1.184000e-01 # as
###################################
## INFORMATION FOR SSMCOUP
###################################
BLOCK SSMCOUP #
      1 1.000000e-01 # madminer
      2 0.000000e+00 # kr
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      6 1.733000e+02 # ymt
###################################
## INFORMATION FOR DECAY
###################################
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)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## 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/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 82
###################################
BLOCK QNUMBERS 82 # 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)
###################################
## INFORMATION FOR QNUMBERS 32
###################################
BLOCK QNUMBERS 32 # zp
      1 0 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 34
###################################
BLOCK QNUMBERS 34 # wp+
      1 3 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.350000e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495200e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 4.170000e-03
#
# PDG Width
DECAY 32 8.959000e+01
#
# PDG Width
DECAY 34 5.214011e-02
# BR NDA ID1 ID2 ...
   3.241072e-01 2 4 -3 # 0.0168989855111
   3.241072e-01 2 2 -1 # 0.0168989855111
   1.080357e-01 2 16 -15 # 0.00563299517038
   1.080357e-01 2 14 -13 # 0.00563299517038
   1.080357e-01 2 12 -11 # 0.00563299517038
   2.767832e-02 2 6 -5 # 0.00144315070656
#
# PDG Width
DECAY 9000002 2.495200e+00
#
# PDG Width
DECAY 9000003 2.085000e+00
#
# PDG Width
DECAY 9000004 2.085000e+00
</slha>
<initrwgt>
<weightgroup name='mg_reweighting' weight_name_strategy='includeIdInWeightName'>
<weight id='bm_wp_0_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_1_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 300.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_1_kl_1'>set param_card decay 34 0.08917674 # orig: 0.05214011
set param_card mass 34 300.0 # orig: 200.0
</weight>
<weight id='bm_wp_2_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 400.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_2_kl_1'>set param_card decay 34 0.1257875 # orig: 0.05214011
set param_card mass 34 400.0 # orig: 200.0
</weight>
<weight id='bm_wp_3_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 500.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_3_kl_1'>set param_card decay 34 0.1614136 # orig: 0.05214011
set param_card mass 34 500.0 # orig: 200.0
</weight>
<weight id='bm_wp_4_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 600.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_4_kl_1'>set param_card decay 34 0.1964585 # orig: 0.05214011
set param_card mass 34 600.0 # orig: 200.0
</weight>
<weight id='bm_wp_5_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 700.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_5_kl_1'>set param_card decay 34 0.2311548 # orig: 0.05214011
set param_card mass 34 700.0 # orig: 200.0
</weight>
<weight id='bm_wp_6_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 800.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_6_kl_1'>set param_card decay 34 0.2656292 # orig: 0.05214011
set param_card mass 34 800.0 # orig: 200.0
</weight>
<weight id='bm_wp_7_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 900.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_7_kl_1'>set param_card decay 34 0.2999543 # orig: 0.05214011
set param_card mass 34 900.0 # orig: 200.0
</weight>
<weight id='bm_wp_8_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1000.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_8_kl_1'>set param_card decay 34 0.3341744 # orig: 0.05214011
set param_card mass 34 1000.0 # orig: 200.0
</weight>
<weight id='bm_wp_9_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1100.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_9_kl_1'>set param_card decay 34 0.368318 # orig: 0.05214011
set param_card mass 34 1100.0 # orig: 200.0
</weight>
<weight id='bm_wp_10_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1200.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_10_kl_1'>set param_card decay 34 0.4024041 # orig: 0.05214011
set param_card mass 34 1200.0 # orig: 200.0
</weight>
<weight id='bm_wp_11_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1300.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_11_kl_1'>set param_card decay 34 0.4364459 # orig: 0.05214011
set param_card mass 34 1300.0 # orig: 200.0
</weight>
<weight id='bm_wp_12_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1400.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_12_kl_1'>set param_card decay 34 0.4704529 # orig: 0.05214011
set param_card mass 34 1400.0 # orig: 200.0
</weight>
<weight id='bm_wp_13_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1500.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_13_kl_1'>set param_card decay 34 0.5044321 # orig: 0.05214011
set param_card mass 34 1500.0 # orig: 200.0
</weight>
<weight id='bm_wp_14_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1600.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_14_kl_1'>set param_card decay 34 0.5383886 # orig: 0.05214011
set param_card mass 34 1600.0 # orig: 200.0
</weight>
<weight id='bm_wp_15_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1700.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_15_kl_1'>set param_card decay 34 0.5723265 # orig: 0.05214011
set param_card mass 34 1700.0 # orig: 200.0
</weight>
<weight id='bm_wp_16_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1800.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_16_kl_1'>set param_card decay 34 0.6062488 # orig: 0.05214011
set param_card mass 34 1800.0 # orig: 200.0
</weight>
<weight id='bm_wp_17_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 1900.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_17_kl_1'>set param_card decay 34 0.6401581 # orig: 0.05214011
set param_card mass 34 1900.0 # orig: 200.0
</weight>
<weight id='bm_wp_18_kl_0'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 2000.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
<weight id='bm_wp_18_kl_1'>set param_card decay 34 0.6740562 # orig: 0.05214011
set param_card mass 34 2000.0 # orig: 200.0
</weight>
<weight id='sm'>set param_card decay 34 0.0 # orig: 0.05214011
set param_card mass 34 0.0 # orig: 200.0
set param_card ssmcoup 1 0.0 # orig: 0.1
</weight>
</weightgroup>
</initrwgt>
<MGGenerationInfo>
# Number of Events : 100
# Integrated weight (pb) : 14.206661267
</MGGenerationInfo>
</header>
<init>
2212 2212 4.000000e+03 4.000000e+03 0 0 261000 261000 -4 3
5.006986e+00 5.571032e-02 1.423070e+01 1
5.790800e+00 1.140000e-01 1.423070e+01 0
3.408875e+00 4.128236e-02 1.423070e+01 2
<generator name='MadGraph5_aMC@NLO' version='2.8.0'>please cite 1405.0301 </generator>
</init>
<event>
 7 2 +1.4230697e+01 2.11223000e+02 7.81616400e-03 1.04745400e-01
        2 -1 0 0 503 0 +0.0000000000e+00 +0.0000000000e+00 +2.0012541118e+01 2.0012541118e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       -1 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00 -1.2035386873e+03 1.2035386873e+03 0.0000000000e+00 0.0000e+00 1.0000e+00
       34 2 1 2 0 0 +1.6829404362e+01 +6.5822578687e+01 -6.1693771192e+02 6.5209455804e+02 1.9999832880e+02 0.0000e+00 0.0000e+00
      -13 1 3 3 0 0 +4.7425607682e+01 -1.9251122775e+01 -5.0780451581e+02 5.1037752720e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       14 1 3 3 0 0 -3.0596203320e+01 +8.5073701462e+01 -1.0913319611e+02 1.4171703084e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 1 1 2 502 501 -1.9628822946e+01 -1.9107179167e+01 -2.3184442362e+02 2.3345708759e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       21 1 1 2 503 502 +2.7994185836e+00 -4.6715399520e+01 -3.3474401063e+02 3.3799958277e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000" pt_clust_6="23.26150" pt_clust_7="74.43301"></scales>
<clustering>
<clus scale=" 23.262"> 6 7 -1 -1</clus>
<clus scale=" 74.433"> 2 6 1 -1</clus>
<clus scale=" 199.998"> 4 5 -1 -1</clus>
<clus scale=" 211.223"> 1 4 2 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +6.2152602e-07 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +1.0790552e-07 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +3.5234624e-08 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +1.5174320e-08 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +7.6733443e-09 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +4.3162594e-09 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +2.6207694e-09 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +1.6860402e-09 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +1.1351122e-09 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +7.9278335e-10 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +5.7074586e-10 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +4.2151035e-10 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +3.1814580e-10 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +2.4468632e-10 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +1.9130269e-10 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +1.5174377e-10 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +1.2191985e-10 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +9.9087869e-11 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>
<event>
 7 2 +1.4230697e+01 2.04216600e+02 7.81616400e-03 1.05219000e-01
        2 -1 0 0 503 0 -0.0000000000e+00 +0.0000000000e+00 +1.5518175205e+02 1.5518175205e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 -1 0 0 501 502 +0.0000000000e+00 -0.0000000000e+00 -8.4729845901e+02 8.4729845901e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
       34 2 1 2 0 0 -3.4314084803e+01 +4.5688544584e+00 +7.9141782844e+01 2.1796425734e+02 2.0011662809e+02 0.0000e+00 0.0000e+00
      -13 1 3 3 0 0 -1.0792715114e+02 -2.2003919441e+01 +8.5510076857e+01 1.3944323457e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       14 1 3 3 0 0 +7.3613066334e+01 +2.6572773899e+01 -6.3682940135e+00 7.8521022767e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 1 1 2 503 502 +4.9368949984e+01 +4.0313421501e+00 -8.7575764788e+01 1.0061341620e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
        1 1 1 2 501 0 -1.5054865181e+01 -8.6001966085e+00 -6.8368272502e+02 6.8390253752e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000" pt_clust_6="40.71550" pt_clust_7="17.33818"></scales>
<clustering>
<clus scale=" 17.338"> 2 7 1 -1</clus>
<clus scale=" 40.716"> 2 6 1 -1</clus>
<clus scale=" 200.117"> 4 5 -1 -1</clus>
<clus scale=" 17.338"> 2 4 1 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +1.3038869e-05 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +2.2612277e-06 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +7.3811393e-07 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +3.1783058e-07 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +1.6070671e-07 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +9.0392838e-08 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +5.4883293e-08 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +3.5307643e-08 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +2.3770148e-08 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +1.6601291e-08 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +1.1951589e-08 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +8.8264817e-09 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +6.6619714e-09 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +5.1237033e-09 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +4.0058388e-09 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +3.1774716e-09 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +2.5529593e-09 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +2.0748602e-09 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>
<event>
 6 1 +1.4230697e+01 2.01238500e+02 7.81616400e-03 1.05426600e-01
        2 -1 0 0 502 0 +0.0000000000e+00 +0.0000000000e+00 +6.6604724203e+01 6.6604724203e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       -1 -1 0 0 0 501 -0.0000000000e+00 -0.0000000000e+00 -1.7408197191e+02 1.7408197191e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       34 2 1 2 0 0 +6.9366597796e+00 +1.2275799397e+01 -9.9649531897e+01 2.2455952542e+02 2.0074396317e+02 0.0000e+00 0.0000e+00
      -13 1 3 3 0 0 -4.5515867245e+01 -3.1987109618e+01 -1.3596867813e+02 1.4690932845e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       14 1 3 3 0 0 +5.2452527024e+01 +4.4262909015e+01 +3.6319146238e+01 7.7650196968e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 1 1 2 502 501 -6.9366597796e+00 -1.2275799397e+01 -7.8277158101e+00 1.6127170692e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000" pt_clust_6="14.10009"></scales>
<clustering>
<clus scale=" 14.100"> 2 6 1 -1</clus>
<clus scale=" 200.744"> 4 5 -1 -1</clus>
<clus scale=" 201.239"> 1 4 2 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +0.0000000e+00 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>
<event>
 6 1 +1.4230697e+01 2.00752900e+02 7.81616400e-03 1.05460800e-01
       -1 -1 0 0 0 501 -0.0000000000e+00 +0.0000000000e+00 +1.2593499518e+02 1.2593499518e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
        2 -1 0 0 502 0 +0.0000000000e+00 -0.0000000000e+00 -9.6574492046e+01 9.6574492046e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       34 2 1 2 0 0 -1.8125604777e+01 -3.8709671018e-01 +1.9054377224e+01 2.0165514418e+02 1.9993258987e+02 0.0000e+00 0.0000e+00
      -13 1 3 3 0 0 -4.8991841477e+01 +8.3056618277e+01 +4.8500499189e+01 1.0793933848e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       14 1 3 3 0 0 +3.0866236700e+01 -8.3443714987e+01 -2.9446121965e+01 9.3715805698e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       21 1 1 2 502 501 +1.8125604777e+01 +3.8709671018e-01 +1.0306125907e+01 2.0854343039e+01 0.0000000000e+00 0.0000e+00 1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000" pt_clust_6="18.12974"></scales>
<clustering>
<clus scale=" 18.130"> 1 6 2 -1</clus>
<clus scale=" 199.933"> 4 5 -1 -1</clus>
<clus scale=" 200.753"> 2 4 1 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +0.0000000e+00 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>
<event>
 7 2 +1.4230697e+01 2.01764600e+02 7.81616400e-03 1.05389700e-01
       21 -1 0 0 501 502 +0.0000000000e+00 +0.0000000000e+00 +4.7352905079e+02 4.7352905079e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
        1 -1 0 0 503 0 -0.0000000000e+00 -0.0000000000e+00 -5.5849953190e+02 5.5849953190e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
      -34 2 1 2 0 0 -3.2681765762e+00 -2.6965754302e+01 -5.3915576520e+02 5.7569388869e+02 1.9999170342e+02 0.0000e+00 0.0000e+00
       13 1 3 3 0 0 +4.5029882175e+01 -4.2414231170e+01 -5.1643689886e+02 5.2012856852e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
      -14 1 3 3 0 0 -4.8298058751e+01 +1.5448476868e+01 -2.2718866344e+01 5.5565320161e+01 0.0000000000e+00 0.0000e+00 1.0000e+00
       21 1 1 2 503 502 +7.1162127269e+00 +1.0404799454e+01 +4.3294162788e+01 4.5091960113e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
        2 1 1 2 501 0 -3.8480361507e+00 +1.6560954848e+01 +4.1089112130e+02 4.1124273388e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000" pt_clust_6="12.60557" pt_clust_7="26.68826"></scales>
<clustering>
<clus scale=" 12.606"> 1 6 2 -1</clus>
<clus scale=" 26.688"> 1 7 2 -1</clus>
<clus scale=" 199.992"> 4 5 -1 -1</clus>
<clus scale=" 12.606"> 1 4 2 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +6.8159704e-07 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +1.1834200e-07 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +3.8643201e-08 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +1.6642419e-08 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +8.4157725e-09 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +4.7338895e-09 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +2.8743538e-09 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +1.8491833e-09 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +1.2449481e-09 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +8.6949541e-10 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +6.2597327e-10 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +4.6229742e-10 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +3.4893102e-10 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +2.6836335e-10 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +2.0981411e-10 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +1.6642730e-10 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +1.3371748e-10 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +1.0867617e-10 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>
<event>
 5 0 +1.4230697e+01 2.00016800e+02 7.81616400e-03 1.05512900e-01
       -1 -1 0 0 0 501 -0.0000000000e+00 +0.0000000000e+00 +1.2371104561e+02 1.2371104561e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
        2 -1 0 0 501 0 +0.0000000000e+00 -0.0000000000e+00 -8.0847068511e+01 8.0847068511e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
       34 2 1 2 0 0 +0.0000000000e+00 +0.0000000000e+00 +4.2863977103e+01 2.0455811413e+02 2.0001675310e+02 0.0000e+00 0.0000e+00
      -13 1 3 3 0 0 -3.2905743832e+01 +1.9938884667e+01 +1.1583894275e+02 1.2206149171e+02 0.0000000000e+00 0.0000e+00 1.0000e+00
       14 1 3 3 0 0 +3.2905743832e+01 -1.9938884667e+01 -7.2974965652e+01 8.2496622420e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
<scales pt_clust_4="8000.00000" pt_clust_5="8000.00000"></scales>
<clustering>
<clus scale=" 200.017"> 4 5 -1 -1</clus>
<clus scale=" 200.017"> 2 4 1 -1</clus>
</clustering>
<rwgt>
<wgt id='bm_wp_0_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_1_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_2_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_3_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_4_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_5_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_6_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_7_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_8_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_9_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_10_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_11_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_12_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_13_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_14_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_15_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_16_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_17_kl_1'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_0'> +0.0000000e+00 </wgt>
<wgt id='bm_wp_18_kl_1'> +0.0000000e+00 </wgt>
<wgt id='sm'> +0.0000000e+00 </wgt>
</rwgt>
</event>

and so on...

If needed, I can also post the log file...which seems to have nothing out of the usual.

Sincerely,

Simonas

Thanks,

Could you paste the reweight_card in top of that?
and also send me a link to the model (if not a public model you can send it to my email address: <email address hidden>

Cheers,

Olivier

Simonas Drauksas (tehmeh) said : #4

No problem, here is the link to the model:
http://feynrules.irmp.ucl.ac.be/wiki/WZPrimeAtNLO

And here is the reweight_card

# Reweight card generated by MadMiner

# Global setup
change output default
change helicity False

# MadMiner benchmark bm_wp_0_kl_0
launch --rwgt_name=bm_wp_0_kl_0
  set mass 34 200
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_1_kl_0
launch --rwgt_name=bm_wp_1_kl_0
  set mass 34 300
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_1_kl_1
launch --rwgt_name=bm_wp_1_kl_1
  set mass 34 300
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_2_kl_0
launch --rwgt_name=bm_wp_2_kl_0
  set mass 34 400
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_2_kl_1
launch --rwgt_name=bm_wp_2_kl_1
  set mass 34 400
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_3_kl_0
launch --rwgt_name=bm_wp_3_kl_0
  set mass 34 500
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_3_kl_1
launch --rwgt_name=bm_wp_3_kl_1
  set mass 34 500
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_4_kl_0
launch --rwgt_name=bm_wp_4_kl_0
  set mass 34 600
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_4_kl_1
launch --rwgt_name=bm_wp_4_kl_1
  set mass 34 600
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_5_kl_0
launch --rwgt_name=bm_wp_5_kl_0
  set mass 34 700
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_5_kl_1
launch --rwgt_name=bm_wp_5_kl_1
  set mass 34 700
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_6_kl_0
launch --rwgt_name=bm_wp_6_kl_0
  set mass 34 800
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_6_kl_1
launch --rwgt_name=bm_wp_6_kl_1
  set mass 34 800
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_7_kl_0
launch --rwgt_name=bm_wp_7_kl_0
  set mass 34 900
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_7_kl_1
launch --rwgt_name=bm_wp_7_kl_1
  set mass 34 900
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_8_kl_0
launch --rwgt_name=bm_wp_8_kl_0
  set mass 34 1000
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_8_kl_1
launch --rwgt_name=bm_wp_8_kl_1
  set mass 34 1000
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_9_kl_0
launch --rwgt_name=bm_wp_9_kl_0
  set mass 34 1100
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_9_kl_1
launch --rwgt_name=bm_wp_9_kl_1
  set mass 34 1100
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_10_kl_0
launch --rwgt_name=bm_wp_10_kl_0
  set mass 34 1200
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_10_kl_1
launch --rwgt_name=bm_wp_10_kl_1
  set mass 34 1200
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_11_kl_0
launch --rwgt_name=bm_wp_11_kl_0
  set mass 34 1300
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_11_kl_1
launch --rwgt_name=bm_wp_11_kl_1
  set mass 34 1300
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_12_kl_0
launch --rwgt_name=bm_wp_12_kl_0
  set mass 34 1400
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_12_kl_1
launch --rwgt_name=bm_wp_12_kl_1
  set mass 34 1400
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_13_kl_0
launch --rwgt_name=bm_wp_13_kl_0
  set mass 34 1500
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_13_kl_1
launch --rwgt_name=bm_wp_13_kl_1
  set mass 34 1500
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_14_kl_0
launch --rwgt_name=bm_wp_14_kl_0
  set mass 34 1600
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_14_kl_1
launch --rwgt_name=bm_wp_14_kl_1
  set mass 34 1600
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_15_kl_0
launch --rwgt_name=bm_wp_15_kl_0
  set mass 34 1700
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_15_kl_1
launch --rwgt_name=bm_wp_15_kl_1
  set mass 34 1700
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_16_kl_0
launch --rwgt_name=bm_wp_16_kl_0
  set mass 34 1800
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_16_kl_1
launch --rwgt_name=bm_wp_16_kl_1
  set mass 34 1800
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_17_kl_0
launch --rwgt_name=bm_wp_17_kl_0
  set mass 34 1900
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_17_kl_1
launch --rwgt_name=bm_wp_17_kl_1
  set mass 34 1900
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark bm_wp_18_kl_0
launch --rwgt_name=bm_wp_18_kl_0
  set mass 34 2000
  set ssmcoup 1 0.0
  set decay 34 auto

# MadMiner benchmark bm_wp_18_kl_1
launch --rwgt_name=bm_wp_18_kl_1
  set mass 34 2000
  set ssmcoup 1 0.1
  set decay 34 auto

# MadMiner benchmark sm
launch --rwgt_name=sm
  set mass 34 0
  set ssmcoup 1 0.0
  set decay 34 auto

Sincerely,

Simonas

Thanks I can reproduce your run now.

Will report what I can find.

Thanks,

Olivier

Hi,

I think that the issue is related to some optimization with the helicity filtering which are not correctly reset when changing benchmark. I'm working on a patch.

Cheers,

Olivier

Hi,

You will find the patch (and a full tarball) fixing the issue here:
https://bazaar.launchpad.net/~maddevelopers/mg5amcnlo/2.8.3.py2/revision/295

Note that such version is python2 only which should not be a problem for you.

Cheers,

Olivier

Simonas Drauksas (tehmeh) said : #8

Thanks, this works now!