| 1 |
How to make limit plots
|
| 2 |
------------------------
|
| 3 |
|
| 4 |
(1) Make sure you have all the necessary L1 files. You need
|
| 5 |
*ALL* of them unless you know what you are doing.
|
| 6 |
|
| 7 |
(2) Merge the L1 files into L2 files using the script
|
| 8 |
from Marc
|
| 9 |
UserCode/claudioc/ScansICHEP2012/scripts/merge.py
|
| 10 |
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 |
(4) Make sure to have the "right" lookup tables from Nik
|
| 16 |
- LandsOutput_Aprioiri.txt
|
| 17 |
- LandsOutput_Posterior.txt
|
| 18 |
|
| 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 |
CENTERED. |
