examples/filter_events_example.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "a69095ea-7aff-4efa-bcde-928be2559aba",
   "metadata": {},
   "outputs": [],
   "source": [
    "import gzip\n",
    "\n",
    "import hist\n",
    "from skhep_testdata import data_path\n",
    "\n",
    "import pylhe"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "1e6bce3e-b076-4862-8a30-7dd9a08c3284",
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot(data):\n",
    "    lheevents = pylhe.read_lhe_file(data).events\n",
    "    events = pylhe.to_awkward(lheevents)\n",
    "    mass_hist = hist.Hist.new.Reg(30, 50, 150).Weight()\n",
    "    mass_hist.fill(\n",
    "        (events.particles.vector[:, -1] + events.particles.vector[:, -2]).mass,\n",
    "        weight=events.eventinfo.weight,\n",
    "    )\n",
    "    artists = mass_hist.plot()\n",
    "    ax = artists[0].stairs.axes\n",
    "    ax.set_yscale(\"log\")\n",
    "    ax.set_xlabel(\"Mass [GeV]\")\n",
    "    ax.set_ylabel(\"Count\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "7649faa3-9ecc-4761-96ce-857c6ee92868",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "b'<LesHouchesEvents version=\"3.0\">\\n'\n",
      "b'<header>\\n'\n",
      "b'<!--\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                        MadGraph5_aMC@NLO                           *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                           Going Beyond                             *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                   http://madgraph.hep.uiuc.edu                     *\\n'\n",
      "b'#                   http://madgraph.phys.ucl.ac.be                   *\\n'\n",
      "b'#                   http://amcatnlo.cern.ch                          *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                     The MadGraph5_aMC@NLO team                     *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#....................................................................*\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'# This file contains all the information necessary to reproduce      *\\n'\n",
      "b'# the events generated:                                              *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'# 1. software version                                                *\\n'\n",
      "b'# 2. proc_card          : code generation info including model       *\\n'\n",
      "b'# 3. param_card         : model primary parameters in the LH format  *\\n'\n",
      "b'# 4. run_card           : running parameters (collider and cuts)     *\\n'\n",
      "b'# 5. pythia_card        : present only if pythia has been run        *\\n'\n",
      "b'# 6. pgs_card           : present only if pgs has been run           *\\n'\n",
      "b'# 7. delphes_cards      : present only if delphes has been run       *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'-->\\n'\n",
      "b'<MGVersion>\\n'\n",
      "b'3.1.0\\n'\n",
      "b'</MGVersion>\\n'\n",
      "b'<MG5ProcCard>\\n'\n",
      "b'<![CDATA[\\n'\n",
      "b'#************************************************************\\n'\n",
      "b'#*                     MadGraph5_aMC@NLO                    *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*                *                       *                 *\\n'\n",
      "b'#*                  *        * *        *                   *\\n'\n",
      "b'#*                    * * * * 5 * * * *                     *\\n'\n",
      "b'#*                  *        * *        *                   *\\n'\n",
      "b'#*                *                       *                 *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*         VERSION 3.1.0                 2021-03-30         *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*    The MadGraph5_aMC@NLO Development Team - Find us at   *\\n'\n",
      "b'#*    https://server06.fynu.ucl.ac.be/projects/madgraph     *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#************************************************************\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*               Command File for MadGraph5_aMC@NLO         *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#*     run as ./bin/mg5_aMC  filename                       *\\n'\n",
      "b'#*                                                          *\\n'\n",
      "b'#************************************************************\\n'\n",
      "b'set group_subprocesses Auto\\n'\n",
      "b'set ignore_six_quark_processes False\\n'\n",
      "b'set low_mem_multicore_nlo_generation False\\n'\n",
      "b'set complex_mass_scheme False\\n'\n",
      "b'set include_lepton_initiated_processes False\\n'\n",
      "b'set gauge unitary\\n'\n",
      "b'set loop_optimized_output True\\n'\n",
      "b'set loop_color_flows False\\n'\n",
      "b'set max_npoint_for_channel 0\\n'\n",
      "b'set default_unset_couplings 99\\n'\n",
      "b'set max_t_for_channel 99\\n'\n",
      "b'set zerowidth_tchannel True\\n'\n",
      "b'set nlo_mixed_expansion True\\n'\n",
      "b'generate p p > l+ l-\\n'\n",
      "b'define p = g u c d s u~ c~ d~ s~\\n'\n",
      "b'define j = g u c d s u~ c~ d~ s~\\n'\n",
      "b'define l+ = e+ mu+\\n'\n",
      "b'define l- = e- mu-\\n'\n",
      "b'define vl = ve vm vt\\n'\n",
      "b'define vl~ = ve~ vm~ vt~\\n'\n",
      "b'output drell-yan_output\\n'\n",
      "b']]>\\n'\n",
      "b'</MG5ProcCard>\\n'\n",
      "b'<MGProcCard>\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#                        MadGraph/MadEvent                           *\\n'\n",
      "b'#                   http://madgraph.hep.uiuc.edu                     *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                          proc_card.dat                             *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#            This Files is generated by MADGRAPH 5                   *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'# WARNING: This Files is generated for MADEVENT (compatibility issue)*\\n'\n",
      "b'#          This files is NOT a valid MG4 proc_card.dat               *\\n'\n",
      "b'#          Running this in MG4 will NEVER reproduce the result of MG5*\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Process(es) requested : mg2 input                                  *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Begin PROCESS # This is TAG. Do not modify this line\\n'\n",
      "b'p p > l+ l-           #Process\\n'\n",
      "b'# Be carefull the coupling are here in MG5 convention\\n'\n",
      "b'          \\n'\n",
      "b'end_coup           # End the couplings input\\n'\n",
      "b'\\n'\n",
      "b'done               # this tells MG there are no more procs\\n'\n",
      "b'# End PROCESS  # This is TAG. Do not modify this line\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Model information                                                  *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Begin MODEL  # This is TAG. Do not modify this line\\n'\n",
      "b'sm\\n'\n",
      "b'# End   MODEL  # This is TAG. Do not modify this line\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Start multiparticle definitions                                    *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Begin MULTIPARTICLES # This is TAG. Do not modify this line\\n'\n",
      "b'\\n'\n",
      "b'# End  MULTIPARTICLES # This is TAG. Do not modify this line\\n'\n",
      "b'</MGProcCard>\\n'\n",
      "b'<MGRunCard>\\n'\n",
      "b'<![CDATA[\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#                       MadGraph5_aMC@NLO                            *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#                     run_card.dat MadEvent                          *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#  This file is used to set the parameters of the run.               *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#  Some notation/conventions:                                        *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b\"#   Lines starting with a '# ' are info or comments                  *\\n\"\n",
      "b'#                                                                    *\\n'\n",
      "b'#   mind the format:   value    = variable     ! comment             *\\n'\n",
      "b'#                                                                    *\\n'\n",
      "b'#   To display more options, you can type the command:               *\\n'\n",
      "b'#      update full_run_card                                          *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#                                                                    \\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Tag name for the run (one word)                                    *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'  tag_1     = run_tag ! name of the run \\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Number of events and rnd seed                                      *\\n'\n",
      "b'# Warning: Do not generate more than 1M events in a single run       *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'  10000 = nevents ! Number of unweighted events requested \\n'\n",
      "b' 21 = iseed   ! rnd seed (0=assigned automatically=default))\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Collider type and energy                                           *\\n'\n",
      "b'# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton,      *\\n'\n",
      "b'#                3=photon from electron, 4=photon from muon          *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'     1        = lpp1    ! beam 1 type \\n'\n",
      "b'     1        = lpp2    ! beam 2 type\\n'\n",
      "b'     6500.0     = ebeam1  ! beam 1 total energy in GeV\\n'\n",
      "b'     6500.0     = ebeam2  ! beam 2 total energy in GeV\\n'\n",
      "b'# To see polarised beam options: type \"update beam_pol\"\\n'\n",
      "b'\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# PDF CHOICE: this automatically fixes also alpha_s and its evol.    *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'     nn23lo1    = pdlabel     ! PDF set                                     \\n'\n",
      "b'     230000    = lhaid     ! if pdlabel=lhapdf, this is the lhapdf number\\n'\n",
      "b'# To see heavy ion options: type \"update ion_pdf\"\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Renormalization and factorization scales                           *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' False = fixed_ren_scale  ! if .true. use fixed ren scale\\n'\n",
      "b' False        = fixed_fac_scale  ! if .true. use fixed fac scale\\n'\n",
      "b' 91.188  = scale            ! fixed ren scale\\n'\n",
      "b' 91.188  = dsqrt_q2fact1    ! fixed fact scale for pdf1\\n'\n",
      "b' 91.188  = dsqrt_q2fact2    ! fixed fact scale for pdf2\\n'\n",
      "b' -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices\\n'\n",
      "b' 1.0  = scalefact        ! scale factor for event-by-event scales\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Type and output format\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'  False     = gridpack  !True = setting up the grid pack\\n'\n",
      "b'  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)\\n'\n",
      "b'  average =  event_norm       ! average/sum. Normalization of the weight in the LHEF\\n'\n",
      "b'# To see MLM/CKKW  merging options: type \"update MLM\" or \"update CKKW\"\\n'\n",
      "b'\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Phase-Space Optimization strategy (basic options)\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'   0  = nhel          ! using helicities importance sampling or not.\\n'\n",
      "b'                             ! 0: sum over helicity, 1: importance sampling\\n'\n",
      "b'   2  = sde_strategy  ! default integration strategy (hep-ph/2021.xxxxx)\\n'\n",
      "b'                             ! 1 is old strategy (using amp square)\\n'\n",
      "b'\\t\\t\\t     ! 2 is new strategy (using only the denominator)\\n'\n",
      "b'# To see advanced option for Phase-Space optimization: type \"update psoptim\"\\t\\t\\t     \\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Generation bias, check the wiki page below for more information:   *\\n'\n",
      "b\"#  'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *\\n\"\n",
      "b'#*********************************************************************\\n'\n",
      "b' None = bias_module  ! Bias type of bias, [None, ptj_bias, -custom_folder-]\\n'\n",
      "b' {} = bias_parameters ! Specifies the parameters of the module.\\n'\n",
      "b'#\\n'\n",
      "b'#*******************************                                                 \\n'\n",
      "b'# Parton level cuts definition *\\n'\n",
      "b'#*******************************                                     \\n'\n",
      "b'#                                                                    \\n'\n",
      "b'#\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for \"$\" and decay  \\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'  15.0  = bwcutoff      ! (M+/-bwcutoff*Gamma)\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Standard Cuts                                                      *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b\"# Minimum and maximum pt's (for max, -1 means no cut)                *\\n\"\n",
      "b'#*********************************************************************\\n'\n",
      "b' 10.0  = ptl       ! minimum pt for the charged leptons \\n'\n",
      "b' -1.0  = ptlmax    ! maximum pt for the charged leptons\\n'\n",
      "b' {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle\\n'\n",
      "b' {}\\t= pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50}) \\n'\n",
      "b'#\\n'\n",
      "b\"# For display option for energy cut in the partonic center of mass frame type 'update ecut'\\n\"\n",
      "b'#\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Maximum and minimum absolute rapidity (for max, -1 means no cut)   *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' 2.5  = etal    ! max rap for the charged leptons \\n'\n",
      "b' 0.0  = etalmin ! main rap for the charged leptons\\n'\n",
      "b' {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle\\n'\n",
      "b' {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Minimum and maximum DeltaR distance                                *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' 0.4 = drll    ! min distance between leptons \\n'\n",
      "b' -1.0  = drllmax ! max distance between leptons\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Minimum and maximum invariant mass for pairs                       *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' 0.0   = mmll    ! min invariant mass of l+l- (same flavour) lepton pair\\n'\n",
      "b' -1.0  = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair\\n'\n",
      "b' {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})\\n'\n",
      "b\" {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only \\n\"\n",
      "b'                       ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.  \\n'\n",
      "b' #*********************************************************************\\n'\n",
      "b' # Minimum and maximum invariant mass for all letpons                 *\\n'\n",
      "b' #*********************************************************************\\n'\n",
      "b' 0.0   = mmnl    ! min invariant mass for all letpons (l+- and vl) \\n'\n",
      "b' -1.0  = mmnlmax ! max invariant mass for all letpons (l+- and vl) \\n'\n",
      "b' #*********************************************************************\\n'\n",
      "b' # Minimum and maximum pt for 4-momenta sum of leptons / neutrino     *\\n'\n",
      "b' #  for pair of lepton includes only same flavor, opposite charge\\n'\n",
      "b' #*********************************************************************\\n'\n",
      "b' 0.0   = ptllmin  ! Minimum pt for 4-momenta sum of leptons(l and vl)\\n'\n",
      "b' -1.0  = ptllmax  ! Maximum pt for 4-momenta sum of leptons(l and vl)\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Inclusive cuts                                                     *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' 0.0  = xptl ! minimum pt for at least one charged lepton \\n'\n",
      "b' #*********************************************************************\\n'\n",
      "b\" # Control the pt's of leptons sorted by pt                           *\\n\"\n",
      "b' #*********************************************************************\\n'\n",
      "b' 0.0   = ptl1min ! minimum pt for the leading lepton in pt\\n'\n",
      "b' 0.0   = ptl2min ! minimum pt for the second lepton in pt\\n'\n",
      "b' -1.0  = ptl1max ! maximum pt for the leading lepton in pt \\n'\n",
      "b' -1.0  = ptl2max ! maximum pt for the second lepton in pt\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# maximal pdg code for quark to be considered as a light jet         *\\n'\n",
      "b'# (otherwise b cuts are applied)                                     *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b' 4 = maxjetflavor    ! Maximum jet pdg code\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'#\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'# Store info for systematics studies                                 *\\n'\n",
      "b'# WARNING: Do not use for interference type of computation           *\\n'\n",
      "b'#*********************************************************************\\n'\n",
      "b'   True  = use_syst      ! Enable systematics studies\\n'\n",
      "b'#\\n'\n",
      "b'systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]\\n'\n",
      "b\"['--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\\n\"\n",
      "b\"# Syscalc is deprecated but to see the associate options type'update syscalc'\\n\"\n",
      "b']]>\\n'\n",
      "b'</MGRunCard>\\n'\n",
      "b'<slha>\\n'\n",
      "b'######################################################################\\n'\n",
      "b'## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL   ####\\n'\n",
      "b'######################################################################\\n'\n",
      "b'##                                                                  ##\\n'\n",
      "b'##  Width set on Auto will be computed following the information    ##\\n'\n",
      "b'##        present in the decay.py files of the model.               ##\\n'\n",
      "b'##        See  arXiv:1402.1178 for more details.                    ##\\n'\n",
      "b'##                                                                  ##\\n'\n",
      "b'######################################################################\\n'\n",
      "b'\\n'\n",
      "b'###################################\\n'\n",
      "b'## INFORMATION FOR MASS\\n'\n",
      "b'###################################\\n'\n",
      "b'Block mass \\n'\n",
      "b'    5 4.700000e+00 # MB \\n'\n",
      "b'    6 1.730000e+02 # MT \\n'\n",
      "b'   15 1.777000e+00 # MTA \\n'\n",
      "b'   23 9.118800e+01 # MZ \\n'\n",
      "b'   25 1.250000e+02 # MH \\n'\n",
      "b'## Dependent parameters, given by model restrictions.\\n'\n",
      "b'## Those values should be edited following the \\n'\n",
      "b'## analytical expression. MG5 ignores those values \\n'\n",
      "b'## but they are important for interfacing the output of MG5\\n'\n",
      "b'## to external program such as Pythia.\\n'\n",
      "b'  1 0.000000e+00 # d : 0.0 \\n'\n",
      "b'  2 0.000000e+00 # u : 0.0 \\n'\n",
      "b'  3 0.000000e+00 # s : 0.0 \\n'\n",
      "b'  4 0.000000e+00 # c : 0.0 \\n'\n",
      "b'  11 0.000000e+00 # e- : 0.0 \\n'\n",
      "b'  12 0.000000e+00 # ve : 0.0 \\n'\n",
      "b'  13 0.000000e+00 # mu- : 0.0 \\n'\n",
      "b'  14 0.000000e+00 # vm : 0.0 \\n'\n",
      "b'  16 0.000000e+00 # vt : 0.0 \\n'\n",
      "b'  21 0.000000e+00 # g : 0.0 \\n'\n",
      "b'  22 0.000000e+00 # a : 0.0 \\n'\n",
      "b'  24 8.041900e+01 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2))) \\n'\n",
      "b'\\n'\n",
      "b'###################################\\n'\n",
      "b'## INFORMATION FOR SMINPUTS\\n'\n",
      "b'###################################\\n'\n",
      "b'Block sminputs \\n'\n",
      "b'    1 1.325070e+02 # aEWM1 \\n'\n",
      "b'    2 1.166390e-05 # Gf \\n'\n",
      "b'    3 1.180000e-01 # aS (Note that Parameter not used if you use a PDF set) \\n'\n",
      "b'\\n'\n",
      "b'###################################\\n'\n",
      "b'## INFORMATION FOR YUKAWA\\n'\n",
      "b'###################################\\n'\n",
      "b'Block yukawa \\n'\n",
      "b'    5 4.700000e+00 # ymb \\n'\n",
      "b'    6 1.730000e+02 # ymt \\n'\n",
      "b'   15 1.777000e+00 # ymtau \\n'\n",
      "b'\\n'\n",
      "b'###################################\\n'\n",
      "b'## INFORMATION FOR DECAY\\n'\n",
      "b'###################################\\n'\n",
      "b'DECAY   6 1.491500e+00 # WT \\n'\n",
      "b'DECAY  23 2.441404e+00 # WZ \\n'\n",
      "b'DECAY  24 2.047600e+00 # WW \\n'\n",
      "b'DECAY  25 6.382339e-03 # WH \\n'\n",
      "b'## Dependent parameters, given by model restrictions.\\n'\n",
      "b'## Those values should be edited following the \\n'\n",
      "b'## analytical expression. MG5 ignores those values \\n'\n",
      "b'## but they are important for interfacing the output of MG5\\n'\n",
      "b'## to external program such as Pythia.\\n'\n",
      "b'DECAY  1 0.000000e+00 # d : 0.0 \\n'\n",
      "b'DECAY  2 0.000000e+00 # u : 0.0 \\n'\n",
      "b'DECAY  3 0.000000e+00 # s : 0.0 \\n'\n",
      "b'DECAY  4 0.000000e+00 # c : 0.0 \\n'\n",
      "b'DECAY  5 0.000000e+00 # b : 0.0 \\n'\n",
      "b'DECAY  11 0.000000e+00 # e- : 0.0 \\n'\n",
      "b'DECAY  12 0.000000e+00 # ve : 0.0 \\n'\n",
      "b'DECAY  13 0.000000e+00 # mu- : 0.0 \\n'\n",
      "b'DECAY  14 0.000000e+00 # vm : 0.0 \\n'\n",
      "b'DECAY  15 0.000000e+00 # ta- : 0.0 \\n'\n",
      "b'DECAY  16 0.000000e+00 # vt : 0.0 \\n'\n",
      "b'DECAY  21 0.000000e+00 # g : 0.0 \\n'\n",
      "b'DECAY  22 0.000000e+00 # a : 0.0\\n'\n",
      "b'</slha>\\n'\n",
      "b'<initrwgt>\\n'\n",
      "b'<weightgroup name=\"Central scale variation\" combine=\"envelope\">\\n'\n",
      "b'<weight id=\"1\" MUR=\"0.5\"  MUF=\"0.5\"  PDF=\"247000\" > MUR=0.5 MUF=0.5  </weight>\\n'\n",
      "b'<weight id=\"2\" MUR=\"0.5\"  MUF=\"0.5\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=0.5 MUF=0.5 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"3\" MUR=\"0.5\"  MUF=\"0.5\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=0.5 MUF=0.5 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"4\" MUR=\"0.5\"  MUF=\"0.5\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=0.5 MUF=0.5 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"5\" MUR=\"0.5\"  MUF=\"0.5\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=0.5 MUF=0.5 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"6\" MUR=\"0.5\"  MUF=\"1.0\"  PDF=\"247000\" > MUR=0.5  </weight>\\n'\n",
      "b'<weight id=\"7\" MUR=\"0.5\"  MUF=\"1.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=0.5 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"8\" MUR=\"0.5\"  MUF=\"1.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=0.5 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"9\" MUR=\"0.5\"  MUF=\"1.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=0.5 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"10\" MUR=\"0.5\"  MUF=\"1.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=0.5 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"11\" MUR=\"0.5\"  MUF=\"2.0\"  PDF=\"247000\" > MUR=0.5 MUF=2.0  </weight>\\n'\n",
      "b'<weight id=\"12\" MUR=\"0.5\"  MUF=\"2.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=0.5 MUF=2.0 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"13\" MUR=\"0.5\"  MUF=\"2.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=0.5 MUF=2.0 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"14\" MUR=\"0.5\"  MUF=\"2.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=0.5 MUF=2.0 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"15\" MUR=\"0.5\"  MUF=\"2.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=0.5 MUF=2.0 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"16\" MUR=\"1.0\"  MUF=\"0.5\"  PDF=\"247000\" > MUF=0.5  </weight>\\n'\n",
      "b'<weight id=\"17\" MUR=\"1.0\"  MUF=\"0.5\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUF=0.5 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"18\" MUR=\"1.0\"  MUF=\"0.5\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUF=0.5 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"19\" MUR=\"1.0\"  MUF=\"0.5\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUF=0.5 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"20\" MUR=\"1.0\"  MUF=\"0.5\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUF=0.5 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"21\" MUR=\"1.0\"  MUF=\"1.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"22\" MUR=\"1.0\"  MUF=\"1.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"23\" MUR=\"1.0\"  MUF=\"1.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"24\" MUR=\"1.0\"  MUF=\"1.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"25\" MUR=\"1.0\"  MUF=\"2.0\"  PDF=\"247000\" > MUF=2.0  </weight>\\n'\n",
      "b'<weight id=\"26\" MUR=\"1.0\"  MUF=\"2.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUF=2.0 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"27\" MUR=\"1.0\"  MUF=\"2.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUF=2.0 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"28\" MUR=\"1.0\"  MUF=\"2.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUF=2.0 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"29\" MUR=\"1.0\"  MUF=\"2.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUF=2.0 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"30\" MUR=\"2.0\"  MUF=\"0.5\"  PDF=\"247000\" > MUR=2.0 MUF=0.5  </weight>\\n'\n",
      "b'<weight id=\"31\" MUR=\"2.0\"  MUF=\"0.5\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=2.0 MUF=0.5 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"32\" MUR=\"2.0\"  MUF=\"0.5\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=2.0 MUF=0.5 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"33\" MUR=\"2.0\"  MUF=\"0.5\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=2.0 MUF=0.5 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"34\" MUR=\"2.0\"  MUF=\"0.5\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=2.0 MUF=0.5 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"35\" MUR=\"2.0\"  MUF=\"1.0\"  PDF=\"247000\" > MUR=2.0  </weight>\\n'\n",
      "b'<weight id=\"36\" MUR=\"2.0\"  MUF=\"1.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=2.0 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"37\" MUR=\"2.0\"  MUF=\"1.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=2.0 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"38\" MUR=\"2.0\"  MUF=\"1.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=2.0 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"39\" MUR=\"2.0\"  MUF=\"1.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=2.0 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'<weight id=\"40\" MUR=\"2.0\"  MUF=\"2.0\"  PDF=\"247000\" > MUR=2.0 MUF=2.0  </weight>\\n'\n",
      "b'<weight id=\"41\" MUR=\"2.0\"  MUF=\"2.0\"  DYN_SCALE=\"1\"  PDF=\"247000\" > MUR=2.0 MUF=2.0 dyn_scale_choice=sum pt  </weight>\\n'\n",
      "b'<weight id=\"42\" MUR=\"2.0\"  MUF=\"2.0\"  DYN_SCALE=\"2\"  PDF=\"247000\" > MUR=2.0 MUF=2.0 dyn_scale_choice=HT  </weight>\\n'\n",
      "b'<weight id=\"43\" MUR=\"2.0\"  MUF=\"2.0\"  DYN_SCALE=\"3\"  PDF=\"247000\" > MUR=2.0 MUF=2.0 dyn_scale_choice=HT/2  </weight>\\n'\n",
      "b'<weight id=\"44\" MUR=\"2.0\"  MUF=\"2.0\"  DYN_SCALE=\"4\"  PDF=\"247000\" > MUR=2.0 MUF=2.0 dyn_scale_choice=sqrts  </weight>\\n'\n",
      "b'</weightgroup> # scale\\n'\n",
      "b'<weightgroup name=\"NNPDF23_lo_as_0130_qed\" combine=\"replicas\"> # 247000: LO QCD + LO QED PDF set, alphas(MZ)=0.130. mem=0 ; average on replicas; mem=1-100 ; PDF replicas\\n'\n",
      "b'<weight id=\"45\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247000\" >  </weight>\\n'\n",
      "b'<weight id=\"46\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247001\" > PDF=247000 MemberID=1 </weight>\\n'\n",
      "b'<weight id=\"47\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247002\" > PDF=247000 MemberID=2 </weight>\\n'\n",
      "b'<weight id=\"48\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247003\" > PDF=247000 MemberID=3 </weight>\\n'\n",
      "b'<weight id=\"49\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247004\" > PDF=247000 MemberID=4 </weight>\\n'\n",
      "b'<weight id=\"50\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247005\" > PDF=247000 MemberID=5 </weight>\\n'\n",
      "b'<weight id=\"51\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247006\" > PDF=247000 MemberID=6 </weight>\\n'\n",
      "b'<weight id=\"52\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247007\" > PDF=247000 MemberID=7 </weight>\\n'\n",
      "b'<weight id=\"53\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247008\" > PDF=247000 MemberID=8 </weight>\\n'\n",
      "b'<weight id=\"54\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247009\" > PDF=247000 MemberID=9 </weight>\\n'\n",
      "b'<weight id=\"55\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247010\" > PDF=247000 MemberID=10 </weight>\\n'\n",
      "b'<weight id=\"56\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247011\" > PDF=247000 MemberID=11 </weight>\\n'\n",
      "b'<weight id=\"57\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247012\" > PDF=247000 MemberID=12 </weight>\\n'\n",
      "b'<weight id=\"58\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247013\" > PDF=247000 MemberID=13 </weight>\\n'\n",
      "b'<weight id=\"59\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247014\" > PDF=247000 MemberID=14 </weight>\\n'\n",
      "b'<weight id=\"60\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247015\" > PDF=247000 MemberID=15 </weight>\\n'\n",
      "b'<weight id=\"61\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247016\" > PDF=247000 MemberID=16 </weight>\\n'\n",
      "b'<weight id=\"62\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247017\" > PDF=247000 MemberID=17 </weight>\\n'\n",
      "b'<weight id=\"63\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247018\" > PDF=247000 MemberID=18 </weight>\\n'\n",
      "b'<weight id=\"64\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247019\" > PDF=247000 MemberID=19 </weight>\\n'\n",
      "b'<weight id=\"65\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247020\" > PDF=247000 MemberID=20 </weight>\\n'\n",
      "b'<weight id=\"66\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247021\" > PDF=247000 MemberID=21 </weight>\\n'\n",
      "b'<weight id=\"67\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247022\" > PDF=247000 MemberID=22 </weight>\\n'\n",
      "b'<weight id=\"68\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247023\" > PDF=247000 MemberID=23 </weight>\\n'\n",
      "b'<weight id=\"69\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247024\" > PDF=247000 MemberID=24 </weight>\\n'\n",
      "b'<weight id=\"70\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247025\" > PDF=247000 MemberID=25 </weight>\\n'\n",
      "b'<weight id=\"71\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247026\" > PDF=247000 MemberID=26 </weight>\\n'\n",
      "b'<weight id=\"72\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247027\" > PDF=247000 MemberID=27 </weight>\\n'\n",
      "b'<weight id=\"73\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247028\" > PDF=247000 MemberID=28 </weight>\\n'\n",
      "b'<weight id=\"74\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247029\" > PDF=247000 MemberID=29 </weight>\\n'\n",
      "b'<weight id=\"75\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247030\" > PDF=247000 MemberID=30 </weight>\\n'\n",
      "b'<weight id=\"76\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247031\" > PDF=247000 MemberID=31 </weight>\\n'\n",
      "b'<weight id=\"77\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247032\" > PDF=247000 MemberID=32 </weight>\\n'\n",
      "b'<weight id=\"78\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247033\" > PDF=247000 MemberID=33 </weight>\\n'\n",
      "b'<weight id=\"79\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247034\" > PDF=247000 MemberID=34 </weight>\\n'\n",
      "b'<weight id=\"80\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247035\" > PDF=247000 MemberID=35 </weight>\\n'\n",
      "b'<weight id=\"81\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247036\" > PDF=247000 MemberID=36 </weight>\\n'\n",
      "b'<weight id=\"82\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247037\" > PDF=247000 MemberID=37 </weight>\\n'\n",
      "b'<weight id=\"83\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247038\" > PDF=247000 MemberID=38 </weight>\\n'\n",
      "b'<weight id=\"84\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247039\" > PDF=247000 MemberID=39 </weight>\\n'\n",
      "b'<weight id=\"85\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247040\" > PDF=247000 MemberID=40 </weight>\\n'\n",
      "b'<weight id=\"86\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247041\" > PDF=247000 MemberID=41 </weight>\\n'\n",
      "b'<weight id=\"87\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247042\" > PDF=247000 MemberID=42 </weight>\\n'\n",
      "b'<weight id=\"88\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247043\" > PDF=247000 MemberID=43 </weight>\\n'\n",
      "b'<weight id=\"89\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247044\" > PDF=247000 MemberID=44 </weight>\\n'\n",
      "b'<weight id=\"90\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247045\" > PDF=247000 MemberID=45 </weight>\\n'\n",
      "b'<weight id=\"91\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247046\" > PDF=247000 MemberID=46 </weight>\\n'\n",
      "b'<weight id=\"92\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247047\" > PDF=247000 MemberID=47 </weight>\\n'\n",
      "b'<weight id=\"93\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247048\" > PDF=247000 MemberID=48 </weight>\\n'\n",
      "b'<weight id=\"94\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247049\" > PDF=247000 MemberID=49 </weight>\\n'\n",
      "b'<weight id=\"95\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247050\" > PDF=247000 MemberID=50 </weight>\\n'\n",
      "b'<weight id=\"96\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247051\" > PDF=247000 MemberID=51 </weight>\\n'\n",
      "b'<weight id=\"97\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247052\" > PDF=247000 MemberID=52 </weight>\\n'\n",
      "b'<weight id=\"98\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247053\" > PDF=247000 MemberID=53 </weight>\\n'\n",
      "b'<weight id=\"99\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247054\" > PDF=247000 MemberID=54 </weight>\\n'\n",
      "b'<weight id=\"100\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247055\" > PDF=247000 MemberID=55 </weight>\\n'\n",
      "b'<weight id=\"101\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247056\" > PDF=247000 MemberID=56 </weight>\\n'\n",
      "b'<weight id=\"102\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247057\" > PDF=247000 MemberID=57 </weight>\\n'\n",
      "b'<weight id=\"103\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247058\" > PDF=247000 MemberID=58 </weight>\\n'\n",
      "b'<weight id=\"104\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247059\" > PDF=247000 MemberID=59 </weight>\\n'\n",
      "b'<weight id=\"105\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247060\" > PDF=247000 MemberID=60 </weight>\\n'\n",
      "b'<weight id=\"106\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247061\" > PDF=247000 MemberID=61 </weight>\\n'\n",
      "b'<weight id=\"107\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247062\" > PDF=247000 MemberID=62 </weight>\\n'\n",
      "b'<weight id=\"108\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247063\" > PDF=247000 MemberID=63 </weight>\\n'\n",
      "b'<weight id=\"109\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247064\" > PDF=247000 MemberID=64 </weight>\\n'\n",
      "b'<weight id=\"110\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247065\" > PDF=247000 MemberID=65 </weight>\\n'\n",
      "b'<weight id=\"111\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247066\" > PDF=247000 MemberID=66 </weight>\\n'\n",
      "b'<weight id=\"112\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247067\" > PDF=247000 MemberID=67 </weight>\\n'\n",
      "b'<weight id=\"113\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247068\" > PDF=247000 MemberID=68 </weight>\\n'\n",
      "b'<weight id=\"114\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247069\" > PDF=247000 MemberID=69 </weight>\\n'\n",
      "b'<weight id=\"115\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247070\" > PDF=247000 MemberID=70 </weight>\\n'\n",
      "b'<weight id=\"116\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247071\" > PDF=247000 MemberID=71 </weight>\\n'\n",
      "b'<weight id=\"117\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247072\" > PDF=247000 MemberID=72 </weight>\\n'\n",
      "b'<weight id=\"118\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247073\" > PDF=247000 MemberID=73 </weight>\\n'\n",
      "b'<weight id=\"119\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247074\" > PDF=247000 MemberID=74 </weight>\\n'\n",
      "b'<weight id=\"120\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247075\" > PDF=247000 MemberID=75 </weight>\\n'\n",
      "b'<weight id=\"121\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247076\" > PDF=247000 MemberID=76 </weight>\\n'\n",
      "b'<weight id=\"122\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247077\" > PDF=247000 MemberID=77 </weight>\\n'\n",
      "b'<weight id=\"123\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247078\" > PDF=247000 MemberID=78 </weight>\\n'\n",
      "b'<weight id=\"124\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247079\" > PDF=247000 MemberID=79 </weight>\\n'\n",
      "b'<weight id=\"125\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247080\" > PDF=247000 MemberID=80 </weight>\\n'\n",
      "b'<weight id=\"126\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247081\" > PDF=247000 MemberID=81 </weight>\\n'\n",
      "b'<weight id=\"127\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247082\" > PDF=247000 MemberID=82 </weight>\\n'\n",
      "b'<weight id=\"128\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247083\" > PDF=247000 MemberID=83 </weight>\\n'\n",
      "b'<weight id=\"129\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247084\" > PDF=247000 MemberID=84 </weight>\\n'\n",
      "b'<weight id=\"130\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247085\" > PDF=247000 MemberID=85 </weight>\\n'\n",
      "b'<weight id=\"131\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247086\" > PDF=247000 MemberID=86 </weight>\\n'\n",
      "b'<weight id=\"132\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247087\" > PDF=247000 MemberID=87 </weight>\\n'\n",
      "b'<weight id=\"133\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247088\" > PDF=247000 MemberID=88 </weight>\\n'\n",
      "b'<weight id=\"134\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247089\" > PDF=247000 MemberID=89 </weight>\\n'\n",
      "b'<weight id=\"135\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247090\" > PDF=247000 MemberID=90 </weight>\\n'\n",
      "b'<weight id=\"136\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247091\" > PDF=247000 MemberID=91 </weight>\\n'\n",
      "b'<weight id=\"137\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247092\" > PDF=247000 MemberID=92 </weight>\\n'\n",
      "b'<weight id=\"138\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247093\" > PDF=247000 MemberID=93 </weight>\\n'\n",
      "b'<weight id=\"139\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247094\" > PDF=247000 MemberID=94 </weight>\\n'\n",
      "b'<weight id=\"140\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247095\" > PDF=247000 MemberID=95 </weight>\\n'\n",
      "b'<weight id=\"141\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247096\" > PDF=247000 MemberID=96 </weight>\\n'\n",
      "b'<weight id=\"142\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247097\" > PDF=247000 MemberID=97 </weight>\\n'\n",
      "b'<weight id=\"143\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247098\" > PDF=247000 MemberID=98 </weight>\\n'\n",
      "b'<weight id=\"144\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247099\" > PDF=247000 MemberID=99 </weight>\\n'\n",
      "b'<weight id=\"145\" MUR=\"1.0\"  MUF=\"1.0\"  PDF=\"247100\" > PDF=247000 MemberID=100 </weight>\\n'\n",
      "b'</weightgroup>\\n'\n",
      "b'</initrwgt>\\n'\n",
      "b'<MGGenerationInfo>\\n'\n",
      "b'#  Number of Events        :       10000\\n'\n",
      "b'#  Integrated weight (pb)  :       1678.21\\n'\n",
      "b'</MGGenerationInfo>\\n'\n",
      "b'</header>\\n'\n",
      "b'<init>\\n'\n",
      "b'2212 2212 6.500000e+03 6.500000e+03 0 0 247000 247000 -4 1\\n'\n",
      "b'1.678210e+03 5.178922e+00 1.678210e+03 1\\n'\n",
      "b\"<generator name='MadGraph5_aMC@NLO' version='3.1.0'>please cite 1405.0301 </generator>\\n\"\n",
      "b'</init>\\n'\n",
      "b'<event>\\n'\n",
      "b' 4      1 +1.6782100e+03 8.93683800e+01 7.54677100e-03 1.30438800e-01\\n'\n",
      "b'        2 -1    0    0  501    0 +0.0000000000e+00 +0.0000000000e+00 +7.4318030872e+01 7.4318030872e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'       -2 -1    0    0    0  501 -0.0000000000e+00 -0.0000000000e+00 -2.6866655940e+01 2.6866655940e+01 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'      -11  1    1    2    0    0 -6.9396665975e+00 +2.7744220552e+01 -1.5147179771e+01 3.2362598807e+01 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'       11  1    1    2    0    0 +6.9396665975e+00 -2.7744220552e+01 +6.2598554703e+01 6.8822088005e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'<mgrwt>\\n'\n",
      "b'<rscale>  0 0.89368383E+02</rscale>\\n'\n",
      "b'<asrwt>0</asrwt>\\n'\n",
      "b'<pdfrwt beam=\"1\">  1        2 0.11433543E-01 0.89368383E+02</pdfrwt>\\n'\n",
      "b'<pdfrwt beam=\"2\">  1       -2 0.41333317E-02 0.89368383E+02</pdfrwt>\\n'\n",
      "b'<totfact> 0.11893024E+05</totfact>\\n'\n",
      "b'</mgrwt>\\n'\n",
      "b'<rwgt>\\n'\n",
      "b\"<wgt id='1'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='2'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='3'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='4'> +1.1295459e+03 </wgt>\\n\"\n",
      "b\"<wgt id='5'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='6'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='7'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='8'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='9'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='10'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='11'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='12'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='13'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='14'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='15'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='16'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='17'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='18'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='19'> +1.1295459e+03 </wgt>\\n\"\n",
      "b\"<wgt id='20'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='21'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='22'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='23'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='24'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='25'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='26'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='27'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='28'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='29'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='30'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='31'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='32'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='33'> +1.1295459e+03 </wgt>\\n\"\n",
      "b\"<wgt id='34'> +1.4809482e+03 </wgt>\\n\"\n",
      "b\"<wgt id='35'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='36'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='37'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='38'> +1.3469027e+03 </wgt>\\n\"\n",
      "b\"<wgt id='39'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='40'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='41'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='42'> +1.7450852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='43'> +1.5527937e+03 </wgt>\\n\"\n",
      "b\"<wgt id='44'> +1.8613709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='45'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='46'> +1.6553344e+03 </wgt>\\n\"\n",
      "b\"<wgt id='47'> +1.6314370e+03 </wgt>\\n\"\n",
      "b\"<wgt id='48'> +1.6461270e+03 </wgt>\\n\"\n",
      "b\"<wgt id='49'> +1.6199311e+03 </wgt>\\n\"\n",
      "b\"<wgt id='50'> +1.6842943e+03 </wgt>\\n\"\n",
      "b\"<wgt id='51'> +1.6670860e+03 </wgt>\\n\"\n",
      "b\"<wgt id='52'> +1.6093764e+03 </wgt>\\n\"\n",
      "b\"<wgt id='53'> +1.7530965e+03 </wgt>\\n\"\n",
      "b\"<wgt id='54'> +1.6819596e+03 </wgt>\\n\"\n",
      "b\"<wgt id='55'> +1.6515737e+03 </wgt>\\n\"\n",
      "b\"<wgt id='56'> +1.7195013e+03 </wgt>\\n\"\n",
      "b\"<wgt id='57'> +1.6270690e+03 </wgt>\\n\"\n",
      "b\"<wgt id='58'> +1.7513255e+03 </wgt>\\n\"\n",
      "b\"<wgt id='59'> +1.6538338e+03 </wgt>\\n\"\n",
      "b\"<wgt id='60'> +1.6744334e+03 </wgt>\\n\"\n",
      "b\"<wgt id='61'> +1.7554777e+03 </wgt>\\n\"\n",
      "b\"<wgt id='62'> +1.7073818e+03 </wgt>\\n\"\n",
      "b\"<wgt id='63'> +1.6950924e+03 </wgt>\\n\"\n",
      "b\"<wgt id='64'> +1.6950707e+03 </wgt>\\n\"\n",
      "b\"<wgt id='65'> +1.6721252e+03 </wgt>\\n\"\n",
      "b\"<wgt id='66'> +1.7156820e+03 </wgt>\\n\"\n",
      "b\"<wgt id='67'> +1.6861468e+03 </wgt>\\n\"\n",
      "b\"<wgt id='68'> +1.5079316e+03 </wgt>\\n\"\n",
      "b\"<wgt id='69'> +1.5807264e+03 </wgt>\\n\"\n",
      "b\"<wgt id='70'> +1.6787583e+03 </wgt>\\n\"\n",
      "b\"<wgt id='71'> +1.7109969e+03 </wgt>\\n\"\n",
      "b\"<wgt id='72'> +1.7204062e+03 </wgt>\\n\"\n",
      "b\"<wgt id='73'> +1.7409654e+03 </wgt>\\n\"\n",
      "b\"<wgt id='74'> +1.6151974e+03 </wgt>\\n\"\n",
      "b\"<wgt id='75'> +1.5948440e+03 </wgt>\\n\"\n",
      "b\"<wgt id='76'> +1.7636047e+03 </wgt>\\n\"\n",
      "b\"<wgt id='77'> +1.6925935e+03 </wgt>\\n\"\n",
      "b\"<wgt id='78'> +1.7551012e+03 </wgt>\\n\"\n",
      "b\"<wgt id='79'> +1.7148596e+03 </wgt>\\n\"\n",
      "b\"<wgt id='80'> +1.7482177e+03 </wgt>\\n\"\n",
      "b\"<wgt id='81'> +1.6604952e+03 </wgt>\\n\"\n",
      "b\"<wgt id='82'> +1.7277442e+03 </wgt>\\n\"\n",
      "b\"<wgt id='83'> +1.6063368e+03 </wgt>\\n\"\n",
      "b\"<wgt id='84'> +1.6893137e+03 </wgt>\\n\"\n",
      "b\"<wgt id='85'> +1.6910564e+03 </wgt>\\n\"\n",
      "b\"<wgt id='86'> +1.6646386e+03 </wgt>\\n\"\n",
      "b\"<wgt id='87'> +1.7112986e+03 </wgt>\\n\"\n",
      "b\"<wgt id='88'> +1.6774080e+03 </wgt>\\n\"\n",
      "b\"<wgt id='89'> +1.6865639e+03 </wgt>\\n\"\n",
      "b\"<wgt id='90'> +1.7343186e+03 </wgt>\\n\"\n",
      "b\"<wgt id='91'> +1.5960589e+03 </wgt>\\n\"\n",
      "b\"<wgt id='92'> +1.7313950e+03 </wgt>\\n\"\n",
      "b\"<wgt id='93'> +1.6164131e+03 </wgt>\\n\"\n",
      "b\"<wgt id='94'> +1.6859491e+03 </wgt>\\n\"\n",
      "b\"<wgt id='95'> +1.6586816e+03 </wgt>\\n\"\n",
      "b\"<wgt id='96'> +1.6679260e+03 </wgt>\\n\"\n",
      "b\"<wgt id='97'> +1.7226055e+03 </wgt>\\n\"\n",
      "b\"<wgt id='98'> +1.7132696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='99'> +1.7486811e+03 </wgt>\\n\"\n",
      "b\"<wgt id='100'> +1.7295602e+03 </wgt>\\n\"\n",
      "b\"<wgt id='101'> +1.5387633e+03 </wgt>\\n\"\n",
      "b\"<wgt id='102'> +1.5822066e+03 </wgt>\\n\"\n",
      "b\"<wgt id='103'> +1.6420998e+03 </wgt>\\n\"\n",
      "b\"<wgt id='104'> +1.6702810e+03 </wgt>\\n\"\n",
      "b\"<wgt id='105'> +1.7067868e+03 </wgt>\\n\"\n",
      "b\"<wgt id='106'> +1.7466994e+03 </wgt>\\n\"\n",
      "b\"<wgt id='107'> +1.7071306e+03 </wgt>\\n\"\n",
      "b\"<wgt id='108'> +1.7314133e+03 </wgt>\\n\"\n",
      "b\"<wgt id='109'> +1.7322694e+03 </wgt>\\n\"\n",
      "b\"<wgt id='110'> +1.7468692e+03 </wgt>\\n\"\n",
      "b\"<wgt id='111'> +1.7209840e+03 </wgt>\\n\"\n",
      "b\"<wgt id='112'> +1.6840758e+03 </wgt>\\n\"\n",
      "b\"<wgt id='113'> +1.6619731e+03 </wgt>\\n\"\n",
      "b\"<wgt id='114'> +1.7204600e+03 </wgt>\\n\"\n",
      "b\"<wgt id='115'> +1.6914023e+03 </wgt>\\n\"\n",
      "b\"<wgt id='116'> +1.6383799e+03 </wgt>\\n\"\n",
      "b\"<wgt id='117'> +1.6618700e+03 </wgt>\\n\"\n",
      "b\"<wgt id='118'> +1.6116126e+03 </wgt>\\n\"\n",
      "b\"<wgt id='119'> +1.7047949e+03 </wgt>\\n\"\n",
      "b\"<wgt id='120'> +1.6773679e+03 </wgt>\\n\"\n",
      "b\"<wgt id='121'> +1.7120135e+03 </wgt>\\n\"\n",
      "b\"<wgt id='122'> +1.7478342e+03 </wgt>\\n\"\n",
      "b\"<wgt id='123'> +1.7986132e+03 </wgt>\\n\"\n",
      "b\"<wgt id='124'> +1.5957709e+03 </wgt>\\n\"\n",
      "b\"<wgt id='125'> +1.5796458e+03 </wgt>\\n\"\n",
      "b\"<wgt id='126'> +1.6807832e+03 </wgt>\\n\"\n",
      "b\"<wgt id='127'> +1.6826222e+03 </wgt>\\n\"\n",
      "b\"<wgt id='128'> +1.6532672e+03 </wgt>\\n\"\n",
      "b\"<wgt id='129'> +1.6827477e+03 </wgt>\\n\"\n",
      "b\"<wgt id='130'> +1.5799663e+03 </wgt>\\n\"\n",
      "b\"<wgt id='131'> +1.6694572e+03 </wgt>\\n\"\n",
      "b\"<wgt id='132'> +1.5943567e+03 </wgt>\\n\"\n",
      "b\"<wgt id='133'> +1.6317081e+03 </wgt>\\n\"\n",
      "b\"<wgt id='134'> +1.6428968e+03 </wgt>\\n\"\n",
      "b\"<wgt id='135'> +1.6526160e+03 </wgt>\\n\"\n",
      "b\"<wgt id='136'> +1.6842128e+03 </wgt>\\n\"\n",
      "b\"<wgt id='137'> +1.7143852e+03 </wgt>\\n\"\n",
      "b\"<wgt id='138'> +1.6288776e+03 </wgt>\\n\"\n",
      "b\"<wgt id='139'> +1.5835530e+03 </wgt>\\n\"\n",
      "b\"<wgt id='140'> +1.6704753e+03 </wgt>\\n\"\n",
      "b\"<wgt id='141'> +1.6398635e+03 </wgt>\\n\"\n",
      "b\"<wgt id='142'> +1.7582902e+03 </wgt>\\n\"\n",
      "b\"<wgt id='143'> +1.7348300e+03 </wgt>\\n\"\n",
      "b\"<wgt id='144'> +1.7595976e+03 </wgt>\\n\"\n",
      "b\"<wgt id='145'> +1.7022906e+03 </wgt>\\n\"\n",
      "b'</rwgt>\\n'\n",
      "b'</event>\\n'\n",
      "b'<event>\\n'\n",
      "b' 5      1 +1.6782100e+03 9.16732600e+01 7.54677100e-03 1.29885000e-01\\n'\n",
      "b'        2 -1    0    0  501    0 +0.0000000000e+00 +0.0000000000e+00 +2.2528790807e+02 2.2528790807e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'       -2 -1    0    0    0  501 -0.0000000000e+00 -0.0000000000e+00 -9.3258299877e+00 9.3258299877e+00 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'       23  2    1    2    0    0 +0.0000000000e+00 +0.0000000000e+00 +2.1596207808e+02 2.3461373806e+02 9.1673261728e+01 0.0000e+00 0.0000e+00\\n'\n",
      "b'      -13  1    3    3    0    0 -3.2117860394e+01 -4.6108461100e+00 +2.5124142411e+01 4.1037048993e+01 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'       13  1    3    3    0    0 +3.2117860394e+01 +4.6108461100e+00 +1.9083793567e+02 1.9357668907e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'<mgrwt>\\n'\n",
      "b'<rscale>  0 0.91673262E+02</rscale>\\n'\n",
      "b'<asrwt>0</asrwt>\\n'\n",
      "b'<pdfrwt beam=\"1\">  1        2 0.34659674E-01 0.91673262E+02</pdfrwt>\\n'\n",
      "b'<pdfrwt beam=\"2\">  1       -2 0.14347432E-02 0.91673262E+02</pdfrwt>\\n'\n",
      "b'<totfact> 0.15920076E+05</totfact>\\n'\n",
      "b'</mgrwt>\\n'\n",
      "b'<rwgt>\\n'\n",
      "b\"<wgt id='1'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='2'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='3'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='4'> +1.1596394e+03 </wgt>\\n\"\n",
      "b\"<wgt id='5'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='6'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='7'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='8'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='9'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='10'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='11'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='12'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='13'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='14'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='15'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='16'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='17'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='18'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='19'> +1.1596394e+03 </wgt>\\n\"\n",
      "b\"<wgt id='20'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='21'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='22'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='23'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='24'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='25'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='26'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='27'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='28'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='29'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='30'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='31'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='32'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='33'> +1.1596394e+03 </wgt>\\n\"\n",
      "b\"<wgt id='34'> +1.4797373e+03 </wgt>\\n\"\n",
      "b\"<wgt id='35'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='36'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='37'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='38'> +1.3759668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='39'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='40'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='41'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='42'> +1.7725791e+03 </wgt>\\n\"\n",
      "b\"<wgt id='43'> +1.5809763e+03 </wgt>\\n\"\n",
      "b\"<wgt id='44'> +1.8628616e+03 </wgt>\\n\"\n",
      "b\"<wgt id='45'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='46'> +1.6505391e+03 </wgt>\\n\"\n",
      "b\"<wgt id='47'> +1.5966108e+03 </wgt>\\n\"\n",
      "b\"<wgt id='48'> +1.6772701e+03 </wgt>\\n\"\n",
      "b\"<wgt id='49'> +1.6641707e+03 </wgt>\\n\"\n",
      "b\"<wgt id='50'> +1.6234566e+03 </wgt>\\n\"\n",
      "b\"<wgt id='51'> +1.7019117e+03 </wgt>\\n\"\n",
      "b\"<wgt id='52'> +1.5674475e+03 </wgt>\\n\"\n",
      "b\"<wgt id='53'> +1.7235693e+03 </wgt>\\n\"\n",
      "b\"<wgt id='54'> +1.6792830e+03 </wgt>\\n\"\n",
      "b\"<wgt id='55'> +1.6923323e+03 </wgt>\\n\"\n",
      "b\"<wgt id='56'> +1.7390786e+03 </wgt>\\n\"\n",
      "b\"<wgt id='57'> +1.6467857e+03 </wgt>\\n\"\n",
      "b\"<wgt id='58'> +1.7388782e+03 </wgt>\\n\"\n",
      "b\"<wgt id='59'> +1.6653197e+03 </wgt>\\n\"\n",
      "b\"<wgt id='60'> +1.7134374e+03 </wgt>\\n\"\n",
      "b\"<wgt id='61'> +1.7874925e+03 </wgt>\\n\"\n",
      "b\"<wgt id='62'> +1.7295074e+03 </wgt>\\n\"\n",
      "b\"<wgt id='63'> +1.6781383e+03 </wgt>\\n\"\n",
      "b\"<wgt id='64'> +1.7137739e+03 </wgt>\\n\"\n",
      "b\"<wgt id='65'> +1.6368330e+03 </wgt>\\n\"\n",
      "b\"<wgt id='66'> +1.7568244e+03 </wgt>\\n\"\n",
      "b\"<wgt id='67'> +1.6588736e+03 </wgt>\\n\"\n",
      "b\"<wgt id='68'> +1.5083516e+03 </wgt>\\n\"\n",
      "b\"<wgt id='69'> +1.6070216e+03 </wgt>\\n\"\n",
      "b\"<wgt id='70'> +1.6595849e+03 </wgt>\\n\"\n",
      "b\"<wgt id='71'> +1.7594656e+03 </wgt>\\n\"\n",
      "b\"<wgt id='72'> +1.7681309e+03 </wgt>\\n\"\n",
      "b\"<wgt id='73'> +1.7305205e+03 </wgt>\\n\"\n",
      "b\"<wgt id='74'> +1.5991613e+03 </wgt>\\n\"\n",
      "b\"<wgt id='75'> +1.5963959e+03 </wgt>\\n\"\n",
      "b\"<wgt id='76'> +1.7571922e+03 </wgt>\\n\"\n",
      "b\"<wgt id='77'> +1.7147748e+03 </wgt>\\n\"\n",
      "b\"<wgt id='78'> +1.7526111e+03 </wgt>\\n\"\n",
      "b\"<wgt id='79'> +1.7185473e+03 </wgt>\\n\"\n",
      "b\"<wgt id='80'> +1.7406068e+03 </wgt>\\n\"\n",
      "b\"<wgt id='81'> +1.6575406e+03 </wgt>\\n\"\n",
      "b\"<wgt id='82'> +1.7406160e+03 </wgt>\\n\"\n",
      "b\"<wgt id='83'> +1.5878548e+03 </wgt>\\n\"\n",
      "b\"<wgt id='84'> +1.6334682e+03 </wgt>\\n\"\n",
      "b\"<wgt id='85'> +1.6392840e+03 </wgt>\\n\"\n",
      "b\"<wgt id='86'> +1.6191189e+03 </wgt>\\n\"\n",
      "b\"<wgt id='87'> +1.6977026e+03 </wgt>\\n\"\n",
      "b\"<wgt id='88'> +1.6717682e+03 </wgt>\\n\"\n",
      "b\"<wgt id='89'> +1.6228405e+03 </wgt>\\n\"\n",
      "b\"<wgt id='90'> +1.6865946e+03 </wgt>\\n\"\n",
      "b\"<wgt id='91'> +1.6454314e+03 </wgt>\\n\"\n",
      "b\"<wgt id='92'> +1.7344132e+03 </wgt>\\n\"\n",
      "b\"<wgt id='93'> +1.6032781e+03 </wgt>\\n\"\n",
      "b\"<wgt id='94'> +1.6507470e+03 </wgt>\\n\"\n",
      "b\"<wgt id='95'> +1.6654721e+03 </wgt>\\n\"\n",
      "b\"<wgt id='96'> +1.6258657e+03 </wgt>\\n\"\n",
      "b\"<wgt id='97'> +1.7065981e+03 </wgt>\\n\"\n",
      "b\"<wgt id='98'> +1.6883326e+03 </wgt>\\n\"\n",
      "b\"<wgt id='99'> +1.7888981e+03 </wgt>\\n\"\n",
      "b\"<wgt id='100'> +1.8002926e+03 </wgt>\\n\"\n",
      "b\"<wgt id='101'> +1.5502967e+03 </wgt>\\n\"\n",
      "b\"<wgt id='102'> +1.5898092e+03 </wgt>\\n\"\n",
      "b\"<wgt id='103'> +1.5835438e+03 </wgt>\\n\"\n",
      "b\"<wgt id='104'> +1.6640150e+03 </wgt>\\n\"\n",
      "b\"<wgt id='105'> +1.6889330e+03 </wgt>\\n\"\n",
      "b\"<wgt id='106'> +1.7654683e+03 </wgt>\\n\"\n",
      "b\"<wgt id='107'> +1.7128100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='108'> +1.7757611e+03 </wgt>\\n\"\n",
      "b\"<wgt id='109'> +1.7580019e+03 </wgt>\\n\"\n",
      "b\"<wgt id='110'> +1.7506678e+03 </wgt>\\n\"\n",
      "b\"<wgt id='111'> +1.7399736e+03 </wgt>\\n\"\n",
      "b\"<wgt id='112'> +1.7433282e+03 </wgt>\\n\"\n",
      "b\"<wgt id='113'> +1.7098484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='114'> +1.7462820e+03 </wgt>\\n\"\n",
      "b\"<wgt id='115'> +1.6987584e+03 </wgt>\\n\"\n",
      "b\"<wgt id='116'> +1.6599069e+03 </wgt>\\n\"\n",
      "b\"<wgt id='117'> +1.6441628e+03 </wgt>\\n\"\n",
      "b\"<wgt id='118'> +1.6224640e+03 </wgt>\\n\"\n",
      "b\"<wgt id='119'> +1.6878609e+03 </wgt>\\n\"\n",
      "b\"<wgt id='120'> +1.6673198e+03 </wgt>\\n\"\n",
      "b\"<wgt id='121'> +1.7409397e+03 </wgt>\\n\"\n",
      "b\"<wgt id='122'> +1.7359414e+03 </wgt>\\n\"\n",
      "b\"<wgt id='123'> +1.7596464e+03 </wgt>\\n\"\n",
      "b\"<wgt id='124'> +1.5766669e+03 </wgt>\\n\"\n",
      "b\"<wgt id='125'> +1.5483232e+03 </wgt>\\n\"\n",
      "b\"<wgt id='126'> +1.6550804e+03 </wgt>\\n\"\n",
      "b\"<wgt id='127'> +1.6488183e+03 </wgt>\\n\"\n",
      "b\"<wgt id='128'> +1.6913894e+03 </wgt>\\n\"\n",
      "b\"<wgt id='129'> +1.6360688e+03 </wgt>\\n\"\n",
      "b\"<wgt id='130'> +1.6231379e+03 </wgt>\\n\"\n",
      "b\"<wgt id='131'> +1.6348874e+03 </wgt>\\n\"\n",
      "b\"<wgt id='132'> +1.5943627e+03 </wgt>\\n\"\n",
      "b\"<wgt id='133'> +1.6045374e+03 </wgt>\\n\"\n",
      "b\"<wgt id='134'> +1.6687288e+03 </wgt>\\n\"\n",
      "b\"<wgt id='135'> +1.6728018e+03 </wgt>\\n\"\n",
      "b\"<wgt id='136'> +1.6806901e+03 </wgt>\\n\"\n",
      "b\"<wgt id='137'> +1.7128673e+03 </wgt>\\n\"\n",
      "b\"<wgt id='138'> +1.6707171e+03 </wgt>\\n\"\n",
      "b\"<wgt id='139'> +1.5442689e+03 </wgt>\\n\"\n",
      "b\"<wgt id='140'> +1.6675483e+03 </wgt>\\n\"\n",
      "b\"<wgt id='141'> +1.6382448e+03 </wgt>\\n\"\n",
      "b\"<wgt id='142'> +1.8170021e+03 </wgt>\\n\"\n",
      "b\"<wgt id='143'> +1.7308880e+03 </wgt>\\n\"\n",
      "b\"<wgt id='144'> +1.7507165e+03 </wgt>\\n\"\n",
      "b\"<wgt id='145'> +1.6813896e+03 </wgt>\\n\"\n",
      "b'</rwgt>\\n'\n",
      "b'</event>\\n'\n",
      "b'<event>\\n'\n",
      "b' 5      1 +1.6782100e+03 9.33757300e+01 7.54677100e-03 1.29487700e-01\\n'\n",
      "b'       -2 -1    0    0    0  501 -0.0000000000e+00 +0.0000000000e+00 +6.3859886915e+01 6.3859886915e+01 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'        2 -1    0    0  501    0 +0.0000000000e+00 -0.0000000000e+00 -3.4133425006e+01 3.4133425006e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'       23  2    1    2    0    0 +0.0000000000e+00 +0.0000000000e+00 +2.9726461909e+01 9.7993311922e+01 9.3375728344e+01 0.0000e+00 0.0000e+00\\n'\n",
      "b'      -11  1    3    3    0    0 +1.7879950015e+01 -4.3280842501e+00 -3.0169517385e+01 3.5335884099e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00\\n'\n",
      "b'       11  1    3    3    0    0 -1.7879950015e+01 +4.3280842501e+00 +5.9895979294e+01 6.2657427823e+01 0.0000000000e+00 0.0000e+00 1.0000e+00\\n'\n",
      "b'<mgrwt>\\n'\n",
      "b'<rscale>  0 0.93375728E+02</rscale>\\n'\n",
      "b'<asrwt>0</asrwt>\\n'\n",
      "b'<pdfrwt beam=\"2\">  1        2 0.52512961E-02 0.93375728E+02</pdfrwt>\\n'\n",
      "b'<pdfrwt beam=\"1\">  1       -2 0.98245980E-02 0.93375728E+02</pdfrwt>\\n'\n",
      "b'<totfact> 0.85331228E+04</totfact>\\n'\n",
      "b'</mgrwt>\\n'\n",
      "b'<rwgt>\\n'\n",
      "b\"<wgt id='1'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='2'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='3'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='4'> +9.9299050e+02 </wgt>\\n\"\n",
      "b\"<wgt id='5'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='6'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='7'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='8'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='9'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='10'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='11'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='12'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='13'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='14'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='15'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='16'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='17'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='18'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='19'> +9.9299050e+02 </wgt>\\n\"\n",
      "b\"<wgt id='20'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='21'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='22'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='23'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='24'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='25'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='26'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='27'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='28'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='29'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='30'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='31'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='32'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='33'> +9.9299050e+02 </wgt>\\n\"\n",
      "b\"<wgt id='34'> +1.4883985e+03 </wgt>\\n\"\n",
      "b\"<wgt id='35'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='36'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='37'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='38'> +1.2114790e+03 </wgt>\\n\"\n",
      "b\"<wgt id='39'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='40'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='41'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='42'> +1.6146484e+03 </wgt>\\n\"\n",
      "b\"<wgt id='43'> +1.4198331e+03 </wgt>\\n\"\n",
      "b\"<wgt id='44'> +1.8531696e+03 </wgt>\\n\"\n",
      "b\"<wgt id='45'> +1.6782100e+03 </wgt>\\n\"\n",
      "b\"<wgt id='46'> +1.6756355e+03 </wgt>\\n\"\n",
      "b\"<wgt id='47'> +1.6361998e+03 </wgt>\\n\"\n",
      "b\"<wgt id='48'> +1.6303713e+03 </wgt>\\n\"\n",
      "b\"<wgt id='49'> +1.6025788e+03 </wgt>\\n\"\n",
      "b\"<wgt id='50'> +1.6586893e+03 </wgt>\\n\"\n",
      "b\"<wgt id='51'> +1.6492306e+03 </wgt>\\n\"\n",
      "b\"<wgt id='52'> +1.6154644e+03 </wgt>\\n\"\n",
      "b\"<wgt id='53'> +1.7549898e+03 </wgt>\\n\"\n",
      "b\"<wgt id='54'> +1.6698530e+03 </wgt>\\n\"\n",
      "b\"<wgt id='55'> +1.6125796e+03 </wgt>\\n\"\n",
      "b\"<wgt id='56'> +1.7880179e+03 </wgt>\\n\"\n",
      "b\"<wgt id='57'> +1.6069682e+03 </wgt>\\n\"\n",
      "b\"<wgt id='58'> +1.7433293e+03 </wgt>\\n\"\n",
      "b\"<wgt id='59'> +1.6409517e+03 </wgt>\\n\"\n",
      "b\"<wgt id='60'> +1.6602141e+03 </wgt>\\n\"\n",
      "b\"<wgt id='61'> +1.7573271e+03 </wgt>\\n\"\n",
      "b\"<wgt id='62'> +1.7026889e+03 </wgt>\\n\"\n",
      "b\"<wgt id='63'> +1.7251596e+03 </wgt>\\n\"\n",
      "b\"<wgt id='64'> +1.6704608e+03 </wgt>\\n\"\n",
      "b\"<wgt id='65'> +1.6768926e+03 </wgt>\\n\"\n",
      "b\"<wgt id='66'> +1.7026783e+03 </wgt>\\n\"\n",
      "b\"<wgt id='67'> +1.7205159e+03 </wgt>\\n\"\n",
      "b\"<wgt id='68'> +1.4845312e+03 </wgt>\\n\"\n",
      "b\"<wgt id='69'> +1.5523637e+03 </wgt>\\n\"\n",
      "b\"<wgt id='70'> +1.6813983e+03 </wgt>\\n\"\n",
      "b\"<wgt id='71'> +1.6851629e+03 </wgt>\\n\"\n",
      "b\"<wgt id='72'> +1.6971995e+03 </wgt>\\n\"\n",
      "b\"<wgt id='73'> +1.7452855e+03 </wgt>\\n\"\n",
      "b\"<wgt id='74'> +1.5968393e+03 </wgt>\\n\"\n",
      "b\"<wgt id='75'> +1.6141849e+03 </wgt>\\n\"\n",
      "b\"<wgt id='76'> +1.7522634e+03 </wgt>\\n\"\n",
      "b\"<wgt id='77'> +1.6797470e+03 </wgt>\\n\"\n",
      "b\"<wgt id='78'> +1.7517570e+03 </wgt>\\n\"\n",
      "b\"<wgt id='79'> +1.7070047e+03 </wgt>\\n\"\n",
      "b\"<wgt id='80'> +1.7658384e+03 </wgt>\\n\"\n",
      "b\"<wgt id='81'> +1.6692237e+03 </wgt>\\n\"\n",
      "b\"<wgt id='82'> +1.7025974e+03 </wgt>\\n\"\n",
      "b\"<wgt id='83'> +1.6140583e+03 </wgt>\\n\"\n",
      "b\"<wgt id='84'> +1.7053940e+03 </wgt>\\n\"\n",
      "b\"<wgt id='85'> +1.7140653e+03 </wgt>\\n\"\n",
      "b\"<wgt id='86'> +1.6910241e+03 </wgt>\\n\"\n",
      "b\"<wgt id='87'> +1.7080953e+03 </wgt>\\n\"\n",
      "b\"<wgt id='88'> +1.6750669e+03 </wgt>\\n\"\n",
      "b\"<wgt id='89'> +1.7555844e+03 </wgt>\\n\"\n",
      "b\"<wgt id='90'> +1.8098533e+03 </wgt>\\n\"\n",
      "b\"<wgt id='91'> +1.5708833e+03 </wgt>\\n\"\n",
      "b\"<wgt id='92'> +1.7438875e+03 </wgt>\\n\"\n",
      "b\"<wgt id='93'> +1.6318138e+03 </wgt>\\n\"\n",
      "b\"<wgt id='94'> +1.7275750e+03 </wgt>\\n\"\n",
      "b\"<wgt id='95'> +1.6557302e+03 </wgt>\\n\"\n",
      "b\"<wgt id='96'> +1.6843102e+03 </wgt>\\n\"\n",
      "b\"<wgt id='97'> +1.7396353e+03 </wgt>\\n\"\n",
      "b\"<wgt id='98'> +1.7271012e+03 </wgt>\\n\"\n",
      "b\"<wgt id='99'> +1.7353544e+03 </wgt>\\n\"\n",
      "b\"<wgt id='100'> +1.7142687e+03 </wgt>\\n\"\n",
      "b\"<wgt id='101'> +1.5204372e+03 </wgt>\\n\"\n",
      "b\"<wgt id='102'> +1.5724589e+03 </wgt>\\n\"\n",
      "b\"<wgt id='103'> +1.6755699e+03 </wgt>\\n\"\n",
      "b\"<wgt id='104'> +1.6666335e+03 </wgt>\\n\"\n",
      "b\"<wgt id='105'> +1.7045745e+03 </wgt>\\n\"\n",
      "b\"<wgt id='106'> +1.7545989e+03 </wgt>\\n\"\n",
      "b\"<wgt id='107'> +1.6874663e+03 </wgt>\\n\"\n",
      "b\"<wgt id='108'> +1.7353972e+03 </wgt>\\n\"\n",
      "b\"<wgt id='109'> +1.7173414e+03 </wgt>\\n\"\n",
      "b\"<wgt id='110'> +1.7495491e+03 </wgt>\\n\"\n",
      "b\"<wgt id='111'> +1.7122079e+03 </wgt>\\n\"\n",
      "b\"<wgt id='112'> +1.6609661e+03 </wgt>\\n\"\n",
      "b\"<wgt id='113'> +1.6467599e+03 </wgt>\\n\"\n",
      "b\"<wgt id='114'> +1.6984213e+03 </wgt>\\n\"\n",
      "b\"<wgt id='115'> +1.6933318e+03 </wgt>\\n\"\n",
      "b\"<wgt id='116'> +1.6180121e+03 </wgt>\\n\"\n",
      "b\"<wgt id='117'> +1.6804048e+03 </wgt>\\n\"\n",
      "b\"<wgt id='118'> +1.5969053e+03 </wgt>\\n\"\n",
      "b\"<wgt id='119'> +1.7142371e+03 </wgt>\\n\"\n",
      "b\"<wgt id='120'> +1.6758778e+03 </wgt>\\n\"\n",
      "b\"<wgt id='121'> +1.6946474e+03 </wgt>\\n\"\n",
      "b\"<wgt id='122'> +1.7439471e+03 </wgt>\\n\"\n",
      "b\"<wgt id='123'> +1.8105080e+03 </wgt>\\n\"\n",
      "b\"<wgt id='124'> +1.5910868e+03 </wgt>\\n\"\n",
      "b\"<wgt id='125'> +1.5705175e+03 </wgt>\\n\"\n",
      "b\"<wgt id='126'> +1.6980770e+03 </wgt>\\n\"\n",
      "b\"<wgt id='127'> +1.7072316e+03 </wgt>\\n\"\n",
      "b\"<wgt id='128'> +1.6440392e+03 </wgt>\\n\"\n",
      "b\"<wgt id='129'> +1.6651287e+03 </wgt>\\n\"\n",
      "b\"<wgt id='130'> +1.5545586e+03 </wgt>\\n\"\n",
      "b\"<wgt id='131'> +1.6787821e+03 </wgt>\\n\"\n",
      "b\"<wgt id='132'> +1.5967667e+03 </wgt>\\n\"\n",
      "b\"<wgt id='133'> +1.6781010e+03 </wgt>\\n\"\n",
      "b\"<wgt id='134'> +1.6159887e+03 </wgt>\\n\"\n",
      "b\"<wgt id='135'> +1.6230268e+03 </wgt>\\n\"\n",
      "b\"<wgt id='136'> +1.6648668e+03 </wgt>\\n\"\n",
      "b\"<wgt id='137'> +1.7231546e+03 </wgt>\\n\"\n",
      "b\"<wgt id='138'> +1.6322453e+03 </wgt>\\n\"\n"
     ]
    }
   ],
   "source": [
    "lhe_data = data_path(\"pylhe-drell-yan-ll-lhe.gz\")\n",
    "with gzip.open(lhe_data) as f:\n",
    "    for _ in range(1000):\n",
    "        print(f.readline())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1b0954bf-4d67-4ddc-8496-052f508135d5",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<LesHouchesEvents version=\"3.0\">\n",
      "<init>\n",
      "   2212   2212  6.5000000e+03  6.5000000e+03     0     0  247000  247000    -4     1\n",
      " 1.6782100e+03  5.1789220e+00  1.6782100e+03     1\n",
      "<initrwgt>\n",
      "  <weightgroup name=\"Central scale variation\" combine=\"envelope\">\n",
      "    <weight id=\"1\" MUR=\"0.5\" MUF=\"0.5\" PDF=\"247000\">MUR=0.5 MUF=0.5</weight>\n",
      "    <weight id=\"2\" MUR=\"0.5\" MUF=\"0.5\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=0.5 MUF=0.5 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"3\" MUR=\"0.5\" MUF=\"0.5\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=0.5 MUF=0.5 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"4\" MUR=\"0.5\" MUF=\"0.5\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=0.5 MUF=0.5 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"5\" MUR=\"0.5\" MUF=\"0.5\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=0.5 MUF=0.5 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"6\" MUR=\"0.5\" MUF=\"1.0\" PDF=\"247000\">MUR=0.5</weight>\n",
      "    <weight id=\"7\" MUR=\"0.5\" MUF=\"1.0\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=0.5 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"8\" MUR=\"0.5\" MUF=\"1.0\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=0.5 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"9\" MUR=\"0.5\" MUF=\"1.0\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=0.5 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"10\" MUR=\"0.5\" MUF=\"1.0\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=0.5 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"11\" MUR=\"0.5\" MUF=\"2.0\" PDF=\"247000\">MUR=0.5 MUF=2.0</weight>\n",
      "    <weight id=\"12\" MUR=\"0.5\" MUF=\"2.0\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=0.5 MUF=2.0 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"13\" MUR=\"0.5\" MUF=\"2.0\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=0.5 MUF=2.0 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"14\" MUR=\"0.5\" MUF=\"2.0\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=0.5 MUF=2.0 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"15\" MUR=\"0.5\" MUF=\"2.0\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=0.5 MUF=2.0 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"16\" MUR=\"1.0\" MUF=\"0.5\" PDF=\"247000\">MUF=0.5</weight>\n",
      "    <weight id=\"17\" MUR=\"1.0\" MUF=\"0.5\" DYN_SCALE=\"1\" PDF=\"247000\">MUF=0.5 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"18\" MUR=\"1.0\" MUF=\"0.5\" DYN_SCALE=\"2\" PDF=\"247000\">MUF=0.5 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"19\" MUR=\"1.0\" MUF=\"0.5\" DYN_SCALE=\"3\" PDF=\"247000\">MUF=0.5 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"20\" MUR=\"1.0\" MUF=\"0.5\" DYN_SCALE=\"4\" PDF=\"247000\">MUF=0.5 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"21\" MUR=\"1.0\" MUF=\"1.0\" DYN_SCALE=\"1\" PDF=\"247000\">dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"22\" MUR=\"1.0\" MUF=\"1.0\" DYN_SCALE=\"2\" PDF=\"247000\">dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"23\" MUR=\"1.0\" MUF=\"1.0\" DYN_SCALE=\"3\" PDF=\"247000\">dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"24\" MUR=\"1.0\" MUF=\"1.0\" DYN_SCALE=\"4\" PDF=\"247000\">dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"25\" MUR=\"1.0\" MUF=\"2.0\" PDF=\"247000\">MUF=2.0</weight>\n",
      "    <weight id=\"26\" MUR=\"1.0\" MUF=\"2.0\" DYN_SCALE=\"1\" PDF=\"247000\">MUF=2.0 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"27\" MUR=\"1.0\" MUF=\"2.0\" DYN_SCALE=\"2\" PDF=\"247000\">MUF=2.0 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"28\" MUR=\"1.0\" MUF=\"2.0\" DYN_SCALE=\"3\" PDF=\"247000\">MUF=2.0 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"29\" MUR=\"1.0\" MUF=\"2.0\" DYN_SCALE=\"4\" PDF=\"247000\">MUF=2.0 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"30\" MUR=\"2.0\" MUF=\"0.5\" PDF=\"247000\">MUR=2.0 MUF=0.5</weight>\n",
      "    <weight id=\"31\" MUR=\"2.0\" MUF=\"0.5\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=2.0 MUF=0.5 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"32\" MUR=\"2.0\" MUF=\"0.5\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=2.0 MUF=0.5 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"33\" MUR=\"2.0\" MUF=\"0.5\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=2.0 MUF=0.5 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"34\" MUR=\"2.0\" MUF=\"0.5\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=2.0 MUF=0.5 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"35\" MUR=\"2.0\" MUF=\"1.0\" PDF=\"247000\">MUR=2.0</weight>\n",
      "    <weight id=\"36\" MUR=\"2.0\" MUF=\"1.0\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=2.0 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"37\" MUR=\"2.0\" MUF=\"1.0\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=2.0 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"38\" MUR=\"2.0\" MUF=\"1.0\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=2.0 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"39\" MUR=\"2.0\" MUF=\"1.0\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=2.0 dyn_scale_choice=sqrts</weight>\n",
      "    <weight id=\"40\" MUR=\"2.0\" MUF=\"2.0\" PDF=\"247000\">MUR=2.0 MUF=2.0</weight>\n",
      "    <weight id=\"41\" MUR=\"2.0\" MUF=\"2.0\" DYN_SCALE=\"1\" PDF=\"247000\">MUR=2.0 MUF=2.0 dyn_scale_choice=sum pt</weight>\n",
      "    <weight id=\"42\" MUR=\"2.0\" MUF=\"2.0\" DYN_SCALE=\"2\" PDF=\"247000\">MUR=2.0 MUF=2.0 dyn_scale_choice=HT</weight>\n",
      "    <weight id=\"43\" MUR=\"2.0\" MUF=\"2.0\" DYN_SCALE=\"3\" PDF=\"247000\">MUR=2.0 MUF=2.0 dyn_scale_choice=HT/2</weight>\n",
      "    <weight id=\"44\" MUR=\"2.0\" MUF=\"2.0\" DYN_SCALE=\"4\" PDF=\"247000\">MUR=2.0 MUF=2.0 dyn_scale_choice=sqrts</weight>\n",
      "  </weightgroup>\n",
      "  <weightgroup name=\"NNPDF23_lo_as_0130_qed\" combine=\"replicas\">\n",
      "    <weight id=\"45\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247000\" />\n",
      "    <weight id=\"46\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247001\">PDF=247000 MemberID=1</weight>\n",
      "    <weight id=\"47\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247002\">PDF=247000 MemberID=2</weight>\n",
      "    <weight id=\"48\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247003\">PDF=247000 MemberID=3</weight>\n",
      "    <weight id=\"49\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247004\">PDF=247000 MemberID=4</weight>\n",
      "    <weight id=\"50\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247005\">PDF=247000 MemberID=5</weight>\n",
      "    <weight id=\"51\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247006\">PDF=247000 MemberID=6</weight>\n",
      "    <weight id=\"52\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247007\">PDF=247000 MemberID=7</weight>\n",
      "    <weight id=\"53\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247008\">PDF=247000 MemberID=8</weight>\n",
      "    <weight id=\"54\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247009\">PDF=247000 MemberID=9</weight>\n",
      "    <weight id=\"55\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247010\">PDF=247000 MemberID=10</weight>\n",
      "    <weight id=\"56\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247011\">PDF=247000 MemberID=11</weight>\n",
      "    <weight id=\"57\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247012\">PDF=247000 MemberID=12</weight>\n",
      "    <weight id=\"58\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247013\">PDF=247000 MemberID=13</weight>\n",
      "    <weight id=\"59\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247014\">PDF=247000 MemberID=14</weight>\n",
      "    <weight id=\"60\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247015\">PDF=247000 MemberID=15</weight>\n",
      "    <weight id=\"61\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247016\">PDF=247000 MemberID=16</weight>\n",
      "    <weight id=\"62\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247017\">PDF=247000 MemberID=17</weight>\n",
      "    <weight id=\"63\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247018\">PDF=247000 MemberID=18</weight>\n",
      "    <weight id=\"64\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247019\">PDF=247000 MemberID=19</weight>\n",
      "    <weight id=\"65\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247020\">PDF=247000 MemberID=20</weight>\n",
      "    <weight id=\"66\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247021\">PDF=247000 MemberID=21</weight>\n",
      "    <weight id=\"67\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247022\">PDF=247000 MemberID=22</weight>\n",
      "    <weight id=\"68\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247023\">PDF=247000 MemberID=23</weight>\n",
      "    <weight id=\"69\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247024\">PDF=247000 MemberID=24</weight>\n",
      "    <weight id=\"70\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247025\">PDF=247000 MemberID=25</weight>\n",
      "    <weight id=\"71\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247026\">PDF=247000 MemberID=26</weight>\n",
      "    <weight id=\"72\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247027\">PDF=247000 MemberID=27</weight>\n",
      "    <weight id=\"73\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247028\">PDF=247000 MemberID=28</weight>\n",
      "    <weight id=\"74\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247029\">PDF=247000 MemberID=29</weight>\n",
      "    <weight id=\"75\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247030\">PDF=247000 MemberID=30</weight>\n",
      "    <weight id=\"76\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247031\">PDF=247000 MemberID=31</weight>\n",
      "    <weight id=\"77\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247032\">PDF=247000 MemberID=32</weight>\n",
      "    <weight id=\"78\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247033\">PDF=247000 MemberID=33</weight>\n",
      "    <weight id=\"79\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247034\">PDF=247000 MemberID=34</weight>\n",
      "    <weight id=\"80\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247035\">PDF=247000 MemberID=35</weight>\n",
      "    <weight id=\"81\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247036\">PDF=247000 MemberID=36</weight>\n",
      "    <weight id=\"82\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247037\">PDF=247000 MemberID=37</weight>\n",
      "    <weight id=\"83\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247038\">PDF=247000 MemberID=38</weight>\n",
      "    <weight id=\"84\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247039\">PDF=247000 MemberID=39</weight>\n",
      "    <weight id=\"85\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247040\">PDF=247000 MemberID=40</weight>\n",
      "    <weight id=\"86\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247041\">PDF=247000 MemberID=41</weight>\n",
      "    <weight id=\"87\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247042\">PDF=247000 MemberID=42</weight>\n",
      "    <weight id=\"88\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247043\">PDF=247000 MemberID=43</weight>\n",
      "    <weight id=\"89\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247044\">PDF=247000 MemberID=44</weight>\n",
      "    <weight id=\"90\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247045\">PDF=247000 MemberID=45</weight>\n",
      "    <weight id=\"91\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247046\">PDF=247000 MemberID=46</weight>\n",
      "    <weight id=\"92\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247047\">PDF=247000 MemberID=47</weight>\n",
      "    <weight id=\"93\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247048\">PDF=247000 MemberID=48</weight>\n",
      "    <weight id=\"94\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247049\">PDF=247000 MemberID=49</weight>\n",
      "    <weight id=\"95\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247050\">PDF=247000 MemberID=50</weight>\n",
      "    <weight id=\"96\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247051\">PDF=247000 MemberID=51</weight>\n",
      "    <weight id=\"97\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247052\">PDF=247000 MemberID=52</weight>\n",
      "    <weight id=\"98\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247053\">PDF=247000 MemberID=53</weight>\n",
      "    <weight id=\"99\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247054\">PDF=247000 MemberID=54</weight>\n",
      "    <weight id=\"100\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247055\">PDF=247000 MemberID=55</weight>\n",
      "    <weight id=\"101\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247056\">PDF=247000 MemberID=56</weight>\n",
      "    <weight id=\"102\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247057\">PDF=247000 MemberID=57</weight>\n",
      "    <weight id=\"103\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247058\">PDF=247000 MemberID=58</weight>\n",
      "    <weight id=\"104\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247059\">PDF=247000 MemberID=59</weight>\n",
      "    <weight id=\"105\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247060\">PDF=247000 MemberID=60</weight>\n",
      "    <weight id=\"106\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247061\">PDF=247000 MemberID=61</weight>\n",
      "    <weight id=\"107\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247062\">PDF=247000 MemberID=62</weight>\n",
      "    <weight id=\"108\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247063\">PDF=247000 MemberID=63</weight>\n",
      "    <weight id=\"109\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247064\">PDF=247000 MemberID=64</weight>\n",
      "    <weight id=\"110\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247065\">PDF=247000 MemberID=65</weight>\n",
      "    <weight id=\"111\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247066\">PDF=247000 MemberID=66</weight>\n",
      "    <weight id=\"112\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247067\">PDF=247000 MemberID=67</weight>\n",
      "    <weight id=\"113\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247068\">PDF=247000 MemberID=68</weight>\n",
      "    <weight id=\"114\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247069\">PDF=247000 MemberID=69</weight>\n",
      "    <weight id=\"115\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247070\">PDF=247000 MemberID=70</weight>\n",
      "    <weight id=\"116\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247071\">PDF=247000 MemberID=71</weight>\n",
      "    <weight id=\"117\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247072\">PDF=247000 MemberID=72</weight>\n",
      "    <weight id=\"118\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247073\">PDF=247000 MemberID=73</weight>\n",
      "    <weight id=\"119\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247074\">PDF=247000 MemberID=74</weight>\n",
      "    <weight id=\"120\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247075\">PDF=247000 MemberID=75</weight>\n",
      "    <weight id=\"121\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247076\">PDF=247000 MemberID=76</weight>\n",
      "    <weight id=\"122\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247077\">PDF=247000 MemberID=77</weight>\n",
      "    <weight id=\"123\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247078\">PDF=247000 MemberID=78</weight>\n",
      "    <weight id=\"124\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247079\">PDF=247000 MemberID=79</weight>\n",
      "    <weight id=\"125\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247080\">PDF=247000 MemberID=80</weight>\n",
      "    <weight id=\"126\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247081\">PDF=247000 MemberID=81</weight>\n",
      "    <weight id=\"127\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247082\">PDF=247000 MemberID=82</weight>\n",
      "    <weight id=\"128\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247083\">PDF=247000 MemberID=83</weight>\n",
      "    <weight id=\"129\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247084\">PDF=247000 MemberID=84</weight>\n",
      "    <weight id=\"130\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247085\">PDF=247000 MemberID=85</weight>\n",
      "    <weight id=\"131\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247086\">PDF=247000 MemberID=86</weight>\n",
      "    <weight id=\"132\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247087\">PDF=247000 MemberID=87</weight>\n",
      "    <weight id=\"133\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247088\">PDF=247000 MemberID=88</weight>\n",
      "    <weight id=\"134\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247089\">PDF=247000 MemberID=89</weight>\n",
      "    <weight id=\"135\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247090\">PDF=247000 MemberID=90</weight>\n",
      "    <weight id=\"136\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247091\">PDF=247000 MemberID=91</weight>\n",
      "    <weight id=\"137\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247092\">PDF=247000 MemberID=92</weight>\n",
      "    <weight id=\"138\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247093\">PDF=247000 MemberID=93</weight>\n",
      "    <weight id=\"139\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247094\">PDF=247000 MemberID=94</weight>\n",
      "    <weight id=\"140\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247095\">PDF=247000 MemberID=95</weight>\n",
      "    <weight id=\"141\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247096\">PDF=247000 MemberID=96</weight>\n",
      "    <weight id=\"142\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247097\">PDF=247000 MemberID=97</weight>\n",
      "    <weight id=\"143\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247098\">PDF=247000 MemberID=98</weight>\n",
      "    <weight id=\"144\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247099\">PDF=247000 MemberID=99</weight>\n",
      "    <weight id=\"145\" MUR=\"1.0\" MUF=\"1.0\" PDF=\"247100\">PDF=247000 MemberID=100</weight>\n",
      "  </weightgroup>\n",
      "</initrwgt>\n",
      "</init>\n",
      "<event>\n",
      "  4      1  1.6782100000e+03  8.9368380000e+01  7.5467710000e-03  1.3043880000e-01\n",
      "    2  -1   0   0 501   0  0.00000000e+00  0.00000000e+00  7.43180309e+01  7.43180309e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   -2  -1   0   0   0 501 -0.00000000e+00 -0.00000000e+00 -2.68666559e+01  2.68666559e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "  -11   1   1   2   0   0 -6.93966660e+00  2.77442206e+01 -1.51471798e+01  3.23625988e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "   11   1   1   2   0   0  6.93966660e+00 -2.77442206e+01  6.25985547e+01  6.88220880e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "<rwgt>\n",
      " <wgt id='1'> 1.4809e+03</wgt>\n",
      " <wgt id='2'> 1.3469e+03</wgt>\n",
      " <wgt id='3'> 1.3469e+03</wgt>\n",
      " <wgt id='4'> 1.1295e+03</wgt>\n",
      " <wgt id='5'> 1.4809e+03</wgt>\n",
      " <wgt id='6'> 1.6782e+03</wgt>\n",
      " <wgt id='7'> 1.5528e+03</wgt>\n",
      " <wgt id='8'> 1.5528e+03</wgt>\n",
      " <wgt id='9'> 1.3469e+03</wgt>\n",
      " <wgt id='10'> 1.6782e+03</wgt>\n",
      " <wgt id='11'> 1.8614e+03</wgt>\n",
      " <wgt id='12'> 1.7451e+03</wgt>\n",
      " <wgt id='13'> 1.7451e+03</wgt>\n",
      " <wgt id='14'> 1.5528e+03</wgt>\n",
      " <wgt id='15'> 1.8614e+03</wgt>\n",
      " <wgt id='16'> 1.4809e+03</wgt>\n",
      " <wgt id='17'> 1.3469e+03</wgt>\n",
      " <wgt id='18'> 1.3469e+03</wgt>\n",
      " <wgt id='19'> 1.1295e+03</wgt>\n",
      " <wgt id='20'> 1.4809e+03</wgt>\n",
      " <wgt id='21'> 1.5528e+03</wgt>\n",
      " <wgt id='22'> 1.5528e+03</wgt>\n",
      " <wgt id='23'> 1.3469e+03</wgt>\n",
      " <wgt id='24'> 1.6782e+03</wgt>\n",
      " <wgt id='25'> 1.8614e+03</wgt>\n",
      " <wgt id='26'> 1.7451e+03</wgt>\n",
      " <wgt id='27'> 1.7451e+03</wgt>\n",
      " <wgt id='28'> 1.5528e+03</wgt>\n",
      " <wgt id='29'> 1.8614e+03</wgt>\n",
      " <wgt id='30'> 1.4809e+03</wgt>\n",
      " <wgt id='31'> 1.3469e+03</wgt>\n",
      " <wgt id='32'> 1.3469e+03</wgt>\n",
      " <wgt id='33'> 1.1295e+03</wgt>\n",
      " <wgt id='34'> 1.4809e+03</wgt>\n",
      " <wgt id='35'> 1.6782e+03</wgt>\n",
      " <wgt id='36'> 1.5528e+03</wgt>\n",
      " <wgt id='37'> 1.5528e+03</wgt>\n",
      " <wgt id='38'> 1.3469e+03</wgt>\n",
      " <wgt id='39'> 1.6782e+03</wgt>\n",
      " <wgt id='40'> 1.8614e+03</wgt>\n",
      " <wgt id='41'> 1.7451e+03</wgt>\n",
      " <wgt id='42'> 1.7451e+03</wgt>\n",
      " <wgt id='43'> 1.5528e+03</wgt>\n",
      " <wgt id='44'> 1.8614e+03</wgt>\n",
      " <wgt id='45'> 1.6782e+03</wgt>\n",
      " <wgt id='46'> 1.6553e+03</wgt>\n",
      " <wgt id='47'> 1.6314e+03</wgt>\n",
      " <wgt id='48'> 1.6461e+03</wgt>\n",
      " <wgt id='49'> 1.6199e+03</wgt>\n",
      " <wgt id='50'> 1.6843e+03</wgt>\n",
      " <wgt id='51'> 1.6671e+03</wgt>\n",
      " <wgt id='52'> 1.6094e+03</wgt>\n",
      " <wgt id='53'> 1.7531e+03</wgt>\n",
      " <wgt id='54'> 1.6820e+03</wgt>\n",
      " <wgt id='55'> 1.6516e+03</wgt>\n",
      " <wgt id='56'> 1.7195e+03</wgt>\n",
      " <wgt id='57'> 1.6271e+03</wgt>\n",
      " <wgt id='58'> 1.7513e+03</wgt>\n",
      " <wgt id='59'> 1.6538e+03</wgt>\n",
      " <wgt id='60'> 1.6744e+03</wgt>\n",
      " <wgt id='61'> 1.7555e+03</wgt>\n",
      " <wgt id='62'> 1.7074e+03</wgt>\n",
      " <wgt id='63'> 1.6951e+03</wgt>\n",
      " <wgt id='64'> 1.6951e+03</wgt>\n",
      " <wgt id='65'> 1.6721e+03</wgt>\n",
      " <wgt id='66'> 1.7157e+03</wgt>\n",
      " <wgt id='67'> 1.6861e+03</wgt>\n",
      " <wgt id='68'> 1.5079e+03</wgt>\n",
      " <wgt id='69'> 1.5807e+03</wgt>\n",
      " <wgt id='70'> 1.6788e+03</wgt>\n",
      " <wgt id='71'> 1.7110e+03</wgt>\n",
      " <wgt id='72'> 1.7204e+03</wgt>\n",
      " <wgt id='73'> 1.7410e+03</wgt>\n",
      " <wgt id='74'> 1.6152e+03</wgt>\n",
      " <wgt id='75'> 1.5948e+03</wgt>\n",
      " <wgt id='76'> 1.7636e+03</wgt>\n",
      " <wgt id='77'> 1.6926e+03</wgt>\n",
      " <wgt id='78'> 1.7551e+03</wgt>\n",
      " <wgt id='79'> 1.7149e+03</wgt>\n",
      " <wgt id='80'> 1.7482e+03</wgt>\n",
      " <wgt id='81'> 1.6605e+03</wgt>\n",
      " <wgt id='82'> 1.7277e+03</wgt>\n",
      " <wgt id='83'> 1.6063e+03</wgt>\n",
      " <wgt id='84'> 1.6893e+03</wgt>\n",
      " <wgt id='85'> 1.6911e+03</wgt>\n",
      " <wgt id='86'> 1.6646e+03</wgt>\n",
      " <wgt id='87'> 1.7113e+03</wgt>\n",
      " <wgt id='88'> 1.6774e+03</wgt>\n",
      " <wgt id='89'> 1.6866e+03</wgt>\n",
      " <wgt id='90'> 1.7343e+03</wgt>\n",
      " <wgt id='91'> 1.5961e+03</wgt>\n",
      " <wgt id='92'> 1.7314e+03</wgt>\n",
      " <wgt id='93'> 1.6164e+03</wgt>\n",
      " <wgt id='94'> 1.6859e+03</wgt>\n",
      " <wgt id='95'> 1.6587e+03</wgt>\n",
      " <wgt id='96'> 1.6679e+03</wgt>\n",
      " <wgt id='97'> 1.7226e+03</wgt>\n",
      " <wgt id='98'> 1.7133e+03</wgt>\n",
      " <wgt id='99'> 1.7487e+03</wgt>\n",
      " <wgt id='100'> 1.7296e+03</wgt>\n",
      " <wgt id='101'> 1.5388e+03</wgt>\n",
      " <wgt id='102'> 1.5822e+03</wgt>\n",
      " <wgt id='103'> 1.6421e+03</wgt>\n",
      " <wgt id='104'> 1.6703e+03</wgt>\n",
      " <wgt id='105'> 1.7068e+03</wgt>\n",
      " <wgt id='106'> 1.7467e+03</wgt>\n",
      " <wgt id='107'> 1.7071e+03</wgt>\n",
      " <wgt id='108'> 1.7314e+03</wgt>\n",
      " <wgt id='109'> 1.7323e+03</wgt>\n",
      " <wgt id='110'> 1.7469e+03</wgt>\n",
      " <wgt id='111'> 1.7210e+03</wgt>\n",
      " <wgt id='112'> 1.6841e+03</wgt>\n",
      " <wgt id='113'> 1.6620e+03</wgt>\n",
      " <wgt id='114'> 1.7205e+03</wgt>\n",
      " <wgt id='115'> 1.6914e+03</wgt>\n",
      " <wgt id='116'> 1.6384e+03</wgt>\n",
      " <wgt id='117'> 1.6619e+03</wgt>\n",
      " <wgt id='118'> 1.6116e+03</wgt>\n",
      " <wgt id='119'> 1.7048e+03</wgt>\n",
      " <wgt id='120'> 1.6774e+03</wgt>\n",
      " <wgt id='121'> 1.7120e+03</wgt>\n",
      " <wgt id='122'> 1.7478e+03</wgt>\n",
      " <wgt id='123'> 1.7986e+03</wgt>\n",
      " <wgt id='124'> 1.5958e+03</wgt>\n",
      " <wgt id='125'> 1.5796e+03</wgt>\n",
      " <wgt id='126'> 1.6808e+03</wgt>\n",
      " <wgt id='127'> 1.6826e+03</wgt>\n",
      " <wgt id='128'> 1.6533e+03</wgt>\n",
      " <wgt id='129'> 1.6827e+03</wgt>\n",
      " <wgt id='130'> 1.5800e+03</wgt>\n",
      " <wgt id='131'> 1.6695e+03</wgt>\n",
      " <wgt id='132'> 1.5944e+03</wgt>\n",
      " <wgt id='133'> 1.6317e+03</wgt>\n",
      " <wgt id='134'> 1.6429e+03</wgt>\n",
      " <wgt id='135'> 1.6526e+03</wgt>\n",
      " <wgt id='136'> 1.6842e+03</wgt>\n",
      " <wgt id='137'> 1.7144e+03</wgt>\n",
      " <wgt id='138'> 1.6289e+03</wgt>\n",
      " <wgt id='139'> 1.5836e+03</wgt>\n",
      " <wgt id='140'> 1.6705e+03</wgt>\n",
      " <wgt id='141'> 1.6399e+03</wgt>\n",
      " <wgt id='142'> 1.7583e+03</wgt>\n",
      " <wgt id='143'> 1.7348e+03</wgt>\n",
      " <wgt id='144'> 1.7596e+03</wgt>\n",
      " <wgt id='145'> 1.7023e+03</wgt>\n",
      "</rwgt>\n",
      "</event>\n",
      "<event>\n",
      "  5      1  1.6782100000e+03  9.1673260000e+01  7.5467710000e-03  1.2988500000e-01\n",
      "    2  -1   0   0 501   0  0.00000000e+00  0.00000000e+00  2.25287908e+02  2.25287908e+02  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   -2  -1   0   0   0 501 -0.00000000e+00 -0.00000000e+00 -9.32582999e+00  9.32582999e+00  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "   23   2   1   2   0   0  0.00000000e+00  0.00000000e+00  2.15962078e+02  2.34613738e+02  9.16732617e+01  0.0000e+00  0.0000e+00\n",
      "  -13   1   3   3   0   0 -3.21178604e+01 -4.61084611e+00  2.51241424e+01  4.10370490e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "   13   1   3   3   0   0  3.21178604e+01  4.61084611e+00  1.90837936e+02  1.93576689e+02  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "<rwgt>\n",
      " <wgt id='1'> 1.4797e+03</wgt>\n",
      " <wgt id='2'> 1.3760e+03</wgt>\n",
      " <wgt id='3'> 1.3760e+03</wgt>\n",
      " <wgt id='4'> 1.1596e+03</wgt>\n",
      " <wgt id='5'> 1.4797e+03</wgt>\n",
      " <wgt id='6'> 1.6782e+03</wgt>\n",
      " <wgt id='7'> 1.5810e+03</wgt>\n",
      " <wgt id='8'> 1.5810e+03</wgt>\n",
      " <wgt id='9'> 1.3760e+03</wgt>\n",
      " <wgt id='10'> 1.6782e+03</wgt>\n",
      " <wgt id='11'> 1.8629e+03</wgt>\n",
      " <wgt id='12'> 1.7726e+03</wgt>\n",
      " <wgt id='13'> 1.7726e+03</wgt>\n",
      " <wgt id='14'> 1.5810e+03</wgt>\n",
      " <wgt id='15'> 1.8629e+03</wgt>\n",
      " <wgt id='16'> 1.4797e+03</wgt>\n",
      " <wgt id='17'> 1.3760e+03</wgt>\n",
      " <wgt id='18'> 1.3760e+03</wgt>\n",
      " <wgt id='19'> 1.1596e+03</wgt>\n",
      " <wgt id='20'> 1.4797e+03</wgt>\n",
      " <wgt id='21'> 1.5810e+03</wgt>\n",
      " <wgt id='22'> 1.5810e+03</wgt>\n",
      " <wgt id='23'> 1.3760e+03</wgt>\n",
      " <wgt id='24'> 1.6782e+03</wgt>\n",
      " <wgt id='25'> 1.8629e+03</wgt>\n",
      " <wgt id='26'> 1.7726e+03</wgt>\n",
      " <wgt id='27'> 1.7726e+03</wgt>\n",
      " <wgt id='28'> 1.5810e+03</wgt>\n",
      " <wgt id='29'> 1.8629e+03</wgt>\n",
      " <wgt id='30'> 1.4797e+03</wgt>\n",
      " <wgt id='31'> 1.3760e+03</wgt>\n",
      " <wgt id='32'> 1.3760e+03</wgt>\n",
      " <wgt id='33'> 1.1596e+03</wgt>\n",
      " <wgt id='34'> 1.4797e+03</wgt>\n",
      " <wgt id='35'> 1.6782e+03</wgt>\n",
      " <wgt id='36'> 1.5810e+03</wgt>\n",
      " <wgt id='37'> 1.5810e+03</wgt>\n",
      " <wgt id='38'> 1.3760e+03</wgt>\n",
      " <wgt id='39'> 1.6782e+03</wgt>\n",
      " <wgt id='40'> 1.8629e+03</wgt>\n",
      " <wgt id='41'> 1.7726e+03</wgt>\n",
      " <wgt id='42'> 1.7726e+03</wgt>\n",
      " <wgt id='43'> 1.5810e+03</wgt>\n",
      " <wgt id='44'> 1.8629e+03</wgt>\n",
      " <wgt id='45'> 1.6782e+03</wgt>\n",
      " <wgt id='46'> 1.6505e+03</wgt>\n",
      " <wgt id='47'> 1.5966e+03</wgt>\n",
      " <wgt id='48'> 1.6773e+03</wgt>\n",
      " <wgt id='49'> 1.6642e+03</wgt>\n",
      " <wgt id='50'> 1.6235e+03</wgt>\n",
      " <wgt id='51'> 1.7019e+03</wgt>\n",
      " <wgt id='52'> 1.5674e+03</wgt>\n",
      " <wgt id='53'> 1.7236e+03</wgt>\n",
      " <wgt id='54'> 1.6793e+03</wgt>\n",
      " <wgt id='55'> 1.6923e+03</wgt>\n",
      " <wgt id='56'> 1.7391e+03</wgt>\n",
      " <wgt id='57'> 1.6468e+03</wgt>\n",
      " <wgt id='58'> 1.7389e+03</wgt>\n",
      " <wgt id='59'> 1.6653e+03</wgt>\n",
      " <wgt id='60'> 1.7134e+03</wgt>\n",
      " <wgt id='61'> 1.7875e+03</wgt>\n",
      " <wgt id='62'> 1.7295e+03</wgt>\n",
      " <wgt id='63'> 1.6781e+03</wgt>\n",
      " <wgt id='64'> 1.7138e+03</wgt>\n",
      " <wgt id='65'> 1.6368e+03</wgt>\n",
      " <wgt id='66'> 1.7568e+03</wgt>\n",
      " <wgt id='67'> 1.6589e+03</wgt>\n",
      " <wgt id='68'> 1.5084e+03</wgt>\n",
      " <wgt id='69'> 1.6070e+03</wgt>\n",
      " <wgt id='70'> 1.6596e+03</wgt>\n",
      " <wgt id='71'> 1.7595e+03</wgt>\n",
      " <wgt id='72'> 1.7681e+03</wgt>\n",
      " <wgt id='73'> 1.7305e+03</wgt>\n",
      " <wgt id='74'> 1.5992e+03</wgt>\n",
      " <wgt id='75'> 1.5964e+03</wgt>\n",
      " <wgt id='76'> 1.7572e+03</wgt>\n",
      " <wgt id='77'> 1.7148e+03</wgt>\n",
      " <wgt id='78'> 1.7526e+03</wgt>\n",
      " <wgt id='79'> 1.7185e+03</wgt>\n",
      " <wgt id='80'> 1.7406e+03</wgt>\n",
      " <wgt id='81'> 1.6575e+03</wgt>\n",
      " <wgt id='82'> 1.7406e+03</wgt>\n",
      " <wgt id='83'> 1.5879e+03</wgt>\n",
      " <wgt id='84'> 1.6335e+03</wgt>\n",
      " <wgt id='85'> 1.6393e+03</wgt>\n",
      " <wgt id='86'> 1.6191e+03</wgt>\n",
      " <wgt id='87'> 1.6977e+03</wgt>\n",
      " <wgt id='88'> 1.6718e+03</wgt>\n",
      " <wgt id='89'> 1.6228e+03</wgt>\n",
      " <wgt id='90'> 1.6866e+03</wgt>\n",
      " <wgt id='91'> 1.6454e+03</wgt>\n",
      " <wgt id='92'> 1.7344e+03</wgt>\n",
      " <wgt id='93'> 1.6033e+03</wgt>\n",
      " <wgt id='94'> 1.6507e+03</wgt>\n",
      " <wgt id='95'> 1.6655e+03</wgt>\n",
      " <wgt id='96'> 1.6259e+03</wgt>\n",
      " <wgt id='97'> 1.7066e+03</wgt>\n",
      " <wgt id='98'> 1.6883e+03</wgt>\n",
      " <wgt id='99'> 1.7889e+03</wgt>\n",
      " <wgt id='100'> 1.8003e+03</wgt>\n",
      " <wgt id='101'> 1.5503e+03</wgt>\n",
      " <wgt id='102'> 1.5898e+03</wgt>\n",
      " <wgt id='103'> 1.5835e+03</wgt>\n",
      " <wgt id='104'> 1.6640e+03</wgt>\n",
      " <wgt id='105'> 1.6889e+03</wgt>\n",
      " <wgt id='106'> 1.7655e+03</wgt>\n",
      " <wgt id='107'> 1.7128e+03</wgt>\n",
      " <wgt id='108'> 1.7758e+03</wgt>\n",
      " <wgt id='109'> 1.7580e+03</wgt>\n",
      " <wgt id='110'> 1.7507e+03</wgt>\n",
      " <wgt id='111'> 1.7400e+03</wgt>\n",
      " <wgt id='112'> 1.7433e+03</wgt>\n",
      " <wgt id='113'> 1.7098e+03</wgt>\n",
      " <wgt id='114'> 1.7463e+03</wgt>\n",
      " <wgt id='115'> 1.6988e+03</wgt>\n",
      " <wgt id='116'> 1.6599e+03</wgt>\n",
      " <wgt id='117'> 1.6442e+03</wgt>\n",
      " <wgt id='118'> 1.6225e+03</wgt>\n",
      " <wgt id='119'> 1.6879e+03</wgt>\n",
      " <wgt id='120'> 1.6673e+03</wgt>\n",
      " <wgt id='121'> 1.7409e+03</wgt>\n",
      " <wgt id='122'> 1.7359e+03</wgt>\n",
      " <wgt id='123'> 1.7596e+03</wgt>\n",
      " <wgt id='124'> 1.5767e+03</wgt>\n",
      " <wgt id='125'> 1.5483e+03</wgt>\n",
      " <wgt id='126'> 1.6551e+03</wgt>\n",
      " <wgt id='127'> 1.6488e+03</wgt>\n",
      " <wgt id='128'> 1.6914e+03</wgt>\n",
      " <wgt id='129'> 1.6361e+03</wgt>\n",
      " <wgt id='130'> 1.6231e+03</wgt>\n",
      " <wgt id='131'> 1.6349e+03</wgt>\n",
      " <wgt id='132'> 1.5944e+03</wgt>\n",
      " <wgt id='133'> 1.6045e+03</wgt>\n",
      " <wgt id='134'> 1.6687e+03</wgt>\n",
      " <wgt id='135'> 1.6728e+03</wgt>\n",
      " <wgt id='136'> 1.6807e+03</wgt>\n",
      " <wgt id='137'> 1.7129e+03</wgt>\n",
      " <wgt id='138'> 1.6707e+03</wgt>\n",
      " <wgt id='139'> 1.5443e+03</wgt>\n",
      " <wgt id='140'> 1.6675e+03</wgt>\n",
      " <wgt id='141'> 1.6382e+03</wgt>\n",
      " <wgt id='142'> 1.8170e+03</wgt>\n",
      " <wgt id='143'> 1.7309e+03</wgt>\n",
      " <wgt id='144'> 1.7507e+03</wgt>\n",
      " <wgt id='145'> 1.6814e+03</wgt>\n",
      "</rwgt>\n",
      "</event>\n",
      "<event>\n",
      "  5      1  1.6782100000e+03  9.3375730000e+01  7.5467710000e-03  1.2948770000e-01\n",
      "   -2  -1   0   0   0 501 -0.00000000e+00  0.00000000e+00  6.38598869e+01  6.38598869e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "    2  -1   0   0 501   0  0.00000000e+00 -0.00000000e+00 -3.41334250e+01  3.41334250e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   23   2   1   2   0   0  0.00000000e+00  0.00000000e+00  2.97264619e+01  9.79933119e+01  9.33757283e+01  0.0000e+00  0.0000e+00\n",
      "  -11   1   3   3   0   0  1.78799500e+01 -4.32808425e+00 -3.01695174e+01  3.53358841e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   11   1   3   3   0   0 -1.78799500e+01  4.32808425e+00  5.98959793e+01  6.26574278e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "<rwgt>\n",
      " <wgt id='1'> 1.4884e+03</wgt>\n",
      " <wgt id='2'> 1.2115e+03</wgt>\n",
      " <wgt id='3'> 1.2115e+03</wgt>\n",
      " <wgt id='4'> 9.9299e+02</wgt>\n",
      " <wgt id='5'> 1.4884e+03</wgt>\n",
      " <wgt id='6'> 1.6782e+03</wgt>\n",
      " <wgt id='7'> 1.4198e+03</wgt>\n",
      " <wgt id='8'> 1.4198e+03</wgt>\n",
      " <wgt id='9'> 1.2115e+03</wgt>\n",
      " <wgt id='10'> 1.6782e+03</wgt>\n",
      " <wgt id='11'> 1.8532e+03</wgt>\n",
      " <wgt id='12'> 1.6146e+03</wgt>\n",
      " <wgt id='13'> 1.6146e+03</wgt>\n",
      " <wgt id='14'> 1.4198e+03</wgt>\n",
      " <wgt id='15'> 1.8532e+03</wgt>\n",
      " <wgt id='16'> 1.4884e+03</wgt>\n",
      " <wgt id='17'> 1.2115e+03</wgt>\n",
      " <wgt id='18'> 1.2115e+03</wgt>\n",
      " <wgt id='19'> 9.9299e+02</wgt>\n",
      " <wgt id='20'> 1.4884e+03</wgt>\n",
      " <wgt id='21'> 1.4198e+03</wgt>\n",
      " <wgt id='22'> 1.4198e+03</wgt>\n",
      " <wgt id='23'> 1.2115e+03</wgt>\n",
      " <wgt id='24'> 1.6782e+03</wgt>\n",
      " <wgt id='25'> 1.8532e+03</wgt>\n",
      " <wgt id='26'> 1.6146e+03</wgt>\n",
      " <wgt id='27'> 1.6146e+03</wgt>\n",
      " <wgt id='28'> 1.4198e+03</wgt>\n",
      " <wgt id='29'> 1.8532e+03</wgt>\n",
      " <wgt id='30'> 1.4884e+03</wgt>\n",
      " <wgt id='31'> 1.2115e+03</wgt>\n",
      " <wgt id='32'> 1.2115e+03</wgt>\n",
      " <wgt id='33'> 9.9299e+02</wgt>\n",
      " <wgt id='34'> 1.4884e+03</wgt>\n",
      " <wgt id='35'> 1.6782e+03</wgt>\n",
      " <wgt id='36'> 1.4198e+03</wgt>\n",
      " <wgt id='37'> 1.4198e+03</wgt>\n",
      " <wgt id='38'> 1.2115e+03</wgt>\n",
      " <wgt id='39'> 1.6782e+03</wgt>\n",
      " <wgt id='40'> 1.8532e+03</wgt>\n",
      " <wgt id='41'> 1.6146e+03</wgt>\n",
      " <wgt id='42'> 1.6146e+03</wgt>\n",
      " <wgt id='43'> 1.4198e+03</wgt>\n",
      " <wgt id='44'> 1.8532e+03</wgt>\n",
      " <wgt id='45'> 1.6782e+03</wgt>\n",
      " <wgt id='46'> 1.6756e+03</wgt>\n",
      " <wgt id='47'> 1.6362e+03</wgt>\n",
      " <wgt id='48'> 1.6304e+03</wgt>\n",
      " <wgt id='49'> 1.6026e+03</wgt>\n",
      " <wgt id='50'> 1.6587e+03</wgt>\n",
      " <wgt id='51'> 1.6492e+03</wgt>\n",
      " <wgt id='52'> 1.6155e+03</wgt>\n",
      " <wgt id='53'> 1.7550e+03</wgt>\n",
      " <wgt id='54'> 1.6699e+03</wgt>\n",
      " <wgt id='55'> 1.6126e+03</wgt>\n",
      " <wgt id='56'> 1.7880e+03</wgt>\n",
      " <wgt id='57'> 1.6070e+03</wgt>\n",
      " <wgt id='58'> 1.7433e+03</wgt>\n",
      " <wgt id='59'> 1.6410e+03</wgt>\n",
      " <wgt id='60'> 1.6602e+03</wgt>\n",
      " <wgt id='61'> 1.7573e+03</wgt>\n",
      " <wgt id='62'> 1.7027e+03</wgt>\n",
      " <wgt id='63'> 1.7252e+03</wgt>\n",
      " <wgt id='64'> 1.6705e+03</wgt>\n",
      " <wgt id='65'> 1.6769e+03</wgt>\n",
      " <wgt id='66'> 1.7027e+03</wgt>\n",
      " <wgt id='67'> 1.7205e+03</wgt>\n",
      " <wgt id='68'> 1.4845e+03</wgt>\n",
      " <wgt id='69'> 1.5524e+03</wgt>\n",
      " <wgt id='70'> 1.6814e+03</wgt>\n",
      " <wgt id='71'> 1.6852e+03</wgt>\n",
      " <wgt id='72'> 1.6972e+03</wgt>\n",
      " <wgt id='73'> 1.7453e+03</wgt>\n",
      " <wgt id='74'> 1.5968e+03</wgt>\n",
      " <wgt id='75'> 1.6142e+03</wgt>\n",
      " <wgt id='76'> 1.7523e+03</wgt>\n",
      " <wgt id='77'> 1.6797e+03</wgt>\n",
      " <wgt id='78'> 1.7518e+03</wgt>\n",
      " <wgt id='79'> 1.7070e+03</wgt>\n",
      " <wgt id='80'> 1.7658e+03</wgt>\n",
      " <wgt id='81'> 1.6692e+03</wgt>\n",
      " <wgt id='82'> 1.7026e+03</wgt>\n",
      " <wgt id='83'> 1.6141e+03</wgt>\n",
      " <wgt id='84'> 1.7054e+03</wgt>\n",
      " <wgt id='85'> 1.7141e+03</wgt>\n",
      " <wgt id='86'> 1.6910e+03</wgt>\n",
      " <wgt id='87'> 1.7081e+03</wgt>\n",
      " <wgt id='88'> 1.6751e+03</wgt>\n",
      " <wgt id='89'> 1.7556e+03</wgt>\n",
      " <wgt id='90'> 1.8099e+03</wgt>\n",
      " <wgt id='91'> 1.5709e+03</wgt>\n",
      " <wgt id='92'> 1.7439e+03</wgt>\n",
      " <wgt id='93'> 1.6318e+03</wgt>\n",
      " <wgt id='94'> 1.7276e+03</wgt>\n",
      " <wgt id='95'> 1.6557e+03</wgt>\n",
      " <wgt id='96'> 1.6843e+03</wgt>\n",
      " <wgt id='97'> 1.7396e+03</wgt>\n",
      " <wgt id='98'> 1.7271e+03</wgt>\n",
      " <wgt id='99'> 1.7354e+03</wgt>\n",
      " <wgt id='100'> 1.7143e+03</wgt>\n",
      " <wgt id='101'> 1.5204e+03</wgt>\n",
      " <wgt id='102'> 1.5725e+03</wgt>\n",
      " <wgt id='103'> 1.6756e+03</wgt>\n",
      " <wgt id='104'> 1.6666e+03</wgt>\n",
      " <wgt id='105'> 1.7046e+03</wgt>\n",
      " <wgt id='106'> 1.7546e+03</wgt>\n",
      " <wgt id='107'> 1.6875e+03</wgt>\n",
      " <wgt id='108'> 1.7354e+03</wgt>\n",
      " <wgt id='109'> 1.7173e+03</wgt>\n",
      " <wgt id='110'> 1.7495e+03</wgt>\n",
      " <wgt id='111'> 1.7122e+03</wgt>\n",
      " <wgt id='112'> 1.6610e+03</wgt>\n",
      " <wgt id='113'> 1.6468e+03</wgt>\n",
      " <wgt id='114'> 1.6984e+03</wgt>\n",
      " <wgt id='115'> 1.6933e+03</wgt>\n",
      " <wgt id='116'> 1.6180e+03</wgt>\n",
      " <wgt id='117'> 1.6804e+03</wgt>\n",
      " <wgt id='118'> 1.5969e+03</wgt>\n",
      " <wgt id='119'> 1.7142e+03</wgt>\n",
      " <wgt id='120'> 1.6759e+03</wgt>\n",
      " <wgt id='121'> 1.6946e+03</wgt>\n",
      " <wgt id='122'> 1.7439e+03</wgt>\n",
      " <wgt id='123'> 1.8105e+03</wgt>\n",
      " <wgt id='124'> 1.5911e+03</wgt>\n",
      " <wgt id='125'> 1.5705e+03</wgt>\n",
      " <wgt id='126'> 1.6981e+03</wgt>\n",
      " <wgt id='127'> 1.7072e+03</wgt>\n",
      " <wgt id='128'> 1.6440e+03</wgt>\n",
      " <wgt id='129'> 1.6651e+03</wgt>\n",
      " <wgt id='130'> 1.5546e+03</wgt>\n",
      " <wgt id='131'> 1.6788e+03</wgt>\n",
      " <wgt id='132'> 1.5968e+03</wgt>\n",
      " <wgt id='133'> 1.6781e+03</wgt>\n",
      " <wgt id='134'> 1.6160e+03</wgt>\n",
      " <wgt id='135'> 1.6230e+03</wgt>\n",
      " <wgt id='136'> 1.6649e+03</wgt>\n",
      " <wgt id='137'> 1.7232e+03</wgt>\n",
      " <wgt id='138'> 1.6322e+03</wgt>\n",
      " <wgt id='139'> 1.6590e+03</wgt>\n",
      " <wgt id='140'> 1.6470e+03</wgt>\n",
      " <wgt id='141'> 1.6213e+03</wgt>\n",
      " <wgt id='142'> 1.7461e+03</wgt>\n",
      " <wgt id='143'> 1.7365e+03</wgt>\n",
      " <wgt id='144'> 1.7717e+03</wgt>\n",
      " <wgt id='145'> 1.7176e+03</wgt>\n",
      "</rwgt>\n",
      "</event>\n",
      "<event>\n",
      "  4      1  1.6782100000e+03  9.1638730000e+01  7.5467710000e-03  1.2989310000e-01\n",
      "    1  -1   0   0 501   0  0.00000000e+00  0.00000000e+00  4.46125101e+01  4.46125101e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   -1  -1   0   0   0 501 -0.00000000e+00 -0.00000000e+00 -4.70588647e+01  4.70588647e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "  -11   1   1   2   0   0  1.28771854e+01  4.38209466e+01 -4.87369922e+00  4.59331058e+01  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "   11   1   1   2   0   0 -1.28771854e+01 -4.38209466e+01  2.42734462e+00  4.57382691e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "<rwgt>\n",
      " <wgt id='1'> 1.4856e+03</wgt>\n",
      " <wgt id='2'> 1.4847e+03</wgt>\n",
      " <wgt id='3'> 1.4847e+03</wgt>\n",
      " <wgt id='4'> 1.2781e+03</wgt>\n",
      " <wgt id='5'> 1.4856e+03</wgt>\n",
      " <wgt id='6'> 1.6782e+03</wgt>\n",
      " <wgt id='7'> 1.6774e+03</wgt>\n",
      " <wgt id='8'> 1.6774e+03</wgt>\n",
      " <wgt id='9'> 1.4847e+03</wgt>\n",
      " <wgt id='10'> 1.6782e+03</wgt>\n",
      " <wgt id='11'> 1.8560e+03</wgt>\n",
      " <wgt id='12'> 1.8552e+03</wgt>\n",
      " <wgt id='13'> 1.8552e+03</wgt>\n",
      " <wgt id='14'> 1.6774e+03</wgt>\n",
      " <wgt id='15'> 1.8560e+03</wgt>\n",
      " <wgt id='16'> 1.4856e+03</wgt>\n",
      " <wgt id='17'> 1.4847e+03</wgt>\n",
      " <wgt id='18'> 1.4847e+03</wgt>\n",
      " <wgt id='19'> 1.2781e+03</wgt>\n",
      " <wgt id='20'> 1.4856e+03</wgt>\n",
      " <wgt id='21'> 1.6774e+03</wgt>\n",
      " <wgt id='22'> 1.6774e+03</wgt>\n",
      " <wgt id='23'> 1.4847e+03</wgt>\n",
      " <wgt id='24'> 1.6782e+03</wgt>\n",
      " <wgt id='25'> 1.8560e+03</wgt>\n",
      " <wgt id='26'> 1.8552e+03</wgt>\n",
      " <wgt id='27'> 1.8552e+03</wgt>\n",
      " <wgt id='28'> 1.6774e+03</wgt>\n",
      " <wgt id='29'> 1.8560e+03</wgt>\n",
      " <wgt id='30'> 1.4856e+03</wgt>\n",
      " <wgt id='31'> 1.4847e+03</wgt>\n",
      " <wgt id='32'> 1.4847e+03</wgt>\n",
      " <wgt id='33'> 1.2781e+03</wgt>\n",
      " <wgt id='34'> 1.4856e+03</wgt>\n",
      " <wgt id='35'> 1.6782e+03</wgt>\n",
      " <wgt id='36'> 1.6774e+03</wgt>\n",
      " <wgt id='37'> 1.6774e+03</wgt>\n",
      " <wgt id='38'> 1.4847e+03</wgt>\n",
      " <wgt id='39'> 1.6782e+03</wgt>\n",
      " <wgt id='40'> 1.8560e+03</wgt>\n",
      " <wgt id='41'> 1.8552e+03</wgt>\n",
      " <wgt id='42'> 1.8552e+03</wgt>\n",
      " <wgt id='43'> 1.6774e+03</wgt>\n",
      " <wgt id='44'> 1.8560e+03</wgt>\n",
      " <wgt id='45'> 1.6782e+03</wgt>\n",
      " <wgt id='46'> 1.6429e+03</wgt>\n",
      " <wgt id='47'> 1.6468e+03</wgt>\n",
      " <wgt id='48'> 1.6162e+03</wgt>\n",
      " <wgt id='49'> 1.6130e+03</wgt>\n",
      " <wgt id='50'> 1.6856e+03</wgt>\n",
      " <wgt id='51'> 1.6179e+03</wgt>\n",
      " <wgt id='52'> 1.6689e+03</wgt>\n",
      " <wgt id='53'> 1.7848e+03</wgt>\n",
      " <wgt id='54'> 1.6704e+03</wgt>\n",
      " <wgt id='55'> 1.6646e+03</wgt>\n",
      " <wgt id='56'> 1.8328e+03</wgt>\n",
      " <wgt id='57'> 1.5622e+03</wgt>\n",
      " <wgt id='58'> 1.7942e+03</wgt>\n",
      " <wgt id='59'> 1.5703e+03</wgt>\n",
      " <wgt id='60'> 1.6517e+03</wgt>\n",
      " <wgt id='61'> 1.7496e+03</wgt>\n",
      " <wgt id='62'> 1.6967e+03</wgt>\n",
      " <wgt id='63'> 1.7457e+03</wgt>\n",
      " <wgt id='64'> 1.7115e+03</wgt>\n",
      " <wgt id='65'> 1.7416e+03</wgt>\n",
      " <wgt id='66'> 1.6420e+03</wgt>\n",
      " <wgt id='67'> 1.7201e+03</wgt>\n",
      " <wgt id='68'> 1.5384e+03</wgt>\n",
      " <wgt id='69'> 1.5973e+03</wgt>\n",
      " <wgt id='70'> 1.7240e+03</wgt>\n",
      " <wgt id='71'> 1.6539e+03</wgt>\n",
      " <wgt id='72'> 1.6980e+03</wgt>\n",
      " <wgt id='73'> 1.7167e+03</wgt>\n",
      " <wgt id='74'> 1.6108e+03</wgt>\n",
      " <wgt id='75'> 1.6938e+03</wgt>\n",
      " <wgt id='76'> 1.7726e+03</wgt>\n",
      " <wgt id='77'> 1.7171e+03</wgt>\n",
      " <wgt id='78'> 1.7419e+03</wgt>\n",
      " <wgt id='79'> 1.7235e+03</wgt>\n",
      " <wgt id='80'> 1.7931e+03</wgt>\n",
      " <wgt id='81'> 1.7104e+03</wgt>\n",
      " <wgt id='82'> 1.7226e+03</wgt>\n",
      " <wgt id='83'> 1.5645e+03</wgt>\n",
      " <wgt id='84'> 1.6474e+03</wgt>\n",
      " <wgt id='85'> 1.7408e+03</wgt>\n",
      " <wgt id='86'> 1.7017e+03</wgt>\n",
      " <wgt id='87'> 1.7299e+03</wgt>\n",
      " <wgt id='88'> 1.6401e+03</wgt>\n",
      " <wgt id='89'> 1.6868e+03</wgt>\n",
      " <wgt id='90'> 1.8030e+03</wgt>\n",
      " <wgt id='91'> 1.5299e+03</wgt>\n",
      " <wgt id='92'> 1.6931e+03</wgt>\n",
      " <wgt id='93'> 1.6550e+03</wgt>\n",
      " <wgt id='94'> 1.7837e+03</wgt>\n",
      " <wgt id='95'> 1.6653e+03</wgt>\n",
      " <wgt id='96'> 1.6500e+03</wgt>\n",
      " <wgt id='97'> 1.7732e+03</wgt>\n",
      " <wgt id='98'> 1.6635e+03</wgt>\n",
      " <wgt id='99'> 1.7852e+03</wgt>\n",
      " <wgt id='100'> 1.6770e+03</wgt>\n",
      " <wgt id='101'> 1.5004e+03</wgt>\n",
      " <wgt id='102'> 1.5982e+03</wgt>\n",
      " <wgt id='103'> 1.6108e+03</wgt>\n",
      " <wgt id='104'> 1.5672e+03</wgt>\n",
      " <wgt id='105'> 1.7089e+03</wgt>\n",
      " <wgt id='106'> 1.7370e+03</wgt>\n",
      " <wgt id='107'> 1.7035e+03</wgt>\n",
      " <wgt id='108'> 1.7673e+03</wgt>\n",
      " <wgt id='109'> 1.7298e+03</wgt>\n",
      " <wgt id='110'> 1.7556e+03</wgt>\n",
      " <wgt id='111'> 1.7264e+03</wgt>\n",
      " <wgt id='112'> 1.6840e+03</wgt>\n",
      " <wgt id='113'> 1.6705e+03</wgt>\n",
      " <wgt id='114'> 1.6938e+03</wgt>\n",
      " <wgt id='115'> 1.6859e+03</wgt>\n",
      " <wgt id='116'> 1.5958e+03</wgt>\n",
      " <wgt id='117'> 1.6275e+03</wgt>\n",
      " <wgt id='118'> 1.6093e+03</wgt>\n",
      " <wgt id='119'> 1.7016e+03</wgt>\n",
      " <wgt id='120'> 1.6733e+03</wgt>\n",
      " <wgt id='121'> 1.6883e+03</wgt>\n",
      " <wgt id='122'> 1.6874e+03</wgt>\n",
      " <wgt id='123'> 1.7999e+03</wgt>\n",
      " <wgt id='124'> 1.5355e+03</wgt>\n",
      " <wgt id='125'> 1.5223e+03</wgt>\n",
      " <wgt id='126'> 1.7132e+03</wgt>\n",
      " <wgt id='127'> 1.7204e+03</wgt>\n",
      " <wgt id='128'> 1.6315e+03</wgt>\n",
      " <wgt id='129'> 1.6846e+03</wgt>\n",
      " <wgt id='130'> 1.5239e+03</wgt>\n",
      " <wgt id='131'> 1.7031e+03</wgt>\n",
      " <wgt id='132'> 1.5787e+03</wgt>\n",
      " <wgt id='133'> 1.7015e+03</wgt>\n",
      " <wgt id='134'> 1.6418e+03</wgt>\n",
      " <wgt id='135'> 1.6290e+03</wgt>\n",
      " <wgt id='136'> 1.7311e+03</wgt>\n",
      " <wgt id='137'> 1.6709e+03</wgt>\n",
      " <wgt id='138'> 1.5711e+03</wgt>\n",
      " <wgt id='139'> 1.6856e+03</wgt>\n",
      " <wgt id='140'> 1.6782e+03</wgt>\n",
      " <wgt id='141'> 1.6651e+03</wgt>\n",
      " <wgt id='142'> 1.7950e+03</wgt>\n",
      " <wgt id='143'> 1.7211e+03</wgt>\n",
      " <wgt id='144'> 1.7342e+03</wgt>\n",
      " <wgt id='145'> 1.6896e+03</wgt>\n",
      "</rwgt>\n",
      "</event>\n",
      "<event>\n",
      "  4      1  1.6782100000e+03  9.4862500000e+01  7.5467710000e-03  1.2914860000e-01\n",
      "    2  -1   0   0 501   0  0.00000000e+00  0.00000000e+00  1.87598540e+02  1.87598540e+02  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "   -2  -1   0   0   0 501 -0.00000000e+00 -0.00000000e+00 -1.19922238e+01  1.19922238e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "  -11   1   1   2   0   0 -1.95431649e+01 -3.60928948e+01  3.77875936e+01  5.57900940e+01  0.00000000e+00  0.0000e+00 -1.0000e+00\n",
      "   11   1   1   2   0   0  1.95431649e+01  3.60928948e+01  1.37818723e+02  1.43800670e+02  0.00000000e+00  0.0000e+00  1.0000e+00\n",
      "<rwgt>\n",
      " <wgt id='1'> 1.4853e+03</wgt>\n",
      " <wgt id='2'> 1.4434e+03</wgt>\n",
      " <wgt id='3'> 1.4434e+03</wgt>\n",
      " <wgt id='4'> 1.2349e+03</wgt>\n",
      " <wgt id='5'> 1.4853e+03</wgt>\n",
      " <wgt id='6'> 1.6782e+03</wgt>\n",
      " <wgt id='7'> 1.6391e+03</wgt>\n",
      " <wgt id='8'> 1.6391e+03</wgt>\n",
      " <wgt id='9'> 1.4434e+03</wgt>\n",
      " <wgt id='10'> 1.6782e+03</wgt>\n",
      " <wgt id='11'> 1.8570e+03</wgt>\n",
      " <wgt id='12'> 1.8209e+03</wgt>\n",
      " <wgt id='13'> 1.8209e+03</wgt>\n",
      " <wgt id='14'> 1.6391e+03</wgt>\n",
      " <wgt id='15'> 1.8570e+03</wgt>\n",
      " <wgt id='16'> 1.4853e+03</wgt>\n",
      " <wgt id='17'> 1.4434e+03</wgt>\n",
      " <wgt id='18'> 1.4434e+03</wgt>\n",
      " <wgt id='19'> 1.2349e+03</wgt>\n",
      " <wgt id='20'> 1.4853e+03</wgt>\n",
      " <wgt id='21'> 1.6391e+03</wgt>\n",
      " <wgt id='22'> 1.6391e+03</wgt>\n",
      " <wgt id='23'> 1.4434e+03</wgt>\n",
      " <wgt id='24'> 1.6782e+03</wgt>\n",
      " <wgt id='25'> 1.8570e+03</wgt>\n",
      " <wgt id='26'> 1.8209e+03</wgt>\n",
      " <wgt id='27'> 1.8209e+03</wgt>\n",
      " <wgt id='28'> 1.6391e+03</wgt>\n",
      " <wgt id='29'> 1.8570e+03</wgt>\n",
      " <wgt id='30'> 1.4853e+03</wgt>\n",
      " <wgt id='31'> 1.4434e+03</wgt>\n",
      " <wgt id='32'> 1.4434e+03</wgt>\n",
      " <wgt id='33'> 1.2349e+03</wgt>\n",
      " <wgt id='34'> 1.4853e+03</wgt>\n",
      " <wgt id='35'> 1.6782e+03</wgt>\n",
      " <wgt id='36'> 1.6391e+03</wgt>\n",
      " <wgt id='37'> 1.6391e+03</wgt>\n",
      " <wgt id='38'> 1.4434e+03</wgt>\n",
      " <wgt id='39'> 1.6782e+03</wgt>\n",
      " <wgt id='40'> 1.8570e+03</wgt>\n",
      " <wgt id='41'> 1.8209e+03</wgt>\n",
      " <wgt id='42'> 1.8209e+03</wgt>\n",
      " <wgt id='43'> 1.6391e+03</wgt>\n",
      " <wgt id='44'> 1.8570e+03</wgt>\n",
      " <wgt id='45'> 1.6782e+03</wgt>\n",
      " <wgt id='46'> 1.6499e+03</wgt>\n",
      " <wgt id='47'> 1.6066e+03</wgt>\n",
      " <wgt id='48'> 1.6708e+03</wgt>\n",
      " <wgt id='49'> 1.6574e+03</wgt>\n",
      " <wgt id='50'> 1.6404e+03</wgt>\n",
      " <wgt id='51'> 1.6968e+03</wgt>\n",
      " <wgt id='52'> 1.5766e+03</wgt>\n",
      " <wgt id='53'> 1.7301e+03</wgt>\n",
      " <wgt id='54'> 1.6813e+03</wgt>\n",
      " <wgt id='55'> 1.6864e+03</wgt>\n",
      " <wgt id='56'> 1.7276e+03</wgt>\n",
      " <wgt id='57'> 1.6457e+03</wgt>\n",
      " <wgt id='58'> 1.7422e+03</wgt>\n",
      " <wgt id='59'> 1.6643e+03</wgt>\n",
      " <wgt id='60'> 1.7046e+03</wgt>\n",
      " <wgt id='61'> 1.7812e+03</wgt>\n",
      " <wgt id='62'> 1.7269e+03</wgt>\n",
      " <wgt id='63'> 1.6793e+03</wgt>\n",
      " <wgt id='64'> 1.7117e+03</wgt>\n",
      " <wgt id='65'> 1.6415e+03</wgt>\n",
      " <wgt id='66'> 1.7508e+03</wgt>\n",
      " <wgt id='67'> 1.6613e+03</wgt>\n",
      " <wgt id='68'> 1.5098e+03</wgt>\n",
      " <wgt id='69'> 1.6033e+03</wgt>\n",
      " <wgt id='70'> 1.6635e+03</wgt>\n",
      " <wgt id='71'> 1.7519e+03</wgt>\n",
      " <wgt id='72'> 1.7578e+03</wgt>\n",
      " <wgt id='73'> 1.7314e+03</wgt>\n",
      " <wgt id='74'> 1.6045e+03</wgt>\n",
      " <wgt id='75'> 1.5908e+03</wgt>\n",
      " <wgt id='76'> 1.7609e+03</wgt>\n",
      " <wgt id='77'> 1.7125e+03</wgt>\n",
      " <wgt id='78'> 1.7523e+03</wgt>\n",
      " <wgt id='79'> 1.7188e+03</wgt>\n",
      " <wgt id='80'> 1.7395e+03</wgt>\n",
      " <wgt id='81'> 1.6576e+03</wgt>\n",
      " <wgt id='82'> 1.7401e+03</wgt>\n",
      " <wgt id='83'> 1.5894e+03</wgt>\n",
      " <wgt id='84'> 1.6494e+03</wgt>\n",
      " <wgt id='85'> 1.6464e+03</wgt>\n",
      " <wgt id='86'> 1.6246e+03</wgt>\n",
      " <wgt id='87'> 1.7022e+03</wgt>\n",
      " <wgt id='88'> 1.6742e+03</wgt>\n",
      " <wgt id='89'> 1.6311e+03</wgt>\n",
      " <wgt id='90'> 1.6913e+03</wgt>\n",
      " <wgt id='91'> 1.6360e+03</wgt>\n",
      " <wgt id='92'> 1.7300e+03</wgt>\n",
      " <wgt id='93'> 1.6052e+03</wgt>\n",
      " <wgt id='94'> 1.6524e+03</wgt>\n",
      " <wgt id='95'> 1.6642e+03</wgt>\n",
      " <wgt id='96'> 1.6335e+03</wgt>\n",
      " <wgt id='97'> 1.7074e+03</wgt>\n",
      " <wgt id='98'> 1.6903e+03</wgt>\n",
      " <wgt id='99'> 1.7843e+03</wgt>\n",
      " <wgt id='100'> 1.7860e+03</wgt>\n",
      " <wgt id='101'> 1.5502e+03</wgt>\n",
      " <wgt id='102'> 1.5893e+03</wgt>\n",
      " <wgt id='103'> 1.5945e+03</wgt>\n",
      " <wgt id='104'> 1.6666e+03</wgt>\n",
      " <wgt id='105'> 1.6919e+03</wgt>\n",
      " <wgt id='106'> 1.7579e+03</wgt>\n",
      " <wgt id='107'> 1.7130e+03</wgt>\n",
      " <wgt id='108'> 1.7634e+03</wgt>\n",
      " <wgt id='109'> 1.7523e+03</wgt>\n",
      " <wgt id='110'> 1.7510e+03</wgt>\n",
      " <wgt id='111'> 1.7362e+03</wgt>\n",
      " <wgt id='112'> 1.7321e+03</wgt>\n",
      " <wgt id='113'> 1.6999e+03</wgt>\n",
      " <wgt id='114'> 1.7430e+03</wgt>\n",
      " <wgt id='115'> 1.6965e+03</wgt>\n",
      " <wgt id='116'> 1.6589e+03</wgt>\n",
      " <wgt id='117'> 1.6458e+03</wgt>\n",
      " <wgt id='118'> 1.6210e+03</wgt>\n",
      " <wgt id='119'> 1.6908e+03</wgt>\n",
      " <wgt id='120'> 1.6696e+03</wgt>\n",
      " <wgt id='121'> 1.7365e+03</wgt>\n",
      " <wgt id='122'> 1.7383e+03</wgt>\n",
      " <wgt id='123'> 1.7671e+03</wgt>\n",
      " <wgt id='124'> 1.5845e+03</wgt>\n",
      " <wgt id='125'> 1.5555e+03</wgt>\n",
      " <wgt id='126'> 1.6596e+03</wgt>\n",
      " <wgt id='127'> 1.6536e+03</wgt>\n",
      " <wgt id='128'> 1.6854e+03</wgt>\n",
      " <wgt id='129'> 1.6476e+03</wgt>\n",
      " <wgt id='130'> 1.6182e+03</wgt>\n",
      " <wgt id='131'> 1.6414e+03</wgt>\n",
      " <wgt id='132'> 1.5937e+03</wgt>\n",
      " <wgt id='133'> 1.6078e+03</wgt>\n",
      " <wgt id='134'> 1.6649e+03</wgt>\n",
      " <wgt id='135'> 1.6722e+03</wgt>\n",
      " <wgt id='136'> 1.6840e+03</wgt>\n",
      " <wgt id='137'> 1.7125e+03</wgt>\n",
      " <wgt id='138'> 1.6588e+03</wgt>\n",
      " <wgt id='139'> 1.5446e+03</wgt>\n",
      " <wgt id='140'> 1.6702e+03</wgt>\n",
      " <wgt id='141'> 1.6419e+03</wgt>\n",
      " <wgt id='142'> 1.8054e+03</wgt>\n",
      " <wgt id='143'> 1.7315e+03</wgt>\n",
      " <wgt id='144'> 1.7515e+03</wgt>\n",
      " <wgt id='145'> 1.6859e+03</wgt>\n",
      "</rwgt>\n",
      "</event>\n",
      "</LesHouchesEvents>\n"
     ]
    }
   ],
   "source": [
    "lhe_file = pylhe.read_lhe_file(lhe_data)\n",
    "print(\n",
    "    pylhe.write_lhe_string(\n",
    "        lhe_file.init, [lhe_file.events.__next__() for i in range(5)]\n",
    "    )\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b2a72c22-5caa-4505-8c6e-708f9114900f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAGzCAYAAAAmH71NAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAxQElEQVR4nO3de3xTdZ7/8XdIL7QyLUK1UOkFVkViuUhbXVEQhClbLo4yO6KzIiq4skQRO15AULGiqDiIMwQUZhEf4yqMKzK6smLBC4y42haqSBmB2WJRCywXG6HQlPT8/uDXjKH3Nu1JTl7PxyOPBznfk9NPvtbmne/3e86xGYZhCAAAwII6mV0AAABAeyHoAAAAyyLoAAAAyyLoAAAAyyLoAAAAyyLoAAAAyyLoAAAAyyLoAAAAyyLoAAAAyyLoAAAAy4owuwCz1dTU6Pvvv9fPfvYz2Ww2s8sBAADNYBiGfvzxRyUlJalTp4bHbUI+6Hz99deaOHGi3/PXX39d119/fbNe//333ys5ObmdqgMAAO1p//796tWrV4PtNivd1PP48eNKS0vTN998o3POOadZr6moqFDXrl21f/9+xcXFtXOFAAAgENxut5KTk/XDDz8oPj6+wf1CfkTnp95++22NHDmy2SFHkm+6Ki4ujqADAECIaWrZiemLkTdv3qzx48crKSlJNptN69atq7PP0qVL1bt3b3Xu3FkZGRnasmVLvcf605/+5DeNBQAAwpvpQefEiRMaOHCglixZUm/7mjVrNHPmTM2ZM0fbt2/X0KFDlZOTo7KyMr/93G63PvnkE40ZM6bRn1dVVSW32+33AAAA1hRUa3RsNpveeustv4XEV1xxhQYPHqxly5b5tvXr10/XX3+9FixY4Nv2xz/+URs2bNCrr77a6M+YN2+eHn/88TrbKyoqmLoCACBEuN1uxcfHN/n5bfqITmM8Ho+KioqUnZ3ttz07O1tbt27129bcaavZs2eroqLC99i/f39AawYAAMEjqBcjHz58WF6vV4mJiX7bExMTdeDAAd/ziooKff7553rzzTebPGZ0dLSio6MDXisAAAg+QR10ap29otowDL9t8fHxOnjwYEeXBQAAglxQT10lJCTIbrf7jd5I0qFDh+qM8rSUy+WSw+FQVlZWm44DAACCV1AHnaioKGVkZCg/P99ve35+voYMGdKmYzudTpWUlKigoKBNxwEAAMHL9Kmr48ePa+/evb7npaWlKi4uVrdu3ZSSkqLc3FxNmjRJmZmZuvLKK7V8+XKVlZVp2rRpJlYNAABCgelBp7CwUCNGjPA9z83NlSRNnjxZq1at0sSJE3XkyBHl5eWpvLxc6enpWr9+vVJTU80qGQAAhIiguo6OGZp7Hj4AAAgelriOTntiMTIAANbHiA4jOgAAhJzmfn6bvkYHQPAzDEMnq71tPk5MpL3JOw0DQCARdAA06WS1V45HN7T5OCV5oxUbxZ8dAB0nbNfoAAAA6wvbr1Yul0sul0teb9uH44FwUjh3lGKj7M3ev9LjVeb8je1YEQA0LGyDjtPplNPp9C1mAtA8sVF2pp8AhAymrgAAgGURdAAAgGURdAAAgGURdAAAgGWFbdDhFhAAAFhf2AYdp9OpkpISFRQUmF0KAABoJ2EbdAAAgPURdAAAgGURdAAAgGURdAAAgGURdAAAgGWFbdDh9HIAAKwvbIMOp5cDAGB9YRt0AACA9RF0AACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZYVt0OGCgQAAWF/YBh0uGAgAgPWFbdABAADWR9ABAACWRdABAACWRdABAACWRdABAACWRdABAACWRdABAACWRdABAACWRdABAACWFbZBh1tAAABgfWEbdLgFBAAA1he2QQcAAFgfQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFhWhNkFAOgYhmHoZLW3Va+t9LTudQBgNoIOECZOVnvleHSD2WUAQIdi6goAAFgWIzpAGCqcO0qxUfZWvTYmsnWvAwAzEHSAMBQbZVdsFP/7A7A+pq4AAIBlhW3QcblccjgcysrKMrsUAADQTsI26DidTpWUlKigoMDsUgAAQDsJ26ADAACsj6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsi6ADAAAsK8LsAgCEj0qPt9WvjYm0y2azBbAaAOGAoAOgw2TO39jq15bkjVZsFH+yALQMU1cAAMCy+HoEoF3FRNpVkje6Va+t9HjbNAoEAAQdAO3KZrMx5QTANJaYuiotLdWIESPkcDjUv39/nThxwuySAABAELDE16zbbrtN8+fP19ChQ3X06FFFR0ebXRIAAAgCIR90du7cqcjISA0dOlSS1K1bN5MrAgAAwcL0qavNmzdr/PjxSkpKks1m07p16+rss3TpUvXu3VudO3dWRkaGtmzZ4mvbs2ePunTpouuuu06DBw/WU0891YHVAwCAYGZ60Dlx4oQGDhyoJUuW1Nu+Zs0azZw5U3PmzNH27ds1dOhQ5eTkqKysTJJUXV2tLVu2yOVy6dNPP1V+fr7y8/M78i0AAIAgZXrQycnJ0fz58zVhwoR62xctWqQpU6Zo6tSp6tevnxYvXqzk5GQtW7ZMktSrVy9lZWUpOTlZ0dHRGjNmjIqLixv8eVVVVXK73X4PAABgTaYHncZ4PB4VFRUpOzvbb3t2dra2bt0qScrKytLBgwd17Ngx1dTUaPPmzerXr1+Dx1ywYIHi4+N9j+Tk5HZ9DwAAwDxBHXQOHz4sr9erxMREv+2JiYk6cOCAJCkiIkJPPfWUhg0bpgEDBuiiiy7SuHHjGjzm7NmzVVFR4Xvs37+/Xd8DAAAwT0icdXX2jfwMw/DblpOTo5ycnGYdKzo6mtPPAQAIE0E9opOQkCC73e4bval16NChOqM8AAAAZwvqoBMVFaWMjIw6Z1Hl5+dryJAhbTq2y+WSw+FQVlZWm44DAACCl+lTV8ePH9fevXt9z0tLS1VcXKxu3bopJSVFubm5mjRpkjIzM3XllVdq+fLlKisr07Rp09r0c51Op5xOp9xut+Lj49v6NgAAQBAyPegUFhZqxIgRvue5ubmSpMmTJ2vVqlWaOHGijhw5ory8PJWXlys9PV3r169XamqqWSUDAIAQYXrQGT58uAzDaHSf6dOna/r06R1UEQAAsIqgXqMDAADQFmEbdFiMDACA9YVt0HE6nSopKVFBQYHZpQAAgHYStkEHAABYH0EHAABYFkEHAABYVtgGHRYjAwBgfWEbdFiMDACA9YVt0AEAANZH0AEAAJZF0AEAAJZF0AEAAJYVtkGHs64AALC+sA06nHUFAID1hW3QAQAA1kfQAQAAlkXQAQAAlkXQAQAAlhVhdgEAms8wDJ2s9tbbVuk5rcz5myRJhXNHKjYq4qz2+l8HAFYWtkHH5XLJ5XLJ6+WPP0LHyWqvHI9uaHK/2sADAOEubKeuOL0cAADrC9sRHSDUFc4dpdgou+95U1NXPxUTaW+wDQCshKADhKjYKLtfmImNitC+p8eaWBEABJ+wnboCAADWR9ABAACWRdABAACWRdABAACWFbZBx+VyyeFwKCsry+xSAABAOwnboMN1dAAAsL6wDToAAMD6CDoAAMCyCDoAAMCyCDoAAMCyCDoAAMCyCDoAAMCyCDoAQlKl57TSZr2rtFnvqtJz2uxyAAQp7l4OICRUerxnPT9d77/rExNpl81ma5e6AAQ3gg6AkJA5f2MjbZsafW1J3mjFRvHnDghHYTt1xS0gAACwvrD9iuN0OuV0OuV2uxUfH292OQDqERNpV0ne6HrbKj2nfSM5hXNH1hmxqfR4Gx0FAhAewjboAAh+NputWVNOsVERTE0BqBd/GQCEpNioCO17eqzZZQAIcmG7RgcAAFgfQQcAAFgWQQcIIlwEDwACi6ADAAAsi8XIQAcyDEMnq70Ntjd1td+zrw4MAGgcQQfoQCervXI8uqFZ+zZ1tV8AQNOYugIAAJbFiA5gksK5oxQbZffb1tTVfn8qJtLeYBsA4AyCDmCS2Ch7o0GGq/0CQNvxVxQIIlztFwACizU6AADAsgg6AMIOF2YEwkfYBh2XyyWHw6GsrCyzSwEAAO0kbNfoOJ1OOZ1Oud1uxcfHm10OgHZ09oUWm7ow40/FRNpls9napS4A7S9sgw6A8JE5f2MjbY1fmLEkbzRnvwEhLGynrgAAgPXxNQWAJcVE2lWSN7retqYuzFjp8TY6CgQgdBB0AFiSzWZr1pQTF2YErI2pKwAAYFl8jQEQdrgCNRA+GNEBAoyL0QFA8CDoAAAAyyLoAAAAy2pV0OnTp4+OHDlSZ/sPP/ygPn36tLkoIJgZhqFKz+lGH7XqtnkbOTIAINBatRh537598nrr/sGuqqrSd9991+aigGB2storx6MbmrVvU1fdBQC0rxYFnbffftv37w0bNvjdI8rr9WrTpk1KS0sLWHEAAABt0aKgc/3110s6cyGuyZMn+7VFRkYqLS1Nv/3tbwNWHBDsCueOUmyU3W9bU1fdrRUTaa93OwAgcFoUdGpqaiRJvXv3VkFBgRISEtqlKKur9Jz2TX1ww8DQFhtlb/S/H1fdBQBzteovcGlpaaDrAAAACLhWf9XctGmTNm3apEOHDvlGemqtXLmyzYUBoYqr7gJA8GhV0Hn88ceVl5enzMxM9ezZUzabLdB1hTTDMHSyuuHTiM8+/bghMZF2+hYAgDZoVdB58cUXtWrVKk2aNCnQ9VhCoE4/Zv0OAABt06oLBno8Hg0ZMiTQtQBBg/tVAYA1tGq4YOrUqXrttdf0yCOPBLqeVomIiFB6erokKTMzU3/4wx9MrujvWnr6caXHq8z5Gzu0RgAArKpVQefUqVNavny5Nm7cqAEDBigyMtKvfdGiRQEprrm6du2q4uLiDv2ZzcXpxwAAmKdVn7BffvmlBg0aJEn66quv/NpYPItQ0JYF49yvCgBCR6uCzocffhiwAjZv3qyFCxeqqKhI5eXleuutt3xXYK61dOlSLVy4UOXl5br00ku1ePFiDR061NfudruVkZGhmJgYPfnkk7rmmmsCVp+VcKHCv+N+VQAQHkz/pDtx4oQGDhyo22+/Xb/85S/rtK9Zs0YzZ87U0qVLddVVV+mll15STk6OSkpKlJKSIunMTUaTkpL01VdfaezYsdqxY4fi4uLq/XlVVVWqqqryPXe73e3zxhrRluusEFYAAGi+Vn1KjhgxotEpqg8++KDZx8rJyVFOTk6D7YsWLdKUKVM0depUSdLixYu1YcMGLVu2TAsWLJAkJSUlSZLS09PlcDi0e/duZWZm1nu8BQsW6PHHH292fbA+7lcFANbVqqBTuz6nVnV1tYqLi/XVV1/VudlnW3g8HhUVFWnWrFl+27Ozs7V161ZJ0rFjxxQbG6vo6Gh9++23KikpUZ8+fRo85uzZs5Wbm+t77na7lZycHLCaEXpYMA4A1tWqv97PP/98vdvnzZun48ePt6mgnzp8+LC8Xq8SExP9ticmJurAgQOSpF27dumuu+5Sp06dZLPZ9MILL6hbt24NHjM6OlrR0dEBqzGYtOcVmZkyAwCEooB+Wt1yyy26/PLL9dxzzwXysHU+dA3D8G0bMmSIduzYEdCfFyzqO7unqbOBmnsNHq7IDLQOoR8ILQH9P/TTTz9V586dA3a8hIQE2e123+hNrUOHDtUZ5bGipkILZwM1ri0fSNyYEwCsoVVBZ8KECX7PDcNQeXm5CgsLA3q15KioKGVkZCg/P1833HCDb3t+fr5+8YtftOnYLpdLLpdLXq81r4lilSsy8+0ZANAWrfrUiI+P93veqVMn9e3bV3l5ecrOzm7RsY4fP669e/f6npeWlqq4uFjdunVTSkqKcnNzNWnSJGVmZurKK6/U8uXLVVZWpmnTprWmdB+n0ymn0ym3213n/ZgpJtKukrzRDba35Gygxs6Ma2yBbUunzFrycwEA6EitCjovv/xywAooLCzUiBEjfM9rz4iaPHmyVq1apYkTJ+rIkSPKy8tTeXm50tPTtX79eqWmpgashmBis9maPWrRXmcDtWXKjFEXWM3Zwb+5oV+qG/wZoQQ6Xpv+LysqKtKuXbtks9nkcDh02WWXtfgYw4cPl2EYje4zffp0TZ8+vbVlwqK4jQM6QmPBv6l1coQZwHyt+j/w0KFDuummm/TRRx+pa9euMgxDFRUVGjFihFavXq3zzjsv0HUGnNXX6NSnsQW2bZkya8v6nraGlUCcZQYAsK5WBZ177rlHbrdbO3fuVL9+/SRJJSUlmjx5smbMmKHXX389oEW2h2Bdo2MWs6bMuOcUglFjwb+pdXLBurAfCFet+rR67733tHHjRl/IkSSHwyGXy9XixchoGU57rh+3cUAgNTf4NxX6A7m+B0DrtCro1NTUKDIyss72yMhI1dTUtLkoWEtzF2C2Nay09iwzoL2wvgcwX6fWvOjaa6/Vvffeq++//9637bvvvtN9992nkSNHBqw4hJfae06d/fh7e9222gfffAEA9WnV14UlS5boF7/4hdLS0pScnCybzaaysjL1799fr776aqBrbBfhuBi5LZgyA5qH9T1AcGlV0ElOTta2bduUn5+vv/71rzIMQw6HQ6NGjQp0fe2GxcgAWqOp0B+o9T0AAqNFU1cffPCBHA6H3G63JOnnP/+57rnnHs2YMUNZWVm69NJLtWXLlnYpFAAAoKVa9HVi8eLFuvPOOxUXF1enLT4+XnfddZcWLVqkoUOHBqxAhJbW3HG9MW2ZMmO6DQDQoqDzxRdf6JlnnmmwPTs7W88991ybi0Lo4o7rAIBg0qKgc/DgwXpPK/cdLCJC//d//9fmogDAihhlBDpei4LOBRdcoB07dujCCy+st/3LL79Uz549A1JYe+Osq8AJ5B3XAQAIpBYtRh4zZoweffRRnTp1qk7byZMn9dhjj2ncuHEBK649OZ1OlZSUqKCgwOxSQl7tWSaNPWpxLRwAQEdq0YjO3LlztXbtWl188cW6++671bdvX9lsNu3atcs3OjJnzpz2qhUAAKBFWhR0EhMTtXXrVv3bv/2bZs+eLcMwJJ35Rj969GgtXbpUiYmJ7VIoAABAS7X4alWpqalav369jh07pr1798owDF100UU699xz26M+AACAVmv1ZTnPPfdcZWVlBbIWAACAgOL642h3nFILADBLq+5ebgUul0sOh4NRKQAALCxsgw6nlwMIF5We00qb9a7SZr1b5zYsgNWFbdABAADWR9ABAACWRdABAACWRdABAACWRdABgBDAgmKgdQg6AADAsrhgIABYgGEYOlntrbftpyNATY0GxUTaZbPZAlobYKawDToul8t3x3UACHUnq71yPLqhyf0y529qtL0kb7Rio8L2owEWFLZTV1wwEAAA6yO2A4DFFM4dpdgou+95pee0bySncO7IOiM2lR6vMudv7NAagY5C0AEAi4mNsjc4/RQbFcHUFMIKv+0AEIQqPd6znje+oPjs/QGcQdABgCDU2FRSUwuKAfxd2C5GBgAA1seIDgAEiZhIu0ryRtfb1tSC4rOPEyiVntO+09Y59RyhiN9YAAgSNputWUGCBcVA8zF1BQAALIuvBABgcbFREdr39FizywBMEbYjOi6XSw6HQ1lZWWaXAgCWxB3XEQzCNuhwCwgAAKwvbIMOAACwPoIOAMBSmDLDTxF0AACAZXHWFQCEgI46c6ql99j6qZhIu2w2W7vUBbQWQQcA4NOWe2xx5WQEI6auAACAZRG9ASDMteUeW5Ueb6OjQIDZCDoAEOa4xxasjN9YAEBAsJAZwYigAwAICBYyIxixGBkAAFgW8RkA0GpmLWQ2DEMnq731tjFlhp8i6AAAWs2shcwnq71yPLqhyf2YMkPYTl25XC45HA5lZWWZXQoA4CzheL+qcHzPHSFsY6zT6ZTT6ZTb7VZ8fLzZ5QAAWqlw7ijFRtl9z7n2D34qbIMOAMAaYqPsDU4/ce0f8F8fANCgjrqZKNBewnaNDgAAsD6CDgAAsCyCDgAAsCyCDgAAsCyCDgAAsCzOugIAtAvO2EIwYEQHAABYFkEHAABYFkEHAABYFmt0AACWwtog/BRBBwBgukqP96znp+v9d0P7B66O03I8ukGSVJI3ukX3yWrqtYZh6GR1w3U39Z5rxUTaZbPZml1XuCPoAABM19jdxGvvRB7qTlZ7fUGoKY2955YGsHDHGh0AAGBZREIAgCliIu0qyRtdb1ul57RvVKNw7shGRzBiIu3tUl97Kpw7SrFR/nU39p4rPd5GR73QMIIOAMAUNputWVMwsVERlpuqiY2yN/qerPiezcLUFQAAsCyCDgAAIa7Sc1pps95V2qx3Gz1jKxxZJuhUVlYqNTVV999/v9mlAACAIGGZoPPkk0/qiiuuMLsMAAAQRCwRdPbs2aO//vWvGjNmjNmlAACAIGJ60Nm8ebPGjx+vpKQk2Ww2rVu3rs4+S5cuVe/evdW5c2dlZGRoy5Ytfu3333+/FixY0EEVAwCAUGF60Dlx4oQGDhyoJUuW1Nu+Zs0azZw5U3PmzNH27ds1dOhQ5eTkqKysTJL05z//WRdffLEuvvjijiwbAACEANNP0s/JyVFOTk6D7YsWLdKUKVM0depUSdLixYu1YcMGLVu2TAsWLND//M//aPXq1XrjjTd0/PhxVVdXKy4uTo8++mi9x6uqqlJVVZXvudvtDuwbAgAAQcP0EZ3GeDweFRUVKTs72297dna2tm7dKklasGCB9u/fr3379um5557TnXfe2WDIqd0/Pj7e90hOTm7X9wAAQDCz+qnpQR10Dh8+LK/Xq8TERL/tiYmJOnDgQKuOOXv2bFVUVPge+/fvD0SpAAAgCJk+ddUcZ9+O3jCMem9Rf9tttzV5rOjoaEVHRweqNAAAEMSCekQnISFBdru9zujNoUOH6ozyAAAAnC2og05UVJQyMjKUn5/vtz0/P19Dhgxp07FdLpccDoeysrLadBwAABC8TJ+6On78uPbu3et7XlpaquLiYnXr1k0pKSnKzc3VpEmTlJmZqSuvvFLLly9XWVmZpk2b1qaf63Q65XQ65Xa7FR8f39a3AQAIoNioCO17emy7/5xKj/es56fr/XfzjtX4a8/+WWfrqPccbkwPOoWFhRoxYoTveW5uriRp8uTJWrVqlSZOnKgjR44oLy9P5eXlSk9P1/r165WammpWyQAAi8icv7GRtk1tOG7rX4vAMj3oDB8+XIZhNLrP9OnTNX369A6qCAAAWIXpQQcAgI4UE2lXSd7oetsqPad9ozGFc0cqNqr5H5MteW1MpL0FFaMtwjbouFwuuVwueb2Nz5kCAKzFZrM1K8DERkW0KOgE6rUIrKA+66o9OZ1OlZSUqKCgwOxSAAAISaFwVeWwDToAAMD6CDoAAMCymEAEACCE1Hc9nuZe/ycm0l7vLZSsLGyDDouRAQChqLFr/5xpb/gaPiV5o8NukXTYTl2xGBkAAOsLr1gHAEAIauzaP1Lj1/Cp9HibHAWyMoIOAABBrrnX/pG4hs/ZwnbqCgAAWB9BBwAAWFbYBh2XyyWHw6GsrCyzSwEAAO0kbIMOZ10BAGB9YRt0AACA9RF0AACAZXH+GQAAYSIcbx9B0AEAIEyE4+0jmLoCAACWFXrRLEC4qScAIByE++0jwjboOJ1OOZ1Oud1uxcfHm10OAADtItxvH2GtdwMAQBvERkVo39NjO/y1aD+s0QEAAJZF0AEAAJZF0AEAAJZF0AEAAJZF0AEAAJYVtkHH5XLJ4XAoKyvL7FIAAEA7Cdug43Q6VVJSooKCArNLAQAA7SRsgw4AALA+gg4AALAsgg4AALAsbgEBAECI64jbT1R66t4Eu9Jzut5/ny0m0i6bzdYudTWFoAMAAJrU1F3Ma++AXp+SvNGm3SyUqSsAAGBZjOgAAIB6xUTaVZI3usH2Ss9p30hO4dyRfqM2lR5vk6NAHYGgAwAA6mWz2Zo95RQbFWHa9FRjmLoCAACWFbZBh1tAAABgfWEbdLgFBAAA1he2QQcAAFgfQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQQcAAFhW8N19CwAAdJjYqAjte3qs2WW0G0Z0AACAZRF0AACAZRF0AACAZRF0AACAZYVt0HG5XHI4HMrKyjK7FAAA0E7CNug4nU6VlJSooKDA7FIAAEA7CdugAwAArI+gAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALIugAwAALCvC7AIAAEBoio2K0L6nx5pdRqNCfkTnxx9/VFZWlgYNGqT+/ftrxYoVZpcEAACCRMiP6MTGxurjjz9WbGysKisrlZ6ergkTJqh79+5mlwYAAEwW8iM6drtdsbGxkqRTp07J6/XKMAyTqwIAAMHA9KCzefNmjR8/XklJSbLZbFq3bl2dfZYuXarevXurc+fOysjI0JYtW/zaf/jhBw0cOFC9evXSgw8+qISEhA6qHgAABDPTg86JEyc0cOBALVmypN72NWvWaObMmZozZ462b9+uoUOHKicnR2VlZb59unbtqi+++EKlpaV67bXXdPDgwY4qHwAABDHTg05OTo7mz5+vCRMm1Nu+aNEiTZkyRVOnTlW/fv20ePFiJScna9myZXX2TUxM1IABA7R58+YGf15VVZXcbrffAwAAWJPpQacxHo9HRUVFys7O9tuenZ2trVu3SpIOHjzoCytut1ubN29W3759GzzmggULFB8f73skJye33xsAAACmCuqgc/jwYXm9XiUmJvptT0xM1IEDByRJ3377rYYNG6aBAwfq6quv1t13360BAwY0eMzZs2eroqLC99i/f3+7vgcAAGCekDi93Gaz+T03DMO3LSMjQ8XFxc0+VnR0tKKjowNZHgAACFJBPaKTkJAgu93uG72pdejQoTqjPAAAAGcL6qATFRWljIwM5efn+23Pz8/XkCFD2nRsl8slh8OhrKysNh0HAAAEL9Onro4fP669e/f6npeWlqq4uFjdunVTSkqKcnNzNWnSJGVmZurKK6/U8uXLVVZWpmnTprXp5zqdTjmdTrndbsXHx7f1bdThPfmjqsp2yDCym94ZAAALMgxDJ3d/qqNHsxTb43xTajA96BQWFmrEiBG+57m5uZKkyZMna9WqVZo4caKOHDmivLw8lZeXKz09XevXr1dqaqpZJTfq6NGjeva53+rA8sU6fapSO3f+sy4ffJnZZQEA0OGqD3+j/1v3lBybfq97771X9913n7p169ahNZg+dTV8+HAZhlHnsWrVKt8+06dP1759+1RVVaWioiINGzbMvIIbcPToUT3yyCNKS0vT719YrF+MzZEkVVdXm1wZAADmMLynJZ25LMyiRYuUlpamRx55REePHu2wGkwPOmYJ1BqdnwacRYsWadq0aSotLdXDDz8coEoBAAhtDz/8sEpLSzVt2rQODzw2I8zvgFlRUaGuXbtq//79iouLa9Hrfv/73+vFF1+U1+uV0+nU/fffr/PPPzMHuW3bNmVkZGjDxg/0j1kZ7VU+AABBqdJzWoPuXaGDr81SUVGRBg8eLOnMmdPPPfecXC6X7Ha7pk2bpnvuuafF62XdbreSk5P1ww8/NP5aI8zt37/fkNTqx3nnnWfs2LGjznGLioradFwePHjw4MHDKo+ioqI6n5M7duwwzjvvvDYfe//+/Y1+zpu+GNlsSUlJ2r9/v372s5/VuTBhY346onPFFVfUGdGp9fHHH2vQoEEBrrp91Kbjlo5uoWXo545DX3cM+rljhGI/FxcX65prrvHbdujQIS1cuFBLly6V3W7XAw880KoRHcMw9OOPPyopKanR/cI+6HTq1Em9evVq8evi4uL07LPPatasWXr++ef1wgsvyOVy+QJPrS5duoTML2StuLi4kKs5FNHPHYe+7hj0c8cIpX7u0qWL799nB5zc3Nw2n4XVnHAUtouRA6Vbt2564okntG/fPuXm5urFF19U79699dRTT5ldGgAAQeGpp55S79699dJLLyk3N1f79u3TE0880SGnmhN0AuTswPP+++9LkiIjI02uDAAAc9R+Br7//vsdHnBqhf3UVaDVBp777rtPH330kdLT080uqdmio6P12GOPcdPTdkY/dxz6umPQzx0jFPs5PT1db775poYPH97hFwqsFfanlwMAAOti6goAAFgWQQcAAFgWQQcAAFgWQQcAAFgWQScMfffdd7rlllvUvXt3xcbGatCgQSoqKvK1G4ahefPmKSkpSTExMRo+fLh27txpYsWh6fTp05o7d6569+6tmJgY9enTR3l5eaqpqfHtQ1+33ObNmzV+/HglJSXJZrNp3bp1fu3N6dOqqirdc889SkhI0DnnnKPrrrtO3377bQe+i+DXWD9XV1froYceUv/+/XXOOecoKSlJt956q77//nu/Y9DPTWvq9/mn7rrrLtlsNi1evNhvO/3cOIJOmDl27JiuuuoqRUZG6r//+79VUlKi3/72t+ratatvn2effVaLFi3SkiVLVFBQoB49eujnP/+5fvzxR/MKD0HPPPOMXnzxRS1ZskS7du3Ss88+q4ULF+r3v/+9bx/6uuVOnDihgQMHasmSJfW2N6dPZ86cqbfeekurV6/WX/7yFx0/flzjxo2T1+vtqLcR9Brr58rKSm3btk2PPPKItm3bprVr12r37t267rrr/Pajn5vW1O9zrXXr1umzzz6r93YH9HMTAnBfTISQhx56yLj66qsbbK+pqTF69OhhPP30075tp06dMuLj440XX3yxI0q0jLFjxxp33HGH37YJEyYYt9xyi2EY9HUgSDLeeust3/Pm9OkPP/xgREZGGqtXr/bt89133xmdOnUy3nvvvQ6rPZSc3c/1+fzzzw1JxjfffGMYBv3cGg3187fffmtccMEFxldffWWkpqYazz//vK+Nfm4aIzph5u2331ZmZqZ+9atf6fzzz9dll12mFStW+NpLS0t14MABZWdn+7ZFR0frmmuu0datW80oOWRdffXV2rRpk3bv3i1J+uKLL/SXv/xFY8aMkURft4fm9GlRUZGqq6v99klKSlJ6ejr93gYVFRWy2Wy+0WH6OTBqamo0adIkPfDAA7r00kvrtNPPTePKyGHmf//3f7Vs2TLl5ubq4Ycf1ueff64ZM2YoOjpat956qw4cOCBJSkxM9HtdYmKivvnmGzNKDlkPPfSQKioqdMkll8hut8vr9erJJ5/UzTffLEn0dTtoTp8eOHBAUVFROvfcc+vsU/t6tMypU6c0a9Ys/frXv/bdbJJ+DoxnnnlGERERmjFjRr3t9HPTCDphpqamRpmZmb6bjl522WXauXOnli1bpltvvdW3n81m83udYRh1tqFxa9as0auvvqrXXntNl156qYqLizVz5kwlJSVp8uTJvv3o68BrTZ/S761TXV2tm266STU1NVq6dGmT+9PPzVdUVKQXXnhB27Zta3Gf0c9/x9RVmOnZs6ccDofftn79+qmsrEyS1KNHD0mq803g0KFDdb4lo3EPPPCAZs2apZtuukn9+/fXpEmTdN9992nBggWS6Ov20Jw+7dGjhzwej44dO9bgPmie6upq3XjjjSotLVV+fr5vNEeinwNhy5YtOnTokFJSUhQREaGIiAh98803+s1vfqO0tDRJ9HNzEHTCzFVXXaWvv/7ab9vu3buVmpoqSerdu7d69Oih/Px8X7vH49HHH3+sIUOGdGitoa6yslKdOvn/L2a3232nl9PXgdecPs3IyFBkZKTfPuXl5frqq6/o9xaoDTl79uzRxo0b1b17d792+rntJk2apC+//FLFxcW+R1JSkh544AFt2LBBEv3cLGauhEbH+/zzz42IiAjjySefNPbs2WP8x3/8hxEbG2u8+uqrvn2efvppIz4+3li7dq2xY8cO4+abbzZ69uxpuN1uEysPPZMnTzYuuOAC47/+67+M0tJSY+3atUZCQoLx4IMP+vahr1vuxx9/NLZv325s377dkGQsWrTI2L59u+9sn+b06bRp04xevXoZGzduNLZt22Zce+21xsCBA43Tp0+b9baCTmP9XF1dbVx33XVGr169jOLiYqO8vNz3qKqq8h2Dfm5aU7/PZzv7rCvDoJ+bQtAJQ++8846Rnp5uREdHG5dccomxfPlyv/aamhrjscceM3r06GFER0cbw4YNM3bs2GFStaHL7XYb9957r5GSkmJ07tzZ6NOnjzFnzhy/DwL6uuU+/PBDQ1Kdx+TJkw3DaF6fnjx50rj77ruNbt26GTExMca4ceOMsrIyE95N8Gqsn0tLS+ttk2R8+OGHvmPQz01r6vf5bPUFHfq5cTbDMIyOGz8CAADoOKzRAQAAlkXQAQAAlkXQAQAAlkXQAQAAlkXQAQAAlkXQAQAAlkXQAQAAlkXQAQAAlkXQAWBp8+bNk81mk81m0+LFi02rY9++fb46Bg0aZFodQLgh6ABos9tuu002m03Tpk2r0zZ9+nTZbDbddtttHV/Y/3fppZeqvLxc//qv/+q3ffv27Zo4caJ69uyp6Ohopaamaty4cXrnnXfU3IvGjx8/XqNGjaq37dNPP5XNZtO2bduUnJys8vJy/eY3v2nz+wHQfAQdAAGRnJys1atX6+TJk75tp06d0uuvv66UlBQTK5MiIiLUo0cPxcbG+rb9+c9/1j/+4z/q+PHjeuWVV1RSUqI33nhD119/vebOnauKiopmHXvKlCn64IMP9M0339RpW7lypQYNGqTBgwfLbrerR48e6tKlS8DeF4CmEXQABMTgwYOVkpKitWvX+ratXbtWycnJuuyyy/z2fe+993T11Vera9eu6t69u8aNG6e//e1vvnaPx6O7775bPXv2VOfOnZWWlqYFCxb42ufNm6eUlBRFR0crKSlJM2bMaFGtJ06c0JQpUzR27Fi9++67ys7O1j/8wz/o8ssv19SpU/XFF18oPj7et39JSYnGjBmjLl26KDExUZMmTdLhw4clSePGjdP555+vVatW+f2MyspKrVmzRlOmTGlRbQACi6ADIGBuv/12vfzyy77nK1eu1B133FFnvxMnTig3N1cFBQXatGmTOnXqpBtuuEE1NTWSpN/97nd6++239ac//Ulff/21Xn31VaWlpUmS/vM//1PPP/+8XnrpJe3Zs0fr1q1T//79W1Tn+++/ryNHjujBBx9scB+bzSZJKi8v1zXXXKNBgwapsLBQ7733ng4ePKgbb7xR0pnRoltvvVWrVq3ym+5644035PF49C//8i8tqg1AYEWYXQAA65g0aZJmz57tW3j7ySefaPXq1froo4/89vvlL3/p9/zf//3fdf7556ukpETp6ekqKyvTRRddpKuvvlo2m02pqam+fcvKytSjRw+NGjVKkZGRSklJ0eWXX96iOnfv3i1J6tu3r29bQUGBRowY4Xu+evVqjRs3TsuWLdPgwYP11FNP+dpWrlyp5ORk7d69WxdffLHuuOMOLVy4UB999JHvGCtXrtSECRN07rnntqg2AIHFiA6AgElISNDYsWP1yiuv6OWXX9bYsWOVkJBQZ7+//e1v+vWvf60+ffooLi5OvXv3lnQmxEhnFjcXFxerb9++mjFjht5//33fa3/1q1/p5MmT6tOnj+6880699dZbOn36dJtrHzBggIqLi1VcXKwTJ074jllUVKQPP/xQXbp08T0uueQS3/uQpEsuuURDhgzRypUrfdu3bNlS72gWgI5F0AEQUHfccYdWrVqlV155pcEP+vHjx+vIkSNasWKFPvvsM3322WeSzqzNkc6s9yktLdUTTzyhkydP6sYbb9Q///M/Szqz6Pnrr7+Wy+VSTEyMpk+frmHDhqm6urrZNV500UWSpK+//tq3LTo6WhdeeKEuvPBCv31ramo0fvx4XwiqfezZs0fDhg3z7TdlyhS9+eabcrvdevnll5WamqqRI0c2uyYA7YOgAyCg/umf/kkej0cej0ejR4+u037kyBHt2rVLc+fO1ciRI9WvXz8dO3aszn5xcXGaOHGiVqxYoTVr1ujNN9/U0aNHJUkxMTG67rrr9Lvf/U4fffSRPv30U+3YsaPZNWZnZ6tbt2565plnmtx38ODB2rlzp9LS0nxBqPZxzjnn+Pa78cYbZbfb9dprr+mVV17R7bff7lvnA8A8rNEBEFB2u127du3y/fts5557rrp3767ly5erZ8+eKisr06xZs/z2ef7559WzZ08NGjRInTp10htvvKEePXqoa9euWrVqlbxer6644grFxsbqj3/8o2JiYvzW8TSlS5cu+sMf/qCJEydq7NixmjFjhi666CIdP35c7733nl/tTqdTK1as0M0336wHHnhACQkJ2rt3r1avXq0VK1b49uvSpYsmTpyohx9+WBUVFaZeNwjA3zGiAyDg4uLiFBcXV29bp06dtHr1ahUVFSk9PV333XefFi5c6LdPly5d9MwzzygzM1NZWVnat2+f1q9fr06dOqlr165asWKFrrrqKg0YMECbNm3SO++8o+7du7eoxhtuuEFbt25VbGysbr31VvXt21fXXnutPvjgA99CZElKSkrSJ598Iq/Xq9GjRys9PV333nuv4uPj1amT/5/QKVOm6NixYxo1apTp1w4CcIbNaO7lPwEgBM2bN0/r1q1TcXGx2aVICr56AKtjRAeA5e3YsUNdunTR0qVLTauhrKxMXbp08TtNHUD7Y0QHgKUdPXrUt4j5vPPO87vicUc6ffq09u3bJ+nMGV7Jycmm1AGEG4IOAACwLKauAACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZRF0AACAZf0/4fC3hMHDAZgAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(lhe_data)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "e81691f4-31b1-41a2-a503-1393e3019faa",
   "metadata": {},
   "outputs": [],
   "source": [
    "# events were consumed so reload here\n",
    "lhe_file = pylhe.read_lhe_file(lhe_data)\n",
    "\n",
    "\n",
    "def filtered(events):\n",
    "    for e in events:\n",
    "        for p in e.particles:\n",
    "            # only keep Z bosons\n",
    "            if p.id == 23:\n",
    "                yield e\n",
    "\n",
    "\n",
    "pylhe.write_lhe_file(lhe_file.init, filtered(lhe_file.events), \"filtered.lhe.gz\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "ec72748b-b1c6-425b-b2e6-e485e3b2e39f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjoAAAGwCAYAAACgi8/jAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8hTgPZAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAp00lEQVR4nO3de3xU9Z3/8fdkcjExJtxqQiQXqIoEFJTEruANwdBwe1htRbsCCriyjEWktYJSpVkUKBZpy4DCFvFR14V1C1QXFgqoiy1ekkAUGCrQDQYkyIo0ARIyyeT8/vDH1CH3ZJKZfM/r+XjkQc75njnzme8jM/Pmey5fh2VZlgAAAAwUEeoCAAAA2gtBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWJGhLiDUamtrdfz4cV122WVyOByhLgcAADSDZVk6c+aMUlJSFBHR8LiN7YPO8ePHlZqaGuoyAABAKxw9elS9evVqsN32Qeeyyy6T9HVHJSQkhLgaAADQHOXl5UpNTfV/jzfE9kHnwuGqhIQEgg4AAJ1MU6edcDIyAAAwlm2DjtvtVmZmprKzs0NdCgAAaCcOu89eXl5ersTERJWVlXHoCgCATqK539+2HdEBAADmI+gAAABjEXQAAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMaybdBhCggAAMzHFBBMAQE0ybIsVVb72ryf2ChnkzMNA0BzNPf7O7IDawLQSVVW+5T5zNY278eTN1Jx0XzsAOg4tj10BQAAzMd/rQC0SMHcEYqLdjZ7+wqvT1nzt7djRQDQMIIOgBaJi3Zy+AlAp8GhKwAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsWwbdJjUEwAA89k26LhcLnk8HuXn54e6FAAA0E5sG3QAAID5CDoAAMBYBB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWLYNOm63W5mZmcrOzg51KQAAoJ3YNui4XC55PB7l5+eHuhQAANBObBt0AACA+Qg6AADAWAQdAABgLIIOAAAwFkEHAAAYi6ADAACMRdABAADGIugAAABjEXQAAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWJGhLiAYiouLNXnyZH3xxRdyOp364IMPdOmll4a6LCCsWJalympfvW0V3hplzd8hSSqYO1xx0ZEXtdf/OAAId0YEnQcffFDz58/XLbfcoq+++koxMTGhLgkIO5XVPmU+s7XJ7S4EHgAwQacPOvv371dUVJRuueUWSVK3bt1CXBEAAAgXIQ86O3fu1OLFi1VYWKjS0lJt2LBBd911V8A2y5cv1+LFi1VaWqr+/ftr6dKl/mBz6NAhxcfHa9y4cTp27Ji+//3v66mnngrBKwE6j4K5IxQX7fQvN3Xo6ptio5wNtgFAuAl50Dl37pwGDhyohx56SPfcc0+d9nXr1mnmzJlavny5hg4dqpdfflm5ubnyeDxKS0tTdXW13nvvPRUVFenyyy/Xd7/7XWVnZ+vOO++s9/mqqqpUVVXlXy4vL2+31waEq7hoZ0CYiYuO1JGFo0NYEQC0j5BfdZWbm6v58+fr7rvvrrd9yZIlmjJliqZOnap+/fpp6dKlSk1N1YoVKyRJvXr1UnZ2tlJTUxUTE6NRo0apqKiowedbsGCBEhMT/T+pqant8bIAAEAYCHnQaYzX61VhYaFycnIC1ufk5GjXrl2SpOzsbH3xxRc6ffq0amtrtXPnTvXr16/Bfc6ZM0dlZWX+n6NHj7brawAAAKET8kNXjfnyyy/l8/mUlJQUsD4pKUknTpyQJEVGRur555/XrbfeKsuylJOTozFjxjS4z5iYGK7KAgDAJsI66FzgcDgCli3LCliXm5ur3Nzcji4LAACEubA+dNWjRw85nU7/6M0FJ0+erDPKAwAAcLGwDjrR0dEaPHiwtm3bFrB+27ZtGjJkSJv27Xa7lZmZqezs7DbtBwAAhK+QH7o6e/asDh8+7F8uLi5WUVGRunXrprS0NM2aNUsTJkxQVlaWbrrpJq1cuVIlJSWaNm1am57X5XLJ5XKpvLxciYmJbX0ZAAAgDIU86BQUFGjYsGH+5VmzZkmSJk2apDVr1mj8+PE6deqU8vLyVFpaqgEDBmjz5s1KT08PVckAAKCTCHnQuf3222VZVqPbTJ8+XdOnT++gigAAgCnC+hwdAACAtiDoAAAAY9k26HDVFQAA5rNt0HG5XPJ4PMrPzw91KQAAoJ3YNugAAADzEXQAAICxCDoAAMBYBB0AAGAs2wYdrroCAMB8tg06XHUFAID5bBt0AACA+Qg6AEKmwlujjNmblDF7kyq8NaEuB4CBQj6pJwD7qPD6Llquqff3+sRGOeVwONqlLgDmIugA6DBZ87c30raj0cd68kYqLpqPLAAtw6ErAABgLNv+98jtdsvtdsvn8zW9MYBWi41yypM3st62Cm+NfySnYO7wOiM2FV5fo6NAANAU2wYdl8sll8ul8vJyJSYmhrocwFgOh6PBQ05x0ZE6snB0B1cEwE44dAUAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMJZtg47b7VZmZqays7NDXQoAAGgntg06LpdLHo9H+fn5oS4FAAC0E9sGHQAAYD6CDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLNsGHaaAAADAfLYNOkwBAQCA+WwbdAAAgPkIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLNsGHSb1BADAfLYNOkzqCQCA+WwbdAAAgPkIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAIao8NYoY/YmZczepApvTajLAYCwQNABAADGigx1AQCaz7IsVVb76m375ihOfSM6Fd76HwcAJiPoAJ1IZbVPmc9sbXK7rPk7OqAaAAh/HLoCAADGYkQH6KQK5o5QXLTTv1zhrfGP5BTMHa646Ibf3rFRzgbbAMAkBB2gk4qLdjYYZuKiIxsNOgBgFxy6AgAAxiLoAAAAYzG2DRgiLjpSRxaODnUZABBWGNEBAADGsm3QcbvdyszMVHZ2dqhLAQAA7cS2Qcflcsnj8Sg/Pz/UpQAAgHZi26ADAADMR9ABAADGIugAAABjEXQAAICxCDoAOqUKb40yZm9SxuxNqvDWhLocAGGKoAMAAIxF0AEAAMZqVdDp06ePTp06VWf93/72N/Xp06fNRQEAAARDq+a6OnLkiHw+X531VVVV+vzzz9tcFABcrMLru2i5pt7f6xMb5ZTD4WiXugCEtxYFnTfffNP/+9atW5WYmOhf9vl82rFjhzIyMoJWHABckDV/eyNtOxp9rCdvpOKimcMYsKMWvfPvuusuSZLD4dCkSZMC2qKiopSRkaFf/vKXQSsOAACgLVoUdGprayVJvXv3Vn5+vnr06NEuRQGA9PUhJ0/eyHrbKrw1/pGcgrnD64zYVHh9jY4CAbCHVo3lFhcXB7sOAKjD4XA065BTXHQkh6YA1KvVnww7duzQjh07dPLkSf9IzwWrV69uc2GAHVV4a5T5zFZJnFcCAMHQqk/Rn//858rLy1NWVpZ69uzJ1QwAACAstSrovPTSS1qzZo0mTJgQ7HoAoFnioiN1ZOHoUJcBIMy16oaBXq9XQ4YMCXYtAAAAQdWqoDN16lS9/vrrwa4FMJ5lWarw1jT6c0H97XVv1AkAaFirDl2dP39eK1eu1Pbt23XdddcpKioqoH3JkiVBKQ4wTWW1z3+ycVOaugkeAKBprQo6n3zyiQYNGiRJ2rdvX0AbJyYDAIBw0aqg88477wS7DsB2CuaOUFy0M2BdUzfB+6bYKGeDbQCAr3GTDiBE4qKdjQYZboIHAG3Xqk/RYcOGNXqI6u233251QQAAAMHSqqBz4fycC6qrq1VUVKR9+/bVmewTAAAgVFoVdF588cV618+bN09nz55tU0EAAADBEtQTAB544AHdeOONeuGFF4K5W8A2uNsvAARXq24Y2JD3339fl1xySTB32SyRkZEaNGiQBg0apKlTp3b48wMAgPDUqhGdu+++O2DZsiyVlpaqoKBAP/vZz4JSWEt06dJFRUVFHf68AAAgvLUq6CQmJgYsR0REqG/fvsrLy1NOTk5QCgMAAGirVgWdV155JWgF7Ny5U4sXL1ZhYaFKS0u1YcMG3XXXXQHbLF++XIsXL1Zpaan69++vpUuX6pZbbvG3l5eXa/DgwYqNjdVzzz2n2267LWj1AQCAzqtNJyMXFhbqwIEDcjgcyszM1PXXX9/ifZw7d04DBw7UQw89pHvuuadO+7p16zRz5kwtX75cQ4cO1csvv6zc3Fx5PB6lpaVJko4cOaKUlBTt27dPo0eP1t69e5WQkFDv81VVVamqqsq/XF5e3uKaAQBA59CqoHPy5Endd999evfdd9WlSxdZlqWysjINGzZMa9eu1be+9a1m7ys3N1e5ubkNti9ZskRTpkzxn2S8dOlSbd26VStWrNCCBQskSSkpKZKkAQMGKDMzUwcPHlRWVla9+1uwYIF+/vOfN7s+AADQebXqqqsf/ehHKi8v1/79+/XVV1/p9OnT2rdvn8rLyzVjxoygFef1elVYWFjnvJ+cnBzt2rVLknT69Gn/CM2xY8fk8XjUp0+fBvc5Z84clZWV+X+OHj0atHoBdA4V3hplzN6kjNmbVOGtCXU5ANpRq0Z0tmzZou3bt6tfv37+dZmZmXK73UE9GfnLL7+Uz+dTUlJSwPqkpCSdOHFCknTgwAE98sgjioiIkMPh0K9+9St169atwX3GxMQoJiYmaDUCAIDw1aqgU1tbq6ioqDrro6KiVFtb2+aiLnbxvFqWZfnXDRkyRHv37g36cwIAgM6vVYeu7rjjDj322GM6fvy4f93nn3+uxx9/XMOHDw9acT169JDT6fSP3lxw8uTJOqM8AAAAF2tV0Fm2bJnOnDmjjIwMffvb39aVV16p3r1768yZM/rNb34TtOKio6M1ePBgbdu2LWD9tm3bNGTIkDbt2+12KzMzU9nZ2W3aDwAACF+tOnSVmpqq3bt3a9u2bfrLX/4iy7KUmZmpESNGtHhfZ8+e1eHDh/3LxcXFKioqUrdu3ZSWlqZZs2ZpwoQJysrK0k033aSVK1eqpKRE06ZNa03pfi6XSy6XS+Xl5XVugAjALBVe30XLNfX+Xp/YKGedw+cAOo8WBZ23335bjz76qD744AMlJCTozjvv1J133ilJKisrU//+/fXSSy8F3MyvKQUFBRo2bJh/edasWZKkSZMmac2aNRo/frxOnTqlvLw8lZaWasCAAdq8ebPS09NbUjoAG8uav72Rth2NPtaTN1Jx0UGd/xhAB2rRu3fp0qV6+OGH670ZX2Jioh555BEtWbKkRUHn9ttvl2VZjW4zffp0TZ8+vSWlAiFT4a1R5jNbJfElCQCh1qJP4I8//liLFi1qsD0nJ0cvvPBCm4sCgLaKjXLKkzey3rYKb41/JKdg7vA6YbTC62t0FAhA59GioPPFF1/Ue1m5f2eRkfq///u/NhfVEdxut9xut3w+X9MbA+h0HA5Hs0bT4qIjGXUDDNaiq66uuOKKRu9Z88knn6hnz55tLqojuFwueTwe5efnh7oUAADQTloUdEaNGqVnnnlG58+fr9NWWVmpZ599VmPGjAlacQAAAG3RovHauXPnav369br66qv16KOPqm/fvnI4HDpw4ID/MNDTTz/dXrUCYcGyLFVWN3zIs7FLly++zBkA0L5aFHSSkpK0a9cu/fM//7PmzJnjv1rK4XBo5MiRWr58OXcshvEqq33+q6qa0tSlywCA9tXiM/DS09O1efNmnT59WocPH5ZlWbrqqqvUtWvX9qgPAIIuLjpSRxaODnUZADpAqy816Nq1K9MnwPYK5o5QXLQzYF1Tly5fEBvlrHc9ACB4bHtNJZeXIxjiop2NXprMpcsAEFqtmtTTBFxeDgCA+WwbdAAAgPkIOgAAwFgEHQAAYCyCDgAAMBaXgwBBxj1aACB82HZEx+12KzMzk3sBAQBgMNsGHS4vBwDAfLYNOugcKrw1ypi9SRmzN9WZIBMAgKYQdAAAgLEIOjAWo0EAAIIOAAAwFkEHAAAYi6ADAACMRdABAADGsm3Q4YaBaAwnMgOAGWw7BYTL5ZLL5VJ5ebkSExNDXY5tWZalympfg+3fDBmNBY7YKKccDkdQawMAdH62DToID5XVPmU+s7VZ22bN39FgmydvpOKi+XMGAATimwHNUuGt8QeScAwVFd66o0KMBgEAwuvbCrZWMHeE4qKdAesqvDX+kZyCucMDAlaF16es+dslyf9vQxgNAgB74tMdYSMu2tlo4IiLjiSQAABahG8NdFqxUU558kY22N7c0SAAgLkIOui0HA5Hs0d4GA1CsLTlfLVwP9cNMBHvMtheS09krm97AEB4IujA9tpyIjMAILzZ9s7IAADAfLYd0XG73XK73fL5OAxhR205kfni/QAAwpdtg44dp4Bo7ERIu03FwInMAGAPfHpDUvhOxRAXHakjC0cHbX8AAHvhHB20O2YCBwCECiM6qKMtUzGEE0aDAAAEHdTBVAwAAFPwbQUAjbj4BpHNPTG//n01/7Gd5cR+INwRdACgEY0dlm3LzSSbeixTRADBwcnIAADAWPx3AQAu0tgNJZt7M8nWPDZcT+wHOjOCDlAPrtiyt+beULItJ+ZzUj/QMXiXGaQtdzdmRm4AgIkIOgYJ1t2NAQAwhW1PRna73crMzFR2dnaoSwEAAO3EtiM6nXVSz8Ym5vymlt7d+JvaMiN3fYfAOGQGtF1z3/vBfizQ2fHXbqhg3924uSfnNnXFCIfMAAAdybaHrgAAgPkY0UGbNXbPEaljDpkBAFAfgg7arLn3HJG4dwgAoGNx6AoAABiL/1oDQAu05a7Z3HEb6HiM6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZXXdkIV3wAAOyGoBNmLMtSZXXDE10yQSZgDxe/nxt77ze9r+Y/NjbKKYfD0aL9A+GMoBNmKqt9/lmGm8IEmYC5Gpsgty3v/aYey+zmMA1/zWh3HDIDAIQKQSeMFcwdobjowIkumSATMFdjE+Q2973fmsdWeH2NjiABnRlBJ4zFRTsb/TBjgkzALM2dILct730+N2A3XF4OAACMZdug43a7lZmZqezs7FCXAgAA2oltg47L5ZLH41F+fn6oSwEAAO3EtkEHAACYj6ADAACMRdABAADGIugAAABjEXQAAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABiLmd06mbjoSB1ZODrUZQDoYG157/O5ATtjRAcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAADarw1ihj9iZlzN6kCm9Np3gs8E0EHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwljFBp6KiQunp6frJT34S6lIAAECYMCboPPfcc/rOd74T6jIAAEAYMSLoHDp0SH/5y180atSoUJcCAADCSMiDzs6dOzV27FilpKTI4XBo48aNdbZZvny5evfurUsuuUSDBw/We++9F9D+k5/8RAsWLOigigEAQGcR8qBz7tw5DRw4UMuWLau3fd26dZo5c6aefvpp7dmzR7fccotyc3NVUlIiSfrDH/6gq6++WldffXVHlg0AADqByFAXkJubq9zc3AbblyxZoilTpmjq1KmSpKVLl2rr1q1asWKFFixYoA8++EBr167VG2+8obNnz6q6uloJCQl65pln6t1fVVWVqqqq/Mvl5eXBfUEAACBshHxEpzFer1eFhYXKyckJWJ+Tk6Ndu3ZJkhYsWKCjR4/qyJEjeuGFF/Twww83GHIubJ+YmOj/SU1NbdfXAAAAQiesg86XX34pn8+npKSkgPVJSUk6ceJEq/Y5Z84clZWV+X+OHj0ajFIBAEAYCvmhq+ZwOBwBy5Zl1VknSQ8++GCT+4qJiVFMTEywSgMAAGEsrEd0evToIafTWWf05uTJk3VGeQAAAC4W1kEnOjpagwcP1rZt2wLWb9u2TUOGDAlRVQAAoLMI+aGrs2fP6vDhw/7l4uJiFRUVqVu3bkpLS9OsWbM0YcIEZWVl6aabbtLKlStVUlKiadOmtel53W633G63fD5fW18CAAAIUyEPOgUFBRo2bJh/edasWZKkSZMmac2aNRo/frxOnTqlvLw8lZaWasCAAdq8ebPS09Pb9Lwul0sul0vl5eVKTExs074AwBQVXt9FyzX1/t68fXXMY2OjnPWetwlIksOyLCvURYTShaBTVlamhISEUJejCm+NMp/ZKkny5I1UXHTIsygAw33zc6cz4rPSnpr7/R3W5+gAAAC0BREYAGwuNsopT97IetsqvDXKmr9DklQwd3iLRk7a87EVXp+y5m9v9v5gX7YNOpyMDABfczgczQohcdGRrT5EFKrHArY9dOVyueTxeJSfnx/qUgAAQDuxbdABAJipwlujjNmblDF7U4uv9oJ5CDoAAMBYBB0AAGAsgg4AADAWQQcAABjLtkHH7XYrMzNT2dnZoS4FAAC0E9sGHS4vBwDAfLYNOgAAwHwEHQAAYCyCDgAAMBZBBwAAGIugAwAAjGXboMPl5QAAmM+2QYfLywEAMJ9tgw4AADAfQQcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFi2DTrcMBAAAPPZNuhww0AAAMwXGeoCAADhKy46UkcWju5UjwW+ybYjOgAAwHwEHQAAYCyCDgAAMBZBBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsWwbdJgCAgAA89k26DAFBAAA5rNt0AEAAOYj6AAAAGMRdAAAgLEIOgAAwFiRoS4AAIC2qPD6Llquqff3+sRGOeVwONqlLoQHgg4AoFPLmr+9kbYdjT7WkzdScdF8FZqMQ1cAAMBYxFgAQKcTG+WUJ29kvW0V3hr/SE7B3OF1RmwqvL5GR4FgFoIOAKDTcTgczTrkFBcdyaEpm+PQFQAAMBZBJwQqvDXKmL1JGbM3NXlFAACg47Tl87mtn+2heu5QvuaOYNugw6SeAACYz7ZBh0k9AQAwn22DDgAAMB9BBwAAGIugAwAAjEXQAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdAABgLIIOAAAwlu2ndLUsS5JUXl7eYc9Z4a1RbVWF/3lrvjGzbmNtAICmNfU52l6fwW39/A7Vc4fyNbfFhe/tC9/jDXFYTW1huGPHjik1NTXUZQAAgFY4evSoevXq1WC77YNObW2tjh8/rssuu0wOhyPU5YRUeXm5UlNTdfToUSUkJIS6HGPRzx2Hvu4Y9HPHoJ8DWZalM2fOKCUlRRERDZ+JY/vjIhEREY0mQTtKSEjgTdQB6OeOQ193DPq5Y9DPf5eYmNjkNpyMDAAAjEXQAQAAxiLowC8mJkbPPvusYmJiQl2K0ejnjkNfdwz6uWPQz61j+5ORAQCAuRjRAQAAxiLoAAAAYxF0AACAsQg6AADAWAQdG/r888/1wAMPqHv37oqLi9OgQYNUWFjob7csS/PmzVNKSopiY2N1++23a//+/SGsuHOqqanR3Llz1bt3b8XGxqpPnz7Ky8tTbW2tfxv6uuV27typsWPHKiUlRQ6HQxs3bgxob06fVlVV6Uc/+pF69OihSy+9VOPGjdOxY8c68FWEv8b6ubq6Wk8++aSuvfZaXXrppUpJSdHEiRN1/PjxgH3Qz01r6u/5mx555BE5HA4tXbo0YD393DiCjs2cPn1aQ4cOVVRUlP77v/9bHo9Hv/zlL9WlSxf/Nr/4xS+0ZMkSLVu2TPn5+UpOTtadd96pM2fOhK7wTmjRokV66aWXtGzZMh04cEC/+MUvtHjxYv3mN7/xb0Nft9y5c+c0cOBALVu2rN725vTpzJkztWHDBq1du1Z/+tOfdPbsWY0ZM0Y+n6+jXkbYa6yfKyoqtHv3bv3sZz/T7t27tX79eh08eFDjxo0L2I5+blpTf88XbNy4UR9++KFSUlLqtNHPTbBgK08++aR18803N9heW1trJScnWwsXLvSvO3/+vJWYmGi99NJLHVGiMUaPHm1Nnjw5YN3dd99tPfDAA5Zl0dfBIMnasGGDf7k5ffq3v/3NioqKstauXevf5vPPP7ciIiKsLVu2dFjtncnF/Vyfjz76yJJkffbZZ5Zl0c+t0VA/Hzt2zLriiiusffv2Wenp6daLL77ob6Ofm8aIjs28+eabysrK0g9+8ANdfvnluv7667Vq1Sp/e3FxsU6cOKGcnBz/upiYGN12223atWtXKErutG6++Wbt2LFDBw8elCR9/PHH+tOf/qRRo0ZJoq/bQ3P6tLCwUNXV1QHbpKSkaMCAAfR7G5SVlcnhcPhHh+nn4KitrdWECRP0xBNPqH///nXa6eem2X5ST7v53//9X61YsUKzZs3SU089pY8++kgzZsxQTEyMJk6cqBMnTkiSkpKSAh6XlJSkzz77LBQld1pPPvmkysrKdM0118jpdMrn8+m5557T/fffL0n0dTtoTp+eOHFC0dHR6tq1a51tLjweLXP+/HnNnj1bP/zhD/2TTdLPwbFo0SJFRkZqxowZ9bbTz00j6NhMbW2tsrKy9Pzzz0uSrr/+eu3fv18rVqzQxIkT/ds5HI6Ax1mWVWcdGrdu3Tq99tprev3119W/f38VFRVp5syZSklJ0aRJk/zb0dfB15o+pd9bp7q6Wvfdd59qa2u1fPnyJrenn5uvsLBQv/rVr7R79+4W9xn9/HccurKZnj17KjMzM2Bdv379VFJSIklKTk6WpDr/Ezh58mSd/yWjcU888YRmz56t++67T9dee60mTJigxx9/XAsWLJBEX7eH5vRpcnKyvF6vTp8+3eA2aJ7q6mrde++9Ki4u1rZt2/yjORL9HAzvvfeeTp48qbS0NEVGRioyMlKfffaZfvzjHysjI0MS/dwcBB2bGTp0qD799NOAdQcPHlR6erokqXfv3kpOTta2bdv87V6vV//zP/+jIUOGdGitnV1FRYUiIgLfYk6n0395OX0dfM3p08GDBysqKipgm9LSUu3bt49+b4ELIefQoUPavn27unfvHtBOP7fdhAkT9Mknn6ioqMj/k5KSoieeeEJbt26VRD83SyjPhEbH++ijj6zIyEjrueeesw4dOmT927/9mxUXF2e99tpr/m0WLlxoJSYmWuvXr7f27t1r3X///VbPnj2t8vLyEFbe+UyaNMm64oorrP/6r/+yiouLrfXr11s9evSwfvrTn/q3oa9b7syZM9aePXusPXv2WJKsJUuWWHv27PFf7dOcPp02bZrVq1cva/v27dbu3butO+64wxo4cKBVU1MTqpcVdhrr5+rqamvcuHFWr169rKKiIqu0tNT/U1VV5d8H/dy0pv6eL3bxVVeWRT83haBjQ2+99ZY1YMAAKyYmxrrmmmuslStXBrTX1tZazz77rJWcnGzFxMRYt956q7V3794QVdt5lZeXW4899piVlpZmXXLJJVafPn2sp59+OuCLgL5uuXfeeceSVOdn0qRJlmU1r08rKyutRx991OrWrZsVGxtrjRkzxiopKQnBqwlfjfVzcXFxvW2SrHfeece/D/q5aU39PV+svqBDPzfOYVmW1XHjRwAAAB2Hc3QAAICxCDoAAMBYBB0AAGAsgg4AADAWQQcAABiLoAMAAIxF0AEAAMYi6AAAAGMRdAAYbd68eXI4HHI4HFq6dGnI6jhy5Ii/jkGDBoWsDsBuCDoA2uzBBx+Uw+HQtGnT6rRNnz5dDodDDz74YMcX9v/1799fpaWl+qd/+qeA9Xv27NH48ePVs2dPxcTEKD09XWPGjNFbb72l5t40fuzYsRoxYkS9be+//74cDod2796t1NRUlZaW6sc//nGbXw+A5iPoAAiK1NRUrV27VpWVlf5158+f17//+78rLS0thJVJkZGRSk5OVlxcnH/dH/7wB/3DP/yDzp49q1dffVUej0dvvPGG7rrrLs2dO1dlZWXN2veUKVP09ttv67PPPqvTtnr1ag0aNEg33HCDnE6nkpOTFR8fH7TXBaBpBB0AQXHDDTcoLS1N69ev969bv369UlNTdf311wdsu2XLFt18883q0qWLunfvrjFjxuivf/2rv93r9erRRx9Vz549dckllygjI0MLFizwt8+bN09paWmKiYlRSkqKZsyY0aJaz507pylTpmj06NHatGmTcnJy9O1vf1s33nijpk6dqo8//liJiYn+7T0ej0aNGqX4+HglJSVpwoQJ+vLLLyVJY8aM0eWXX641a9YEPEdFRYXWrVunKVOmtKg2AMFF0AEQNA899JBeeeUV//Lq1as1efLkOtudO3dOs2bNUn5+vnbs2KGIiAh973vfU21trSTp17/+td588039x3/8hz799FO99tprysjIkCT953/+p1588UW9/PLLOnTokDZu3Khrr722RXX+8Y9/1KlTp/TTn/60wW0cDockqbS0VLfddpsGDRqkgoICbdmyRV988YXuvfdeSV+PFk2cOFFr1qwJONz1xhtvyOv16h//8R9bVBuA4IoMdQEAzDFhwgTNmTPHf+Ltn//8Z61du1bvvvtuwHb33HNPwPJvf/tbXX755fJ4PBowYIBKSkp01VVX6eabb5bD4VB6erp/25KSEiUnJ2vEiBGKiopSWlqabrzxxhbVefDgQUlS3759/evy8/M1bNgw//LatWs1ZswYrVixQjfccIOef/55f9vq1auVmpqqgwcP6uqrr9bkyZO1ePFivfvuu/59rF69Wnfffbe6du3aotoABBcjOgCCpkePHho9erReffVVvfLKKxo9erR69OhRZ7u//vWv+uEPf6g+ffooISFBvXv3lvR1iJG+Prm5qKhIffv21YwZM/THP/7R/9gf/OAHqqysVJ8+ffTwww9rw4YNqqmpaXPt1113nYqKilRUVKRz587591lYWKh33nlH8fHx/p9rrrnG/zok6ZprrtGQIUO0evVq//r33nuv3tEsAB2LoAMgqCZPnqw1a9bo1VdfbfCLfuzYsTp16pRWrVqlDz/8UB9++KGkr8/Nkb4+36e4uFj/8i//osrKSt177736/ve/L+nrk54//fRTud1uxcbGavr06br11ltVXV3d7BqvuuoqSdKnn37qXxcTE6Mrr7xSV155ZcC2tbW1Gjt2rD8EXfg5dOiQbr31Vv92U6ZM0e9//3uVl5frlVdeUXp6uoYPH97smgC0D4IOgKD67ne/K6/XK6/Xq5EjR9ZpP3XqlA4cOKC5c+dq+PDh6tevn06fPl1nu4SEBI0fP16rVq3SunXr9Pvf/15fffWVJCk2Nlbjxo3Tr3/9a7377rt6//33tXfv3mbXmJOTo27dumnRokVNbnvDDTdo//79ysjI8AehCz+XXnqpf7t7771XTqdTr7/+ul599VU99NBD/vN8AIQO5+gACCqn06kDBw74f79Y165d1b17d61cuVI9e/ZUSUmJZs+eHbDNiy++qJ49e2rQoEGKiIjQG2+8oeTkZHXp0kVr1qyRz+fTd77zHcXFxel3v/udYmNjA87jaUp8fLz+9V//VePHj9fo0aM1Y8YMXXXVVTp79qy2bNkSULvL5dKqVat0//3364knnlCPHj10+PBhrV27VqtWrfJvFx8fr/Hjx+upp55SWVlZSO8bBODvGNEBEHQJCQlKSEioty0iIkJr165VYWGhBgwYoMcff1yLFy8O2CY+Pl6LFi1SVlaWsrOzdeTIEW3evFkRERHq0qWLVq1apaFDh+q6667Tjh079NZbb6l79+4tqvF73/uedu3apbi4OE2cOFF9+/bVHXfcobffftt/IrIkpaSk6M9//rN8Pp9GjhypAQMG6LHHHlNiYqIiIgI/QqdMmaLTp09rxIgRIb93EICvOazm3v4TADqhefPmaePGjSoqKgp1KZLCrx7AdIzoADDe3r17FR8fr+XLl4eshpKSEsXHxwdcpg6g/TGiA8BoX331lf8k5m9961sBdzzuSDU1NTpy5Iikr6/wSk1NDUkdgN0QdAAAgLE4dAUAAIxF0AEAAMYi6AAAAGMRdAAAgLEIOgAAwFgEHQAAYCyCDgAAMBZBBwAAGOv/AdBFlDqFjIIPAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# continue with modified lhe, rivet, pythia ,...\n",
    "plot(\"filtered.lhe.gz\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "67370f77-f471-4901-8daf-9f8979fdb4bf",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}