Changes between Initial Version and Version 1 of MadWeightAdvancedOpttion


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

--

Legend:

Unmodified
Added
Removed
Modified
  • MadWeightAdvancedOpttion

    v1 v1  
     1
     2
     3
     4
     5Note that all available option (and default value) are defined in Source/MadWeight_File/Python/MW_param_default.inc Please read that files in order to know all the possible options. The main options are explained here. But not all of them.
     6
     7==== Block Run ====
     8
     9Define the cluster and some constraint on the run.
     10||  '''Options for block MW_run''' |||| '''option for job schedulling (step=2,5,6)''' ||
     11||  '''tag''' |||| '''Effect'''  ||
     12||  1 || || cluster choice (0: single machine// 1: condor // 2: SGE)  ||
     13|| || 11  || condor requirements (only for condor) ||
     14||  2 || || number of data to analyse ||
     15|| || 21  || avoid to consider the X first event of the file input.lhco ||
     16|| || 22  || if True (T) add new data directory to those already existing -use only for create new dir then launch job with relaunch/refine- (since 2.1.12) ||
     17||  3 || || control efficiency block by block ||
     18||  4 || || normalizes Weight by cross-section ||
     19||  5 || || number of points in MadGraph integral ||
     20||  6 || || number of points in MadWeight integral ||
     21||  7 || || [0.01] requested final precision for each weight  ||
     22|| || 710  || [0.5] minimal requested precision for the first refinment criteria (see below)  ||
     23|| || 711  || [0.02] first refinment criteria (see below)  ||
     24|| || 720  || [0.1] requested precision for the second refinment criteria (see below)  ||
     25|| || 721  || [2d-3] second refinment criteria  ||
     26|| || 731  || [2] number of iteration in permutation selection -first loop-  ||
     27|| || 732  || [2] number of iteration in configuration search (not used anymore -since MW2.4-)  ||
     28|| || 733  || [5] number of iteration in permutation selection -second loop-  ||
     29|| || 734  || [1] number of iteration in plotting diagram mode (not used for the moment)  ||
     30|| || 735  || [10] maximal number of iteration in final loop  ||
     31|| || 801  || [0] stop the control step when the number of idle job is strictly lower than ... ||
     32|| || 802  || [0] stop the control step when the total number of idle and running job is lower or equal to ... ||
     33||  9 || || [F] use the cut defined in run_card.dat for weight computation  ||
     34|| || 91  || [F] use the cut defined in run_card.dat for cross-section computation ||
     35||  weight_format |||| ['%(card)s.%(event)\t%(value)s\t%(error)s'] string formatting the weighted output (python string formatting--see below--) ||
     36||  cross_format |||| ['%(card)s\t%(value)s\t%(error)s'] string formatting the cross-section output (python string formating--see below--) ||
     37
     38In these Block you can define some (trully) '''advanced user''' option
     39
     40'''How are computed the weights''' If you want to change the defaults for the integration treatment, this explains how the computation is performed. The integral to compute is in a fact cut in different pieces computed one after the other (permutation and different channel of integration) The computation is done in the following way
     41   1. Launching a survey of 2 iterations with '''MW_RUN_6''' /2[=5000] points in each integrals.
     42   1. Check the survey, if the max of the survey has a precision less than '''MW_RUN_710''' [50%] then going back to step one but with 5 times more points launch for each integrals.
     43   1. We organize the contribution in the decreasing contribution order.
     44   1. We refine the contribution in this order, up to the points where the contribution are pointless: when the contribution is smaller than the error on the most important contribution time '''MW_RUN_711''' [0.01].
     45      1. This refinement is stopped when the error is smaller than the error of the main contribution time '''MW_RUN_721''' [0.1]
     46      1. This refinement is done in '''MW_RUN_735''' iterations, on the '''MW_RUN_733''' first iteration the grid is adapted and we launch '''MW_RUN_6''' points by iteration. Reducing '''MW_RUN_733''' can solve some problems of grid instabilities
     47      1. For the following iteration, the grid is not refined anymore and the number of points is computed in order to achieve the requested precision. In consequence changing '''MW_RUN_735''' doesn't influences the program except if you reach the maximal number of authorized points by iteration
     48
     49Concerning the output formatting: The syntax is available here: [http://docs.python.org/release/2.5.2/lib/typesseq-strings.html http://docs.python.org/release/2.5.2/lib/typesseq-strings.html]
     50
     51The available information are
     52   1. card: The card number
     53   1. number: The number associated to the MW counting (number of events selected) --not available for cross-section--
     54   1. lhco_number: number of the event in the LHCO file --not available for cross-section--
     55   1. trigger: number of the trigger in the LHCO file --not available for cross-section--
     56   1. value: result
     57   1. error: integration error on the result
     58
     59==== Block Perm ====
     60
     61Define how to treat the permutation in MadWeight.
     62
     63||  '''Option for block MW_perm''' |||| '''option for treating permutation (step=6)''' ||
     64||  '''tag''' |||| '''Effect'''  ||
     65||  1 || || makes permutation ||
     66||  2 || || [T] permutates jet with bjet for weight computation  ||
     67|| || 21 || [follows rule 2] don't consider b-tagging for the selection of events (change since 2.2.1)  ||
     68|| 3 || || authorize one ISR jet (not implemented yet) ||
     69==== Block Parameter ====
     70
     71MadWeight can run on different parameter value on the same run. For that you have three different choice (differentiate by the value of the entry with tag 1)
     72   * value '''0''': use the param_card defined by the user. Those one must be enter in the Card directory with name: param_card_1.dat, param_card_2.dat,...
     73   * value '''1''': Apply from the param_card.dat the modification define in MadWeight_card.dat (see below)
     74   * value '''2''': Same than one but all parameter modification are correlated (see below)
     75
     76||  '''Option for block MW_parameter''' ||||
     77||  '''tag''' || '''Effect'''  ||
     78||  1 || how create all the param card (see above for value description)  ||
     79||  2 || [F] cumulative mode add those card to the previous one (and create new event directory for those new card if step 5 already runs)  ||
     80||  n*10+1 || name of the block to modify for variable n (branching ratio is special see below) ||
     81||  n*10+2 || tag of the value to modify (if they are a double tag use two times this line (in the correct order) ||
     82||  n*10+3 || different value to insert (so you can use this tag more than once)  ||
     83
     84
     85==== Block Generation ====
     86
     87This block parametrize how MadWeightAnalyzer will find the approriate phase-space generator. More information on those option in MadWeightAnalyzer page.
     88
     89||  '''Option for block MW_Gen''' ||||
     90||  '''tag''' || '''Effect'''  ||
     91||  1 || chooses which propagator to allign in ECS: 0-> thinest/ 1-> chooses more local possibility  ||
     92||  2 || chooses which propagator to allign in blob: 0-> thinest/ 1-> chooses more local possibility  ||
     93||  3 || maximal multiplicity for ECS sector  ||
     94||  4 || uses blob solutions which optimizes Breit-Wigner integration  ||
     95||  5 || uses blob solutions which optimizes TransferFunction integration  ||
     96||  6 || use s balanced solutions (more than one true is possible)  ||
     97||  10 || allows ECS A (0 neut) for change of variables  ||
     98||  11 || allows ECS B (1 neut) for change of variables  ||
     99||  12 || allows ECS C (1 neut) for change of variables  ||
     100||  13 || allows ECS D (2 neut ttbar) for change of variables  ||
     101||  14 || allows ECS E (2 neut HWW) for change of variables  ||
     102||  15 || allows ECS F (2 neut gamma-gamma) for change of variables  ||
     103
     104==== Block !EventSelection ====
     105
     106This Block define is defined in version 2.1.13 for the first time. This block authorize some cuts selection on the particle present in the lhco file. With that blok you can easily constraint the the events used in MadWeight whitout using a complex selection. This blok has a special syntax like the following
     107{{{
     108Block EventSelection
     109#   TAG      VALUE                     UTILITY
     110     PartName_variable          min           #    first entry imposes a minimal value
     111     PartName_variable          max           #    second entries imposes a maximum value         
     112}}}
     113
     114!PartName must belong to
     115
     116   * begin (constraint on the first line of the lhco file)
     117   * jet
     118   * bjet (if btag used -MW_perm tag 2 and 21-)
     119   * electron
     120   * positron
     121   * photon
     122   * muon
     123   * amuon
     124   * tau (lhcotype=5)
     125   * atau
     126   * miss
     127   * init (lhcotype=7)
     128
     129!variable correspond to the different column of the lhco file and those name are the following one:
     130   * card
     131   * lhcoid
     132   * eta
     133   * phi
     134   * pt
     135   * jmass
     136   * ntrk
     137   * btag
     138   * hadem
     139   * dum1
     140   * dum2
     141
     142special cut on the energy is also availble
     143   * E
     144
     145Be carefull that you can change with that the definition between a jet/bjet, electron/positron,....
     146
     147-- Main.OlivierMattelaer - 02 Dec 2008
     148   * [attachment:generate_min_events.py.txt generate_min_events.py.txt]: Generate for a set of param_card
     149
     150