pythia fails showering photon initiated processes

Can you try with the 2.4.0 beta version?
(available here:https://launchpad.net/mg5amcnlo)
It should normally fix your problem,

Cheers,

Olivier

Hi Olivier,
Thank you very much for your fast reply.
It fixed my problem. now the standard model process is generated and goes through pythia successfully.

but, it seems that in beta version (what is beta version for?) new physics models can not be included. I also test with default models in models directory of madgraph.
when I want to import model it gives errors (in the following).
how can I test new physics model which contribute to photon initiated process.
Thanks a lot,
Reza

************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.4.0.beta 2016-04-06 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* and *
* http://amcatnlo.web.cern.ch/amcatnlo/ *
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
************************************************************
load MG5 configuration from input/mg5_configuration.txt
fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
 MG5_aMC> set fastjet /PATH/TO/fastjet-config

set lhapdf to lhapdf-config
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
Loading default model: sm
INFO: Restrict model sm with file models/sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Defined multiparticle p = g u c d s u~ c~ d~ s~
Defined multiparticle j = g u c d s u~ c~ d~ s~
Defined multiparticle l+ = e+ mu+
Defined multiparticle l- = e- mu-
Defined multiparticle vl = ve vm vt
Defined multiparticle vl~ = ve~ vm~ vt~
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
MG5_aMC>import model 2HDM
Command "import model 2HDM" interrupted with error:
AttributeError : 'module' object has no attribute 'load_model'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

Hi,

version 2.4.0.beta is perfectly able to import BSM model.
The reason why it crash for you is likely because you have wrongly untar your BSM model in the wrong directory.
and at the same time overwrite one or more of the python file of mg5.
This is typically the case for the __init__.py file,

You can either re-download/re-unpack the full package and/or just copy the __init__.py from another version of MG5 by doing something like
cp 2.3.3/models/__init__.py 2.4.0/models/__init__.py

Cheers,

Olivier

> On Apr 6, 2016, at 19:37, reza <email address hidden> wrote:
>
> Question #290074 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290074
>
> reza posted a new comment:
> Hi Olivier,
> Thank you very much for your fast reply.
> It fixed my problem. now the standard model process is generated and goes through pythia successfully.
>
> but, it seems that in beta version (what is beta version for?) new physics models can not be included. I also test with default models in models directory of madgraph.
> when I want to import model it gives errors (in the following).
> how can I test new physics model which contribute to photon initiated process.
> Thanks a lot,
> Reza
>
> ************************************************************
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 2.4.0.beta 2016-04-06 *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * and *
> * http://amcatnlo.web.cern.ch/amcatnlo/ *
> * *
> * Type 'help' for in-line help. *
> * Type 'tutorial' to learn how MG5 works *
> * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
> * Type 'tutorial MadLoop' to learn how MadLoop works *
> * *
> ************************************************************
> load MG5 configuration from input/mg5_configuration.txt
> fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
> Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
> MG5_aMC> set fastjet /PATH/TO/fastjet-config
>
> set lhapdf to lhapdf-config
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> Loading default model: sm
> INFO: Restrict model sm with file models/sm/restrict_default.dat .
> INFO: Run "set stdout_level DEBUG" before import for more information.
> INFO: Change particles name to pass to MG5 convention
> Defined multiparticle p = g u c d s u~ c~ d~ s~
> Defined multiparticle j = g u c d s u~ c~ d~ s~
> Defined multiparticle l+ = e+ mu+
> Defined multiparticle l- = e- mu-
> Defined multiparticle vl = ve vm vt
> Defined multiparticle vl~ = ve~ vm~ vt~
> Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h w- ta- ta+
> MG5_aMC>import model 2HDM
> Command "import model 2HDM" interrupted with error:
> AttributeError : 'module' object has no attribute 'load_model'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

When I re-unpacked the beta version BSM models are imported.

But I have another problem, pythia run with no error and every thing seems OK but the out put file is empty and (Sorry I did not look at the pythia output LHE file before).
I am confused!
I test the standard model [ a a > l+ l-] process.
at the end of tag_1_pythia.log file I see the following error

Thanks,
Reza

Starting event loop

     Error type 2 has occured after 1 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 1 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 2 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 2 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 3 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 3 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 4 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 4 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 5 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system

     Error type 2 has occured after 5 PYEXEC calls:
     (PYSTRF:) failed to reconstruct jet system
  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 = 21111 *************
 ********* Total number of errors, including junctions = 21111 *************
 ********* Total number of warnings = 0 *************
 ********* Fraction of events that fail fragmentation cuts = 1.00000 *********

 Cross section (pb): 0.0000000000000000

Hi,

This is a question for the pythia author,
I guess that this is related to the fact that you do not have any QCD charged particles.
So the QCD shower is clearly irrelevant in this case. I’m not actually even sure that PY6 can handle such type of processes.

Cheers,

Olivier

> On Apr 6, 2016, at 22:22, reza <email address hidden> wrote:
>
> Question #290074 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/290074
>
> reza posted a new comment:
> Hi Olivier,
>
> When I re-unpacked the beta version BSM models are imported.
>
> But I have another problem, pythia run with no error and every thing seems OK but the out put file is empty and (Sorry I did not look at the pythia output LHE file before).
> I am confused!
> I test the standard model [ a a > l+ l-] process.
> at the end of tag_1_pythia.log file I see the following error
>
> Thanks,
> Reza
>
> Starting event loop
>
> Error type 2 has occured after 1 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 1 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 2 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 2 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 3 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 3 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 4 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 4 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 5 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
>
> Error type 2 has occured after 5 PYEXEC calls:
> (PYSTRF:) failed to reconstruct jet system
> 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 = 21111 *************
> ********* Total number of errors, including junctions = 21111 *************
> ********* Total number of warnings = 0 *************
> ********* Fraction of events that fail fragmentation cuts = 1.00000 *********
>
> Cross section (pb): 0.0000000000000000
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Dear experts of MG5

I am studying proton - photon collision going to two jets: p a > j j, with MG5_3_4_1. I get events at parton level (an lhe file), but the HEPMC file after Pythia8 is empty. No error messages are printed in the terminal screen.

Is it possible to pass events form this processes through Pythia8?

Thanks,
Leandro