Changes between Initial Version and Version 1 of MW_development


Ignore:
Timestamp:
03/20/12 16:17:50 (7 years ago)
Author:
trac
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • MW_development

    v1 v1  
     1
     2
     3== Low mass SM Higgs search ==
     4
     5LowMassSMHiggs
     6
     7== Ressource Material ==
     8
     9   1. CVS Version of the code: [http://cp3wks05.fynu.ucl.ac.be/cgi-bin/ViewVC/viewvc.cgi/MG_ME/?pathrev=MW2  version of the code]
     10   1. Ongoing Paper: [http://www.fynu.ucl.ac.be/viewvc/omatt/article/?root=cvsfynu  paper]
     11
     12== Code V2.5 ==
     13
     14After a new discussion (September 09) between Olivier and Pierre and inspired by a discussion to Fabio. We want to include the multi-channel techniques in MW. In the previous code (2.0), an overconstrained systems lead to unaligned peak, but also to different change of variable, different because they are not aligning the same peak. In consequence, the multi-channel techniques is possible to performed.
     15
     16In principle, the multi-channel is pointless because the correct area of the PS should be selected by any of those change of variables (the constraints should be consistent) but MW computes also weights for different permutation/background where those constraints are in conflicts, causing slow integration. For those events/permutations, the multi-channel can contibutes to fastenize the integral (a lot?)
     17
     18=== Development of Code ===
     19
     20In a first step, (in progress), we will keep the different changes of variables provides by the python supra-code and run in multi-channel in those cases. Some script are added in the Python supra-script in order to compute the weights associated to each channel. The 'MadEvent trick' is used in order to choose wisely the weights of each channel.
     21
     22A second step, will be to check to modify the supra-script in order to create additional change of variables in order to have all peaks aligned in at least one change of variables
     23
     24We could(should?) think to fuse the supra code with Madgraph5 in order to have a full multi-channel program (multi-channel on the square matrix element and on the peak position). But this will be MG5 dependant and has some drawback (impossibility to modify Transfer Function after the diagram generation)
     25
     26=== Test of V2.5 ===
     27
     28Test of V2.5 against 2.0:
     29||  '''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 ||
     31||  ^ || W at 100 || 15.80653 +-8.8e-05 ||  15.80328 +-2.0-05 || OK ||
     32||  ^ || 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 ||
     34|| ^ || 100  ||  46.64555+- 8e-05 || 46.632452+- 6.1e-05 || OK ||
     35|| ^ || 150  ||  56.563823+- 6.5e-05 || 56.575452+-6.9e-05 || OK ||
     36|| semi-lept tt || Mt=150 || 813.6807+-0.0015 || 813.6377+-0.0013 || OK ||
     37|| ^ ||  152.5 || 813.3573+-0.0343 || 812.5640 +-0.0032 || OK ||
     38|| ^ ||  155 || 812.6200+-0.0015 || 812.1788+-0.0024 || OK ||
     39|| ^ ||  157.5 || 812.8469+-0.0013 || 812.8812+-0.0015 || OK ||
     40|| ^ ||  160 || 814.1727+-0.0015 || 814.2527+-0.0047 || OK ||
     41|| ^ ||  162.5 || 815.5738+-0.0013 || 816.0374+-0.0037 || OK ||
     42|| ^ ||  165 || 819.1933+-0.0019 || 818.2402+-0.0014 || more confident in 2.5 ||
     43|| ^ ||  167.5 || 822.4768+-0.0026 || 821.9353+-0.0038 || OK ||
     44|| ^ ||  170 || 827.7806+-0.0013 || 825.9346+-0.0019 || more confident in 2.5  ||
     45
     46=== Check of version 2.5 ===
     47
     48Phase-Space Check
     49
     50|| '''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 || ||||
     62|| || 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 ||||
     70|| 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 ||||
     77
     78Cross-section Check
     79|| '''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 ||
     82|| SMU SMU || F || 1.020 +- 0.005 ||
     83|| tt> fully leptonique || D || 1.00056 +- 0.025 ||
     84
     85== Code V2.0 ==
     86
     87After a discussion (18 jan 08) between Olivier and Pierre, we have planned to write a new version of the code, with particular emphase on the generic and modular characters of the algorithm. In this new code, a supra-code (python script) will produce a phase-space generator (in fortran) that is specific to a certain topology. This version is currently on CVS (see link above).
     88
     89=== Advantages of V2.0 ===
     90
     91The main advantage of this new organization of the code resides in the factorization of two difficulties:
     92
     93   * Given a topology, find a phase space generator that is appropriate to ME method
     94   * write a code that select '''automatically''' the adequate phase-space parametrization for any topology
     95
     96In the new code, the main code (which is ''almost'' the only topology-dependent file) has a very simple structure, that can be handled easily, even written explicitly for a given topology. Only the way to generate it automatically is intricate. If we choose the supra code to generate this file, it is still very easy to check its contents.
     97
     98Also in this new version, it will be really easy to implement new blocks or new classes of ECS in the future.
     99
     100=== The supra-code (python script) ===
     101
     102This is for sure the most difficult part. The supra code will
     103
     104   * read the topology and record it in an appropriate structure: MadWeightAnalyzer MadWeightDiagramClass
     105   * determine the class of ECS as well as the block structure of each blob, and organize the blocks in the adequate order
     106   * write the main code in fortran (see below) that is the topology-dependent part of the phase-space generator
     107
     108=== The phase-space generator ===
     109
     110The associated fortran code will partially result from the supra-code. It will contain several things: PhaseSpaceGenerator
     111
     112==== I. Generic subroutines ====
     113
     114One subroutine for each block and each class of ECS. As discussed in the paper, a block (or a class of ECS) is nothing else than a specific local change of phase-space parametrization. Given the inputs (parameters+dynamical variables), the subroutine associated to a block will solve the change of variables and compute the jacobian. These subroutines are generic so that there is no need to generate them with a supra code.
     115
     116==== II. The main code ====
     117
     118This 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:
     119
     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.
     121
     122   * 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'''
     123
     124   * generation/resolution of the kinematics associated to the blobs (+ computation of the associated jacobian). These are made of blocks, in the sense discussed in the paper. The kinematics of any blob will be generated/resolved by calling (a certain number of) block subroutines defined in I.
     125
     126   * resolution of the change of variables corresponding to the ECS at work. This will be done by the call of a generic subroutine (see I.)
     127
     128So the main code will be made of external calls to several blocks and one class of ECS. The number of calls, their ordering, the inputs that are passed through the subroutines depend on the topology. These issues will be handled in the supra-code.
     129
     130==== III. Status and checks ====
     131
     132Each block or ECS has been checked at least once by computing the phase space volume with a phase-space parametrization involving this block.
     133
     134||  '''Block_name''' || '''Toplology'''  || '''PS volume'''  ||  '''check''' ||
     135|| Block A  ||  ECS B + bl A (massless)  ||  3.78 +/- 0.008 || 3.89  ||
     136|| Block B  || ECS B + bl B (massless)  || 0.0166 +/- 0.0001  ||  0.0166  ||
     137|| Block B  || ECS B + bl B (massive)  || 0.001407 +/- 9e-6  ||  0.001401  ||
     138|| Block C  || ECS B + bl C (massless)  || (6.28+/-0.014)e-5  || 6.30e-5  ||
     139|| Block C  || ECS B + bl C (massive)  || (9.75+/-0.02)e-6  || 9.74e-6  ||
     140|| Block D  || ECS B + bl D (massless)  || (6.27+/-0.02)e-5  || 6.30e-5  ||
     141|| Block E  || ECS B + bl E (massless)  || 0.0165 +/-0.0001  ||  0.0166  ||
     142|| ECS A  || ECS A (3FP, massless)  || (6.301+/-0.01)e-5  || 6.299e-5  ||
     143|| ECS A  || ECS A (3FP, massive)  || (9.741+/-0.01)e-6  || 9.755e-6  ||
     144|| ECS B  || ECS B (3FP, massive)  || (9.740+/-?)e-6  || 9.755e-6  ||
     145|| ECS C  || ECS C (3FP, massless)  || (6.293+/-0.003)e-5  || 6.299e-5  ||
     146|| ECS C  || ECS C (3FP, massive)  || (1.716+/-0.001)e-5  || 1.715e-5  ||
     147|| ECS D  || ECS D (massless)  || 43.73+/-12.41  || 44.60  ||
     148|| ECS D  || ECS D ( massive)  || 692.29+/-17.88  || 694.28  ||
     149|| ECS E  || ECS E (massless)  || 0.01662+/-?  || 0.01662 ||
     150|| ECS E  || ECS E ( massive)  || 0.001627  || 0.001631  ||
     151|| ECS F2  || ECS E (massless)  || 0.01664+/-2e-5  || 0.01662  ||
     152|| ECS F2  || ECS E (massive)  || 0.001631+/-3e-6  || 0.001631  ||
     153
     154==== IV. Differential weight ====
     155
     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.
     157
     158-
     159
     160--- MadWeight: the general code
     161
     162==== Status and checks ====
     163
     164Checked processes (after the fuse of MW in MG trunk):
     165
     166||  '''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 ||||
     177
     178==== Things to do for the next Release ====
     179   1. Grid initialization.
     180   1. relative weight between solution in Blob/ECS
     181
     182==== 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.
     189
     190-- Main.OlivierMattelaer - 24 Nov 2008
     191
     192
     193