Very few events survive after pythia8 showering

Hi,

Using "all" here within your syntax means that you generate a lot of process with very different kind of phsycis/ double counting/....
This is going to be a mess and I'm not even sure that even the diagram generated is going to be what you would expect.

Additionally, MLM/CKKW will only remove the double counting related to QCD.
But the default setup of pythia8_card should be wrong for that case, since you only have QCD
radiation from the "3all" and therefore have to specify the lowest number of QCD jet in the lowest sample.
Which will not be detected automatically here.

As far as I know we/Pythia8 does not have the equivalent of MLM for photon radiation. How where you planning to handle that issue?

Your syntax has also some other weird double counting/symmetry issue
indeed I would expect the first line to contain
e+ e- > z z, z > mu+ mu-
which is a four muon final state
but you also have
e+ e- > z mu+ mu-, z > mu+ mu-
Which has also four muon and which also includes the double onshell Z but where we do have break the symmetry between the muon and therefore you need to implement a cut that forbid to have the double onshell contribution here.

So this is likely not going to fly correctly with such syntax. You need to split this into smaller simulation and carefully be sure that everything is under control.

Cheers,

Olivier

> On 28 Oct 2022, at 03:15, Nikhil Bachhawat <email address hidden> wrote:
>
> New question #703625 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/703625
>
> Hi,
>
> I am trying to generate
>
> e+ e- > z all, z > mu+ mu-
> add process e+ e- > z all all, z > mu+ mu-
> add process e+ e- > z all all all, z > mu+ mu-
> add process e+ e- > z all all all all, z > mu+ mu-
> Switch on pythia8
> 0
> 2
> 300 GeV e+ and 300 GeV e- beams
> 300000 events
>
> and then i run it.
>
> I get this error at the end
>
> Command "shower pythia8 run_01 -f" interrupted with error:
> InvalidCmd : Pythia8 shower interrupted with return code 139.
> You can find more information in this log file:
> /mnt/e/mg5_amc_V3_4_1/eezjets_test2/Events/run_01/tag_2_pythia8.log
>
> I have read a bunch of the other forum posts and from what I gather I have a segmentation error
> in Pythia8. I tried multi-core and single-core (the error above is from single-core run) runs but it
> didn't help. Multi core is worse. I have to manually run the shower in -i mode and then i get some
> events after pythia showering. In single core I get more events in the end after Pythia8
> showering. I also tried reducing xqcut in the run card to 10 but it didn't help either.
>
> In the best case I get ~2k out of 100k events requested in the final output file. I gather from
> another thread that after removing double counted events, if 2k events are left then this
> error message would make sense? Is this the only problem or could there be a bug in pythia8
> or a parameter I should change?
>
> I have attached the pythia log file and my run card below.
>
> Thank you very much for your help,
> Nikhil
>
> #*********************************************************************
> # MadGraph5_aMC@NLO *
> # *
> # run_card.dat MadEvent *
> # *
> # This file is used to set the parameters of the run. *
> # *
> # Some notation/conventions: *
> # *
> # Lines starting with a '# ' are info or comments *
> # *
> # mind the format: value = variable ! comment *
> # *
> # To display more options, you can type the command: *
> # update to_full *
> #*********************************************************************
> #
> #*********************************************************************
> # Tag name for the run (one word) *
> #*********************************************************************
> tag_3 = run_tag ! name of the run
> #*********************************************************************
> # Number of events and rnd seed *
> # Warning: Do not generate more than 1M events in a single run *
> #*********************************************************************
> 300000 = nevents ! Number of unweighted events requested
> 0 = iseed ! rnd seed (0=assigned automatically=default))
> #*********************************************************************
> # Collider type and energy *
> # lpp: 0=No PDF, 1=proton, -1=antiproton, 2=elastic photon of proton,*
> # +/-3=PDF of electron/positron beam *
> # +/-4=PDF of muon/antimuon beam *
> #*********************************************************************
> 0 = lpp1 ! beam 1 type
> 0 = lpp2 ! beam 2 type
> 300.0 = ebeam1 ! beam 1 total energy in GeV
> 300.0 = ebeam2 ! beam 2 total energy in GeV
> #*********************************************************************
> # Beam polarization from -100 (left-handed) to 100 (right-handed) *
> #*********************************************************************
> 0.0 = polbeam1 ! beam polarization for beam 1
> 0.0 = polbeam2 ! beam polarization for beam 2
>
>
> #*********************************************************************
> # PDF CHOICE: this automatically fixes alpha_s and its evol. *
> # pdlabel: lhapdf=LHAPDF (installation needed) [1412.7420] *
> # iww=Improved Weizsaecker-Williams Approx.[hep-ph/9310350] *
> # eva=Effective W/Z/A Approx. [21yy.zzzzz] *
> # none=No PDF, same as lhapdf with lppx=0 *
> #*********************************************************************
> none = pdlabel1 ! PDF type for beam #1
> none = pdlabel2 ! PDF type for beam #2230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
> # To see heavy ion options: type "update ion_pdf"
> #*********************************************************************
> # Renormalization and factorization scales *
> #*********************************************************************
> False = fixed_ren_scale ! if .true. use fixed ren scale
> False = fixed_fac_scale ! if .true. use fixed fac scale
> 91.188 = scale ! fixed ren scale
> 91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
> 91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
> -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
> 1.0 = scalefact ! scale factor for event-by-event scales
>
>
> #*********************************************************************
> # Type and output format
> #*********************************************************************
> False = gridpack !True = setting up the grid pack
> -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
> average = event_norm ! average/sum. Normalization of the weight in the LHEF
> #*********************************************************************
> # Matching parameter (MLM only)
> #*********************************************************************
> 1 = ickkw ! 0 no matching, 1 MLM
> 1.0 = alpsfact ! scale factor for QCD emission vx
> False = chcluster ! cluster only according to channel diag
> 5 = asrwgtflavor ! highest quark flavor for a_s reweight
> True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
> ! (turn off for VBF and single top processes)
> 30.0 = xqcut ! minimum kt jet measure between partons
>
> #***********************************************************************
> # Turn on either the ktdurham or ptlund cut to activate *
> # CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
> #***********************************************************************
> -1.0 = ktdurham
> 0.4 = dparameter
> -1.0 = ptlund
> 1, 2, 3, 4, 5, 6, 21 = pdgs_for_merging_cut ! PDGs for two cuts above
>
> #*********************************************************************
> #
> #*********************************************************************
> # Phase-Space Optimization strategy (basic options)
> #*********************************************************************
> 0 = nhel ! using helicities importance sampling or not.
> ! 0: sum over helicity, 1: importance sampling
> 1 = sde_strategy ! default integration strategy (hep-ph/2021.00773)
> ! 1 is old strategy (using amp square)
> ! 2 is new strategy (using only the denominator)
> # To see advanced option for Phase-Space optimization: type "update psoptim"
> #*********************************************************************
> # Generation bias, check the wiki page below for more information: *
> # 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
> #*********************************************************************
> None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
> {} = bias_parameters ! Specifies the parameters of the module.
> #
> #*******************************
> # Parton level cuts definition *
> #*******************************
> 0.0 = dsqrt_shat ! minimal shat for full process
> #
> #
> #*********************************************************************
> # BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
> #*********************************************************************
> 15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
> #*********************************************************************
> # Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
> # (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
> #*********************************************************************
> False = cut_decays ! Cut decay products
> #*********************************************************************
> # Standard Cuts *
> #*********************************************************************
> # Minimum and maximum pt's (for max, -1 means no cut) *
> #*********************************************************************
> 20.0 = ptj ! minimum pt for the jets
> 0.0 = ptb ! minimum pt for the b
> 10.0 = pta ! minimum pt for the photons
> 10.0 = ptl ! minimum pt for the charged leptons
> 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
> -1.0 = ptjmax ! maximum pt for the jets
> -1.0 = ptbmax ! maximum pt for the b
> -1.0 = ptamax ! maximum pt for the photons
> -1.0 = ptlmax ! maximum pt for the charged leptons
> -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
> {} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
> #*********************************************************************
> # Minimum and maximum E's (in the center of mass frame) *
> #*********************************************************************
> 0.0 = ej ! minimum E for the jets
> 0.0 = eb ! minimum E for the b
> 0.0 = ea ! minimum E for the photons
> 0.0 = el ! minimum E for the charged leptons
> -1.0 = ejmax ! maximum E for the jets
> -1.0 = ebmax ! maximum E for the b
> -1.0 = eamax ! maximum E for the photons
> -1.0 = elmax ! maximum E for the charged leptons
> {} = e_min_pdg ! E cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = e_max_pdg ! E cut for other particles (syntax e.g. {6: 100, 25: 50})
>
> #*********************************************************************
> # Maximum and minimum absolute rapidity (for max, -1 means no cut) *
> #*********************************************************************
> 5.0 = etaj ! max rap for the jets
> -1.0 = etab ! max rap for the b
> 2.5 = etaa ! max rap for the photons
> 2.5 = etal ! max rap for the charged leptons
> 0.0 = etajmin ! min rap for the jets
> 0.0 = etabmin ! min rap for the b
> 0.0 = etaamin ! min rap for the photons
> 0.0 = etalmin ! main rap for the charged leptons
> {} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
> {} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
> #*********************************************************************
> # Minimum and maximum DeltaR distance *
> #*********************************************************************
> 0.0 = drjj ! min distance between jets
> 0.0 = drbb ! min distance between b's
> 0.4 = drll ! min distance between leptons
> 0.4 = draa ! min distance between gammas
> 0.0 = drbj ! min distance between b and jet
> 0.4 = draj ! min distance between gamma and jet
> 0.0 = drjl ! min distance between jet and lepton
> 0.0 = drab ! min distance between gamma and b
> 0.0 = drbl ! min distance between b and lepton
> 0.4 = dral ! min distance between gamma and lepton
> -1.0 = drjjmax ! max distance between jets
> -1.0 = drbbmax ! max distance between b's
> -1.0 = drllmax ! max distance between leptons
> -1.0 = draamax ! max distance between gammas
> -1.0 = drbjmax ! max distance between b and jet
> -1.0 = drajmax ! max distance between gamma and jet
> -1.0 = drjlmax ! max distance between jet and lepton
> -1.0 = drabmax ! max distance between gamma and b
> -1.0 = drblmax ! max distance between b and lepton
> -1.0 = dralmax ! maxdistance between gamma and lepton
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> # WARNING: for four lepton final state mmll cut require to have *
> # different lepton masses for each flavor! *
> #*********************************************************************
> 0.0 = mmjj ! min invariant mass of a jet pair
> 0.0 = mmbb ! min invariant mass of a b pair
> 0.0 = mmaa ! min invariant mass of gamma gamma pair
> 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> -1.0 = mmjjmax ! max invariant mass of a jet pair
> -1.0 = mmbbmax ! max invariant mass of a b pair
> -1.0 = mmaamax ! max invariant mass of gamma gamma pair
> -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> {} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
> {'default': False} = mxx_only_part_antipart ! if True the invariant mass is applied only
> ! to pairs of particle/antiparticle and not to pairs of the same pdg codes.
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
> -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
> #*********************************************************************
> # Minimum and maximum pt for 4-momenta sum of leptons / neutrino *
> # for pair of lepton includes only same flavor, opposite charge
> #*********************************************************************
> 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
> -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
> #*********************************************************************
> # Inclusive cuts *
> #*********************************************************************
> 0.0 = ptheavy ! minimum pt for at least one heavy final state
> 0.0 = xptj ! minimum pt for at least one jet
> 0.0 = xptb ! minimum pt for at least one b
> 0.0 = xpta ! minimum pt for at least one photon
> 0.0 = xptl ! minimum pt for at least one charged lepton
> #*********************************************************************
> # Control the pt's of the jets sorted by pt *
> #*********************************************************************
> 0.0 = ptj1min ! minimum pt for the leading jet in pt
> 0.0 = ptj2min ! minimum pt for the second jet in pt
> 0.0 = ptj3min ! minimum pt for the third jet in pt
> 0.0 = ptj4min ! minimum pt for the fourth jet in pt
> -1.0 = ptj1max ! maximum pt for the leading jet in pt
> -1.0 = ptj2max ! maximum pt for the second jet in pt
> -1.0 = ptj3max ! maximum pt for the third jet in pt
> -1.0 = ptj4max ! maximum pt for the fourth jet in pt
> 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
> #*********************************************************************
> # Control the pt's of leptons sorted by pt *
> #*********************************************************************
> 0.0 = ptl1min ! minimum pt for the leading lepton in pt
> 0.0 = ptl2min ! minimum pt for the second lepton in pt
> 0.0 = ptl3min ! minimum pt for the third lepton in pt
> 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
> -1.0 = ptl1max ! maximum pt for the leading lepton in pt
> -1.0 = ptl2max ! maximum pt for the second lepton in pt
> -1.0 = ptl3max ! maximum pt for the third lepton in pt
> -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
> #*********************************************************************
> # Control the Ht(k)=Sum of k leading jets *
> #*********************************************************************
> 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
> -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
> 0.0 = ihtmin !inclusive Ht for all partons (including b)
> -1.0 = ihtmax !inclusive Ht for all partons (including b)
> 0.0 = ht2min ! minimum Ht for the two leading jets
> 0.0 = ht3min ! minimum Ht for the three leading jets
> 0.0 = ht4min ! minimum Ht for the four leading jets
> -1.0 = ht2max ! maximum Ht for the two leading jets
> -1.0 = ht3max ! maximum Ht for the three leading jets
> -1.0 = ht4max ! maximum Ht for the four leading jets
> #***********************************************************************
> # Photon-isolation cuts, according to hep-ph/9801442 *
> # When ptgmin=0, all the other parameters are ignored *
> # When ptgmin>0, pta and draj are not going to be used *
> #***********************************************************************
> 0.0 = ptgmin ! Min photon transverse momentum
> 0.4 = r0gamma ! Radius of isolation code
> 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> True = isoem ! isolate photons from EM energy (photons and leptons)
> #*********************************************************************
> # WBF cuts *
> #*********************************************************************
> 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
> 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
> #*********************************************************************
> # maximal pdg code for quark to be considered as a light jet *
> # (otherwise b cuts are applied) *
> #*********************************************************************
> 4 = maxjetflavor ! Maximum jet pdg code
> #*********************************************************************
> #
> #*********************************************************************
> # Store info for systematics studies *
> # WARNING: Do not use for interference type of computation *
> #*********************************************************************
> False = use_syst ! Enable systematics studies
> #
> systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
> ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset', '--alps=0.5,1,2'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
> #*********************************************************************
> # Additional hidden parameters
> #*********************************************************************
> 230000 = lhaid # hidden_parameter
>
>
>
>
>
>
> *------------------------------------------------------------------------------------*
> | |
> | *------------------------------------------------------------------------------* |
> | | | |
> | | | |
> | | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | |
> | | P P Y Y T H H I A A This is PYTHIA version 8.306 | |
> | | PPP Y T HHHHH I AAAAA Last date of change: 28 Jun 2021 | |
> | | P Y T H H I A A | |
> | | P Y T H H III A A Now is 27 Oct 2022 at 19:18:30 | |
> | | | |
> | | Program documentation and an archive of historic versions is found on: | |
> | | | |
> | | https://pythia.org/ | |
> | | | |
> | | PYTHIA is authored by a collaboration consisting of: | |
> | | | |
> | | Christian Bierlich, Nishita Desai, Leif Gellersen, Ilkka Helenius, Philip | |
> | | Ilten, Leif Lönnblad, Stephen Mrenna, Stefan Prestel, Christian Preuss, | |
> | | Torbjörn Sjöstrand, Peter Skands, Marius Utheim and Rob Verheyen. | |
> | | | |
> | | The complete list of authors, including contact information and | |
> | | affiliations, can be found on https://pythia.org/. | |
> | | Problems or bugs should be reported on email at <email address hidden>. | |
> | | | |
> | | The main program reference is 'An Introduction to PYTHIA 8.2', | |
> | | T. Sjöstrand et al, Comput. Phys. Commun. 191 (2015) 159 | |
> | | [arXiv:1410.3012 [hep-ph]] | |
> | | | |
> | | The main physics reference is the 'PYTHIA 6.4 Physics and Manual', | |
> | | T. Sjöstrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 [hep-ph/0603175] | |
> | | | |
> | | PYTHIA is released under the GNU General Public Licence version 2 or later.| |
> | | Please respect the MCnet Guidelines for Event Generator Authors and Users. | |
> | | | |
> | | Disclaimer: this program comes without any guarantees. | |
> | | Beware of errors and use common sense when interpreting results. | |
> | | | |
> | | Copyright (C) 2021 Torbjörn Sjöstrand | |
> | | | |
> | | | |
> | *------------------------------------------------------------------------------* |
> | |
> *------------------------------------------------------------------------------------*
>
>
> *-------- Madgraph parameters --------*
> *---------------------------------------*
>
> *----- Madgraph matching parameters -----*
> | qCut | 45.000 |
> | nQmatch | 4 |
> | clFact | 1.000 |
> | Jet algorithm | kT |
> | eTjetMin | 20.000 |
> | etaJetMax | 1000.000 |
> | jetAllow | 1 |
> | Mode | exclusive |
> *-----------------------------------------*
> *----------------------- SusyLesHouches SUSY/BSM Interface ------------------------*
> | Last Change 12 Apr 2017 - P. Skands
> | Parsing: unweighted_events.lhe.gz
> | (SLHA::readFile) line 473 - storing non-SLHA(2) block: yukawa
> *------------------------------------------------------------------------------------*
> PYTHIA Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring MASS entries for id = {1,2,3,4,5,11,12,13,14,15,16,21,22,23,24} (m0 < SLHA:minMassSM)
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring DECAY tables for id = {23,24,1,2,3,4,5,11,12,13,14,15,16,21,22} (m0 < SLHA:minMassSM)
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables for id = {6,25} (total width provided but no Branching Ratios)
>
> *------- PYTHIA Process Initialization --------------------------*
> | |
> | We collide e+ with e- at a CM energy of 6.000e+02 GeV |
> | |
> |------------------------------------------------------------------|
> | | |
> | Subprocess Code | Estimated |
> | | max (mb) |
> | | |
> |------------------------------------------------------------------|
> | | |
> | Les Houches User Process(es) 9999 | 4.000e-09 |
> | |
> *------- End PYTHIA Process Initialization -----------------------*
>
> *------- PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings (changes only) ------------------*
> | |
> | Name | Now | Default Min Max |
> | | | |
> | Beams:frameType | 4 | 1 1 5 |
> | Beams:LHEF | unweighted_events.lhe.gz | void |
> | Beams:setProductionScalesFromLHEF | on | off |
> | Check:epTolErr | 0.0100000 | 1.0000e-04 |
> | HEPMCoutput:file | tag_1_pythia8_events.hepmc | void |
> | HEPMCoutput:scaling | 1.0000e+09 | 1.00000 0.0 |
> | JetMatching:coneRadius | 1.00000 | 0.70000 0.10000 |
> | JetMatching:etaJetMax | 1000.000 | 2.50000 0.10000 |
> | JetMatching:merge | on | off |
> | JetMatching:nJetMax | 2 | -1 -1 |
> | JetMatching:nQmatch | 4 | 5 3 6 |
> | JetMatching:qCut | 45.00000 | 10.00000 0.0 |
> | JetMatching:setMad | off | on |
> | LHEFInputs:nSubruns | 1 | 0 0 |
> | Main:numberOfEvents | 0 | 1000 0 |
> | Main:subrun | 0 | -999 0 |
> | SysCalc:qWeed | 30.00000 | 0.0 0.0 |
> | |
> *------- End PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings -----------------------------*
>
> -------- PYTHIA Particle Data Table (changed only) ------------------------------------------------------------------------------
>
> id name antiName spn chg col m0 mWidth mMin mMax tau0 res dec ext vis wid
> no onMode bRatio meMode products
>
> 6 t tbar 2 2 1 171.00000 1.30203 163.54250 178.45750 1.51554e-13 1 1 0 1 0
> 0 1 0.0000664 0 24 1
> 1 1 0.0017362 0 24 3
> 2 1 0.9981973 0 24 5
> 3 0 0.0000000 0 37 5
>
> 25 h0 1 0 0 125.00000 0.00408 124.96809 125.03191 4.83238e-11 1 1 0 0 0
> 0 1 0.0000009 0 1 -1
> 1 1 0.0000002 0 2 -2
> 2 1 0.0002163 0 3 -3
> 3 1 0.0288493 0 4 -4
> 4 1 0.5768873 0 5 -5
> 5 1 0.0000000 0 6 -6
> 6 1 0.0000000 0 11 -11
> 7 1 0.0002185 0 13 -13
> 8 1 0.0628569 0 15 -15
> 9 1 0.0850383 0 21 21
> 10 1 0.0022965 0 22 22
> 11 1 0.0016155 0 22 23
> 12 1 0.0261758 0 23 23
> 13 1 0.2158445 0 24 -24
> 14 1 0.0000000 103 1000022 1000022
> 15 1 0.0000000 103 1000023 1000022
> 16 1 0.0000000 103 1000023 1000023
> 17 1 0.0000000 103 1000025 1000022
> 18 1 0.0000000 103 1000025 1000023
> 19 1 0.0000000 103 1000025 1000025
> 20 1 0.0000000 103 1000035 1000022
> 21 1 0.0000000 103 1000035 1000023
> 22 1 0.0000000 103 1000035 1000025
> 23 1 0.0000000 103 1000035 1000035
> 24 1 0.0000000 103 1000024 -1000024
> 25 1 0.0000000 103 1000024 -1000037
> 26 1 0.0000000 103 1000037 -1000024
> 27 1 0.0000000 103 1000037 -1000037
> 28 1 0.0000000 103 1000001 -1000001
> 29 1 0.0000000 103 2000001 -2000001
> 30 1 0.0000000 103 1000001 -2000001
> 31 1 0.0000000 103 -1000001 2000001
> 32 1 0.0000000 103 1000002 -1000002
> 33 1 0.0000000 103 2000002 -2000002
> 34 1 0.0000000 103 1000002 -2000002
> 35 1 0.0000000 103 -1000002 2000002
> 36 1 0.0000000 103 1000003 -1000003
> 37 1 0.0000000 103 2000003 -2000003
> 38 1 0.0000000 103 1000003 -2000003
> 39 1 0.0000000 103 -1000003 2000003
> 40 1 0.0000000 103 1000004 -1000004
> 41 1 0.0000000 103 2000004 -2000004
> 42 1 0.0000000 103 1000004 -2000004
> 43 1 0.0000000 103 -1000004 2000004
> 44 1 0.0000000 103 1000005 -1000005
> 45 1 0.0000000 103 2000005 -2000005
> 46 1 0.0000000 103 1000005 -2000005
> 47 1 0.0000000 103 -1000005 2000005
> 48 1 0.0000000 103 1000006 -1000006
> 49 1 0.0000000 103 2000006 -2000006
> 50 1 0.0000000 103 1000006 -2000006
> 51 1 0.0000000 103 -1000006 2000006
> 52 1 0.0000000 103 1000011 -1000011
> 53 1 0.0000000 103 2000011 -2000011
> 54 1 0.0000000 103 1000011 -2000011
> 55 1 0.0000000 103 -1000011 2000011
> 56 1 0.0000000 103 1000012 -1000012
> 57 1 0.0000000 103 2000012 -2000012
> 58 1 0.0000000 103 1000012 -2000012
> 59 1 0.0000000 103 -1000012 2000012
> 60 1 0.0000000 103 1000013 -1000013
> 61 1 0.0000000 103 2000013 -2000013
> 62 1 0.0000000 103 1000013 -2000013
> 63 1 0.0000000 103 -1000013 2000013
> 64 1 0.0000000 103 1000014 -1000014
> 65 1 0.0000000 103 2000014 -2000014
> 66 1 0.0000000 103 1000014 -2000014
> 67 1 0.0000000 103 -1000014 2000014
> 68 1 0.0000000 103 1000015 -1000015
> 69 1 0.0000000 103 2000015 -2000015
> 70 1 0.0000000 103 1000015 -2000015
> 71 1 0.0000000 103 -1000015 2000015
> 72 1 0.0000000 103 1000016 -1000016
> 73 1 0.0000000 103 2000016 -2000016
> 74 1 0.0000000 103 1000016 -2000016
> 75 1 0.0000000 103 -1000016 2000016
>
> -------- End PYTHIA Particle Data Table -----------------------------------------------------------------------------------------
>
> Counting events in input LHEF, please wait.
> Found 300000 events in input LHEF.
>
>
>
> Start generating events
> PYTHIA Info from Pythia::init: be aware that successive calls to init() do not clear previous settings.
>
> *-------- Madgraph parameters --------*
> *---------------------------------------*
>
> *----- Madgraph matching parameters -----*
> | qCut | 45.0000000 |
> | nQmatch | 4 |
> | clFact | 1.0000000 |
> | Jet algorithm | kT |
> | eTjetMin | 20.0000000 |
> | etaJetMax | 1000.0000000 |
> | jetAllow | 1 |
> | Mode | exclusive |
> *-----------------------------------------*
> *----------------------- SusyLesHouches SUSY/BSM Interface ------------------------*
> | Last Change 12 Apr 2017 - P. Skands
> | Parsing: unweighted_events.lhe.gz
> | (SLHA::readFile) line 473 - storing non-SLHA(2) block: yukawa
> *------------------------------------------------------------------------------------*
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring MASS entries for id = {1,2,3,4,5,11,12,13,14,15,16,21,22,23,24} (m0 < SLHA:minMassSM)
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring DECAY tables for id = {23,24,1,2,3,4,5,11,12,13,14,15,16,21,22} (m0 < SLHA:minMassSM)
> PYTHIA Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables for id = {6,25} (total width provided but no Branching Ratios)
>
> *------- PYTHIA Process Initialization --------------------------*
> | |
> | We collide e+ with e- at a CM energy of 6.000e+02 GeV |
> | |
> |------------------------------------------------------------------|
> | | |
> | Subprocess Code | Estimated |
> | | max (mb) |
> | | |
> |------------------------------------------------------------------|
> | | |
> | Les Houches User Process(es) 9999 | 4.000e-09 |
> | |
> *------- End PYTHIA Process Initialization -----------------------*
>
> *------- PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings (changes only) ------------------*
> | |
> | Name | Now | Default Min Max |
> | | | |
> | Beams:frameType | 4 | 1 1 5 |
> | Beams:LHEF | unweighted_events.lhe.gz | void |
> | Beams:setProductionScalesFromLHEF | on | off |
> | Check:epTolErr | 0.0100000 | 1.0000e-04 |
> | HEPMCoutput:file | tag_1_pythia8_events.hepmc | void |
> | HEPMCoutput:scaling | 1.0000e+09 | 1.00000 0.0 |
> | JetMatching:coneRadius | 1.00000 | 0.70000 0.10000 |
> | JetMatching:etaJetMax | 1000.000 | 2.50000 0.10000 |
> | JetMatching:merge | on | off |
> | JetMatching:nJetMax | 2 | -1 -1 |
> | JetMatching:nQmatch | 4 | 5 3 6 |
> | JetMatching:qCut | 45.00000 | 10.00000 0.0 |
> | JetMatching:setMad | off | on |
> | LHEFInputs:nSubruns | 1 | 0 0 |
> | Main:numberOfEvents | 0 | 1000 0 |
> | Main:subrun | 0 | -999 0 |
> | SysCalc:qWeed | 30.00000 | 0.0 0.0 |
> | |
> *------- End PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings -----------------------------*
>
> -------- PYTHIA Particle Data Table (changed only) ------------------------------------------------------------------------------
>
> id name antiName spn chg col m0 mWidth mMin mMax tau0 res dec ext vis wid
> no onMode bRatio meMode products
>
> 6 t tbar 2 2 1 171.00000 1.30203 163.54250 178.45750 1.51554e-13 1 1 0 1 0
> 0 1 0.0000664 0 24 1
> 1 1 0.0017362 0 24 3
> 2 1 0.9981973 0 24 5
> 3 0 0.0000000 0 37 5
>
> 25 h0 1 0 0 125.00000 0.00408 124.96809 125.03191 4.83238e-11 1 1 0 0 0
> 0 1 0.0000009 0 1 -1
> 1 1 0.0000002 0 2 -2
> 2 1 0.0002163 0 3 -3
> 3 1 0.0288493 0 4 -4
> 4 1 0.5768873 0 5 -5
> 5 1 0.0000000 0 6 -6
> 6 1 0.0000000 0 11 -11
> 7 1 0.0002185 0 13 -13
> 8 1 0.0628569 0 15 -15
> 9 1 0.0850383 0 21 21
> 10 1 0.0022965 0 22 22
> 11 1 0.0016155 0 22 23
> 12 1 0.0261758 0 23 23
> 13 1 0.2158445 0 24 -24
> 14 1 0.0000000 103 1000022 1000022
> 15 1 0.0000000 103 1000023 1000022
> 16 1 0.0000000 103 1000023 1000023
> 17 1 0.0000000 103 1000025 1000022
> 18 1 0.0000000 103 1000025 1000023
> 19 1 0.0000000 103 1000025 1000025
> 20 1 0.0000000 103 1000035 1000022
> 21 1 0.0000000 103 1000035 1000023
> 22 1 0.0000000 103 1000035 1000025
> 23 1 0.0000000 103 1000035 1000035
> 24 1 0.0000000 103 1000024 -1000024
> 25 1 0.0000000 103 1000024 -1000037
> 26 1 0.0000000 103 1000037 -1000024
> 27 1 0.0000000 103 1000037 -1000037
> 28 1 0.0000000 103 1000001 -1000001
> 29 1 0.0000000 103 2000001 -2000001
> 30 1 0.0000000 103 1000001 -2000001
> 31 1 0.0000000 103 -1000001 2000001
> 32 1 0.0000000 103 1000002 -1000002
> 33 1 0.0000000 103 2000002 -2000002
> 34 1 0.0000000 103 1000002 -2000002
> 35 1 0.0000000 103 -1000002 2000002
> 36 1 0.0000000 103 1000003 -1000003
> 37 1 0.0000000 103 2000003 -2000003
> 38 1 0.0000000 103 1000003 -2000003
> 39 1 0.0000000 103 -1000003 2000003
> 40 1 0.0000000 103 1000004 -1000004
> 41 1 0.0000000 103 2000004 -2000004
> 42 1 0.0000000 103 1000004 -2000004
> 43 1 0.0000000 103 -1000004 2000004
> 44 1 0.0000000 103 1000005 -1000005
> 45 1 0.0000000 103 2000005 -2000005
> 46 1 0.0000000 103 1000005 -2000005
> 47 1 0.0000000 103 -1000005 2000005
> 48 1 0.0000000 103 1000006 -1000006
> 49 1 0.0000000 103 2000006 -2000006
> 50 1 0.0000000 103 1000006 -2000006
> 51 1 0.0000000 103 -1000006 2000006
> 52 1 0.0000000 103 1000011 -1000011
> 53 1 0.0000000 103 2000011 -2000011
> 54 1 0.0000000 103 1000011 -2000011
> 55 1 0.0000000 103 -1000011 2000011
> 56 1 0.0000000 103 1000012 -1000012
> 57 1 0.0000000 103 2000012 -2000012
> 58 1 0.0000000 103 1000012 -2000012
> 59 1 0.0000000 103 -1000012 2000012
> 60 1 0.0000000 103 1000013 -1000013
> 61 1 0.0000000 103 2000013 -2000013
> 62 1 0.0000000 103 1000013 -2000013
> 63 1 0.0000000 103 -1000013 2000013
> 64 1 0.0000000 103 1000014 -1000014
> 65 1 0.0000000 103 2000014 -2000014
> 66 1 0.0000000 103 1000014 -2000014
> 67 1 0.0000000 103 -1000014 2000014
> 68 1 0.0000000 103 1000015 -1000015
> 69 1 0.0000000 103 2000015 -2000015
> 70 1 0.0000000 103 1000015 -2000015
> 71 1 0.0000000 103 -1000015 2000015
> 72 1 0.0000000 103 1000016 -1000016
> 73 1 0.0000000 103 2000016 -2000016
> 74 1 0.0000000 103 1000016 -2000016
> 75 1 0.0000000 103 -1000016 2000016
>
> -------- End PYTHIA Particle Data Table -----------------------------------------------------------------------------------------
>
>
> -------- LHA event information and listing ----------------------------------------------------------------------
>
> process = 2 weight = 8.2822e-02 scale = 0.0000e+00 (GeV)
> alpha_em = 7.5468e-03 alpha_strong = 9.9989e-02
>
> Participating Particles
> no id stat mothers colours p_x p_y p_z e m tau spin
> 1 -11 -1 0 0 0 0 0.000 0.000 300.000 300.000 0.000 0.000 1.000
> 2 11 -1 0 0 0 0 -0.000 -0.000 -300.000 300.000 0.000 0.000 -1.000
> 3 23 2 1 2 0 0 62.316 -192.850 -202.949 302.415 95.875 0.000 0.000
> 4 23 2 1 2 0 0 -62.316 192.850 202.949 297.585 79.334 0.000 0.000
> 5 -13 1 4 4 0 0 -23.641 5.956 10.496 26.543 0.000 0.000 -1.000
> 6 13 1 4 4 0 0 -38.675 186.895 192.454 271.042 0.000 0.000 1.000
> 7 2 1 3 3 501 0 24.007 -177.333 -184.894 257.311 0.000 0.000 -1.000
> 8 -2 1 3 3 0 501 38.309 -15.518 -18.055 45.104 0.000 0.000 1.000
>
> -------- End LHA event information and listing ------------------------------------------------------------------
>
> -------- PYTHIA Info Listing ----------------------------------------
>
> Beam A: id = -11, pz = 3.000e+02, e = 3.000e+02, m = 5.110e-04.
> Beam B: id = 11, pz = -3.000e+02, e = 3.000e+02, m = 5.110e-04.
>
> In 1: id = -11, x = 1.000e+00, pdf = 0.000e+00 at Q2 = 0.000e+00.
> In 2: id = 11, x = 1.000e+00, pdf = 0.000e+00 at same Q2.
>
> Process Les Houches User Process(es) with code 9999 is 2 -> 2.
> Subprocess user process 2 with code 2 is 2 -> 2.
> It has sHat = 3.600e+05, tHat = -2.940e+05, uHat = -5.049e+04,
> pTHat = 2.027e+02, m3Hat = 9.587e+01, m4Hat = 7.933e+01,
> thetaHat = 2.357e+00, phiHat = -1.258e+00.
> alphaEM = 7.547e-03, alphaS = 9.999e-02 at Q2 = 0.000e+00.
>
> Impact parameter b = 0.000e+00 gives enhancement factor = 1.000e+00.
> Max pT scale for MPI = 6.000e+02, ISR = 6.000e+02, FSR = 6.000e+02.
> Number of MPI = 1, ISR = 0, FSRproc = 0, FSRreson = 12.
>
> -------- End PYTHIA Info Listing ------------------------------------
>
> -------- PYTHIA Event Listing (hard process) -----------------------------------------------------------------------------------
>
> no id name status mothers daughters colours p_x p_y p_z e m
> 0 90 (system) -11 0 0 0 0 0 0 0.000 0.000 0.000 600.000 600.000
> 1 -11 (e+) -12 0 0 3 0 0 0 0.000 0.000 300.000 300.000 0.001
> 2 11 (e-) -12 0 0 4 0 0 0 0.000 0.000 -300.000 300.000 0.001
> 3 -11 (e+) -21 1 0 5 6 0 0 0.000 0.000 300.000 300.000 0.000
> 4 11 (e-) -21 2 0 5 6 0 0 0.000 0.000 -300.000 300.000 0.000
> 5 23 (Z0) -22 3 4 7 8 0 0 62.316 -192.850 -202.949 302.415 95.875
> 6 23 (Z0) -22 3 4 9 10 0 0 -62.316 192.850 202.949 297.585 79.334
> 7 2 u 23 5 0 0 0 501 0 24.007 -177.333 -184.894 257.311 0.000
> 8 -2 ubar 23 5 0 0 0 0 501 38.309 -15.518 -18.055 45.104 0.000
> 9 -13 mu+ 23 6 0 0 0 0 0 -23.641 5.956 10.496 26.543 0.106
> 10 13 mu- 23 6 0 0 0 0 0 -38.675 186.894 192.453 271.042 0.106
> Charge sum: 0.000 Momentum sum: 0.000 0.000 -0.000 600.000 600.000
>
> -------- End PYTHIA Event Listing -----------------------------------------------------------------------------------------------
>
> -------- PYTHIA Event Listing (complete event) ---------------------------------------------------------------------------------
>
> no id name status mothers daughters colours p_x p_y p_z e m
> 0 90 (system) -11 0 0 0 0 0 0 0.000 0.000 0.000 600.000 600.000
> 1 -11 (e+) -12 0 0 3 0 0 0 0.000 0.000 300.000 300.000 0.001
> 2 11 (e-) -12 0 0 4 0 0 0 0.000 0.000 -300.000 300.000 0.001
> 3 -11 (e+) -21 1 0 5 6 0 0 0.000 0.000 300.000 300.000 0.000
> 4 11 (e-) -21 2 0 5 6 0 0 0.000 0.000 -300.000 300.000 0.000
> 5 23 (Z0) -22 3 4 7 8 0 0 62.316 -192.850 -202.949 302.415 95.875
> 6 23 (Z0) -22 3 4 42 43 0 0 -62.316 192.850 202.949 297.585 79.334
> 7 2 (u) -23 5 0 11 11 501 0 24.007 -177.333 -184.894 257.311 0.000
> 8 -2 (ubar) -23 5 0 9 10 0 501 38.309 -15.518 -18.055 45.104 0.000
> 9 -2 (ubar) -51 8 0 14 14 0 502 21.113 -3.387 -9.281 23.310 0.000
> 10 21 (g) -51 8 0 12 13 502 501 17.575 -14.931 -11.694 25.857 0.000
> 11 2 (u) -52 7 7 17 17 501 0 23.628 -174.533 -181.975 253.249 0.000
> 12 21 (g) -51 10 0 15 16 503 501 10.087 -5.785 -3.900 12.265 0.000
> 13 21 (g) -51 10 0 18 19 502 503 10.483 -9.626 -9.111 16.899 0.000
> 14 -2 (ubar) -52 9 9 26 26 0 502 18.118 -2.906 -7.964 20.003 0.000
> 15 21 (g) -51 12 0 20 20 503 504 8.604 -3.606 -2.876 9.762 0.000
> 16 21 (g) -51 12 0 21 22 504 501 2.120 -6.885 -5.930 9.331 0.000
> 17 2 (u) -52 11 11 29 29 501 0 22.991 -169.827 -177.069 246.421 0.000
> 18 21 (g) -51 13 0 24 25 502 505 7.838 -7.257 -7.779 13.214 0.000
> 19 21 (g) -51 13 0 38 38 505 503 6.799 -4.110 -2.720 8.397 0.000
> 20 21 (g) -52 15 15 23 23 503 504 4.450 -1.865 -1.488 5.049 0.000
> 21 21 (g) -51 16 0 27 28 506 501 1.952 -4.590 -5.181 7.192 0.000
> 22 21 (g) -51 16 0 52 52 504 506 0.781 -2.551 -0.953 2.833 0.000
> 23 21 (g) -52 20 20 53 53 503 504 3.838 -1.609 -1.283 4.355 0.000
> 24 21 (g) -51 18 0 36 37 507 505 3.260 -3.836 -3.908 6.373 0.000
> 25 21 (g) -51 18 0 32 32 502 507 5.160 -3.514 -4.128 7.484 0.000
> 26 -2 (ubar) -52 14 14 30 31 0 502 17.535 -2.813 -7.708 19.360 0.000
> 27 21 (g) -51 21 0 33 34 506 508 2.364 -6.888 -6.821 9.978 0.000
> 28 21 (g) -51 21 0 35 35 508 501 0.813 -6.756 -7.800 10.352 0.000
> 29 2 (u) -52 17 17 39 40 501 0 21.765 -160.773 -167.629 233.283 0.000
> 30 -2 (ubar) -51 26 0 59 59 0 509 16.393 -2.589 -7.751 18.317 0.000
> 31 21 (g) -51 26 0 58 58 509 502 1.674 -0.586 -0.383 1.814 0.000
> 32 21 (g) -52 25 25 57 57 502 507 4.628 -3.152 -3.702 6.713 0.000
> 33 21 (g) -51 27 0 51 51 506 510 0.092 -1.782 -2.222 2.849 0.000
> 34 21 (g) -51 27 0 50 50 510 508 2.681 -8.499 -8.516 12.327 0.000
> 35 21 (g) -52 28 28 41 41 508 501 0.405 -3.363 -3.883 5.153 0.000
> 36 21 (g) -51 24 0 56 56 507 511 2.776 -2.330 -2.998 4.703 0.000
> 37 21 (g) -51 24 0 55 55 511 505 1.379 -2.047 -1.268 2.775 0.000
> 38 21 (g) -52 19 19 54 54 505 503 5.904 -3.569 -2.362 7.292 0.000
> 39 2 (u) -51 29 0 47 47 512 0 15.817 -113.164 -118.026 164.276 0.000
> 40 21 (g) -51 29 0 48 48 501 512 6.101 -48.876 -51.065 70.949 0.000
> 41 21 (g) -52 35 35 49 49 508 501 0.252 -2.097 -2.421 3.212 0.000
> 42 -13 (mu+) -23 6 0 46 46 0 0 -23.641 5.956 10.496 26.543 0.106
> 43 13 (mu-) -23 6 0 44 45 0 0 -38.675 186.894 192.453 271.042 0.106
> 44 13 mu- 51 43 0 0 0 0 0 -38.674 186.894 192.453 271.041 0.106
> 45 22 gamma 51 43 0 0 0 0 0 -0.000 0.001 0.001 0.001 0.000
> 46 -13 mu+ 52 42 42 0 0 0 0 -23.641 5.956 10.496 26.543 0.106
> 47 2 (u) -71 39 39 60 82 512 0 15.817 -113.164 -118.026 164.276 0.000
> 48 21 (g) -71 40 40 60 82 501 512 6.101 -48.876 -51.065 70.949 0.000
> 49 21 (g) -71 41 41 60 82 508 501 0.252 -2.097 -2.421 3.212 0.000
> 50 21 (g) -71 34 34 60 82 510 508 2.681 -8.499 -8.516 12.327 0.000
> 51 21 (g) -71 33 33 60 82 506 510 0.092 -1.782 -2.222 2.849 0.000
> 52 21 (g) -71 22 22 60 82 504 506 0.781 -2.551 -0.953 2.833 0.000
> 53 21 (g) -71 23 23 60 82 503 504 3.838 -1.609 -1.283 4.355 0.000
> 54 21 (g) -71 38 38 60 82 505 503 5.904 -3.569 -2.362 7.292 0.000
> 55 21 (g) -71 37 37 60 82 511 505 1.379 -2.047 -1.268 2.775 0.000
> 56 21 (g) -71 36 36 60 82 507 511 2.776 -2.330 -2.998 4.703 0.000
> 57 21 (g) -71 32 32 60 82 502 507 4.628 -3.152 -3.702 6.713 0.000
> 58 21 (g) -71 31 31 60 82 509 502 1.674 -0.586 -0.383 1.814 0.000
> 59 -2 (ubar) -71 30 30 60 82 0 509 16.393 -2.589 -7.751 18.317 0.000
> 60 223 (omega) -83 47 59 91 93 0 0 10.840 -76.360 -79.500 110.766 0.788
> 61 211 pi+ 83 47 59 0 0 0 0 4.501 -34.269 -35.544 49.578 0.140
> 62 -211 pi- 83 47 59 0 0 0 0 3.908 -26.783 -28.242 39.118 0.140
> 63 211 pi+ 83 47 59 0 0 0 0 0.736 -8.724 -9.584 12.982 0.140
> 64 -211 pi- 83 47 59 0 0 0 0 0.790 -6.143 -6.360 8.878 0.140
> 65 211 pi+ 83 47 59 0 0 0 0 1.017 -6.780 -6.419 9.393 0.140
> 66 -213 (rho-) -83 47 59 83 84 0 0 1.710 -6.284 -6.707 9.378 0.741
> 67 211 pi+ 83 47 59 0 0 0 0 0.428 -5.789 -5.982 8.337 0.140
> 68 -211 pi- 83 47 59 0 0 0 0 0.652 -0.892 -0.756 1.346 0.140
> 69 221 (eta) -83 47 59 94 95 0 0 0.580 -1.762 -1.229 2.292 0.548
> 70 113 (rho0) -83 47 59 85 86 0 0 0.897 -1.583 -2.027 2.915 1.040
> 71 211 pi+ 83 47 59 0 0 0 0 0.048 -0.593 -0.282 0.673 0.140
> 72 223 (omega) -84 47 59 96 98 0 0 1.533 -1.064 -1.059 2.285 0.785
> 73 1114 (Delta-) -84 47 59 87 88 0 0 2.708 -2.521 -0.788 3.991 1.274
> 74 211 pi+ 84 47 59 0 0 0 0 1.459 -1.028 -1.155 2.131 0.140
> 75 -3122 (Lambdabar0) -84 47 59 99 100 0 0 3.115 -2.145 -2.022 4.431 1.116
> 76 -321 K- 84 47 59 0 0 0 0 3.796 -2.570 -2.433 5.214 0.494
> 77 211 pi+ 84 47 59 0 0 0 0 0.797 -0.635 -0.265 1.062 0.140
> 78 221 (eta) -84 47 59 101 102 0 0 1.750 -1.802 -2.322 3.464 0.548
> 79 111 (pi0) -84 47 59 103 104 0 0 5.564 -1.943 -2.834 6.541 0.135
> 80 -211 pi- 84 47 59 0 0 0 0 0.563 -0.054 0.021 0.583 0.140
> 81 213 (rho+) -84 47 59 89 90 0 0 9.509 -2.375 -5.204 11.141 0.985
> 82 -211 pi- 84 47 59 0 0 0 0 5.414 -0.751 -2.258 5.915 0.140
> 83 -211 pi- 91 66 0 0 0 0 0 0.264 -1.794 -1.572 2.405 0.140
> 84 111 (pi0) -91 66 0 105 106 0 0 1.446 -4.489 -5.135 6.973 0.135
> 85 211 pi+ 91 70 0 0 0 0 0 0.667 -1.239 -0.826 1.638 0.140
> 86 -211 pi- 91 70 0 0 0 0 0 0.230 -0.344 -1.201 1.278 0.140
> 87 2112 n0 91 73 0 0 0 0 0 2.350 -2.395 -0.565 3.530 0.940
> 88 -211 pi- 91 73 0 0 0 0 0 0.358 -0.126 -0.222 0.461 0.140
> 89 211 pi+ 91 81 0 0 0 0 0 1.160 -0.119 -0.367 1.230 0.140
> 90 111 (pi0) -91 81 0 107 108 0 0 8.349 -2.256 -4.838 9.911 0.135
> 91 211 pi+ 91 60 0 0 0 0 0 5.216 -36.159 -37.505 52.358 0.140
> 92 -211 pi- 91 60 0 0 0 0 0 5.145 -36.342 -38.020 52.846 0.140
> 93 111 (pi0) -91 60 0 109 110 0 0 0.480 -3.859 -3.975 5.562 0.135
> 94 22 gamma 91 69 0 0 0 0 0 0.496 -1.666 -1.276 2.156 0.000
> 95 22 gamma 91 69 0 0 0 0 0 0.084 -0.096 0.046 0.136 0.000
> 96 111 (pi0) -91 72 0 111 112 0 0 0.093 -0.079 0.082 0.200 0.135
> 97 11 e- 91 72 0 0 0 0 0 1.225 -0.840 -0.969 1.773 0.001
> 98 -11 e+ 91 72 0 0 0 0 0 0.216 -0.146 -0.172 0.312 0.001
> 99 -2112 nbar0 91 75 0 0 0 0 0 2.650 -1.856 -1.841 3.839 0.940
> 100 111 (pi0) -91 75 0 113 114 0 0 0.465 -0.290 -0.181 0.593 0.135
> 101 22 gamma 91 78 0 0 0 0 0 0.238 -0.545 -0.673 0.898 0.000
> 102 22 gamma 91 78 0 0 0 0 0 1.512 -1.256 -1.649 2.565 0.000
> 103 22 gamma 91 79 0 0 0 0 0 1.907 -0.613 -0.919 2.204 0.000
> 104 22 gamma 91 79 0 0 0 0 0 3.657 -1.331 -1.915 4.337 0.000
> 105 22 gamma 91 84 0 0 0 0 0 1.214 -3.729 -4.196 5.744 0.000
> 106 22 gamma 91 84 0 0 0 0 0 0.232 -0.760 -0.939 1.230 0.000
> 107 22 gamma 91 90 0 0 0 0 0 4.970 -1.290 -2.823 5.859 0.000
> 108 22 gamma 91 90 0 0 0 0 0 3.380 -0.966 -2.015 4.051 0.000
> 109 22 gamma 91 93 0 0 0 0 0 0.440 -3.127 -3.246 4.528 0.000
> 110 22 gamma 91 93 0 0 0 0 0 0.040 -0.732 -0.729 1.034 0.000
> 111 22 gamma 91 96 0 0 0 0 0 0.072 0.019 0.008 0.075 0.000
> 112 22 gamma 91 96 0 0 0 0 0 0.020 -0.098 0.074 0.124 0.000
> 113 22 gamma 91 100 0 0 0 0 0 0.412 -0.195 -0.144 0.478 0.000
> 114 22 gamma 91 100 0 0 0 0 0 0.054 -0.095 -0.037 0.115 0.000
> Charge sum: 0.000 Momentum sum: 0.000 0.000 -0.000 600.000 600.000
>
> -------- End PYTHIA Event Listing -----------------------------------------------------------------------------------------------
>
> Pythia::next(): 1000 events have been generated
>
> Pythia::next(): 2000 events have been generated
>
> Pythia::next(): 3000 events have been generated
> /mnt/e/mg5_amc_V3_4_1/eezjets_test2/Events/run_01/run_shower.sh: line 2: 62704 Segmentation fault LD_LIBRARY_PATH=/mnt/e/mg5_amc_V3_4_1/HEPTools/lib:$LD_LIBRARY_PATH /mnt/e/mg5_amc_V3_4_1/HEPTools/MG5aMC_PY8_interface/MG5aMC_PY8_interface tag_1_pythia8.cmd
>
>
>
>
>
>
>
>
>
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

I am studying the double higgs diagram, especially in this case the channel,

e+ e- > Z H H QED = 6, Z > mu+ mu-, H > b b~

The sample in the question posted is intended to be a full backgraound sample of the form

e+ e- > z + n jets where we want n up to 4. with the z decaying to mu+ mu- on-shell.

Ideally I wanted to compare against a full background and so I was using 'all' but depending on
what might cause issues and how much they contribute to the background I can try narrowing
down that definition.

If you have any suggestions to simulate such a background, I would greatly appreciate it. There is a
paper that looks at only certain diagrams that they believe would pose a problem to the signal but
I was hoping for a more complete approach than that if possible.

Thanks,
Nikhil

> e+ e- > z + n jets where we want n up to 4. with the z decaying to mu+
> mu- on-shell.

So this would rather be this syntax:
import model sm-no_b_mass # I guess you want to be in five flavor scheme here.
generate e+ e- > z 2j , z > mu+ mu-
add process e+ e- > z 3j , z > mu+ mu-
add process e+ e- > z 4j , z > mu+ mu-

Given the signal that you have, it might make sense to not include the 2j process since it will likey be vetoed by your analysis.
import model sm-no_b_mass # I guess you want to be in five flavor scheme here.
generate e+ e- > z 3j , z > mu+ mu-
add process e+ e- > z 4j , z > mu+ mu-

Cheers,

Olivier
> On 28 Oct 2022, at 04:50, Nikhil Bachhawat <email address hidden> wrote:
>
> Question #703625 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/703625
>
> Status: Answered => Open
>
> Nikhil Bachhawat is still having a problem:
> Hi,
>
> I am studying the double higgs diagram, especially in this case the
> channel,
>
> e+ e- > Z H H QED = 6, Z > mu+ mu-, H > b b~
>
> The sample in the question posted is intended to be a full backgraound
> sample of the form
>
> e+ e- > z + n jets where we want n up to 4. with the z decaying to mu+
> mu- on-shell.
>
> Ideally I wanted to compare against a full background and so I was using 'all' but depending on
> what might cause issues and how much they contribute to the background I can try narrowing
> down that definition.
>
> If you have any suggestions to simulate such a background, I would greatly appreciate it. There is a
> paper that looks at only certain diagrams that they believe would pose a problem to the signal but
> I was hoping for a more complete approach than that if possible.
>
> Thanks,
> Nikhil
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

import model sm-no_b_mass or even using 4 flavor sm which is what I was doing earlier both give an error message like the one below.

generate e+ e- > z 2j , z > mu+ mu-
add process e+ e- > z 3j , z > mu+ mu-
add process e+ e- > z 4j , z > mu+ mu-

************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.4.1 2022-09-01 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /mnt/e/mg5_amc_V3_4_1/eezjets_test3/Cards/me5_configuration.txt
INFO: load configuration from /mnt/e/mg5_amc_V3_4_1/input/mg5_configuration.txt
INFO: load configuration from /mnt/e/mg5_amc_V3_4_1/eezjets_test3/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
No valid eps viewer found. Please set in ./input/mg5_configuration.txt
No valid web browser found. Please set in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/================ Description =================|=========== values ===========|====== other options ======\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=========================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>1
The following switches determine which programs are run:
/================ Description =================|=========== values ===========|====== other options ======\
| 1. Choose the shower/hadronization program | shower = Pythia8 | OFF |
| 2. Choose the detector simulation program | detector = Not Avail. | Please install module |
| 3. Choose an analysis package (plot/convert) | analysis = Not Avail. | Please install module |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = Not Avail. | Please install module |
\=========================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=OFF' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, pythia8, enter path][90s to answer]
>2
open /mnt/e/mg5_amc_V3_4_1/eezjets_test3/Cards/run_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. pythia8 : pythia8_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, pythia8, enter path]
>0
INFO: Update the dependent parameter of the param_card.dat
Generating 300000 events with run name run_01
survey run_01
INFO: compile directory
Not able to open file /mnt/e/mg5_amc_V3_4_1/eezjets_test3/crossx.html since no program configured.Please set one in ./input/mg5_configuration.txt
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: Compiling for process 1/4.
INFO: P3_ll_zggqq_z_ll
INFO: Compiling for process 2/4.
INFO: P3_ll_zqqqq_z_ll
INFO: Compiling for process 3/4.
INFO: P2_ll_zgqq_z_ll
INFO: Compiling for process 4/4.
INFO: P1_ll_zqq_z_ll
INFO: P3_ll_zggqq_z_ll
INFO: P3_ll_zqqqq_z_ll
INFO: P2_ll_zgqq_z_ll
INFO: P1_ll_zqq_z_ll
INFO: Idle: 1, Running: 17, Completed: 38 [ current time: 11h00 ]
INFO: Idle: 0, Running: 18, Completed: 38 [ 0.027s ]
INFO: Idle: 0, Running: 7, Completed: 49 [ 8.2s ]
INFO: Idle: 0, Running: 3, Completed: 53 [ 16.1s ]
INFO: Idle: 0, Running: 0, Completed: 56 [ 34.7s ]
INFO: Idle: 0, Running: 0, Completed: 56 [ 34.7s ]
INFO: End survey
refine 300000
Creating Jobs
INFO: Refine results to 300000
INFO: Generating 300000.0 unweighted events.
sum of cpu time of last step: 5m48s
INFO: Effective Luminosity 48517340.022755675 pb^-1
INFO: need to improve 24 channels
- Current estimate of cross-section: 0.0074200275578 +- 0.00010011084868982622
    P3_ll_zggqq_z_ll
    P3_ll_zqqqq_z_ll
    P2_ll_zgqq_z_ll
    P1_ll_zqq_z_ll
INFO: Idle: 226, Running: 64, Completed: 35 [ current time: 11h01 ]
INFO: Idle: 225, Running: 64, Completed: 36 [ 0.74s ]
INFO: Idle: 186, Running: 64, Completed: 75 [ 3.9s ]
INFO: Idle: 151, Running: 64, Completed: 110 [ 6.9s ]
INFO: Idle: 111, Running: 64, Completed: 150 [ 9.9s ]
INFO: Idle: 66, Running: 64, Completed: 195 [ 13.3s ]
INFO: Idle: 25, Running: 64, Completed: 236 [ 16.6s ]
INFO: Idle: 0, Running: 55, Completed: 270 [ 19.6s ]
INFO: Idle: 0, Running: 17, Completed: 308 [ 22.6s ]
INFO: Idle: 0, Running: 8, Completed: 317 [ 26.2s ]
INFO: Idle: 0, Running: 7, Completed: 318 [ 44.8s ]
INFO: Idle: 0, Running: 6, Completed: 319 [ 2m 54s ]
INFO: Idle: 0, Running: 5, Completed: 320 [ 3m 0s ]
INFO: Idle: 0, Running: 4, Completed: 321 [ 5m 0s ]
INFO: Idle: 0, Running: 3, Completed: 322 [ 5m 16s ]
INFO: Idle: 0, Running: 2, Completed: 323 [ 5m 23s ]
INFO: Idle: 0, Running: 0, Completed: 325 [ 5m 28s ]
INFO: Idle: 0, Running: 0, Completed: 325 [ 5m 28s ]
INFO: Combining runs
sum of cpu time of last step: 34m23s
INFO: finish refine
refine 300000 --treshold=0.9
No need for second refine due to stability of cross-section
INFO: Combining Events
  === Results Summary for run: run_01 tag: tag_1 ===

     Cross-section : 0.007788 +- 3.828e-06 pb
     Nb of events : 300000

store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
decay_events -from_cards
INFO: Running Pythia8 [arXiv:1410.3012]
No user-defined value for Pythia8 parameter 'JetMatching:nJetMax'. Setting it automatically to 2.
Splitting .lhe event file for PY8 parallelization...
Submitting Pythia8 jobs...
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_7/run_PY8.sh: line 2: 4820 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
WARNING: program /mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_7/run_PY8.sh 7 launch ends with non zero status: 139. Stop all computation
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_11/run_PY8.sh: line 2: 4828 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_15/run_PY8.sh: line 2: 4836 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_6/run_PY8.sh: line 2: 4818 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_18/run_PY8.sh: line 2: 4842 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Terminated
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_16/run_PY8.sh: line 2: 4838 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_21/run_PY8.sh: line 2: 4848 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_9/run_PY8.sh: line 2: 4823 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Terminated
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_17/run_PY8.sh: line 2: 4840 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_8/run_PY8.sh: line 2: 4822 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Terminated
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_20/run_PY8.sh: line 2: 4846 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_12/run_PY8.sh: line 2: 4830 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Terminated
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_10/run_PY8.sh: line 2: 4826 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_22/run_PY8.sh: line 2: 4850 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Terminated
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_26/run_PY8.sh: line 2: 4858 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_5/run_PY8.sh: line 2: 4816 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_13/run_PY8.sh: line 2: 4832 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_24/run_PY8.sh: line 2: 4854 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_28/run_PY8.sh: line 2: 4862 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_32/run_PY8.sh: line 2: 4870 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_25/run_PY8.sh: line 2: 4856 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_29/run_PY8.sh: line 2: 4864 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_27/run_PY8.sh: line 2: 4860 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_31/run_PY8.sh: line 2: 4868 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_34/run_PY8.sh: line 2: 4874 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_35/run_PY8.sh: line 2: 4876 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_23/run_PY8.sh: line 2: 4852 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 1 Idle, 32 Running, 31 Done [1m42s]
Pythia8 shower jobs: 0 Idle, 33 Running, 31 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_36/run_PY8.sh: line 2: 4878 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 32 Running, 32 Done [1m42s]
Terminated
Pythia8 shower jobs: 0 Idle, 31 Running, 33 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_33/run_PY8.sh: line 2: 4872 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 30 Running, 34 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_37/run_PY8.sh: line 2: 4880 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 29 Running, 35 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_19/run_PY8.sh: line 2: 4844 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 28 Running, 36 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_38/run_PY8.sh: line 2: 4882 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 27 Running, 37 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_40/run_PY8.sh: line 2: 4886 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 26 Running, 38 Done [1m42s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_39/run_PY8.sh: line 2: 4884 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 25 Running, 39 Done [1m43s]
Terminated
Pythia8 shower jobs: 0 Idle, 24 Running, 40 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_45/run_PY8.sh: line 2: 4896 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 23 Running, 41 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_43/run_PY8.sh: line 2: 4892 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 22 Running, 42 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_42/run_PY8.sh: line 2: 4890 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 21 Running, 43 Done [1m43s]
Terminated
Pythia8 shower jobs: 0 Idle, 20 Running, 44 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_47/run_PY8.sh: line 2: 4903 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 19 Running, 45 Done [1m43s]
Terminated
Pythia8 shower jobs: 0 Idle, 18 Running, 46 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_46/run_PY8.sh: line 2: 4898 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 17 Running, 47 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_50/run_PY8.sh: line 2: 4914 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 16 Running, 48 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_51/run_PY8.sh: line 2: 4920 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 15 Running, 49 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_48/run_PY8.sh: line 2: 4906 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 14 Running, 50 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_49/run_PY8.sh: line 2: 4911 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 13 Running, 51 Done [1m43s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_53/run_PY8.sh: line 2: 4927 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 12 Running, 52 Done [1m44s]
Terminated
Pythia8 shower jobs: 0 Idle, 11 Running, 53 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_54/run_PY8.sh: line 2: 4932 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 10 Running, 54 Done [1m44s]
Terminated
Pythia8 shower jobs: 0 Idle, 9 Running, 55 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_59/run_PY8.sh: line 2: 4953 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 8 Running, 56 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_57/run_PY8.sh: line 2: 4945 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 7 Running, 57 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_61/run_PY8.sh: line 2: 4963 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 6 Running, 58 Done [1m44s]
Terminated
Pythia8 shower jobs: 0 Idle, 5 Running, 59 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_60/run_PY8.sh: line 2: 4958 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 4 Running, 60 Done [1m44s]
/mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_63/run_PY8.sh: line 2: 4972 Segmentation fault ./MG5aMC_PY8_interface PY8Card.dat &> PY8_log.txt
Pythia8 shower jobs: 0 Idle, 3 Running, 61 Done [1m44s]
Terminated
Pythia8 shower jobs: 0 Idle, 2 Running, 62 Done [1m44s]
Terminated
Pythia8 shower jobs: 0 Idle, 1 Running, 63 Done [1m44s]
Pythia8 shower jobs: 0 Idle, 0 Running, 64 Done [1m44s]
Command "generate_events run_01" interrupted with error:
Exception : program /mnt/e/mg5_amc_V3_4_1/eezjets_test3/Events/run_01/PY8_parallelization/split_7/run_PY8.sh 7 launch ends with non zero status: 139. Stop all computation
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/mnt/e/mg5_amc_V3_4_1/eezjets_test3/run_01_tag_1_debug.log'.
Please attach this file to your report.
INFO: storing files of previous run
INFO: Done
INFO:
quit

Even in a single core run I get a similar error and segmentation fault and even with requesting 300k events and running pythia8 manually, i get very few events remaining. I am trying only the 3j and 4j process in single core and will update if I get an error there too.

Thanks,
Nikhil

Hi,

Segfaut in pythia8 typically occurs when pythia8 fails to correctly identify the basic production process.
But since I'm not an author of Pythia8 I can not comment more on this and on how fix the pythia8 setup to fix such issue.

Cheers,

Olivier

> On 30 Oct 2022, at 17:35, Nikhil Bachhawat <email address hidden> wrote:
>
> : 300