LHEF for Matching

From Wiki Les Houches 09

(Difference between revisions)
Jump to: navigation, search
Line 42: Line 42:
definition one-way, and the information can be made more easily
definition one-way, and the information can be made more easily
provided. The suggested attributes for the <tt>xsecinfo</tt> tag are as follows
provided. The suggested attributes for the <tt>xsecinfo</tt> tag are as follows
-
(R means the attribute is required. 'D=' means a default is assumed if
+
(R means the attribute is required. ''D='' means a default is assumed if
the attribute is absent).
the attribute is absent).
Line 50: Line 50:
* <tt>meanweight</tt> (D=1): The average weight of the events in the file (same unit as maxweight.
* <tt>meanweight</tt> (D=1): The average weight of the events in the file (same unit as maxweight.
* <tt>negweights</tt> (D=no): If yes, the file may contain events with negative weights.
* <tt>negweights</tt> (D=no): If yes, the file may contain events with negative weights.
-
* <tt>varweights</tt> (D=no): If yes, the file may contain varying weights (if no, all events are weighted with maxweight=meanweight).
+
* <tt>varweights</tt> (D=no): If yes, the file may contain varying weights (if no, all events are weighted with <tt>maxweight=meanweight</tt>).
-
The information given per process in the HEPRUP common block can then
+
The information given per process in the <tt>HEPRUP</tt> common block can then
safely be ignored and the cross section for each process is assumed to
safely be ignored and the cross section for each process is assumed to
-
be ''sumprocw*totxsec/maxweight/neve''.
+
be <tt>sumprocw*totxsec/maxweight/neve</tt>.
-
== The ''cutsinfo'' tag (optional) ==
+
== The <tt>cutsinfo</tt> tag (optional) ==
Here the MEG should supply information about the cuts used in the
Here the MEG should supply information about the cuts used in the
-
generation. Inside this tag any number of ''cut'' tags can be supplied.
+
generation. Inside this tag any number of <tt>cut</tt> tags can be supplied.
-
=== The cut tag ===
+
=== The <tt>cut</tt> tag ===
The cut tag has information of a particular cut used. The information
The cut tag has information of a particular cut used. The information
Line 68: Line 68:
the maximum and/or minimum value of that parameter. The attributes are
the maximum and/or minimum value of that parameter. The attributes are
-
* type (R): This defines the variable which is cut. The following values are allowed (the lab frame is assumed in all cases):
+
* <tt>type<tt> (R): This defines the variable which is cut. The following values are allowed (the lab frame is assumed in all cases):
** "invmass" the invariant mass of two particles (or if only one particle type is given, the rest mass of that particle).
** "invmass" the invariant mass of two particles (or if only one particle type is given, the rest mass of that particle).
** "kt": the transverse momenta of a particle (in GeV)
** "kt": the transverse momenta of a particle (in GeV)
Line 77: Line 77:
** "ETmiss": the sum of the pt of the neutrinos in the process
** "ETmiss": the sum of the pt of the neutrinos in the process
** ...
** ...
-
* p1 (D=any): The particle for which this cut applies. This can either be a number corresponding to a given particle PDG code or one of the following
+
* <tt>p1</tt> (D=any): The particle for which this cut applies. This can either be a number corresponding to a given particle PDG code or one of the following
** "any": any particle
** "any": any particle
** "jet": any coloured particle
** "jet": any coloured particle
Line 85: Line 85:
** "l" ("l+" or "l-"): any (positively or negatively) charged lepton
** "l" ("l+" or "l-"): any (positively or negatively) charged lepton
** ...
** ...
-
* p2 (D=0): If cut is between pairs of particles, this attribute should be non-zero. Allowed values are the same as for ''p1''
+
* </tt><tt>p2 (D=0): If cut is between pairs of particles, this attribute should be non-zero. Allowed values are the same as for <tt>p1</tt>
-
* min (D=-inf): The minimum value of this variable
+
* <tt>min</tt> (D=-inf): The minimum value of this variable
-
* max (D=+inf): The maximum value of this variable
+
* <tt>max</tt> (D=+inf): The maximum value of this variable
-
== The ''procinfo'' tag (optional) ==
+
== The <tt>procinfo</tt> tag (optional) ==
-
For each process number used in the NPRUP variable in the HEPEUP
+
For each process number used in the <tt>NPRUP</tt> variable in the <tt>HEPEUP</tt>
common block we can have additional information given in the following
common block we can have additional information given in the following
-
attributeds:
+
attributes:
-
* nprup (D=0): The process number for which the information is given. 0 means all processes.
+
* <tt>nprup</tt> (D=0): The process number for which the information is given. 0 means all processes.
-
* loops (D=-1): The number of loops used in calculating this process.
+
* <tt>loops</tt> (D=-1): The number of loops used in calculating this process.
-
* qcdorder (D=-1): The number of QCD vertices used in calculating this process.
+
* <tt>qcdorder</tt> (D=-1): The number of QCD vertices used in calculating this process.
-
* eworder (D=-1): The number of electro-weak vertices used in calculating this process.
+
* <tt>eworder</tt> (D=-1): The number of electro-weak vertices used in calculating this process.
In all cases above ''-1'' means that no information is given.
In all cases above ''-1'' means that no information is given.
-
== The ''mergeinfo'' tag (optional) ==
+
== The <tt>mergeinfo</tt> tag (optional) ==
For some merging schemes (eg. for CKKW) it is possible to reweight the
For some merging schemes (eg. for CKKW) it is possible to reweight the
the events with Sudakov form factors already in the MEG. This can be
the events with Sudakov form factors already in the MEG. This can be
described in the following attributes.
described in the following attributes.
-
* nprup (D=0): The process number for which the information is given. "0" means all processes.
+
* <tt>nprup</tt> (D=0): The process number for which the information is given. "0" means all processes.
-
* reweightscheme (R): valid values are
+
* <tt>reweightscheme</tt> (R): valid values are
** "None": no reweighting done
** "None": no reweighting done
** "CKKW-NLL": events are reweighted by the standard CKKW NLL Sudakov form factor
** "CKKW-NLL": events are reweighted by the standard CKKW NLL Sudakov form factor
** "CKKW-NLL1": same as above but weight is required to be less than 1
** "CKKW-NLL1": same as above but weight is required to be less than 1
** "VetoSudakov": events are reweighted with 0 or 1 with a probability given by a (non-analytic) form factor.
** "VetoSudakov": events are reweighted with 0 or 1 with a probability given by a (non-analytic) form factor.
-
* reconstruction (R): The clustering used to find the scales in the Sudakov form factors (and alpha_S). The following values are valid
+
* <tt>reconstruction</tt> (R): The clustering used to find the scales in the Sudakov form factors (and alpha_S). The following values are valid
** "ktalg": the kt-algorithm used in the original CKKW scheme
** "ktalg": the kt-algorithm used in the original CKKW scheme
** ...
** ...
-
* maxmult (D=no): If yes the corresponding process is reweighted as if it is the maximum multiplicity process (ie. the Sudakov for the last step down to the merging scale is not included.
+
* <tt>maxmult</tt> (D=no): If yes the corresponding process is reweighted as if it is the maximum multiplicity process (ie. the Sudakov for the last step down to the merging scale is not included.
-
== The ''nloinfo'' tag (optional) ==
+
== The <tt>nloinfo</tt> tag (optional) ==
If the events were produced by a NLO generator the following
If the events were produced by a NLO generator the following
information can be supplied as attributes:
information can be supplied as attributes:
-
* scheme (R): The NLO (subtraction) scheme used. The following values are valid:
+
* <tt>scheme</tt> (R): The NLO (subtraction) scheme used. The following values are valid:
** "CSdipole": Catani-Seymour subtraction.  
** "CSdipole": Catani-Seymour subtraction.  
** "POWHEG": Nasons scheme for giving one extra jet with Sudakov form factors correct to NLO.
** "POWHEG": Nasons scheme for giving one extra jet with Sudakov form factors correct to NLO.
** "MC@NLO": "parton-shower subtracted"
** "MC@NLO": "parton-shower subtracted"
** "Exclusive": The 1-loop events are generated according to the exclusive NLO cross section (according to the cuts), while 0-loop events are pure tree-level.
** "Exclusive": The 1-loop events are generated according to the exclusive NLO cross section (according to the cuts), while 0-loop events are pure tree-level.
-
* rfscheme (D=MSbar): The factorization and renormalization scheme used (MSbar, DIS)
+
* <tt>rfscheme</tt> (D=MSbar): The factorization and renormalization scheme used (MSbar, DIS)
Note that for NLO generators (except for the Exclusive scheme), events
Note that for NLO generators (except for the Exclusive scheme), events
with real emission corrections will be produced which will not respect
with real emission corrections will be produced which will not respect
the cuts. For the POWHEG scheme there is an additional cut imposed on
the cuts. For the POWHEG scheme there is an additional cut imposed on
-
the emission, which can be specified by "cut" tags (see the "cutsinfo"
+
the emission, which can be specified by "cut" tags (see the <tt>cutsinfo</tt>
-
tag above) inside the the nloinfo tag.
+
tag above) inside the the <tt>nloinfo</tt> tag.
-
== The ''pdferrors'' tag (optional) ==
+
== The <tt>pdferrors</tt> tag (optional) ==
It is sometimes efficient to supply different weights for the same
It is sometimes efficient to supply different weights for the same
Line 143: Line 143:
have this information and the only attribute in this tag specifies the
have this information and the only attribute in this tag specifies the
number of sets used.
number of sets used.
-
* nsets (R): the number of error sets used.
+
* <tt>nsets</tt> (R): the number of error sets used.
= Event information =
= Event information =
This is information about a particular event to be included inside the
This is information about a particular event to be included inside the
-
"event" tag after the compulsory block of parameters corresponding to
+
<tt>event</tt> tag after the compulsory block of parameters corresponding to
-
the HEPEUP common block.
+
the <tt>HEPEUP</tt> common block.
-
== The ''weightinfo'' tag (optional) ==
+
== The <tt>weightinfo</tt> tag (optional) ==
Here some information can be supplied to help the PSG to reweight
Here some information can be supplied to help the PSG to reweight
events according to its needs. The following attributes may be given
events according to its needs. The following attributes may be given
-
* weight (D=XWGTUP): the weight of the event as described in the xsecinfo tag. The default is the XWGTUP variable of the HEPEUP common block.
+
* <tt>weight</tt> (D=XWGTUP): the weight of the event as described in the <tt>xsecinfo</tt> tag. The default is the <tt>XWGTUP</tt> variable of the <tt>HEPEUP</tt> common block.
-
* born (D=0): the relative size of the tree-level cross section for this event. Summing ''born*weight'' (times ''totxsec/neve/maxweight'') will give the born cross section for this process. If 0 this information is not given.
+
* <tt>born</tt> (D=0): the relative size of the tree-level cross section for this event. Summing <tt>born*weight</tt> and multiplying with <tt>totxsec/neve/maxweight</tt> will give the born cross section for this process. If 0 this information is not given.
-
* sudakov (D=0): The sudakov form factor used to weight this event (see the "reweightscheme" attribute of the "mergeinfo" tag).
+
* <tt>sudakov</tt> (D=0): The sudakov form factor used to weight this event (see the <tt>reweightscheme</tt> attribute of the <tt>mergeinfo</tt> tag).
-
* subtraction (D=0): The relative size of the subtraction term used for a 0-loop event (is typically negative)
+
* <tt>subtraction</tt> (D=0): The relative size of the subtraction term used for a 0-loop event (is typically negative)
-
The relationships between these terms are typically ''born*sudakov=1'' and
+
The relationships between these terms are typically <tt>born*sudakov=1</tt> and
-
''born+subtraction=1''.
+
<tt>born+subtraction=1</tt>.
-
== The ''alphas'' tag (optional) ==
+
== The <tt>alphas</tt> tag (optional) ==
-
Information about the alpha_S used in generating this
+
Information about the <math>\alpha_S</math> used in generating this
event. Attributes are
event. Attributes are
-
* scale (D=SCALUP): The scale in GeV used in the alpha_S (default is the SCALUP variable in the HEPEUP common block).
+
* <tt>scale</tt> (D=SCALUP): The scale in GeV used in the alpha_S (default is the <tt>SCALUP</tt> variable in the <tt>HEPEUP</tt> common block).
-
* value (D=AQCDUP): The value of alpha_S used (default is the AQCDUP variable in the HEPEUP common block).
+
* <tt>value</tt> (D=AQCDUP): The value of alpha_S used (default is the <tt>AQCDUP</tt> variable in the <tt>HEPEUP</tt> common block).
If different values of alpha_s has been used in different vertices in
If different values of alpha_s has been used in different vertices in
this even a white-space separated list of values can be given in this
this even a white-space separated list of values can be given in this
tag. The order of the values is dictated by the "reconstruction"
tag. The order of the values is dictated by the "reconstruction"
-
attribute in the "mergeinfo" tag above (falling scales) if present.
+
attribute in the <tt>mergeinfo</tt> tag above (falling scales) if present.
-
== The ''pdfinfo'' tag (optional) ==
+
== The <tt>pdfinfo</tt> tag (optional) ==
Here the MEG can supply information about the PDFs used for this
Here the MEG can supply information about the PDFs used for this

Revision as of 13:13, 10 June 2009

Contents

Introduction

The Les Houches Event File accord (LHEF) was introduced in 2006 and was meant as a standard way of representing in an ASCII file the information in the old Les Houches common block accord from 2001. LHEF was constructed using XML tags to be easy to extend (although some additional structure is assumed inside some tags which is not formulated in XML). The standard has been extremely usefull, and has been used a lot to interface matrix element generators and parton shower programs.

The standards are found here:

As the matching and merging of tree-level matrix elements and parton showers are now becoming the state-of-art, it is reasonable to let this be reflected in an updated file format to standardise how relevant information should be given by the matrix element generators (MEGs) in a usable fashion for the parton shower programs (PSGs). Furthermore, looking ahead towards the time when NLO matching and merging with parton showers hopefully will become the standard, it might be worth thinking about which additional information is needed for the parton showers in this case.

Below you will find a rough draft of additional XML tags suggsted to be included in the new standard (version 2.0).

Run information

This is information about all events in the file to be included inside the "init" tag after the compulsary block of parameters corresponding to the HEPRUP common block.

The xsecinfo tag (required)

The information in the HEPRUP common block is in principle sufficient to figure out the cross sections of the processes involved. However, the way things are specified is a bit confusing and complicated since it was assumed to be used to pass information between the MEG and PSG in both directions. For the event file, the communication is per definition one-way, and the information can be made more easily provided. The suggested attributes for the xsecinfo tag are as follows (R means the attribute is required. D= means a default is assumed if the attribute is absent).

  • neve (R): the number of events in the file
  • totxsec (R): the total cross section (in units of pb) of all processes in the file
  • maxweight (D=1): the maximum weight of any event in the file (in an arbitrary unit)
  • meanweight (D=1): The average weight of the events in the file (same unit as maxweight.
  • negweights (D=no): If yes, the file may contain events with negative weights.
  • varweights (D=no): If yes, the file may contain varying weights (if no, all events are weighted with maxweight=meanweight).

The information given per process in the HEPRUP common block can then safely be ignored and the cross section for each process is assumed to be sumprocw*totxsec/maxweight/neve.


The cutsinfo tag (optional)

Here the MEG should supply information about the cuts used in the generation. Inside this tag any number of cut tags can be supplied.

The cut tag

The cut tag has information of a particular cut used. The information needed is which particle(s) are affected, what variable is used and the maximum and/or minimum value of that parameter. The attributes are

  • type<tt> (R): This defines the variable which is cut. The following values are allowed (the lab frame is assumed in all cases):
    • "invmass" the invariant mass of two particles (or if only one particle type is given, the rest mass of that particle).
    • "kt": the transverse momenta of a particle (in GeV)
    • "eta": the pseudorapidity of a particle
    • "y": the true rapidity of a particle
    • "deltaR": the pseudorapidity--azimuthal-angle difference between two particles.
    • "E": the energy of a particle
    • "ETmiss": the sum of the pt of the neutrinos in the process
    • ...
  • <tt>p1 (D=any): The particle for which this cut applies. This can either be a number corresponding to a given particle PDG code or one of the following
    • "any": any particle
    • "jet": any coloured particle
    • "uds": ("udsc") any light quark
    • "guds": ("gudsc") any light quark or gluon
    • "b" ("c"): any b or bbar (c or cbar) quark
    • "l" ("l+" or "l-"): any (positively or negatively) charged lepton
    • ...
  • </tt>p2 (D=0): If cut is between pairs of particles, this attribute should be non-zero. Allowed values are the same as for <tt>p1
  • min (D=-inf): The minimum value of this variable
  • max (D=+inf): The maximum value of this variable


The procinfo tag (optional)

For each process number used in the NPRUP variable in the HEPEUP common block we can have additional information given in the following attributes:

  • nprup (D=0): The process number for which the information is given. 0 means all processes.
  • loops (D=-1): The number of loops used in calculating this process.
  • qcdorder (D=-1): The number of QCD vertices used in calculating this process.
  • eworder (D=-1): The number of electro-weak vertices used in calculating this process.

In all cases above -1 means that no information is given.

The mergeinfo tag (optional)

For some merging schemes (eg. for CKKW) it is possible to reweight the the events with Sudakov form factors already in the MEG. This can be described in the following attributes.

  • nprup (D=0): The process number for which the information is given. "0" means all processes.
  • reweightscheme (R): valid values are
    • "None": no reweighting done
    • "CKKW-NLL": events are reweighted by the standard CKKW NLL Sudakov form factor
    • "CKKW-NLL1": same as above but weight is required to be less than 1
    • "VetoSudakov": events are reweighted with 0 or 1 with a probability given by a (non-analytic) form factor.
  • reconstruction (R): The clustering used to find the scales in the Sudakov form factors (and alpha_S). The following values are valid
    • "ktalg": the kt-algorithm used in the original CKKW scheme
    • ...
  • maxmult (D=no): If yes the corresponding process is reweighted as if it is the maximum multiplicity process (ie. the Sudakov for the last step down to the merging scale is not included.


The nloinfo tag (optional)

If the events were produced by a NLO generator the following information can be supplied as attributes:

  • scheme (R): The NLO (subtraction) scheme used. The following values are valid:
    • "CSdipole": Catani-Seymour subtraction.
    • "POWHEG": Nasons scheme for giving one extra jet with Sudakov form factors correct to NLO.
    • "MC@NLO": "parton-shower subtracted"
    • "Exclusive": The 1-loop events are generated according to the exclusive NLO cross section (according to the cuts), while 0-loop events are pure tree-level.
  • rfscheme (D=MSbar): The factorization and renormalization scheme used (MSbar, DIS)

Note that for NLO generators (except for the Exclusive scheme), events with real emission corrections will be produced which will not respect the cuts. For the POWHEG scheme there is an additional cut imposed on the emission, which can be specified by "cut" tags (see the cutsinfo tag above) inside the the nloinfo tag.

The pdferrors tag (optional)

It is sometimes efficient to supply different weights for the same event using a sets of error PDFs. If this is done, each event must have this information and the only attribute in this tag specifies the number of sets used.

  • nsets (R): the number of error sets used.

Event information

This is information about a particular event to be included inside the event tag after the compulsory block of parameters corresponding to the HEPEUP common block.

The weightinfo tag (optional)

Here some information can be supplied to help the PSG to reweight events according to its needs. The following attributes may be given

  • weight (D=XWGTUP): the weight of the event as described in the xsecinfo tag. The default is the XWGTUP variable of the HEPEUP common block.
  • born (D=0): the relative size of the tree-level cross section for this event. Summing born*weight and multiplying with totxsec/neve/maxweight will give the born cross section for this process. If 0 this information is not given.
  • sudakov (D=0): The sudakov form factor used to weight this event (see the reweightscheme attribute of the mergeinfo tag).
  • subtraction (D=0): The relative size of the subtraction term used for a 0-loop event (is typically negative)

The relationships between these terms are typically born*sudakov=1 and born+subtraction=1.


The alphas tag (optional)

Information about the αS used in generating this event. Attributes are

  • scale (D=SCALUP): The scale in GeV used in the alpha_S (default is the SCALUP variable in the HEPEUP common block).
  • value (D=AQCDUP): The value of alpha_S used (default is the AQCDUP variable in the HEPEUP common block).

If different values of alpha_s has been used in different vertices in this even a white-space separated list of values can be given in this tag. The order of the values is dictated by the "reconstruction" attribute in the mergeinfo tag above (falling scales) if present.

The pdfinfo tag (optional)

Here the MEG can supply information about the PDFs used for this event. The following attributes are available

  • p1 (R): the type of the incoming particle 1.
  • p2 (R): the type of the incoming particle 2.
  • x1 (R): the x-value used for the incoming particle 1.
  • x2 (R): the x-value used for the incoming particle 2.
  • xf1 (R): the value of the pdf for the incoming particle 1.
  • xf2 (R): the value of the pdf for the incoming particle 2.
  • scale (D=SCALUP): The scale in GeV used in the PDFs (default is the SCALUP variable in the HEPEUP common block).

The pdfinfo tag may contain a white-space separated list of additional weights corresponding to reweighting this event according to the error sets of the LHAPDF set used.

Personal tools