Update on the HZZ 250 GeV +80/-30 250/fb Analysis for the
...
5 September 2013 SiD Meeting:
Table of Contents | ||
---|---|---|
|
...
Analysis:
- divide into 4 jet / 6 jet topologies
- apply preselection depending on topology
- train/apply TMVA's
- validate with cut table
- check remaining backgrounds
^variables_id_c1.png|||||\ ^variables_id_c2.png|||||\ ^variables_id_c3.png|||||\
...
Preselection variables:
Code Block |
---|
Evis<140.0: (4 Jet category) # nTrks>10 && # EvisJETS>60.0 && # PTvisJETS>5.0 && # hmass>50. && hmass<140. && # fabs(fabs(dph)-3.14159)>0.15 && # ej1>40. && # jetthrust<0.99 Evis>140: (6 jet category) # nTrks>20 && # EvisJETS<220.0 && # PTvisJETS>5.0 && # hmass>100. && hmass<150. && # fabs(fabs(dph)-3.14159)>0.05 && # ej1<130. && # jetthrust<0.99 |
...
Remaining backgrounds:
Code Block |
---|
$ awk '{print sed 's/\./,/g' zzhpassing.dat | awk '{if (NF==18) print $3,$7,$15}' zzhpassing.dat | sort -n | uniq -c | awk '{if (NF==4) aa+=($1*$4);print $1*$2*$4,"\t"$0}' | sort -k END{print "Sum of weights = "aa}' | sort -k 1,1 -n 1278744 Sum of weights = 762,5 12,5 1 106562 4f_sz 12.,500000 1278792 12,5 1 106566 4f_sw 12.,500000 1278864 12,5 1 106572 4f_sz 12.,500000 1278900 12,5 1 106575 4f_zz 12.,500000 1279284 12,5 1 106607 2f_z_ 12.,500000 2557224 25 2 106551 4f_ww 12.,500000 2557272 25 2 106553 4f_zz 12.,500000 2557704 25 2 106571 4f_sz 12.,500000 2557752 25 2 106573 4f_zz 12.,500000 3836664 37,5 3 106574 4f_zz 12.,500000 3836808 37,5 3 106578 4f_ww 12.,500000 12792960 125 10 106608 2f_z_ 12.,500000 16625856 162,5 13 106576 4f_zz 12.,500000 24299556 237,5 19 106577 4f_ww 12.,500000 |
Variable distributions without any preselection:
^variables_id_c1-nopresel.png|||||\ ^variables_id_c2-nopresel.png|||||\ ^variables_id_c3-nopresel.png|||||\
nu nu H only:
^variables_id_c1-4jonly.png|||||\
Code Block |
---|
--- Classifier ( #signal, #backgr.) Optimal-cut S/sqrt(S+B) NSig NBkg EffSig EffBkg
--- --------------------------------------------------------------------------------------------------
--- Cuts: (172.67554, 53762.5) 0.0350 2.42407 5.876106 0 0.03403 0
--- Likelihood: (172.67554, 53762.5) 0.7786 0.910341 161.7188 31396.51 0.9365 0.584
--- Fisher: (172.67554, 53762.5) 0.0942 5.84493 34.16319 0 0.1978 0
--- BDTG: (172.67554, 53762.5) 0.4529 4.48342 20.10101 0 0.1164 0
--- BDT: (172.67554, 53762.5) 0.1717 5.33182 28.42833 0 0.1646 0
--- --------------------------------------------------------------------------------------------------
|
qqH,llH only:
|
Background contributions for Evis>140:
Code Block |
---|
[neal@localhost weights]$ sed 's/\./,/g' zzhpassing.dat | awk '{if ($18>=140. && NF==18) print $3,$7,$15}' | sort -n | uniq -c | awk '{if (NF==4) aa+=($1*$4);print $1*$4"\t"$0} END{print "Sum of weights = "aa}' | sort -k 1,1 -n
Sum of weights = 187,5
12,5 1 106551 4f_ww 12,500000
12,5 1 106577 4f_ww 12,500000
25 2 106553 4f_zz 12,500000
25 2 106573 4f_zz 12,500000
25 2 106576 4f_zz 12,500000
37,5 3 106574 4f_zz 12,500000
50 4 106608 2f_z_ 12,500000
|
Background contributions for Evis<140:
Code Block |
---|
[neal@localhost weights]$ sed 's/\./,/g' zzhpassing.dat | awk '{if ($18<140. && NF==18) print $3,$7,$15}' | sort -n | uniq -c | awk '{if (NF==4) aa+=($1*$4);print $1*$4"\t"$0} END{print "Sum of weights = "aa}' | sort -k 1,1 -n
Sum of weights = 575
12,5 1 106551 4f_ww 12,500000
12,5 1 106562 4f_sz 12,500000
12,5 1 106566 4f_sw 12,500000
12,5 1 106572 4f_sz 12,500000
12,5 1 106575 4f_zz 12,500000
12,5 1 106607 2f_z_ 12,500000
25 2 106571 4f_sz 12,500000
37,5 3 106578 4f_ww 12,500000
75 6 106608 2f_z_ 12,500000
137,5 11 106576 4f_zz 12,500000
225 18 106577 4f_ww 12,500000 <***************
|
Code Block |
---|
--- Classifier ( #signal, #backgr.) Optimal-cut S/sqrt(S+B) NSig NBkg EffSig EffBkg
--- ------ |
Code Block |
--- Classifier ( #signal, #backgr.) Optimal-cut S/sqrt(S+B) NSig NBkg EffSig EffBkg --- -------------------------------------------------------------------------------------------------- --- Cuts: (491.39883, 1258440) -0.0050 0 0 0 0 0 --- Likelihood: (491.39883, 1258440) 1.0000 0.575508 479.389 693382.9 0.9756 0.551 --- Fisher: (491.39883, 1258440) 0.0274 1.59132 2.532298 0 0.005153 0 --- BDTG: (491.39883, 1258440) -0.9906 0.918097 174.4497 35930.22 0.355 0.02855 --- BDT: (491.39883, 1258440) 0.0076 2.53129 6.407406 0 0.01304 0 --- -------------------------------------------------------------------------------------------------- |
Background contributions for Evis<140:
Code Block |
---|
[neal@localhost weights]$ sed 's/\./,/g' zzhpassing.dat | awk '{if ($18<140. && NF==18) print $3,$7,$15}' | sort -n | uniq -c | awk '{if (NF==4) aa+=($1*$4);print $1*$4"\t"$0} END{print "Sum of weights = "aa}' | sort -k 1,1 -n
Sum of weights = 575
12,5 1 106551 4f_ww 12,500000
12,5 1 106562 4f_sz 12,500000
12,5 1 106566 4f_sw 12,500000
12,5 1 106572 4f_sz 12,500000
12,5 1 106575 4f_zz 12,500000
12,5 1 106607 2f_z_ 12,500000
25 2 106571 4f_sz 12,500000
37,5 3 106578 4f_ww 12,500000
75 6 106608 2f_z_ 12,500000
137,5 11 106576 4f_zz 12,500000
225 18 106577 4f_ww 12,500000
|
Background contributions for Evis>140:
...
----------
--- Cuts: (309.89453, 1128865) -0.0050 0 0 0 0 0
--- Likelihood: (309.89453, 1128865) 1.0000 0.378385 299.2573 625191.7 0.9657 0.5538
--- Fisher: (309.89453, 1128865) 0.0088 1.63217 150.0501 8301.587 0.4842 0.007354
--- BDTG: (309.89453, 1128865) -0.9950 1.40141 134.4158 9065.234 0.4337 0.00803
--- BDT: (309.89453, 1128865) 0.0177 3.63521 38.5486 73.90136 0.1244 6.547e-05
--- --------------------------------------------------------------------------------------------------
after raising reco higgs mass cut from 50 GeV to 100 GeV:
--- Classifier ( #signal, #backgr.) Optimal-cut S/sqrt(S+B) NSig NBkg EffSig EffBkg
--- --------------------------------------------------------------------------------------------------
--- Cuts: (209.30833, 51800) -0.0050 0 0 0 0 0
--- Likelihood: (209.30833, 51800) 0.9105 1.08931 191.9554 30860.64 0.9171 0.5958
--- Fisher: (209.30833, 51800) 0.0176 1.78558 116.8125 4162.945 0.5581 0.08037
--- BDTG: (209.30833, 51800) -0.5098 1.00598 204.0195 40926.24 0.9747 0.7901
--- BDT: (209.30833, 51800) 0.0361 2.92564 72.9263 548.4119 0.3484 0.01059
--- --------------------------------------------------------------------------------------------------
|
Evis<140 GeV and only idrup = 106577 for the background:
Potential tools for reducing the 4f_WW background:
- higher cut on number of charged tracks (unlikely to help S)
- extra cuts but keep any event that has high B-likeness
- redoing n-tuples adding variables for jet combos consistent with W decays (IN PROGRESS)
Note concerning ZZ background: wider cut on acoplanarity should help