March 25

Some test numbers

ZZ->2l2nu WZ->2l2nu WW->2l2nu ttbar Z+jets BG Observed
single lepton 767.03 3772.97 20404.73 80733.47 11501041.70 11606719.91 3405157
HLT selection/weights 767.03 3772.97 20404.73 80733.47 11501041.70 11606719.91 2991924
Good PV 758.30 3747.46 20321.69 79989.98 11335481.19 11440298.62 2917298
Data quality 758.30 3747.46 20321.69 79989.98 11335481.19 11440298.62 2917298
Cosmics rejection 758.24 3747.14 20320.74 79987.04 11331315.70 11436128.86 2913076
Dilepton 138.67 340.90 832.77 3791.25 1293789.66 1298893.24 1322224
Z mass 128.48 260.97 194.88 827.36 1198804.75 1200216.45 1207984

February 16

ADC studies

Some things that would be useful to implement

If I am sampling, at least, at 25 ns/sample, this will require I collect 40 MS/s. This is well within the capacity of the ADC (500 MS/s). This means I could collect ~11000 orbits, process the data,

  • Modify sample size: It would be useful to increase the sample size (currently 2 ns) to be closer to the size of a bunch (25 ns). This means the

  • Zero length encoding: only 14 control words...

January 27

BCM1F Lumi

Look through lumi talks. For the talk next week...

HF

  • 2 slide description of system
  • how it is used in measurement
  • summarize statistical/systematic errors
  • how it fails at max lumi

Rates

  • total
  • per channel
  • gated rates (8 ns from bptx) -- describe reduction in rate
  • lumi versus mu (per bunch crossing rate)
  • (get gated rate data from David and, if needed, TDC data from Roberval)

Plans

  • Define variable that is proportional to lumi (gated rate diamond-by-diamond (1 or more, 2 or more, etc.), gated OR, zero counting...)
  • understand bunch-by-bunch statistics (what is publishing rate that achieves 10% error)
  • systematics
    • deadtime (front-end saturation, time to return to baseline)
    • efficiency
    • bcm1f acceptance
  • sketch out data flow (gated rates to "tdc" data to software)
    • bin requirements (bunch clock)
    • timebase orbit trigger
    • time over threshold

January 23

ADC tests

Testing ADC to be used for BCM1F lumi measurement. Basic idea is to use the ADC as a TDC, i.e. reconstruct the bunch structure to obtain bunch-by-bunch lumi measurement. Test setup consists of the following:

  • LeCroy scope (400MHz)
  • Stanford pulse generator
  • Tektroniks arbitrary signal generator
  • VME crate with VME-CONET bridge and V1721 digitizer (ADC). Also plugged in, but not operational (currently) are two TDCs.

With this setup, I am recreating an approximation to the beam bunch structure. I use the stanford to simulate the orbit clock (f~11,111 Hz) and consequently as the trigger source for the signal generator, the scope, and the ADC. The current setting for the signal generator is:

  • period = 100 ns
  • amplitude = 1 Vpp
  • pulse/train = 40
  • pulse width = 10 ns
  • delay wrt trigger = 25 μs

The triggering of the signal generator and the ADC are in time so it is likely that there is a ~5 ns offset in the readout due to travel time from the signal generator to the ADC.

Dec 13

OF yields

WZ WW ttbar tW Z+jets BG Observed
single lepton 3688.56 21050.72 83996.98 78945.42 10560081.67 10747763.35 40315234
HLT selection/weights 3688.56 21050.72 83996.98 78945.42 10560081.67 10747763.35 27266446
Good PV 3686.02 21030.11 83992.51 78942.27 10544636.08 10732287.00 27202567
Data quality 3686.02 21030.11 83992.51 78942.27 10544636.08 10732287.00 27135665
Cosmics rejection 3685.56 21028.66 83988.54 78941.68 10538281.48 10725925.92 27135634
Dilepton 215.65 2375.43 8928.20 974.37 9358.36 21852.00 21885
Z mass 45.56 555.57 1987.64 214.12 807.73 3610.63 3787
Z pT 21.45 162.26 1090.57 117.84 55.76 1447.88 1484
3rd lepton veto 9.19 161.74 1031.12 113.45 38.41 1353.90 1396
b-jet veto 9.19 161.74 1031.12 113.45 38.41 1353.90 1396

WZ WW ttbar tW Z+jets BG Observed
0 b-jet 8.33 ± 0.21 ± 1.00 149.94 ± 2.04 ± 0.03 164.32 ± 1.38 36.24 ± 1.59 33.42 ± 4.50 392.24 ± 5.46 425
1 b-jet 0.84 ± 0.07 ± 0.10 11.16 ± 0.57 ± 0.01 495.35 ± 2.40 64.79 ± 2.13 4.95 ± 1.77 577.10 ± 3.71 583
2 b-jet 0.02 ± 0.01 ± 0.00 0.61 ± 0.13 ± 0.00 348.56 ± 2.01 11.99 ± 0.91 0.04 ± 0.04 361.22 ± 2.21 362
3 b-jet 0.00 ± 0.00 ± 0.00 0.03 ± 0.03 ± 0.00 21.40 ± 0.50 0.43 ± 0.17 0.00 ± 0.00 21.86 ± 0.52 26

Dec 8

OF yields for 2011

ZZ WZ WW ttbar tW Z+jets BG Observed
single lepton 729.43 3797.37 21671.72 86474.91 81274.33 10871605.79 11065553.55 40315234
HLT selection/weights 729.43 3797.37 21671.72 86474.91 81274.33 10871605.79 11065553.55 18053416
Good PV 728.64 3794.76 21650.50 86470.31 81271.08 10855704.56 11049619.84 18042209
Data quality 728.64 3794.76 21650.50 86470.31 81271.08 10855704.56 11049619.84 18002888
Cosmics rejection 728.54 3794.28 21649.01 86466.22 81270.47 10849162.49 11043071.02 18002869
Dilepton 1.52 222.01 2445.50 9191.58 1003.11 9634.43 22498.16 21533
Z mass 0.11 46.91 571.96 2046.28 220.44 831.56 3717.26 3730
Z pT 0.05 22.08 167.05 1122.74 121.32 57.40 1490.64 1462
3rd lepton veto 0.04 9.42 166.51 1061.56 116.88 39.54 1393.95 1374
b-jet veto 0.04 8.53 154.36 169.19 37.31 34.40 403.84 418

Dec 5

photonPt_Combined.png

Dec 1

2011 combined numbers (Muons)

ZZ WZ WW ttbar tW Z+jets BG Observed
single lepton 666.01 3441.60 19641.33 78368.64 73659.86 9853061.19 10028838.63 13955484
HLT selection/weights 666.01 3441.60 19641.33 78368.64 73659.86 9853061.19 10028838.63 9087316
Good PV 665.37 3439.24 19622.10 78364.47 73656.91 9838649.72 10014397.80 9034579
Data quality 665.37 3439.24 19622.10 78364.47 73656.91 9838649.72 10014397.80 9015963
Cosmics rejection 665.26 3438.80 19620.75 78360.77 73656.36 9832720.57 10008462.52 9009947
Dilepton 190.47 462.23 1230.89 4695.02 513.56 1742771.87 1749864.05 1877609
Z mass 176.00 354.29 287.75 1053.81 113.87 1603297.18 1605282.91 1691815
Z pT 81.74 168.60 84.47 570.96 63.95 93285.26 94254.97 98498
3rd lepton veto 81.67 96.17 84.34 567.58 63.59 93151.84 94045.18 98228
b-jet veto 77.82 86.99 78.44 88.27 20.50 81095.20 81447.23 84388

ZZ WZ WW ttbar tW Z+jets BG Signal Observed
250 23.16 ± 0.25 ± 2.55 16.72 ± 0.29 ± 2.01 18.24 ± 0.69 ± 0.03 19.29 ± 0.46 5.77 ± 0.62 16.48 ± 3.25 99.66 ± 4.73 15.50 ± 2.13 100
300 15.01 ± 0.20 ± 1.65 8.63 ± 0.21 ± 1.04 3.80 ± 0.31 ± 0.03 5.40 ± 0.24 1.74 ± 0.36 1.56 ± 0.82 36.14 ± 2.20 13.66 ± 2.12 45
350 10.00 ± 0.17 ± 1.10 5.06 ± 0.16 ± 0.61 0.93 ± 0.16 ± 0.03 1.08 ± 0.11 0.46 ± 0.19 1.18 ± 0.83 18.72 ± 1.55 13.28 ± 2.50 18
400 7.73 ± 0.15 ± 0.85 3.30 ± 0.13 ± 0.40 0.20 ± 0.08 ± 0.04 0.30 ± 0.06 0.23 ± 0.13 0.59 ± 0.59 12.36 ± 1.14 11.19 ± 1.82 14
500 4.75 ± 0.12 ± 0.52 1.43 ± 0.08 ± 0.17 0.00 ± 0.00 ± 0.00 0.03 ± 0.02 0.00 ± 0.00 0.00 ± 0.00 6.21 ± 0.57 4.70 ± 0.83 9
600 2.77 ± 0.10 ± 0.31 0.90 ± 0.07 ± 0.11 0.00 ± 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 1.17 ± 0.92 4.85 ± 0.98 2.12 ± 0.40 3

Nov23

Muons

  • 2011B

ZZ WZ WW ttbar tW Z+jets (data) BG Observed
0 jet 2.12 ± 0.05 ± 0.23 4.50 ± 0.11 ± 0.54 3.15 ± 0.23 ± 0.00 163.02 ± 1.29 18.98 ± 0.85 5884.45 ± 74.81 6057.21 ± 74.82 6601
1 jet 0.11 ± 0.01 ± 0.01 0.29 ± 0.03 ± 0.04 0.21 ± 0.06 ± 0.00 108.17 ± 0.92 3.56 ± 0.38 394.69 ± 19.04 503.47 ± 19.07 507
2 jet 0.01 ± 0.00 ± 0.00 0.02 ± 0.01 ± 0.00 0.00 ± 0.00 ± 0.00 7.02 ± 0.23 0.21 ± 0.08 17.94 ± 2.87 24.98 ± 2.88 26
3 jet 0.00 ± 0.00 ± 0.00 0.00 ± 0.00 ± 0.00 0.00 ± 0.00 ± 0.00 0.51 ± 0.06 0.00 ± 0.00 0.61 ± 0.15 1.12 ± 0.16 3

  • 2011A

ZZ WZ WW ttbar tW Z+jets (data) BG Observed
0 jet 1.88 ± 0.04 ± 0.21 3.90 ± 0.09 ± 0.47 2.42 ± 0.16 ± 0.01 142.01 ± 1.03 17.36 ± 0.70 5119.47 ± 27.39 5269.65 ± 27.42 5799
1 jet 0.10 ± 0.01 ± 0.01 0.26 ± 0.02 ± 0.03 0.13 ± 0.04 ± 0.00 93.76 ± 0.71 2.89 ± 0.28 276.65 ± 5.70 370.90 ± 5.74 435
2 jet 0.00 ± 0.00 ± 0.00 0.02 ± 0.01 ± 0.00 0.00 ± 0.00 ± 0.00 5.96 ± 0.18 0.21 ± 0.08 12.68 ± 1.38 18.67 ± 1.39 27
3 jet 0.00 ± 0.00 ± 0.00 0.00 ± 0.00 ± 0.00 0.00 ± 0.00 ± 0.00 0.40 ± 0.04 0.00 ± 0.00 0.37 ± 0.13 0.78 ± 0.14 5

Nov 17

Systematics for HZZ2l2nu

  • Everything from them except,

    • top background - Done by Radek
    • WW - take values from HWWlnulnu
    • pile-up - check effect of using low PU versus high PU events
    • third letpon veto - !?

Test table

ZZ WZ WW top Z+jets (data) BG Observed Signal
250 14.93 ± 0.25 ± 1.64 12.23 ± 0.33 ± 1.47 10.05 ± 0.75 ± 0.30 17.94 ± 4.84 ± 2.35 13.53 ± 2.06 ± 2.71 67.81 ± 5.87 76.0 9.18 ± 0.10
300 7.34 ± 0.18 ± 0.81 4.40 ± 0.20 ± 0.53 1.31 ± 0.28 ± 0.33 3.37 ± 0.92 ± 0.44 8.62 ± 1.72 ± 2.09 24.89 ± 2.78 19.0 6.91 ± 0.08
350 5.14 ± 0.15 ± 0.57 2.73 ± 0.16 ± 0.33 0.04 ± 0.01 ± 0.50 0.76 ± 0.21 ± 0.10 3.19 ± 1.03 ± 0.91 11.83 ± 1.86 4.0 6.91 ± 0.07
400 3.81 ± 0.13 ± 0.42 1.78 ± 0.13 ± 0.21 0.01 ± 0.01 ± 0.50 0.21 ± 0.06 ± 0.03 2.75 ± 1.05 ± 0.90 8.55 ± 1.79 7.0 5.92 ± 0.06
500 2.27 ± 0.10 ± 0.25 0.92 ± 0.09 ± 0.11 0.01 ± 0.01 ± 0.00 0.02 ± 0.01 ± 0.00 1.98 ± 0.90 ± 0.82 5.21 ± 1.41 6.0 2.63 ± 0.03
600 1.23 ± 0.08 ± 0.14 0.40 ± 0.06 ± 0.05 0.01 ± 0.01 ± 0.00 0.00 ± 0.00 ± 0.00 0.45 ± 0.45 ± 0.23 2.11 ± 0.79 2.0 1.03 ± 0.01

Nov 6

We want to simulate individual leptons for the Z+jet from photon data background estimate. The current approach being pursued by Andrey is to save the lepton pairs from dilepton events passing our preselection and then do some mapping when the photon sample is run over. There are a few issues with this:

  • You need to save a separate ntuple and then loop over it before you analyzer the photon sample. This means we have spend more time running code and take up more space with some extra file.

  • To make the association you bin the lepton pairs by the qt of the Z candidate and then randomly select pairs of leptons from the bin based on the photon pt. This means you'll get the right qt, but the phi and eta of the Z won't be the same as that of the photon. This is an issue if you want to make angular correlations between the leptons and the MET or jets. This, I believe, is a serious limitation.

The first issue is purely technical whereas the second will mess up our physics and calls the whole procedure into question. So I propose an alternate method that should simplify the association and take care of any issues we might have with angular correlations. The idea is to, in essence, randomly generate a lead lepton (as a Lorentz vector) based on information taken from the double muon/electron data and the photon 4-momentum. Since you know the 4-momentum of the both the photon and the lead lepton you get the trailing lepton from the difference. Naively, I would do the following:

  • Get the pt, eta, and phi distributions of the lead lepton from the analysis sample

  • Randomly sample these histograms (in the same fashion we do to assign a photon mass) and generate the 4-momentum for the lead lepton.

  • Take the difference of the photon P4 and the lead lepton P4 to derive the trailing lepton P4.

This takes care of the first issue, but suffers from the second issue. Additionally, we now ignore correlations between the lead lepton pt and the Z pt. To address these I do the following:

  • To include the pt correlations, I parameterize the lead lepton pt as a function of the Z pt. Technically, this means that I store a 2D histogram and randomly sample from the profile slice corresponding to the photon/Z pt. (Alternatively, I can use what I do for angular correlations.)

  • To include angular correlations, I randomly sample the Δφ and Δη between the lead lepton and the Z instead just the φ and η of the lead lepton. Whatever value I get is then just added to photon η/φ to get the value for the lepton.

Now to test it...

UPDATE They had actually already taken care of the issue with the angular correlations of the lepton and the Z being preserved after mapping to the photon. This is done by rotating the dilepton to match the phi of the Z. I had realized this was an option, but had not implemented it.

Nov ?

tW ttbar GluGluWW WW WZ ZZ PhotonJets BG DATA
single lepton 35983.60 38283.89 303.20 8680.89 1236.79 360.75 31621020.50 31705869.62 11985727
HLT selection/weights 35983.60 38283.89 303.20 8680.89 1236.79 360.75 4806886.00 4891735.12 7642710
Good PV 35982.09 38281.85 303.15 8670.75 1235.97 360.28 4758828.00 4843662.10 7586132
Data quality 35982.09 38281.85 303.15 8670.75 1235.97 360.28 4741392.00 4826226.10 7571606
Cosmics rejection 35981.82 38280.04 303.14 8670.01 1235.79 360.24 4446965.00 4531796.05 7568873
Dilepton 249.24 2280.10 20.32 547.55 216.89 94.09 165427.48 168835.67 873144
Z mass 55.03 511.12 4.27 130.50 177.10 85.24 165209.92 166173.18 783204
Z pT 48.19 446.97 3.50 89.94 140.82 69.65 165209.92 166008.99 168094
3rd lepton veto 48.04 444.32 3.47 89.74 82.87 69.60 164995.18 165733.22 167784
b-jet veto 15.30 70.19 3.40 87.88 79.20 67.83 155515.42 155839.21 156400

tW ttbar GluGluWW WW WZ ZZ PhotonJets GluGluHiggs VBFHiggs BG DATA Signal S/B Significance
250 3.64 ± 0.32 15.79 ± 0.28 0.88 ± 0.06 ± 0.44 9.17 ± 0.74 12.23 ± 0.33 14.93 ± 0.25 11.16 ± 1.97 8.28 ± 0.10 0.90 ± 0.02 67.81 ± 2.19 76 9.18 ± 0.10 0.14 ± 0.00 1.05 ± 0.00
300 0.61 ± 0.13 3.13 ± 0.12 0.14 ± 0.02 ± 0.07 1.16 ± 0.28 4.40 ± 0.20 7.34 ± 0.18 7.49 ± 1.80 6.26 ± 0.08 0.65 ± 0.01 24.28 ± 1.85 19 6.91 ± 0.08 0.28 ± 0.02 1.24 ± 0.01
350 0.17 ± 0.07 0.69 ± 0.06 0.04 ± 0.01 ± 0.02 0.01 ± 0.00 2.73 ± 0.16 5.14 ± 0.15 2.58 ± 1.16 6.38 ± 0.07 0.53 ± 0.01 11.35 ± 1.18 4 6.91 ± 0.07 0.61 ± 0.06 1.62 ± 0.02
400 0.11 ± 0.06 0.20 ± 0.03 0.01 ± 0.01 ± 0.00 0.00 ± 0.00 1.78 ± 0.13 3.81 ± 0.13 3.78 ± 1.55 5.54 ± 0.06 0.38 ± 0.01 9.68 ± 1.56 7 5.92 ± 0.06 0.61 ± 0.10 1.50 ± 0.04
450 0.01 ± 0.01 0.04 ± 0.02 0.00 ± 0.00 ± 0.00 0.01 ± 0.01 1.35 ± 0.11 2.94 ± 0.12 2.69 ± 1.35 3.83 ± 0.04 0.29 ± 0.00 7.06 ± 1.36 6 4.12 ± 0.04 0.58 ± 0.11 1.23 ± 0.04
500 0.00 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 ± 0.00 0.01 ± 0.01 0.92 ± 0.09 2.27 ± 0.10 2.69 ± 1.35 2.41 ± 0.03 0.22 ± 0.00 5.91 ± 1.35 6 2.63 ± 0.03 0.45 ± 0.10 0.90 ± 0.05
550 0.00 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 ± 0.00 0.01 ± 0.01 0.62 ± 0.08 1.74 ± 0.09 2.02 ± 1.17 1.48 ± 0.02 0.19 ± 0.00 4.40 ± 1.17 2 1.67 ± 0.02 0.38 ± 0.10 0.68 ± 0.05
600 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 ± 0.00 0.01 ± 0.01 0.40 ± 0.06 1.23 ± 0.08 0.00 ± 0.00 0.91 ± 0.01 0.12 ± 0.00 1.65 ± 0.10 2 1.03 ± 0.01 0.62 ± 0.04 0.63 ± 0.02

27 October

tW ttbar WW WZ ZZ DYToTauTau BG DATA
M = 250 11.31 ± 0.49 49.69 ± 0.42 38.37 ± 1.33 2.65 ± 0.13 0.02 ± 0.01 1.77 ± 1.77 103.81 ± 2.31 115

26 October

Does this even make any sense,

tW ttbar WZ WW ZZ DYToMuMu DYToTauTau BG DATA
0 jet 36.75 ± 0.88 246.06 ± 0.93 7.39 ± 0.22 90.42 ± 2.04 0.09 ± 0.01 3.64 ± 0.66 31.90 ± 7.52 416.26 ± 7.93 538
1 jet 67.55 ± 1.19 753.78 ± 1.62 0.43 ± 0.05 5.32 ± 0.50 0.01 ± 0.00 0.12 ± 0.12 3.54 ± 2.51 830.75 ± 3.25 760
2 jet 16.68 ± 0.59 535.30 ± 1.37 0.03 ± 0.01 0.37 ± 0.13 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 552.38 ± 1.50 471
3 jet 0.94 ± 0.14 41.17 ± 0.38 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 42.11 ± 0.40 40

23 October

Made some new ntuples (details) and new Higgs numbers. I uncovered a change in the name of the double electron analysis

19 October

More b-jet numbers

ZZ WZ WW ttbar tW DYToMuMu BG DATA
0 jet 21.62 ± 0.30 14.62 ± 0.35 8.66 ± 0.73 24.40 ± 0.33 4.37 ± 0.35 19.54 ± 1.76 93.21 ± 2.01 114
1 jet 0.73 ± 0.06 0.74 ± 0.08 0.55 ± 0.18 73.44 ± 0.58 7.72 ± 0.47 4.39 ± 0.81 87.57 ± 1.12 100
2 jet 0.03 ± 0.01 0.02 ± 0.01 0.00 ± 0.00 50.70 ± 0.48 1.57 ± 0.21 0.09 ± 0.08 52.41 ± 0.53 60
3 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 3.70 ± 0.13 0.14 ± 0.06 0.00 ± 0.00 3.85 ± 0.14 2

17 October

Test tables

tW ttbar WZ WW ZZ DYToMuMu BG DATA
single lepton 34224.46 33918.34 1176.33 8256.50 343.11 2533100.27 2611019.02 11985727
HLT selection/weights 34224.46 33918.34 1176.33 8256.50 343.11 2533100.27 2611019.02 7642710
Good PV 34223.03 33917.01 1175.55 8246.86 342.67 2529517.00 2607422.12 7586132
Data quality 34223.03 33917.01 1175.55 8246.86 342.67 2529517.00 2607422.12 7586132
Cosmics rejection 34222.77 33915.29 1175.38 8246.16 342.63 2526951.77 2604854.00 7583398
Dilepton 237.28 2170.58 207.13 523.74 91.07 823911.20 827141.00 874587
Z mass 52.37 486.48 169.05 124.81 82.49 748347.51 749262.72 784504
Z pT 45.86 425.35 134.39 86.01 67.59 156487.77 157246.97 168379
3rd lepton veto 45.71 423.45 79.08 85.81 67.54 156307.66 157009.26 168069
b-jet veto 15.20 69.54 75.62 84.04 65.83 148032.33 148342.56 156855

tW ttbar WZ WW ZZ DYToMuMu GluGlu VBF BG Signal S/B Significance DATA
250 3.59 ± 0.31 15.50 ± 0.27 11.76 ± 0.31 8.90 ± 0.72 14.47 ± 0.24 5.28 ± 0.86 7.91 ± 0.10 0.76 ± 0.01 59.50 ± 1.26 8.68 ± 0.10 0.15 ± 0.00 1.05 ± 0.00 78
300 0.62 ± 0.13 3.01 ± 0.12 4.23 ± 0.19 1.11 ± 0.27 7.23 ± 0.17 1.13 ± 0.40 5.96 ± 0.07 0.56 ± 0.01 17.33 ± 0.57 6.52 ± 0.07 0.38 ± 0.01 1.33 ± 0.01 20
350 0.16 ± 0.07 0.65 ± 0.05 2.61 ± 0.15 0.01 ± 0.00 5.12 ± 0.15 0.57 ± 0.32 6.07 ± 0.07 0.45 ± 0.01 9.12 ± 0.39 6.52 ± 0.07 0.72 ± 0.03 1.65 ± 0.01 4
400 0.11 ± 0.06 0.18 ± 0.03 1.71 ± 0.12 0.00 ± 0.00 3.84 ± 0.13 0.69 ± 0.33 5.27 ± 0.06 0.32 ± 0.01 6.52 ± 0.38 5.59 ± 0.06 0.86 ± 0.05 1.61 ± 0.02 7
450 0.01 ± 0.01 0.04 ± 0.01 1.29 ± 0.11 0.01 ± 0.01 3.00 ± 0.12 0.14 ± 0.09 3.64 ± 0.04 0.24 ± 0.00 4.49 ± 0.19 3.88 ± 0.04 0.86 ± 0.04 1.34 ± 0.01 6
500 0.00 ± 0.00 0.01 ± 0.01 0.87 ± 0.09 0.01 ± 0.01 2.35 ± 0.11 0.17 ± 0.10 2.29 ± 0.02 0.19 ± 0.00 3.43 ± 0.17 2.48 ± 0.02 0.72 ± 0.04 1.02 ± 0.01 6
550 0.00 ± 0.00 0.01 ± 0.01 0.59 ± 0.07 0.01 ± 0.01 1.83 ± 0.10 0.27 ± 0.19 1.41 ± 0.02 0.16 ± 0.00 2.71 ± 0.23 1.57 ± 0.02 0.58 ± 0.05 0.76 ± 0.02 2
600 0.00 ± 0.00 0.00 ± 0.00 0.38 ± 0.06 0.01 ± 0.01 1.33 ± 0.09 0.25 ± 0.19 0.86 ± 0.01 0.17 ± 0.00 1.98 ± 0.21 1.03 ± 0.01 0.52 ± 0.06 0.60 ± 0.02 2

13 October

For Radek

  • Muons

ZZ WZ WW ttbar tW DYToMuMu DATA BG
0 jet 22.33 ± 0.29 14.77 ± 0.33 11.24 ± 0.77 19.07 ± 0.28 4.19 ± 0.32 25.39 ± 1.87 107 97.00 ± 2.11
1 jet 0.84 ± 0.06 0.86 ± 0.08 0.85 ± 0.22 64.73 ± 0.51 7.70 ± 0.44 4.92 ± 0.79 87 79.89 ± 1.07
2 jet 0.02 ± 0.01 0.03 ± 0.01 0.01 ± 0.01 52.93 ± 0.46 1.87 ± 0.22 0.25 ± 0.18 63 55.11 ± 0.54
3 jet 0.01 ± 0.01 0.00 ± 0.00 0.00 ± 0.00 4.44 ± 0.13 0.15 ± 0.06 0.00 ± 0.00 4 4.60 ± 0.15

  • Electrons

ZZ WZ WW ttbar DYToEE tW DATA BG
0 jet 16.11 ± 0.25 11.22 ± 0.29 8.73 ± 0.68 12.98 ± 0.23 13.02 ± 1.30 3.08 ± 0.28 79 65.13 ± 1.55
1 jet 0.63 ± 0.05 0.56 ± 0.06 0.76 ± 0.21 45.30 ± 0.43 4.85 ± 0.84 5.38 ± 0.37 59 57.48 ± 1.04
2 jet 0.04 ± 0.01 0.03 ± 0.02 0.03 ± 0.02 36.81 ± 0.38 0.40 ± 0.18 1.61 ± 0.20 33 38.93 ± 0.47
3 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 2.98 ± 0.11 0.00 ± 0.00 0.17 ± 0.07 2 3.15 ± 0.13

11 October

For Radek

ZZ WZ WW ttbar tW DYToMuMu DATA BG
0 jet 18.62 ± 0.26 10.98 ± 0.28 6.90 ± 0.60 10.07 ± 0.20 2.50 ± 0.25 29.71 ± 1.96 61 78.78 ± 2.11
1 jet 1.05 ± 0.06 1.57 ± 0.11 0.82 ± 0.22 48.17 ± 0.44 6.30 ± 0.40 5.94 ± 0.85 68 63.85 ± 1.07
2 jet 0.07 ± 0.02 0.10 ± 0.03 0.01 ± 0.01 52.81 ± 0.46 2.05 ± 0.23 0.40 ± 0.24 66 55.43 ± 0.57
3 jet 0.01 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 6.55 ± 0.16 0.25 ± 0.08 0.00 ± 0.00 8 6.81 ± 0.18

tW ttbar WZ WW ZZ PhotonJets GluGlu VBF DATA BG Signal S/B Significance
250 1.49 ± 0.19 6.46 ± 0.18 5.57 ± 0.20 4.10 ± 0.46 7.85 ± 0.16 2.32 ± 0.48 3.69 ± 0.06 0.52 ± 0.01 22 27.80 ± 0.76 4.21 ± 0.06 0.15 ± 0.00 0.74 ± 0.00
300 0.59 ± 0.12 2.87 ± 0.12 3.17 ± 0.15 1.00 ± 0.22 6.21 ± 0.15 0.75 ± 0.28 4.23 ± 0.06 0.46 ± 0.01 10 14.59 ± 0.45 4.70 ± 0.06 0.32 ± 0.01 1.07 ± 0.01
350 0.19 ± 0.07 0.92 ± 0.07 1.77 ± 0.11 0.16 ± 0.10 4.47 ± 0.13 0.71 ± 0.34 4.93 ± 0.06 0.37 ± 0.01 7 8.23 ± 0.41 5.30 ± 0.06 0.64 ± 0.03 1.44 ± 0.01
400 0.11 ± 0.06 0.37 ± 0.04 1.67 ± 0.11 0.04 ± 0.02 4.39 ± 0.13 1.17 ± 0.39 4.26 ± 0.05 0.29 ± 0.00 4 7.76 ± 0.43 4.55 ± 0.05 0.59 ± 0.03 1.30 ± 0.01
500 0.06 ± 0.04 0.19 ± 0.03 1.66 ± 0.11 0.01 ± 0.01 4.49 ± 0.13 1.29 ± 0.41 1.72 ± 0.02 0.19 ± 0.00 5 7.70 ± 0.45 1.90 ± 0.02 0.25 ± 0.01 0.61 ± 0.01
600 0.03 ± 0.03 0.04 ± 0.02 1.13 ± 0.09 0.01 ± 0.01 3.51 ± 0.12 0.61 ± 0.27 0.77 ± 0.01 0.16 ± 0.00 3 5.33 ± 0.31 0.92 ± 0.01 0.17 ± 0.01 0.37 ± 0.01

6 October

HZZllnunu

  • Include Radek's numbers into table.
  • Add signal (gg plus VBF) over BG and significance to table.
  • Redo yields with Andrey's MET cut for masses 250-300.

  • test table:

tW ttbar WZ WW ZZ PhotonJets Signal VBF DATA BG S/B Significance
250 2.15 ± 0.23 7.87 ± 0.20 8.00 ± 0.24 7.31 ± 0.62 11.27 ± 0.20 6.31 ± 2.65 4.77 ± 0.07 0.60 ± 0.01 40 42.91 ± 2.76 0.13 ± 0.01 0.77 ± 0.01
300 0.87 ± 0.14 3.88 ± 0.14 5.13 ± 0.19 2.17 ± 0.34 8.97 ± 0.18 2.18 ± 0.65 5.09 ± 0.06 0.54 ± 0.01 21 23.21 ± 0.80 0.24 ± 0.01 1.05 ± 0.01
350 0.34 ± 0.09 1.42 ± 0.08 3.26 ± 0.15 0.67 ± 0.19 7.35 ± 0.17 1.59 ± 0.47 6.15 ± 0.07 0.46 ± 0.01 12 14.62 ± 0.57 0.45 ± 0.02 1.43 ± 0.01
400 0.13 ± 0.06 0.48 ± 0.05 2.38 ± 0.13 0.11 ± 0.08 5.77 ± 0.15 2.39 ± 0.55 4.72 ± 0.05 0.33 ± 0.01 5 11.26 ± 0.60 0.45 ± 0.02 1.25 ± 0.01
500 0.06 ± 0.04 0.19 ± 0.03 1.66 ± 0.11 0.01 ± 0.01 4.49 ± 0.13 1.98 ± 0.52 1.72 ± 0.02 0.19 ± 0.00 5 8.38 ± 0.56 0.23 ± 0.02 0.59 ± 0.01
600 0.03 ± 0.03 0.04 ± 0.02 1.13 ± 0.09 0.01 ± 0.01 3.51 ± 0.12 1.23 ± 0.44 0.77 ± 0.01 0.16 ± 0.00 3 5.95 ± 0.47 0.16 ± 0.01 0.35 ± 0.01

4 October

HZZllnunu

tW ttbar WZ WW ZZ DYToEE DATA BG
Generated 1624374 10340000 205000 211000 220000 29500000    
single lepton 28762.61 34486.83 965.56 6777.18 281.64 2055759.25 26884304 2127033.07
HLT selection/weights 28762.61 34486.83 965.56 6777.18 281.64 2055759.25 12042889 2127033.07
Good PV 28761.41 34485.48 964.55 6766.58 281.21 2051712.82 12041675 2122972.05
Data quality 28761.41 34485.48 964.55 6766.58 281.21 2051712.82 12041675 2122972.05
Cosmics rejection 28761.19 34483.73 964.41 6766.01 281.18 2051712.82 12041675 2122969.34
Dilepton 140.47 1545.34 149.04 292.21 51.59 455634.63 389670 457813.29
Z mass 31.71 342.84 107.12 71.01 47.04 418438.63 352422 419038.35
Z pT 28.07 302.62 86.29 49.74 38.68 89233.27 76989 89738.66
3rd lepton veto 25.68 267.52 45.10 49.48 38.55 88506.44 76206 88932.77
b-jet veto 5.81 23.23 40.11 46.66 36.65 78657.28 65315 78809.74

tW ttbar WZ WW ZZ DYToEE Signal VBF DATA BG
250 1.18 ± 0.17 4.94 ± 0.15 5.05 ± 0.19 5.50 ± 0.52 7.46 ± 0.16 4.61 ± 0.69 3.28 ± 0.06 0.42 ± 0.01 28 28.75 ± 0.92
300 0.72 ± 0.13 2.32 ± 0.11 3.60 ± 0.16 1.59 ± 0.27 6.00 ± 0.14 0.59 ± 0.25 3.60 ± 0.05 0.38 ± 0.01 22 14.81 ± 0.46
350 0.24 ± 0.07 0.92 ± 0.07 2.66 ± 0.14 0.37 ± 0.13 4.90 ± 0.13 0.47 ± 0.27 4.40 ± 0.05 0.34 ± 0.01 11 9.57 ± 0.37
400 0.07 ± 0.03 0.33 ± 0.04 1.79 ± 0.11 0.16 ± 0.08 4.06 ± 0.12 0.75 ± 0.33 3.42 ± 0.04 0.24 ± 0.00 10 7.15 ± 0.38
500 0.03 ± 0.02 0.13 ± 0.03 1.32 ± 0.10 0.11 ± 0.08 3.23 ± 0.11 0.16 ± 0.12 1.31 ± 0.02 0.14 ± 0.00 6 4.99 ± 0.21
600 0.03 ± 0.02 0.05 ± 0.02 0.85 ± 0.08 0.11 ± 0.08 2.55 ± 0.10 0.22 ± 0.14 0.57 ± 0.01 0.12 ± 0.00 2 3.81 ± 0.21

Muons

tW ttbar WZ WW ZZ DYToMuMu DATA BG
Generated 1624374 10340000 205000 211000 220000 29700000    
single lepton 28762.61 25716.18 965.56 6777.18 281.64 1561825.18 8302580 1624328.35
HLT selection/weights 28762.61 25716.18 965.56 6777.18 281.64 1561825.18 5036176 1624328.35
Good PV 28761.41 25715.15 964.55 6766.58 281.21 1559373.99 4967397 1621862.90
Data quality 28761.41 25715.15 964.55 6766.58 281.21 1559373.99 4967397 1621862.90
Cosmics rejection 28761.19 25713.79 964.41 6766.01 281.18 1557790.94 4965668 1620277.52
Dilepton 212.02 1705.18 181.88 460.80 79.90 532070.05 538940 534709.82
Z mass 46.75 381.88 148.32 109.68 72.36 483284.81 483438 484043.81
Z pT 40.93 333.80 117.97 75.66 59.32 101080.74 103847 101708.42
3rd lepton veto 37.58 294.63 63.13 74.99 59.05 100125.48 102631 100654.86
b-jet veto 9.19 25.87 55.27 71.44 56.25 89114.96 87835 89332.98

tW ttbar WZ WW ZZ DYToMuMu Signal VBF DATA BG
250 2.15 ± 0.23 5.70 ± 0.17 8.00 ± 0.24 7.31 ± 0.62 11.27 ± 0.20 4.94 ± 0.73 4.77 ± 0.07 0.60 ± 0.01 32 39.37 ± 1.04
300 0.87 ± 0.14 2.74 ± 0.11 5.13 ± 0.19 2.17 ± 0.34 8.97 ± 0.18 1.79 ± 0.46 5.09 ± 0.06 0.54 ± 0.01 17 21.68 ± 0.66
350 0.34 ± 0.09 0.97 ± 0.07 3.26 ± 0.15 0.67 ± 0.19 7.35 ± 0.17 1.48 ± 0.47 6.15 ± 0.07 0.46 ± 0.01 9 14.06 ± 0.57
400 0.13 ± 0.06 0.32 ± 0.04 2.38 ± 0.13 0.11 ± 0.08 5.77 ± 0.15 2.19 ± 0.54 4.72 ± 0.05 0.33 ± 0.01 4 10.91 ± 0.59
500 0.06 ± 0.04 0.12 ± 0.02 1.66 ± 0.11 0.01 ± 0.01 4.49 ± 0.13 1.22 ± 0.39 1.72 ± 0.02 0.19 ± 0.00 4 7.56 ± 0.43
600 0.03 ± 0.03 0.03 ± 0.01 1.13 ± 0.09 0.01 ± 0.01 3.51 ± 0.12 0.56 ± 0.26 0.77 ± 0.01 0.16 ± 0.00 3 5.27 ± 0.30

3 October

HZZllnunu

Need to put together numbers for tomorrow. I'm going to need:

  • Include efficiencies for triggers (with η-pT parametrization) and lepton reconstruction.
  • Calculate VBF signal cross-sections. Also, include samples in analysis.
  • Send Radek the input for his ttbar estimation, e.g., b-jet multiplicity with MET > 70 GeV.

Muons

DYToMuMu ZZ WZ WW ttbar tW DATA BG
0 jet 50.82 ± 2.37 20.47 ± 0.26 12.89 ± 0.29 10.45 ± 0.71 14.04 ± 0.25 3.20 ± 0.26 87 111.87 ± 2.53
1 jet 12.50 ± 1.19 1.10 ± 0.06 1.75 ± 0.11 0.97 ± 0.23 66.01 ± 0.55 7.11 ± 0.40 96 89.45 ± 1.39
2 jet 0.87 ± 0.35 0.06 ± 0.01 0.13 ± 0.03 0.01 ± 0.01 72.66 ± 0.57 2.27 ± 0.23 76 76.00 ± 0.71
3 jet 0.00 ± 0.00 0.01 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 8.43 ± 0.19 0.27 ± 0.08 8 8.70 ± 0.21

Electrons

DYToEE ZZ WZ WW ttbar tW DATA BG
0 jet 23.74 ± 1.44 12.78 ± 0.20 8.44 ± 0.22 7.37 ± 0.56 6.77 ± 0.15 1.96 ± 0.19 67 61.06 ± 1.59
1 jet 7.88 ± 0.90 0.72 ± 0.05 1.10 ± 0.08 0.69 ± 0.17 32.43 ± 0.32 4.54 ± 0.30 59 47.35 ± 1.02
2 jet 0.90 ± 0.30 0.03 ± 0.01 0.06 ± 0.02 0.02 ± 0.02 35.92 ± 0.34 1.68 ± 0.18 45 38.62 ± 0.49
3 jet 0.09 ± 0.09 0.01 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 4.27 ± 0.12 0.13 ± 0.05 4 4.51 ± 0.16

2 October

HZZllnunu

Need to put together numbers for tomorrow. I'm going to need:

  • both electrons and muons with the Z+jets background from the photon data
  • The ttbar sample should be rescaled so the the b-jet multiplicity matches for the 2-jet bin.
  • Include efficiencies for triggers (with η-pT parametrization) and lepton reconstruction.
  • Update selection information on twiki and submit revisions of the code for Anton.

BRM

Look at some recent TDC data, and fit (linear?) the albedo before some non-colliding bunches.

30 September

More test tables (ZJets from data)

tW ttbar WZ WW ZZ ZJets DATA BG
single lepton 28762.61 27841.15 965.56 6777.18 281.64 33443762.00 9883487 33508390.14
HLT selection/weights 28762.61 27841.15 965.56 6777.18 281.64 1801594.00 6158706 1866222.14
Good PV 28761.41 27840.06 964.55 6766.58 281.21 1798165.00 6079793 1862778.81
Data quality 28761.41 27840.06 964.55 6766.58 281.21 1798165.00 6079793 1862778.81
Cosmics rejection 28761.19 27838.64 964.41 6766.01 281.18 1798165.00 6077699 1862776.44
Dilepton 192.58 1720.71 165.26 419.14 72.58 186820.80 663854 189391.07
Z mass 42.46 385.27 134.76 99.77 65.73 186603.28 595395 187331.26
Z pT 37.16 336.86 107.13 68.77 53.86 131129.88 127819 131733.66
3rd lepton veto 34.13 297.47 57.33 68.17 53.61 118121.15 126314 118631.85
b-jet veto 10.19 40.33 54.31 66.17 51.96 105934.36 115776 106157.32

tW ttbar WZ WW ZZ ZJets (Data) Signal DATA BG
HZZ250 2.02 ± 0.21 7.03 ± 0.16 6.78 ± 0.21 6.05 ± 0.53 9.76 ± 0.17 16.69 ± 9.85 4.43 ± 0.07 40 48.33 ± 9.87
HZZ300 0.89 ± 0.14 3.75 ± 0.12 4.73 ± 0.18 1.83 ± 0.30 7.98 ± 0.16 3.99 ± 1.48 4.96 ± 0.06 19 23.17 ± 1.54
HZZ350 0.36 ± 0.09 1.46 ± 0.07 3.21 ± 0.15 0.49 ± 0.16 6.69 ± 0.15 2.54 ± 0.98 6.13 ± 0.07 13 14.74 ± 1.02
HZZ400 0.17 ± 0.06 0.52 ± 0.04 2.40 ± 0.13 0.10 ± 0.07 5.34 ± 0.14 3.64 ± 1.22 4.88 ± 0.05 5 12.17 ± 1.24
HZZ450 0.10 ± 0.05 0.20 ± 0.03 1.69 ± 0.11 0.01 ± 0.01 4.14 ± 0.12 2.58 ± 1.06 3.70 ± 0.04 5 8.72 ± 1.08
HZZ500 0.03 ± 0.03 0.06 ± 0.02 1.16 ± 0.09 0.01 ± 0.01 3.25 ± 0.11 0.92 ± 0.66 1.79 ± 0.02 3 5.43 ± 0.67
HZZ550 0.00 ± 0.00 0.02 ± 0.01 0.85 ± 0.08 0.01 ± 0.01 2.56 ± 0.10 0.92 ± 0.66 1.38 ± 0.01 4 4.38 ± 0.67
HZZ600 0.00 ± 0.00 0.01 ± 0.01 0.59 ± 0.06 0.01 ± 0.01 2.05 ± 0.09 0.92 ± 0.66 0.88 ± 0.01 4 3.59 ± 0.66

tW ttbar WZ WW ZZ ZJets (MC) Signal DATA BG
HZZ250 2.02 ± 0.21 7.03 ± 0.16 6.78 ± 0.21 6.05 ± 0.53 9.76 ± 0.17 4.00 ± 0.66 4.43 ± 0.07 40 35.64 ± 0.93
HZZ300 0.89 ± 0.14 3.75 ± 0.12 4.73 ± 0.18 1.83 ± 0.30 7.98 ± 0.16 1.91 ± 0.47 4.96 ± 0.06 19 21.10 ± 0.63
HZZ350 0.36 ± 0.09 1.46 ± 0.07 3.21 ± 0.15 0.49 ± 0.16 6.69 ± 0.15 1.95 ± 0.51 6.13 ± 0.07 13 14.16 ± 0.58
HZZ400 0.17 ± 0.06 0.52 ± 0.04 2.40 ± 0.13 0.10 ± 0.07 5.34 ± 0.14 2.08 ± 0.50 4.88 ± 0.05 5 10.61 ± 0.55
HZZ450 0.10 ± 0.05 0.20 ± 0.03 1.69 ± 0.11 0.01 ± 0.01 4.14 ± 0.12 1.17 ± 0.37 3.70 ± 0.04 5 7.31 ± 0.41
HZZ500 0.03 ± 0.03 0.06 ± 0.02 1.16 ± 0.09 0.01 ± 0.01 3.25 ± 0.11 0.55 ± 0.24 1.79 ± 0.02 3 5.07 ± 0.28
HZZ550 0.00 ± 0.00 0.02 ± 0.01 0.85 ± 0.08 0.01 ± 0.01 2.56 ± 0.10 0.55 ± 0.24 1.38 ± 0.01 4 4.01 ± 0.27
HZZ600 0.00 ± 0.00 0.01 ± 0.01 0.59 ± 0.06 0.01 ± 0.01 2.05 ± 0.09 0.28 ± 0.16 0.88 ± 0.01 4 2.94 ± 0.20

29 September

BCM1F signal transfer

Andrey and I tested the lines and determined the mapping for the ethernet connection between the BCM1F rack and the gated rack. We found that one of the ports (second from the top) in the LVNIM send module has 2 faulty channels (5 and 6). Tomorrow I will check the third port. This will require that I remove the card and set the port to send mode if it is not already.

22 September

Test tables with signal MC

tW ttbar WZ WW ZZ ZMuMuJets DATA BG
Generated 1624374 10340000 205000 211000 220000 29700000    
single lepton 31958.46 30934.61 1072.85 7530.19 312.93 2310268.35 9883487 2382077.39
HLT selection/weights 31958.46 30934.61 1072.85 7530.19 312.93 2310268.35 6158706 2382077.39
Good PV 31957.12 30933.40 1071.72 7518.43 312.46 2306661.37 6079793 2378454.50
Data quality 31957.12 30933.40 1071.72 7518.43 312.46 2306661.37 6079793 2378454.50
Cosmics rejection 31956.88 30931.83 1071.57 7517.79 312.42 2304321.90 6077699 2376112.39
Dilepton 213.98 1911.90 183.62 465.71 80.64 731131.58 663854 733987.43
Z mass 47.18 428.07 149.73 110.86 73.03 663994.73 595395 664803.60
Z pT 41.29 374.29 119.03 76.41 59.84 138777.23 127819 139448.10
3rd lepton veto 37.92 330.52 63.70 75.74 59.56 137472.12 126314 138039.55
b-jet veto 11.32 44.81 60.35 73.52 57.73 128264.96 115776 128512.69

tW ttbar WZ WW ZZ ZMuMuJets Signal DATA BG
M250 2.64 ± 0.25 9.19 ± 0.19 7.81 ± 0.24 6.42 ± 0.57 10.87 ± 0.20 1.40 ± 0.42 5.33 ± 0.08 34 38.34 ± 0.84
M300 1.03 ± 0.16 4.25 ± 0.13 5.23 ± 0.20 1.58 ± 0.29 8.82 ± 0.18 1.01 ± 0.33 5.81 ± 0.07 23 21.91 ± 0.55
M350 0.42 ± 0.10 1.42 ± 0.08 3.47 ± 0.16 0.37 ± 0.14 7.32 ± 0.17 1.02 ± 0.40 6.99 ± 0.07 13 14.02 ± 0.50
M400 0.20 ± 0.07 0.50 ± 0.05 2.47 ± 0.13 0.03 ± 0.02 5.75 ± 0.15 1.88 ± 0.50 5.38 ± 0.06 5 10.82 ± 0.55
M450 0.08 ± 0.04 0.17 ± 0.03 1.77 ± 0.12 0.01 ± 0.01 4.45 ± 0.14 0.74 ± 0.29 4.08 ± 0.04 5 7.21 ± 0.35
M500 0.03 ± 0.03 0.05 ± 0.01 1.16 ± 0.09 0.01 ± 0.01 3.34 ± 0.12 0.59 ± 0.25 1.97 ± 0.02 4 5.19 ± 0.30
M550 0.00 ± 0.00 0.01 ± 0.01 0.80 ± 0.08 0.01 ± 0.01 2.71 ± 0.11 0.24 ± 0.17 1.52 ± 0.02 4 3.77 ± 0.21
M600 0.00 ± 0.00 0.01 ± 0.01 0.58 ± 0.07 0.01 ± 0.01 2.16 ± 0.10 0.08 ± 0.05 0.97 ± 0.01 4 2.84 ± 0.13

21 September

Numbers for Radek

b-jet pT > 30 GeV

bJetMultPostVeto.png

ZMuMuJets ZZ WZ WW ttbar tW DATA BG
0 jet 24.57 ± 1.86 22.12 ± 0.29 14.54 ± 0.33 10.39 ± 0.73 23.17 ± 0.31 4.43 ± 0.33 110 99.22 ± 2.09
1 jet 4.16 ± 0.72 0.71 ± 0.05 0.68 ± 0.07 0.47 ± 0.15 67.28 ± 0.52 7.68 ± 0.44 82 80.97 ± 1.01
2 jet 0.09 ± 0.09 0.02 ± 0.01 0.02 ± 0.01 0.00 ± 0.00 45.48 ± 0.43 1.34 ± 0.18 58 46.95 ± 0.48
3 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 3.17 ± 0.11 0.11 ± 0.05 2 3.28 ± 0.12

b-jet pT > 20 GeV

ZMuMuJets ZZ WZ WW ttbar tW DATA BG
0 jet 24.14 ± 1.84 22.00 ± 0.29 14.34 ± 0.32 10.10 ± 0.71 18.81 ± 0.28 3.94 ± 0.31 102 93.34 ± 2.06
1 jet 4.58 ± 0.76 0.83 ± 0.06 0.86 ± 0.08 0.76 ± 0.20 63.74 ± 0.51 7.59 ± 0.44 86 78.35 ± 1.04
2 jet 0.09 ± 0.09 0.02 ± 0.01 0.03 ± 0.01 0.01 ± 0.01 52.05 ± 0.46 1.89 ± 0.22 60 54.08 ± 0.52
3 jet 0.00 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 0.00 ± 0.00 4.39 ± 0.13 0.12 ± 0.05 4 4.52 ± 0.14

15 September

Run meeting

  • Implemented albedo corrections for BKGD 2.

Test numbers for comparison w/ Andrey

ZJets ZZ WZ WW ttbar tW DATA Background
Generated 29500000 220000 205000 220000 10340000 795000    
Single PV + single lepton 2325924.99 312.93 1072.85 7221.45 30934.61 10106.78 9883487 2375573.62
HLT selection/weights 2325924.99 312.93 1072.85 7221.45 30934.61 10106.78 6158706 2375573.62
Good PV 2322293.57 312.46 1071.72 7210.17 30933.40 10106.38 6079793 2371927.70
Data quality 2322293.57 312.46 1071.72 7210.17 30933.40 10106.38 6068779 2371927.70
Cosmics rejection 2319938.25 312.42 1071.57 7209.56 30931.83 10106.25 6066685 2369569.87
Dilepton 735607.19 80.96 184.66 444.97 1907.18 105.08 662856 738330.03
Z mass 668003.38 73.29 150.62 107.19 426.79 23.94 594497 668785.21
Z pT 139611.98 60.07 119.62 73.84 373.01 20.85 127624 140259.37
3rd lepton veto 138206.66 59.73 64.38 73.07 325.87 19.44 126036 138749.15
b-jet veto 128986.18 57.91 60.95 70.98 44.01 6.15 115538 129226.17
M250 5.27 ± 0.88 11.64 ± 0.21 8.45 ± 0.27 7.09 ± 0.64 9.35 ± 0.21 1.00 ± 0.17 48 42.81
M300 1.92 ± 0.50 9.19 ± 0.20 5.50 ± 0.21 2.07 ± 0.35 4.54 ± 0.14 0.42 ± 0.11 24 23.64
M350 1.97 ± 0.55 7.58 ± 0.18 3.70 ± 0.18 0.61 ± 0.20 1.65 ± 0.09 0.16 ± 0.07 14 15.69
M400 2.22 ± 0.56 5.99 ± 0.16 2.71 ± 0.15 0.09 ± 0.08 0.58 ± 0.05 0.07 ± 0.04 5 11.66
M450 1.26 ± 0.41 4.66 ± 0.15 1.95 ± 0.13 0.00 ± 0.00 0.23 ± 0.03 0.00 ± 0.00 5 8.10
M500 0.55 ± 0.24 3.62 ± 0.13 1.36 ± 0.11 0.00 ± 0.00 0.08 ± 0.02 0.00 ± 0.00 3 5.60
M550 0.55 ± 0.24 2.85 ± 0.12 1.01 ± 0.09 0.00 ± 0.00 0.03 ± 0.01 0.00 ± 0.00 4 4.45
M600 0.30 ± 0.18 2.29 ± 0.11 0.70 ± 0.08 0.00 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 4 3.30

14 September

  • Preselection + MET > 70 - b-jet veto

ZJets ZZ WZ WW ttbar tW DATA BG
Generated 29500000 220000 205000 220000 10340000 795000    
Single lepton 2325924.99 312.93 1072.85 7221.45 30934.61 10106.78 9547834 2375573.62
HLT selection/weights 2325924.99 312.93 1072.85 7221.45 30934.61 10106.78 5901429 2375573.62
Good PV 2322293.57 312.46 1071.72 7210.17 30933.40 10106.38 5846314 2371927.70
Data quality 2322293.57 312.46 1071.72 7210.17 30933.40 10106.38 5834995 2371927.70
Cosmics 2319938.25 312.42 1071.57 7209.56 30931.83 10106.25 5833026 2369569.87
Dilepton 735607.19 80.96 184.66 444.97 1907.18 105.08 635413 738330.03
Z mass 668003.38 73.29 150.62 107.19 426.79 23.94 569957 668785.21
Z pT 139611.98 60.07 119.62 73.84 373.01 20.85 122268 140259.37
3rd lepton veto 138206.66 59.73 64.38 73.07 325.87 19.44 120758 138749.15
MET 56.91 24.14 16.60 12.18 143.96 6.94 232 260.72

29 August

  • With mass cuts

ZMuMuJets ZZ WZ WW ttbar tW DATA Background
Generated 29500000 220000 205000 220000 10340000 795000  
Single PV + single lepton 2320104.93 ± 510.31 312.42 ± 0.78 1071.54 ± 2.49 7051.31 ± 16.95 30852.62 ± 9.95 10104.50 ± 14.10 9883487 2369497.31 ± 510.89
HLT selection/weights 2320104.93 ± 510.31 312.42 ± 0.78 1071.54 ± 2.49 7051.31 ± 16.95 30852.62 ± 9.95 10104.50 ± 14.10 6158706 2369497.31 ± 510.89
Data quality 2306865.33 ± 508.85 311.20 ± 0.78 1068.52 ± 2.48 7022.79 ± 16.92 30847.54 ± 9.95 10103.47 ± 14.10 6136701 2356218.85 ± 509.43
Dilepton 740846.47 ± 345.31 81.69 ± 0.56 186.11 ± 1.25 438.56 ± 5.07 1941.00 ± 2.99 106.91 ± 1.74 681195 743600.74 ± 345.36
Z mass 669988.96 ± 328.36 73.64 ± 0.53 151.12 ± 1.12 106.05 ± 2.50 435.32 ± 1.41 24.34 ± 0.83 608393 670779.43 ± 328.38
Z pT 140156.61 ± 150.18 60.34 ± 0.49 119.91 ± 1.00 72.94 ± 2.07 380.07 ± 1.32 21.20 ± 0.78 130770 140811.07 ± 150.21
3rd lepton veto 131754.34 ± 145.60 58.66 ± 0.48 63.49 ± 0.73 70.65 ± 2.04 292.33 ± 1.16 18.13 ± 0.72 122263 132257.60 ± 145.62
b-jet veto 123176.16 ± 140.77 56.99 ± 0.47 60.34 ± 0.71 68.73 ± 2.01 41.57 ± 0.44 5.75 ± 0.41 112505 123409.54 ± 140.79
M250 1.42 ± 0.48 10.75 ± 0.21 8.01 ± 0.26 5.97 ± 0.57 8.57 ± 0.20 0.85 ± 0.15 31 35.56 ± 0.86
M300 0.44 ± 0.16 8.66 ± 0.19 5.20 ± 0.21 1.33 ± 0.27 3.99 ± 0.13 0.35 ± 0.10 24 19.96 ± 0.45
M350 0.82 ± 0.36 7.15 ± 0.18 3.48 ± 0.17 0.33 ± 0.14 1.36 ± 0.08 0.11 ± 0.06 13 13.25 ± 0.47
M400 0.10 ± 0.10 3.79 ± 0.13 1.45 ± 0.11 0.00 ± 0.00 0.05 ± 0.01 0.00 ± 0.00 2 5.40 ± 0.20
M450 0.10 ± 0.10 2.89 ± 0.12 0.94 ± 0.09 0.00 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 1 3.95 ± 0.18
M500 0.10 ± 0.10 2.16 ± 0.10 0.58 ± 0.07 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 2 2.85 ± 0.16
M550 0.00 ± 0.00 1.78 ± 0.09 0.42 ± 0.06 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 2 2.20 ± 0.11
M600 0.00 ± 0.00 1.41 ± 0.08 0.28 ± 0.05 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 2 1.68 ± 0.10

26 August

  • preselction numbers

ZMuMuJets ZZ WZ WW ttbar tW DATA Background
Generated 29500000 220000 205000 220000 10340000 795000  
Single PV + single lepton 2320104.93 ± 510.31 312.42 ± 0.78 1071.54 ± 2.49 7051.31 ± 16.95 30852.62 ± 9.95 10104.50 ± 14.10 9883487 2369497.31 ± 510.89
HLT selection/weights 2320104.93 ± 510.31 312.42 ± 0.78 1071.54 ± 2.49 7051.31 ± 16.95 30852.62 ± 9.95 10104.50 ± 14.10 6158706 2369497.31 ± 510.89
Data quality 2306865.33 ± 508.85 311.20 ± 0.78 1068.52 ± 2.48 7022.79 ± 16.92 30847.54 ± 9.95 10103.47 ± 14.10 6136701 2356218.85 ± 509.43
Dilepton 740846.47 ± 345.31 81.69 ± 0.56 186.11 ± 1.25 438.56 ± 5.07 1941.00 ± 2.99 106.91 ± 1.74 681195 743600.74 ± 345.36
Z mass 669988.96 ± 328.36 73.64 ± 0.53 151.12 ± 1.12 106.05 ± 2.50 435.32 ± 1.41 24.34 ± 0.83 608393 670779.43 ± 328.38
Z pT 140156.61 ± 150.18 60.34 ± 0.49 119.91 ± 1.00 72.94 ± 2.07 380.07 ± 1.32 21.20 ± 0.78 130770 140811.07 ± 150.21
3rd lepton veto 131754.34 ± 145.60 58.66 ± 0.48 63.49 ± 0.73 70.65 ± 2.04 292.33 ± 1.16 18.13 ± 0.72 122263 132257.60 ± 145.62
b-jet veto 123176.16 ± 140.77 56.99 ± 0.47 60.34 ± 0.71 68.73 ± 2.01 41.57 ± 0.44 5.75 ± 0.41 112505 123409.54 ± 140.79

  • Multiplicity for b-jets for M300 higgs

Z+Jets ZZ WZ WW ttbar tW DATA Background
0 jet 0.44 ± 0.16 8.83 ± 0.19 5.32 ± 0.21 1.47 ± 0.29 4.06 ± 0.14 0.35 ± 0.10 24 20.45 ± 0.47
1 jet 0.13 ± 0.13 0.18 ± 0.03 0.18 ± 0.04 0.11 ± 0.08 11.90 ± 0.24 0.75 ± 0.15 15 13.26 ± 0.32
2 jet 0.00 ± 0.00 0.01 ± 0.01 0.01 ± 0.01 0.00 ± 0.00 8.48 ± 0.20 0.14 ± 0.07 8 8.64 ± 0.21
3 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.51 ± 0.05 0.04 ± 0.03 1 0.55 ± 0.05

25 August

Z+Jets ZZ WZ WW ttbar tW Background
0 jet 0.06 ± 0.06 0.37 ± 0.02 0.35 ± 0.03 0.16 ± 0.06 0.29 ± 0.02 0.04 ± 0.02 1.27 ± 0.09
1 jet 0.00 ± 0.00 0.01 ± 0.00 0.01 ± 0.01 0.00 ± 0.00 0.76 ± 0.03 0.05 ± 0.02 0.83 ± 0.04
2 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.66 ± 0.03 0.00 ± 0.00 0.66 ± 0.03
3 jet 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.00 ± 0.00 0.03 ± 0.01 0.00 ± 0.00 0.03 ± 0.01

24 August

Just a test...

Z+Jets ZZ WZ WW ttbar tW Signal Data Background
Generated 29500000 220000 205000 220000 10340000 795000 96990  
Single PV + single lepton 2290731.00 ± 1513.52 18659.00 ± 136.60 21083.00 ± 145.20 19604.00 ± 140.01 1111636.00 ± 1054.34 59888.00 ± 244.72 10329.00 ± 101.63 169186 0.0
HLT selection/weights 828189.00 ± 910.05 3950.00 ± 62.85 4040.00 ± 63.56 1393.00 ± 37.32 79862.00 ± 282.60 719.00 ± 26.81 2830.00 ± 53.20 18306 31309607.0
Data quality 748913.00 ± 865.40 3551.00 ± 59.59 3319.00 ± 57.61 349.00 ± 18.68 17736.00 ± 133.18 166.00 ± 12.88 2622.00 ± 51.21 16350 31309607.0
Dilepton 156366.00 ± 395.43 2870.00 ± 53.57 2637.00 ± 51.35 250.00 ± 15.81 15483.00 ± 124.43 150.00 ± 12.25 2558.00 ± 50.58 3522 3521601.0
Z mass 155164.00 ± 393.91 2855.00 ± 53.43 1570.00 ± 39.62 249.00 ± 15.78 14207.00 ± 119.19 145.00 ± 12.04 2531.00 ± 50.31 3486 918153.0
Z pT 155164.00 ± 393.91 2855.00 ± 53.43 1570.00 ± 39.62 249.00 ± 15.78 14207.00 ± 119.19 145.00 ± 12.04 2531.00 ± 50.31 3486 774034.0
3rd lepton veto 134652.00 ± 366.95 2784.00 ± 52.76 1482.00 ± 38.50 242.00 ± 15.56 12343.00 ± 111.10 133.00 ± 11.53 2382.00 ± 48.81 3083 177756.0
b-jet veto 6.00 ± 2.45 652.00 ± 25.53 201.00 ± 14.18 18.00 ± 4.24 3597.00 ± 59.97 32.00 ± 5.66 1722.00 ± 41.50 3 174190.0
Δφ(MET, jet) 1.00 ± 1.00 388.00 ± 19.70 126.00 ± 11.22 8.00 ± 2.83 1075.00 ± 32.79 9.00 ± 3.00 1543.00 ± 39.28 0 174190.0
MET 1.00 ± 1.00 388.00 ± 19.70 126.00 ± 11.22 8.00 ± 2.83 1075.00 ± 32.79 9.00 ± 3.00 1543.00 ± 39.28 0 151636.0
Final (raw) 1 388 126 8 1075 9 1543 0  

17 August

Z+Jets ZZ WZ WW ttbar tW Signal Data Background
Single PV + single lepton 1273705.61 171.80 589.19 3964.91 16996.50 5552.97 30.42 7954347 1300981.0
HLT selection/weights 1273705.61 171.80 589.19 3964.91 16996.50 5552.97 30.42 4283530 1300981.0
Data quality 1267464.49 171.27 587.80 3951.68 16993.73 5552.44 30.42 4276240 1294721.4
Dilepton 408159.32 44.95 102.41 246.53 1068.63 58.76 7.36 446227 409680.6
Z mass 369161.17 40.53 83.17 59.53 239.75 13.37 6.87 398618 369597.5
3rd lepton veto 368313.42 40.38 47.87 59.25 212.22 12.49 6.80 397479 368685.6
b-jet veto 368313.42 40.38 47.87 59.25 212.22 12.49 6.80 397479 368685.6
Z pT 76339.86 33.09 37.73 40.75 185.34 10.89 6.64 84616 76647.7
Δφ(MET, jet) 68056.65 32.15 35.39 39.70 160.01 9.87 6.19 74681 68333.8
MET 1.88 9.03 5.36 2.38 49.31 2.30 4.46 86 70.3
MT 0.49 5.00 3.18 0.90 15.42 0.78 3.85 35 25.8

12 August

Changes to talk

  • Underline name in title page (DONE)
  • only my name in footer (uh...pretty much)
  • specify which cuts are affected by our changes selection slides (uh...pretty much)
  • add jet-vertex association for b-jets (uh...)
  • remove unchanged parameters from 3rd lepton veto slide (DONE)
  • include slide before PU with changes to yields for each change; add S/B (DONE)
  • add muon isolation slide; point is, too sensitive to PU. (DONE)
  • add nJet < 2 cut for final cut (uh...)

11 August

Z+Jets ZZ WZ WW ttbar tW Signal Data Background
μμ 0.63 4.99 3.22 1.14 3.00 0.34 3.74 13 13.3
ee 0.23 3.87 2.65 0.87 2.18 0.15 3.01 19 10.0

9 August


| Two oppositely charged leptons:    |	22693	|24220.4	|
| Z mass window:                     |	20968	|22380.5	|
| Third lepton veto:                 |	20739	|22136.5	|
| b-jet veto:                        |	18842	|20049.1	|
| Z Qt :                             |	18618	|19807.6	|
| deltaPhi(MET, jet):                |	18618	|19807.6	|
| MT:                                |	14339	|15260.7	|
| MET:                               |	13480	|14358.4	|

Friday update

  • reproduce EPS results with 300 GeV Higgs

Changes

  • lower threshold on b-jets
  • tighter third lepton vetor
  • jet-vertex association
  • change to muon isolation
  • electron endcap H/E
  • MET corrections for PU

ntuples

For next ntuple production...

Data

For each dataset there are three periods: May 10th rereco, prompt v4, and the Aug 5th rereco. We may also include prompt v6.

MC

Signal

For the signal we have 3 types of samples (ZZ*, VBFZZ, and WW) and 8 mass bins (140*, 200*, 250, 300*, 350, 400*, 450, 500, 550, 600). These need to be run over the grid.

3 August

Test yields for M=300 Higgs

MUONS

Z+Jets ZZ WZ WW ttbar tW Signal Data Background
Generated 29500000 220000 205000 220000 10340000 795000 96990 -- 41376990
Ntuple 20653393.00 161732.00 185585.00 177114.00 9631195.00 513701.00 84130.00 7954347 31322720.0
HLT selection/weights 1323652.61 126.04 589.19 3964.91 16996.50 7451.27 18.56 4277626 1352780.5
data quality bits 1322626.36 131.41 628.02 4162.60 17656.30 7773.52 19.49 4276240 1352978.2
Dilepton 514641.44 27.85 119.41 291.01 1264.69 91.92 5.27 446119 516436.3
Z mass 440521.28 25.14 96.95 69.45 284.27 20.81 4.91 398565 441017.9
3rd lepton veto 424153.34 24.45 65.00 67.26 205.18 17.39 4.60 390175 424532.6
b-jet veto 417632.82 24.20 63.74 66.73 54.61 8.04 4.48 382104 417850.1
Z pT 64527.93 14.39 36.25 29.21 38.58 5.90 4.21 72838 64652.3
Δφ(MET, jet) 57160.75 14.39 36.25 29.21 38.58 5.90 4.21 62897 57285.1
MT 6708.60 3.40 6.83 2.04 6.44 0.77 3.10 247 6728.1
MET 2.82 2.86 4.72 1.31 4.52 0.47 2.68 10 16.7
MET/QT 2.82 2.86 4.72 1.31 4.52 0.47 2.68 10 16.7
MET + QT 2.82 2.86 4.72 1.31 4.52 0.47 2.68 10 16.7
nJets < 2 2.69 2.56 3.75 1.16 1.02 0.16 1.91 3 11.3

Comparing to the EPS selection

Muon selection
| CUT DESCRIPTION                    |  May10   |   Prompt  |
| Initial number of events:          |  2878369 |   5075978 |
| Pass HLT selection:                |  957013  |   3326517 |
| Cosmics veto:                      |  955777  |   3321849 |
| Two oppositely charged leptons:    |  92949   |   536503  |
| Z mass window:                     |  61761   |   373238  |
| lepton pt cut :                 |  56372   |   341701  |
| b-jet veto:                        |  55406   |   335918  |
| Z Qt :                             |  11178   |   68005   |
| deltaPhi(MET, jet):                |  10341   |   63068   |
| MT:                                |  50      |   252     |
| MET:                               |  2       |   11      |

| CUT DESCRIPTION                    |  May10   |   Prompt   |
| Initial number of events:          |  7001065 |   16172294 |
| Pass HLT selection:                |  1785871 |   8294984  |
| Cosmics veto:                      |  1785310 |   8292508  |
| Two oppositely charged leptons:    |  64449   |   301050   |
| Z mass window:                     |  53376   |   249979   |
| lepton pt cut  :                 |  50315   |   235409   |
| b-jet veto:                        |  49423   |   231450   |
| Z Qt :                             |  10222   |   47737    |
| deltaPhi(MET, jet):                |  9478    |   44110    |
| MT:                                |  47      |   249      |
| MET:                               |  3       |   17       |


1 August

Some test yields

Trying to recreate the results for the 300 GeV Higgs.

Z+Jets ZZ WZ WW ttbar tW Signal (M=200 GeV) Data (DoubleMu) Background
Initial 1225733.976 802.779 336.140 2549.545 16132.754 7451.271 55.711 7954347 1253006.5
HLT selection 1225733.976 802.779 336.140 2549.545 16132.754 7451.271 55.711 4301175 1253006.5
data quality bits 1281581.022 836.992 358.290 2676.665 16759.023 7773.519 59.170 4294302 1309985.5
Dilepton 461733.431 177.385 68.125 187.128 1200.419 91.917 15.249 443289 463458.4
Z mass 417631.602 160.159 55.309 44.659 269.821 20.811 14.525 395688 418182.4
3rd lepton veto (> 10 GeV) 410220.868 155.722 37.083 43.248 194.749 17.395 13.815 387675 410669.1
b-jet veto (> 20 GeV) 406970.796 154.166 36.365 42.907 51.833 8.037 13.521 384006 407264.1
Z pT (> 25 GeV) 406970.796 154.166 36.365 42.907 51.833 8.037 13.521 384006 407264.1
Δφ(MET, jet) 406970.796 154.166 36.365 42.907 51.833 8.037 13.521 384006 407264.1
MT > 0 GeV 406970.796 154.166 36.365 42.907 51.833 8.037 13.521 384006 407264.1
MET 406970.796 154.166 36.365 42.907 51.833 8.037 13.521 384006 407264.1
MET/QT 406970.010 154.166 36.365 42.907 51.833 8.037 13.521 384005 407263.3
MET + QT 406970.010 154.166 36.365 42.907 51.833 8.037 13.521 384005 407263.3
nJets < 2 406970.010 154.166 36.365 42.907 51.833 8.037 13.521 384005 407263.3

Analysis code instructions

To run the analysis code check out the code along with the necessary libraries,

cvs co UserCode/NatesAnalyzer
cd UserCode/NatesAnalyzer/HiggsAnalyzer

In this directory are several scripts which, depending on how you run, perform different functions. The core of the package is the higgsAnalyzer.C and higgsAnalyzer.h. Various parameters can be specified (trigger selection, lepton selection, analysis cuts, etc.) at the top of higgsAnalyzer.C. The two main selection choices are "electron" and "muon" (there are photon selections under development).

Now you need to specify a dataset to run over. If you are running on lpc, it should run out of the box; otherwise, you must specify the root file you want to run over in runLocal.C. This should be done with

fChain->Add()

Once this is done, just type the command,

root -l runLocal.C

Altenatively, you can run using the higgsLocal.csh script. In this case you are using the higgsAnalyzer_Template.C file as the body of the code and the datasets are specified in a text file in the sourceFiles directory. To run,

./higgsLocal.csh SUFFIX TRIGGERS DATASET SELECTION

Basically, you are now specifying some of the more important parameters that were specified in the code in the command line. The dataset field should contain the name of the text file that lists the paths to the datasets. For running over many datasets at once there is a script that will submit multiple jobs to the lpc batch system. This is done with the multiBatch.csh script,

./multiBatch.csh BATCHCONFIG

The one field you specify should let the script know about a text file that has a list of the arguments you would give to the higgsLocal.csh script.

25 July

Using the higgs analyzer on lxplus

Log in to lxplus and then, if it is available, go to your scratch area. We need to set up a release of CMSSW

cd scratch0
cmsrel CMSSW_4_2_4
cd CMSSW_4_2_4/src
cmsenv

The checkout the necessary code from CVS,

cvs co -r V00-00-01 /UserCode/NatesAnalyzer/HiggsAnalyzer
cvs co /UserCode/NatesAnalyzer/src
cd /UserCode/NatesAnalyzer/HiggsAnalyzer

The important files here are higgsAnalyzer.C, higgsAnalyzer.h, and runLocal.C. In the analyzer, you make plots, apply selections, and do any calculations that need to be done. You specify what dataset you want to run over in runLocal.C. The path of the dataset is saved in a text file in the sourceFiles directory. For now, it points to a file called TEST.txt which gives the path to data from the double muon stream stored in my castor space. This file can be copied from my local space on lxplus, ~naodell/public/TEST.txt. Make sure it is place in the sourceFiles directory. Finally, to run, just enter

root -l runLocal.C

19 July

Draft of ntuple production instructions

From lpc

   cmsrel CMSSW_4_2_4
   cd CMSSW_4_2_4/src

   cvs co -d Higgs/ntupleProducer UserCode/NWU/ntupleProducer
   cvs co -r V110523_BE -d DataFormats/AnomalousEcalDataFormats UserCode/csander/DataFormats/AnomalousEcalDataFormats
   cvs co -r V110523_BE -d PhysicsTools/EcalAnomalousEventFilter UserCode/csander/PhysicsTools/EcalAnomalousEventFilter
   cvs co -r V110523_BE -d Sandbox UserCode/csander/Sandbox
   cvs co -r V19MAY2011_v3 JetMETAnalysis/ecalDeadCellTools
   cvs co -r V04-04-04 JetMETCorrections/Type1MET

   scram b

   cd Higgs/ntupleProducer/test
   cmsenv
   crabenv
   
   

Then modify crabNtuples_*.cfg (Data or MC depending on the sample) to your liking. By default it will save to dcache (which I think is preferable), but save it wherever you want. The config file (ntupleProducer_cfg.py) also needs to be modified depending on whether you are running over MC or data. It should be sufficient to change the tag name.

15 July

Questions for HWW lepton fake team

  • The fake rates are calculated using a combination of data from the Jet and Photon streams. What is the relative contribution from each stream? I see that the vast majority of events come from the Photon stream.

  • How do you calculate the fake rate from the W+Jets sample? Do you do it by MC truth matching or do you use the same method as for the data, but exclude the electron coming from the W? If it is the former, do you consider a fake to be a reconstructed electron not coming from a W or Z?

7 July

Acceptances from HZZ analysis

The order of Andrey's cuts are:

  • some preselection

  • 76 < Mll < 106

  • cut on pT,l1 and pT,l2

  • veto on b-jets

  • cut on MET

The details of the preselection need to be clarified, but it should apply some standard event quality cuts plus a requirement that the two leptons pass ID, isolation, pT, and eta cuts.

Dataset number of events passing μ (e) % of events μ (e)
ttbar 17541 (11974) 2 (1)
DY->ll 461527 (302982) 75 (53)
ggH200 22501 (15610) 40 (28)

19 June

Lepton fakes

To do

  • Determine lead jet pT cut: this needs to be applied to improve agreement between the pT spectrum in the calibration sample and W+jets MC. It also allows for the determination of the associated systematic uncertainty.

  • Get the fake rates from the calibration samples (Photon and Jet streams combined) and MC samples (W+jets and QCD).

13 June

BRM

Remote connection to access DIP cernts.cern.ch --> cerntscms01.

Fake studies

This is going to be somewhat useful for both the VBF analysis and the H->ZZ->llnunu analyses. First step in setting up the analysis is to measure the fake rates parameterized by pT and η (and maybe φ). This I could do now with the PromptReco jet triggered sample. At first I'll have to do this without a trigger selection since I'm not quite sure which one to choose. Prescales probably won't matter, but I will want to maximize the number of objects.

10 June

BCM1F

Some rough calibration factors, shifting everything to 60 ns

Channel pilot peak  
1/1 60.5 -0.05
1/2 62.1 -2.1
1/3 51.35 +8.65
1/4 54.25 +5.75
2/1 60.3 -0.3
2/2 62.0 -2.0
2/3 54.9 +5.1
2/4 54.1 +5.9

7 June

BCM1F

Some results from the TDC tests. We ran a pulse generator directly into the fan in/fan out in with varying parameters to try to isolate the source of the noise observered in channels 1/2, 1/3, 2/2, and 2/3. I put the relevant links here on the calibration twiki, https://twiki.cern.ch/twiki/bin/view/Main/BCM1FCalibrationStudies.

VBF lepton fakes

I hope this is still useful. I've made a ntuple producer just for the fakes. I don't save jets in it for now which might be a problem if I need to include a VBF selection or want to parameterize the fake rate in terms of so jet quantities. For now I save the 4-momentum of the tracks for the electron and muon denominators. The basic procedure is outlined in AN-2009/120. For the case of the electrons, I match gsfTracks to super clusters and then apply some loose cuts. What I need to determine is if I will need some extra information in the case that I later want to apply cuts for electron ID and isolation. This goes for the loosest denominator for the muons which is just an isolated track (the other two are muon objects so I just save them as I would any muon). By the end of today I should produce some ntuples and get some basic numbers.

1 June

BCM1F

To access the 30 minute BMC1F bunch data log onto brmbcmctrl3, and they are located in /d1/brmpro/data/Transfer/cache. Yay!

30 May

BCM1F

I'm going to need to determine the correct offset factors for each channel. The best way to determine this (I think) will be to look at non-colliding bunches, fit a gaussian, and compare the mean of the fit to the center of the bunch. Just from eyeing the centroid of the first colliding bunch for Fill 1815 I would say the offsets are approximately,

Channel Offset
11 -8
12 +5
13 -4
14 -1
21 -7
22 -3
23 -2
24 +6

27 May

BCM1F

Apparently the BCM1F has not been set correctly for the past year or so, so we are now analyzing the effect of changing the high voltage. I made a twiki, HERE!. From Anne:

1) non-colliding bunch associated to beam 1
2) non-colliding where associated to beam 2
3) the pilot
4) a normal colliding bunch.

We need to have these samples for:
100 V old thresholds
250 V old thresholds (the sample you are using is fine)
250 V thresholds of last night

Right now, in the twiki, there are fits for colliding bunches at 250V with the old thresholds. So I will need to do 12 sets of fits, where in the case of colliding bunches I have to also do a series of double Gaussian fits.

25 May

Back to the grind after three weeks of nonsense and travelling. I need to get my thoughts together on what to do now. Some things,

  • I'm meeting with Anne later in the afternoon; she already gave an outline of what should be done in an email, but we'll likely work out what the deal is going to be with shifts.

  • For the VBF stuff I have two things I should be looking at: (1) make overlays of some of the standard distributions I've been looking at, and (2) start getting the objects to look at the fakes.

  • τ triggers: catch up on this shit... I had a couple of items that were requested, but that was so long ago that I'm not quite sure if they still need to be done...

18 May

Candidacy passed!

τ triggers

A lot to do here...

From Michal:

make new evaluation of e+tau trigger rates
as a function of el and tau pT with following conditions:
   1. Electron as it is (CalIdVT_CalIsoT_TrkIdT_TrkIsoT), while tau with
tightIso instead of looseIso
   2. Electron identified as current one
(CalIdVT_CalIsoT_TrkIdT_TrkIsoT), but with |eta|<2.1
(L1_SingleEG12_eta1p17), tau as it is (looseIso)
   3. Electron with looser CaloIso, namely CaloIsoVL
(https://twiki.cern.ch/twiki/bin/viewauth/CMS/EgammaWorkingPointsv3),
rest without changes 

From Artur:

 
   * can you send me updated plot  of the rate vs cut (tau Et, mu Pt w/o
mu isolation) for 1E33?
 (page 22 of attached pdf)

   * can you check the rate table, and fill the missing numbers?
 (page 23 of attached pdf)
   

From Simone:

   Nate, could you provide rate numbers for 2E33 for the following triggers:

   HLT_Ele18(with usual EleID and Isolation)_LooseIsoTau20
   HLT_Ele18(with usual EleID and Isolation)_TightIsoTau20
   

Trigger name L1 seed Prescale Rate(Hz)
Electon
OpenHLT_Ele18_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_LooseIsoPFTau20 L1_SingleEG12 1 6.2 +- 0.04
OpenHLT_Ele18_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_TightIsoPFTau20 L1_SingleEG12 1 4.2 +- 0.04

1 May

Higgs

After speaking with Pietro, we decided I would look at fake rates for central objects in the VBF topology. I guess we have determined that the QCD BG will be well modelled based on what I showed before. Now we want to look at central leptons (e and μ). This will likely rely on a data-driven technique that makes use of 'fakeable' objects; in the case of electrons these will be (gsf)tracks associated to superclusters. So I will need to add the latter two things to my ntuple to carry out the study. For muons I may need some extra information from the muon chambers. For reference,

29 April

τ trigger update

For the meeting today:

  • Some plots showing the rate dependence on τ and/or &muon; pT.

  • Show table with changes to rates from Josh's modifications with current numbers.

Trigger name L1 seed Prescale New changes (1e33) 1e33 2e33 Run 163334 (3.5e32)
Electon
OpenHLT_Ele15_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_LooseIsoPFTau15 L1_SingleEG12 1 11.9 ± 0.1 10.6 ± 0.1 21.2 ± 0.2 6.02
OpenHLT_Ele15_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_LooseIsoPFTau20 L1_SingleEG12 1 7.71 ± 0.05 6.6 ± 0.1 13.3 ± 0.1 4.1
OpenHLT_Ele18_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_LooseIsoPFTau15 L1_SingleEG15 1 9.23 ± 0.05 --- 16.8 ± 0.1 ---
OpenHLT_Ele18_CaloIdVT_TrkIdT_CaloIsoT_TrkIsoT_LooseIsoPFTau20 L1_SingleEG15 1 6.2 ± 0.04 --- 11.04 ± 0.12 ---
Mu
OpenHLT_IsoMu12_LooseIsoPFTau10 L1_SingleMu10 1 10.85 ± 0.07 7.243 ± 0.1 --- 4.3
OpenHLT_IsoMu12_LooseIsoPFTau15 L1_SingleMu10 1 5.98 ± 0.06 --- 11.3 ± 0.1 ---
OpenHLT_IsoMu15_LooseIsoPFTau15 L1_SingleMu10 1 3.13 ± 0.04 --- 6.1 ± 0.1 ---
OpenHLT_Mu15_LooseIsoPFTau20 L1_SingleMu10 1 13.68 ± 0.08 15.2 ± 0.1 30.4 ± 0.2 5.9

19 April

DPS

I think I've taken care of the b-tag issue. There were at least two bugs in the filling of the b-tag discriminator. First, Steve used a box cut to associate jets using b-tags to the PFJets. For the Δφ cut he didn't take the absolute value so that the comparison could conceivably have allowed an association between jets thought their Δφ separation was large. This was unlikely to happen, though, because the Δη cut was so tight (< 0.005).

The bigger problem is that the b-tags were made using some CaloJet collection (or possibly some jet collection using a different clustering algorithm) so that jet axis was sometimes rather far apart and jets were not being associated to the b-tag seed value. This meant certain jets were not given any value for the b-tag whereas they should all be given some value. This could have been addressed by loosening the dR of the jet association so that the b-tag jet lied within the cone of the PFJet. This still could have problems; instead I changed the collection that is used to generate the b-tags to the collection I save in my ntuples (ak5PF jets) and check that dR == 0. This seems to work; just need to check the results.

8 April

DPS

I changed the ordering of the if statement to fill the jet containers (b-tag/no-tag) and included a requirement that no-tag jets have a discriminator value less than 0.

ptHat NInit (106) Nfinal,b-tag Nfinal,noTag(103) Normalization (nb-1)
[15,30] -- -- 0.01 6.68
[30,50] 3.26 36 15.3 61.37
[50,80] 3.19 1142 158 501.8
[80, 120] 3.2 6415 428 4090.3
[120,170] 2.74 13929 572 26499

τ trigger

Looks like the dz variables are a go. Some distributions without any cuts on the twiki .

7 April

DPS

May have fixed it

6 April

DPS

The issue with the b-tag ptHat distribution needs to be addressed. Before and after selection of events with b-tagging on and off,

ptHat_preBranch.png ptHat_2Bjets.png

I'm working under the assumption that I have a bug and it probably comes from when I fill the b-jet container. Possibly I am overcounting bjets...

ptHat NInit (106) Nfinal,b-tag Nfinal,noTag(103) Normalization (nb-1)
[15,30] -- -- 0.01 6.68
[30,50] 3.26 19 15.3 61.37
[50,80] 3.19 1545 158 501.8
[80, 120] 3.2 13229 428 4090.3
[120,170] 2.74 36118 572 26499

3 April

Weekend trigger shifts cancelled! I need to get my to-do in order here now that this nonsense is out of the way.

BRM

Meeting with Anne & co. on Monday to discuss what I have so far and how to proceed. Previous entry has where things currently stand...

Training to become a BRM expert will begin this week as well.

τ trigger

I'll check the vertex z distributions for taus, leptons, and vertices tonight. Then right a simple analyzer that makes Δz distributions. Should have something tomorrow or Tuesday so I can begin calculating the effect on the rates and PU enhancement.

DPS

A couple things here,

  • More MC! I'll run for 2 million events with ptHat 15/10.

30 March

Trigger training shift...

BRM

See the twiki. At the bottom is the to-do...

The first thing I should figure out is how to get a file that has the bunch-by-bunch intensities. It basically just needs two columns, one for each beam. I think the distinction between colliding, non-colliding, ghost/parasitic, and empty could be done simply by checking the intensity and whether both beams are filled. So for one, I need to find a source of such a file; kind of like this, but in a format that is easily read by a script, like a txt file. Then once I have this, I identify the bunch types and store this in an array. Then as I loop over the data I look up the status of the bunch each time I have hit. Once this is done, I fill the relevant histogram. Simple enough.

In the case of non-colliding bunches, I'll take the surrounding (+-1) bunches and include them in the rate (for now). Also, I will need to subtract off the pedestal. To do this, I subtract off the average of the sum of bunches [-2,-4] and [2,4].

28 March

BRM

Since I have to give an update on what I've done for the BCM1F discriminator (uhh...well...), I need to make something today. Currently I can get the rates, the bunch structure plots, and some miscellaneous other stuff. This should, more or less, be sufficient for what Nicola asked for. Because I have only one day to prepare something for the meeting now, I think the simplest thing to do is look at the rate for colliding and non-colliding bunches, compare the rates to each other and compare them to the bunch scheme available on cmswbm.

27 March

DPS

So...I developed a jet-parton association scheme starting about a week or so ago. The first attempt was not quite successful. In that case I did the following:

  • Get genParticles associated to jet.

  • From HepMC::GenEvent get particles (partons) that are the products of the hardest scattering (status = 23).

  • Trace ancestors/descendents of both collections to status = 71 particles.

  • Look for overlap.

The problem with this approach is that there is a huge amount of overlap and doing something simple, like the fraction of matched from one parton to the total matched, doesn't really help. So after trying the complicated thing, I am now going to use a simple dR association (with a possible dpT selection as well). Just from looking at the printouts, this looks very promisiing. So the first step is make some histograms,

  • dR between jet and nearest parton

  • dR between jet and second nearest parton

  • dpT between jet and nearest parton

  • number of matched jets (identify parton), not matched, multiple matches. Do this with a tight cut (dR < 0.5), something a bit looser (dR < 0.8 of 0.9) and just taking the nearest.

Some of the plots are here.

jetPartonDR.png

18 March

τ trigger

There was a bug in the oHLT ntuple code (something with track seeding from the pixels) so I am recalculating the rates for the usual trigger paths as well as the PU enhancement factors. Michal Bluj also requested that I look at the effect of lowering the threshold for the CaloId on the electron x-triggers.

DPS

Because we have been plagued by useless fit results, I decided we would need to clean up our signal sample so that it is in fact a signal sample. This requires that the jets which are paired to produce our DV distributions are 'tagged' so that we can say that they actually come from seperate hard scatterings. Digging around in the MC information has revealed that this is no trivial task. Roughly, what I do, or plan to do, is 1) associate recoJets to genJets, 2) retrieve collection of genParticles used to construct genJet, 3) ?.

10 March

So much for being 'daily'. Anyway, I'll summarize what I've done over the last couple of days.

DPS

A couple days ago I made a number of fits using our isolated dijet scheme, here. Smin and ΔS seem consistent, but the other two are certainly not. I would explect, though, that Smin and SpT would agree because they are basically the same, and Sphi and ΔS would be closer to agreeing because they both measure angular correlations. Realizing that this fitting method is sensitive to low statistics, I restricted the ranges to only include non-zero bin ranges. The result is here. The only DV that changes is Smin; it agrees quite well with SpT now as I anticipated. The reason for this is likely due to the difficulty when the bin 9-9.5 and (maybe) 6.5-7 are included.

The issue is that the fitting method is obviously sensitive to the fitting range in a somewhat unpredictable way as can be seen here. What I have done here is restrict the bin range to [1-9], [1-8], [1-7], and [1-6], respectively. The results are definitely sensitive to the selected fit range and are quite unstable. The binning used is also going to be somewhat of a concern. I am repeating this in a more systematic way.

Assuming this effect arises due to sensitivity to statistics, I consider our modified standard selection (lower pT on jets 3 and 4) which has a much higher event yield. Here are the fit results using the entire range for each,

DV pT,3 + pT,4 DPS forced Pythia8 QCD Chi2/nDOF Probablility
SMin 20 0.683315 +/- 0.0453319 0.316677 +/- 0.0452781 1.14289 0.301551
  30 0.236942 +/- 0.0684679 0.763077 +/- 0.0695081 1.23038 0.225534
  40 0.257489 +/- 0.0928911 0.742523 +/- 0.0942959 0.558875 0.929939
  50 0.0677737 +/- 0.190128 0.9324 +/- 0.192548 0.598892 0.903456
deltaS 20 -0.0511974 +/- 0.0419507 1.05119 +/- 0.0447697 0.44778 0.969959
  30 -0.107051 +/- 0.0454348 1.10705 +/- 0.049374 0.254022 0.998788
  40 -0.148065 +/- 0.0636131 1.14806 +/- 0.0692001 0.241436 0.999119
  50 -0.244766 +/- 0.102777 1.24488 +/- 0.111892 0.219256 0.999525
Spt 20 0.0244053 +/- 0.0460903 0.975627 +/- 0.0465043 1.93703 0.000467081
  30 0.130362 +/- 0.0477762 0.869618 +/- 0.0487975 0.932268 0.589078
  40 0.179387 +/- 0.0945634 0.820804 +/- 0.0948816 0.570859 0.984436
  50 0.102581 +/- 0.104497 0.897379 +/- 0.106065 0.386395 0.999766
Sphi 20 0.274443 +/- 0.0376487 0.725558 +/- 0.0378575 0.949458 0.551997
  30 0.186563 +/- 0.0480685 0.813483 +/- 0.0483528 0.753585 0.848784
  40 0.181946 +/- 0.134658 0.817969 +/- 0.135743 0.447943 0.997716
  50 -0.169413 +/- 0.162867 1.16949 +/- 0.162076 0.378076 0.99962

3 March

DPS

Now I add the topology selection and get the following yield,

Initial number of events: 801803
Pass HLT selection 754542
One and only one primary vertex: 552959
Vertex displacement in z < cm: 552746
Vertex ndof >= 4: 552415
abs(η) < 2.4 : 323067
Miscellaneous jet quality cuts: 313393
p_{T,1/2} > 20 GeV and p_{T,3/4} > 10: 43176
p_{T,1} + p_{T,2} > 75 GeV: 15712
isolated dijets 313

If I make no restriction on the run range, this becomes

Initial number of events: 7814238
Pass HLT selection 2197572
One and only one primary vertex: 1122399
Vertex displacement in z < cm: 1121557
Vertex ndof >= 4: 1120807
abs(η) < 2.4 : 657983
Miscellaneous jet quality cuts: 638001
p_{T,1/2} > 20 GeV and p_{T,3/4} > 10: 87281
p_{T,1} + p_{T,2} > 75 GeV: 31744
isolated dijets 627

2 March

DPS

Testing out JetMETTau for Run2010A (runNumber < 137028) with residual corrections and extending range. Without residual corrections the event yield is,

Initial number of events: 801803
Pass HLT selection 754542
One and only one primary vertex: 552959
Vertex displacement in z < cm: 552746
Vertex ndof >= 4: 552415
abs(η) < 2.4 : 323067
Miscellaneous jet quality cuts: 313393
p_{T,1/2} > 20 GeV and p_{T,3/4} > 20: 3906
p_{T,1} + p_{T,2} > 75 GeV: 2644

If I lower the pT threshold on the jets in the second pair I get

p_{T,1/2} > 20 GeV and p_{T,3/4} > 10: 43176
p_{T,1} + p_{T,2} > 75 GeV: 15712

I discovered that the residual corrections are being stored incorrectly; I was applying the L2 and L3 corrections a second time while retrieving the residual corrections. This will have to be corrected in the next round of ntuples.

1 March

HCAL Conditions

From Anton:


Tag: HcalRespCorrs_v3.10_offline,
use the etaPhiCorr_030111 files.
IOV:
"period1" -> 1-148655
"period2" -> 148656...eot

"Both phi and eta-dependent scale corrections are included. The phi
corrections are identical to v2.10 and scale corrections are derived on top
them. The effect of cluster contamination is not accounted for yet.
The IOV separation is to account for hardware changes in HF (see comments
for v2.01)."


Similarly, make
Tag: HcalRespCorrs_v2.10_offline,
use the phiCorr_030111 files.
IOV:
"period1" -> 1-148655
"period2" -> 148656...eot

"Only phi corrections are included. The phi corrections in HBM
(ieta=[-1,-15], iphi=6), (ieta=[-1,-7], iphi=32) are set to 1.
Channels HB(-1,36), HB(-5,36), HB(-9,36) are adjusted down relative to v2.01
to correct for behavior seen in Dec re-reco.
The IOV separation is to account for hardware changes in HF (see comments
for v2.01)."

   

Done and done.

DPS

We tentatively had decided to switch to running over MinBias stream data from the Nov 4th ReReco. Initially it looked as though we might have had a sufficiently large number of events to proceed, but once you include the HLT selection HLT_L1_BscMinBiasOR_BptxPlusORMinus the event yield is negligible,

Initial number of events: 14552305
Pass HLT selection 1489011
One and only one primary vertex: 627543
Vertex displacement in z < cm: 627078
Vertex ndof >= 4: 626643
abs(η) < 2.4 : 132112
Miscellaneous jet quality cuts: 127465
p_{T,1/2} > 20 GeV and p_{T,3/4} > 10: 1514
p_{T,1} + p_{T,2} > 75 GeV: 344

On the other hand if you run without the trigger selection, you get a high event yield, but funky distributions, [[][here]]. Radek is looking into what triggers have all of these events, but it probably will just make this more complicated. So, back to comissioning and Jet(MET(Tau)) data. We want to increase our statistics so I am going to look at using an alternative trigger which does not become prescaled as quickly as the Jet_15U trigger.

23 Feb

HCAL conditions

Anton sent a request for me to load test tags for three new types of conditions objects: RecoParams, LongRecoParams, and MCParams. They need to be run in CMSSW_4_2_0_pre5 (instructions for setting up the environment are here. There are a couple extra packages that need to be checked out on top of the release,

   cvs co -r V03-07-01 CondFormats/HcalObjects
   cvs co -r V06-14-13 CondFormats/DataRecord
   cvs co -r V03-04-00 CondCore/HcalPlugins
   cvs co -r V02-19-00 CalibCalorimetry/HcalAlgos
   cvs co -r V02-13-00 CalibCalorimetry/HcalPlugins
   cvs co -r V02-05-00 CondTools/Hcal

and the test tags are Salavat's lxplus space, /afs/cern.ch/user/a/abdullin/public/Conditions420. I was able to load all of the tags to a local database, but, alas, when I attempted to load them to Orcon by way of dropbox I was unsuccessful. Andrey took care of the rest.

τ trigger

Pile-up effects in the cross triggers seem to be under control compared to previous number (see here), but there is a problem: I don't know why these numbers are different than the previous numbers. As far as I can tell nothing changed. I'm starting from "scratch" to recreate the rates now.

DPS

New ntuples are in the oven. I messed up the trigger mapping and the typecasting of certain variables in the previous iteration so I'm redoing them now. I also have to run over Run2010A MinBias data since Mayda wants to switch streams (or change horses mid-stream). I suspect the motivation is, given we are already up against 2011 data-taking, that using comissioning data will not appear appealing to the ARC; they'l consider it a waste of time. Alas, the time already spent.

18 Feb

τ trigger

I reran the rates for some of the electron x-triggers and renamed them so they conform to the eGamma POG naming conventions (which changed again). I also ran over some of the pile-up samples (finally), and saw that there is a significant enhancement in the rates when pile-up is present especially for the e-τ triggers. There are some considerations to make when reviewing these results: the statistics are low so making a definitive statement is out of the question, and there are certain results which are counter-intuitive to the extent that they are either wrong or corrupted by the low statistics. All of these results are up on the twiki and were presented in the τ trigger meeting this afternoon.

14 Feb

DPS

First order of business: Look at events where we can isolate dijet pairs. This is done by making two additional requirements : pairing the highest two pT jets, then the next two, require that the first pair be located on one side of the detector, in z/η, and the other pair in the other side. On top of this require that Δφ of the pairs be greater than 2. The distributions are here

11 Feb

DPS

It's not, at this point, exactly clear where this analysis is going. There are systematic studies that need to be done, first of which is the effect of the JES. One useful thing would be to look at the correction factors in η and pT. I need to get the overall uncertainty for the current corrections and then see how this affects the fits like before.

The meeting with Rick Field was useful. He had three useful suggestions,

  • Look at jets with a tighter cone size, say 0.2 rather than 0.5. This seems like it would be non-trivial. Since these collections aren't available in the official samples it would require that I rerun the jet collection producer which could be time time consuming. I'll hold off on this until some compelling reason to pursue this presents itself. Currently available collections are,

      vector               "ak5PFJets"
      vector               "ak7PFJets"        
      vector               "iterativeCone5PFJets"     
      vector               "kt4PFJets"       
      vector               "kt6PFJets" 
      

  • Consider the effects of a scan of the pT cut on the second jet pair. That is, require the first pair to be hard (pT > 30 GeV), and allow the second pair to go as low at 10 GeV and raise it to the threshold of the first two. This is straight forward enough and will allow us to see how strongly dependent DPS is on our energy scale.

  • Finally, or firstly, look for the signal. Since we assume DPS will consist of two uncorrelated, superimposed dijet events, look for exactly this. In order to avoid jets "eating each other up" I should include a Δη restriction. This value should be scanned.

I also ran a fit over what I have right now, even with the funky signal sample. The results are

      Fit param:		deltaS:			S_min			sPt			sPhi

      SIG fraction: -0.0773653 +/- 0.0795894	-0.0189561 +/- 0.197025	0.158164 +/- 0.0880736	-0.0881456 +/- 0.177555
      BG fraction:  1.07735 +/- 0.0855524	1.01902 +/- 0.198393	0.841907 +/- 0.0882447	1.08816 +/- 0.178004

      Chi^2/NDF:		0.171545			0.87513			1.01845			0.594215
      Probablilty:		1			0.609827		0.433988		0.970685

      

The priority for now is on regenerating the signal sample. It should be a large quantity so I should get the filter that Radek uses, given that mine doesn't work. Next I should pursue the last two items from Rick.

7 Feb

τ trigger

Simone had asked that I look at the effects of pile-up by running with our triggers over some MC samples. The result is that none of the events pass the triggers which is not surprising because the one sample is MinBias and one is QCD with ptHat > 15 GeV.

4 Feb

HCAL conditions

Updates need to be made for reprocessing of HI data. The only tags for which this is the case are Gains and RespCorrs.

  • Gains : The current tag is usable since the current IOV in the tag, HcalGains_v2.32_offline, begins with a run the precedes the beginning of the HI run. The need changes are: to include the ZDC calibration constants from Taylan for IOV 148655...151249; then propagate the changes from the HLT tag to the offline tag.

  • Respcorrs : Append the current IOV in the HLT tag.

3 Feb

τ trigger

Was seeing a discrepancy between my numbers and Aruna's. The source was the track isolation which I forgot to change the isolation variables from ints to floats (again). The result was that all of the rates were underestimated by some factor. It is now all under control, but the ntuples do not have the proper photon isolation. Comparison for validation.

DPS

Some new distributions with out the single primary vertex requirement. For these I require that there be at least one primary vertex and that all of the tag jets are associated to this vertex. There aren't really any significant changes, but two things are worth noting. The number of events increases by 50% and the dip in S is less pronounced.

1 Feb

τ trigger

Doh! Submitted some bad code. Set the photon isolation to the track isolation. The effect is probably not very significant, but it could do something depending on whether the other cuts have a tendency to enhance the effect of high energy photons in the isolation annulus.

27 Jan

DPS

Had to rerun the data ntuples to include a low threshold (minBias) trigger for the trigger study. Also included technical triggers this time though their effect is completely trivial in my case (< 10 events). I need to implement a script that will allow me to run over multiple datasets, binned QCD samples for instance, in one batch submission. It would also be useful to split the submissions by individual root files.

τ trigger

25 Jan

τ trigger

Further investigations of τ-lepton x-trigger for 1e33 menu. Want to verify Aruna's numbers and add some extra options. See twiki for details.

20 Jan

Focus has been on DPS and τ trigger development.

DPS

Recreated trigger study to determine 2-jet HT cut that makes us fully (or ~99%) efficient.

17 Jan

τ trigger

Current results are on twiki and in a talk given on Friday, here. Simone suggested that Pedrame, Aruna, Anton, and myself meet today to discuss where we should be going from here. Either way there are two things that need to be done ASAP:

  • Update ntuple code with changes to track/photon isolation in the repository.

  • Run μ-τ x-trigger analysis for e-τ x-trigger.

VBF

There a couple of things that will need to be done:

  • Clean up leading jet η distribution. Contact QCD/JMT people for fix.

  • Look at μ and e distributions. Will need to get appropriate MC and run over e/γ sample for electrons.

  • Unfold distributions using bin-by-bin corrections.

MPI

Update all 4-jet DPS plots to include 38X MC and rereco. Going to need to do the same for the b-jets, but first need to fix bug with luminosity normalisation of distributions.

BRM

Just run the damn code! Need to have something in this department by the end of the week. Can just be something preliminary, but something.

12 Jan

τ trigger

Some updates on the twiki. To summarize, We have looked at the effect on the rates of varying the pT of the leading track and the track(s) in the isolation annulus. Cutting on the lead track pT, but this is probably just because the variation in the value of the threshold is much greater compared to that of the isolation track. We also looked at the effect of isolating the muon and varying the η threshold of the muon. In the former case, the effect is quite significant, in the latter, not so much.

Simone's assessment is that we use the isolated muon. When ntuples are ready I am to do the same for electrons.

10 Jan

VBF

Met with Pietro this morning to discuss some plans. Need to give a talk on Friday to the forward physics people. The contents will be much of the material that I had already presented to the VBF folks with the following additions:

  • Clean up bad runs/ excess of events in -2.5 < |η| < -2.0.

  • Include plots from discussion with Pietro today: Δφ vs. Δη. This should be done in two ways: one is the normal way; the other is to do a scan of the minimum Δη value allowed and take the sum of all events for a given Δφ. The latter needs some thought and may be slightly involved.

  • Investigate in particular the case that you have two jets in the central region.

  • Include leptons in the central region. This will require ntuplizing the mu and e/gamma stream data and the relevant MC.

  • Unfold spectra to go from detector level distributions to hadron level.

τ trigger

Meeting with Simone at 2:30. Need to work out bug. Edit: done.

Conditions

Loaded a test tag for respCorrs in response to a request by Anton. Dropbox seems to be down.

09 Jan

Back to the grind...

Tau trigger

Today I am doing a scan of the trigger rates for the tau-muon cross-section. The four relevant parameters in the trigger are the ET of the seed calotower of the tau; the leading track to reconstruct the tau; the maximum allowed track pT in the isolation annulus; and the maximum allowed photon pT in the isolation annulus. The latter two parameters, i.e., the isolation parameters, were originally just saved as the number of tracks/photons in the isolation region. Since the goal is to optimise the tau portion of the trigger, which will be somewhat involved, I've made a twiki to post my results to here.

There is a bug, perhaps, that is preventing the isolation cut from having any effect. Need to work this out ASAP.

-- NateOdell - 09-Jan-2011

Topic attachments
I Attachment History Action Size Date Who Comment
PNGpng photonPt_Combined.png r1 manage 8.5 K 2011-12-05 - 16:22 NateOdell  
Edit | Attach | Watch | Print version | History: r130 < r129 < r128 < r127 < r126 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r130 - 2012-03-25 - NateOdell
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Main All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2023 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
or Ideas, requests, problems regarding TWiki? use Discourse or Send feedback