| 1 | |
| 2 | |
| 3 | === Definition === |
| 4 | |
| 5 | === Practical Computations === |
| 6 | |
| 7 | === How to run === |
| 8 | |
| 9 | This is module is include in MW version XXXX |
| 10 | |
| 11 | In order to use this module You need to modify three files |
| 12 | 1. dbook.inc: changes the definitions accordingly of what you want to do.<br />[attachment:dbook.inc dbook.inc]: Example of dbook.inc in order to plot three different variables. |
| 13 | 1. topgraph.f: defines the value to plot<br />[attachment:topgraph.f topgraph.f]: Example of topgraph.f in order to have Mtt, Pt of the top quark and angular distribution of the top quark. |
| 14 | 1. Madweight_card.dat: |
| 15 | |
| 16 | You should add the following option in the MadWeight _card |
| 17 | || BLOCK || variable name || example of value || explanation || |
| 18 | || mw_run || histo || T || # creates histogram for value defined in topgraph.f || |
| 19 | || mw_run || histo_by_events || T || # creates histogram for each events (Optional) False by default || |
| 20 | |
| 21 | === Kinematical fit === |
| 22 | |
| 23 | ==== What is done ==== |
| 24 | * matrix element to one |
| 25 | * mass fix at pole mass |
| 26 | * integration on the transfer function |
| 27 | |
| 28 | ==== How to run ==== |
| 29 | * use the KIN_TEMP in ingrid: /home/fynu/omatt/scratch/MG_ME_DMTT/ |
| 30 | * copy that Template |
| 31 | * ./bin/PassToMadWeight |
| 32 | * configure the dbook.inc/topgraph.f (see above) |
| 33 | * configure all the card |
| 34 | * launch |
| 35 | |
| 36 | === Estimating error and deviation between theory experimental result === |
| 37 | |
| 38 | ==== The principle ==== |
| 39 | |
| 40 | Usual estimation of the error on each bin based on the square root of the number of entry is not valid in the current scope due to the strong correlation between the bin entry. We then need a more evolve way to estimate the standard deviation. |
| 41 | |
| 42 | Additionally, we want to be able to claim the validity of an experimental curves compare to it's theoretical couterpart. Due to the correlation, the usual chi-square method is not valid and we need to use an extended version based on the correlation matrix. |
| 43 | |
| 44 | In consequence, in this section we will present two tools. The first one computes the standard deviation and the correlation matrix. it's base on the computation of the expected value/standard deviation/correlation of the value of one bin for an experiment containing a random number of experimental event following of poisson of parameter one. (this method was validate using some simple case) |
| 45 | |
| 46 | ==== Compute the error/sigma/correlation ==== |
| 47 | |
| 48 | requirement: python (at least 2.5) with numpy installed |
| 49 | |
| 50 | program: [attachment:correlation.py.txt correlation.py.txt]: correlation.py.txt (note that you have to remove the .txt) |
| 51 | |
| 52 | after the ''''''main'''''' == '''''name''''': |
| 53 | |
| 54 | edit the init_dir (position of data) |
| 55 | |
| 56 | edit the output_dir (where to write the output) then write the following line: |
| 57 | |
| 58 | MW_correlation_mult(XXX, YYY, bin(VVV, WWW)) where |
| 59 | 1. XXX is the number of events in your directory |
| 60 | 1. YYY is the variable number |
| 61 | 1. VVV is the value of the first bin |
| 62 | 1. WWW is the width of a bin (we suppose that we have 50 bins) then python correlation.py |
| 63 | |
| 64 | ==== Estimate the chi-square ==== |
| 65 | |
| 66 | requirement: python (at least 2.5) with numpy installed |
| 67 | |
| 68 | program: [attachment:dist.py.txt dist.py.txt]: dist.py.txt (note that you have to remove the .txt) |
| 69 | |
| 70 | all modification to perform, set the end of the file |
| 71 | |
| 72 | <verbatim> |
| 73 | |
| 74 | <div id="''mcePaste">if ''''''main''''''=='''''name''''':</div> <div id="''mcePaste"> obj2=Chi_carre('MW_mean','MW_error','MW_correlation',Nb_event)</div> <div id="''mcePaste"> chi''carre=obj2.compute_from_file('MW_SM_event',mode=5)</div> <div id="''mcePaste"> print chi''carre</verbatim></div> |
| 75 | |
| 76 | 'MW_mean','MW_error','MW_correlation' are the value computed for one event. |
| 77 | |
| 78 | Nb_event is the number of event if your analyzed histogram. |
| 79 | |
| 80 | 'MW_SM_event' is the histogram files. |
| 81 | |
| 82 | -- Main.OlivierMattelaer - 13 Oct 2009 |
| 83 | |
| 84 | |
| 85 | |
| 86 | |
| 87 | |
| 88 | |