Changes between Version 1 and Version 2 of MW_development


Ignore:
Timestamp:
04/06/12 16:33:02 (8 years ago)
Author:
trac
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • MW_development

    v1 v2  
    2828Test of V2.5 against 2.0:
    2929||  '''Process''' || param || '''2.0'''  ||  '''2.5''' || Comment ||
    30 ||  W production (->e-) || W at 80.4  ||  14.36273+-7.3e-05 || 14.36218 +-1.5e-05 || OK ||
     30||  W production (->e-) || W at 80.4  ||  14.36273+-7.3e-05 || 14.36218 +-1.5e-05 || OK ||
    3131||  ^ || W at 100 || 15.80653 +-8.8e-05 ||  15.80328 +-2.0-05 || OK ||
    3232||  ^ || W at 150 || 17.8815 +- 1.0e-04 ||  17.87489 +-3.0e-05 || OK ||
    33 || (W>lnu)(W>jj) || 80.4  ||  48.27716+- 7e-05 || 48.255909+-6.6e-05 || OK ||
     33|| (W>lnu)(W>jj) || 80.4  ||  48.27716+- 7e-05 || 48.255909+-6.6e-05 || OK ||
    3434|| ^ || 100  ||  46.64555+- 8e-05 || 46.632452+- 6.1e-05 || OK ||
    3535|| ^ || 150  ||  56.563823+- 6.5e-05 || 56.575452+-6.9e-05 || OK ||
     
    4949
    5050|| '''process''' || '''block''' || V2.5 || true value || remarque ||||
    51 || pp>e+e-u || '''A'''  || (6.2952+-0.0017)E-5  || 6.299E-5 ||||||
    52 || bg>(t>b(W>jj))jj  || A+E || 149301+-282  || 149748 || ||||
    53 || jj(W>jj)jj  || A+D  || 3.2503+- 0.0046  || 3.254 ||||||
    54 || jj(W>eve)  || A+C  || 3.244+- 0.008  || 3.254 ||||||
    55 || jj(t>b(W>eve))  || A+B  || 150800+-1900  || 149748 || ||||
    56 || jj(W>bt) with (t>b(W>eve)) || A+A || (5.233+-0.004)E+19  || 5.22E9 || ||||
    57 || bg>(t>b(W>jj))jj  || multi-channel  || 148500+- 1300  || 149748 || (A+DD/A+E/A+11) (5e-17/1e-16/148500) || ||
    58 || bg>(t>b(W>jj))jj  || multi-channel  || || 149748 || (A+DD/A+E)  ||||
    59 || (W>lve)  || '''B'''  || (1.32 +-0.08)E-8 || 1.28E-8 || instability in the jacobian (if pz(l)=E(l)) ||||
    60 || (W>lve)j  || B  || (6.2927 +- 0.0039)E-5  || 6.299E-5 || non sensitive to the jac=0 area ||||
    61 ||  uu~>(W+>b~(t>b(W+>e+ve)))(W->e-ve~) || B+A  ||  5.201 +-0.046 E9 || 5.22E9 || ||||
     51|| pp>e+e-u || '''A'''  || (6.2952+-0.0017)E-5  || 6.299E-5 ||||||
     52|| bg>(t>b(W>jj))jj  || A+E || 149301+-282  || 149748 || ||||
     53|| jj(W>jj)jj  || A+D  || 3.2503+- 0.0046  || 3.254 ||||||
     54|| jj(W>eve)  || A+C  || 3.244+- 0.008  || 3.254 ||||||
     55|| jj(t>b(W>eve))  || A+B  || 150800+-1900  || 149748 || ||||
     56|| jj(W>bt) with (t>b(W>eve)) || A+A || (5.233+-0.004)E+19  || 5.22E9 || ||||
     57|| bg>(t>b(W>jj))jj  || multi-channel  || 148500+- 1300  || 149748 || (A+DD/A+E/A+11) (5e-17/1e-16/148500) || ||
     58|| bg>(t>b(W>jj))jj  || multi-channel  || || 149748 || (A+DD/A+E)  ||||
     59|| (W>lve)  || '''B'''  || (1.32 +-0.08)E-8 || 1.28E-8 || instability in the jacobian (if pz(l)=E(l)) ||||
     60|| (W>lve)j  || B  || (6.2927 +- 0.0039)E-5  || 6.299E-5 || non sensitive to the jac=0 area ||||
     61||  uu~>(W+>b~(t>b(W+>e+ve)))(W->e-ve~) || B+A  ||  5.201 +-0.046 E9 || 5.22E9 || ||||
    6262|| || B+B  || 149000+-500 || 149748 || ||||
    63 || (W>lv)(W>lv) || B+C  || 3.238+- 0.007 || 3.254 ||||||
    64 || (W>lv)(Z>mu+mu-) || B+D || 3.24 +- 0.05  || 3.254 ||||||
    65 || bg>(t>b(W+>jj))(W->e-ve~) || B+E || 148100+- 1100 || 149748  ||||||
    66 || (W>lv)(Z>mu+mu-) || B+fuse || 3.24 +- 0.04 || 3.254 ||||||
    67 ||  (W>lv)(Z>mu+mu-) || multi-channel  || 3.254+-0.003 || 3.254  ||||||
    68 || (W>svt~(sta->ta-n1)) || B+C (both massif)  || (5.959+-0.006)E-5 || 5.967E-5 ||||||
    69 || pp>(t>b(W+>e+ve)) || '''C'''  || 1.28E-1  || 6.2982E-005 || WRONG ||||
     63|| (W>lv)(W>lv) || B+C  || 3.238+- 0.007 || 3.254 ||||||
     64|| (W>lv)(Z>mu+mu-) || B+D || 3.24 +- 0.05  || 3.254 ||||||
     65|| bg>(t>b(W+>jj))(W->e-ve~) || B+E || 148100+- 1100 || 149748  ||||||
     66|| (W>lv)(Z>mu+mu-) || B+fuse || 3.24 +- 0.04 || 3.254 ||||||
     67||  (W>lv)(Z>mu+mu-) || multi-channel  || 3.254+-0.003 || 3.254  ||||||
     68|| (W>svt~(sta->ta-n1)) || B+C (both massif)  || (5.959+-0.006)E-5 || 5.967E-5 ||||||
     69|| pp>(t>b(W+>e+ve)) || '''C'''  || 1.28E-1  || 6.2982E-005 || WRONG ||||
    7070|| tt~ || '''D'''  ||  5.247 +- 0.021 E+9  || 5.227E+09 || ||||
    71 || H>WW>e+e-veve~ || '''E'''  || 531.567756+-121.365377  || 3.254 || instability in jac ???? ||||
    72 || H>WWj>e+e-veve~j || E  || || 149753 || ||||
    73 || WW>e+e-veve~ || '''F'''  || 565.237354+-123.048905  || 3.254 || instability in jac ????? ||||
    74 || WWj>e+e-veve~j || E  || || 149753 || ||||
    75 || bu > td; t>benu || C || 0.0468 || 3.257 || in progress ||||
    76 || bu > td; t>benu || B || 3.25 || 3.257 || in progress ||||
     71|| H>WW>e+e-veve~ || '''E'''  || 531.567756+-121.365377  || 3.254 || instability in jac ???? ||||
     72|| H>WWj>e+e-veve~j || E  || || 149753 || ||||
     73|| WW>e+e-veve~ || '''F'''  || 565.237354+-123.048905  || 3.254 || instability in jac ????? ||||
     74|| WWj>e+e-veve~j || E  || || 149753 || ||||
     75|| bu > td; t>benu || C || 0.0468 || 3.257 || in progress ||||
     76|| bu > td; t>benu || B || 3.25 || 3.257 || in progress ||||
    7777
    7878Cross-section Check
    7979|| '''process''' || '''block''' || ratio ME/MW ||
    80 || (W>lv)  || B  || 0.9991 +- 0.0014 ||
    81 || (W>svt~(sta->ta-n1)) || B+C || 1.003 +- 0.006 ||
     80|| (W>lv)  || B  || 0.9991 +- 0.0014 ||
     81|| (W>svt~(sta->ta-n1)) || B+C || 1.003 +- 0.006 ||
    8282|| SMU SMU || F || 1.020 +- 0.005 ||
    83 || tt> fully leptonique || D || 1.00056 +- 0.025 ||
     83|| tt> fully leptonique || D || 1.00056 +- 0.025 ||
    8484
    8585== Code V2.0 ==
     
    118118This code will be specific to the topology at work, and generated by the supra-code. The aim of the main code is to generate a phase-space point and compute the associated phase-space weight and jacobian. This will be done according to the following steps:
    119119
    120    * generation of invariant masses linked to factorized propagators. This generation will occur in an ascendant way so that we avoid unphysical regions as much as possible. This is particularly useful for topologies like H>WW*, where at least one of the W is off-shell.
     120   * generation of invariant masses linked to factorized propagators. This generation will occur in an ascendant way so that we avoid unphysical regions as much as possible. This is particularly useful for topologies like H>WW*, where at least one of the W is off-shell.
    121121
    122122   * generation of momenta that are factorized (according to the decision made by the supra-code) according to the transfer function. '''We should generate points according to Gaussian distributions'''
     
    154154==== IV. Differential weight ====
    155155
    156 A histrogram package (dbook) has been set up in MadWeight, in order to draw differential distributions of the weight. One first application is the reconstruction of top quark pair invariant mass in tt>2l events.
     156A histrogram package (dbook) has been set up in MadWeight, in order to draw differential distributions of the weight. One first application is the reconstruction of top quark pair invariant mass in tt>2l events.
    157157
    158158-
     
    165165
    166166||  '''Process''' || '''status'''  ||||
    167 || gg> (h>(w+>e+ ve )(w- > mu- vm~ ) )  ||  OK! || Olivier ||
    168 || gg<script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/themes/advanced/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikibuttons/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikiimage/langs/en.js"></script>8;gt;( h&gt;(w+&gt;e+ ve )(w- &gt; e- ve~ ) )  || -&gt; event  || Olivier ||
    169 || gg&gt;( h&gt;(w+&gt;e+ ve (w- &gt; e- ve~ ) )  ||  -&gt; event  || Olivier ||
    170 || pp&gt;(er+&gt;e+ n1) (er-&gt;e- n1)  ||  OK! || Pierre ||
    171 ||  pp&gt;(mur+&gt;mu+ n1) (mur-&gt;mu- <script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/themes/advanced/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikibuttons/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikiimage/langs/en.js"></script>n1) )  ||  not checked yet || Pierre ||
    172 || pp&gt;(t&gt;b(W+&gt;e+ve))(t~&gt;b~(W-&gt;mu-vm~)) ||  Ok! || Pierre ||
    173 || pp&gt;(t&gt;b(W+&gt;mu+vm)) (t~&gt;b~(W-&gt;mu-vm~)) || not checked yet ||||
    174 || pp&gt;(t&gt;b(W+&gt;mu+vm)) (t~&gt;b~(W-&gt;mu-vm~)) ||  not checked yet ||||
    175 || pp&gt;(t&gt;b(W+&gt;mu+vm)) (t~&gt;b~(W-&gt;jj)) ||  not checked yet ||||
    176 || pp&gt;(t&gt;b(W+&gt;jj)) (t~&gt;b~(W-&gt;jj)) || not checked yet ||||
     167|| gg> (h>(w+>e+ ve )(w- > mu- vm~ ) )  ||  OK! || Olivier ||
     168|| gg<script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/themes/advanced/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikibuttons/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikiimage/langs/en.js"></script>8;gt;( h>(w+>e+ ve )(w- > e- ve~ ) )  || -> event  || Olivier ||
     169|| gg>( h>(w+>e+ ve (w- > e- ve~ ) )  ||  -> event  || Olivier ||
     170|| pp>(er+>e+ n1) (er->e- n1)  ||  OK! || Pierre ||
     171||  pp>(mur+>mu+ n1) (mur->mu- <script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/themes/advanced/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikibuttons/langs/en.js"></script><script type="text/javascript" src="http://cp3wks05.fynu.ucl.ac.be/twiki/pub/TWiki/TinyMCEPlugin/tinymce/jscripts/tiny_mce/plugins/twikiimage/langs/en.js"></script>n1) )  ||  not checked yet || Pierre ||
     172|| pp>(t>b(W+>e+ve))(t~>b~(W->mu-vm~)) ||  Ok! || Pierre ||
     173|| pp>(t>b(W+>mu+vm)) (t~>b~(W->mu-vm~)) || not checked yet ||||
     174|| pp>(t>b(W+>mu+vm)) (t~>b~(W->mu-vm~)) ||  not checked yet ||||
     175|| pp>(t>b(W+>mu+vm)) (t~>b~(W->jj)) ||  not checked yet ||||
     176|| pp>(t>b(W+>jj)) (t~>b~(W->jj)) || not checked yet ||||
    177177
    178178==== Things to do for the next Release ====
     
    181181
    182182==== Different Idea ====
    183    1. '''Grid initialization''': <br /> The transfer functions coud(should?) have their own grid this could be achieved by integrating the one-dimensional transfer functions. If this is needed (probably only for some hard cases), this could be done on an event-by-event basis. In other case a typical events could do the job once.
    184    1. '''Tau decay''';<br /> it could be possible to have a special change of variable for the tau&gt;lepton decay, using the coliniarity of the decay. A starting grid for the neutrinos angles could be used if we want to stay completely general.
    185    1. '''Transfer function''':<br /> we could use a supra code trick for the transfer functions. This will give gain of time (~1%) but the main goal is to simplify the way to add a transfer functions in special case. -&gt; '''Done in Version 2.5'''
    186    1. '''Permutation''':<br /> we could/should develop a python code treating the permutation (supressing the identical permutation). This could be achieved with a MG stand alone test.
    187    1. '''grid for solution choices''': <br />For the moment each solution (in change of variables) has the same probability to choose. Adding an adaptative grid on this could be interesting (-&gt;include those variable in VEGAS). An other idea will be to use the jacobian in order to improve our choices
    188    1. '''Adaptative precision''':<br />For the moment the weight in each SubProcess the weights are computed with the same precision, in a normalized mode or in order to have refine procedure, a adaptative precision could be cut the processing time by a significative factor.
     183   1. '''Grid initialization''': [[br]] The transfer functions coud(should?) have their own grid this could be achieved by integrating the one-dimensional transfer functions. If this is needed (probably only for some hard cases), this could be done on an event-by-event basis. In other case a typical events could do the job once.
     184   1. '''Tau decay''';[[br]] it could be possible to have a special change of variable for the tau>lepton decay, using the coliniarity of the decay. A starting grid for the neutrinos angles could be used if we want to stay completely general.
     185   1. '''Transfer function''':[[br]] we could use a supra code trick for the transfer functions. This will give gain of time (~1%) but the main goal is to simplify the way to add a transfer functions in special case. -> '''Done in Version 2.5'''
     186   1. '''Permutation''':[[br]] we could/should develop a python code treating the permutation (supressing the identical permutation). This could be achieved with a MG stand alone test.
     187   1. '''grid for solution choices''': [[br]]For the moment each solution (in change of variables) has the same probability to choose. Adding an adaptative grid on this could be interesting (->include those variable in VEGAS). An other idea will be to use the jacobian in order to improve our choices
     188   1. '''Adaptative precision''':[[br]]For the moment the weight in each SubProcess the weights are computed with the same precision, in a normalized mode or in order to have refine procedure, a adaptative precision could be cut the processing time by a significative factor.
    189189
    190190-- Main.OlivierMattelaer - 24 Nov 2008
     
    192192
    193193
     194