Turn off PDF systematic weights in LO mode

Asked by Josh McFayden

Hi all,

I want to run SysCalc in leading order mode but only generate the scale variations and not the PDF (we will do this with an external tool later), is this possible?

I tried just setting sys_pdf to a dummy value, e.g.:
T = use_syst
0.5 1 2 = sys_scalefact
0.5 1 2 = sys_alpsfact
30 50 = sys_matchscale
F = sys_pdf

But that gives the error:
"SysCalc Failed. Please read the associate log to see the reason. Did you install the associate PDF set?

I tried with v2_3_2.

Cheers,

Josh.

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Josh McFayden
Solved:
Last query:
Last reply:
Revision history for this message
Alexis Kalogeropoulos (alkaloge) said :
#1

Hi Josh

You dont need to add the line

 F = sys_pdf

SysCalc will complains, as cannot find "F" as PDF set. To NOT run SysCalc for PDFs you simply have to remove this line

Regards

Alexis

Revision history for this message
Alexis Kalogeropoulos (alkaloge) said :
#2

Discussing with Olivier, currently if the "sys_pdf" is missing, MG will pass by default the CT10nlo.grid to SysCalc - So, I guess for now you can maybe live with that and in due time the behaviour will change if this line is missing not passing anything to SysCalc at all.

Regards

Alexis

Revision history for this message
Josh McFayden (mcfayden) said :
#3

Thanks Alexis.

Yes, it would be great if we could turn it off completely.

Could this change in behaviour go into the next release?

Best,

Josh.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#4

To be more clear,

Alexis claim that SysCalc can do that if run in standalone.
I will make a patch (on Friday) to create the correct syscalc output to support this situation.

Cheers,

Olivier

On 07 Oct 2015, at 15:27, Alexis Kalogeropoulos <email address hidden> wrote:

> Question #272136 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/272136
>
> Alexis Kalogeropoulos posted a new comment:
> Discussing with Olivier, currently if the "sys_pdf" is missing, MG will
> pass by default the CT10nlo.grid to SysCalc - So, I guess for now you
> can maybe live with that and in due time the behaviour will change if
> this line is missing not passing anything to SysCalc at all.
>
> Regards
>
> Alexis
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Josh McFayden (mcfayden) said :
#5

Hi Olivier,

Is there a patch for this yet?

Best,

Josh.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#6

Sorry,

doing that right now

Olivier
On 12 Oct 2015, at 11:08, Josh McFayden <email address hidden> wrote:

> Question #272136 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/272136
>
> Status: Answered => Open
>
> Josh McFayden is still having a problem:
> Hi Olivier,
>
> Is there a patch for this yet?
>
> Best,
>
> Josh.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#7

Here is the patch:

=== modified file 'madgraph/interface/madevent_interface.py'
--- old/madgraph/interface/madevent_interface.py 2015-10-11 19:18:37 +0000
+++ new/madgraph/interface/madevent_interface.py 2015-10-12 11:13:31 +0000
@@ -4160,6 +4160,12 @@
         card = pjoin(self.me_dir, 'bin','internal', 'syscalc_card.dat')
         template = open(pjoin(self.me_dir, 'bin','internal', 'syscalc_template.dat')).read()
         self.run_card['sys_pdf'] = self.run_card['sys_pdf'].split('#',1)[0].replace('&&',' \n ')
+
+ if self.run_card['sys_pdf'].lower() in ['', 'f', 'false', 'none', '.false.']:
+ self.run_card['sys_pdf'] = ''
+ if self.run_card['sys_alpsfact'].lower() in ['', 'f', 'false', 'none','.false.']:
+ self.run_card['sys_alpsfact'] = ''
+
         # check if the scalecorrelation parameter is define:
         if not 'sys_scalecorrelation' in self.run_card:
             self.run_card['sys_scalecorrelation'] = -1

=== modified file 'madgraph/various/banner.py'
--- old/madgraph/various/banner.py 2015-10-09 10:39:28 +0000
+++ new/madgraph/various/banner.py 2015-10-12 10:40:46 +0000
@@ -1627,7 +1627,7 @@
         self.add_param("xqcut", 0.0, cut=True)
         self.add_param("use_syst", True)
         self.add_param("sys_scalefact", "0.5 1 2", include=False)
- self.add_param("sys_alpsfact", "0.5 1 2", include=False)
+ self.add_param("sys_alpsfact", "None", include=False)
         self.add_param("sys_matchscale", "30 50", include=False)
         self.add_param("sys_pdf", "Ct10nlo.LHgrid", include=False)
         self.add_param("sys_scalecorrelation", -1, include=False)
@@ -1796,6 +1796,7 @@
                 self['use_syst'] = False
                 self['drjj'] = 0
                 self['drjl'] = 0
+ self['sys_alpsfact'] = "0.5 1 2"

Cheers,

Olivier

Revision history for this message
Josh McFayden (mcfayden) said :
#8

Thanks Olivier!