Changes between Initial Version and Version 1 of MadWeightFAQ

03/20/12 16:16:17 (8 years ago)



  • MadWeightFAQ

    v1 v1  
     3== MadWeight Frequently Asked Questions ==
     7=== General ===
     9==== Why did you choose to have only analytical change of variables? ====
     11Simply because it's much faster. If we computing time for the class d for analytical/numeric resolution of equation, we have a factor 2.
     13==== How do we inform that a 'non standard' particle is not visible in the detector. ====
     15If a particle is consider by default as visible in the detector, when in your case she is not (decaying in the detector in invisible particle exclusively).
     16you can add her PDG number to the invisible list present at the begining of this file: Source/Madeight_File/Python/
     19=== LHCO event file ===
     21==== Do I have to perform selection on the events? ====
     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-&gt;WW-&gt;(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>
     25MadWeight (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.
     27If you have three jet and want anyway to perform the Matrix Element Method (MEM) on this event, then yes you have to made a pre-selection.
     29==== I have two (or more) neutrino (invisible particle) in my events how can I disentangle the neutrinos implusion ====
     31That is not a problem at all. The information for the missing energy is not use from the events for MW integration. But this information can be use in the transfer functions. Note that all events in lhco files must have a missing ET line (line starting with 6). In fact in practise, MW compute the missing ET and disantagle the neutrinos impulsion for you (and this in fact point by point in the phase space)
     33==== Which information are mandatory/not use in the LHCO file? ====
     35'''MW''' needs two type of information for visible particle
     36   1. '''identification of the particle''' which mean the identification tag but also depending of the type the charge of the particle and/or the tagging
     37   1. '''kinematical information''' (that is obvious)
     39All the rest (in particular the missing ET line) can be put to zero. Note that for format compatibility you need to have those number written in the file (in particular the missing ET line)
     41=== Output file ===
     43==== How can I extract the individual event weights to create my own Likelihood function ? ====
     45You 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)
     47The 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}}}
     50'''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)
     52'''nb_event''' is the event number. (you can read this event in ./SubProcesses/MW_P_XXXXXXXXXXXXX/MY_RUN_NAME/MY_RUN_NAME_1_NB_EVENT/verif.lhco)
     54'''value''' is the (un)normalized weight
     56'''error''' is the error on this value
     58=== Possible causes of failure ===
     60==== failure of ./bin/madweight param (step: 1) ====
     61   * {{{
     62IOError: [Errno 2] No such file or directory: './Cards/ident_mw_card.dat'
     65This indicates that your transfer functions is not loaded properly. Did you perform the generation step $&gt;./bin/newprocess ? If you Downloaded your directory from the web did you assign your transfer functions with $&gt;./bin/ ?
     67   * 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==== failure of ./bin/madweight analyzer (step: 2) ====
     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$&gt; ./bin/newprocess<br /> Indeed the proc_card.dat must be the same between this moment and the time where you launch madweight analyzer.
     75   * {{{
     77                      you probably put bad option in your card
     80This error could arrive for two different reason
     81   1. You don't use the DecayChain format in your proc_card.dat (-&gt; change your proc_card.dat)
     82   1. You apply too strict restriction on available change of variable in MadWeight_card.dat
     84==== zero weight numbers ====
     85   * check the cuts in the run_card.dat. If cuts are applied on rapidity and you assume delta-like function for all visible angles, all phase-space points might be cut.
     87In this specific case, you can type "./bin/madweight relaunch" just after having suppressed the cuts in the run_card.dat. Note that you can bypass all the cut check via the MadWeight _Card.dat
     89   * An other current possibility is that your decay is forbidden (same coupling are at zero)
     91==== result files in directory Event are empty ====
     93   * make sure that you have collected the results:
     95$&gt;./bin/ collect plot
     97-- Main.OlivierMattelaer - 19 Mar 2009