| 1 |
peller |
1.1 |
VHbb Analysis
|
| 2 |
|
|
Philipp Eller, 01.10.2012
|
| 3 |
|
|
|
| 4 |
|
|
|
| 5 |
|
|
EXPLANATION:
|
| 6 |
|
|
___________
|
| 7 |
|
|
|
| 8 |
|
|
The Analysis is split up in several steps. Having the step2 Ntupels the code does the following:
|
| 9 |
|
|
|
| 10 |
|
|
- prepare the environement (skim trees, write sample.info): This uses the info about your MC samples and data stored in *samples*.cfg and creates a new subdirectory 'env/'
|
| 11 |
|
|
- write systematics and the regression: This creates a new subdirectory 'sys/'
|
| 12 |
|
|
(optional: - do a new training)
|
| 13 |
|
|
- apply all MVAs to the trees. this creates the subdirectory 'MVAout/'
|
| 14 |
|
|
- write the RooFit Workspaces and the DataCards
|
| 15 |
|
|
|
| 16 |
|
|
|
| 17 |
|
|
EXAMPLE OF STEPS TO RUN:
|
| 18 |
|
|
_______________________
|
| 19 |
|
|
*** the analysis is designed to be run on the T3 cluster, you can run locally by not using the 'submitThem.py' script
|
| 20 |
|
|
*** make sure you understand the config*TeV, where everything can be set for the Analysis. If something needs to be changed it should be done via this config, NOT by changing the code!
|
| 21 |
|
|
|
| 22 |
|
|
- Edit pathConfig*TeV and make sure the directories for plots and Limits exist, point the right path to the code and your directory where you have the trees)
|
| 23 |
|
|
- Edit *samples*.cfg according to your MC samples and Data
|
| 24 |
|
|
(if you need to change something later there is an updater script updateinfo.py)
|
| 25 |
|
|
|
| 26 |
|
|
- run ./submitThem.py -T 8TeV -J prep
|
| 27 |
|
|
|
| 28 |
|
|
check with ./showinfo.py PATH/env/ that everything is correct
|
| 29 |
|
|
|
| 30 |
|
|
- run ./submitThem.py -T 8TeV -J sys
|
| 31 |
|
|
- run ./submitThem.py -T 8TeV -J eval
|
| 32 |
|
|
|
| 33 |
|
|
(if you want to test or run individual samples, add the option -S SAMPLE, e.g. '-S ZH125')
|
| 34 |
|
|
|
| 35 |
|
|
make some plots (at the latest here) with
|
| 36 |
|
|
- run ./submitThem.py -T 8TeV -J plot -M 125 (makes plot with ZH125 as Signal)
|
| 37 |
|
|
|
| 38 |
|
|
make the datacards:
|
| 39 |
|
|
-run ./submitThem.py -T 8TeV -J dc -M 125 (or -M all for all Masses)
|
| 40 |
|
|
|
| 41 |
|
|
you should now have DC, TH and WS files for every Mass and Region in the Limit Directory.
|
| 42 |
|
|
|
| 43 |
|
|
From there, using CMSSW:
|
| 44 |
|
|
|
| 45 |
|
|
to combine channels
|
| 46 |
|
|
- combineCards.py vhbb_DC_*1.txt vhbb_DC_2*.txt ... >> vhbb_combinedDC.txt
|
| 47 |
|
|
|
| 48 |
|
|
then calculate the Limit, for example Asymptotic
|
| 49 |
|
|
- combine -M Asymptotic vhbb_combinedDC.txt
|
| 50 |
|
|
|
| 51 |
|
|
or with Toys the option -t # (# number of toys)
|
