Produce vector boson strictly on-shell

Asked by Harikrishnan Ramani

Hello authors,
I am interested in a pp->ww->eevv process.
hence i generate p p > w+ w- , w+ > e+ ve , w- > e- ve~ as suggested.
I understand that toggling with bw parameter in runcard can reduce the bw width M+- bw x gamma. However for my application I need the Intermediate W masses at strictly mW. Is there anyway i could do this inside madgraph? setting bw=0 gives me zero cross section because madgraph merely trims the tail and 0 phase space is left.

I understand I can do p p > w+ w- and then shower/decay it with pythia. but I would prefer to do the whole thing inside madgraph

Thank You,
Harikrishnan

Question information

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

Hi Harikrishnan,

I would propose to use Madspin for the decay where you can switch off the BW effect (require MG5 v2.0.0 or newer).
So you generate
p p > w+ w-
then ask Madspin for the decay
and edit the madspin_card to have only the requested electron decay
and in that file you set the following options:
set BW_effect False

Another option, is to do the same with the DECAY package but this is less precise.

Cheers,

Olivier

On Mar 6, 2014, at 1:51 AM, Harikrishnan Ramani <email address hidden> wrote:

> New question #245046 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/245046
>
> Hello authors,
> I am interested in a pp->ww->eevv process.
> hence i generate p p > w+ w- , w+ > e+ ve , w- > e- ve~ as suggested.
> I understand that toggling with bw parameter in runcard can reduce the bw width M+- bw x gamma. However for my application I need the Intermediate W masses at strictly mW. Is there anyway i could do this inside madgraph? setting bw=0 gives me zero cross section because madgraph merely trims the tail and 0 phase space is left.
>
> I understand I can do p p > w+ w- and then shower/decay it with pythia. but I would prefer to do the whole thing inside madgraph
>
> Thank You,
> Harikrishnan
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Harikrishnan Ramani (harikrishramani) said :
#2

Hi Olivier,
Thank you for your proposal. However, it still produces Ws off-shell.
Copied the lhe file including first event for your perusal.
Thank you.
<LesHouchesEvents version="1.0">
<init>
  2212 2212 0.40000E+04 0.40000E+04 0 0 0 0 3 1
 0.45190E+09 0.00000E+00 0.10000E+01 100
</init>
<!--
# Warning! Never use this file for detector studies!
-->
<!--

<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.1.0
</mgversion>
<mg5proccard>
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.1.0 2014-02-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 loop_optimized_output True
set gauge unitary
set complex_mass_scheme False
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~
generate p p > w+ w-
output likecms
</mg5proccard>
<mgproccard>
#*********************************************************************
# 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 > w+ w- #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
sm
# 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>
#*********************************************************************
# 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 *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Run to generate the grid pack *
#*********************************************************************
  .false. = gridpack !True = setting up the grid pack
#*********************************************************************
# 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. *
#*********************************************************************
     11 = nevents ! Number of unweighted events requested
 189 = 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 = ebeam1 ! beam 1 total energy in GeV
     4000 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
        0 = polbeam1 ! beam polarization for beam 1
        0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
 'cteq6l1' = pdlabel ! PDF set
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 T = fixed_ren_scale ! if .true. use fixed ren scale
 T = fixed_fac_scale ! if .true. use fixed fac scale
 160.770 = scale ! fixed ren scale
 160.770 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 160.770 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 1 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Matching - Warning! ickkw > 1 is still beta
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
 1 = highestmult ! for ickkw=2, highest mult group
 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
 1 = alpsfact ! scale factor for QCD emission vx
 F = chcluster ! cluster only according to channel diag
 T = pdfwgt ! for ickkw=1, perform pdf reweighting
 5 = asrwgtflavor ! highest quark flavor for a_s reweight
 T = clusinfo ! include clustering tag in output
#*********************************************************************
#**********************************************************
#
#**********************************************************
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#**********************************************************
   T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#**********************************************************
#
#**********************************
# BW cutoff (M+/-bwcutoff*Gamma)
#**********************************
  15 = bwcutoff ! (M+/-bwcutoff*Gamma)
#**********************************************************
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#**********************************************************
   T = cut_decays ! Cut decay products
#*************************************************************
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#*************************************************************
   0 = nhel ! Number of helicities used per event
#*******************
# Standard Cuts
#*******************
#
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
  0 = ptj ! minimum pt for the jets
  0 = ptb ! minimum pt for the b
  0 = pta ! minimum pt for the photons
  0 = ptl ! minimum pt for the charged leptons
  0 = misset ! minimum missing Et (sum of neutrino's momenta)
  0 = ptheavy ! minimum pt for one heavy final state
  0 = ptonium ! minimum pt for the quarkonium states
 -1 = ptjmax ! maximum pt for the jets
 -1 = ptbmax ! maximum pt for the b
 -1 = ptamax ! maximum pt for the photons
 -1 = ptlmax ! maximum pt for the charged leptons
 -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#*********************************************************************
# Minimum and maximum E's (in the lab frame) *
#*********************************************************************
  0 = ej ! minimum E for the jets
  0 = eb ! minimum E for the b
  0 = ea ! minimum E for the photons
  0 = el ! minimum E for the charged leptons
 -1 = ejmax ! maximum E for the jets
 -1 = ebmax ! maximum E for the b
 -1 = eamax ! maximum E for the photons
 -1 = elmax ! maximum E for the charged leptons
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
  -1 = etaj ! max rap for the jets
  -1 = etab ! max rap for the b
  -1 = etaa ! max rap for the photons
  -1 = etal ! max rap for the charged leptons
  -1 = etaonium ! max rap for the quarkonium states
   0 = etajmin ! min rap for the jets
   0 = etabmin ! min rap for the b
   0 = etaamin ! min rap for the photons
   0 = etalmin ! main rap for the charged leptons
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.4 = drjj ! min distance between jets
 0 = drbb ! min distance between b's
 0.4 = drll ! min distance between leptons
 0.4 = draa ! min distance between gammas
 0 = 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 = drab ! min distance between gamma and b
 0 = drbl ! min distance between b and lepton
 0.4 = dral ! min distance between gamma and lepton
 -1 = drjjmax ! max distance between jets
 -1 = drbbmax ! max distance between b's
 -1 = drllmax ! max distance between leptons
 -1 = draamax ! max distance between gammas
 -1 = drbjmax ! max distance between b and jet
 -1 = drajmax ! max distance between gamma and jet
 -1 = drjlmax ! max distance between jet and lepton
 -1 = drabmax ! max distance between gamma and b
 -1 = drblmax ! max distance between b and lepton
 -1 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 0 = mmjj ! min invariant mass of a jet pair
 0 = mmbb ! min invariant mass of a b pair
 0 = mmaa ! min invariant mass of gamma gamma pair
 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1 = mmjjmax ! max invariant mass of a jet pair
 -1 = mmbbmax ! max invariant mass of a b pair
 -1 = mmaamax ! max invariant mass of gamma gamma pair
 -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0 = xptj ! minimum pt for at least one jet
 0 = xptb ! minimum pt for at least one b
 0 = xpta ! minimum pt for at least one photon
 0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0 = ptj1min ! minimum pt for the leading jet in pt
 0 = ptj2min ! minimum pt for the second jet in pt
 0 = ptj3min ! minimum pt for the third jet in pt
 0 = ptj4min ! minimum pt for the fourth jet in pt
 -1 = ptj1max ! maximum pt for the leading jet in pt
 -1 = ptj2max ! maximum pt for the second jet in pt
 -1 = ptj3max ! maximum pt for the third jet in pt
 -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
 0 = ptl2min ! minimum pt for the second lepton in pt
 0 = ptl3min ! minimum pt for the third lepton in pt
 0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1 = ptl1max ! maximum pt for the leading lepton in pt
 -1 = ptl2max ! maximum pt for the second lepton in pt
 -1 = ptl3max ! maximum pt for the third lepton in pt
 -1 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1 = htjmax ! maximum jet HT=Sum(jet pt)
 0 = ihtmin !inclusive Ht for all partons (including b)
 -1 = ihtmax !inclusive Ht for all partons (including b)
 0 = ht2min ! minimum Ht for the two leading jets
 0 = ht3min ! minimum Ht for the three leading jets
 0 = ht4min ! minimum Ht for the four leading jets
 -1 = ht2max ! maximum Ht for the two leading jets
 -1 = ht3max ! maximum Ht for the three leading jets
 -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
 0 = deltaeta ! minimum rapidity for two jets in the WBF case
#*********************************************************************
# KT DURHAM CUT *
#*********************************************************************
 -1 = ktdurham
 0.4 = dparameter
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
# Jet measure cuts *
#*********************************************************************
 0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: If use_syst is T, matched Pythia output is *
# meaningful ONLY if plotted taking matchscale *
# reweighting into account! *
#*********************************************************************
   F = use_syst ! Enable systematics studies
#
#**************************************
# Parameter of the systematics study
# will be use by SysCalc (if install)
#**************************************
#
0.5 1 2 = sys_scalefact # Central scale factors
0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
30 50 = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
CT10nlo.LHgrid = sys_pdf # matching scales
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf
</mgruncard>
<slha>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.1187600e+01 # MZ
   25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 80.38500 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.322330e+02 # aEWM1
    2 1.166380e-05 # Gf
    3 1.180000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.700000e+00 # ymb
    6 1.730000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
</slha>
<MGGenerationInfo>
# Number of Events : 11
# Integrated weight (pb) : 35.7846
</MGGenerationInfo>
<mgpythiacard>
!...Parton showering on or off
      MSTP(61)=1
      MSTP(71)=1

!...Fragmentation/hadronization on or off
      MSTJ(1)=1

!...Multiple interactions on or off
      MSTP(81)=20

!...Don't stop execution after 10 errors
      MSTU(21)=1

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets

      LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
</mgpythiacard>
<madspin>
#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
 set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set BW_effect False
#
# specify the decay for the final state particles
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > e+ ve
decay w- > e- ve~
decay z > all all
# running the actual code
launch
</madspin>
</header>
</LesHouchesEvents>

-->
<event>
 11 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
      -2 -1 0 0 0 0 -0.54179131E+01 -0.25883280E+01 0.50677325E+02 0.51031799E+02 0.00000000E+00 0. 0.
       2 -1 0 0 0 0 -0.27751624E+00 0.35994440E+01 -0.16756743E+03 0.16760631E+03 0.00000000E+00 0. 0.
      24 2 1 1 0 0 0.35950721E+02 0.25561351E+02 -0.81317312E+02 0.11809097E+03 0.73397111E+02 0. 0.
     -24 2 1 1 0 0 -0.41646150E+02 -0.24550236E+02 -0.35572794E+02 0.10054714E+03 0.80667143E+02 0. 0.
     -11 2 3 3 0 0 -0.12415025E+02 -0.16554797E+02 -0.27197544E+02 0.34174560E+02 0.51000000E-03 0. 0.
      12 2 3 3 0 0 0.48365745E+02 0.42116148E+02 -0.54119768E+02 0.83916414E+02 0.00000000E+00 0. 0.
      11 2 4 4 0 0 -0.46063140E+02 -0.65992536E-01 -0.52963448E+02 0.70192194E+02 0.51000000E-03 0. 0.
     -12 2 4 4 0 0 0.44169898E+01 -0.24484243E+02 0.17390654E+02 0.30354947E+02 0.00000000E+00 0. 0.
     -11 1 0 0 0 0 -0.12921121E+02 -0.14762599E+02 -0.26326985E+02 0.32832908E+02 0.51000000E-03 0. 0.
      11 1 0 0 0 0 -0.46027223E+02 -0.65895362E-01 -0.52922065E+02 0.70137398E+02 0.51000000E-03 0. 0.
      12 1 0 0 0 0 0.52279052E+02 0.17193306E+02 0.00000000E+00 0.55033708E+02 0.00000000E+00 0. 0.
</event>

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

I’m adding Pierre in cc.
It might better know those options than me.
The second possibility is probably to have the line
 set BW_cut 0

In the madspin_card.

Cheers,

Olivier
On Mar 6, 2014, at 8:31 PM, Harikrishnan Ramani <email address hidden> wrote:

> Question #245046 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/245046
>
> Status: Answered => Open
>
> Harikrishnan Ramani is still having a problem:
> Hi Olivier,
> Thank you for your proposal. However, it still produces Ws off-shell.
> Copied the lhe file including first event for your perusal.
> Thank you.
> <LesHouchesEvents version="1.0">
> <init>
> 2212 2212 0.40000E+04 0.40000E+04 0 0 0 0 3 1
> 0.45190E+09 0.00000E+00 0.10000E+01 100
> </init>
> <!--
> # Warning! Never use this file for detector studies!
> -->
> <!--
>
> <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.1.0
> </mgversion>
> <mg5proccard>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.1.0 2014-02-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 loop_optimized_output True
> set gauge unitary
> set complex_mass_scheme False
> 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~
> generate p p > w+ w-
> output likecms
> </mg5proccard>
> <mgproccard>
> #*********************************************************************
> # 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 > w+ w- #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
> sm
> # 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>
> #*********************************************************************
> # 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 *
> #*********************************************************************
> #
> #*******************
> # Running parameters
> #*******************
> #
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_1 = run_tag ! name of the run
> #*********************************************************************
> # Run to generate the grid pack *
> #*********************************************************************
> .false. = gridpack !True = setting up the grid pack
> #*********************************************************************
> # 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. *
> #*********************************************************************
> 11 = nevents ! Number of unweighted events requested
> 189 = 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 = ebeam1 ! beam 1 total energy in GeV
> 4000 = ebeam2 ! beam 2 total energy in GeV
> #*********************************************************************
> # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> #*********************************************************************
> 0 = polbeam1 ! beam polarization for beam 1
> 0 = polbeam2 ! beam polarization for beam 2
> #*********************************************************************
> # PDF CHOICE: this automatically fixes also alpha_s and its evol. *
> #*********************************************************************
> 'cteq6l1' = pdlabel ! PDF set
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> T = fixed_ren_scale ! if .true. use fixed ren scale
> T = fixed_fac_scale ! if .true. use fixed fac scale
> 160.770 = scale ! fixed ren scale
> 160.770 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 160.770 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> 1 = scalefact ! scale factor for event-by-event scales
> #*********************************************************************
> # Matching - Warning! ickkw > 1 is still beta
> #*********************************************************************
> 0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
> 1 = highestmult ! for ickkw=2, highest mult group
> 1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
> 1 = alpsfact ! scale factor for QCD emission vx
> F = chcluster ! cluster only according to channel diag
> T = pdfwgt ! for ickkw=1, perform pdf reweighting
> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
> T = clusinfo ! include clustering tag in output
> #*********************************************************************
> #**********************************************************
> #
> #**********************************************************
> # Automatic ptj and mjj cuts if xqcut > 0
> # (turn off for VBF and single top processes)
> #**********************************************************
> T = auto_ptj_mjj ! Automatic setting of ptj and mjj
> #**********************************************************
> #
> #**********************************
> # BW cutoff (M+/-bwcutoff*Gamma)
> #**********************************
> 15 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #**********************************************************
> # Apply pt/E/eta/dr/mij cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #**********************************************************
> T = cut_decays ! Cut decay products
> #*************************************************************
> # Number of helicities to sum per event (0 = all helicities)
> # 0 gives more stable result, but longer run time (needed for
> # long decay chains e.g.).
> # Use >=2 if most helicities contribute, e.g. pure QCD.
> #*************************************************************
> 0 = nhel ! Number of helicities used per event
> #*******************
> # Standard Cuts
> #*******************
> #
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 0 = ptj ! minimum pt for the jets
> 0 = ptb ! minimum pt for the b
> 0 = pta ! minimum pt for the photons
> 0 = ptl ! minimum pt for the charged leptons
> 0 = misset ! minimum missing Et (sum of neutrino's momenta)
> 0 = ptheavy ! minimum pt for one heavy final state
> 0 = ptonium ! minimum pt for the quarkonium states
> -1 = ptjmax ! maximum pt for the jets
> -1 = ptbmax ! maximum pt for the b
> -1 = ptamax ! maximum pt for the photons
> -1 = ptlmax ! maximum pt for the charged leptons
> -1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> #*********************************************************************
> # Minimum and maximum E's (in the lab frame) *
> #*********************************************************************
> 0 = ej ! minimum E for the jets
> 0 = eb ! minimum E for the b
> 0 = ea ! minimum E for the photons
> 0 = el ! minimum E for the charged leptons
> -1 = ejmax ! maximum E for the jets
> -1 = ebmax ! maximum E for the b
> -1 = eamax ! maximum E for the photons
> -1 = elmax ! maximum E for the charged leptons
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> -1 = etaj ! max rap for the jets
> -1 = etab ! max rap for the b
> -1 = etaa ! max rap for the photons
> -1 = etal ! max rap for the charged leptons
> -1 = etaonium ! max rap for the quarkonium states
> 0 = etajmin ! min rap for the jets
> 0 = etabmin ! min rap for the b
> 0 = etaamin ! min rap for the photons
> 0 = etalmin ! main rap for the charged leptons
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.4 = drjj ! min distance between jets
> 0 = drbb ! min distance between b's
> 0.4 = drll ! min distance between leptons
> 0.4 = draa ! min distance between gammas
> 0 = 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 = drab ! min distance between gamma and b
> 0 = drbl ! min distance between b and lepton
> 0.4 = dral ! min distance between gamma and lepton
> -1 = drjjmax ! max distance between jets
> -1 = drbbmax ! max distance between b's
> -1 = drllmax ! max distance between leptons
> -1 = draamax ! max distance between gammas
> -1 = drbjmax ! max distance between b and jet
> -1 = drajmax ! max distance between gamma and jet
> -1 = drjlmax ! max distance between jet and lepton
> -1 = drabmax ! max distance between gamma and b
> -1 = drblmax ! max distance between b and lepton
> -1 = dralmax ! maxdistance between gamma and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> 0 = mmjj ! min invariant mass of a jet pair
> 0 = mmbb ! min invariant mass of a b pair
> 0 = mmaa ! min invariant mass of gamma gamma pair
> 0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1 = mmjjmax ! max invariant mass of a jet pair
> -1 = mmbbmax ! max invariant mass of a b pair
> -1 = mmaamax ! max invariant mass of gamma gamma pair
> -1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons *
> #*********************************************************************
> 0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0 = xptj ! minimum pt for at least one jet
> 0 = xptb ! minimum pt for at least one b
> 0 = xpta ! minimum pt for at least one photon
> 0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0 = ptj1min ! minimum pt for the leading jet in pt
> 0 = ptj2min ! minimum pt for the second jet in pt
> 0 = ptj3min ! minimum pt for the third jet in pt
> 0 = ptj4min ! minimum pt for the fourth jet in pt
> -1 = ptj1max ! maximum pt for the leading jet in pt
> -1 = ptj2max ! maximum pt for the second jet in pt
> -1 = ptj3max ! maximum pt for the third jet in pt
> -1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
> 0 = ptl2min ! minimum pt for the second lepton in pt
> 0 = ptl3min ! minimum pt for the third lepton in pt
> 0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1 = ptl1max ! maximum pt for the leading lepton in pt
> -1 = ptl2max ! maximum pt for the second lepton in pt
> -1 = ptl3max ! maximum pt for the third lepton in pt
> -1 = ptl4max ! maximum pt for the fourth lepton in pt
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1 = htjmax ! maximum jet HT=Sum(jet pt)
> 0 = ihtmin !inclusive Ht for all partons (including b)
> -1 = ihtmax !inclusive Ht for all partons (including b)
> 0 = ht2min ! minimum Ht for the two leading jets
> 0 = ht3min ! minimum Ht for the three leading jets
> 0 = ht4min ! minimum Ht for the four leading jets
> -1 = ht2max ! maximum Ht for the two leading jets
> -1 = ht3max ! maximum Ht for the three leading jets
> -1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
> 0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #*********************************************************************
> # KT DURHAM CUT *
> #*********************************************************************
> -1 = ktdurham
> 0.4 = dparameter
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> # Jet measure cuts *
> #*********************************************************************
> 0 = xqcut ! minimum kt jet measure between partons
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: If use_syst is T, matched Pythia output is *
> # meaningful ONLY if plotted taking matchscale *
> # reweighting into account! *
> #*********************************************************************
> F = use_syst ! Enable systematics studies
> #
> #**************************************
> # Parameter of the systematics study
> # will be use by SysCalc (if install)
> #**************************************
> #
> 0.5 1 2 = sys_scalefact # Central scale factors
> 0.5 1 2 = sys_alpsfact # \alpha_s emission scale factors
> 30 50 = sys_matchscale # variation of merging scale
> # PDF sets and number of members (0 or none for all members).
> CT10nlo.LHgrid = sys_pdf # matching scales
> # MSTW2008nlo68cl.LHgrid 1 = sys_pdf
> </mgruncard>
> <slha>
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.1187600e+01 # MZ
> 25 1.250000e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.000000 # d : 0.0
> 2 0.000000 # u : 0.0
> 3 0.000000 # s : 0.0
> 4 0.000000 # c : 0.0
> 11 0.000000 # e- : 0.0
> 12 0.000000 # ve : 0.0
> 13 0.000000 # mu- : 0.0
> 14 0.000000 # vm : 0.0
> 16 0.000000 # vt : 0.0
> 21 0.000000 # g : 0.0
> 22 0.000000 # a : 0.0
> 24 80.38500 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 1.322330e+02 # aEWM1
> 2 1.166380e-05 # Gf
> 3 1.180000e-01 # aS
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.000000 # d : 0.0
> DECAY 2 0.000000 # u : 0.0
> DECAY 3 0.000000 # s : 0.0
> DECAY 4 0.000000 # c : 0.0
> DECAY 5 0.000000 # b : 0.0
> DECAY 11 0.000000 # e- : 0.0
> DECAY 12 0.000000 # ve : 0.0
> DECAY 13 0.000000 # mu- : 0.0
> DECAY 14 0.000000 # vm : 0.0
> DECAY 15 0.000000 # ta- : 0.0
> DECAY 16 0.000000 # vt : 0.0
> DECAY 21 0.000000 # g : 0.0
> DECAY 22 0.000000 # a : 0.0
> </slha>
> <MGGenerationInfo>
> # Number of Events : 11
> # Integrated weight (pb) : 35.7846
> </MGGenerationInfo>
> <mgpythiacard>
> !...Parton showering on or off
> MSTP(61)=1
> MSTP(71)=1
>
> !...Fragmentation/hadronization on or off
> MSTJ(1)=1
>
> !...Multiple interactions on or off
> MSTP(81)=20
>
> !...Don't stop execution after 10 errors
> MSTU(21)=1
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
>
> LHAPATH=/Users/harikrishnan/Downloads/MG5_aMC_v2_1_0/pythia-pgs/src/PDFsets
> </mgpythiacard>
> <madspin>
> #************************************************************
> #* MadSpin *
> #* *
> #* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
> #* *
> #* Part of the MadGraph5_aMC@NLO Framework: *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #Some options (uncomment to apply)
> #
> # set seed 1
> # set Nevents_for_max_weigth 75 # number of events for the estimate of the max. weight
> # set BW_cut 15 # cut on how far the particle can be off-shell
> set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
> set BW_effect False
> #
> # specify the decay for the final state particles
> decay t > w+ b, w+ > all all
> decay t~ > w- b~, w- > all all
> decay w+ > e+ ve
> decay w- > e- ve~
> decay z > all all
> # running the actual code
> launch
> </madspin>
> </header>
> </LesHouchesEvents>
>
> -->
> <event>
> 11 100 0.10000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00
> -2 -1 0 0 0 0 -0.54179131E+01 -0.25883280E+01 0.50677325E+02 0.51031799E+02 0.00000000E+00 0. 0.
> 2 -1 0 0 0 0 -0.27751624E+00 0.35994440E+01 -0.16756743E+03 0.16760631E+03 0.00000000E+00 0. 0.
> 24 2 1 1 0 0 0.35950721E+02 0.25561351E+02 -0.81317312E+02 0.11809097E+03 0.73397111E+02 0. 0.
> -24 2 1 1 0 0 -0.41646150E+02 -0.24550236E+02 -0.35572794E+02 0.10054714E+03 0.80667143E+02 0. 0.
> -11 2 3 3 0 0 -0.12415025E+02 -0.16554797E+02 -0.27197544E+02 0.34174560E+02 0.51000000E-03 0. 0.
> 12 2 3 3 0 0 0.48365745E+02 0.42116148E+02 -0.54119768E+02 0.83916414E+02 0.00000000E+00 0. 0.
> 11 2 4 4 0 0 -0.46063140E+02 -0.65992536E-01 -0.52963448E+02 0.70192194E+02 0.51000000E-03 0. 0.
> -12 2 4 4 0 0 0.44169898E+01 -0.24484243E+02 0.17390654E+02 0.30354947E+02 0.00000000E+00 0. 0.
> -11 1 0 0 0 0 -0.12921121E+02 -0.14762599E+02 -0.26326985E+02 0.32832908E+02 0.51000000E-03 0. 0.
> 11 1 0 0 0 0 -0.46027223E+02 -0.65895362E-01 -0.52922065E+02 0.70137398E+02 0.51000000E-03 0. 0.
> 12 1 0 0 0 0 0.52279052E+02 0.17193306E+02 0.00000000E+00 0.55033708E+02 0.00000000E+00 0. 0.
> </event>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

Provide an answer of your own, or ask Harikrishnan Ramani for more information if necessary.

To post a message you must log in.