No systematic variations calculated, weights missing in hepmc file

Asked by Früd Braren

Hi,

I've got the problem that I don't get systematic uncertainties and the corresponding weights in the hepmc file if I use [QCD] at the end of 'generate', i.e. loop corrections. I'm using Pythia8 to shower the events, and the Madgraph version is 5.2.6.5.

As an example: If I do 'p p > t t~', I get systematic uncertainties and all the corresponding weights in the .hepmc file. The 'sm'-model is used.

If I do 'p p > t t~ [QCD]', I don't get systematic variations, and the weights are present in the .lhe file, but not in the .hepmc file. The 'loop-sm'-model is used automatically.
A further difference is that the run_card.dat looks quite different in the two cases (longer / more options in the first case).
Some pertinent details from the run_card in the [QCD] case:
lhapdf = pdlabel
 91700 = lhaid
 True = fixed_ren_scale
 True = fixed_fac_scale
 91.118 = muR_ref_fixed
 91.118 = muF_ref_fixed
 -1 = dynamical_scale_choice
 1.0 = muR_over_ref
 1.0 = muF_over_ref
 1.0, 2.0, 0.5 = rw_rscale
 1.0, 2.0, 0.5 = rw_fscale
 True = reweight_scale
 True = reweight_PDF
 True = store_rwgt_info

The summary I get in the first case (no [QCD]) looks like this:
# original cross-section: 611.966
# scale variation: +31.6% -22.5%
# central scheme variation: + 0% -45.1%
# PDF variation: +1.88% -1.88%
#
# dynamical scheme # 1 : 437.959 +28.3% -20.7% # \sum ET
# dynamical scheme # 2 : 356.783 +26.1% -19.5% # \sum\sqrt{m^2+pt^2}
# dynamical scheme # 3 : 449.901 +28.3% -20.7% # 0.5 \sum\sqrt{m^2+pt^2}
# dynamical scheme # 4 : 335.91 +25.4% -19.1% # \sqrt{\hat s}

The summary I get in the second case looks like this:
   --------------------------------------------------------------
      Summary:
      Process p p > t t~ [QCD]
      Run at p-p collider (6500.0 + 6500.0 GeV)
      Number of events generated: 100
      Total cross section: 8.034e+02 +- 5.8e+00 pb
   --------------------------------------------------------------

Any help with getting the systematics variations right in the loop-corrections case would be much appreciated.

Cheers,
Früd

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Paolo Torrielli Edit question
Last query:
Last reply:
Revision history for this message
Paolo Torrielli (paolotorriell) said :
#1

Dear Früd,

you can also see the discussion at https://answers.launchpad.net/mg5amcnlo/+question/678584,
which is related to the problem you're having.

The short answer is that currently, if you use the shower interface
that is distributed with the code, the weights are simply not passed.
We are working with some Pythia authors in order to introduce this
feature, which will become available in one of the future releases.

Best regards.
Paolo

Can you help with this problem?

Provide an answer of your own, or ask Früd Braren for more information if necessary.

To post a message you must log in.