Running MG 2.6.5 (also 2.6.6) to generate events, shower and simulate detector, no pythia or delphes event is generated. The only output is the parton-level lhe events.

Asked by HH on 2019-09-01

Running MG 2.6.5 (also 2.6.6) to generate events, shower and simulate detector, no pythia or delphes event is generated. The only output is the parton-level lhe events. why?

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2019-09-02
Last query:
2019-09-02
Last reply:
2019-09-02

Do you look at the Pythia log and/or the Delphi’s log file. I assume that you have successfully generate your Parton event.

Get Outlook for iOS<https://aka.ms/o0ukef>
________________________________
From: <email address hidden> <email address hidden> on behalf of jack <email address hidden>
Sent: Sunday, September 1, 2019 4:57:52 PM
To: Olivier Mattelaer <email address hidden>
Subject: [Question #683567]: running ./bin/generate_event with --laststep=delphes option, no delphes event is stored in output root file

New question #683567 on MadGraph5_aMC@NLO:
https://answers.launchpad.net/mg5amcnlo/+question/683567

Running ./bin/generate_event and also running event generator in CLI, I get neither pythia event root file nor delphes event root file.
I also used ./bin/madEvent generator with option --laststep=delphes with no result. Why is this? I use python 2.6.6 and MadGraph 2.6.6. Is my madgraph works properly with this version of python ?

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

HH (hh-madgrapher) said : #3

Yes the parton level events are fully generated, however pythia and delphes events are NOT.
Looking at the pythia log file I found:

*------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 20 Abort from Pythia::next: processLevel failed; giving up |
 | 1 Abort from Pythia::next: reached end of Les Houches Events File |
 | 100 Error in ResonanceDecays::next: failed to find workable decay channel |
 | 2 Error in ResonanceWidths::init: resonance mass too small |
 | 2 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 2 Info from SLHAinterface::initSLHA: importing DECAY tables |
 | 2 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 2 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 2 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring QNUMBERS |
 | 2 Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

does this give you any insight to the problem ?

Thanks in advance

I guess that it means that madgraph refuses to activate Pythia or Delphi’s. This typically happens when the installation of those codes fails. Can you check Pythia alone first to check if it works in that case. If not then can you re install Pythia?

Then when Pythia is working do the same for Delphes

Get Outlook for iOS<https://aka.ms/o0ukef>
________________________________
From: <email address hidden> <email address hidden> on behalf of jack <email address hidden>
Sent: Sunday, September 1, 2019 5:27:31 PM
To: Olivier Mattelaer <email address hidden>
Subject: Re: [Question #683567]: running ./bin/generate_event with --laststep=delphes option, no delphes event is stored in output root file

Question #683567 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/683567

jack posted a new comment:
It does not even give me a pythia or delphes log. There are only 2 files
in the run_01 directory which are unweighted-event files containing
parton level events. Parton level events are generated properly but the
problem is shower and detector !!!

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

Looks like Pythia can not handle a decay. Likely an error on your side. I do not have enough information to comment more

Get Outlook for iOS<https://aka.ms/o0ukef>
________________________________
From: <email address hidden> <email address hidden> on behalf of jack <email address hidden>
Sent: Sunday, September 1, 2019 5:46:51 PM
To: Olivier Mattelaer <email address hidden>
Subject: Re: [Question #683567]: running ./bin/generate_event with --laststep=delphes option, no delphes event is stored in output root file

Question #683567 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/683567

jack posted a new comment:
Yes the parton level events are fully generated, however pythia and delphes events are NOT.
Looking at the pythia log file I found:

*------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 20 Abort from Pythia::next: processLevel failed; giving up |
 | 1 Abort from Pythia::next: reached end of Les Houches Events File |
 | 100 Error in ResonanceDecays::next: failed to find workable decay channel |
 | 2 Error in ResonanceWidths::init: resonance mass too small |
 | 2 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 2 Info from SLHAinterface::initSLHA: importing DECAY tables |
 | 2 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 2 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 2 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring QNUMBERS |
 | 2 Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

does this give you any insight to the problem ?

Thanks in advance

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

HH (hh-madgrapher) said : #6

what information do you need to guess the problem?

Thanks for your consideration and time

If you have a debug file, it typically contains the information that I need.
Otherwise the banner associated to the run.

But I guess that:
https://answers.launchpad.net/mg5amcnlo/+question/683569
is your issue no?

Olivier

> On 2 Sep 2019, at 08:22, jack <email address hidden> wrote:
>
> Question #683567 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/683567
>
> Status: Answered => Open
>
> jack is still having a problem:
> what information do you need to guess the problem?
>
> Thanks for your consideration and time
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

HH (hh-madgrapher) said : #8

Dear Dr. Olivier,

I want MG to generate parton-level events, then shower using Pythia8 and then simulate detector using Delphes. The problem is that no pythia events and no delphes events are generated. MG just generates parton-level events.

I use the command
./bin/mg5_aMC 4proc.txt
where 4proc.txt contains

import model ALP_linear_UFO/
define jb = g u c d s u~ c~ d~ s~ b b~
generate p p > jb jb ax
output test2

Then I use
cd test2
./bin/madevent ../4gen.txt
where ../4gen.txt contains

multi_run 1
shower = Pythia8
detector = Delphes
analysis = OFF
madspin = OFF
reweight = OFF
0
set cgtil 0.5
set cwtil 0
set cbtil 0
set caphi 0
set width 9000005 Auto
set nevents 100
set ebeam1 6500
set ebeam2 6500
0

Here is the banner content:

<!--
#*********************************************************************
# *
# MadGraph5_aMC@NLO *
# *
# Going Beyond *
# *
# http://madgraph.hep.uiuc.edu *
# http://madgraph.phys.ucl.ac.be *
# http://amcatnlo.cern.ch *
# *
# The MadGraph5_aMC@NLO team *
# *
#....................................................................*
# *
# This file contains all the information necessary to reproduce *
# the events generated: *
# *
# 1. software version *
# 2. proc_card : code generation info including model *
# 3. param_card : model primary parameters in the LH format *
# 4. run_card : running parameters (collider and cuts) *
# 5. pythia_card : present only if pythia has been run *
# 6. pgs_card : present only if pgs has been run *
# 7. delphes_cards : present only if delphes has been run *
# *
# *
#*********************************************************************
-->
<MGVersion>
#2.6.5
</MGVersion>
<MG5ProcCard>
<![CDATA[
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model /home/haghighat/sw/MG5_aMC_v2_6_5/models/ALP_linear_UFO
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define jb = g u c d s u~ c~ d~ s~ b b~
generate p p > jb jb ax
output test2
]]>
</MG5ProcCard>
<MGProcCard>
#*********************************************************************
# MadGraph/MadEvent *
# http://madgraph.hep.uiuc.edu *
# *
# proc_card.dat *
#*********************************************************************
# *
# This Files is generated by MADGRAPH 5 *
# *
# WARNING: This Files is generated for MADEVENT (compatibility issue)*
# This files is NOT a valid MG4 proc_card.dat *
# Running this in MG4 will NEVER reproduce the result of MG5*
# *
#*********************************************************************
#*********************************************************************
# Process(es) requested : mg2 input *
#*********************************************************************
# Begin PROCESS # This is TAG. Do not modify this line
p p > jb jb ax #Process
# Be carefull the coupling are here in MG5 convention

end_coup # End the couplings input

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

# End MULTIPARTICLES # This is TAG. Do not modify this line
</MGProcCard>
<MGRunCard>
<![CDATA[
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# To display more options, you can type the command: *
# update full_run_card *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  tag_1 = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  100 = nevents ! Number of unweighted events requested
 21 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
  1 = lpp1 ! beam 1 type
  1 = lpp2 ! beam 2 type
  6500.0 = ebeam1 ! beam 1 total energy in GeV
  6500.0 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
  nn23lo1 = pdlabel ! PDF set
  230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# To see heavy ion options: type "update ion_pdf"
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
  False = fixed_ren_scale ! if .true. use fixed ren scale
  False = fixed_fac_scale ! if .true. use fixed fac scale
  91.188 = scale ! fixed ren scale
  91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
  91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
  -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
  1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  3.0 = lhe_version ! Change the way clustering information pass to shower.
  True = clusinfo ! include clustering tag in output
  average = event_norm ! average/sum. Normalization of the weight in the LHEF

#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
  0 = ickkw ! 0 no matching, 1 MLM
  1.0 = alpsfact ! scale factor for QCD emission vx
  False = chcluster ! cluster only according to channel diag
  4 = asrwgtflavor ! highest quark flavor for a_s reweight
  False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
                                   ! (turn off for VBF and single top processes)
  0.0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
  0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
  None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
  {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
  False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
  20.0 = ptj ! minimum pt for the jets
  0.0 = ptb ! minimum pt for the b
  10.0 = 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.4 = 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.4 = 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 *
#*********************************************************************
  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
#***********************************************************************
# 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
#*********************************************************************
# 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 *
#*********************************************************************
  True = use_syst ! Enable systematics studies
#
  systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments ! see: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule
# Syscalc is deprecated but to see the associate options type'update syscalc'#*********************************************************************
# Additional hidden parameters
#*********************************************************************
  ['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_arguments # Choose the argment to pass to the systematics command. like --mur=0.25,1,4. Look at the help of the systematics function for more details.
]]>
</MGRunCard>
<slha>
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR ALPPARS
###################################
BLOCK ALPPARS #
      1 1.000000e+03 # fa
      2 5.000000e-01 # cgtil
      3 0.000000e+00 # cwtil
      4 0.000000e+00 # cbtil
      5 0.000000e+00 # caphi
###################################
## INFORMATION FOR CKMBLOCK
###################################
BLOCK CKMBLOCK #
      1 2.277360e-01 # cabi
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      1 5.040000e-03 # md
      2 2.550000e-03 # mu
      3 1.010000e-01 # ms
      4 1.270000e+00 # mc
      5 4.700000e+00 # mb
      6 1.720000e+02 # mt
      11 5.110000e-04 # me
      13 1.056600e-01 # mmu
      15 1.777000e+00 # mta
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh
      9000005 1.000000e-03 # ma
      12 0.000000e+00 # ve : 0.0
      14 0.000000e+00 # vm : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 7.982436e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.279000e+02 # aewm1
      2 1.166370e-05 # gf
      3 1.184000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      1 5.040000e-03 # ymdo
      2 2.550000e-03 # ymup
      3 1.010000e-01 # yms
      4 1.270000e+00 # ymc
      5 4.700000e+00 # ymb
      6 1.720000e+02 # ymt
      11 5.110000e-04 # yme
      13 1.056600e-01 # ymm
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 9000005 # ax
      1 0 # 3 times electric charge
      2 1 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.508336e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495200e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 4.070000e-03
#
# PDG Width
DECAY 9000005 1.591549e-16
# BR NDA ID1 ID2 ...
   1.000000e+00 2 21 21 # 1.59154943092e-16
</slha>
<MGGenerationInfo>
# Number of Events : 100
# Integrated weight (pb) : 3.89547E+04
</MGGenerationInfo>

Here is the tail of tag_1_pythia8.log:

*------- PYTHIA Event and Cross Section Statistics -------------------------------------------------------------*
 | |
 | Subprocess Code | Number of events | sigma +- delta |
 | | Tried Selected Accepted | (estimated) (mb) |
 | | | |
 |-----------------------------------------------------------------------------------------------------------------|
 | | | |
 | Les Houches User Process(es) 9999 | 100 100 0 | 0.000e+00 0.000e+00 |
 | ... whereof user classification code 1 | 100 100 0 | |
 | | | |
 | sum | 100 100 0 | 0.000e+00 0.000e+00 |
 | |
 *------- End PYTHIA Event and Cross Section Statistics ----------------------------------------------------------*

 *------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------*
 | |
 | times message |
 | |
 | 20 Abort from Pythia::next: processLevel failed; giving up |
 | 1 Abort from Pythia::next: reached end of Les Houches Events File |
 | 100 Error in ResonanceDecays::next: failed to find workable decay channel |
 | 2 Error in ResonanceWidths::init: resonance mass too small |
 | 2 Info from SLHAinterface::initSLHA: No MODSEL found, keeping internal SUSY switched off |
 | 2 Info from SLHAinterface::initSLHA: importing DECAY tables |
 | 2 Info from SLHAinterface::initSLHA: importing MASS entries |
 | 2 Warning in SLHAinterface::initSLHA: ignoring DECAY tables |
 | 2 Warning in SLHAinterface::initSLHA: ignoring MASS entries |
 | 1 Warning in SLHAinterface::initSLHA: ignoring QNUMBERS |
 | 2 Warning in SLHAinterface::initSLHA: ignoring empty DECAY tables |
 | |
 *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------*

 Contribution of sample 0 to the inclusive cross section : 0.00000000e+00 +- 0.00000000e+00

Inclusive cross section: 0.00000000e+00 +- 0.00000000e+00 mb

Any help would be appreciated

Hi,

You might want to contact pythia author for dedicated help on this here.
My guess is that the width is so small that you have issue with the hadronization model.
(indeed your particles decay into two gluon according to your param_card.)
But the width is so small that this is not physically possible. Your particle should decay directly into meson/hadron. But this is something that I do not know how to compute since this is not perturbative.

Note that for Delphes, you also have to worried about displaced vertex here.

Cheers,

Olivier

HH (hh-madgrapher) said : #10

Dear Dr. Olivier,

Thanks for your answer. I ran MG using the Standard Model instead of the ALP and it worked well. I got a delphes root file containing simulated events as I needed. However, following the same procedure for the ALP model, I gain no result. Is it possible that the ALP model file is the source of problem ?

My ALP model files contains these files:

ALP_linear_UFO.log decays.pyo model.pkl propagators.py
coupling_orders.py function_library.py object_library.py propagators.pyo
coupling_orders.pyo function_library.pyo object_library.pyo vertices.py
couplings.py __init__.py parameters.py vertices.pyo
couplings.pyo __init__.pyo parameters.pyo write_param_card.py
CT_couplings.py lorentz.py particles.py write_param_card.pyo
decays.py lorentz.pyo particles.pyo

It was a tarball. I just used tar -xvf command to untar it.

HH (hh-madgrapher) said : #11

For completeness, I should say that I want the quarks and gluons to produce jets, so that I can study the di-jet process p p > jb jb ax. The ax doesn't decay in the detector.

HH (hh-madgrapher) said : #12

I just produced the p p > g g events in ALP model with no problem. delphes output root file contains events as I wanted. the difference between this process and the previous process is the ax (axion-like) particle. So the problem is now reduced to the ax particle. But how come ?!!

I see the below warning when I include ax in the process:

WARNING: Particle 9000005 with very small width (1.59155e-16): Learn about special handling here: https://answers.launchpad.net/mg5amcnlo/+faq/3053

I read the document. states that it uses a fake width and I can change the threshold using command
 1e-12 = small_width_treatment

does changing the threshold make anything better?

HH (hh-madgrapher) said : #13

I reduced the threshold of small_width_treatment to 1e-20 to prevent using fake width and I got the same result. no delphes or pythia event in output.

One solution might be to forbid the decay of ax within pythia8 (please read py8 manual to learn how to do that).
Note that Delphes can then crash on that later. And you likely need to change the configuration of Delphes to consider it as missing energy.
Another idea is to change the PDG code to put it as a DM candidate such that Delphes (and maybe even PY8) will handle that particle correctly. (Now this idea might not work).

Cheers,

Olivier

HH (hh-madgrapher) said : #15

Thanks Olivier Mattelaer, that solved my question.

HH (hh-madgrapher) said : #16

Thanks Dr Olivier for your advice.

In previous runs, the ax decay width was set to Auto. Now, I just left it to stay zero and ran the program. delphes events are fully generated with no problem.

This is actually what I wanted since the ax particle is assumed to be invisible (no decay in detector) in this study.

again thank you