pythia-pgs and extension of the MSSM

Hi Chiara,

I will forward this to the FR author, in order to see if they are aware of such problem. Otherwise, you might need to either report to pythia author (or to use pythia8).
The pythia author may suggest to use pythia8 since they want to drop pythia6 anyway.

Cheers,

Olivier

On Jul 11, 2013, at 3:36 PM, Chiara Arina <email address hidden> wrote:

> New question #232276 on MadGraph5:
> https://answers.launchpad.net/madgraph5/+question/232276
>
> Dear MadGraph team,
>
> I have implemented an extension of the MSSM with feynrules and analyzed some processes specific to the model with madgraph. Since at parton level everything is working fine, I would like to go one step further and simulate the events with pythia-psg. However the pythia-pgs part does not work properly because the model is not implement in pythia. I tried the same process within the MSSM and everything is fine. The same process within the nmssm implemented in madgraph produces the same error as the model I implemented (see below). Is there a way to make pythia-pgs or pythia working with extensions of the mssm?
>
> Thanks,
>
> best
>
> Chiara
>
>
> Here you can find the the log from the pythia run in the case of the nmssm (the process is direct stau production):
>
> 1
> ******************************************************************************
> ******************************************************************************
> ** **
> ** **
> ** *......* Welcome to the Lund Monte Carlo! **
> ** *:::!!:::::::::::* **
> ** *::::::!!::::::::::::::* PPP Y Y TTTTT H H III A **
> ** *::::::::!!::::::::::::::::* P P Y Y T H H I A A **
> ** *:::::::::!!:::::::::::::::::* PPP Y T HHHHH I AAAAA **
> ** *:::::::::!!:::::::::::::::::* P Y T H H I A A **
> ** *::::::::!!::::::::::::::::*! P Y T H H III A A **
> ** *::::::!!::::::::::::::* !! **
> ** !! *:::!!:::::::::::* !! This is PYTHIA version 6.426 **
> ** !! !* -><- * !! Last date of change: 16 Nov 2011 **
> ** !! !! !! **
> ** !! !! !! Now is 0 Jan 2000 at 0:00:00 **
> ** !! !! **
> ** !! lh !! Disclaimer: this program comes **
> ** !! !! without any guarantees. Beware **
> ** !! hh !! of errors and use common sense **
> ** !! ll !! when interpreting results. **
> ** !! !! **
> ** !! Copyright T. Sjostrand (2011) **
> ** **
> ** An archive of program versions and documentation is found on the web: **
> ** http://www.thep.lu.se/~torbjorn/Pythia.html **
> ** **
> ** When you cite this program, the official reference is to the 6.4 manual: **
> ** T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 **
> ** (LU TP 06-13, FERMILAB-PUB-06-052-CD-T) [hep-ph/0603175]. **
> ** **
> ** Also remember that the program, to a large extent, represents original **
> ** physics research. Other publications of special relevance to your **
> ** studies may therefore deserve separate mention. **
> ** **
> ** Main author: Torbjorn Sjostrand; Department of Theoretical Physics, **
> ** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; **
> ** phone: + 46 - 46 - 222 48 16; e-mail: <email address hidden> **
> ** Author: Stephen Mrenna; Computing Division, GDS Group, **
> ** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; **
> ** phone: + 1 - 630 - 840 - 2556; e-mail: <email address hidden> **
> ** Author: Peter Skands; CERN/PH-TH, CH-1211 Geneva, Switzerland **
> ** phone: + 41 - 22 - 767 24 47; e-mail: <email address hidden> **
> ** **
> ** **
> ******************************************************************************
> ******************************************************************************
> MSTU(11) changed from 6 to 6
> MDCY(C15,1) changed from 1 to 0
> MSEL changed from 1 to 0
>
> Reading pythia input card
> mstp(61) changed from 2 to 1
> mstp(71) changed from 1 to 1
> mstj(1) changed from 1 to 1
> mstp(81) changed from 1 to 20
> mstu(21) changed from 2 to 1
> Set LHAPATH to /Users/chiara/Documents/software_linux/MadGraph5_v1_5_9/pythia-pgs/src/PDFsets
> Opened file unweighted_events.lhe
> Running as inclusive sample (highest mult.)
>
> Initiating pythia with external process
> 1****************** PYINIT: initialization of PYTHIA routines *****************
> Reading parameter sminputs 1 127.92000000000000
> Reading parameter sminputs 3 0.11720000000000000
> Reading model: nmssm
> IMSS(21) changed from 0 to 24
> IMSS(22) changed from 0 to 24
> Initializing PYR with random seed 21
> MSTP(52) changed from 1 to 2
> Found parameter ickkw 0
> MSTP(98) changed from 0 to 1
> ==== PYTHIA WILL USE LHAPDF ====
> *************************************
> * LHAPDF Version 4.0 *
> *************************************
>
>>>>>>> PDF description: <<<<<<
> CTEQ6L1 - LO with LO alpha_s
> Reference:
> J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky,
> W.K. Tung
> hep-ph/0201195
>>>>>>> <<<<<<
>
> Parametrization: CTEQ6
>
> warning empty CtLhQCDSET call: NAME=IHDN VALUE= 1.0000000000000000
> warning empty CtLhQCDSET call: NAME=QMAX VALUE= 100000.00000000000
> warning empty CtLhQCDSET call: NAME=XMIN VALUE= 9.9999999999999995E-007
> warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
> warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
> ==============================================
> PDFset name /Users/chiara/Documents/software_linux/MadGraph5_v1_5_9/pythia-pgs/src/PDFsets/c
> with 1 members
> ==== initialized. ===========================
> warning empty CtLhQCDSET call: NAME=NX VALUE= 95.000000000000000
> warning empty CtLhQCDSET call: NAME=NT VALUE= 18.000000000000000
> Strong coupling at Mz for PDF is: 0.12978
> ****************** PYSLHA v1.14: SUSY/BSM SPECTRUM INTERFACE *****************
> * (PYSLHA:) Last Change 10 Jun 2010 - P.Z. Skands
> * (PYSLHA:) Reading QNUMBERS for KF = 1000022
> * (PYSLHA:) KF = 1000022 already exists: ~chi_10 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000023
> * (PYSLHA:) KF = 1000023 already exists: ~chi_20 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000025
> * (PYSLHA:) KF = 1000025 already exists: ~chi_30 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000035
> * (PYSLHA:) KF = 1000035 already exists: ~chi_40 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000045
> * (PYSLHA:) Reading QNUMBERS for KF = 1000024
> * (PYSLHA:) KF = 1000024 already exists: ~chi_1+ . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000037
> * (PYSLHA:) KF = 1000037 already exists: ~chi_2+ . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000021
> * (PYSLHA:) KF = 1000021 already exists: ~g . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 35
> * (PYSLHA:) KF = 35 already exists: H0 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 45
> * (PYSLHA:) KF = 45 already exists: . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 36
> * (PYSLHA:) KF = 36 already exists: A0 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 46
> * (PYSLHA:) KF = 46 already exists: . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 37
> * (PYSLHA:) KF = 37 already exists: H+ . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000012
> * (PYSLHA:) KF = 1000012 already exists: ~nu_eL . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000014
> * (PYSLHA:) KF = 1000014 already exists: ~nu_muL . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000016
> * (PYSLHA:) KF = 1000016 already exists: ~nu_tauL . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000011
> * (PYSLHA:) KF = 1000011 already exists: ~e_L- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000013
> * (PYSLHA:) KF = 1000013 already exists: ~mu_L- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000015
> * (PYSLHA:) KF = 1000015 already exists: ~tau_1- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000011
> * (PYSLHA:) KF = 2000011 already exists: ~e_R- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000013
> * (PYSLHA:) KF = 2000013 already exists: ~mu_R- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000015
> * (PYSLHA:) KF = 2000015 already exists: ~tau_2- . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000002
> * (PYSLHA:) KF = 1000002 already exists: ~u_L . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000004
> * (PYSLHA:) KF = 1000004 already exists: ~c_L . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000006
> * (PYSLHA:) KF = 1000006 already exists: ~t_1 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000002
> * (PYSLHA:) KF = 2000002 already exists: ~u_R . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000004
> * (PYSLHA:) KF = 2000004 already exists: ~c_R . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000006
> * (PYSLHA:) KF = 2000006 already exists: ~t_2 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000001
> * (PYSLHA:) KF = 1000001 already exists: ~d_L . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000003
> * (PYSLHA:) KF = 1000003 already exists: ~s_L . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 1000005
> * (PYSLHA:) KF = 1000005 already exists: ~b_1 . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000001
> * (PYSLHA:) KF = 2000001 already exists: ~d_R . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000003
> * (PYSLHA:) KF = 2000003 already exists: ~s_R . Entry ignored.
> * (PYSLHA:) Reading QNUMBERS for KF = 2000005
> * (PYSLHA:) KF = 2000005 already exists: ~b_2 . Entry ignored.
> * (PYSLHA:) Ignoring MASS entry for KF = 5 (SLHA read-in not allowed)
> * (PYSLHA:) Reading MASS entry for KF = 6, pole mass = 171.400
> * (PYSLHA:) Ignoring MASS entry for KF = 15 (SLHA read-in not allowed)
> * (PYSLHA:) Reading MASS entry for KF = 23, pole mass = 91.187
> * (PYSLHA:) Reading MASS entry for KF = 24, pole mass = 80.939
> * (PYSLHA:) Reading MASS entry for KF = 25, pole mass = 119.164
> * (PYSLHA:) Reading MASS entry for KF = 35, pole mass = 1016.551
> * (PYSLHA:) Reading MASS entry for KF = 36, pole mass = 40.398
> * (PYSLHA:) Reading MASS entry for KF = 37, pole mass = 1022.472
> * (PYSLHA:) Reading MASS entry for KF = 45, pole mass = 2112.576
> * (PYSLHA:) Reading MASS entry for KF = 46, pole mass = 1021.042
> * (PYSLHA:) Reading MASS entry for KF = 1000001, pole mass = 1066.519
> * (PYSLHA:) Reading MASS entry for KF = 1000002, pole mass = 1063.635
> * (PYSLHA:) Reading MASS entry for KF = 1000005, pole mass = 866.365
> * (PYSLHA:) Reading MASS entry for KF = 1000006, pole mass = 499.229
> * (PYSLHA:) Reading MASS entry for KF = 1000011, pole mass = 380.176
> * (PYSLHA:) Reading MASS entry for KF = 1000012, pole mass = 372.121
> * (PYSLHA:) Reading MASS entry for KF = 1000015, pole mass = 214.740
> * (PYSLHA:) Reading MASS entry for KF = 1000016, pole mass = 360.341
> * (PYSLHA:) Reading MASS entry for KF = 1000021, pole mass = 1151.543
> * (PYSLHA:) Reading MASS entry for KF = 1000022, pole mass = 208.142
> * (PYSLHA:) Reading MASS entry for KF = 1000023, pole mass = 397.851
> * (PYSLHA:) Reading MASS entry for KF = 1000024, pole mass = 397.829
> * (PYSLHA:) Reading MASS entry for KF = 1000025, pole mass = -963.981
> * (PYSLHA:) Reading MASS entry for KF = 1000035, pole mass = 969.594
> * (PYSLHA:) Reading MASS entry for KF = 1000037, pole mass = 970.137
> * (PYSLHA:) Reading MASS entry for KF = 1000045, pole mass = 2094.274
> * (PYSLHA:) Reading MASS entry for KF = 2000001, pole mass = 1023.754
> * (PYSLHA:) Reading MASS entry for KF = 2000002, pole mass = 1027.259
> * (PYSLHA:) Reading MASS entry for KF = 2000005, pole mass = 992.138
> * (PYSLHA:) Reading MASS entry for KF = 2000006, pole mass = 935.527
> * (PYSLHA:) Reading MASS entry for KF = 2000011, pole mass = 261.024
> * (PYSLHA:) Reading MASS entry for KF = 2000015, pole mass = 374.857
> * (PYSLHA:) Ignoring MASS entry for KF = 1 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 2 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 3 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 4 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 11 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 12 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 13 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 14 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 16 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 21 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 22 (SLHA read-in not allowed)
> * (PYSLHA:) Reading MASS entry for KF = 1000014, pole mass = 372.121
> * (PYSLHA:) Reading MASS entry for KF = 1000013, pole mass = 380.176
> * (PYSLHA:) Reading MASS entry for KF = 2000013, pole mass = 261.024
> * (PYSLHA:) Reading MASS entry for KF = 1000004, pole mass = 1063.635
> * (PYSLHA:) Reading MASS entry for KF = 2000004, pole mass = 1027.259
> * (PYSLHA:) Reading MASS entry for KF = 1000003, pole mass = 1066.519
> * (PYSLHA:) Reading MASS entry for KF = 2000003, pole mass = 1023.754
> * (PYSLHA:) Allowing DECAY table for KF = 6 but this is NOT recommended.
> * (PYSLHA:) Allowing DECAY table for KF = 23 but this is NOT recommended.
> * (PYSLHA:) Allowing DECAY table for KF = 24 but this is NOT recommended.
> * (PYSLHA:) Ignoring DECAY table for KF = 1 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 2 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 3 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 4 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 5 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 11 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 12 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 13 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 14 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 15 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 16 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 21 (SLHA read-in not allowed)
>
>
> ********************************************************
> STDHEP version 5.04.01 - Aug. 29, 2005
> ********************************************************
>
>
> STDXWOPEN WARNING: I/O is initialized for stdhep only
> * (PYSLHA:) Ignoring DECAY table for KF = 22 (SLHA read-in not allowed)
> * (PYSLHA:) Reading SLHA stable particle KF = 1000022, ~chi_10
>
> ==============================================================================
> I I
> I PYTHIA will be initialized for p+ on p+ user configuration I
> I with 7000.000 GeV on 7000.000 GeV beam energies I
> I I
> I corresponding to 14000.000 GeV center-of-mass energy I
> I I
> ==============================================================================
>
> ******** PYMAXI: summary of differential cross-section maximum search ********
>
> ==========================================================
> I I I
> I ISUB Subprocess name I Maximum value I
> I I I
> ==========================================================
> I I I
> I 4 User process 1 I 5.5237D-12 I
> I I I
> ==========================================================
>
> ********************** PYINIT: initialization completed **********************
> Initialisation done
>
> Generating all events
> stdhep initialized
>
>
>
> ***************************************************************************
> * *****TAUOLA LIBRARY: VERSION 2.5 ****** *
> * ***********JUNE 1994*************** *
> * **AUTHORS: S.JADACH, Z.WAS************* *
> * **R. DECKER, M. JEZABEK, J.H.KUEHN***** *
> * **AVAILABLE FROM: WASM AT CERNVM ****** *
> * ***** PUBLISHED IN COMP. PHYS. COMM.*** *
> * *******CERN-TH-5856 SEPTEMBER 1990***** *
> * *******CERN-TH-6195 SEPTEMBER 1991***** *
> * *******CERN-TH-6793 NOVEMBER 1992***** *
> * **5 or more pi dec.: precision limited *
> * ******DEXAY ROUTINE: INITIALIZATION**** *
> * 0 JAK1 = DECAY MODE FERMION1 (TAU+) *
> * 0 JAK2 = DECAY MODE FERMION2 (TAU-) *
> ***************************************************************************
>
>
> ***************************************************************************
>
> _______________________________________________________________________________
> TAUOLA VERSION 2.1 PARAMETER SET 0
> _______________________________________________________________________________
> Printing end run statistics.
> GENP bank used for input.
> Mult. taus in event have same chirality.
> Taus polarized in W/Z rest frame.
> Helicity: 1.00000000
> Lifetime: 295.600006 fs
> _______________________________________________________________________________
> STANDARD MODEL PARAMETERS
> Gfermi = 1.16639203E-05
> Gv = 1.00000000
> Ga = -1.00000000
> Cos(Cabibbo_angle) = 0.975000024
> Sin(Cabibbo_angle) = 0.222204715
> Gamma_e = 4.04914356E-13
> 1/alpha = 137.03598956100001
> 1/(alpha*pi) = 2.3228196271906998E-003
> xk0 = 0.10000000149011612
> _______________________________________________________________________________
> KAON/A1 DECAY PARAMETERS
> bra1 = 0.500000000
> brk0 = 0.500000000
> brk0b = 0.500000000
> brks = 0.666700006
> _______________________________________________________________________________
> MASSES
> e = 5.10998943E-04 nu_e = 0.00000000
> mu = 0.105658367 nu_mu = 0.00000000
> tau = 1.77699006 nu_tau= 9.99999978E-03
> pi = 0.139570177
> pi0 = 0.134976596
> rho = 0.775799990 rho wid.= 0.150299996
> A1 = 1.23000002 A1 wid. = 0.598999977
> K = 0.493676990
> K0 = 0.497648001
> K* = 0.891659975 K* wid. = 5.07999994E-02
> _______________________________________________________________________________
> TAU BRANCHING RATIOS
> nchan = 22
> Tau --> e nu nu 1.00000000
> Tau --> mu nu nu 0.979799986
> Tau --> pi nu 0.649559975
> Tau --> Rho (--> 2pi) nu 1.39901996
> Tau --> A1 (--> 3pi) nu 1.00055003
> Tau --> K nu 3.71199995E-02
> Tau --> K* nu 8.02399963E-02
> Tau --> 3pi pi0 nu 0.312229991
> Tau --> pi 3pi0 nu 0.00000000
> Tau --> 3pi 2pi0 nu 0.00000000
> Tau --> 5pi nu 0.00000000
> Tau --> 5pi pi0 nu 0.00000000
> Tau --> 3pi 3pi0 nu 0.00000000
> Tau --> 2K pi nu 0.00000000
> Tau --> pi K0 K0bar nu 0.00000000
> Tau --> K K0 pi0 nu 0.00000000
> Tau --> K 2pi0 nu 0.00000000
> Tau --> K 2pi nu 0.00000000
> Tau --> pi K0bar pi0 nu 0.00000000
> Tau --> eta pi pi0 nu 0.00000000
> Tau --> pi pi0 gam nu 0.00000000
> Tau --> K K0 nu 0.00000000
> _______________________________________________________________________________
> Starting event loop
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Advisory warning type 3 given after 0 PYEXEC calls:
> (PYRESD:) Failed to decay particle -1000015 with mass 214.740
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
>
> Error type 9 has occured after 0 PYEXEC calls:
> (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
> Failed to read LHEF event information,
> assume end of file has been reached.
> 1********* PYSTAT: Statistics on Number of Events and Cross-sections *********
>
> ==============================================================================
> I I I I
> I Subprocess I Number of points I Sigma I
> I I I I
> I----------------------------------I----------------------------I (mb) I
> I I I I
> I N:o Type I Generated Tried I I
> I I I I
> ==============================================================================
> I I I I
> I 0 All included subprocesses I 0 10000 I 0.000D+00 I
> I 4 User process 1 I 0 10000 I 0.000D+00 I
> I I I I
> ==============================================================================
>
> ********* Total number of errors, excluding junctions = 1111 *************
> ********* Total number of errors, including junctions = 1113 *************
> ********* Total number of warnings = 10001 *************
> ********* Fraction of events that fail fragmentation cuts = 1.00000 *********
>
> Cross section (pb): 0.0000000000000000
>
>
>
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Dear Olivier,

we came up with a solution that works for what we want to do. As you say, there is the possibility of using via madgraph pythia8, which we also tried and it is working. The drawback is that it is more difficult to do the analysis with pgs. As far as it concerns us, we found the answer to the question.

Thanks,

best

Chiara

On 17 Jul 2013, at 22:11, Olivier Mattelaer wrote:

Your question #232276 on MadGraph5 changed:
https://answers.launchpad.net/madgraph5/+question/232276

   Status: Open => Answered

Olivier Mattelaer proposed the following answer:
Hi Chiara,

I will forward this to the FR author, in order to see if they are aware of such problem. Otherwise, you might need to either report to pythia author (or to use pythia8).
The pythia author may suggest to use pythia8 since they want to drop pythia6 anyway.

Cheers,

Olivier

On Jul 11, 2013, at 3:36 PM, Chiara Arina <email address hidden> wrote:

New question #232276 on MadGraph5:
https://answers.launchpad.net/madgraph5/+question/232276

Dear MadGraph team,

I have implemented an extension of the MSSM with feynrules and analyzed some processes specific to the model with madgraph. Since at parton level everything is working fine, I would like to go one step further and simulate the events with pythia-psg. However the pythia-pgs part does not work properly because the model is not implement in pythia. I tried the same process within the MSSM and everything is fine. The same process within the nmssm implemented in madgraph produces the same error as the model I implemented (see below). Is there a way to make pythia-pgs or pythia working with extensions of the mssm?

Thanks,

best

Chiara

Here you can find the the log from the pythia run in the case of the nmssm (the process is direct stau production):

1
******************************************************************************
******************************************************************************
** **
** **
** *......* Welcome to the Lund Monte Carlo! **
** *:::!!:::::::::::* **
** *::::::!!::::::::::::::* PPP Y Y TTTTT H H III A **
** *::::::::!!::::::::::::::::* P P Y Y T H H I A A **
** *:::::::::!!:::::::::::::::::* PPP Y T HHHHH I AAAAA **
** *:::::::::!!:::::::::::::::::* P Y T H H I A A **
** *::::::::!!::::::::::::::::*! P Y T H H III A A **
** *::::::!!::::::::::::::* !! **
** !! *:::!!:::::::::::* !! This is PYTHIA version 6.426 **
** !! !* -><- * !! Last date of change: 16 Nov 2011 **
** !! !! !! **
** !! !! !! Now is 0 Jan 2000 at 0:00:00 **
** !! !! **
** !! lh !! Disclaimer: this program comes **
** !! !! without any guarantees. Beware **
** !! hh !! of errors and use common sense **
** !! ll !! when interpreting results. **
** !! !! **
** !! Copyright T. Sjostrand (2011) **
** **
** An archive of program versions and documentation is found on the web: **
** http://www.thep.lu.se/~torbjorn/Pythia.html **
** **
** When you cite this program, the official reference is to the 6.4 manual: **
** T. Sjostrand, S. Mrenna and P. Skands, JHEP05 (2006) 026 **
** (LU TP 06-13, FERMILAB-PUB-06-052-CD-T) [hep-ph/0603175]. **
** **
** Also remember that the program, to a large extent, represents original **
** physics research. Other publications of special relevance to your **
** studies may therefore deserve separate mention. **
** **
** Main author: Torbjorn Sjostrand; Department of Theoretical Physics, **
** Lund University, Solvegatan 14A, S-223 62 Lund, Sweden; **
** phone: + 46 - 46 - 222 48 16; e-mail: <email address hidden> **
** Author: Stephen Mrenna; Computing Division, GDS Group, **
** Fermi National Accelerator Laboratory, MS 234, Batavia, IL 60510, USA; **
** phone: + 1 - 630 - 840 - 2556; e-mail: <email address hidden> **
** Author: Peter Skands; CERN/PH-TH, CH-1211 Geneva, Switzerland **
** phone: + 41 - 22 - 767 24 47; e-mail: <email address hidden> **
** **
** **
******************************************************************************
******************************************************************************
   MSTU(11) changed from 6 to 6
   MDCY(C15,1) changed from 1 to 0
   MSEL changed from 1 to 0

Reading pythia input card
   mstp(61) changed from 2 to 1
   mstp(71) changed from 1 to 1
   mstj(1) changed from 1 to 1
   mstp(81) changed from 1 to 20
   mstu(21) changed from 2 to 1
Set LHAPATH to /Users/chiara/Documents/software_linux/MadGraph5_v1_5_9/pythia-pgs/src/PDFsets
Opened file unweighted_events.lhe
Running as inclusive sample (highest mult.)

Initiating pythia with external process
1****************** PYINIT: initialization of PYTHIA routines *****************
Reading parameter sminputs 1 127.92000000000000
Reading parameter sminputs 3 0.11720000000000000
Reading model: nmssm
   IMSS(21) changed from 0 to 24
   IMSS(22) changed from 0 to 24
Initializing PYR with random seed 21
   MSTP(52) changed from 1 to 2
Found parameter ickkw 0
   MSTP(98) changed from 0 to 1
==== PYTHIA WILL USE LHAPDF ====
*************************************
* LHAPDF Version 4.0 *
*************************************

PDF description: <<<<<<
CTEQ6L1 - LO with LO alpha_s
Reference:
J. Pumplin, D.R. Stump, J. Huston, H.L. Lai, P. Nadolsky,
W.K. Tung
hep-ph/0201195
                <<<<<<

Parametrization: CTEQ6

warning empty CtLhQCDSET call: NAME=IHDN VALUE= 1.0000000000000000
warning empty CtLhQCDSET call: NAME=QMAX VALUE= 100000.00000000000
warning empty CtLhQCDSET call: NAME=XMIN VALUE= 9.9999999999999995E-007
warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
==============================================
PDFset name /Users/chiara/Documents/software_linux/MadGraph5_v1_5_9/pythia-pgs/src/PDFsets/c
with 1 members
==== initialized. ===========================
warning empty CtLhQCDSET call: NAME=NX VALUE= 95.000000000000000
warning empty CtLhQCDSET call: NAME=NT VALUE= 18.000000000000000
Strong coupling at Mz for PDF is: 0.12978
****************** PYSLHA v1.14: SUSY/BSM SPECTRUM INTERFACE *****************
* (PYSLHA:) Last Change 10 Jun 2010 - P.Z. Skands
* (PYSLHA:) Reading QNUMBERS for KF = 1000022
* (PYSLHA:) KF = 1000022 already exists: ~chi_10 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000023
* (PYSLHA:) KF = 1000023 already exists: ~chi_20 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000025
* (PYSLHA:) KF = 1000025 already exists: ~chi_30 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000035
* (PYSLHA:) KF = 1000035 already exists: ~chi_40 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000045
* (PYSLHA:) Reading QNUMBERS for KF = 1000024
* (PYSLHA:) KF = 1000024 already exists: ~chi_1+ . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000037
* (PYSLHA:) KF = 1000037 already exists: ~chi_2+ . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000021
* (PYSLHA:) KF = 1000021 already exists: ~g . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 35
* (PYSLHA:) KF = 35 already exists: H0 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 45
* (PYSLHA:) KF = 45 already exists: . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 36
* (PYSLHA:) KF = 36 already exists: A0 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 46
* (PYSLHA:) KF = 46 already exists: . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 37
* (PYSLHA:) KF = 37 already exists: H+ . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000012
* (PYSLHA:) KF = 1000012 already exists: ~nu_eL . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000014
* (PYSLHA:) KF = 1000014 already exists: ~nu_muL . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000016
* (PYSLHA:) KF = 1000016 already exists: ~nu_tauL . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000011
* (PYSLHA:) KF = 1000011 already exists: ~e_L- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000013
* (PYSLHA:) KF = 1000013 already exists: ~mu_L- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000015
* (PYSLHA:) KF = 1000015 already exists: ~tau_1- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000011
* (PYSLHA:) KF = 2000011 already exists: ~e_R- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000013
* (PYSLHA:) KF = 2000013 already exists: ~mu_R- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000015
* (PYSLHA:) KF = 2000015 already exists: ~tau_2- . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000002
* (PYSLHA:) KF = 1000002 already exists: ~u_L . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000004
* (PYSLHA:) KF = 1000004 already exists: ~c_L . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000006
* (PYSLHA:) KF = 1000006 already exists: ~t_1 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000002
* (PYSLHA:) KF = 2000002 already exists: ~u_R . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000004
* (PYSLHA:) KF = 2000004 already exists: ~c_R . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000006
* (PYSLHA:) KF = 2000006 already exists: ~t_2 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000001
* (PYSLHA:) KF = 1000001 already exists: ~d_L . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000003
* (PYSLHA:) KF = 1000003 already exists: ~s_L . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 1000005
* (PYSLHA:) KF = 1000005 already exists: ~b_1 . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000001
* (PYSLHA:) KF = 2000001 already exists: ~d_R . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000003
* (PYSLHA:) KF = 2000003 already exists: ~s_R . Entry ignored.
* (PYSLHA:) Reading QNUMBERS for KF = 2000005
* (PYSLHA:) KF = 2000005 already exists: ~b_2 . Entry ignored.
* (PYSLHA:) Ignoring MASS entry for KF = 5 (SLHA read-in not allowed)
* (PYSLHA:) Reading MASS entry for KF = 6, pole mass = 171.400
* (PYSLHA:) Ignoring MASS entry for KF = 15 (SLHA read-in not allowed)
* (PYSLHA:) Reading MASS entry for KF = 23, pole mass = 91.187
* (PYSLHA:) Reading MASS entry for KF = 24, pole mass = 80.939
* (PYSLHA:) Reading MASS entry for KF = 25, pole mass = 119.164
* (PYSLHA:) Reading MASS entry for KF = 35, pole mass = 1016.551
* (PYSLHA:) Reading MASS entry for KF = 36, pole mass = 40.398
* (PYSLHA:) Reading MASS entry for KF = 37, pole mass = 1022.472
* (PYSLHA:) Reading MASS entry for KF = 45, pole mass = 2112.576
* (PYSLHA:) Reading MASS entry for KF = 46, pole mass = 1021.042
* (PYSLHA:) Reading MASS entry for KF = 1000001, pole mass = 1066.519
* (PYSLHA:) Reading MASS entry for KF = 1000002, pole mass = 1063.635
* (PYSLHA:) Reading MASS entry for KF = 1000005, pole mass = 866.365
* (PYSLHA:) Reading MASS entry for KF = 1000006, pole mass = 499.229
* (PYSLHA:) Reading MASS entry for KF = 1000011, pole mass = 380.176
* (PYSLHA:) Reading MASS entry for KF = 1000012, pole mass = 372.121
* (PYSLHA:) Reading MASS entry for KF = 1000015, pole mass = 214.740
* (PYSLHA:) Reading MASS entry for KF = 1000016, pole mass = 360.341
* (PYSLHA:) Reading MASS entry for KF = 1000021, pole mass = 1151.543
* (PYSLHA:) Reading MASS entry for KF = 1000022, pole mass = 208.142
* (PYSLHA:) Reading MASS entry for KF = 1000023, pole mass = 397.851
* (PYSLHA:) Reading MASS entry for KF = 1000024, pole mass = 397.829
* (PYSLHA:) Reading MASS entry for KF = 1000025, pole mass = -963.981
* (PYSLHA:) Reading MASS entry for KF = 1000035, pole mass = 969.594
* (PYSLHA:) Reading MASS entry for KF = 1000037, pole mass = 970.137
* (PYSLHA:) Reading MASS entry for KF = 1000045, pole mass = 2094.274
* (PYSLHA:) Reading MASS entry for KF = 2000001, pole mass = 1023.754
* (PYSLHA:) Reading MASS entry for KF = 2000002, pole mass = 1027.259
* (PYSLHA:) Reading MASS entry for KF = 2000005, pole mass = 992.138
* (PYSLHA:) Reading MASS entry for KF = 2000006, pole mass = 935.527
* (PYSLHA:) Reading MASS entry for KF = 2000011, pole mass = 261.024
* (PYSLHA:) Reading MASS entry for KF = 2000015, pole mass = 374.857
* (PYSLHA:) Ignoring MASS entry for KF = 1 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 2 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 3 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 4 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 11 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 12 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 13 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 14 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 16 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 21 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring MASS entry for KF = 22 (SLHA read-in not allowed)
* (PYSLHA:) Reading MASS entry for KF = 1000014, pole mass = 372.121
* (PYSLHA:) Reading MASS entry for KF = 1000013, pole mass = 380.176
* (PYSLHA:) Reading MASS entry for KF = 2000013, pole mass = 261.024
* (PYSLHA:) Reading MASS entry for KF = 1000004, pole mass = 1063.635
* (PYSLHA:) Reading MASS entry for KF = 2000004, pole mass = 1027.259
* (PYSLHA:) Reading MASS entry for KF = 1000003, pole mass = 1066.519
* (PYSLHA:) Reading MASS entry for KF = 2000003, pole mass = 1023.754
* (PYSLHA:) Allowing DECAY table for KF = 6 but this is NOT recommended.
* (PYSLHA:) Allowing DECAY table for KF = 23 but this is NOT recommended.
* (PYSLHA:) Allowing DECAY table for KF = 24 but this is NOT recommended.
* (PYSLHA:) Ignoring DECAY table for KF = 1 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 2 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 3 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 4 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 5 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 11 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 12 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 13 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 14 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 15 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 16 (SLHA read-in not allowed)
* (PYSLHA:) Ignoring DECAY table for KF = 21 (SLHA read-in not allowed)

        ********************************************************
                STDHEP version 5.04.01 - Aug. 29, 2005
        ********************************************************

STDXWOPEN WARNING: I/O is initialized for stdhep only
* (PYSLHA:) Ignoring DECAY table for KF = 22 (SLHA read-in not allowed)
* (PYSLHA:) Reading SLHA stable particle KF = 1000022, ~chi_10

==============================================================================
I I
I PYTHIA will be initialized for p+ on p+ user configuration I
I with 7000.000 GeV on 7000.000 GeV beam energies I
I I
I corresponding to 14000.000 GeV center-of-mass energy I
I I
==============================================================================

******** PYMAXI: summary of differential cross-section maximum search ********

         ==========================================================
         I I I
         I ISUB Subprocess name I Maximum value I
         I I I
         ==========================================================
         I I I
         I 4 User process 1 I 5.5237D-12 I
         I I I
         ==========================================================

********************** PYINIT: initialization completed **********************
Initialisation done

Generating all events
stdhep initialized

***************************************************************************
* *****TAUOLA LIBRARY: VERSION 2.5 ****** *
* ***********JUNE 1994*************** *
* **AUTHORS: S.JADACH, Z.WAS************* *
* **R. DECKER, M. JEZABEK, J.H.KUEHN***** *
* **AVAILABLE FROM: WASM AT CERNVM ****** *
* ***** PUBLISHED IN COMP. PHYS. COMM.*** *
* *******CERN-TH-5856 SEPTEMBER 1990***** *
* *******CERN-TH-6195 SEPTEMBER 1991***** *
* *******CERN-TH-6793 NOVEMBER 1992***** *
* **5 or more pi dec.: precision limited *
* ******DEXAY ROUTINE: INITIALIZATION**** *
* 0 JAK1 = DECAY MODE FERMION1 (TAU+) *
* 0 JAK2 = DECAY MODE FERMION2 (TAU-) *
***************************************************************************

***************************************************************************

_______________________________________________________________________________
TAUOLA VERSION 2.1 PARAMETER SET 0
_______________________________________________________________________________
Printing end run statistics.
GENP bank used for input.
Mult. taus in event have same chirality.
Taus polarized in W/Z rest frame.
Helicity: 1.00000000
Lifetime: 295.600006 fs
_______________________________________________________________________________
STANDARD MODEL PARAMETERS
Gfermi = 1.16639203E-05
Gv = 1.00000000
Ga = -1.00000000
Cos(Cabibbo_angle) = 0.975000024
Sin(Cabibbo_angle) = 0.222204715
Gamma_e = 4.04914356E-13
1/alpha = 137.03598956100001
1/(alpha*pi) = 2.3228196271906998E-003
xk0 = 0.10000000149011612
_______________________________________________________________________________
KAON/A1 DECAY PARAMETERS
bra1 = 0.500000000
brk0 = 0.500000000
brk0b = 0.500000000
brks = 0.666700006
_______________________________________________________________________________
MASSES
e = 5.10998943E-04 nu_e = 0.00000000
mu = 0.105658367 nu_mu = 0.00000000
tau = 1.77699006 nu_tau= 9.99999978E-03
pi = 0.139570177
pi0 = 0.134976596
rho = 0.775799990 rho wid.= 0.150299996
A1 = 1.23000002 A1 wid. = 0.598999977
K = 0.493676990
K0 = 0.497648001
K* = 0.891659975 K* wid. = 5.07999994E-02
_______________________________________________________________________________
TAU BRANCHING RATIOS
nchan = 22
Tau --> e nu nu 1.00000000
Tau --> mu nu nu 0.979799986
Tau --> pi nu 0.649559975
Tau --> Rho (--> 2pi) nu 1.39901996
Tau --> A1 (--> 3pi) nu 1.00055003
Tau --> K nu 3.71199995E-02
Tau --> K* nu 8.02399963E-02
Tau --> 3pi pi0 nu 0.312229991
Tau --> pi 3pi0 nu 0.00000000
Tau --> 3pi 2pi0 nu 0.00000000
Tau --> 5pi nu 0.00000000
Tau --> 5pi pi0 nu 0.00000000
Tau --> 3pi 3pi0 nu 0.00000000
Tau --> 2K pi nu 0.00000000
Tau --> pi K0 K0bar nu 0.00000000
Tau --> K K0 pi0 nu 0.00000000
Tau --> K 2pi0 nu 0.00000000
Tau --> K 2pi nu 0.00000000
Tau --> pi K0bar pi0 nu 0.00000000
Tau --> eta pi pi0 nu 0.00000000
Tau --> pi pi0 gam nu 0.00000000
Tau --> K K0 nu 0.00000000
_______________________________________________________________________________
Starting event loop

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Advisory warning type 3 given after 0 PYEXEC calls:
   (PYRESD:) Failed to decay particle -1000015 with mass 214.740

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.

   Error type 9 has occured after 0 PYEXEC calls:
   (PYEVNW:) failed to evolve shower or multiple interactions. Returning.
Failed to read LHEF event information,
assume end of file has been reached.
1********* PYSTAT: Statistics on Number of Events and Cross-sections *********

==============================================================================
I I I I
I Subprocess I Number of points I Sigma I
I I I I
I----------------------------------I----------------------------I (mb) I
I I I I
I N:o Type I Generated Tried I I
I I I I
==============================================================================
I I I I
I 0 All included subprocesses I 0 10000 I 0.000D+00 I
I 4 User process 1 I 0 10000 I 0.000D+00 I
I I I I
==============================================================================

********* Total number of errors, excluding junctions = 1111 *************
********* Total number of errors, including junctions = 1113 *************
********* Total number of warnings = 10001 *************
********* Fraction of events that fail fragmentation cuts = 1.00000 *********

Cross section (pb): 0.0000000000000000

--
You received this question notification because you are a member of
MadTeam, which is an answer contact for MadGraph5.

--
If this answers your question, please go to the following page to let us
know that it is solved:
https://answers.launchpad.net/madgraph5/+question/232276/+confirm?answer_id=0

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/madgraph5/+question/232276

You received this question notification because you asked the question.

Hi Chiara,

Could you describe your solution? This might help me to answer future question and/or help other people with similar question.

Thanks,

Olivier

Thanks Benjamin Fuks, that solved my question.

Hi Olivier,

we could use pythia8, as suggested, however we would loose some nice features of madgraph of having more than two particles in the final state and also the pgs part (or at least it would more complicated to get the same level of analysis). For our analysis, the whole matrix element is computed at parton level with madgraph, after which we are left with SM particles plus a stable particle which is supersymmetric (however not the neutralino but the sneutrino). Even if we extend the MSSM with right-handed neutrinos, the sneutrino is already present in the MSSM. Therefore to do hadronization of the SM particles only, we forced pythia to load our model as it was the MSSM and to recognize the sneutrino as stable, by modifying a file in the src part of pythia. If you want to know more details we can send you the modified file.

best

Chiara

------------------------------------------
Chiara Arina
Postdoc researcher
Institute for Theoretical Physics
University of Amsterdam

On 17 Jul 2013, at 23:21, Olivier Mattelaer wrote:

Your question #232276 on MadGraph5 changed:
https://answers.launchpad.net/madgraph5/+question/232276

Olivier Mattelaer posted a new comment:
Hi Chiara,

Could you describe your solution? This might help me to answer future
question and/or help other people with similar question.

Thanks,

Olivier

--
You received this question notification because you asked the question.

Hi Chiara,

My two cents: maybe could you use delphes for the fastsim part? It can read the HepMC output of pythia8 out of the box.

Cheers,

Benjamin

Hi Benjamin,

thanks for your advise; in this particular case we found an alternative solution within pythia-pgs, as I told Olivier. However we weren't aware of the possibility you mention and we will keep it in mind for the future.

best

Chiara

On 18 Jul 2013, at 17:06, Benjamin Fuks wrote:

Your question #232276 on MadGraph5 changed:
https://answers.launchpad.net/madgraph5/+question/232276

Benjamin Fuks posted a new comment:
Hi Chiara,

My two cents: maybe could you use delphes for the fastsim part? It can
read the HepMC output of pythia8 out of the box.

Cheers,

Benjamin

--
You received this question notification because you asked the question.

Chiara,

You say that you had madgraph calculate the "whole" matrix element so that pythia only showers on SM particles. I am having a very similar error when generating color octet channels for quarkonium. How did you modify proc_card.dat (specifically the process and model you used) to get Madgraph to leave only SM particles at the end but have intermediate states of non-SM particles?

Color octet channels supposedly are used as follows:

gg>gcc~[3S18to9900443]

But I get a similar "Failed to decay..." message

Thank you.