Table of Contents

outline	true

Update on the HZZ 250 GeV +80/-30 250/fb Analysis for 17-19 September 2013:

Table of Contents
outline true

...

Analysis Procedure:

divide into 4 jet / 6 jet topologies
- this is ZH with H->ZZ* , Z-> nn,ll,qq
apply preselection depending on topology
train/apply TMVA'sTMVAs
validate with cut table
check remaining backgrounds

Distributions before preselection with only a cut on the reconstructed Higgs mass:

...

Preselection:

Simplified and making only clean cuts ...

Code Block

Evis<140.0: (4 Jet category)
#y34>0.0#hmass>95. && hmass<140.
#PTvisJETS>25.0 && PTvisJETS<70.0
#nTrks>5.
#ej1<120.

Evis>140: (6 jet category)
#y34>0.0
#hmass>90. && hmass<160.
#ej1<120.
#jetthrust<0.98
#nTrks>5

Distributions after preselection:

aaa

(Note: the lepton pair mass plot has the entries at less than 10 GeV randomly set to 10 and 130 GeV so that the signal peak can be seen.)

The TMVA variables:

Evis JETS
Mvis JETS
PTvis JETS
B Likeness
C likeness
Number of HE Electrons
Higgs Mass
Z energy
Zstr energy
Z cos theta
Zstr cos theta
Z mass
Zstr mass
azimuthal separation of Z jets
thrust
Number Charged Tracks
Number of electrons
Number of muons

----NEW VARIABLES:

y34
y56
lepton pair (PDG ID1 = -ID2) mass closest to Mz
jet pair mass closest to mW
jet pair1 - mW
+
jet pair2 - mW

Performance of different MVA options:

BEFORE

Some signals are just not reasonable to try to select:

1) decays like l+l-H where H->ffnunu
Even in cases where the ff are leptons the initial l+l- are more energetic
and those were being selected. Image Added
Image Added
Image Added

eeH, H -> cc ee
Image Added

nu nu H, H -> nu nu nu nu
Image Added

		HZZ generated sample stats

		events	weight
I108056.001.stdhep	10000	20000	1.29E-004	e1e1h_zz.Gwhizard-
I108055.001.stdhep	10000	20000	4.36E-003	e1e1h_zz.Gwhizard-
I108054.001.stdhep	10000	20000	4.11E-004	e1e1h_zz.Gwhizard-
I108053.001.stdhep	10000	20000	2.70E-005	e1e1h_zz.Gwhizard-
I108058.001.stdhep	10000	20000	4.29E-003	e2e2h_zz.Gwhizard-
I108057.001.stdhep	10000	20000	4.00E-004	e2e2h_zz.Gwhizard-
I108060.001.stdhep	10000	20000	4.28E-003	e3e3h_zz.Gwhizard-
I108059.001.stdhep	10000	20000	4.00E-004	e3e3h_zz.Gwhizard-
		Sum =	1.43E-002	10.40%
I108062.001.stdhep	10000	20000	2.54E-002	nnh_zz.Gwhizard-1_
I108061.001.stdhep	10000	20000	3.01E-003	nnh_zz.Gwhizard-1_
		Sum =	2.84E-002	20.68%
I108064.001.stdhep	10000	20000	8.67E-002	qqh_zz.Gwhizard-1_
I108063.001.stdhep	10000	20000	8.08E-003	qqh_zz.Gwhizard-1_
		Sum =	9.47E-002	68.93%


		All =	1.37E-001

Mass of lepton pair with mass closest to mZ:

Image Added

The following also shows the jet pair mass for those that have mass closest to mW:

Image Added

Note: Compared to last week plots and tables there was also a simple programming error which I've already fixed.

...

Performance of different MVA options:

BEFORE

Code Block


 [neal@localhost weights]$ grep -A 8 -i optimal- ../hzzv14-presel-v4-new-common0410results.txt
 --- Classifier   (  #signal, #backgr.)  Optimal-cut  S/sqrt(S+B)      NSig      NBkg   EffSig   EffBkg
 -- ---------------------

Code Block


 [neal@localhost weights]$ grep -A 8 -i optimal- ../hzzv14-presel-v4-new-common0410results.txt
 --- Classifier   (  #signal, #backgr.)  Optimal-cut  S/sqrt(S+B)      NSig      NBkg   EffSig   EffBkg
 -- ---------------------------------------------------------------------------------------------
 ---       Cuts:  ( 664.0744,1312202.5)      -0.0050            0         0         0        0        0
 --- Likelihood:  ( 664.0744,1312202.5)       1.0000     0.728991  644.0848  779980.1   0.9699   0.5944
 ---     Fisher:  ( 664.0744,1312202.5)       0.0120      1.70335  258.2099  22721.13   0.3888  0.01732
 ---       BDTG:  ( 664.0744,1312202.5)      -0.9823      1.17754  538.5353  208621.3    0.811    0.159
 ---        BDT:  ( 664.0744,1312202.5)      -0.0462      3.42709  100.6751  762.2891   0.1516 0.0005809
 -- ------------------------------------------------------------------------
 ---------------------

NOW:

Code Block


     if ((fp = fopen("HZZ-TMVA-vars-ffh_zz_all_SM_background_+80e-_-30e+-v6-mini.txt","r"))==NULL) exit(0);
[neal@localhost test]$ source storeresultsHZZv4.sh hzz-withlep-withmW-n2-Mhiggspresel-plusptntrkp8g16y34nvtx


--- ==================================================================================================
--- Classifier   (  #signal, #backgr.)  Optimal-cut  S/sqrt(S+B)      NSig      NBkg   EffSig   EffBkg
---       Cuts:  ( 664.0744,1312202.5)      -0.0050            0         0         0        0        0
 --- Likelihood:  ( 664.0744,1312202.5)       1.0000     0.728991  644.0848  779980.1   0.9699   0.5944
 ---     Fisher:  ( 664.0744,1312202.5)       0.0120      1.70335  258.2099  22721.13   0.3888  0.01732
 ---       BDTG:  ( 664.0744,1312202.5)      -0.9823      1.17754  538.5353  208621.3    0.811    0.159
 ---        BDT:  ( 664.0744,1312202.5)      -0.0462      3.42709  100.6751  762.2891   0.1516 0.0005809
 -- ---------------------------------------------------------------------------------------------

NOW:

Code Block


-----
---     if ((fp Cuts:=  (1035.8644,  1196975)      -0.0050            0         0         0        0        0
--- Likelihood:  (1035.8644,  1196975)       1.0000      1.66886  315.8824  35510.98   0.3049  0.02967
---     Fisher:  (1035.8644,  1196975fopen("HZZ-TMVA-vars-ffh_zz_all_SM_background_+80e-_-30e+-v6-mini.txt","r"))==NULL) exit(0);
[neal@localhost test]$ source storeresultsHZZv4.sh hzz-withlep-withmW-n2-Mhiggspresel-plusptntrkp8g16y34nvtx


--- ==================================================================================================
--- Classifier   (  #signal, #backgr.)  Optimal-cut  S/sqrt(S+B)      NSig 0.0047     NBkg 2.38122  586.8862EffSig  60157.76   0.5666  0.05026 EffBkg
---       BDTG:  (1035.8644,  1196975)      -0.9795        2.625  293.8675  12238.81   0.2837  0.01022
---        BDT:  (1035.8644,  1196975)      -0.0427       4.9889  146.0711  711.2005    0.141 0.0005942
--- ----------------------------------------------------------------------------------------------------------------------------

Image Removed

Image Removed
Image Removed
Image Removed Image Removed

Spent a lot of time playing with different preselection cuts and then settled for only making cuts where
I was cutting almost no signal away and then letting the MVA's do the rest.

Some signals are just not reasonable to try to select:

1) decays like l+l-H where H->ffnunu
Even in cases where the ff are leptons the initial l+l- are more energetic
and those were being selected. Image Removed
Image Removed
Image Removed

Image Removed
Image Removed
Image Removed
Changed to checking all opposite signed lepton pairs for consistency with coming
from a Z.

Image Removed

		HZZ generated sample stats

		events	weight
I108056.001.stdhep	10000	20000	1.29E-004	e1e1h_zz.Gwhizard-
I108055.001.stdhep	10000	20000	4.36E-003	e1e1h_zz.Gwhizard-
I108054.001.stdhep	10000	20000	4.11E-004	e1e1h_zz.Gwhizard-
I108053.001.stdhep	10000	20000	2.70E-005	e1e1h_zz.Gwhizard-
I108058.001.stdhep	10000	20000	4.29E-003	e2e2h_zz.Gwhizard-
I108057.001.stdhep	10000	20000	4.00E-004	e2e2h_zz.Gwhizard-
I108060.001.stdhep	10000	20000	4.28E-003	e3e3h_zz.Gwhizard-
I108059.001.stdhep	10000	20000	4.00E-004	e3e3h_zz.Gwhizard-
		Sum =	1.43E-002	10.40%
I108062.001.stdhep	10000	20000	2.54E-002	nnh_zz.Gwhizard-1_
I108061.001.stdhep	10000	20000	3.01E-003	nnh_zz.Gwhizard-1_
		Sum =	2.84E-002	20.68%
I108064.001.stdhep	10000	20000	8.67E-002	qqh_zz.Gwhizard-1_
I108063.001.stdhep	10000	20000	8.08E-003	qqh_zz.Gwhizard-1_
		Sum =	9.47E-002	68.93%


		All =	1.37E-001

2) There was also a simple programming error which I've already fixed.

Cut table for BDT: (NEW)NEW:

Code Block


 cuts__ (Cut Name                          ):       all   others        2f       4f       6f       aa       1f       3f       5f        q      e+mu      tau       nu  Signif.
 cut #0 (all                               ): 7.769e+07   625000  14507062  1019938        0 54073784        0  7464875        0       948       91       45      262      0.15  Delta(sig*BR) = 6.548594 +/-  0.046879
 cut #1 (y34>0.y34>0.                      ): 7.705e+07   625000  14449800  1018312        0 53568784        0  7387500        0       947       91       45      262      0.15  Delta(sig*BR) = 6.523157 +/-  0.046724
 cut #2 (95.<hmass<140. 95.<hmass<140.     ): 1.257e+07    150000 10264438   825238        0   387992        0   938000         0      901       86       43      213      0.35  Delta(sig*BR) =  2.854636 +/- 0.021075
 cut #3 (25.<PTvisJETS<70 PTvisJETS no cut ): 5.637e+06    25000   4545812   744212        0     8750        0   312125        0       866       84       37      180      0.49  Delta(sig*BR) = 2.033046 +/-  0.015436
 cut #4 (nTrks>5 nTrks>5                   ): 3.826e+06        0   2957175   679838        0     8750        0   179500        0       866       81       36      176      0.59  Delta(sig*BR) = 1.688536 +/-  0.012688
 cut #5 (jetthrust no cut jetthrust < 0.98 ): 3.826e+06        0   2957175   679838        0     8750        0   179500        0       866       81       36      176      0.59  Delta(sig*BR) = 1.688536 +/-  0.012688
 cut #6 (ej1<120. ej1<120.                 ): 2.144e+06        0   1480162   539312        0     8750        0   114500        0       844       78       35      176      0.77  Delta(sig*BR) = 1.292336 +/-  0.009952
 cut #7 (MVA                               ): 9.252e+02        0        50      712        0        0        0        0        0       79        27        3       53      5.35  Delta(sig*BR) = 0.186991 +/- 0.010382

Remaining backgrounds: (NEW)

Code Block


[neal@localhost weights|../../../../../../../../../]$ sed 's/\./,/g' zzhpassing.dat \| awk '{if (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 = 762,5
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106561 4f_sz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106564 4f_sw 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106573 4f_zz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106580 4f_zz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106607 2f_z\_ 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106551 4f_ww 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106562 4f_sz 12,500000
37,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106575 4f_zz 12,500000
37,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106608 2f_z\_ 12,500000
50&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4 106578 4f_ww 12,500000
62,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5 106572 4f_sz 12,500000
100&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 8 106574 4f_zz 12,500000
125&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10 106577 4f_ww 12,500000
237,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19 106576 4f_zz 12,500000

Cut table for BDT: (OLD):

Code Block


isigall= 1345.984602 isigpresel= 1460.101057 igood= 114.116455 ibad= 762.500000
intot= 77692005.984714 ipre= 876.616455
cuts_\_ (Cut Name&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ):&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; all&nbsp;&nbsp; others&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 6f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; aa&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5f&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; q&nbsp;&nbsp;&nbsp;&nbsp; e+mu&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; tau&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; nu&nbsp; Signif.
cut #0 (all&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 7.769e+07&nbsp;&nbsp; 625000 14507062&nbsp; 1019938&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0 54073784&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 7464875&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 948&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 91&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 45&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 262&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.15&nbsp; Delta(sig*BR) = 6.548594 \+/\- 0.046879
cut #1 (20.<PTvisJETS&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 1.283e+07&nbsp;&nbsp;&nbsp; 37500&nbsp; 5525512&nbsp;&nbsp; 746188&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 2968250&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 3546875&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 806&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 69&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 44&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 259&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.33&nbsp; Delta(sig*BR) = 3.042092 \+/\- 0.024172
cut #2 (EvisJETS<220.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 4.637e+06&nbsp;&nbsp;&nbsp; 25000&nbsp; 3287550&nbsp;&nbsp; 514400&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 130196&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 678875&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 608&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 42&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 41&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 219&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.42&nbsp; Delta(sig*BR) = 2.365748 \+/\- 0.020248
cut #3 (50.<hmass<140.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 3.306e+06&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 2542562&nbsp;&nbsp; 378675&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 60426&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 323875&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 536&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 38&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 35&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 206&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.45&nbsp; Delta(sig*BR) = 2.233245 \+/\- 0.019503
cut #4 (\|\|dph\|-3.14159\|>0.15 ): 2.772e+06&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 2083912&nbsp;&nbsp; 364150&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 57097&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp; 266000&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 514&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 36&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 33&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 202&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.47&nbsp; Delta(sig*BR) = 2.120439 \+/- 0.018965
cut #5 (nTrks>10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 1.426e+06&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 1136375&nbsp;&nbsp; 169088&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 30283&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 89875&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 466&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 177&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.56&nbsp; Delta(sig*BR) = 1.777889 \+/\- 0.017994
cut #6 (jetthrust<0.99&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 1.325e+06&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 1044700&nbsp;&nbsp; 165588&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 30033&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 84125&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 463&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 173&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.58&nbsp; Delta(sig*BR) = 1.730704 \+/\- 0.017728
cut #7 (ej1>40.&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 1.313e+06&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp; 1035112&nbsp;&nbsp; 164625&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 29215&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp; 83250&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 462&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 173&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.58&nbsp; Delta(sig*BR) = 1.725415 \+/\- 0.017704
cut #8 (MVA&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ): 8.766e+02&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 138&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 625&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 51&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 62&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.85&nbsp; Delta(sig*BR) = 0.259452 \+/\- 0.015456

Remaining backgrounds: (OLD)

Code Block


$ sed 's/\./,/g' zzhpassing.dat \| awk '{if (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 = 762,5
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106562 4f_sz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106566 4f_sw 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106572 4f_sz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106575 4f_zz 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106607 2f_z\_ 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106551 4f_ww 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106553 4f_zz 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106571 4f_sz 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106573 4f_zz 12,500000
37,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106574 4f_zz 12,500000
37,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106578 4f_ww 12,500000
125&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10 106608 2f_z\_ 12,500000
162,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 13 106576 4f_zz 12,500000
237,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 19 106577 4f_ww 12,500000

Background contributions for Evis>140:

Code Block


$ 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&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106551 4f_ww 12,500000
12,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106577 4f_ww 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106553 4f_zz 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106573 4f_zz 12,500000
25&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106576 4f_zz 12,500000
37,5&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106574 4f_zz 12,500000
50&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 4 106608 2f_z\_ 12,500000

Background contributions for Evis<140: (OLD)

Image Removed

Code Block


 []$ 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---------------------
---       Cuts:  (1035.8644,  1196975)      -0.0050            0         0         0        0        0
--- Likelihood:  (1035.8644,  1196975)       1.0000      1.66886  315.8824  35510.98   0.3049  0.02967
---     Fisher:  (1035.8644,  1196975)       0.0047      2.38122  586.8862  60157.76   0.5666  0.05026
---       BDTG:  (1035.8644,  1196975)      -0.9795        2.625  293.8675  12238.81   0.2837  0.01022
---        BDT:  (1035.8644,  1196975)      -0.0427       4.9889  146.0711  711.2005    0.141 0.0005942
--- --------------------------------------------------------------------------------------------------
              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 ^^^^^^

BDT plots:

Image Added

Image Added
Image Added
Image Added Image Added

...

Cut table for BDT: (NEW)NEW:

Code Block


 cuts__ (Cut Name                        1 106607 2f_z\_ 12,500000
25            2 106571 4f_sz 12,500000
37,5    ):       all   others        2f       4f       6f       aa       1f       3f       5f        q      e+mu      tau       nu  Signif.
 cut #0 (all                       3 106578 4f_ww 12,500000
75     ): 7.769e+07   625000  14507062  1019938        0 654073784        106608 2f_z\_ 12,500000
137,5        11 106576 4f_zz 12,500000
225   0  7464875        0       948       91       45      262      0.15  Delta(sig*BR) = 6.548594 +/-  0.046879
 cut #1 (y34>0.y34>0.             18 106577 4f_ww 12,500000 <**************\*

Code Block


--- Classifier&nbsp;&nbsp; (&nbsp; #signal, #backgr.)&nbsp; Optimal-cut&nbsp; S/sqrt(S+B)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NSig&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NBkg&nbsp;&nbsp; EffSig&nbsp;&nbsp; EffBkg
--- \-------------------------------------------------------------------------------------------------\-
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Cuts:&nbsp; (309.89453,&nbsp; 1128865)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.0050&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0
--- Likelihood:&nbsp; (309.89453,&nbsp; 1128865)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0000&nbsp;&nbsp;&nbsp;&nbsp; 0.378385&nbsp; 299.2573&nbsp; 625191.7&nbsp;&nbsp; 0.9657&nbsp;&nbsp; 0.5538
\---&nbsp;&nbsp;&nbsp;&nbsp; Fisher:&nbsp; (309.89453,&nbsp; 1128865)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0088&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.63217&nbsp; 150.0501&nbsp; 8301.587&nbsp;&nbsp; 0.4842 0.007354
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDTG:&nbsp; (309.89453,&nbsp; 1128865)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.9950&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.40141&nbsp; 134.4158&nbsp; 9065.234&nbsp;&nbsp; 0.4337&nbsp; 0.00803
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDT:&nbsp; (309.89453,&nbsp; 1128865)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0177&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.63521&nbsp;&nbsp; 38.5486&nbsp; 73.90136&nbsp;&nbsp; 0.1244 6.547e-05
--- \-------------------------------------------------------------------------------------------------\-

after raising reco higgs mass cut from 50 GeV to 100 GeV:

&nbsp;--\- Classifier&nbsp;&nbsp; (&nbsp; #signal, #backgr.)&nbsp; Optimal-cut&nbsp; S/sqrt(S+B)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NSig&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NBkg&nbsp;&nbsp; EffSig&nbsp;&nbsp; EffBkg
--- \-------------------------------------------------------------------------------------------------\-
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Cuts:&nbsp; (209.30833,&nbsp;&nbsp;&nbsp; 51800)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.0050&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0
--- Likelihood:&nbsp; (209.30833,&nbsp;&nbsp;&nbsp; 51800)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.9105&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.08931&nbsp; 191.9554&nbsp; 30860.64&nbsp;&nbsp; 0.9171&nbsp;&nbsp; 0.5958
\---&nbsp;&nbsp;&nbsp;&nbsp; Fisher:&nbsp; (209.30833,&nbsp;&nbsp;&nbsp; 51800)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0176&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.78558&nbsp; 116.8125&nbsp; 4162.945&nbsp;&nbsp; 0.5581&nbsp; 0.08037
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDTG:&nbsp; (209.30833,&nbsp;&nbsp;&nbsp; 51800)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.5098&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.00598&nbsp; 204.0195&nbsp; 40926.24&nbsp;&nbsp; 0.9747&nbsp;&nbsp; 0.7901
\---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDT:&nbsp; (209.30833,&nbsp;&nbsp;&nbsp; 51800)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0361&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2.92564&nbsp;&nbsp; 72.9263&nbsp; 548.4119&nbsp;&nbsp; 0.3484&nbsp; 0.01059
--- \-------------------------------------------------------------------------------------------------\-

...

Evis<140 GeV and only idrup = 106577 for the background: (OLD)

Image Removed
Image Removed
Image Removed

       ): 7.705e+07   625000  14449800  1018312        0 53568784        0  7387500        0       947       91       45      262      0.15  Delta(sig*BR) = 6.523157 +/-  0.046724
 cut #2 (95.<hmass<140. 95.<hmass<140.     ): 1.257e+07    150000 10264438   825238        0   387992        0   938000         0      901       86       43      213      0.35  Delta(sig*BR) =  2.854636 +/- 0.021075
 cut #3 (25.<PTvisJETS<70 PTvisJETS no cut ): 5.637e+06    25000   4545812   744212        0     8750        0   312125        0       866       84       37      180      0.49  Delta(sig*BR) = 2.033046 +/-  0.015436
 cut #4 (nTrks>5 nTrks>5                   ): 3.826e+06        0   2957175   679838        0     8750        0   179500        0       866       81       36      176      0.59  Delta(sig*BR) = 1.688536 +/-  0.012688
 cut #5 (jetthrust no cut jetthrust < 0.98 ): 3.826e+06        0   2957175   679838        0     8750        0   179500        0       866       81       36      176      0.59  Delta(sig*BR) = 1.688536 +/-  0.012688
 cut #6 (ej1<120. ej1<120.                 ): 2.144e+06        0   1480162   539312        0     8750        0   114500        0       844       78       35      176      0.77  Delta(sig*BR) = 1.292336 +/-  0.009952
 cut #7 (MVA                               ): 9.252e+02        0        50      712        0        0        0        0        0       79        27        3       53      5.35  Delta(sig*BR) = 0.186991 +/- 0.010382

Remaining backgrounds: (NEW)

Code Block


$ sed 's/\./,/g' zzhpassing.dat | awk '{if (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 = 762,5
12,5          1 106561 4f_sz 12,500000
12,5          1 106564 4f_sw 12,500000
12,5          1 106573 4f_zz 12,500000
12,5          1 106580 4f_zz 12,500000
12,5          1 106607 2f_z_ 12,500000
25            2 106551 4f_ww 12,500000
25            2 106562 4f_sz 12,500000
37,5          3 106575 4f_zz 12,500000
37,5          3 106608 2f_z_ 12,500000
50            4 106578 4f_ww 12,500000
62,5          5 106572 4f_sz 12,500000
100           8 106574 4f_zz 12,500000
125          10 106577 4f_ww 12,500000
237,5        19 106576 4f_zz 12,500000  <----

Plans:

just noticed that the preselection nTrks cut is on the PFO count and not the #charged tracks distribution shown so I will switch to
using the charged tracks count in the preselection
may still be able to do better with the leptonic Z decay selection
- high evis events with well reconstructed leptonic Z decays and a clear 6 jet topology (using y56 cut)
write report

...

Child pages

Versions Compared

Old Version 9

New Version Current

Key

Update on the HZZ 250 GeV +80/-30 250/fb Analysis for 17-19 September 2013:

Table of Contents
outline true

Analysis Procedure:

Distributions before preselection with only a cut on the reconstructed Higgs mass:

Preselection:

Distributions after preselection:

The TMVA variables:

Performance of different MVA options:

Some signals are just not reasonable to try to select:

Performance of different MVA options:

Some signals are just not reasonable to try to select:

Cut table for BDT: (NEW)NEW:

Remaining backgrounds: (NEW)

Cut table for BDT: (OLD):

Remaining backgrounds: (OLD)

Background contributions for Evis>140:

Background contributions for Evis<140: (OLD)

Cut table for BDT: (NEW)NEW:

Evis<140 GeV and only idrup = 106577 for the background: (OLD)

Remaining backgrounds: (NEW)

Plans:

Child pages

Page History

Versions Compared

Old Version 9

New Version Current

Key

Update on the HZZ 250 GeV +80/-30 250/fb Analysis for 17-19 September 2013:

Table of Contentsoutlinetrue

Analysis Procedure:

Distributions before preselection with only a cut on the reconstructed Higgs mass:

Preselection:

Distributions after preselection:

The TMVA variables:

Performance of different MVA options:

Some signals are just not reasonable to try to select:

Performance of different MVA options:

Some signals are just not reasonable to try to select:

Cut table for BDT: (NEW)NEW:

Remaining backgrounds: (NEW)

Cut table for BDT: (OLD):

Remaining backgrounds: (OLD)

Background contributions for Evis>140:

Background contributions for Evis<140: (OLD)

Cut table for BDT: (NEW)NEW:

Evis<140 GeV and only idrup = 106577 for the background: (OLD)

Remaining backgrounds: (NEW)

Plans:

Table of Contents
outline true