Generation of many events

Do 4 generation of 50k
(you have the multi_run option to launch them in a easy way)

./bin/madevent
then type
multi_run 4

Cheers,

Olivier
> On Jan 26, 2016, at 23:07, sgurruti <email address hidden> wrote:
>
> New question #283732 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/283732
>
> Hello,
> I want to generate 200k events in Madgraph, and run them on Pythia but the file "run_card.dat" says:
>
> << If you want to run Pythia, avoid more than 50k events in a run >>
>
> My question is: How can I generate my required number of events?
>
> Best regards,
> Sebastian
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

I am using ./bin/mg5_aMC GetCrossSection.dat to run the following card:

  set group_subprocesses Auto
  set ignore_six_quark_processes False
  set gauge unitary
  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 u~ c~ d~ s~
  define l+ = e+ mu+
  define l- = e- mu-
  define vl = ve vm vt
  define vl~ = ve~ vm~ vt~
  import model InertDoublet_UFO
  generate p p > w+ h2 h2, (w+ > j j)
  add process p p > w- h2 h2, (w- > j j)
  output CrossSection_3_100
  launch
  1
  2
  3
  0
  set lam2 3.743240
  set lamL -0.494199
  set mmH0 407.606971
  set mmA0 457.952477
  set mmHch 471.009235
  set Wh 0.004070
  set mmh 125
  set WH0 0
  set WA0 Auto
  set WHch Auto
  set ebeam1 7000
  set ebeam2 7000
  0
  launch CrossSection_3_100 -i
  print_results --path=./CrossSection_3_100/cs.txt --format=short

I don't know how do in this context... can you help me?

Best regards,
Sebastian

Hi,

./bin/madevent is the equivalent of launch -i, so you can do the following:

 set group_subprocesses Auto
  set ignore_six_quark_processes False
  set gauge unitary
  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 u~ c~ d~ s~
  define l+ = e+ mu+
  define l- = e- mu-
  define vl = ve vm vt
  define vl~ = ve~ vm~ vt~
  import model InertDoublet_UFO
  generate p p > w+ h2 h2, (w+ > j j)
  add process p p > w- h2 h2, (w- > j j)
  output CrossSection_3_100
  launch CrossSection_3_100 -i
  multi_run X
 pgs=ON
 delphes=ON
  set lam2 3.743240
  set lamL -0.494199
  set mmH0 407.606971
  set mmA0 457.952477
  set mmHch 471.009235
  set Wh 0.004070
  set mmh 125
  set WH0 0
  set WA0 Auto
  set WHch Auto
  set ebeam1 7000
  set ebeam2 7000
  print_results --path=./CrossSection_3_100/cs.txt --format=short

or you can put in a given file:

  multi_run X
 pgs=ON
 delphes=ON
  set lam2 3.743240
  set lamL -0.494199
  set mmH0 407.606971
  set mmA0 457.952477
  set mmHch 471.009235
  set Wh 0.004070
  set mmh 125
  set WH0 0
  set WA0 Auto
  set WHch Auto
  set ebeam1 7000
  set ebeam2 7000
  print_results --path=./CrossSection_3_100/cs.txt --format=short

and run it like ./bin/madevent PATH.

One comment here. You should never use the number "1","2","3" when scripting since their number will not be kept retro-compatible in future version. It is better to use the associate flag.

Cheers,

Olivier

Thank you so much!

Hi Olivier,
                         I have a similar type question. Suppose I have generated the events using multi_run option with pythia=ON. It gives .hep files for each run. But how to combine all the .hep files ?

Regards,
Shiba

You do not combine them, you just use all of them (load all of them) in your analysis framework

Cheers,

Olivier
> On Jan 28, 2016, at 05:37, shibasipu <email address hidden> wrote:
>
> Question #283732 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/283732
>
> shibasipu posted a new comment:
> Hi Olivier,
> I have a similar type question. Suppose I have generated the events using multi_run option with pythia=ON. It gives .hep files for each run. But how to combine all the .hep files ?
>
>
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.