Changes between Version 1 and Version 2 of MadWeightFAQ
- Timestamp:
- Apr 6, 2012, 4:33:02 PM (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
MadWeightFAQ
v1 v2 21 21 ==== Do I have to perform selection on the events? ==== 22 22 23 <em> <strong>Question from Andre Sznajder</strong><br /> The matrix element method uses just outgoing the hard scattering partons/leptons. As a concrete example let's take the Higgs to WW semi-leptonic mode H->WW->(mu nu)(qq'). As I am using CMS reconstructed event as an input, do I have to make a pre-selection of the reconstructed muons, jets and missing energy that corresponds to the parton level objects ? Or I can put in the LHCO all reconstructed jets, and leptons and it will somehow choose the correct ones ? </em>23 <em> **g>Question from Andre Sznajder**[[br]] The matrix element method uses just outgoing the hard scattering partons/leptons. As a concrete example let's take the Higgs to WW semi-leptonic mode H->WW->(mu nu)(qq'). As I am using CMS reconstructed event as an input, do I have to make a pre-selection of the reconstructed muons, jets and missing energy that corresponds to the parton level objects ? Or I can put in the LHCO all reconstructed jets, and leptons and it will somehow choose the correct ones ? </em> 24 24 25 25 MadWeight (MW) can only process event with the correct number of reconstructed events. So MW will first reject all the "wrong" events. So in practise, if you have two muons or three jet for an event, this one will not be accepted. … … 43 43 ==== How can I extract the individual event weights to create my own Likelihood function ? ==== 44 44 45 You will have files with that information: <br /> ./MYPROCESS/Events/MY_RUN_NAME/MY_RUN_NAME_weights.out<br /> ./MYPROCESS/Events/MY_RUN_NAME/MY_RUN_NAME_norm_weights.out<br />(by default MY_RUN_NAME is fermi)45 You will have files with that information:[[br]] ./MYPROCESS/Events/MY_RUN_NAME/MY_RUN_NAME_weights.out[[br]] ./MYPROCESS/Events/MY_RUN_NAME/MY_RUN_NAME_norm_weights.out[[br]] (by default MY_RUN_NAME is fermi) 46 46 47 47 The first file contains all the unormalized weight. The second want contains the weight normalized by the cross section. (Usually people are interrested in this one) The format for each file is the following: 48 {{{ nb_card.nb_event value error}}} 48 {{{ 49 nb_card.nb_event value error 50 }}} 49 51 50 52 '''nb_card''' is the card number (the XXX in Cards/param_card_XXX.dat). a correspondance with your modify parameter is show in ./MYPROCESS/Events/MY_RUN_NAME/likelihood_value.out or in ./Cards/mapping_card.dat) … … 63 65 }}} 64 66 65 This indicates that your transfer functions is not loaded properly. Did you perform the generation step $ >./bin/newprocess ? If you Downloaded your directory from the web did you assign your transfer functions with $>./bin/change_tf.py ?67 This indicates that your transfer functions is not loaded properly. Did you perform the generation step $>./bin/newprocess ? If you Downloaded your directory from the web did you assign your transfer functions with $>./bin/change_tf.py ? 66 68 67 69 * Did you enter an appropriate MadWeight_card.dat ? If you ask to change the value of a parameter which does not exist in the model you are working in, the code will not be able to produce the cards correctly. … … 69 71 ==== failure of ./bin/madweight analyzer (step: 2) ==== 70 72 71 * {{{ bad proc_card information: missing request process tag }}} Did you modify the proc_card.dat after the launch of ./bin/newprocess ? Then you have to type again 73 * {{{ 74 bad proc_card information: missing request process tag 75 }}} Did you modify the proc_card.dat after the launch of ./bin/newprocess ? Then you have to type again 72 76 73 $ > ./bin/newprocess<br />Indeed the proc_card.dat must be the same between this moment and the time where you launch madweight analyzer.77 $> ./bin/newprocess[[br]] Indeed the proc_card.dat must be the same between this moment and the time where you launch madweight analyzer. 74 78 75 79 * {{{ … … 79 83 80 84 This error could arrive for two different reason 81 1. You don't use the DecayChain format in your proc_card.dat (- >change your proc_card.dat)85 1. You don't use the DecayChain format in your proc_card.dat (-> change your proc_card.dat) 82 86 1. You apply too strict restriction on available change of variable in MadWeight_card.dat 83 87 … … 93 97 * make sure that you have collected the results: 94 98 95 $ >./bin/madweight.py collect plot99 $>./bin/madweight.py collect plot 96 100 97 101 -- Main.OlivierMattelaer - 19 Mar 2009
