Versions Compared

Old Version 2

changes.mady.by.user Homer Neal

Saved on

compared with

New Version Current

changes.mady.by.user Homer Neal

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
Comment: Migration of unmigrated content due to installation of a new plugin

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

Image Added
Image Added
Image Added

...

TMVA driver routine and cut-table code:

sidHZZtmva-v4.C

zzhcuttabmini_v4.c

...

Preselection variables:

Code Block

Evis<140.0: (4 Jet category)
# nTrks>10 &&
# EvisJETS>60.0 &&&nbsp;
# PTvisJETS>5.0 &&&nbsp;
# hmass>50.  && hmass<140. &&&nbsp;
# ]




[^sidHZZtmva-v4.C]



factory->AddVariable( "evisjets", "Evis JETS", "GeV", 'F' );
&nbsp;&nbsp; factory->AddVariable( "mvisjets", "Mvis JETS", "GeV", 'F' );
&nbsp;&nbsp; factory->AddVariable( "ptvisjets", "PTvis JETS", "GeV/c", 'F' );
&nbsp;&nbsp; factory->AddVariable( "blike", "B Likeness", "", 'F' );
&nbsp;&nbsp; factory->AddVariable( "clike", "C likeness", "", 'F' );
&nbsp;&nbsp; factory->AddVariable( "nelec", "Number of HE Electrons", 'F' );
&nbsp;&nbsp; factory->AddVariable( "hmass", "Higgs Mass","GeV/c2", 'F' );
&nbsp;&nbsp; factory->AddVariable( "ejet1", "Z energy","GeV", 'F' );
&nbsp;&nbsp; factory->AddVariable( "ejet2",&nbsp; "Zstr energy", "GeV", 'F' );
&nbsp;&nbsp; factory->AddVariable( "cjet1",&nbsp; "Z cos theta","", 'F' );
&nbsp;&nbsp; factory->AddVariable( "cjet2",&nbsp; "Zstr cos theta", "", 'F' );
&nbsp;&nbsp; factory->AddVariable( "mjet1",&nbsp; "Z mass", "GeV/c2", 'F' );
&nbsp;&nbsp; factory->AddVariable( "mjet2",&nbsp; "Zstr mass", "GeV/c2", 'F' );
&nbsp;&nbsp; factory->AddVariable( "dphi",&nbsp; "azimuthal separation of Z jets", "", 'F' );
&nbsp;&nbsp; factory->AddVariable( "jetthrust",&nbsp; "thrust", "1", 'F' );
&nbsp;&nbsp; factory->AddVariable( "nchtrks",&nbsp; "Number Charged Tracks", "", 'F' );

&nbsp;&nbsp; factory->AddVariable( "tjpart1",&nbsp; "Number of electrons", "", 'F' );
&nbsp;&nbsp; factory->AddVariable( "tjpart2",&nbsp; "Number of muons", "", 'F' );

if ( EvisJETS<140.0) {
&nbsp;&nbsp; &nbsp; if (nTrks>10 && EvisJETS>60.0 && PTvisJETS>5.0 && hmass>50. && hmass<140. && fabs(fabs(dph)-3.14159)>0.15  &&&nbsp;
# ej1>40. &&&nbsp;
# jetthrust<0.99) presel=1; // #4 LOW EVIS
&nbsp;&nbsp;&nbsp;&nbsp; } else {
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; if (

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) presel=1; // #10
&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; }&nbsp; // common0410

\[neal@localhost weights\

The TMVA variables:

  1. Evis JETS
  2. Mvis JETS
  3. PTvis JETS
  4. B Likeness
  5. C likeness
  6. Number of HE Electrons
  7. Higgs Mass
  8. Z energy
  9. Zstr energy
  10. Z cos theta
  11. Zstr cos theta
  12. Z mass
  13. Zstr mass
  14. azimuthal separation of Z jets
  15. thrust
  16. Number Charged Tracks
  17. Number of electrons
  18. Number of muons

    Performance of different MVA options:

Code Block

 [neal@localhost weights]$ grep \-A 8 \-i optimal\- ../hzzv14-presel-v4-new-common0410results.txt
\ --\- Classifier&nbsp;&nbsp;Classifier   (&nbsp;  #signal, #backgr.)&nbsp;  Optimal-cut&nbsp;cut  S/sqrt(S+B)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NSig&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; NBkg&nbsp;&nbsp; EffSig&nbsp;&nbsp; EffBkg
\--\- \--\--\      NSig      NBkg   EffSig   EffBkg
 -- -------------------------------------------------------------------------------------------------\-
\ ---&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
\--\- Likelihood:&nbsp;      -0.0050            0         0         0        0        0
 --- Likelihood:  ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       1.0000&nbsp;&nbsp;&nbsp;&nbsp;0000     0.728991&nbsp;728991  644.0848&nbsp;0848  779980.1&nbsp;&nbsp;1   0.9699&nbsp;&nbsp;9699   0.5944
\ ---&nbsp;&nbsp;&nbsp;&nbsp;     Fisher:&nbsp;  ( 664.0744,1312202.5)&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;       0.0120&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;0120      1.70335&nbsp;70335  258.2099&nbsp;2099  22721.13&nbsp;&nbsp;13   0.3888&nbsp;3888  0.01732
\ ---&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.9823      1.17754  538.5353  208621.3    0.811    0.159
\ ---&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.0462      3.42709  100.6751  762.2891   0.1516 0.0005809
\ --\- \--\--\-----------------------------------------------------------------------------------------

...

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

Image Added

Image Added

Image Added

Cut table for BDT:

Code Block
---\-

{code}
isigall= 1345.984602 isigpresel= 14651460.839489101057 igood= 119114.854887116455 ibad= 862762.500000
 intot= 77692005.984714 ipre= 982876.354887616455
 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.11 15  Delta(sig*BR) = 96.302456548594 +/- 0.088343046879
 cut #1 (20.<PTvisJETS       ): 1.283e+07    37500  5525512    746188        0  2968250        0  3546875        0      806        69       44      259      0.23 33  Delta(sig*BR) = 43.441477042092 +/- 0.046833024172
 cut #2 (EvisJETS<220.       ): 4.637e+06    25000  3287550    514400        0   130196        0   678875        0      608        42       41      219      0.28 42  Delta(sig*BR) = 32.541462365748 +/- 0.041839020248
 cut #3 (50.<hmass<140.      ): 3.306e+06        0  2542562    378675        0    60426        0   323875        0      536        38       35      206      0.29 45  Delta(sig*BR) = 32.393216233245 +/- 0.042005019503
 cut #4 (||dph|-3.14159|>0.15 ): 2.772e+06        0  2083912    364150        0    57097        0   266000        0      514        36       33      202      0.31 47  Delta(sig*BR) = 32.238194120439 +/- 0.041071018965
 cut #5 (nTrks>10            ): 1.426e+06        0  1136375    169088        0    30283        0    89875        0      466        10       19      177      0.39 56  Delta(sig*BR) = 21.565081777889 +/- 0.034652017994
 cut #6 (jetthrust<0.99      ): 1.325e+06        0  1044700    165588        0    30033        0    84125        0      463        10       19      173      0.40 58  Delta(sig*BR) = 21.488716730704 +/- 0.033852017728
 cut #7 (ej1>40.             ): 1.313e+06        0  1035112    164625        0    29215        0    83250        0      462        10       19      173      0.40 58  Delta(sig*BR) = 21.479266725415 +/- 0.033752017704
 cut #8 (MVA                 ): 98.824e766e+02        0      162     138      625        700        0        0        0        0        0       54       51         0        1       65     62      13.71 85  Delta(sig*BR) = 0.583134259452 +/- 0.038515
{code}

[^TrackSubdetectorHitNumbersDriver.java]015456

Remaining backgrounds:

Code Block

$ 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       1 106562 4f_sz 12.500000
1278792       1 106566 4f_sw 12.500000
1278864       1 106572 4f_sz 12.500000
1278900       1 106575 4f_zz 12.500000
1279284       1 106607 2f_z_ 12.500000
2557224       2 106551 4f_ww 12.500000
2557272       2 106553 4f_zz 12.500000
2557704       2 106571 4f_sz 12.500000
2557752       2 106573 4f_zz 12.500000
3836664       3 106574 4f_zz 12.500000
3836808       3 106578 4f_ww 12.500000
12792960      10 106608 2f_z_ 12.500000
16625856      13 106576 4f_zz 12.500000
24299556      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