MadGraph fails to generate the requested number of events

Asked by HH on 2020-08-19

Hi

I'm using MG version 2.7.3 and I need to generate some events for a BSM process. I want MG to generate 100k events but it generates only about 20 events. The model is a UFO BSM model and no error arises while importing it to MG. Below you can see MG's output lhe event file when nevents=100. The more weird thing is that sometimes MG generates successfully all the requested events! However, repeating the generation, it fails again and generates about 10-20 events even for nevents=100k.
Your help is highly appreciated.

#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model /home/MyFiles/sw/MG5_aMC_v2_7_3/models/ALP_line\
ar_LFV_UFO
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 jb = j b b~
generate e- e+ > z z weighted<=99, z > jb jb , (z > mu- mu+, mu+ > e+ \
ax)
output zz2llax13-rh
]]>
</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
e- e+ > z z , z > jb jb , (z > mu- mu+ , mu+ > e+ ax) #Process
# Be carefull the coupling are here in MG5 convention
weighted<=99

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
ALP_linear_LFV_UFO
# 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 *
#*********************************************************************
  0 = lpp1 ! beam 1 type
  0 = lpp2 ! beam 2 type
  500.0 = ebeam1 ! beam 1 total energy in GeV
  500.0 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
  0.0 = polbeam1 ! beam polarization for beam 1
  0.0 = polbeam2 ! beam polarization for beam 2

#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
  nn23lo1 = pdlabel ! PDF set
  230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
  False = fixed_ren_scale ! if .true. use fixed ren scale
  False = fixed_fac_scale ! if .true. use fixed fac scale
  91.188 = scale ! fixed ren scale
  91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
  91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
  -1 = dynamical_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
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"

#*********************************************************************
#
#*********************************************************************
# 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
  0.0 = ptb ! minimum pt for the b
  10.0 = ptl ! minimum pt for the charged leptons
  -1.0 = ptjmax ! maximum pt for the jets
  -1.0 = ptbmax ! maximum pt for the b
  -1.0 = ptlmax ! maximum pt for the charged leptons
  {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
  {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0.0 = ej ! minimum E for the jets
  0.0 = eb ! minimum E for the b
  0.0 = ea ! minimum E for the photons
  0.0 = el ! minimum E for the charged leptons
  -1.0 = ejmax ! maximum E for the jets
  -1.0 = ebmax ! maximum E for the b
  -1.0 = eamax ! maximum E for the photons
  -1.0 = elmax ! maximum E for the charged leptons
  {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
  {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})

#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
  5.0 = etaj ! max rap for the jets
  -1.0 = etab ! max rap for the b
  2.5 = etal ! max rap for the charged leptons
  0.0 = etabmin ! min rap for the b
  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.0 = drbb ! min distance between b's
  0.4 = drll ! min distance between leptons
  0.0 = drbj ! min distance between b and jet
  0.4 = drjl ! min distance between jet and lepton
  0.0 = drbl ! min distance between b and lepton
  -1.0 = drjjmax ! max distance between jets
  -1.0 = drbbmax ! max distance between b's
  -1.0 = drllmax ! max distance between leptons
  -1.0 = drbjmax ! max distance between b and jet
  -1.0 = drjlmax ! max distance between jet and lepton
  -1.0 = drblmax ! max distance between b and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
  0.0 = mmjj ! min invariant mass of a jet pair
  0.0 = mmbb ! min invariant mass of a b pair
  0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
  -1.0 = mmjjmax ! max invariant mass of a jet pair
  -1.0 = mmbbmax ! max invariant mass of a b pair
  -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
  {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
  {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
 #*********************************************************************
 # Minimum and maximum invariant mass for all letpons *
 #*********************************************************************
  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 / neutrino *
 # for pair of lepton includes only same flavor, opposite charge
 #*********************************************************************
  0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
  -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
  0.0 = ptheavy ! minimum pt for at least one heavy final state
  0.0 = xptj ! minimum pt for at least one jet
  0.0 = xptb ! minimum pt for at least one b
  0.0 = xptl ! minimum pt for at least one charged lepton
 #*********************************************************************
 # Control the pt's of the jets sorted by pt *
 #*********************************************************************
  0.0 = ptj1min ! minimum pt for the leading jet in pt
  0.0 = ptj2min ! minimum pt for the second jet in pt
  -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 pt's of leptons sorted by pt *
 #*********************************************************************
  0.0 = ptl1min ! minimum pt for the leading lepton in pt
  0.0 = ptl2min ! minimum pt for the second lepton in pt
  -1.0 = ptl1max ! maximum pt for the leading lepton in pt
  -1.0 = ptl2max ! maximum pt for the second lepton in pt
 #*********************************************************************
 # Control the Ht(k)=Sum of k leading jets *
 #*********************************************************************
  0.0 = htjmin ! minimum jet HT=Sum(jet pt)
  -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
  0.0 = ihtmin !inclusive Ht for all partons (including b)
  -1.0 = ihtmax !inclusive Ht for all partons (including b)
 #*********************************************************************
 # 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) *
#*********************************************************************
  4 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
  False = use_syst ! Enable systematics studies
#
  systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'#*********************************************************************
# Additional hidden parameters
#*********************************************************************
  ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments # Choose the argment to pass to the systematics command. like --mur=0.25,1,4. Look at the help of the systematics function for more details.
]]>
</MGRunCard>
<slha>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR ALPPARS
###################################
BLOCK ALPPARS #
      1 1.000000e+03 # fa
      2 0.000000e+00 # cgtil
      3 0.000000e+00 # cwtil
      4 0.000000e+00 # cbtil
      5 0.000000e+00 # caphi
      6 1.000000e+00 # cuni
      7 1.000000e+00 # cvem
      8 1.000000e+00 # cvet
      9 1.000000e+00 # cvmt
      10 1.000000e+00 # caem
      11 1.000000e+00 # caet
      12 1.000000e+00 # camt
###################################
## INFORMATION FOR CKMBLOCK
###################################
BLOCK CKMBLOCK #
      1 2.277360e-01 # cabi
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      1 5.040000e-03 # md
      2 2.550000e-03 # mu
      3 1.010000e-01 # ms
      4 1.270000e+00 # mc
      5 4.700000e+00 # mb
      6 1.720000e+02 # mt
      11 5.110000e-04 # me
      13 1.056600e-01 # mmu
      15 1.777000e+00 # mta
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh
      9000005 1.000000e-03 # ma
      12 0.000000e+00 # ve : 0.0
      14 0.000000e+00 # vm : 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.982436e+01 # 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.279000e+02 # aewm1
      2 1.166370e-05 # gf
      3 1.184000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      1 5.040000e-03 # ymdo
      2 2.550000e-03 # ymup
      3 1.010000e-01 # yms
      4 1.270000e+00 # ymc
      5 4.700000e+00 # ymb
      6 1.720000e+02 # ymt
      11 5.110000e-04 # yme
      13 1.056600e-01 # ymm
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.508336e+00 # wt
DECAY 23 2.495200e+00 # wz
DECAY 24 2.085000e+00 # ww
DECAY 25 4.070000e-03 # wh
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 9000005 0.000000e+00 # ax : 0.0
###################################
## INFORMATION FOR QNUMBERS 9000005
###################################
BLOCK QNUMBERS 9000005 # ax
      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 0 # particle/antiparticle distinction (0=own anti)
</slha>
<MGGenerationInfo>
# Number of Events : 100
# Integrated weight (pb) : 2.77229e-08
</MGGenerationInfo>
</header>
<init>
11 -11 5.000000e+02 5.000000e+02 0 0 247000 247000 -4 1
2.772290e-08 5.351424e-10 2.772290e-08 1
<generator name='MadGraph5_aMC@NLO' version='2.7.3'>please cite 1405.0301 </generator>
</init>
<event>
 9 1 +2.7722900e-08 5.01452400e+02 7.81860800e-03 9.42167800e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 -8.5687059939e+00 +2.4081084389e+01 -4.8954220926e+02 4.9854762783e+02 9.0801108040e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +8.5687059939e+00 -2.4081084389e+01 +4.8954220926e+02 5.0145237217e+02 1.0559159799e+02 0.0000e+00 0.0000e+00
        5 1 4 4 501 0 -1.2761664477e+00 +3.9981873604e+01 +1.3041924828e+02 1.3649706642e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 4 4 0 501 +9.8448724416e+00 -6.4062957994e+01 +3.5912296097e+02 3.6495530575e+02 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +3.0455815943e+01 +6.3364173267e+00 -6.9813173711e+01 7.6430342776e+01 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 -2.6822313581e+01 +1.8473700585e+01 -3.1354676602e+02 3.1523370473e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 -1.2202208356e+01 -7.2903352232e-01 -1.0618226953e+02 1.0688358032e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.01131900e+02 7.81860800e-03 9.42239700e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +1.4516094565e+02 -2.8794425467e+01 -4.6778492081e+02 4.9886807130e+02 9.0254093201e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.4516094565e+02 +2.8794425467e+01 +4.6778492081e+02 5.0113192870e+02 1.0202773519e+02 0.0000e+00 0.0000e+00
        1 1 3 3 501 0 +4.1794586826e+01 +3.1831652057e+01 -1.6787494482e+02 1.7590349253e+02 5.0400000000e-03 0.0000e+00 -1.0000e+00
       -1 1 3 3 0 501 +1.0336635882e+02 -6.0626077524e+01 -2.9990997599e+02 3.2296457877e+02 5.0400000000e-03 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 -4.6340389500e+01 +1.7513793412e+01 +6.4402917338e+01 8.1252148183e+01 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -7.9227143622e+01 +3.4028816626e+01 +2.8568203481e+02 2.9841100124e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -1.9593412523e+01 -2.2748184570e+01 +1.1769996866e+02 1.2146877928e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00005700e+02 7.81860800e-03 9.42492700e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 -8.7819883540e+01 +2.2552074667e+02 +4.2784105272e+02 5.0000570584e+02 9.1574015727e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +8.7819883540e+01 -2.2552074667e+02 -4.2784105272e+02 4.9999429416e+02 9.1511685971e+01 0.0000e+00 0.0000e+00
        1 1 4 4 501 0 +1.0222810354e+02 -2.1881489840e+02 -4.1189641840e+02 4.7748215090e+02 5.0400000000e-03 0.0000e+00 -1.0000e+00
       -1 1 4 4 0 501 -1.4408220003e+01 -6.7058482704e+00 -1.5944634325e+01 2.2512143253e+01 5.0400000000e-03 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 -5.5239294100e+01 +1.6211551513e+02 +2.7152204551e+02 3.2102500248e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 -1.4741546114e+00 +1.9511407459e+00 +6.2730719965e+01 6.2778366573e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 -3.1106434828e+01 +6.1454090799e+01 +9.3588287250e+01 1.1620233679e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00096600e+02 7.81860800e-03 9.42472300e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 -6.7263780694e+00 -1.6290300255e+01 -4.9128635666e+02 5.0009661036e+02 9.1780811727e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +6.7263780694e+00 +1.6290300255e+01 +4.9128635666e+02 4.9990338964e+02 9.0722084901e+01 0.0000e+00 0.0000e+00
        5 1 4 4 501 0 +2.7564281876e+01 +4.6794380823e+01 +2.6179694253e+02 2.6741210300e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 4 4 0 501 -2.0837903807e+01 -3.0504080569e+01 +2.2948941413e+02 2.3249128664e+02 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 -1.3844453256e+01 +9.9711137956e+00 -3.2475511136e+02 3.2520299125e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 -4.7406362024e+00 -3.1820413660e+01 -6.5050341738e+01 7.2571063914e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +1.1858711389e+01 +5.5589996103e+00 -1.0148090356e+02 1.0232255520e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00335700e+02 7.81860800e-03 9.42418500e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 -4.6104202019e+01 -9.6935656470e+00 -4.8945745702e+02 4.9966430302e+02 8.8748244021e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +4.6104202019e+01 +9.6935656470e+00 +4.8945745702e+02 5.0033569698e+02 9.2453473594e+01 0.0000e+00 0.0000e+00
        4 1 4 4 501 0 +4.5410331706e+01 -4.3770365928e+00 +4.8235241443e+02 4.8450667826e+02 1.2700000000e+00 0.0000e+00 -1.0000e+00
       -4 1 4 4 0 501 +6.9387031284e-01 +1.4070602240e+01 +7.1050425969e+00 1.5829018722e+01 1.2700000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 -4.6939926509e+01 -3.4294992471e+00 -4.6209816346e+02 4.6448879642e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +6.8128812980e+00 -7.9305370797e+00 -4.3152040729e+00 1.1310603715e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 -5.9771568077e+00 +1.6664706798e+00 -2.3044089491e+01 2.3864902881e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00063600e+02 7.81860800e-03 9.42479700e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 +5.7960736542e+01 +9.8610747947e+01 +4.7784974042e+02 4.9993638732e+02 9.2263158148e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -5.7960736542e+01 -9.8610747947e+01 -4.7784974042e+02 5.0006361268e+02 9.2950071095e+01 0.0000e+00 0.0000e+00
        4 1 3 3 501 0 +5.5928333709e+01 +8.5248265274e+00 +1.7032914094e+02 1.7948337063e+02 1.2700000000e+00 0.0000e+00 -1.0000e+00
       -4 1 3 3 0 501 +2.0324028333e+00 +9.0085921420e+01 +3.0752059948e+02 3.2045301669e+02 1.2700000000e+00 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 -1.5900092058e+01 -2.9050919562e+00 -2.7950239978e+01 3.2287452138e+01 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -3.1679517596e+01 -1.0082813330e+02 -3.1246396993e+02 3.2985396284e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -1.0381126888e+01 +5.1224773119e+00 -1.3743553052e+02 1.3792219770e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.01542400e+02 7.81860800e-03 9.42147700e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +1.4551366235e+02 +4.7106678566e+02 -8.1178482839e+00 5.0154240168e+02 9.1655543979e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.4551366235e+02 -4.7106678566e+02 +8.1178482839e+00 4.9845759832e+02 7.2910461350e+01 0.0000e+00 0.0000e+00
        1 1 3 3 501 0 +8.0329178701e+00 +8.8453765086e+00 -1.5102873424e+01 1.9257862450e+01 5.0400000000e-03 0.0000e+00 -1.0000e+00
       -1 1 3 3 0 501 +1.3748074448e+02 +4.6222140915e+02 +6.9850251404e+00 4.8228453923e+02 5.0400000000e-03 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 -6.1388361899e+01 -1.2274889280e+02 +1.0550328339e+01 1.3764861878e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -4.5707425269e+01 -1.8929045801e+02 +1.9665558308e+01 1.9572118026e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -3.8417875178e+01 -1.5902743486e+02 -2.2098038364e+01 1.6508779928e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00689100e+02 7.81860800e-03 9.42339100e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 +1.1234074203e+02 +3.8688622397e+02 +2.8312327030e+02 5.0068913765e+02 9.0716226577e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.1234074203e+02 -3.8688622397e+02 -2.8312327030e+02 4.9931086235e+02 8.2771725022e+01 0.0000e+00 0.0000e+00
        2 1 3 3 501 0 -6.1620975716e-01 -2.4114485849e+00 -3.4583875023e+00 4.2608977178e+00 2.5500000000e-03 0.0000e+00 -1.0000e+00
       -2 1 3 3 0 501 +1.1295695178e+02 +3.8929767255e+02 +2.8658165780e+02 4.9642823993e+02 2.5500000000e-03 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 -7.4815827833e+00 -1.1752953452e+02 -1.1638696401e+02 1.6557506491e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -3.3756695480e+01 -8.3751028413e+01 -5.4757278338e+01 1.0560354531e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -7.1102463762e+01 -1.8560566103e+02 -1.1197902795e+02 2.2813225213e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00852300e+02 7.81860800e-03 9.42302500e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +9.0973757160e+00 -9.0790032205e+01 -4.8384296810e+02 5.0085229412e+02 9.1779140036e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -9.0973757160e+00 +9.0790032205e+01 +4.8384296810e+02 4.9914770588e+02 8.1968422649e+01 0.0000e+00 0.0000e+00
        3 1 4 4 501 0 -2.7719582600e+01 +7.6053979976e+01 +4.5323835473e+02 4.6041025133e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 1 4 4 0 501 +1.8622206884e+01 +1.4736052228e+01 +3.0604613376e+01 3.8737454555e+01 1.0100000000e-01 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +1.1021365264e+01 -8.9442743283e+01 -4.8117941247e+02 4.8954582315e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 -5.1819703608e+00 -3.7341167060e+00 +1.4218510051e-01 6.3887918465e+00 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +3.2579808124e+00 +2.3868277844e+00 -2.8057407398e+00 4.9176791218e+00 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00053700e+02 7.81860800e-03 9.42481900e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 -2.4552899045e+01 +1.2079071880e+02 -4.7605371625e+02 4.9994627502e+02 9.0143743851e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +2.4552899045e+01 -1.2079071880e+02 +4.7605371625e+02 5.0005372498e+02 9.0737778869e+01 0.0000e+00 0.0000e+00
        5 1 3 3 501 0 -1.2602969622e+01 +1.2663346689e+02 -4.6471266928e+02 4.8184522906e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 3 3 0 501 -1.1949929424e+01 -5.8427480913e+00 -1.1341046963e+01 1.8101045956e+01 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 -1.4534605911e+01 -8.4531548739e+01 +3.5521621724e+02 3.6542497131e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 +3.4988925687e+01 -1.7380194742e+01 +9.1621747704e+01 9.9603417324e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 +4.0985792686e+00 -1.8878975316e+01 +2.9215751301e+01 3.5025336345e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.01609800e+02 7.81860800e-03 9.42132600e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 -2.8224192250e+02 +7.9055193001e+01 +3.9278616542e+02 4.9839019549e+02 9.0562623936e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +2.8224192250e+02 -7.9055193001e+01 -3.9278616542e+02 5.0160980451e+02 1.0687000459e+02 0.0000e+00 0.0000e+00
        5 1 4 4 501 0 +7.3124730290e+01 +6.2630020413e+00 -4.2918427826e+01 8.5150060609e+01 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 4 4 0 501 +2.0911719221e+02 -8.5318195042e+01 -3.4986773759e+02 4.1645974390e+02 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 -1.2084849871e-01 +4.6933931285e+00 +2.3272010658e+01 2.3741107546e+01 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 -9.0125678109e+01 +4.9624928520e+01 +1.5101174363e+02 1.8272881026e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 -1.9199539589e+02 +2.4736871352e+01 +2.1850241114e+02 2.9192027768e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00999800e+02 7.81860800e-03 9.42269300e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 +2.3093144421e+02 -2.6741600640e+02 +3.4029803585e+02 5.0099983917e+02 1.0177147557e+02 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -2.3093144421e+02 +2.6741600640e+02 -3.4029803585e+02 4.9900016083e+02 9.1420757504e+01 0.0000e+00 0.0000e+00
        5 1 4 4 501 0 -1.9866301595e+02 +2.1748757868e+02 -3.2242214931e+02 4.3674474599e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 4 4 0 501 -3.2268428263e+01 +4.9928427719e+01 -1.7875886536e+01 6.2255414837e+01 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +1.9431409165e+02 -1.8778804690e+02 +2.2768475516e+02 3.5335913129e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +1.6129496853e+01 -6.2329849838e+01 +9.7563146390e+01 1.1689199452e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +2.0487855709e+01 -1.7298109660e+01 +1.5050134301e+01 3.0748713352e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00420400e+02 7.81860800e-03 9.42399500e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 -7.0752902142e+00 +2.0322340825e+01 -4.9114154653e+02 4.9957958273e+02 8.8863285394e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +7.0752902142e+00 -2.0322340825e+01 +4.9114154653e+02 5.0042041727e+02 9.3474691892e+01 0.0000e+00 0.0000e+00
        5 1 3 3 501 0 +1.6315294630e+01 +1.2124181472e+01 -4.5514610696e+02 4.5562402625e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 3 3 0 501 -2.3390584844e+01 +8.1981593529e+00 -3.5995439565e+01 4.3955556481e+01 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 +2.3563808635e+01 +3.1598674486e+01 +2.5834843603e+02 2.6133820017e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -7.9792002516e+00 -2.5191465848e+01 +1.1270681984e+02 1.1576314105e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -8.5093181693e+00 -2.6729549463e+01 +1.2008629066e+02 1.2331907604e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00349400e+02 7.81860800e-03 9.42415400e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 -3.8294145079e+01 +6.4824552854e+00 -4.9031914735e+02 5.0034941073e+02 9.1805243780e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +3.8294145079e+01 -6.4824552854e+00 +4.9031914735e+02 4.9965058927e+02 8.7916900079e+01 0.0000e+00 0.0000e+00
        2 1 3 3 501 0 +1.6354233049e+01 -3.0259107366e+00 -3.3517298617e+02 3.3558538070e+02 2.5500000000e-03 0.0000e+00 1.0000e+00
       -2 1 3 3 0 501 -5.4648378128e+01 +9.5083660220e+00 -1.5514616118e+02 1.6476403003e+02 2.5500000000e-03 0.0000e+00 -1.0000e+00
       13 1 4 4 0 0 +5.3579433878e+01 +1.4287430380e+01 +2.1294835784e+02 2.2004976862e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -1.8807594931e+01 -6.2255266708e+00 +1.7997043677e+02 1.8105756245e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 +3.5223061322e+00 -1.4544358995e+01 +9.7400352742e+01 9.8543258188e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00074900e+02 7.81860800e-03 9.42477200e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +1.0975694230e+02 -1.7608762782e+02 -4.4555472838e+02 4.9992507863e+02 9.1392719658e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.0975694230e+02 +1.7608762782e+02 +4.4555472838e+02 5.0007492137e+02 9.2208849594e+01 0.0000e+00 0.0000e+00
        1 1 4 4 501 0 -9.2699189476e+01 +6.2283367475e+01 +2.0651655013e+02 2.3477956278e+02 5.0400000000e-03 0.0000e+00 1.0000e+00
       -1 1 4 4 0 501 -1.7057752825e+01 +1.1380426034e+02 +2.3903817825e+02 2.6529535859e+02 5.0400000000e-03 0.0000e+00 -1.0000e+00
       13 1 3 3 0 0 -1.7625395352e+00 -1.1116953158e+01 -4.6815467408e+00 1.2191031009e+01 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +2.2178840192e+00 +5.3558295265e+00 -1.8334557234e+01 1.9229142172e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +1.0930159782e+02 -1.7032650419e+02 -4.2253862441e+02 4.6850490545e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.01104300e+02 7.81860800e-03 9.42245900e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 +1.1943581696e+02 +3.3076130676e+02 +3.4166546210e+02 4.9889571211e+02 9.2161201972e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.1943581696e+02 -3.3076130676e+02 -3.4166546210e+02 5.0110428789e+02 1.0345174201e+02 0.0000e+00 0.0000e+00
        3 1 3 3 501 0 +5.7779716583e+01 +6.2596791409e+01 +1.1046661490e+02 1.3949816182e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 1 3 3 0 501 +6.1656100381e+01 +2.6816451536e+02 +2.3119884720e+02 3.5939755029e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
       13 1 4 4 0 0 +6.3672683948e+00 -1.1647809624e+02 -1.2768095511e+02 1.7294544363e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 4 4 0 0 -1.0000692568e+02 -1.7582250028e+02 -1.7925604268e+02 2.7027331653e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 4 4 0 0 -2.5796159681e+01 -3.8460710245e+01 -3.4728464305e+01 5.7885527732e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00021900e+02 7.81860800e-03 9.42489100e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 -1.8294216674e+01 +4.8060819032e+02 -1.0408161523e+02 4.9997806975e+02 8.8465680515e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 +1.8294216674e+01 -4.8060819032e+02 +1.0408161523e+02 5.0002193025e+02 8.8713229779e+01 0.0000e+00 0.0000e+00
        3 1 4 4 501 0 +2.6060322989e+01 -4.5664873908e+02 +1.1641975745e+02 4.7197540345e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 1 4 4 0 501 -7.7661063151e+00 -2.3959451241e+01 -1.2338142223e+01 2.8046526800e+01 1.0100000000e-01 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 -3.7088229508e+01 +1.8490654064e+02 -4.4534513602e+01 1.9377641653e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +1.8389203810e+01 +3.7393179231e+01 -3.0161086237e+00 4.1779296084e+01 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +4.0480902328e-01 +2.5830847045e+02 -5.6530993005e+01 2.6442235713e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00674100e+02 7.81860800e-03 9.42342500e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +1.6925322991e+01 -1.3548633623e+02 +4.7189822486e+02 4.9932591867e+02 8.9417137950e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.6925322991e+01 +1.3548633623e+02 -4.7189822486e+02 5.0067408133e+02 9.6662232677e+01 0.0000e+00 0.0000e+00
        3 1 4 4 501 0 -7.3363913477e+00 -1.2910052829e+01 -1.9814947504e+01 2.4761551802e+01 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 1 4 4 0 501 -9.5889316431e+00 +1.4839638906e+02 -4.5208327736e+02 4.7591252953e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +1.0102847513e+01 -1.9268583633e+01 +1.4961086650e+02 1.5118455074e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +2.6754462386e+01 -4.8381520609e+01 +1.9325292745e+02 2.0100563864e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 -1.9931986909e+01 -6.7836231992e+01 +1.2903443090e+02 1.4713572929e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00115100e+02 7.81860800e-03 9.42468100e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
       23 2 1 2 0 0 +1.7879825361e+01 -1.7169818129e+01 -4.9128571986e+02 4.9988494729e+02 8.8931497071e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.7879825361e+01 +1.7169818129e+01 +4.9128571986e+02 5.0011505271e+02 9.0215944260e+01 0.0000e+00 0.0000e+00
        3 1 4 4 501 0 -5.3086113521e+00 -3.5329105247e+01 +1.3004536861e+02 1.3486339433e+02 1.0100000000e-01 0.0000e+00 -1.0000e+00
       -3 1 4 4 0 501 -1.2571214009e+01 +5.2498923376e+01 +3.6124035126e+02 3.6525165839e+02 1.0100000000e-01 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +1.3778881228e+01 -3.7435271789e+01 -1.6589317013e+02 1.7062183976e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +2.9832692881e+00 +3.2918111932e+01 -2.8525768816e+02 2.8716624217e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +1.1176748447e+00 -1.2652658271e+01 -4.0134861575e+01 4.2096865359e+01 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
<event>
 9 1 +2.7722900e-08 5.00121800e+02 7.81860800e-03 9.42466600e-02
       11 -1 0 0 0 0 +0.0000000000e+00 +0.0000000000e+00 +5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
      -11 -1 0 0 0 0 -0.0000000000e+00 -0.0000000000e+00 -5.0000000000e+02 5.0000000000e+02 5.1100000000e-04 0.0000e+00 1.0000e+00
       23 2 1 2 0 0 +1.1709481882e+02 -2.6171522587e+01 -4.7681985543e+02 4.9987815074e+02 9.0137926904e+01 0.0000e+00 0.0000e+00
       23 2 1 2 0 0 -1.1709481882e+02 +2.6171522587e+01 +4.7681985543e+02 5.0012184926e+02 9.1479748547e+01 0.0000e+00 0.0000e+00
        5 1 4 4 501 0 -6.4781625088e+01 +2.9800198390e+01 +4.0394906520e+02 4.1022146220e+02 4.7000000000e+00 0.0000e+00 -1.0000e+00
       -5 1 4 4 0 501 -5.2313193728e+01 -3.6286758026e+00 +7.2870790226e+01 8.9900387064e+01 4.7000000000e+00 0.0000e+00 1.0000e+00
       13 1 3 3 0 0 +4.3967164476e+01 +9.2351886676e+00 -2.2921267916e+02 2.3357410754e+02 1.0566000000e-01 0.0000e+00 1.0000e+00
      -11 1 3 3 0 0 +3.8681927428e+01 -4.5697940987e+01 -1.2612541731e+02 1.3961452006e+02 5.1100000000e-04 0.0000e+00 -1.0000e+00
  9000005 1 3 3 0 0 +3.4445726913e+01 +1.0291229732e+01 -1.2148175895e+02 1.2668952313e+02 1.0000000000e-03 0.0000e+00 0.0000e+00
</event>
</LesHouchesEvents>

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2020-08-22
Last query:
2020-08-22
Last reply:
2020-08-21

This cross-section should be infinite.

In some cases the phase-space integrator does not find the singularity and returns wrongly a finite cross-section and the events that you requuest) and in some cases it hits close enough of the singularity and the cross-section blows up (as it should be) which prevents the event generation
(In order to generate events, you need the matrix-element to be finite over the full phase-space.

So please do not use any of those events even if you have the correct number of events generated.

The issue is that you can not ask the muon to decay via an onshell decay
(as you ask with this syntax: z > mu- mu+, mu+ > e+ ax)
and at the same time tell madgraph that muon are stable:
(as you ask with this line DECAY 13 0.000000e+00 # mu- : 0.0)

Note that the way your model is written the width of the muon is fixed to zero from your model.
you need to change the UFO model if you want to have MG5aMC to use a value different of zero.

This being said, the width of the muon is very small and including such width can lead to numerical accuracy problem. We have a security and a special handling fo such case but you need to be aware of their impact on distribution: https://answers.launchpad.net/mg5amcnlo/+faq/3053 <https://answers.launchpad.net/mg5amcnlo/+faq/3053>

Cheers,

Olivier

HH (hh-madgrapher) said : #2

Dear Olivier,

Thanks a lot for your prompt answer. I modified my UFO model so that the decay width of the muon is automatically computed based on the given couplings. This gives a width of O(10^-11). My event generation problem was solved. Thanks again. As a complementary question: I need to perform a scan over a coupling range to find xsection as a function of that coupling. The decay width in this range is of the same order as O(10^-11) or even smaller. Will I encounter any problem in obtaining true xsections as the decay width is so small? In the link you sent it is said that the total xsection is not altered by the small width treatment procedure. Actually I didn't know what you mean by the total xsection. I need the xsection of this process: e- e+ > z z , z > j j , (z > mu- mu+ , mu+ > e+ ax) with exactly the specified final state and decays.

total cross-section is the sum of the cross-section computed by MG5aMC (since you can have more than one process at the time).

In principle your scan should work but this is your job to perform validation and be sure that this is indeed the case.

Cheers,

Olivier

> On 21 Aug 2020, at 09:20, HH <email address hidden> wrote:
>
> Question #692457 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692457
>
> Status: Answered => Open
>
> HH is still having a problem:
> Dear Olivier,
>
> Thanks a lot for your prompt answer. I modified my UFO model so that the
> decay width of the muon is automatically computed based on the given
> couplings. This gives a width of O(10^-11). My event generation problem
> was solved. Thanks again. As a complementary question: I need to perform
> a scan over a coupling range to find xsection as a function of that
> coupling. The decay width in this range is of the same order as
> O(10^-11) or even smaller. Will I encounter any problem in obtaining
> true xsections as the decay width is so small? In the link you sent it
> is said that the total xsection is not altered by the small width
> treatment procedure. Actually I didn't know what you mean by the total
> xsection. I need the xsection of this process: e- e+ > z z , z > j j ,
> (z > mu- mu+ , mu+ > e+ ax) with exactly the specified final state and
> decays.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

HH (hh-madgrapher) said : #4

Dear Olivier,

Thanks a lot for you help.

HH (hh-madgrapher) said : #5

Thanks Olivier Mattelaer, that solved my question.