Changes between Version 16 and Version 17 of Reweight


Ignore:
Timestamp:
Feb 1, 2016, 11:32:55 AM (9 years ago)
Author:
Olivier Mattelaer
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • Reweight

    v16 v17  
    88= Reweighting of NLO sample
    99
    10  As in 2.3.3, we only offer "LO" reweighting of NLO sample (also call "kamikaze" reweighting). This consists to reweight the S events by the born matrix element and the H events by the real matrix element. We have shown that for EW EFT theory this is a correct approach. The way to run the code is the same as for the LO (see below)
     10 As in 2.3.3, we only offer Leading Order type reweighting of NLO sample (also call "kamikaze" reweighting). This consists to reweight the S events by the born matrix element and the H events by the real matrix element. We have shown that for EW EFT theory this is a correct approach. The way to run the code is the same as for the LO (see below)
    1111 
    1212 The correct NLO reweighting is only available on alpha version so far and can be requested by email to olivier Mattelaer. This version will be available in 2.4.0
    … …  
    1414 In order to perform the correct NLO reweighting, you need to generate the original samples with the following options:
    1515  True = keep_rwgt_info ! keep the information for reweighting in the lhe file (huge impact on size file)
    16  After the reweighting will automatically be the NLO one and everything follows the LO implementation.
    17  In order to run the "kamikaze" reweighting on this sample, you can include in the reweighting_card the following command: "change mode LO"
     16 After the reweighting will automatically be the NLO and, for the user point of view, the interface to run the re-weighting is the same as the LO one described below.
     17 In order to run the "kamikaze" reweighting on this sample, you can include in the reweighting_card the following command (at the beginning of the file):
     18"""
     19    change mode LO
     20"""
    1821
    1922
    2023= Technical details
    2124 == Limitation
    22   1. The two theoretical hypothesis should differ ONLY by difference in the model parameters. So this should be the same pdf/cut (you can change model and process definition).
    23   2. The scale choice for the reweighting is done events by events and is given by the value in the LHE file. The same scale is used for both Hypothesis.
    24   3. In presence of decay chain the order of the particle in the events file is important. This is important if you want to use this tools with LHE events not produced by MadGraph/MadEvent.
    25 
    26  == How to get the code
    27  The code is released with any Madgraph version (since 2.0)
     25  1.
     26
     27
     28The two theoretical hypothesis should differ ONLY by difference in the model parameters. So this should be the same pdf/cut (you can change model and process definition).
     29  2. Note that the functional form of the hard scale can not be changed.
     30  3. In presence of decay chain the order of the particle in the events file is important. This is important if you want to use this tools with LHE events not produced by MadGraph5_aMC@NLO.
    2831
    2932 == Do I need special package
    30  Starting from MG5.2.3.2, you need to have f2py on your machine. The easiest way to install f2py is to install numpy (if not already done).
     33 Starting from MadGraph5_aMC@NLO.2.3.2, you need to have f2py on your machine. The easiest way to install f2py is to install numpy (if not already done).
    3134 
    3235 == How to use the code on the flight.
     36
     37--> Add comment: Same for the NLO. Here some examples.
    3338
    3439 When running the generation of events (./bin/generate_events from the local directory or "launch" via the mg5 interface)
    … …  
    7782 The default file is basically empty and will make the re-weighting to crash since both theoretical hypothesis are identical.
    7883
    79  == How to use the code afterwards.
     84 == How to use the code after the generation of events as been completed.
     85
     86--> Add for the NLO
    8087
    8188 You can also use the madevent interface as explained below:
    8289 
    8390 1. go to the process directory.
    84  2. launch the '''./bin/madevent''' script
     91 2. launch the '''./bin/madevent''' script  
    8592 3. type '''reweight RUN_NAME'''
    8693 4. then you will see the following question:
    … …  
    9299 [0, done, 1, reweight, enter path][60s to answer]
    93100}}}
    94     If you didn't define the content of this file already. You need to do it. The syntax is explained inside the file, and you can see example below (validation section). The '''important''' point, is that the first line should be launch, and then you specify which parameter, you want to modify. This is the exact same syntax has for scripting a scan over parameter.
     101    If you didn't define the content of reweight_card.dat already. You need to do it. The syntax is explained inside the file, and you can see example below (validation section). The '''important''' point, is that the first line should be launch, and then you specify which parameter, you want to modify. This is the exact same syntax has for scripting a scan over parameter.
    95102 5. exit the file and you are done (the script will run).
    96103
    … …  
    102109
    103110 = Options
    104   Note all the following options will be available in MG5 version 2.3.2
     111  Note all the following options will be available in Madgraph5_aMC@NLO since version 2.3.2
    105112   1. "change model <XXX>" performed the reweighting within a new model (you then need to profide a full param_card and not a difference)
    106113   2. "change process <XXX>" change the process definition of the process.
    … …  
    109116   5. "change helicity <True|False>": perform reweighting for the given helicity (True --default--) or do the sum over helicity.
    110117   6. "change rwgt_dir <PATH>": change directory where the computation is performed. This can be use to avoid to recreate/recompile the fortran executable if pointing to a previously existing directory.
     118   7. change mode LO:
    111119
    112120
    113121 = Input/Output format
    114122
    115  1. the output format follows the Leshouches agreement (not yet published): http://phystev.in2p3.fr/wiki/2013:groups:tools_lheextension#the_proposal
     123 1. the output format follows the Leshouches agreement version 3 [link arxiv]
    116124
    117125as an example the header looks like:
    … …  
    181189}}}
    182190
    183 4. An example of the full log is:
    184 {{{
    185 ************************************************************
    186 *                                                          *
    187 *           W E L C O M E  to  M A D G R A P H  5          *
    188 *                      M A D E V E N T                     *
    189 *                                                          *
    190 *                 *                       *                *
    191 *                   *        * *        *                  *
    192 *                     * * * * 5 * * * *                    *
    193 *                   *        * *        *                  *
    194 *                 *                       *                *
    195 *                                                          *
    196 *         VERSION 2.0.0.beta3           2013-02-14         *
    197 *                                                          *
    198 *    The MadGraph Development Team - Please visit us at    *
    199 *    https://server06.fynu.ucl.ac.be/projects/madgraph     *
    200 *                                                          *
    201 *               Type 'help' for in-line help.              *
    202 *                                                          *
    203 ************************************************************
    204 INFO: load configuration from /Users/omatt/.mg5/mg5_configuration.txt 
    205 INFO: load configuration from /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Cards/me5_configuration.txt 
    206 INFO: load configuration from /Users/omatt/Documents/eclipse/reweight_module/input/mg5_configuration.txt 
    207 INFO: load configuration from /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Cards/me5_configuration.txt 
    208 MGME5>reweight run_01
    209 Do you want to edit one cards (press enter to bypass editing)?
    210   1 / reweight  : reweight_card.dat
    211  you can also
    212    - enter the path to a valid card.
    213  [0, done, 1, reweight, enter path][60s to answer]
    214 >1
    215 open /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Cards/reweight_card.dat
    216 Do you want to edit one cards (press enter to bypass editing)?
    217   1 / reweight  : reweight_card.dat
    218  you can also
    219    - enter the path to a valid card.
    220  [0, done, 1, reweight, enter path]
    221 >
    222 INFO: Extracting the banner ...
    223 INFO: process: p p > e+ e-
    224 INFO: options: 
    225 INFO: Running Reweight
    226 launch
    227 INFO: detected model: sm. Loading...
    228 INFO: generating the square matrix element for reweighting
    229 INFO: generate p p > e+ e- ;
    230 INFO: Done 1.197
    231 set aewm1 100
    232 INFO: modify param_card information BLOCK sminputs with id (1,) set to 100.0
    233 Do you want to edit one cards (press enter to bypass editing)?
    234   1 / param     : param_card.dat
    235   9 / plot      : plot_card.dat
    236  you can also
    237    - enter the path to a valid card or banner.
    238    - use the 'set' command to modify a parameter directly.
    239      The set option works only for param_card and run_card.
    240      Type 'help set' for more information on this command.
    241  [0, done, 1, param, 9, plot, enter path]
    242 The answer to the previous question is not set in your input file
    243 Use 0 value
    244 INFO: starts to compute weight for events with the following modification to the param_card:
    245 INFO: set param_card sminputs 1 100.0
    246  
    247 INFO: Event nb 0 current time: 23h39
    248 INFO: Event nb 1000  0.61s 
    249 INFO: Event nb 2000  1.1s 
    250 INFO: Event nb 3000  1.5s 
    251 INFO: Event nb 4000  2.1s 
    252 INFO: Event nb 5000  2.5s 
    253 INFO: Event nb 6000  2.9s 
    254 INFO: Event nb 7000  3.3s 
    255 INFO: Event nb 8000  3.8s 
    256 INFO: Event nb 9000  4.2s 
    257 INFO: All event done  (nb_event: 10000)  4.6s 
    258 INFO: Event /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Events/run_01/unweighted_events.lhe have now the additional weight
    259 INFO: new cross-section is : 1135.25 pb
    260 launch
    261 set aewm1 200
    262 INFO: modify param_card information BLOCK sminputs with id (1,) set to 200.0
    263 Do you want to edit one cards (press enter to bypass editing)?
    264   1 / param     : param_card.dat
    265   9 / plot      : plot_card.dat
    266  you can also
    267    - enter the path to a valid card or banner.
    268    - use the 'set' command to modify a parameter directly.
    269      The set option works only for param_card and run_card.
    270      Type 'help set' for more information on this command.
    271  [0, done, 1, param, 9, plot, enter path]
    272 The answer to the previous question is not set in your input file
    273 Use 0 value
    274 INFO: starts to compute weight for events with the following modification to the param_card:
    275 INFO: set param_card sminputs 1 200.0
    276  
    277 INFO: Event nb 0 current time: 23h39
    278 INFO: Event nb 1000  0.48s 
    279 INFO: Event nb 2000  1s 
    280 INFO: Event nb 3000  1.5s 
    281 INFO: Event nb 4000  2s 
    282 INFO: Event nb 5000  2.4s 
    283 INFO: Event nb 6000  2.9s 
    284 INFO: Event nb 7000  3.3s 
    285 INFO: Event nb 8000  3.8s 
    286 INFO: Event nb 9000  4.2s 
    287 INFO: All event done  (nb_event: 10000)  4.7s 
    288 INFO: Event /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Events/run_01/unweighted_events.lhe have now the additional weight
    289 INFO: new cross-section is : 1095.28 pb
    290 launch
    291 set aewm1 300
    292 INFO: modify param_card information BLOCK sminputs with id (1,) set to 300.0
    293 Do you want to edit one cards (press enter to bypass editing)?
    294   1 / param     : param_card.dat
    295   9 / plot      : plot_card.dat
    296  you can also
    297    - enter the path to a valid card or banner.
    298    - use the 'set' command to modify a parameter directly.
    299      The set option works only for param_card and run_card.
    300      Type 'help set' for more information on this command.
    301  [0, done, 1, param, 9, plot, enter path]
    302 The answer to the previous question is not set in your input file
    303 Use 0 value
    304 INFO: starts to compute weight for events with the following modification to the param_card:
    305 INFO: set param_card sminputs 1 300.0
    306  
    307 INFO: Event nb 0 current time: 23h39
    308 INFO: Event nb 1000  0.62s 
    309 INFO: Event nb 2000  1.1s 
    310 INFO: Event nb 3000  1.5s 
    311 INFO: Event nb 4000  2s 
    312 INFO: Event nb 5000  2.5s 
    313 INFO: Event nb 6000  2.9s 
    314 INFO: Event nb 7000  3.4s 
    315 INFO: Event nb 8000  3.8s 
    316 INFO: Event nb 9000  4.3s 
    317 INFO: All event done  (nb_event: 10000)  4.7s 
    318 INFO: Event /Users/omatt/Documents/eclipse/reweight_module/PROC_sm_11/Events/run_01/unweighted_events.lhe have now the additional weight
    319 INFO: new cross-section is : 1329.52 pb
    320 INFO: Original cross-section: 949.58665707 +- 2.2589158361 pb
    321 INFO: Computed cross-section:
    322 INFO: 1 : 1135.24526333
    323 INFO: 2 : 1095.27584154
    324 INFO: 3 : 1329.51775579
    325 }}}
    326191 = Validation
    327192