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.
Edit | Attach | Watch | Print version | History: r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r1 - 2015-09-10 - unknown
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright &© 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback