Scripts for 2013 VH tau Moriond PAS
Install the HiggsToTauTau limit package,
cmsrel CMSSW_5_2_5
cd CMSSW_5_2_5/src/
cmsenv
cvs co -r V02-02-08 HiggsAnalysis/CombinedLimit
cvs co -r V00-03-10 HiggsAnalysis/HiggsToTauTau
python HiggsAnalysis/HiggsToTauTau/scripts/init.py --tag V00-01-00-noModels
# Remove any leftover shape files
rm HiggsAnalysis/HiggsToTauTau/setup/vhtt/*root
scram b -j 4; rehashthen check out this package:
git clone https://github.com/ekfriis/HIG-12-053.gitand get the latest VH datacards:
cd auxiliaries/datacards/collected/
cvs up -A -d vhttAll the tricks to build the results are contained in the Makefile. The important commands are:
cd HIG-12-053
# Run the post fit and make all the final mass distribution plots
make massplots
# Compute all the limits
make limits
# Plot the limits (they show up in limits/*pdf)
make plotlimits
# Make vh_table.tex for the PAS
make plots/vh_table.texor you can just run "make all", which runs everything.
Do the following to make the signal injected plots.
# Run a bunch of jobs with different pseudoexperiments
./inject_step1.sh
# Wait for those jobs to all finish, then collect all the results
./inject_step2.sh
# Plot the results
./inject_step3.shThis example is for the UW cluster. Omit the --condor option in the scripts to run on LXBatch.
https://twiki.cern.ch/twiki/bin/view/CMS/HiggsCombinationMoriond2013#Inputs
Checkout the vhtt area of the HCG SVN.
svn co svn+ssh://svn.cern.ch/reps/cmshcg/trunk/moriond2013/searches/vhtt hcg/searches/vhtt
svn co svn+ssh://svn.cern.ch/reps/cmshcg/trunk/moriond2013/couplings/vhtt hcg/couplings/vhttGenerate all intermediate masses required by the HCG group. Additionally, generate cards with the "Two Higgs format", where the SM Higgs 125 is considered as a background.
MASSLIST=coupling_masses.txt make limitdir
MASSLIST=searches_masses.txt SECONDHIGGS=1 make limitdirThe cards with second Higgs appear in auxiliaries_2ndhiggs/.../. Now run
./copy_to_hcg_svn.shto copy the appropriate cards into the HCG SVN area, then svn commit. You can double check that all the datacards are OK by running:
./double_check_svn_cards.sh