"set nevents 50000" with error:

Dear Yingu,

without the second script (and the command that you use to run it) I would not be able to help you.

Cheers,

Olivier

set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
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~
import model ZpHiggs_UFO
generate p p > n1 n1 Hs, (Hs > b b~)
add process p p > n1 n1 Hs j, (Hs > b b~)

output ms50_m300
launch
pythia=OFF
set nevents 50000
set DECAY 9000006 Auto # WS
set DECAY 9000007 Auto # wZp
set misset 400.0
set GQ 0.25
set GX 0.23
set MHS 50
set MZP 300
set MDM 100

Hi,

The problem is with the command pythia=OFF
In order for that command to be understood by the program, you have to have pythia-pgs installed.
Without pythia-pgs this command is invalid.

You can replace this command by shower=OFF to fix your problem.
(an alternative is obviously to install pythia-pgs)

Cheers,

Olivier

> On Jan 11, 2017, at 11:33, yingpu xiahou <email address hidden> wrote:
>
> Question #435186 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/435186
>
> yingpu xiahou posted a new comment:
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model sm
> 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~
> import model ZpHiggs_UFO
> generate p p > n1 n1 Hs, (Hs > b b~)
> add process p p > n1 n1 Hs j, (Hs > b b~)
>
> output ms50_m300
> launch
> pythia=OFF
> set nevents 50000
> set DECAY 9000006 Auto # WS
> set DECAY 9000007 Auto # wZp
> set misset 400.0
> set GQ 0.25
> set GX 0.23
> set MHS 50
> set MZP 300
> set MDM 100
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

But getting an error from 'set nevents' is a sign of pythia being misconfigured.
Either something is wrong with the setup or with the input file.I don't think my setup is wrong, so I change the .input file a little bit.Just as you recommend, change "pythia=OFF" into "shower=OFF".

E.g. if you have "pythia=OFF" my understanding is that "set nevents" will not work.
Instead "shower=OFF" should be used.

After I change "pythia=OFF" into "shower=OFF", I tried another .input file,and new error comes:

 INFO: generating the production square matrix element
INFO: generate p p > Zp S --no_warning=duplicate;
INFO: Done 2.429
Command "generate_events run_01" interrupted with error:
Exception : No decay define for process.
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/storage/yx16344/MG5_aMC_v2_5_2/ZpHiggs/run_01_tag_1_debug.log'.
Please attach this file to your report.

Here is my script:

set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s b t u~ c~ d~ s~ b~ t~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model ZpHiggs_UFO --modelname
generate p p > Zp S
#add process p p > chi chi~ j j j

output ZpHiggs
launch
shower=OFF
madspin=ON
#decay Zp > DM DM
#decay S > b b~
set nevents 20000
#set misset 60
#set asrwgtflavor 5
#set Mchi 10
#set Lambda 100

Btw, how can I attach a file for you , I want to attach run_01_tag_1_debug.log.

Thanks,Cheers

Yingpu

Hi,

Could you please tell me what is the difference between 'install pythia-pgs' and 'shower=OFF' in order to solve the 'set nevents' error?

I tried both and both can solve my problem.

But, I am wondering what is the difference between them.Which method is better?

Thanks,Cheers

Yingpu

Hi,

You really solve my problem!

Thank you so much!

Cheers,

Yingpu

Thanks Olivier Mattelaer, that solved my question.