TagProbeSF/calcSF_shift.sh at master · cmantill/TagProbeSF · 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
78
#!/bin/bash
object=$1
algo=$2
wp=$3
whichbit=$4 #For bit version, enter either 04 or 05
scale=$5
scaledat=$6
smear=$7 # by default should be 0.5
smeardat=$8 # by default should be 0.1

declare -a ptranges
ptranges=("incl")

workdir=${object}"_"${algo}
mkdir ${workdir}

## loop over the pt bins
for ptrange in "${ptranges[@]}";
do
    ##make templates
    echo "make templates"
    echo "using bit version ${whichbit}"
    rm -r -f T_n2/
    rm -r -f T_n2_scale${scale}_${scaledat}_smear${smear}_${smeardat}/

    # for bacon skim -> comment this if the templates are already done in coffea
    nbins=20
    xmin=50
    xmax=120
    cmdpass=$(echo 'makeSFTemplates.C("'${object}'","'${algo}'","'${wp}'","'${ptrange}'","'${whichbit}'",true,'${nbins}','${xmin}','${xmax}','${scale}')')
    cmdfail=$(echo 'makeSFTemplates.C("'${object}'","'${algo}'","'${wp}'","'${ptrange}'","'${whichbit}'",false,'${nbins}','${xmin}','${xmax}','${scale}')')
    root -l -q ${cmdpass}
    root -l -q ${cmdfail}

    # pre-templates done

    ## add scale and smear
    inputname=${object}"_"${algo}"_"${wp}"_"${whichbit}"_"${ptrange}
    #python makeSmearShift.py --ifile ${workdir}/${inputname}_pass_pre.root --smear ${7} --scale ${5}
    #python makeSmearShift.py --ifile ${workdir}/${inputname}_fail_pre.root --smear ${7} --scale ${5}

    mv ${workdir}/${inputname}_pass_pre.root ${workdir}/${inputname}_pass.root
    mv ${workdir}/${inputname}_fail_pre.root ${workdir}/${inputname}_fail.root

    ## make datacard
    echo "make datacard"
    echo " "
    cd ${workdir}
    inputname=${object}"_"${algo}"_"${wp}"_"${whichbit}"_"${ptrange}
    cp ../makeSFDatacard.C .
    cmdmakedatacard=$(echo 'makeSFDatacard.C("'${inputname}'",'${scaledat}','${smeardat}')')
    root -l -q ${cmdmakedatacard} > sf.txt
    sed -n -i '3,$ p' sf.txt

    ## do the tag and probe
    echo "run the tag and probe"
    if [ ${object} == "T" ];
    then
	text2workspace.py -m 125 -P HiggsAnalysis.CombinedLimit.TagAndProbeExtended:tagAndProbe sf.txt --PO categories=catp2,catp1
    elif [ ${object} == "W" ];
    then
	text2workspace.py -m 125 -P HiggsAnalysis.CombinedLimit.TagAndProbeExtended:tagAndProbe sf.txt --PO categories=catp2,catp1
    fi
    mv sf.root sf"_"${inputname}".root"
    echo "Do the MultiDimFit"
    combine -M MultiDimFit -m 125 sf"_"${inputname}.root --algo=singles --robustFit=1 --cminDefaultMinimizerTolerance 5.
    echo "Run the FitDiagnostics"
    combine -M FitDiagnostics -m 125 sf"_"${inputname}.root --saveShapes --saveWithUncertainties --robustFit=1 --cminDefaultMinimizerTolerance 5.
    mv fitDiagnostics.root sf"_"fitDiagnostics"_"${inputname}".root"
    mv sf.txt sf"_"datacard"_"${inputname}".txt"
    cd ../

    # plots
    cmdmake=$(echo 'makePlots.C("'${workdir}'","'${inputname}'","'${wp}'","'${ptrange}'","'${whichbit}'",'${scale}','${scaledat}','${smear}','${smeardat}','${xmin}','${xmax}.','${nbins}',"mass")')
    root -l -q ${cmdmake}
done

mv T_n2/ T_n2_scale${scale}gev_${scaledat}_smear${smear}_${smeardat}/