Yield output for search bins

Here the procedure fro obtaining event yields for Rcs/Kappa determination, yield prediction and, later on, limit calculation. NB! Following works with the branch: DESY-CMGTools-from-CMSSW_7_4_7_METcorr in DESY-CMS-SUS The development is ongoing in the folder CMGTools/TTHAnalysis/python/plotter/susy-1lep/RcsDevel. The macros available are:

• makeBinYields.py creates files with yield histograms according to the bin dictionary provided in searchBins.py. Usage:

       ./makeBinYields.py [options]
       

  where the options could be the same as in mcAnalysis.py plus the following:
  • -—mca : provide the mca file with the samples (like for mcAnalysis.py)
  • -—cut : provide the cut file (like for mcAnalysis.py)
  • -v (—verbose) : set verbosity level (0,1,2)
  • -—od : output directory for root files
  • -—grid : make a grid of yields for bkg/data: a TH2F with 2 bins on the y-axis: selected,anti-selected and 3 x-bins: muon,mu+ele,electrons yields.
  • -—signal : make a grid of yields for signal: mGluino vs mLSP with positive mGo for electrons and negative for muons
  • -—chunk : run only the n-th chunk (bin) of the bin dictionary
  • -b (--batch) : run all chunks on the batch system (NAF/DESY for now only) Example:

          ./makeBinYields.py --mca mca-SignalScan.txt --signal -v 1 -l 3 --od lumi3fb/scan -b #for signal
          ./makeBinYields.py --mca mca-Spring15.txt --grid -v 1 -l 3 --od lumi3fb/grid -b # for backgrounds/data
          

• mergeBins.py merges the bins according to the main/side-band and signal/control-region definitions. The output files are stored in the merged sub-folder of the input directory. The output files contain four directories: SR_MB, CR_MB and SR_SB, CR_SB. Usage:

        ./mergeBins.py in/put/dir/pattern
        

  where the input directory pattern should either be an existing directory ending with a / or, a pattern like input/dir/LT1_.

• appendRcsKappa.py appends Rcs and Kappa histograms for all samples in the merged root file: Rcs is obtained from dividing SR/CR, and Kappa comes from dividing the MB_Rcs with SB_Rcs. Usage:

        ./appendRcsKappa.py in/put/dir/pattern
        

  where the input directory pattern should either be an existing directory ending with a / or, a pattern like input/dir/LT1_. The input files should be products of the merge script.

• calcSyst.py creates _syst histograms from the up and down variation histograms stored for each sample. The maximum deviation is taken for each bin. Usage:

        ./calcSyst.py in/put/dir/pattern
        

  where the input directory pattern should either be an existing directory ending with a / or, a pattern like input/dir/LT1_. The input files should be products of the merge script.