Changes between Initial Version and Version 1 of DifferentialWeight


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

--

Legend:

Unmodified
Added
Removed
Modified
  • DifferentialWeight

    v1 v1  
     1
     2
     3=== Definition ===
     4
     5=== Practical Computations ===
     6
     7=== How to run ===
     8
     9This is module is include in MW version XXXX
     10
     11In 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
     16You 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
     40Usual 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
     42Additionally, 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
     44In 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
     48requirement: python (at least 2.5) with numpy installed
     49
     50program: [attachment:correlation.py.txt correlation.py.txt]: correlation.py.txt (note that you have to remove the .txt)
     51
     52after the ''''''main'''''' == '''''name''''':
     53
     54edit the init_dir (position of data)
     55
     56edit the output_dir (where to write the output) then write the following line:
     57
     58MW_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
     66requirement: python (at least 2.5) with numpy installed
     67
     68program: [attachment:dist.py.txt dist.py.txt]: dist.py.txt (note that you have to remove the .txt)
     69
     70all modification to perform, set the end of the file
     71
     72&lt;verbatim&gt;
     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&lt;/verbatim&gt;</div>
     75
     76'MW_mean','MW_error','MW_correlation' are the value computed for one event.
     77
     78Nb_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