2013:groups:tools_lheextension
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision | |||
|
2013:groups:tools_lheextension [2013/09/20 10:07] benjamin.fuks |
2013:groups:tools_lheextension [2013/10/31 15:09] (current) benjamin.fuks |
||
|---|---|---|---|
| Line 12: | Line 12: | ||
| * a reorganization of the weights, collecting them into categories. | * a reorganization of the weights, collecting them into categories. | ||
| - | The current version of the proposal is currently being implemented in the NLO tools, MadGraph 5 as well as in MadAnalysis 5. Changes are of course still possible. **Any comment, suggestion or criticism is welcome. Please update [[#the_proposal|the text]]. Points to be discussed can be found [[#To_be_discussed|below]] the proposal.** | + | The current version of the proposal has been implemented in MadGraph5_aMC@NLO, as well as in MadAnalysis 5. |
| Line 26: | Line 26: | ||
| <initrwgt> | <initrwgt> | ||
| <weight id='1'> This is the original event weight </weight> | <weight id='1'> This is the original event weight </weight> | ||
| - | <weightgroup type='scale_variation' combine='envelope'> | + | <weightgroup name='scale_variation' combine='envelope'> |
| <weight id='2'> muR=2.0 </weight> | <weight id='2'> muR=2.0 </weight> | ||
| <weight id='3'> muR=0.5 </weight> | <weight id='3'> muR=0.5 </weight> | ||
| </weightgroup> | </weightgroup> | ||
| - | <weightgroup type="mrst2008e40" combine="hessian"> | + | <weightgroup name="mrst2008e40" combine="hessian"> |
| <weight id='4'> set01 </weight> | <weight id='4'> set01 </weight> | ||
| <weight id='5'> set02 </weight> | <weight id='5'> set02 </weight> | ||
| ... | ... | ||
| </weightgroup> | </weightgroup> | ||
| - | <weightgroup type='Qmatch_variation' combine='envelope'> | + | <weightgroup name='Qmatch_variation' combine='envelope'> |
| <weight id='44'> Qmatch=20 </weight> | <weight id='44'> Qmatch=20 </weight> | ||
| <weight id='45'> Qmatch=40 </weight> | <weight id='45'> Qmatch=40 </weight> | ||
| Line 46: | Line 46: | ||
| This information in the header should be human-readable and explain what the weights with the corresponding identifiers mean. It can simply contain all the parameters that were used in to generate this weight; or only the ones that were changed compared to the original run; or simply a sentence explaining what this number means. It's up to the user that is doing the analysis to make sure that this information is correctly used (and up to the authors of the codes to make sure that the user has enough information to understand what the weights correspond to). | This information in the header should be human-readable and explain what the weights with the corresponding identifiers mean. It can simply contain all the parameters that were used in to generate this weight; or only the ones that were changed compared to the original run; or simply a sentence explaining what this number means. It's up to the user that is doing the analysis to make sure that this information is correctly used (and up to the authors of the codes to make sure that the user has enough information to understand what the weights correspond to). | ||
| - | The weightgroup tag allows to group several weights together (to have the information about how to combine weights to obtain, //e.g.//, scale variation or pdf uncertainties. The attributes ''combine'' is optional. It indicates how to combine the uncertainties. Possible arguments are ''none'', ''hessian'', ''envelope'' or ''gaussian''. If not specified, the default choice is ''combine='none''', all the curves associated with each weights being kept independent. For ''combine='hessian''', the first weight is the central value and the next weights correspond respectively to the positive and negative variations along a specific direction of the parameter space. This is not very XML-friendly, so any suggestion here is very welcome. | + | The weightgroup tag allows to group several weights together (to have the information about how to combine weights to obtain, //e.g.//, scale variation or pdf uncertainties. The attributes ''combine'' is optional. It indicates how to combine the uncertainties. Possible arguments are ''none'', ''hessian'', ''envelope'' or ''gaussian''. If not specified, the default choice is ''combine='none''', all the curves associated with each weights being kept independent. For ''combine='hessian''', the first weight is the central value and the next weights correspond respectively to the positive and negative variations along a specific direction of the parameter space. |
| ===Second part: within each event=== | ===Second part: within each event=== | ||
| We start with an example: | We start with an example: | ||
| <code> | <code> | ||
| - | <event id='evtid'> | + | <event> |
| 7 100 0.10000000E+01 0.20000000E+00 0.00000000E+00 0.00000000E+00 | 7 100 0.10000000E+01 0.20000000E+00 0.00000000E+00 0.00000000E+00 | ||
| -2 -1 0 0 0 0 0.12699952E+01 0.55429630E+01 0.57634577E+02 0.57914435E+02 0.00000000E+00 0. 0. | -2 -1 0 0 0 0 0.12699952E+01 0.55429630E+01 0.57634577E+02 0.57914435E+02 0.00000000E+00 0. 0. | ||
| Line 61: | Line 61: | ||
| 13 1 0 0 0 0 -0.48048241E+01 -0.14284099E+02 -0.31291705E+02 0.34731884E+02 0.10566000E+00 0. 0. | 13 1 0 0 0 0 -0.48048241E+01 -0.14284099E+02 -0.31291705E+02 0.34731884E+02 0.10566000E+00 0. 0. | ||
| <rwgt> | <rwgt> | ||
| - | <wgt id='1'> 1.001e+00 </wgt> | + | <wgt id='1'> 1.000e+00 </wgt> |
| <wgt id='2'> 0.204e+00 </wgt> | <wgt id='2'> 0.204e+00 </wgt> | ||
| <wgt id='3'> 1.564e+00 </wgt> | <wgt id='3'> 1.564e+00 </wgt> | ||
| Line 72: | Line 72: | ||
| The numbers should be normalized in the same way as the original weight of the event, //i.e.// if the weights sum-up to the total cross section, also the new ''<nw>'' weights should sum up to the total cross section (which is in general slightly different because different parameters were used). | The numbers should be normalized in the same way as the original weight of the event, //i.e.// if the weights sum-up to the total cross section, also the new ''<nw>'' weights should sum up to the total cross section (which is in general slightly different because different parameters were used). | ||
| - | |||
| - | The event ''id'' is important as soon as we have to deal with event files after showering/hadronization (stdHEP or HEPMC event files). This would allow to pass the reweighing information under the price of having multiple file to read. This also avoids to have to extend the standard stdHEP and HEPMC formats. When used jointly with a stdHEP or HEPMC file, the LHE file only contains the reweighting information. | ||
| - | ==== To be discussed ==== | ||
| - | * The ''<clustering>'' tag of 1003.1643 | ||
2013/groups/tools_lheextension.1379664449.txt.gz ยท Last modified: 2013/09/20 10:07 by benjamin.fuks
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: Public Domain
