wrong number of unweighted events

Hi Alexandra,

That's weird.
could you provide to me (omatt __AT__ illinois.edu) yours Cards for the RS model?
(I need the proc_card, param_card and run_card)

Thanks,

Olivier

Hi again,
to the problem above, I don't think this should be the source of error, but anyway
to generate events to a large sample of parameters I'm using simple macros like that ones

http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/acarvalh/test/Macros/vary_inMRLR.sh?revision=1.1&view=markup
(generates samples with rigth number of events)
http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/acarvalh/test/Macros/vary_masses_Gr.sh?revision=1.1&view=markup
(generate samples with wrong number of events)

this error depends on the model, but not of the process chosen...
thanks again for help
Alexandra

Hi Alexandra,

Thanks for the model and the cards.
I have no problem with the generation and the number of events.
So Two questions:
1) Which version of MG5 did you use? (If not 1.5.2 could you try with this one?)
2) How do you get the number of events (did you get this number of the html page, from the line
    Nb of events : 20000 in the terminal or by counting the nb of events in a text file/ root file and if so which file?

Cheers,

Olivier

On Oct 23, 2012, at 3:36 PM, alexandra <email address hidden> wrote:

> Question #212130 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/212130
>
> Status: Needs information => Open
>
> alexandra gave more information on the question:
> Hi again,
> to the problem above, I don't think this should be the source of error, but anyway
> to generate events to a large sample of parameters I'm using simple macros like that ones
>
> http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/acarvalh/test/Macros/vary_inMRLR.sh?revision=1.1&view=markup
> (generates samples with rigth number of events)
> http://cmssw.cvs.cern.ch/cgi-bin/cmssw.cgi/UserCode/acarvalh/test/Macros/vary_masses_Gr.sh?revision=1.1&view=markup
> (generate samples with wrong number of events)
>
> this error depends on the model, but not of the process chosen...
> thanks again for help
> Alexandra
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Hi Olivier

I use MadGraph5_v1_4_5, I'll try the new one now

I get the number of events by counting with grep the .lhe (like on the examples below)
(indeed I only discovered it by counting with a fastjet code, after pass by pythia)
indeed on the LHE header I can see something like

  20000 = nevents ! Number of unweighted events requested
      71 = iseed ! rnd seed (0=assigned automatically=default))

(the iseed was automatic (0) before the generation)
but in below of the .lhe file I see I have

</MGRunCard>
<MGGenerationInfo>
# Number of Events : 79967
# Integrated weight (pb) : .36606E-06
# Max wgt : .25658E-10
# Average wgt : .11980E-10

##################################################
examples

> grep -o "<event>" MGraviton_3000.lhe | wc -w
79967

>grep -o "<event>" MGraviton_1000.lhe | wc -w
40527

 > grep -o "<event>" MGraviton_2000.lhe | wc -w
45262

> grep -o "<event>" MGraviton_500.lhe | wc -w
41914

but even if the new version works, there would be no apparent reason?
cheers
Alexandra

Hi Alexandra,

In the analysis script that you point to me. you do:
gunzip Events/M_$i/events.lhe.gz
mv Events/M_$i/events.lhe Graviton_Parton_on/MGraviton_$i.lhe
instead of
gunzip Events/M_$i/unweighted_events.lhe.gz
mv Events/M_$i/events.lhe Graviton_Parton_on/MGraviton_$i.lhe

(while in the other script you take the correct file)
You always have a larger number of weighted events compare to the physical unweighted events.

Cheers,

Olivier

On Oct 23, 2012, at 7:05 PM, alexandra <email address hidden> wrote:

> Question #212130 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/212130
>
> Status: Answered => Open
>
> alexandra is still having a problem:
> Hi Olivier
>
> I use MadGraph5_v1_4_5, I'll try the new one now
>
> I get the number of events by counting with grep the .lhe (like on the examples below)
> (indeed I only discovered it by counting with a fastjet code, after pass by pythia)
> indeed on the LHE header I can see something like
>
> 20000 = nevents ! Number of unweighted events requested
> 71 = iseed ! rnd seed (0=assigned automatically=default))
>
> (the iseed was automatic (0) before the generation)
> but in below of the .lhe file I see I have
>
> </MGRunCard>
> <MGGenerationInfo>
> # Number of Events : 79967
> # Integrated weight (pb) : .36606E-06
> # Max wgt : .25658E-10
> # Average wgt : .11980E-10
>
> ##################################################
> examples
>
>> grep -o "<event>" MGraviton_3000.lhe | wc -w
> 79967
>
>> grep -o "<event>" MGraviton_1000.lhe | wc -w
> 40527
>
>> grep -o "<event>" MGraviton_2000.lhe | wc -w
> 45262
>
>> grep -o "<event>" MGraviton_500.lhe | wc -w
> 41914
>
> but even if the new version works, there would be no apparent reason?
> cheers
> Alexandra
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Hi Olivier
thanks much!
so sorry for bother you about a so silly typo :-/
Alexandra