This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
mem:mem [2013/07/05 14:28] – decosa | mem:mem [2013/07/05 15:21] (current) – decosa | ||
---|---|---|---|
Line 4: | Line 4: | ||
This is the development branch, the most up to date. | 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. | 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. | 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 a **set of datacards** and an **LHCO file**. | + | </ |
- | 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: | + | |
+ | 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 34: | 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 43: | 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/ | ||
+ |