Update on the HZZ 250 GeV +80/-30 250/fb Analysis for the 29 August 2013 SiD Meeting:

Archive of the preparation and analysis code:

[^subproclst-data8-delay-e250-350Gflav-4jets]

[^flavortag-all-batch-revtx-350-4jets.xml] [^revertex-all-batch.xml]

[^runitallbatch2] [^suballbatch2] [^allSME250p80m30.lst] [^ffhzzE250p80m30.lst]

----
[^FastJetProcessor.cpp]

----
All difference in the LCFIPlus routines were commented out; not including these should make no difference in the output.

[^algoEtc.cc] [^FlavorTag.cc] [^flavtag.cc]

----

Analysis code:

[^MyDBDvvH.java]

[^MyDBDzHzzAnalysis.java]

[^IsolatedHighPElectronIdentifier.java]

----

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]

!variables_id_c1.png|border=1!
!variables_id_c2.png|border=1!
!variables_id_c3.png|border=1!

----

TMVA driver routine and cut-table code:

[^sidHZZtmva-v4.C]

[^zzhcuttabmini_v4.c]

----

Preselection variables:

Evis<140.0: (4 Jet category)
\# nTrks>10 &&
\# EvisJETS>60.0 &&&nbsp;
\# PTvisJETS>5.0 &&&nbsp;
\# hmass>50.&nbsp; && hmass<140. &&&nbsp;
\# fabs(fabs(dph)-3.14159)>0.15&nbsp; &&&nbsp;
\# ej1>40. &&&nbsp;
\# jetthrust<0.99

Evis>140: (6 jet category)

\# nTrks>20 &&
\# EvisJETS<220.0 &&
\# PTvisJETS>5.0&nbsp; &&
\# hmass>100. && hmass<150. &&
\# fabs(fabs(dph)-3.14159)>0.05 &&
\# ej1<130. &&
\# jetthrust<0.99

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
----

Performance of different MVA options:

&nbsp;\[neal@localhost weights\]$ grep \-A 8 \-i optimal\- ../hzzv14-presel-v4-new-common0410results.txt
&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;&nbsp;&nbsp; Cuts:&nbsp; ( 664.0744,1312202.5)&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
&nbsp;--\- Likelihood:&nbsp; ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.0000&nbsp;&nbsp;&nbsp;&nbsp; 0.728991&nbsp; 644.0848&nbsp; 779980.1&nbsp;&nbsp; 0.9699&nbsp;&nbsp; 0.5944
&nbsp;---&nbsp;&nbsp;&nbsp;&nbsp; Fisher:&nbsp; ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.0120&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.70335&nbsp; 258.2099&nbsp; 22721.13&nbsp;&nbsp; 0.3888&nbsp; 0.01732
&nbsp;---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDTG:&nbsp; ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.9823&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1.17754&nbsp; 538.5353&nbsp; 208621.3&nbsp;&nbsp;&nbsp; 0.811&nbsp;&nbsp;&nbsp; 0.159
&nbsp;---&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; BDT:&nbsp; ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; \-0.0462&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3.42709&nbsp; 100.6751&nbsp; 762.2891&nbsp;&nbsp; 0.1516 0.0005809
&nbsp;-\- \--------------------------------------------------------------------------------------------\-

----

Combined 4-jet 6-jet performance using a cut at EvisJETS of 140 GeV to determine the preselection cuts to be used:

!mvaeffs_BDT.png|border=1!

!mvaeffs_Fisher.png|border=1!

!rejBvsS.png|border=1!

Cut table for BDT:

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:

$ awk '{print $3,$7,$15}' zzhpassing.dat \| sort \-n \| uniq \-c \| awk '{if (NF==4) print $1*$2*$4,$0}' \| sort \-k 1,1 \-n
1278744&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106562 4f_sz 12.500000
1278792&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106566 4f_sw 12.500000
1278864&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106572 4f_sz 12.500000
1278900&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106575 4f_zz 12.500000
1279284&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 1 106607 2f_z\_ 12.500000
2557224&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106551 4f_ww 12.500000
2557272&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106553 4f_zz 12.500000
2557704&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106571 4f_sz 12.500000
2557752&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2 106573 4f_zz 12.500000
3836664&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106574 4f_zz 12.500000
3836808&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 3 106578 4f_ww 12.500000
12792960&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 10 106608 2f_z\_ 12.500000
16625856&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 13 106576 4f_zz 12.500000
24299556&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 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]

qqH,llH only:

[^variables_id_c1-no4j.png]

[^TrackSubdetectorHitNumbersDriver.java]

Child pages

HZZ 250 GeV Analysis - Update 1

Update on the HZZ 250 GeV +80/-30 250/fb Analysis for the 29 August 2013 SiD Meeting:

Archive of the preparation and analysis code:

Analysis code:

Analysis:

TMVA driver routine and cut-table code:

Preselection variables:

The TMVA variables:

Performance of different MVA options:

Combined 4-jet 6-jet performance using a cut at EvisJETS of 140 GeV to determine the preselection cuts to be used:

Cut table for BDT:

Remaining backgrounds:

Variable distributions without any preselection:

nu nu H only:

qqH,llH only: