1 |
claudioc |
1.1 |
How to make limit plots
|
2 |
|
|
------------------------
|
3 |
|
|
|
4 |
claudioc |
1.2 |
(1) Make sure you have all the necessary L1 files. You need
|
5 |
|
|
*ALL* of them unless you know what you are doing.
|
6 |
claudioc |
1.1 |
|
7 |
|
|
(2) Merge the L1 files into L2 files using the script
|
8 |
claudioc |
1.2 |
from Marc
|
9 |
|
|
UserCode/claudioc/ScansICHEP2012/scripts/merge.py
|
10 |
claudioc |
1.1 |
inside the directory that contains the L1 files
|
11 |
|
|
|
12 |
|
|
(3) The L2 files coming out of Marc's script need to be sorted, eg,
|
13 |
|
|
sort L2_file_from_Marc.txt >> L2_file_from_Marc_sorted.txt
|
14 |
|
|
|
15 |
claudioc |
1.2 |
(4) Make sure to have the "right" lookup tables from Nik
|
16 |
|
|
- LandsOutput_Aprioiri.txt
|
17 |
|
|
- LandsOutput_Posterior.txt
|
18 |
claudioc |
1.1 |
|
19 |
|
|
(5) Edit makeNtuple.C and change the names of the L2 files
|
20 |
|
|
to be whatever you want
|
21 |
|
|
|
22 |
|
|
(6) The name of the xsection file is also hardwired in makeNtuple.C
|
23 |
|
|
(It defaults to the gluino xsection, if you are doing sbottom
|
24 |
|
|
pairs, change the name of the file!)
|
25 |
|
|
|
26 |
|
|
(7) Dump the L2 files into a flat ntuple
|
27 |
|
|
|
28 |
|
|
unix> root -l
|
29 |
|
|
root> .L Functions.cc
|
30 |
|
|
root> .L makeNtuple()
|
31 |
|
|
root> makeNtuple()
|
32 |
|
|
|
33 |
|
|
The output is ntuple.root
|
34 |
|
|
|
35 |
|
|
(8) To make plots
|
36 |
|
|
|
37 |
|
|
unix> root -l ntuple.root
|
38 |
|
|
root> .L makePlots_T1tttt.C
|
39 |
|
|
root> makePlots()
|
40 |
|
|
|
41 |
|
|
The histogram binning etc is hardwired in the C file,
|
42 |
|
|
you need to change it appropriately for a different model.
|
43 |
|
|
MAKE SURE THAT THE BINNING CORRESPONDS TO THE GRANULARITY
|
44 |
|
|
OF THE L2 FILES. ALSO MAKE SURE THAT THE BINS ARE PROPERLY
|
45 |
claudioc |
1.2 |
CENTERED. |