generate_events taking forever/pythia.hep problem

Hi,

Looks like your cross-section is infinite:
> 25929438.0451 +- 10758000000.4

this sounds to be the case, since you have very large cross-section with very large error. Which is typical when the function that you try to integrate is not bounded.
In such case, the code is indeed very slow and ends up with failing to converge/produce events.

Are you sure that all the width of s-channel particle are set to a physical value?
Did you put the cut required for your process?

Cheers,

Olivier

> On Sep 1, 2016, at 06:22, Hemang <email address hidden> wrote:
>
> New question #376598 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/376598
>
> I am a bit of a beginner with madgraph and I was trying to get data on an event.
>
> MG5_aMC>import model WpZp_MonoXp
> INFO: Change particles name to pass to MG5 convention
> Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
> 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+ zp wp+ h hzero azero h+ w- wp- h- ta- n1 nd ta+ n1~ nd~
> MG5_aMC>generate wp+ zp > t b~ n1 n1~
> INFO: Checking for minimal orders which gives processes.
> INFO: Please specify coupling orders to bypass this step.
> INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS OV1e+OV1tau+2*QED+QCD+OgR+OV1mu
> INFO: Trying process: wp+ zp > t b~ n1 n1~ WEIGHTED<=4 @1
> INFO: Process has 10 diagrams
> 1 processes with 10 diagrams generated in 0.205 s
> Total: 1 processes with 10 diagrams
> MG5_aMC>output proj
> INFO: directory /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj already exists.
> If you continue this directory will be deleted and replaced.
> Do you want to continue? [y, n][60s to answer]
>> y
> INFO: initialize a new directory: proj
> INFO: remove old information in proj
> INFO: Organizing processes into subprocess groups
> INFO: Generating Helas calls for process: wp+ zp > t b~ n1 n1~ WEIGHTED<=4 @1
> INFO: Processing color information for process: wp+ zp > t b~ n1 n1~ @1
> INFO: Creating files in directory P1_wppzp_tbxn1n1x
> INFO: Generating Feynman diagrams for Process: wp+ zp > t b~ n1 n1~ WEIGHTED<=4 @1
> INFO: Finding symmetric diagrams for subprocess group wppzp_tbxn1n1x
> History written to /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj/Cards/proc_card_mg5.dat
> Generated helas calls for 1 subprocesses (10 diagrams) in 0.047 s
> Wrote files for 29 helas calls in 0.086 s
> Export UFO model to MG4 format
> ALOHA: aloha creates FFV3 routines
> ALOHA: aloha creates FFV4 routines
> ALOHA: aloha creates FFV5 routines
> save configuration file to /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj/Cards/me5_configuration.txt
> INFO: Use Fortran compiler gfortran
> INFO: Use c++ compiler clang
> INFO: Generate jpeg diagrams
> nice: gs: No such file or directory
> nice: gs: No such file or directory
> cp: card.jpg: No such file or directory
> INFO: Generate web pages
> Output to directory /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj done.
> Type "launch" to generate events from this process, or see
> /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj/README
> Run "open index.html" to see more information about this process.
> MG5_aMC>quit
>
>
>
>
> Hemangs-MBP:mg5_amc_v2_4_3 hemang$ cd proj
> Hemangs-MBP:proj hemang$ ./bin/generate_events test
> launch in debug mode
> ************************************************************
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * M A D E V E N T *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 5.2.4.3 20xx-xx-xx *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * *
> * Type 'help' for in-line help. *
> * *
> ************************************************************
> INFO: load configuration from /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj/Cards/me5_configuration.txt
> INFO: load configuration from /Users/hemang/Downloads/MG5_aMC_v2_4_3/input/mg5_configuration.txt
> INFO: load configuration from /Users/hemang/Downloads/MG5_aMC_v2_4_3/proj/Cards/me5_configuration.txt
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> generate_events test
> The following switches determine which programs are run:
> 1 Run the pythia shower/hadronization: pythia=OFF
> 2 Run PGS as detector simulator: pgs=OFF
> 3 Run Delphes as detector simulator: delphes=OFF
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weights to the events based on changing model parameters: reweight=Not available (requires NumPy)
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>> 1
> The following switches determine which programs are run:
> 1 Run the pythia shower/hadronization: pythia=ON
> 2 Run PGS as detector simulator: pgs=OFF
> 3 Run Delphes as detector simulator: delphes=OFF
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weights to the events based on changing model parameters: reweight=Not available (requires NumPy)
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>> 2
> The following switches determine which programs are run:
> 1 Run the pythia shower/hadronization: pythia=ON
> 2 Run PGS as detector simulator: pgs=ON
> 3 Run Delphes as detector simulator: delphes=OFF
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weights to the events based on changing model parameters: reweight=Not available (requires NumPy)
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>> 3
> The following switches determine which programs are run:
> 1 Run the pythia shower/hadronization: pythia=ON
> 2 Run PGS as detector simulator: pgs=ON
> 3 Run Delphes as detector simulator: delphes=ON
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weights to the events based on changing model parameters: reweight=Not available (requires NumPy)
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>>
> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> 3 / pythia : pythia_card.dat
> 4 / pgs : pgs_card.dat
> 5 / delphes : delphes_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.
> - call an external program (ASperGE/MadWidth/...).
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, 3, pythia, 4, enter path, ... ][60s to answer]
>> 0
> Generating 10000 events with run name test
> survey test
> INFO: compile directory
> compile Source Directory
> Using random number seed offset = 21
> INFO: Running Survey
> Creating Jobs
> Working on SubProcesses
> INFO: P1_wppzp_tbxn1n1x
> INFO: Idle: 1, Running: 3, Completed: 0 [ current time: 21h43 ]
> INFO: Idle: 0, Running: 3, Completed: 1 [ 13.2s ]
> INFO: Idle: 0, Running: 2, Completed: 2 [ 20.3s ]
> INFO: Idle: 0, Running: 0, Completed: 4 [ 20.6s ]
> INFO: End survey
> refine 10000
> Creating Jobs
> INFO: Refine results to 10000
> INFO: Generating 10000.0 unweigthed events.
> INFO: Effective Luminosity 0.000462794449272 pb^-1
> INFO: need to improve 3 channels
> Current estimate of cross-section: 25929438.0451 +- 10758000000.4
> P1_wppzp_tbxn1n1x
> INFO: Idle: 6, Running: 8, Completed: 0 [ current time: 21h44 ]
> INFO: Idle: 5, Running: 8, Completed: 1 [ 6.1s ]
> INFO: Idle: 4, Running: 8, Completed: 2 [ 22m 28s ]
>
>
>
>
>
>
>
>
> The last time I ran this I ended up waiting 7 + hours for it to finish and it ended up having problems with this
> No file file pythia_events.hep currently available
>
> I'm not sure how to proceed because I don't think this should take 7 hours and I'm not sure how to fix the pythia_events issue either.
> Thanks
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hello,
so I may be a bit in over my head but I'm not sure what function I'm trying to integrate? Is that something that is included when i generate a process and output it ? I'm not really sure what I should edit and if you could explain it to me that would really help a lot.

Hi,

> so I may be a bit in over my head but I’m not sure what function I'm trying to integrate?

you are trying to integrate the matrix element of the process that you specified.

> Is that something that is included when i generate a process and output it ?

Those command correspond to create a matrix-element (i.e. the function that you need to integrate to get the cross-section and based on that function we generate the events according to the associate density of probability)

> I’m not really sure what I should edit and if you could explain it to me that would really help a lot.

The first think you need to do is to understand why your matrix-element can be infinite and for which condition.
If you need help on that you should ask some of your colleague who know the details of the model that you are using.

With that knowledge, you then need to either fix the problem in the param_card (if this is related to a mass/width problem).
If you need to implement a cut to forbid part of the phase-space, you should either include it in the run_card.dat (for simple cut) or implementing yourself
in SubProcesses/cuts.f

Cheers,

Olivier

> On Sep 1, 2016, at 08:57, Hemang <email address hidden> wrote:
>
> Question #376598 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/376598
>
> Status: Answered => Open
>
> Hemang is still having a problem:
> Hello,
> so I may be a bit in over my head but I'm not sure what function I'm trying to integrate? Is that something that is included when i generate a process and output it ? I'm not really sure what I should edit and if you could explain it to me that would really help a lot.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.