HHggTauTau_flashggFinalFit/make_plots.sh at master · fsetti/HHggTauTau_flashggFinalFit · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
#!/usr/bin/env bash

set -x

source /cvmfs/cms.cern.ch/cmsset_default.sh
source /vols/grid/cms/setup.sh

tag=09May2023_nonRes_SM_bkgOnly

cmsenv
source setup.sh

nToys=500

make_toys(){
	pushd Plots
		rm -rf SplusBModels$tag
		python makeToys.py --inputWSFile ../Combine/Datacard_bkgOnly.root --ext $tag --dryRun --nToys $nToys
		#python makeToys.py --inputWSFile ../Combine/Datacard_ggtt_resBkg_syst.root --ext $tag --dryRun --nToys $nToys
		iter=0
		while [ $iter -lt $nToys ]
		do
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+1)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+2)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+3)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+4)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+5)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+6)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+7)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+8)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$(($iter+9)).sh	&
			./SplusBModels${tag}/toys/jobs/sub_toy_$iter.sh
			iter=$(($iter+10))
		done
	popd
}

make_SpB(){
	pushd Plots
		python makeSplusBModelPlot.py --inputWSFile ../Combine/Datacard_bkgOnly.root --cat "all" --doBands --ext $tag --parameterMap r:0
		#python makeSplusBModelPlot.py --inputWSFile ../Combine/Datacard_ggtt_resBkg_syst.root --cat "all" --doBands --ext $tag --doResonantBackground --unblind --parameterMap r:10

		mkdir -p /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag
		mkdir -p /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag/SplusBModels

		cp SplusBModels$tag/SR1_CMS_hgg_mass.png /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag/SplusBModels/SR1_CMS_hgg_mass_bkg.png
		cp SplusBModels$tag/SR2_CMS_hgg_mass.png /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag/SplusBModels/SR2_CMS_hgg_mass_bkg.png
		cp SplusBModels$tag/SR1_CMS_hgg_mass.pdf /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag/SplusBModels/SR1_CMS_hgg_mass_bkg.pdf
		cp SplusBModels$tag/SR2_CMS_hgg_mass.pdf /home/users/fsetti/public_html/HH2ggtautau/flashggFinalFit/Plots/$tag/SplusBModels/SR2_CMS_hgg_mass_bkg.pdf
	popd

}

make_plot_ws(){
	pushd Combine
		#rm -rf Models
		#mkdir -p Models
		#mkdir -p Models/signal
		#mkdir -p Models/background
		#cp ../Signal/outdir_packaged/CMS-HGG*.root ./Models/signal/
		#cp ../Background/outdir_$tag/CMS-HGG*.root ./Models/background/
		cp ../Datacard/Datacard_$tag.txt Datacard.txt
		#echo "r_singleHiggs rateParam * VH_*_hgg 1" >> Datacard.txt
		#echo "r_singleHiggs rateParam * ttH_*_hgg 1" >> Datacard.txt
		#echo "r_singleHiggs rateParam * ggH_*_hgg 1" >> Datacard.txt
		#echo "r_singleHiggs rateParam * VBFH_*_hgg 1" >> Datacard.txt
		#echo "nuisance edit freeze r_singleHiggs" >> Datacard.txt

		python RunText2Workspace.py --mode  plots --dryRun
		./t2w_jobs/t2w_plots.sh
	popd
}


#make_toys
make_SpB
#make_plot_ws