systematics for NLO process

The only information reported is the summary for the cross-section.
the additionaly weight are also present whithin the lhe file such that you can plot the diffferential cross-section with the error.

Cheers,

Olivier

how to plot differential cross-section with errors?
or is their any why to calculate systematics for NLO like LO?
as below
I use LO to compute cross-section which automatically give systematics for lhapdf error sets
run_01
p p
6500.0 x 6500.0 GeV
 tag_1
8.692 ± 0.027 ± systematics

Yes you can run the exact same program as for LO to have the systematics computation at NLO.
See here:
https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Systematics#Systematicspythonmodule1

Cheers,

Olivier

so it means i have to write below to lines in run card of NLO

  systematics = systematics_program

 ['--pdf=central', '--dyn=-1,3', '--together=mur,muf'] = systematics_arguments

right??

And have to setup one additional parameter to include in the chef file the information needed for the systematics program to run.

Cheers,

Olivier

PS: I should check if this is still working in 3.3.0 since they have slightly change the format. So if you have issue please use our LTS version (2.9.6)

> On 23 Nov 2021, at 18:15, Anees khan <email address hidden> wrote:
>
> Question #699600 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/699600
>
> Status: Answered => Open
>
> Anees khan is still having a problem:
> so it means i have to write below to lines in run card of NLO
>
> systematics = systematics_program
>
> ['--pdf=central', '--dyn=-1,3', '--together=mur,muf'] =
> systematics_arguments
>
>
> right??
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

where is this chef file?

Sorry chef was a typo.
This should have been "run_card"

You also have to set the parameter
True = store_reweight_info

I set True = store_reweight_info
but i can't see systimatics for NLO calculations?

Yes that's normal those are hidden parameter.

Cheers,

Olivier

> On 24 Nov 2021, at 08:55, Anees khan <email address hidden> wrote:
>
> Question #699600 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/699600
>
> Status: Answered => Open
>
> Anees khan is still having a problem:
> I set True = store_reweight_info
> but i can't see systimatics for NLO calculations?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Then how can i see those hidden parameters?

can any one help me?

Sure, i have copy paste here the instruction from our webpage:

   In order to use this module you have to set the run_card parameter "store_reweight_info" to True (i.e. having the following line in the run_card)

> True = store_reweight_info
This is not the default value of this parameter. If you want that this module is directly called after the event generation you have to add the following lines in the run_card:

  systematics = systematics_program
parameter of the systematics module can be specified in systematics_argument for example with the following line in the run_card:

> ['--pdf=central', '--dyn=-1,3', '--together=mur,muf'] = systematics_arguments

The list of available options are the same as the one at LO (see above)

As for the LO case, a dedicated function exists in the aMC@NLO interface (available via either launch -i or via the ./bin/aMC@NLO script)

> On 29 Nov 2021, at 07:35, Anees khan <email address hidden> wrote:
>
> Question #699600 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/699600
>
> Anees khan gave more information on the question:
> can any one help me?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

and please note the line
"you have to **add** the following lines in the run_card:"

Cheers,

Olivier

> On 29 Nov 2021, at 07:42, Olivier Mattelaer <email address hidden> wrote:
>
> Sure, i have copy paste here the instruction from our webpage:
>
>
>
> In order to use this module you have to set the run_card parameter "store_reweight_info" to True (i.e. having the following line in the run_card)
>
>> True = store_reweight_info
> This is not the default value of this parameter. If you want that this module is directly called after the event generation you have to add the following lines in the run_card:
>
> systematics = systematics_program
> parameter of the systematics module can be specified in systematics_argument for example with the following line in the run_card:
>
>> ['--pdf=central', '--dyn=-1,3', '--together=mur,muf'] = systematics_arguments
>
> The list of available options are the same as the one at LO (see above)
>
> As for the LO case, a dedicated function exists in the aMC@NLO interface (available via either launch -i or via the ./bin/aMC@NLO script)
>
>
>> On 29 Nov 2021, at 07:35, Anees khan <<email address hidden> <mailto:<email address hidden>>> wrote:
>>
>> Question #699600 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/699600 <https://answers.launchpad.net/mg5amcnlo/+question/699600>
>>
>> Anees khan gave more information on the question:
>> can any one help me?
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>

an error came

INFO: Events generated
systematics run_02 --mur=0.5,1,2 --muf=0.5,1,2 --pdf=errorset
INFO: Running Systematics computation
INFO: Idle: 4, Running: 0, Completed: 0 [ current time: 17h25 ]
WARNING: program /usr/bin/python3 -O /home/anees/MG5_aMC_v3.3.0/MG5_aMC_v3_3_0/anees/bin/internal/systematics.py events.lhe.gz ./tmp_0_events.lhe.gz --mur=0.5,1,2 --muf=0.5,1,2 --pdf=errorset --start_event=0 --stop_event=2500 --result=./log_sys_0.txt --lhapdf_config=/home/anees/MG5_aMC_v3.3.0/MG5_aMC_v3_3_0/HEPTools/lhapdf6_py3/bin/lhapdf-config launch ends with non zero status: 1. Stop all computation
INFO: Running Systematics computation
Command "launch auto " interrupted with error:
KeyError : 'scalefact'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/anees/MG5_aMC_v3.3.0/MG5_aMC_v3_3_0/anees/run_02_tag_1_debug.log'.
Please attach this file to your report.
INFO:
quit
INFO:

Hi,

As i was expecting, this mode is not working in our 3.3.0 release, you should use our long term stable release (2.9.7) for having this mode working. (3.3.1 should have that fixed)

Cheers,

Olivier

yes now systematics are included but with no data in it
https://ibb.co/qpTpv96
how can i solve this now?

and i also need systematics after decay (in madspin card)

yes now systematics are included but with no data in it
https://ibb.co/qpTpv96
how can i solve this now?

The above problem is solved now,but showing systematics before decay.I also need systematics after decay
How can i do this?

You can just multiply by the Branching ratio (if unique) to have the result after madspin.

Cheers,

Olivier

> On 30 Nov 2021, at 07:30, Anees khan <email address hidden> wrote:
>
> Question #699600 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/699600
>
> Anees khan gave more information on the question:
> yes now systematics are included but with no data in it
> https://ibb.co/qpTpv96
> how can i solve this now?
>
> The above problem is solved now,but showing systematics before decay.I also need systematics after decay
> How can i do this?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Thanks Olivier Mattelaer, that solved my question.

Thanks For your help

https://ibb.co/7tw69Rk
when i try to calculate systematics for HERAPDF20_NLO_EIG and ABMP16_4_nlo there is no data (for errorset)stored in parton_systematic.log??why this heppen?

yes me too having same issue with these pdfs ,but showing errorset calculations of NNPDF30_nlo errorsets

any solution?

Did you check if those pdf have an assocciated method to compute the error?

Cheers,

Olivier

HERAPDF20_NLO_EIG ErrorType: hessian
ABMP16_5_nlo ErrorType: symmhessian

Hum weird then, I'm looking into it now.

With this series of command:
generate p p > t t~ [QCD]
output
launch
shower=OFF
set pdlabel=lhapdf
set lhaid 60600
set store_reweight_info T
set systematics_program systematics

I do not reproduce the error.
Actually I did face the error once but never succeed to repeat it.
Maybe an issue if you do not have the pdf set on your system?

With this series of command:
import model sm-no_b_mass
generate p p > z b [QCD]
output
launch
shower=OFF
IN RUN CARD
set pdlabel=lhapdf
set lhaid =61100
True = store_rwgt_info ! Store info for reweighting in LHE file
systematics = systematics_program
  ['--pdf=errorset', '--dyn=-1,3', '--together=mur,muf'] = systematics_arguments

still no data in parton_systematic.log

But here you mix four and five flavor model.
(your model is five flavor but you ask for b jet flavor which only makes sense in four flavour)
So the meaning of such computation at NLO is not clear to me and even more how to do the running (should it be in four or five flavor?)

So if the issue is only when you mix four and five flavor, I would not care about it.

Cheers,

Olivier