This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
mem:mem [2013/07/05 11:51] – decosa | mem:mem [2013/07/05 15:21] (current) – decosa | ||
---|---|---|---|
Line 1: | Line 1: | ||
| | ||
- | Madweight requires as input a **set of datacards** and an **LHCO file**. | + | 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. | ||
+ | |||
+ | < | ||
+ | bzr branch lp: | ||
+ | </ | ||
+ | |||
+ | 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 7: | Line 30: | ||
* More info about the structure of LHCO files can be found at: http:// | * More info about the structure of LHCO files can be found at: http:// | ||
- | 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: |
- | * '' | + | |
- | * '' | + | |
- | 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: | + | < |
+ | 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: | ||
* [[http:// | * [[http:// | ||
- | 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**__: |
+ | |||
+ | ===== 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/ | ||
+ |