This shows you the differences between two versions of the page.
| Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
| mem:mem [2013/07/05 12:55] – decosa | mem:mem [2013/07/05 15:21] (current) – decosa | ||
|---|---|---|---|
| Line 1: | Line 1: | ||
| | | ||
| - | In order to have MadWeight on the machine just check out it via bazaar (bazaar is installed on our tier 3): | + | In order to have MadWeight on the machine just check out it via bazaar (bazaar is installed on our tier 3). |
| - | * '' | + | |
| This is the development branch, the most up to date. | This is the development branch, the most up to date. | ||
| + | < | ||
| + | bzr branch lp: | ||
| + | </ | ||
| - | Madweight requires as input a **set of datacards** and an **LHCO file**. | + | Once MadWeight is on the working machine, the first step is to calculate the matrix elements for the process to consider: the hypothesis HP. |
| + | It is possible to it by launching MadGraph and defining the hypothesis process. | ||
| + | |||
| + | < | ||
| + | cd madweight_mc_perm | ||
| + | ./ | ||
| + | generate p p > t h b~ j $$ w+ w-, h > b b~ , t > b l+ vl # generating the process | ||
| + | add process p p > t~ h b j $$ w+ w-, h > b b~ , t~ > b~ l- vl~ # adding the complementary channel | ||
| + | output madweight tHbb_HP | ||
| + | </ | ||
| + | |||
| + | |||
| + | Madweight requires as input an **LHCO file** and a **set of datacards**. | ||
| + | |||
| + | |||
| + | |||
| + | ===== LHCO files ===== | ||
| A **LHCO** file is a text file collecting basic info on all the events of the sample. | A **LHCO** file is a text file collecting basic info on all the events of the sample. | ||
| Line 13: | Line 31: | ||
| __**LHE to LHCO**__: lhco files can be converted from LHE files. To do this a script is available from MadGraph website. Download it in the working directory: | __**LHE to LHCO**__: lhco files can be converted from LHE files. To do this a script is available from MadGraph website. Download it in the working directory: | ||
| - | * '' | + | |
| - | * '' | + | < |
| + | wget https:// | ||
| + | mv lhe2lhco.py.txt lhe2lhco.py | ||
| + | </ | ||
| __**EDM LHE to LHCO**__: lhco files can converted also from EDM LHE files using a standard tool provided within the CMS framework (ExternalLHEAsciiDumper). Copy it and change the name of the input file: | __**EDM LHE to LHCO**__: lhco files can converted also from EDM LHE files using a standard tool provided within the CMS framework (ExternalLHEAsciiDumper). Copy it and change the name of the input file: | ||
| Line 22: | Line 43: | ||
| __**ROOT to LHCO**__: To convert root files to LHCO files a standard tool does not exist, because it depends on the structure of the user root file. A root/lhco converter script to convert tH analysis root ntuples in lhco file is going to be implemented and it will be posted here as soon as it will be ready. | __**ROOT to LHCO**__: To convert root files to LHCO files a standard tool does not exist, because it depends on the structure of the user root file. A root/lhco converter script to convert tH analysis root ntuples in lhco file is going to be implemented and it will be posted here as soon as it will be ready. | ||
| + | |||
| + | ===== Datacards ===== | ||
| + | |||
| + | Datacards are stored in the sundirectory Cards of the created folder tHbb_HP. | ||
| + | The datacards to customize are three: param_card.dat, | ||
| + | |||
| + | * __**param_card.dat**__: | ||
| + | * __**run_card.dat**__: | ||
| + | * __**Madweight_card.dat**__: | ||
| + | < | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | | ||
| + | </ | ||
| + | |||
| + | |||
| + | |||
| + | |||
| + | ===== Producing weights ===== | ||
| + | |||
| + | In order to apply the matrix element method to a sample S under the hypothesis HP, Madweight has to be run using the LHCO file relative to the sample S and the datacards related to the hypothesis HP. | ||
| + | * Copy the .lhco file in the Events folder with name '' | ||
| + | < | ||
| + | cd tHbb_HP | ||
| + | cp $MYSAMPLE_PATH/ | ||
| + | </ | ||
| + | * Set the transfer function (TF), choosing it among those available: | ||
| + | < | ||
| + | ./ | ||
| + | define_transfer_fct | ||
| + | Please choose your transfer_function between | ||
| + | 0 / dbl_gauss_pt_jet | ||
| + | 1 / uniform | ||
| + | 2 / jet_gaussian | ||
| + | 3 / user | ||
| + | 4 / single_gaussian | ||
| + | 5 / gauss_on_leptons | ||
| + | 6 / all_delta | ||
| + | | ||
| + | </ | ||
| + | |||
| + | * It is possible to define a custom transfer function. Instructions will come soon. | ||
| + | * Now that everything has been set launch MadWeight to compute the probabilities for each events: | ||
| + | < | ||
| + | ./ | ||
| + | </ | ||
| + | At this point MadWeight propose to check the datacards, change again transfer function if desired or specify another lhco file. If everything has been already set properly before, then just type 0 to go ahead. | ||
| + | |||
| + | |||
| + | In order to apply the MEM over different samples, but under the same hypothesis, follow the next steps: | ||
| + | |||
| + | * run_card.dat and param_card.dat have to remain unchanged once the hypothesis is defined. | ||
| + | * The parameter '' | ||
| + | * Also the lhco file has to be updated each time MadWeight has to be run over a new sample: the lhco file related to the specific sample must be placed in the Events directory and named '' | ||
| + | |||
| + | |||
| + | |||
| + | ===== Results ===== | ||
| + | |||
| + | Results will be collecting in the file weights.out in the Events/ | ||