problem with pythia 8: Fail to produce pythia output

Asked by Tania Robens

Hi -

I am trying to run pythia from the command line (choosing pythia=on), and get the following output:

INFO: Running Pythia
Fail to produce pythia output.

with the tag card

 ******************************************************************************
 ******************************************************************************
 ** **
 ** **
 ** *......* 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)=1
     mstp(61) changed from 2 to 1
       mstp(71)=1
     mstp(71) changed from 1 to 1

       mstj(1)=1
     mstj(1) changed from 1 to 1

       mstp(81)=20
     mstp(81) changed from 1 to 20

       mstu(21)=1
     mstu(21) changed from 2 to 1

       lhaid= 10041
 Read LHAID = 10041

 Set LHAPATH to /amnt/remote/pkts08.raid1/users/trobens/programs/madgraph/MG5_aMC_v2_1_2/pythia-pgs/src/PDFsets

 Set LHAPATH to /amnt/remote/pkts08.raid1/users/trobens/programs/madgraph/MG5_aMC_v2_1_2/pythia-pgs/src/PDFsets

 Set LHAPATH to /amnt/remote/pkts08.raid1/users/trobens/programs/madgraph/MG5_aMC_v2_1_2/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 model: InertDou
     IMSS(21) changed from 0 to 24
     IMSS(22) changed from 0 to 24
 Initializing PYR with random seed 48
     MSTP(52) changed from 1 to 2
 Found parameter ickkw 0
     MSTP(98) changed from 0 to 1
 ==== PYTHIA WILL USE LHAPDF ====
 WRONG LHAPDF set number = 200400 given! STOP EXE!
 WILL USE 10041 INSTEAD
At line 467 of file pyslha.f
Fortran runtime error: Bad real number in item 2 of list input
 *************************************
 * LHAPDF Version 4.0 *
 *************************************

 >>>>>> PDF description: <<<<<<
 CTEQ6L - LO with NLO 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 /amnt/remote/pkts08.raid1/users/trobens/programs/madgraph/MG5_aMC_v2_1_2/pythia-
 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.11798
 ****************** PYSLHA v1.14: SUSY/BSM SPECTRUM INTERFACE *****************
 * (PYSLHA:) Last Change 10 Jun 2010 - P.Z. Skands
 * (PYSLHA:) Reading QNUMBERS for KF = 35
 * (PYSLHA:) KF = 35 already exists: H0 . Entry ignored.
 * (PYSLHA:) Reading QNUMBERS for KF = 36
 * (PYSLHA:) KF = 36 already exists: A0 . Entry ignored.
 * (PYSLHA:) Reading QNUMBERS for KF = 37
 * (PYSLHA:) KF = 37 already exists: H+ . Entry ignored.
 * (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 = 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 = 15 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring MASS entry for KF = 16 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring MASS entry for KF = 4 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring MASS entry for KF = 5 (SLHA read-in not allowed)
 * (PYSLHA:) Reading MASS entry for KF = 6, pole mass = 175.000
 * (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 = 23, pole mass = 91.188
 * (PYSLHA:) Reading MASS entry for KF = 24, pole mass = 79.947
 * (PYSLHA:) Reading MASS entry for KF = 25, pole mass = 150.000
 * (PYSLHA:) Reading MASS entry for KF = 35, pole mass = 40.000

I am not so sure what the problem is... as you can see i run it w a 2hdm parameter card, running it per default (eg for w+ w- in the standard model or so) does not produce this, so i dont think its pdf related eg. ??

the parameter card is given by

Block frblock
1 7.757952e-03 # aEWM1
2 1.166370e-05 # Gf
3 9.118760e+01 # MMZ
4 7.994700e+01 # MMW
5 1.200000e+00 # MMC
6 4.230000e+00 # MMB
7 1.750000e+02 # MMT
8 1.172000e-01 # aS
9 1.000000e+02 # QS
10 -0.06875 # lamL
11 0.1 # lam2
12 150 # mmh
13 40 # mmH0
14 140 # mmA0
15 140 # mmHch

Block mass
 1 0.000000e+00 # MD
   2 0.000000e+00 # MU
  3 1.000000e-01 # MS
 11 0.000000e+00 # Me
 12 0.000000e+00 # Mnue
 13 1.057000e-01 # MM
 14 0.000000e+00 # Mnum
 15 1.777000e+00 # MTA
 16 0.000000e+00 # Mnut
 4 1.200000 # c : MMC
 5 4.230000 # b : MMB
 6 175.000000 # t : MMT
 21 0.000000 # g : 0.0
 22 0.000000 # a : 0.0
 23 91.187600 # z : MMZ
 24 79.947000 # w+ : MMW
25 150 # h : mmh
35 40 # h2 : mmH0
36 140# h3 : mmA0
37 140 # h+ : mmHch

# PDG Width
DECAY 6 1.47419196e+00 # top decays
# BR NDA ID1 ID2
       1.00000000e+00 2 24 5
DECAY 35 0.00000000e+00 # h1 decays, lightest CP-even Higgs
# BR NDA ID1 ID2
DECAY 25 8.10412825e-02 # h2 decays, heaviest CP-even Higgs
# BR NDA ID1 ID2
       4.49447644e-05 2 3 -3
       1.40601505e-03 2 4 -4
       3.48524670e-02 2 5 -5
       3.17202348e-10 2 11 -11
       1.35613501e-05 2 13 -13
       3.83201622e-03 2 15 -15
       2.69093035e-04 2 22 22
       1.71565885e-02 2 23 23
       1.41871946e-01 2 24 -24
       4.79299358e-04 2 23 22
       6.73109068e-03 2 21 21
       7.93342978e-01 2 25 25
DECAY 36 3.10965379e-02 # h3 decays, CP-odd Higgs
# BR NDA ID1 ID2
       1.00000000e+00 2 23 25
DECAY 37 1.05271876e-01 # Charged Higgs decays
# BR NDA ID1 ID2
       1.00000000e+00 2 24 25
#

Block QNUMBERS 35 # h2
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)
Block QNUMBERS 36 # h3
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)
Block QNUMBERS 37 # h+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/Antiparticle distinction (0=own anti)

 (i initially thought its a problem about decay specification...)

please get back to me concerning this

thanks best tania

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Rikkert Frederix (frederix) said :
#1

Dear Tania,

Could you please send me a small sample event file and the pythia_card.dat for which you get the failure? If you cannot upload them here, you can also send them to me by e-mail, rikkert.frederix_@_cern.ch (remove the underscores).

Best regards,
Rikkert

Revision history for this message
Rikkert Frederix (frederix) said :
#2

Ok. I probably found it.

You need to have a space between the values for the parameters and the comments in the param_card.dat. In particular the line that specifies the mass of particle 36 (the 'h3') reads:

36 140# h3 : mmA0

but you need a space between the '140' and the hash '#':

36 140 # h3 : mmA0

Please, let me know if this solves the problem.

Cheers,
Rik

Can you help with this problem?

Provide an answer of your own, or ask Tania Robens for more information if necessary.

To post a message you must log in.