Search for doubly charged Higgs bosons in e-μ channel

at CMS

1

Feb. 23, 2013 Jongseok Lee(Sungkyunkwan University)

Contents• 7TeV Samples

• Lepton cut variables

• Data driven background estimation

• Comparison with Mario

• Systematic uncertainties

• Effect of tau

• 8TeV analysis

• Summary and Plan 2

Samples• Used CMSSW version : 4_4_2

• MC signal samples (√s=7TeV) : full simulation data → H++H-- to l+l+l-l-

→ M(H±±) = 150, 170, 200, 250, 300, 400, 500GeV • MC background samples : full simulation data → ZZJetsTo4L, WZJetsTo3LNu TTTo2L2Nu2B, DYJetsToLL

• Data samples : 2011 data, Lint(μμ)= 4.75 fb-1 , Lint(ee)= 4.98 fb-1 ,

Lint (eμ)= 4.98 fb-1

3

calculated by Michael Spira

Mμμ

4

nm>=2, pt>10, |η|<2.4, global, good, mu_hits_val>20, mu_nstation>0, relIso<0.2, 60<Mμμ<120

lumi : 4.983 fb-1

cross section of DY: 3048 fb

lumi : 4.753 fb-1

cross section of DY : 3048 fb

lumi : 4.753 fb-1

cross section of DYfor 60<Mll<120 : 2916 fb

Standard model cross sections for CMS at 7TeV : https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSections

Muon cut variables

5

nm>=2, pt>10, |η|<2.4, global, good, nhits_val>20, nsation>=1, relIso<0.2, 60<Mμμ<120

7TeV, 4.753 fb-1

7TeV, 4.753 fb-1

7TeV, 4.753 fb-1

7TeV, 4.753 fb-1

7TeV, 4.753 fb-1

7TeV, 4.753 fb-1

Mee

6

ne>=2, pt>10, |η|<2.5, charge consistent, eidLoose, relIso<0.2, 60<Mμμ<120

lumi : 4.983 fb-1

cross section of DY: 3048 fb

lumi : 4.983 fb-1

cross section of DYfor 60<Mll<120 : 2916 fb

Standard model cross sections for CMS at 7TeV : https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSections

Electron cut variables

7

ne>=2, pt>10, |η|<2.5, eidLoose, relIso<0.2

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1 7TeV, 4.983

fb-1

7TeV, 4.983 fb-1 7TeV, 4.983

fb-1

7TeV, 4.983 fb-1

Sideband method

8

yuncertaint lstatistica N assumed isit then y,uncertaint lstatistica N If

sidebands. in the events data of # total theis N

data, from estimatedregion signal in the events background of # theis N where

,1 N

1 oferror relative a and

,1 N N

yuncertaint lstatistica N assumed isit then y,uncertaint lstatistica N if

, N assumed isit then 0, N if

MC. by the prediced as regions sideband and signal in the events theare N and N where

,N

N

region. signal in theon contributi background theestimate toused is method sidebandA

BGSRBGSR

DataSB

BGSR

DataSB

DataSBBGSR

SRSR

SRSB

SBSR

SB

SR

# of events at 7TeV & 4.753fb-1

9

mmmmMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 0.879±0.019

1.50±0.18±0.27

0

170 0.855±0.017

1.46±0.18±0.26

0

200 0.486±0.013

0.815±0.098±0.144

0

250 0.296±0.010

0.491±0.059±0.087

0

300 0.118±0.006

0.194±0.024±0.034

0

350 0.109±0.006

0.179±0.023±0.032

0

400 0.072±0.005

0.118±0.016±0.021

0

500 0.025±0.003

0.040±0.007±0.007

0

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10, in signal region (5σ)

± statistical error ± α error ± relative error

# of events at 7TeV & 4.983fb-1

10

eeeeMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 0.451±0.013

0.889±0.052 ±0.216

0

170 0.381±0.012

0.745±0.044 ±0.181

0

200 0.309±0.011

0.600±0.036 ±0.145

0

250 0.072±0.006

0.136±0.013 ±0.033

0

300 0.049±0.004

0.092±0.009 ±0.022

0

350 0.017±0.002

0.031±0.005 ±0.008

0

400 0.017±0.002

0.031±0.005 ± 0.008

0

500 0.0014±0.0007

0.0025±0.0013 ±0.0006

0

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10, in signal region (5σ)

± statistical error ± α error ± relative error

# of events at 7TeV & 4.983fb-1

11

ememMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 1.40±0.02

2.39±0.10 ±0.36

4.0±2.0

170 1.20±0.02

2.09±0.09 ±0.31

3.0±1.7

200 0.642±0.015

1.14±0.05 ±0.168

1.0±1.0

250 0.334±0.011

0.599±0.030 ±0.087

0

300 0.203±0.009

0.363±0.021 ±0.053

0

350 0.114±0.006

0.203±0.013 ±0.030

0

400 0.068±0.005

0.120±0.010 ±0.017

0

500 0.012±0.002

0.022±0.004 ±0.003

0

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10, in signal region (5σ)

± statistical error ± α error ± relative error

ABCD method

12

signal

backgound

A

C

B

D

A B C D : regionsnA nB nC nD : event counts

relIso<0.2 relIso>0.2

M4l>200

M4l<200

Background is independent for relIso & m4 :

The estimated number of background events in A :

This can be done completely without MC information-> Data driven background estimation

D

B

C

A

n

n

n

n

D

CBestA n

nnn

# of events at 7TeV & 4.983fb-1

13

mmmmregion B region C region D region A A-estimated

data 19.0±4.4

17.0±4.1

142.0±11.9

9.0±3.0

2.3±0.8

MC 13.9±1.3

2.56±0.04

49.0±3.6

11.8±0.1

singnal(150GeV)

33.0±1.2

0 0.048±0.048

241±2

nlep+>=2, nlep->=2, m(H±±)>10

± statistical error

# of events at 7TeV & 4.983fb-1

14

eeeeregion B region C region D region A A-estimated

data 5.0±2.2

7.0±2.6

5.0±2.2

5.0±2.2

7.0±5.2

MC 3.8±0.9

1.6±0.4

8.8±3.3

5.81±0.05

singnal(150GeV)

21.3±1.0

0 0 133±2

nlep+>=2, nlep->=2, m(H±±)>10

± statistical error

# of events at 7TeV & 4.983fb-1

15

ememregion B region C region D region A A-estimated

data 16.0±4.0

14.0±3.7

86.0±9.3

21.0±4.6

2.6±1.0

MC 19.3±1.5

3.8±0.4

57.3±4.7

16.7±0.1

singnal(150GeV)

26.6±1.1

0 0 180±2

nlep+>=2, nlep->=2, m(H±±)>10

± statistical error

Comparison with Mario for 4μ channel

at 7TeV & 4.753fb-1

16

4μ & event selection(case 1)

17

• Muon selectionpT>10 GeV, |η|<2.4, global muon, good muon, hits_val>20, nstation>=1

• Preselection : nlep+>=2, nlep->=2, Mll>10 GeV

• Event selection - mass dependent cut- quad-lepton mass : M4l > cut value

- relative isolation : (IsoPttrack + IsoPtcalorimeter)/pt < cut value

- minimum one event for backgrounds

4μ & event selection(case 2)

18

• Muon selectionpT>10 GeV, |η|<2.4, global muon, good muon, hits_val>20, nstation>=1

• Preselection : nlep+>=2, nlep->=2,

pT,1st>20 GeV, pT,2nd>10 GeV,

Mll>12 GeV,

relIsoworst + relIsonext to worst < 0.35

• Event selection - mass dependent cut- sumpt : ∑pT > 0.6mΦ + 130GeV

- mass window : [0.9mΦ ; 1.1mΦ]

- minimum one event for backgrounds

19

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 247±2

4.45±0.16

39.2 2.0±1.4

170GeV 151±1

3.19±0.16

29.9 2.0±1.4

200GeV 76.6±0.7

1.78±0.15

21.0 2±1.4

250GeV 30.4±0.3

1.02±0.15

12.4 0

300GeV 12.8±0.2

0.688±0.153

7.39 0

350GeV 6.01±0.07

0.557±0.153

4.52 0

400GeV 3.02±0.03

0.488±0.153

2.80 0

500GeV 0.849±0.013

0.415±0.153

1.06 0

preselection, sumpt, mass window (case 2)

# of expected eventspreselection, M4l, relIso (case 1)

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 220±2

0.494±0.166

47.5 0

170GeV 132±2

0.477±0.166

35.0 0

200GeV 69.0±0.7

0.457±0.166

23.7 0

250GeV 28.0±0.3

0.425±0.166

13.5 0

300GeV 11.1±0.2

0.419±0.166

7.38 0

350GeV 5.15±0.08

0.416±0.166

4.31 0

400GeV 3.03±0.04

0.413±0.166

2.92 0

500GeV 0.700±0.014

0.411±0.166

0.90 0

Exclusion limit for 4μ at 7TeV

20

Expected Limit : 385 GeVObserved Limit : 399 GeV

preselection, sumpt, mass window (case 2)preselection, M4l, relIso (case 1)

Expected Limit : 375 GeVObserved Limit : 399 GeV

4μ & event selectionby Mario (case 3)

21

• Muon selectionpT>5 GeV, |η|<2.4, global muon, tight muon id

• Preselection : pT,1st>20 GeV, pT,2nd>10 GeV,

Mll>12 GeV,

relIsoworst + relIsonext to worst < 0.35,

sip<4

• Event selection - mass dependent cut- sumpt : ∑pT > 0.6mΦ + 130GeV

- mass window : [0.9mΦ ; 1.1mΦ]

- minimum one event for backgrounds

22

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 247±2

4.45±0.16

39.2 2.0±1.4

170GeV 151±1

3.19±0.16

29.9 2.0±1.4

200GeV 76.6±0.7

1.78±0.15

21.0 2±1.4

250GeV 30.4±0.3

1.02±0.15

12.4 0

300GeV 12.8±0.2

0.688±0.153

7.39 0

350GeV 6.01±0.07

0.557±0.153

4.52 0

400GeV 3.02±0.03

0.488±0.153

2.80 0

500GeV 0.849±0.013

0.415±0.153

1.06 0

Mario’s result at 7TeV & 4.93fb-1 (case 3)

# of expected eventspreselection, M4l, relIso (case 1)

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 215±2

0.53±0.43

46.5 0

170GeV 131±2

0.51±0.43

34.5 0

200GeV 66.9±0.7

0.48±0.43

23.1 0

250GeV 26.0±0.3

0.46±0.43

12.7 0

300GeV 10.7±0.2

0.46±0.43

7.06 0

350GeV 4.87±0.08

0.46±0.43

4.05 0

400GeV 2.41±0.04

0.46±0.43

2.39 0

500GeV 0.65±0.02

0.46±0.43

0.91 0

Exclusion limit for 4μ at 7TeV

23

Mario’s result at 7TeV & 4.93fb-1 (case 3)

Result for pair production only : 395GeV

Expected Limit : 385 GeVObserved Limit : 399 GeV

preselection, M4l, relIso (case 1)

24

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 220±2

0.494±0.166

47.5 0

170GeV 132±2

0.477±0.166

35.0 0

200GeV 69.0±0.7

0.457±0.166

23.7 0

250GeV 28.0±0.3

0.425±0.166

13.5 0

300GeV 11.1±0.2

0.419±0.166

7.38 0

350GeV 5.15±0.08

0.416±0.166

4.31 0

400GeV 3.03±0.04

0.413±0.166

2.92 0

500GeV 0.700±0.014

0.411±0.166

0.90 0

preselection, sumpt, mass window (case 2)

# of expected events

# of expected eventsM(H++) H++H– to

μμμμBkg.(MC) ScL Data

150GeV 215±2

0.53±0.43

46.5 0

170GeV 131±2

0.51±0.43

34.5 0

200GeV 66.9±0.7

0.48±0.43

23.1 0

250GeV 26.0±0.3

0.46±0.43

12.7 0

300GeV 10.7±0.2

0.46±0.43

7.06 0

350GeV 4.87±0.08

0.46±0.43

4.05 0

400GeV 2.41±0.04

0.46±0.43

2.39 0

500GeV 0.65±0.02

0.46±0.43

0.91 0

Mario’s result at 7TeV & 4.93fb-1 (case 3)

Exclusion limit for 4μ at 7TeV

25

Mario’s result at 7TeV & 4.63fb-1 (case 3)

Result for pair production only : 395GeV

preselection, sumpt, mass window (case 2)

Expected Limit : 375 GeVObserved Limit : 399 GeV

Comparison with Mario for 4e channel

at 7TeV & 4.983fb-1

26

4e & event selection (case 1)

27

• Electron selectionpT>10 GeV, |η|<2.5, |dxy|<0.02, eidLoose

• Preselection : nlep+>=2, nlep->=2, Mll>10 GeV

• Event selection - mass dependent cut- quad-lepton mass : M4l > cut value

- relative isolation : (IsoPttrack + IsoPtcalorimeter)/pt < cut value

- minimum one event for backgrounds

4e & event selection(case 2)

28

• Electron selection pT>10 GeV, |η|<2.5, |dxy|<0.02, eidLoose

• Preselection : nlep+>=2, nlep->=2,

pT,1st>20 GeV, pT,2nd>10 GeV,

Mll>12 GeV,

relIsoworst + relIsonext to worst < 0.35

• Event selection - mass dependent cut- sumpt : ∑pT > 0.6mΦ + 130GeV

- mass window : [0.9mΦ ; 1.1mΦ]

- minimum one event for backgrounds

29

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 150±2

2.67±0.18

30.6 0

170GeV 94.0±1.6

1.95±0.18

23.7 0

200GeV 48.3±0.8

1.32±0.18

16.2 0

250GeV 19.7±0.3

0.808±0.174

9.67 0

300GeV 7.88±0.17

0.607±0.174

5.39 0

350GeV 3.57±0.08

0.516±0.174

3.13 0

400GeV 1.74±0.04

0.481±0.174

1.82 0

500GeV 0.476±0.015

0.448±0.174

0.621 0

preselection, sumpt, mass window at 7TeV & 4.983fb-1

(case 2)

# of expected eventspreselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1)

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 132±2

0.474±0.174

35.0 0

170GeV 83.5±1.5

0.461±0.174

26.6 0

200GeV 42.5±0.8

0.455±0.174

17.5 0

250GeV 17.6±0.3

0.445±0.174

9.91 0

300GeV 7.10±0.17

0.435±0.174

5.36 0

350GeV 3.25±0.08

0.433±0.174

3.05 0

400GeV 1.89±0.05

0.431±0.174

2.01 0

500GeV 0.433±0.015

0.431±0.174

0.579 0

Exclusion limit

30

Expected Limit : 346 GeVObserved Limit : 364 GeV

preselection, sumpt, mass window at 7TeV & 4.983fb-1 (case 2)

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1)

Expected Limit : 341 GeVObserved Limit : 356 GeV

4e & event selectionby Mario (case 3)

31

• Electron selectionpT>15 GeV, |η|<2.5, CiCid, ipT <2mm

• Preselection : pT,1st>20 GeV, pT,2nd>10 GeV,

Mll>12 GeV,

relIsoworst + relIsonext to worst < 0.35,

sip<4

• Event selection - mass dependent cut- sumpt : ∑pT > 0.6mΦ + 130GeV

- mass window : [0.9mΦ ; 1.1mΦ]

- minimum one event for backgrounds

32

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 150±2

2.67±0.18

30.6 2

170GeV 94.0±1.6

1.95±0.18

23.7 0

200GeV 48.3±0.8

1.32±0.18

16.2 0

250GeV 19.7±0.3

0.808±0.174

9.67 0

300GeV 7.88±0.17

0.607±0.174

5.39 0

350GeV 3.57±0.08

0.516±0.174

3.13 0

400GeV 1.74±0.04

0.481±0.174

1.82 0

500GeV 0.476±0.015

0.448±0.174

0.621 0

# of expected events

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 144±2

0.50±0.43

36.7 0

170GeV 92.5±1.47

0.48±0.43

28.2 0

200GeV 47.8±0.7

0.47±0.43

18.7 0

250GeV 19.2±0.3

0.46±0.43

10.5 0

300GeV 8.10±0.17

0.46±0.43

5.82 0

350GeV 3.72±0.08

0.46±0.43

3.32 0

400GeV 1.79±0.04

0.46±0.43

1.89 0

500GeV 0.480±0.015

0.46±0.43

0.62 0

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1) Mario’s result at 7TeV & 4.93fb-1 (case 3)

Exclusion limit for 4e at 7TeV

33

Result for pair production only : 384GeV

Expected Limit : 346 GeVObserved Limit : 364 GeV

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1) Mario’s result at 7TeV & 4.93fb-1 (case 3)

34

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 132±2

0.474±0.174

35.0 0

170GeV 83.5±1.5

0.461±0.174

26.6 0

200GeV 42.5±0.8

±0.4550.174

17.5 0

250GeV 17.6±0.3

0.445±0.174

9.91 0

300GeV 7.10±0.17

0.435±0.174

5.36 0

350GeV 3.25±0.08

0.433±0.174

3.05 0

400GeV 1.89±0.05

0.431±0.174

2.01 0

500GeV 0.433±0.015

0.431±0.174

0.579 0

# of expected events

# of expected eventsM(H++) H++H– to

eeeeBkg.(MC) ScL Data

150GeV 144±2

0.50±0.43

36.7 0

170GeV 92.5±1.47

0.48±0.43

28.2 0

200GeV 47.8±0.7

0.47±0.43

18.7 0

250GeV 19.2±0.3

0.46±0.43

10.5 0

300GeV 8.10±0.17

0.46±0.43

5.82 0

350GeV 3.72±0.08

0.46±0.43

3.32 0

400GeV 1.79±0.04

0.46±0.43

1.89 0

500GeV 0.480±0.015

0.46±0.43

0.62 0

Mario’s result at 7TeV & 4.93fb-1 (case 3)preselection, sumpt, mass window at 7TeV & 4.983fb-1

(case 2)

Exclusion limit for 4e at 7TeV

35

Expected Limit : 341 GeVObserved Limit : 356 GeV

Result for pair production only : 384GeV

Mario’s result at 7TeV & 4.93fb-1 (case 3)

preselection, sumpt, mass window at 7TeV & 4.983fb-1 (case 2)

Comparison with Mario for eμeμ channel

at 7TeV & 4.983fb-1

36

37

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 186±2

5.40±0.18

31.6 4.0±2.0

170GeV 116±2

4.02±0.18

24.1 5±2.2

200GeV 61.7±0.8

2.60±0.18

17.0 0

250GeV 25.4±0.4

1.37±0.18

10.4 0

300GeV 10.3±0.2

0.889±0.174

6.02 0

350GeV 4.77±0.08

0.657±0.174

3.66 0

400GeV 2.85±0.05

0.575±0.174

2.55 0

500GeV 0.703±0.015

0.470±0.174

0.861 0

# of expected events

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 173±2

0.550±0.174

40.6 0

170GeV 106±2

0.529±0.174

30.3 0

200GeV 56.4±0.8

0.486±0.174

20.7 0

250GeV 22.7±0.4

0.454±0.174

11.7 0

300GeV 9.34±0.18

0.442±0.174

6.48 0

350GeV 4.31±0.08

0.437±0.174

3.74 0

400GeV 2.47±0.05

0.434±0.174

2.47 0

500GeV 0.605±0.015

0.431±0.174

0.779 0

preselection, sumpt, mass window at 7TeV & 4.983fb-1

(case 2)preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1)

Exclusion limit

38

preselection, sumpt, mass window at 7TeV & 4.983fb-1 (case 2)

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1)

Expected Limit : 363 GeVObserved Limit : 381 GeV

Expected Limit : 370 GeVObserved Limit : 393 GeV

39

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 186±2

5.40±0.18

31.6 4.0±2.0

170GeV 116±2

4.02±0.18

24.1 5±2.2

200GeV 61.7±0.8

2.60±0.18

17.0 0

250GeV 25.4±0.4

1.37±0.18

10.4 0

300GeV 10.3±0.2

0.889±0.174

6.02 0

350GeV 4.77±0.08

0.657±0.174

3.66 0

400GeV 2.85±0.05

0.575±0.174

2.55 0

500GeV 0.703±0.015

0.470±0.174

0.861 0

# of expected events

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 176±2

0.57±0.43

40.9 0

170GeV 109±2

0.55±0.43

30.7 0

200GeV 56.9±0.8

0.50±0.43

20.8 0

250GeV 21.7±0.3

0.47±0.43

11.3 0

300GeV 9.32±0.18

0.46±0.43

6.42 0

350GeV 4.29±0.08

0.46±0.43

3.69 0

400GeV 2.01±0.04

0.46±0.43

2.07 0

500GeV 0.580±0.015

0.46±0.43

0.733 0

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1) Mario’s result at 7TeV & 4.93fb-1 (case 3)

Exclusion limit

40

Result for pair production only : 391GeV

Expected Limit : 370 GeVObserved Limit : 393 GeV

preselection, M4l, relIso at 7TeV & 4.983fb-1 (case 1) Mario’s result at 7TeV & 4.93fb-1 (case 3)

41

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 173±2

0.550±0.174

40.6 0

170GeV 106±2

0.529±0.174

30.3 0

200GeV 56.4±0.8

0.486±0.174

20.7 0

250GeV 22.7±0.4

0.454±0.174

11.7 0

300GeV 9.34±0.18

0.442±0.174

6.48 0

350GeV 4.31±0.08

0.437±0.174

3.74 0

400GeV 2.47±0.05

0.434±0.174

2.47 0

500GeV 0.605±0.015

0.431±0.174

0.779 0

# of expected events

# of expected eventsM(H++) H++H– to

eμeμBkg.(MC) ScL Data

150GeV 176±2

0.57±0.43

40.9 0

170GeV 109±2

0.55±0.43

30.7 0

200GeV 56.9±0.8

0.50±0.43

20.8 0

250GeV 21.7±0.3

0.47±0.43

11.3 0

300GeV 9.32±0.18

0.46±0.43

6.42 0

350GeV 4.29±0.08

0.46±0.43

3.69 0

400GeV 2.01±0.04

0.46±0.43

2.07 0

500GeV 0.580±0.015

0.46±0.43

0.733 0

Mario’s result at 7TeV & 4.93fb-1 (case 3)preselection, sumpt, mass window at 7TeV & 4.983fb-1

(case 2)

Exclusion limit

42

Result for pair production only : 391GeV

Expected Limit : 363 GeVObserved Limit : 381 GeV

Mario’s result at 7TeV & 4.93fb-1 (case 3)

preselection, sumpt, mass window at 7TeV & 4.983fb-1 (case 2)

Systematic uncertainties

• Trigger and primary vertex finding : 1.5%

• Signal cross section : 10~15%

• Luminosity (for signal only) : 4.5%

43

Effect of tau

44

45

# of events VS channel(m200)

Gen2 channels are separated by gen leptons matched 4 final reco leptons.

No cut

nlep+>=2, nlep-

>=2

nlep+>=2, nlep-

>=2nlep+>=2, nlep->=2,|mp-mn|/(mp+mn)<0.1

46

# of events VS channel(ZZJetsTo4L)

nlep+>=2, nlep-

>=2

nlep+>=2, nlep-

>=2

nlep+>=2, nlep->=2,|mp-mn|/(mp+mn)<0.1

Gen2 channels are separated by gen leptons matched 4 final reco leptons.

47

# of events VS channel(TTTo2L2Nu2B)

Gen2 channels are separated by gen leptons matched 4 final reco leptons.

nlep+>=2, nlep-

>=2

nlep+>=2, nlep-

>=2

nlep+>=2, nlep->=2,|mp-mn|/(mp+mn)<0.1

48

Cut variables(ZZJetsTo4L)

nlep+>=2, nlep-

>=2

recoch(eeee) &&gen2ch(eeee)

recoch(eeee) &&gen2ch(etet)

recoch(eeee) &&gen2ch(eeee)

recoch(eeee) &&gen2ch(etet)

49

Cut variables(m200)

nlep+>=2, nlep-

>=2

recoch(eeee) &&gen2ch(eeee)

recoch(eeee) &&gen2ch(eeee)

recoch(eeee) &&(gen2ch(eeet)||gen2ch(eett)||gen2ch(etet)||gen2ch(ettt)||gen2ch(tttt))

recoch(eeee) &&(gen2ch(eeet)||gen2ch(eett)||gen2ch(etet)||gen2ch(ettt)||gen2ch(tttt))

effciencygen2ch recoch

M(H++)presel presel &&

|mp-mn|/(mp+mn)<0.1

presel&&

m4>300

presel && |mp-mn|/

(mp+mn)<0.1&& m4>300

etet

eeee

150GeV 0.727%

0.279% 0.388% 0.170%

300GeV 1.09% 0.239% 1.02% 0.205%

500GeV 1.63% 0.260% 1.63% 0.260%

ZZJetsTo4L

0.085%

0.0155% 0.0067%

0.00084%

eeem

150GeV 1.61% 0.569% 0.788% 0.303%

300GeV 2.42% 0.820% 2.36% 0.820%

500GeV 1.82% 0.390% 1.82% 0.390%

ZZJetsTo4L

0.204%

0.0376% 0.0166%

0.00304%

emem

150GeV 0.933%

0.424% 0.594% 0.303%

300GeV 1.16% 0.307% 1.13% 0.307%

500GeV 2.54% 0.975% 2.54% 0.975%

ZZJetsTo4L

0.113%

0.0212% 0.0079%

0.00076% 50

Cut efficiency

effciencygen2ch recoch

M(H++)presel presel &&

|mp-mn|/(mp+mn)<0.1

presel&&

m4>300

presel && |mp-mn|/

(mp+mn)<0.1&& m4>300

mtmt

mmmm

150GeV 1.11% 0.341% 0.584% 0.195%

300GeV 1.53% 0.477% 1.40% 0.443%

500GeV 2.34% 0.895% 2.34% 0.895%

ZZJetsTo4L

0.134%

0.0244% 0.0116%

0.00143%

emmm

150GeV 2.24% 1.06% 1.17% 0.669%

300GeV 3.37% 1.19% 3.17% 1.19%

500GeV 4.06% 1.10% 4.06% 1.10%

ZZJetsTo4L

0.246%

0.0490% 0.0201%

0.00270%

emem

150GeV 1.08% 0.389% 0.547% 0.182%

300GeV 0.988%

0.409% 0.954% 0.409%

500GeV 1.10% 0.344% 1.10% 0.344%

ZZJetsTo4L

0.113%

0.0199% 0.0109%

0.00135%51

Cut efficiency

Analysis at 8TeV & 5.098 fb-

1

52

Used samples/TTTo2L2Nu2B_8TeV-powheg-pythia6/jslee-v1-Summer12_TTTo2L2Nu2B-f16eeefb6ffa56020dc09a3f8c0db443/USER/ZZJetsTo4L_TuneZ2star_8TeV-madgraph-tauola/jslee-v1-Summer12_ZZJetsTo4L-f16eeefb6ffa56020dc09a3f8c0db443/USER/ZZJetsTo2L2Q_TuneZ2star_8TeV-madgraph-tauola/jslee-v1-Summer12_ZZJetsTo2L2Q-f16eeefb6ffa56020dc09a3f8c0db443/

USER/WZJetsTo3LNu_TuneZ2_8TeV-madgraph-tauola/jslee-v1-Summer12_WZJetsTo3LNu-f16eeefb6ffa56020dc09a3f8c0db443/

USER/WWJetsTo2L2Nu_TuneZ2star_8TeV-madgraph-tauola/jslee-v1-Summer12_WWJetsTo2L2Nu-

f16eeefb6ffa56020dc09a3f8c0db443/USER/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball/jslee-v1-Summer12_DYJetsToLL_M-50-

f16eeefb6ffa56020dc09a3f8c0db443/USER/DYJetsToLL_M-10To50filter_8TeV-madgraph/jslee-v1-Summer12_DYJetsToLL_M-10To50-

f16eeefb6ffa56020dc09a3f8c0db443/USER/T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola/jslee-v1-Summer12_T_tW-channel-

f16eeefb6ffa56020dc09a3f8c0db443/USER/T_s-channel_TuneZ2star_8TeV-powheg-tauola/jslee-v1-Summer12_T_s-channel-f16eeefb6ffa56020dc09a3f8c0db443/USER/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola/jslee-v1-Summer12_Tbar_tW-channel-

f16eeefb6ffa56020dc09a3f8c0db443/USER/Tbar_t-channel_TuneZ2star_8TeV-powheg-tauola/jslee-v1-Summer12_Tbar_t-channel-f16eeefb6ffa56020dc09a3f8c0db443/

USER

/DoubleMu/jslee-v1-Summer12_DoubleMu_Run2012A_PromptReco_v1-0c6dd9aa6d8f92e3c92cf94d77ab8809/USER/DoubleMu/jslee-v1-Summer12_DoubleMu_Run2012B_PromptReco_v1-0c6dd9aa6d8f92e3c92cf94d77ab8809/USER/DoubleElectron/jslee-v1-Summer12_DoubleElectron_Run2012A_PromptReco_v1-cbcf422786c4af5101f9b073742156b8/USER/DoubleElectron/jslee-v1-Summer12_DoubleElectron_Run2012B_PromptReco_v1-cbcf422786c4af5101f9b073742156b8/USER

Cross sections : https://twiki.cern.ch/twiki/bin/viewauth/CMS/HZZSamples8TeV

53

HLT efficiencies(Nhlt/Ngen)

54

• Used paths (8TeV):HLT_Mu17_Mu8HLT_Ele17_Ele8HLT_Mu8_Ele17HLT_Mu17_Ele8

• Used paths (7TeV):HLT_DoubleMu7HLT_Mu13_Mu8HLT_Mu17_Mu8HLT_Ele17_Ele8HLT_Mu8_Ele17HLT_Mu17_Ele8

Skim efficiencies(Nskim/Nhlt)

55

• Skim conditions :firstpt>15 GeV && secondpt>10 GeVfor e, μ

Primary vertex

56

# of vertices8TeV, 5.098 fb-1

normalized chi2

8TeV, 5.098 fb-1

nlep>=2, relIso<0.2

# of vertices7TeV, 4.983 fb-1

normalized chi2

7TeV, 4.983 fb-1

# of events at 8TeV & 5.098fb-1

57

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Channel background(MC)

Data

mmmm 19.4±2.8

58.0±7.6

eeee 68.8±6.4

85.0±9.2

emem 105.3±7.7

145.0±12.0

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10, M4l>100

Channel background(MC)

Data

mmmm 16.5±2.6

25.0±5.0

eeee 60.2±6.0

61.0±7.8

emem 85.0±6.9

77.0±8.8

# of events at 8TeV & 5.098fb-1

58

mmmmMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 3.87±0.03

15.0±0.7

9.0±3

nlep+>=2, nlep->=2, m(H±±)>10, in signal region (relIso<0.2 && M4l >200)

± statistical error ± relative error

mmmmMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 5.76±1.60

22.4±1.0

12.0±3.5

nlep+>=2, nlep->=2, m(H±±)>10, in signal region (relIso<0.2 && |mp-mn|/(mp+mn)<0.2)

mmmmMass of H++

(GeV)In signal region

(MC background)In signal region

by sideband method

In signal region(Data)

150 4.69±1.60

18.4±0.8

5.0±2.2

nlep+>=2, nlep->=2, m(H±±)>10, in signal region (relIso<0.2 && |mp-mn|/(mp+mn)<0.1)

Mμμ

59

nm>=2, pt>10, |η|<2.4, global, good, mu_hits_val>20, mu_nstation>0, relIso<0.2, 60<Mμμ<120

lumi : 5.098 fb-1

cross section of DY: 3533 fb

lumi : 3.496 fb-1

cross section of DY : 3533 fb

lumi : 3.496 fb-1

cross section of DYfor 60<Mll<120 : 3352 fb

Standard model cross sections for CMS at 8TeV : https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV

Muon cut variables

60

nm>=2, pt>10, |η|<2.4, global, good, nhits_val>20, nchi2<10, nsation>=1, relIso<0.2, 60<Mμμ<120

Mee

61

ne>=2, pt>10, |η|<2.5, charge consistent, eidLoose, relIso<0.2, 60<Mμμ<120

lumi : 5.098 fb-1

cross section of DY: 3533 fb

lumi : 1.712 fb-1

cross section of DY : 3533 fb

lumi : 1.712 fb-1

cross section of DYfor 60<Mll<120 : 3352 fb

Standard model cross sections for CMS at 8TeV : https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV

Electron cut variables

62

ne>=2, pt>10, |η|<2.5, charge consistent, eidLoose, relIso<0.2, 60<Mμμ<120

Summary

• No signal excess is observed.

• 95% Confidence Level lower limits are set on the H++ mass of 364GeV in the eeee channel, of 393GeV in the eμeμ, of 399GeV in the μμμμ at 7TeV.

63

64

Plan per month

13.2 13.3 13.4 13.5 13.6 13.7 13.8

2012 data pat sample production

8TeV signal sampleproduction

MC VS 2012 data with Z->ll

Lepton ID study

Pileup study

Systematic Study

Paper writing

Defense

Plan

?

backup

65

• Primary vertex

• Cut variables for event selection

• Gen matching effciency for leptons

Primary vertex

66

z position8TeV, 5.098 fb-1

rho position8TeV, 5.098 fb-1

nlep>=2, relIso<0.2

z position7TeV, 4.983 fb-1

rho position7TeV, 4.983 fb-1

Cut variables for event selection(mmmm)

67

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

68

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Cut variables for event selection(mmmm)

69

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Cut variables for event selection(eeee)

70

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Cut variables for event selection(eeee)

71

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Cut variables for event selection(emem)

72

8TeV, 5.098 fb-1

8TeV, 5.098 fb-1

7TeV, 4.983 fb-1

7TeV, 4.983 fb-1

nlep+>=2, nlep->=2, relIso<0.5, m(H±±)>10

Cut variables for event selection(emem)

Minimum one MC background event after cuts from Mario

• This comes from MC uncertainty.

• Basically it shows the upper bound on what amount of events could have survived if you had done more MC so the number comes from limiting MC.

• You cannot claim that you'd see 0 background.

• You overestimate the background strongly and your observed limit would drop far below -2 sigma.

73

74

eeee

M(H++)gench(eeee)

&& recoch(eeee)

gench(eeee)&&

recoch(eeee)&&gen2ch(eee

e)

effciency

150GeV

3832 3816 99.6%

250GeV

2317 2306 99.5%

350GeV

1495 1490 99.7%

500GeV

718 718 100%

Channel matching at 7TeV

mmmm

M(H++)gench(mmmm

)&&

recoch(mmmm)

gench(mmmm)&&

recoch(mmmm)&&gen2ch(mmm

m)

effciency

150GeV

5717 5715 99.9%

250GeV

3318 3316 99.9%

350GeV

2184 2183 99.9%

500GeV

1180 1179 99.9%emem

M(H++)gench(emem)

&& recoch(emem)

gench(emem)&&

recoch(emem)&&gen2ch(eme

m)

effciency

150GeV

4539 4520 99.6%

250GeV

2599 2593 99.8%

350GeV

1753 1751 99.9%

500GeV

957 956 99.9%

Genmatching effciency for leptons at 8TeV

75

Electron cut variables

76

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5

Electron cut variables

77

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5

Electron cut variables

78

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5

Electron cut variables

79

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5, chargeok, |d0|<0.02, eidLoose, nchi2<10, sip<4, relIso<0.2

Electron cut variables

80

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5, chargeok, |d0|<0.02, eidLoose, nchi2<10, sip<4, relIso<0.2

Electron cut variables

81

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5, chargeok, |d0|<0.02, eidLoose, nchi2<10, sip<4, relIso<0.2

Muon cut variables

82

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5, global muons

Muon cut variables

83

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.5, global muons

Muon cut variables

84

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.4, global muons

Muon cut variables

85

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.4, global muons, good muons, nchi2<10, nhits_val>20, nstation>=1, relIso<0.2

Muon cut variables

86

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.4, global muons, good muons, nchi2<10, nhits_val>20, nstation>=1, relIso<0.2

Muon cut variables

87

without gen matchingwith gen matching(ΔR<0.01, charge)

8TeV, m200 // pt>10, |η|<2.4, global muons, good muons, nchi2<10, nhits_val>20, nstation>=1, relIso<0.2