Error detected in sub-command generate_events

Hi Daniel,
thanks for the feedback, but I would need a bit more information.
What do you mean exactly by 'without any response'?
You do 'generate events', edit the cards, press enter, and what happens?
Are the various jobs running, and then they get stuck at some point, or they
don't even start? If we are in the former situation, what is the stage at which
they get stuck? While setting up grids, computing upper envelope, or generating
events? What is the process / model you are using? How many events or what
precision do you require?
Cheers.
Paolo

Hi Paolo,

ok, when I press enter after to do 'generate events' MG5 calculates the cross-section at parton level, in the terminal reads,
.
.
.
End Parton
Running Pythia

and continue in this step for the eternity! until I stop them with the above response "Error detected in sub-command generate_events"..... The model that I use here is own model that contain new particles (I import the SM model and the Effective_Higgs_Couplings models into the model). This behavior of MG5 ist for some processes (pp > U U~) in particular with the new quarks. The quantity of events is 10000. If you want I can send you the folders. Thanks.

Cheers.

daniel

Hi Daniel,
just to understand better, you are running LO MadGraph5, right?
Have you tried with just 10 events? Some times if you run Pythia
with Underlying Event on, it may take a while to shower the full
sample, while with 10 events it should take few seconds.
If you get the same problem even in that case, it seems that Pythia
either is caught in some loop, or waiting for more input from the
user (I don’t really know). Could you try and use your input Pythia
card for some other more standard process, and see if you get the
same problem or not?
Thanks.
Cheers.
Paolo

On 11 Jun 2014, at 23:21, daniel camargo <email address hidden> wrote:

> Question #250051 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/250051
>
> Status: Needs information => Open
>
> daniel camargo gave more information on the question:
> Hi Paolo,
>
> ok, when I press enter after to do 'generate events' MG5 calculates the cross-section at parton level, in the terminal reads,
> .
> .
> .
> End Parton
> Running Pythia
>
> and continue in this step for the eternity! until I stop them with the
> above response "Error detected in sub-command generate_events"..... The
> model that I use here is own model that contain new particles (I import
> the SM model and the Effective_Higgs_Couplings models into the model).
> This behavior of MG5 ist for some processes (pp > U U~) in particular
> with the new quarks. The quantity of events is 10000. If you want I can
> send you the folders. Thanks.
>
> Cheers.
>
> daniel
>
> --
> You received this question notification because you are a direct
> subscriber of the question.

Hi Paolo,

I done it and continue with the same problem... here the run...

daniel@daniel:~/Documents/MadGraph5/CSU13$ bin/newprocess_mg5
Running...
done
daniel@daniel:~/Documents/MadGraph5/CSU13$ bin/generate_events
No module named madgraph.interface.extended_cmd
No module named madgraph
No module named madgraph
************************************************************
* *
* W E L C O M E to M A D G R A P H 5 *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 5.1.5.5 *
* *
* The MadGraph Development Team - Please visit us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
load configuration from /home/daniel/Documents/MadGraph5/CSU13/Cards/me5_configuration.txt
load configuration from /home/daniel/Documents/MadGraph5/input/mg5_configuration.txt
load configuration from /home/daniel/Documents/MadGraph5/CSU13/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
generate_events
Which programs do you want to run?
  0 / auto : running existing card
  1 / parton : Madevent
  2 / pythia : MadEvent + Pythia.
  3 / pgs : MadEvent + Pythia + PGS.
 [0, 1, 2, 3, auto, parton, pythia, pgs][60s to answer]
>0
Will run in mode pgs
Do you want to edit one cards (press enter to bypass editing)?
  1 / param : param_card.dat (be carefull about parameter consistency, especially widths)
  2 / run : run_card.dat
  3 / pythia : pythia_card.dat
  4 / pgs : pgs_card.dat
  9 / plot : plot_card.dat
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
 [0, done, 1, param, 2, run, 3, pythia, 4, enter path, ... ][60s to answer]
>0
Generating 10 events with run name run_05
survey run_05
compile directory
/home/daniel/Documents/MadGraph5/CSU13/bin/internal/ufomodel/particles.py
write ./param_card.dat
Using random number seed offset = 57
Running Survey
Creating Jobs
Working on SubProcesses
    P1_gg_u1u1x
    P1_uux_u1u1x
    P1_ccx_u1u1x
    P1_ddx_u1u1x
    P1_ssx_u1u1x
 Idle: 0 Running: 0 Finish: 6
End survey
refine 10
Creating Jobs
Refine results to 10
    P1_gg_u1u1x
    P1_uux_u1u1x
    P1_ccx_u1u1x
    P1_ddx_u1u1x
    P1_ssx_u1u1x
Combining runs
finish refine
refine 10
Creating Jobs
Refine results to 10
    P1_gg_u1u1x
    P1_uux_u1u1x
    P1_ccx_u1u1x
    P1_ddx_u1u1x
    P1_ssx_u1u1x
Combining runs
finish refine
combine_events
Combining Events
Creating root files
** Calculating number of events to process. Please wait...
** Input file contains 10 events
[#########################] (100.00%) : 10/10 entries processed
** Exiting...
  === Results Summary for run: run_05 tag: tag_1 ===

     Cross-section : 33.02 +- 0.2264 pb
     Nb of events : 10

Creating Plots for parton level
End Plots for parton level
store_events
Storing parton level results
End Parton
Running Pythia

And continue runing.....I done also the pp> t t~ process and the problem disappear .

Cheers.

daniel

Hi Daniel,
since for pp>tt there is no problem, could you write in detail what
are the shell commands you have used to generate the process? I mean
in particular the model you import, the (multi)particles you define,
the process you generate, and so on.
After that, I’ll try to redo exactly what you do and reproduce your
situation.
Cheers.
Paolo

On 12 Jun 2014, at 00:31, daniel camargo <email address hidden> wrote:

> Question #250051 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/250051
>
> Status: Answered => Open
>
> daniel camargo is still having a problem:
> Hi Paolo,
>
> I done it and continue with the same problem... here the run...
>
> daniel@daniel:~/Documents/MadGraph5/CSU13$ bin/newprocess_mg5
> Running...
> done
> daniel@daniel:~/Documents/MadGraph5/CSU13$ bin/generate_events
> No module named madgraph.interface.extended_cmd
> No module named madgraph
> No module named madgraph
> ************************************************************
> * *
> * W E L C O M E to M A D G R A P H 5 *
> * M A D E V E N T *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 5.1.5.5 *
> * *
> * The MadGraph Development Team - Please visit us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * *
> * Type 'help' for in-line help. *
> * *
> ************************************************************
> load configuration from /home/daniel/Documents/MadGraph5/CSU13/Cards/me5_configuration.txt
> load configuration from /home/daniel/Documents/MadGraph5/input/mg5_configuration.txt
> load configuration from /home/daniel/Documents/MadGraph5/CSU13/Cards/me5_configuration.txt
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
> Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> generate_events
> Which programs do you want to run?
> 0 / auto : running existing card
> 1 / parton : Madevent
> 2 / pythia : MadEvent + Pythia.
> 3 / pgs : MadEvent + Pythia + PGS.
> [0, 1, 2, 3, auto, parton, pythia, pgs][60s to answer]
>> 0
> Will run in mode pgs
> Do you want to edit one cards (press enter to bypass editing)?
> 1 / param : param_card.dat (be carefull about parameter consistency, especially widths)
> 2 / run : run_card.dat
> 3 / pythia : pythia_card.dat
> 4 / pgs : pgs_card.dat
> 9 / plot : plot_card.dat
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> [0, done, 1, param, 2, run, 3, pythia, 4, enter path, ... ][60s to answer]
>> 0
> Generating 10 events with run name run_05
> survey run_05
> compile directory
> /home/daniel/Documents/MadGraph5/CSU13/bin/internal/ufomodel/particles.py
> write ./param_card.dat
> Using random number seed offset = 57
> Running Survey
> Creating Jobs
> Working on SubProcesses
> P1_gg_u1u1x
> P1_uux_u1u1x
> P1_ccx_u1u1x
> P1_ddx_u1u1x
> P1_ssx_u1u1x
> Idle: 0 Running: 0 Finish: 6
> End survey
> refine 10
> Creating Jobs
> Refine results to 10
> P1_gg_u1u1x
> P1_uux_u1u1x
> P1_ccx_u1u1x
> P1_ddx_u1u1x
> P1_ssx_u1u1x
> Combining runs
> finish refine
> refine 10
> Creating Jobs
> Refine results to 10
> P1_gg_u1u1x
> P1_uux_u1u1x
> P1_ccx_u1u1x
> P1_ddx_u1u1x
> P1_ssx_u1u1x
> Combining runs
> finish refine
> combine_events
> Combining Events
> Creating root files
> ** Calculating number of events to process. Please wait...
> ** Input file contains 10 events
> [#########################] (100.00%) : 10/10 entries processed
> ** Exiting...
> === Results Summary for run: run_05 tag: tag_1 ===
>
> Cross-section : 33.02 +- 0.2264 pb
> Nb of events : 10
>
> Creating Plots for parton level
> End Plots for parton level
> store_events
> Storing parton level results
> End Parton
> Running Pythia
>
>
> And continue runing.....I done also the pp> t t~ process and the problem
> disappear .
>
> Cheers.
>
> daniel
>
> --
> You received this question notification because you are a direct
> subscriber of the question.

Hi Paolo,

I create a copy of the "Template" folder and then import my model in the "proc_card_mg5.dat" .

I import a model that contains the SM model, the Higgs_Effective_Couplings_UFO model and my own model Exotics_UFO.
The multiparticles here are:

Defined multiparticle all = g a ve vm vt ve~ vm~ vt~ u c t d s b u1 x1 y1 d1 u~ c~ t~ d~ s~ b~ u1~ x1~ y1~ d1~ z w+ h g0 g+ hg h1 w- g- e- mu- ta- e1 n1 e2 f1 e+ mu+ ta+ e1~ n1~ e2~ f1~

where hg its a sm-like Higgs particle and y1,x1,u1,d1 are quarks some of the with exotic charges.

after to do copies of the pythia_card_default.dat and pgs_card_default.dat I generete the matrix element with the command "bin/newprocess_mg5" after of that I do "bin/generate_events"

Above the copy of the "proc_card_mg5.dat" that I used.

#************************************************************
#* MadGraph 5 *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 1.5.5 2012-11-18 *
#* *
#* The MadGraph Development Team - Please visit us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph 5 *
#* *
#* run as ./bin/mg5 filename *
#* *
#************************************************************

set group_subprocesses Auto
set ignore_six_quark_processes False
set gauge unitary
set complex_mass_scheme False
import model Higgs_Effective_Couplings_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 multiparticle labels
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~
# Specify process(es) to run
generate p p > U1 U1~ @1
# Output processes to MadEvent directory
output -f

Cheers

daniel

Dear Daniel,

could you first try to run completely from the shell, instead of copying the folder?
I mean: go in the folder above Template, do ./bin/mg5_aMC, and then type the
commands within the shell, namely 'import model ...', 'generate ...', 'output proc',
'launch'. Then you would type '1' if you want pythia, and run (maybe better to try
with the default pythia card to be sure).

If this gives the same problem, could you then look (before and after interrupting
the run) at the file proc/Events/run_**/tag_*_pythia.log ('proc' is a fake name: the
true name is what you typed after 'output')? Is it there? What does it say?

Then, after interrupting, could you go in proc, do './bin/madevent', type 'pythia' in
the python shell and see what happens?

Thanks.
Cheers.
Paolo

Hi Paolo,

I done that you said me and yes, the pythia continue with the problem...so before the stopped of the run the "tag_1_pythia.log" file say the following:

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 /home/daniel/Documents/MadGraph5/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.90000000000001
 Reading parameter sminputs 2 1.16636999999999998E-005
 Reading parameter sminputs 3 0.11840000000000001
 Reading model: Higgs_Ef
     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.99999999999999955E-007
 warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
 warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
 ==============================================
 PDFset name /home/daniel/Documents/MadGraph5/pythia-pgs/src/PDFsets/cteq6ll.LHpdf
 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 = 250
 * (PYSLHA:) Reading QNUMBERS for KF = 251
 * (PYSLHA:) Reading QNUMBERS for KF = 9000006
 * (PYSLHA:) Reading QNUMBERS for KF = 9000007
 * (PYSLHA:) Reading QNUMBERS for KF = 9000008
 * (PYSLHA:) Reading QNUMBERS for KF = 9000009
 * (PYSLHA:) Reading QNUMBERS for KF = 9000010
 * (PYSLHA:) Reading QNUMBERS for KF = 9000011
 * (PYSLHA:) Reading QNUMBERS for KF = 9000012
 * (PYSLHA:) Reading QNUMBERS for KF = 9000013
 * (PYSLHA:) Reading QNUMBERS for KF = 9000014
 * (PYSLHA:) Reading QNUMBERS for KF = 9000015
 * (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 = 5 (SLHA read-in not allowed)
 * (PYSLHA:) Reading MASS entry for KF = 6, pole mass = 750.000

     Advisory warning type 9 given after 0 PYEXEC calls:
     (PYSLHA:) Note top quark mass, M = 750.00000GeV
 * (PYSLHA:) Ignoring MASS entry for KF = 11 (SLHA read-in not allowed)

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

 STDXWOPEN WARNING: I/O is initialized for stdhep only
 * (PYSLHA:) Ignoring MASS entry for KF = 13 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring MASS entry for KF = 15 (SLHA read-in not allowed)
 * (PYSLHA:) Reading MASS entry for KF = 23, pole mass = 91.188
 * (PYSLHA:) Reading MASS entry for KF = 25, pole mass = 120.000
 * (PYSLHA:) Reading MASS entry for KF = 9000006, pole mass = 300.000
 * (PYSLHA:) Reading MASS entry for KF = 9000007, pole mass = 3000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000008, pole mass = 3000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000009, pole mass = 3000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000010, pole mass = 1000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000011, pole mass = 1000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000012, pole mass = 1000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000013, pole mass = 1000.000
 * (PYSLHA:) Reading MASS entry for KF = 9000014, pole mass = 400.000
 * (PYSLHA:) Reading MASS entry for KF = 9000015, pole mass = 120.000
 * (PYSLHA:) Ignoring MASS entry for KF = 12 (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 = 24, pole mass = 79.824
 * (PYSLHA:) Ignoring MASS entry for KF = 250 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring MASS entry for KF = 251 (SLHA read-in not allowed)
 * (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)
 * (PYSLHA:) Ignoring DECAY table for KF = 22 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring DECAY table for KF = 250 (SLHA read-in not allowed)
 * (PYSLHA:) Ignoring DECAY table for KF = 251 (SLHA read-in not allowed)

 ==============================================================================
 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 0 I 3.7285D-08 I
           I I I
           ==========================================================

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

 Generating all events
 stdhep initialized

The file "tag_1_pythia.log" remains equal after to stop pythia.
Now I done the other run through the proc file and this happened:

daniel@daniel:~/Documents/MadGraph5$ cd PROC_Higgs_Effective_Couplings_UFO_3/
daniel@daniel:~/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3$ ./bin/madevent
launch in debug mode
************************************************************
* *
* W E L C O M E to M A D G R A P H 5 *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 5.1.5.5 *
* *
* The MadGraph Development Team - Please visit us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
load configuration from /home/daniel/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3/Cards/me5_configuration.txt
load configuration from /home/daniel/Documents/MadGraph5/input/mg5_configuration.txt
load configuration from /home/daniel/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
MGME5>pythia
Which programs do you want to run?
    0 / auto : running existing card
    1 / pythia : Pythia
    2 / pgs : Pythia + PGS
 [0, 1, 2, auto, pythia, pgs][60s to answer]
>1
Will run in mode pythia
Do you want to edit one cards (press enter to bypass editing)?
  1 / pythia : pythia_card.dat
  9 / plot : plot_card.dat
  Path to a valid card.
 [0, done, 1, pythia, 9, plot, enter path][60s to answer]
>0
compile directory
Running Pythia

and remains so for a long time... after to stop the run show me the following;

^CTraceback (most recent call last):
  File "./bin/madevent", line 83, in <module>
    subprocess.call([sys.executable] + ['-O'] + sys.argv)
  File "/usr/lib/python2.7/subprocess.py", line 493, in call
Stop by the user
stopping all operation
            in order to quit madevent please enter exit
pythia Stop by the user
MGME5> return Popen(*popenargs, **kwargs).wait()
  File "/usr/lib/python2.7/subprocess.py", line 1291, in wait
    pid, sts = _eintr_retry_call(os.waitpid, self.pid, 0)
  File "/usr/lib/python2.7/subprocess.py", line 478, in _eintr_retry_call
    return func(*args)
KeyboardInterrupt

thanks Paolo for the patience with this problem!

Cheers

daniel

Hi Daniel,
it then seems to be a problem linked to the pdg codes of your particles.

The last line printed out by the code is ’stdhep initialized’, namely line
286 of file pythia-pgs/src/pythia.f. The code then calls the routine pdgrdtb
at line 289.
The latter routine is in pythia-pgs/libraries/PGS4/src/stdhep-dir/src/stdhep/pdgrdtb.F.

I cannot tell what causes the problem right now, since I don’t know this part
of the code very well. However, I would suggest you to go in that code, insert
some write(*,*) statements at various stages of that routine, and recompile and
run, so to see at which exact point the codes gets caught into an infinite loop,
and why.

In general I may imagine something fishy is happening in stdhep about the
recognition of the pdg code of your non-standard particles: are the conventions
you’re using for these pdg codes accepted by stdhep (see for example the website
http://cd-docdb.fnal.gov/cgi-bin/RetrieveFile?docid=903&version=14&filename=stdhep_505_manual.pdf)?
What happens if you change those pdg codes?

Cheers.
Paolo

On 14 Jun 2014, at 04:47, daniel camargo <email address hidden> wrote:

> Question #250051 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/250051
>
> Status: Answered => Open
>
> daniel camargo is still having a problem:
> Hi Paolo,
>
> I done that you said me and yes, the pythia continue with the
> problem...so before the stopped of the run the "tag_1_pythia.log" file
> say the following:
>
> 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 /home/daniel/Documents/MadGraph5/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.90000000000001
> Reading parameter sminputs 2 1.16636999999999998E-005
> Reading parameter sminputs 3 0.11840000000000001
> Reading model: Higgs_Ef
> 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.99999999999999955E-007
> warning empty CtLhQCDSET call: NAME=QINI VALUE= 1.3000000000000000
> warning empty CtLhQCDSET call: NAME=IKNL VALUE= 1.0000000000000000
> ==============================================
> PDFset name /home/daniel/Documents/MadGraph5/pythia-pgs/src/PDFsets/cteq6ll.LHpdf
> 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 = 250
> * (PYSLHA:) Reading QNUMBERS for KF = 251
> * (PYSLHA:) Reading QNUMBERS for KF = 9000006
> * (PYSLHA:) Reading QNUMBERS for KF = 9000007
> * (PYSLHA:) Reading QNUMBERS for KF = 9000008
> * (PYSLHA:) Reading QNUMBERS for KF = 9000009
> * (PYSLHA:) Reading QNUMBERS for KF = 9000010
> * (PYSLHA:) Reading QNUMBERS for KF = 9000011
> * (PYSLHA:) Reading QNUMBERS for KF = 9000012
> * (PYSLHA:) Reading QNUMBERS for KF = 9000013
> * (PYSLHA:) Reading QNUMBERS for KF = 9000014
> * (PYSLHA:) Reading QNUMBERS for KF = 9000015
> * (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 = 5 (SLHA read-in not allowed)
> * (PYSLHA:) Reading MASS entry for KF = 6, pole mass = 750.000
>
> Advisory warning type 9 given after 0 PYEXEC calls:
> (PYSLHA:) Note top quark mass, M = 750.00000GeV
> * (PYSLHA:) Ignoring MASS entry for KF = 11 (SLHA read-in not allowed)
>
>
> ********************************************************
> STDHEP version 5.04.01 - Aug. 29, 2005
> ********************************************************
>
>
> STDXWOPEN WARNING: I/O is initialized for stdhep only
> * (PYSLHA:) Ignoring MASS entry for KF = 13 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 15 (SLHA read-in not allowed)
> * (PYSLHA:) Reading MASS entry for KF = 23, pole mass = 91.188
> * (PYSLHA:) Reading MASS entry for KF = 25, pole mass = 120.000
> * (PYSLHA:) Reading MASS entry for KF = 9000006, pole mass = 300.000
> * (PYSLHA:) Reading MASS entry for KF = 9000007, pole mass = 3000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000008, pole mass = 3000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000009, pole mass = 3000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000010, pole mass = 1000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000011, pole mass = 1000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000012, pole mass = 1000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000013, pole mass = 1000.000
> * (PYSLHA:) Reading MASS entry for KF = 9000014, pole mass = 400.000
> * (PYSLHA:) Reading MASS entry for KF = 9000015, pole mass = 120.000
> * (PYSLHA:) Ignoring MASS entry for KF = 12 (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 = 24, pole mass = 79.824
> * (PYSLHA:) Ignoring MASS entry for KF = 250 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring MASS entry for KF = 251 (SLHA read-in not allowed)
> * (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)
> * (PYSLHA:) Ignoring DECAY table for KF = 22 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 250 (SLHA read-in not allowed)
> * (PYSLHA:) Ignoring DECAY table for KF = 251 (SLHA read-in not allowed)
>
> ==============================================================================
> 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 0 I 3.7285D-08 I
> I I I
> ==========================================================
>
> ********************** PYINIT: initialization completed **********************
> Initialisation done
>
> Generating all events
> stdhep initialized
>
>
> The file "tag_1_pythia.log" remains equal after to stop pythia.
> Now I done the other run through the proc file and this happened:
>
> daniel@daniel:~/Documents/MadGraph5$ cd PROC_Higgs_Effective_Couplings_UFO_3/
> daniel@daniel:~/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3$ ./bin/madevent
> launch in debug mode
> ************************************************************
> * *
> * W E L C O M E to M A D G R A P H 5 *
> * M A D E V E N T *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 5.1.5.5 *
> * *
> * The MadGraph Development Team - Please visit us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * *
> * Type 'help' for in-line help. *
> * *
> ************************************************************
> load configuration from /home/daniel/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3/Cards/me5_configuration.txt
> load configuration from /home/daniel/Documents/MadGraph5/input/mg5_configuration.txt
> load configuration from /home/daniel/Documents/MadGraph5/PROC_Higgs_Effective_Couplings_UFO_3/Cards/me5_configuration.txt
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> Using default eps viewer "gv". Set another one in ./input/mg5_configuration.txt
> Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> MGME5>pythia
> Which programs do you want to run?
> 0 / auto : running existing card
> 1 / pythia : Pythia
> 2 / pgs : Pythia + PGS
> [0, 1, 2, auto, pythia, pgs][60s to answer]
>> 1
> Will run in mode pythia
> Do you want to edit one cards (press enter to bypass editing)?
> 1 / pythia : pythia_card.dat
> 9 / plot : plot_card.dat
> Path to a valid card.
> [0, done, 1, pythia, 9, plot, enter path][60s to answer]
>> 0
> compile directory
> Running Pythia
>
>
> and remains so for a long time... after to stop the run show me the following;
>
> ^CTraceback (most recent call last):
> File "./bin/madevent", line 83, in <module>
> subprocess.call([sys.executable] + ['-O'] + sys.argv)
> File "/usr/lib/python2.7/subprocess.py", line 493, in call
> Stop by the user
> stopping all operation
> in order to quit madevent please enter exit
> pythia Stop by the user
> MGME5> return Popen(*popenargs, **kwargs).wait()
> File "/usr/lib/python2.7/subprocess.py", line 1291, in wait
> pid, sts = _eintr_retry_call(os.waitpid, self.pid, 0)
> File "/usr/lib/python2.7/subprocess.py", line 478, in _eintr_retry_call
> return func(*args)
> KeyboardInterrupt
>
> thanks Paolo for the patience with this problem!
>
> Cheers
>
> daniel
>
> --
> You received this question notification because you are a direct
> subscriber of the question.

Hi Paolo,

With the luck that I did not need to change the pdg codes or something like that. My colleague Alexandre Alves found the solution in a straightforward way, that is the pythia package need to have the final states of the process SM-like and the process p p > U1 U1~ has a no SM-like final state. If I do decay the U1 particle in SM particles the pythia run without any ṕroblem. Any way I believe that I was not clear with the process and sorry for that!! ....thank you so much Paolo for the help .

daniel

Hi Daniel,
good that you solved. I did not know that particular
feature in Pythia, it is good to know anyhow.
Thanks.
Cheers,
Paolo

On 16 Jun 2014, at 22:16, daniel camargo <email address hidden> wrote:

> Question #250051 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/250051
>
> Status: Answered => Solved
>
> daniel camargo confirmed that the question is solved:
> Hi Paolo,
>
> With the luck that I did not need to change the pdg codes or something
> like that. My colleague Alexandre Alves found the solution in a
> straightforward way, that is the pythia package need to have the final
> states of the process SM-like and the process p p > U1 U1~ has a no SM-
> like final state. If I do decay the U1 particle in SM particles the
> pythia run without any ṕroblem. Any way I believe that I was not clear
> with the process and sorry for that!! ....thank you so much Paolo for
> the help .
>
> daniel
>
> --
> You received this question notification because you are a direct
> subscriber of the question.