problem with pythia 8: Fail to produce pythia output
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! **
** *:::!!:::::::::::* **
** *::::::
** *::::::
** *::::::
** *::::::
** *::::::
** *::::::
** !! *:::!!:::::::::::* !! 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://
** **
** 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-
** **
** 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/
Set LHAPATH to /amnt/remote/
Set LHAPATH to /amnt/remote/
Opened file unweighted_
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.9999999999999
warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
======
PDFset name /amnt/remote/
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
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
DECAY 36 3.10965379e-02 # h3 decays, CP-odd Higgs
# BR NDA ID1 ID2
DECAY 37 1.05271876e-01 # Charged Higgs decays
# BR NDA ID1 ID2
#
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/
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/
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/
(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
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Tania Robens for more information if necessary.