CERN Twiki Users-Lei Li

https://twiki.cern.ch/twiki/bin/view/Main/LeiLi


2014 Summer

VVJJ_HVTsamples

VVJJ_JHUGen

VVJJ_CalcHEP

2012 Fall Semester

ZH2llbb

MC_Tuning

VVanalysis


2012 Summer at CERN

CERN_week1--06/04/2012-06/10/2012

CERN_week2--06/11/2012-06/17/2012

CERN_week3--06/18/2012-06/25/2012

CERN_week4--06/25/2012-07/1/2012

CERN_week5--07/02/2012-07/08/2012

CERN_week6--07/09/2012-07/15/2012

CERN_week7--07/16/2012-07/22/2012

CERN_week8--07/23/2012-07/29/2012


JetMuScanSubstructure_Study(02/15/2012--)

01/30/2012

01/24/2012

01/06/2012

12/16/2011

It turns out that the reason why I got tight more than loose does have something to do with the release. I tried a data11 sample (just one run) switching to realease 17 ( used to be 16) and got the results below which seem better. I will apply release 17 to MC11 samples for the next step.

E_real_Pt.png E_real_eta.png

11/08/2011

Kevin downloaded all of the W->e+nu+jet samples and I used them to calculate the fake rate:

mc10_7TeV.107680.AlpgenJimmyWenuNp0_pt20.merge.NTUP_TOP.e600_s933_s946_r2302_r2300_p572/

mc10_7TeV.107681.AlpgenJimmyWenuNp1_pt20.merge.NTUP_TOP.e600_s933_s946_r2302_r2300_p572/

mc10_7TeV.107682.AlpgenJimmyWenuNp2_pt20.merge.NTUP_TOP.e600_s933_s946_r2302_r2300_p572/

mc10_7TeV.107683.AlpgenJimmyWenuNp3_pt20.merge.NTUP_TOP.e760_s933_s946_r2302_r2300_p572/

mc10_7TeV.107684.AlpgenJimmyWenuNp4_pt20.merge.NTUP_TOP.e760_s933_s946_r2302_r2300_p572/

mc10_7TeV.107685.AlpgenJimmyWenuNp5_pt20.merge.NTUP_TOP.e760_s933_s946_r2302_r2300_p572/

Total number of events processed is 5425861.

To find out how many events are there for each RunNumber, go to: TopD3PDCorrections/TopD3PDCorrections/MCweight.h

This is the formula to calculate the weight:

Wweight=((FE* x-sec)*Kfactor*Ldata)/Nevnts

NMC-low=Wweight*Nmeasured

After running ZFakes.cxx, ZFakes.MC.ttbar.tight.root and ZFakes.MC.ttbar.loose.root are generated for counting Nmeasured , just do like:

root [2] physics_aida->Draw("el_n","(event_weight)*((el_cl_pt[user_topEl_index]>40000)&&(el_cl_pt[user_topEl_index]<50000)&&(el_eta[user_topEl_index]>1.800)&&(el_eta[user_topEl_index]<2.500))")

(Long64_t)4700

root [3] htemp->Integral()

(const Double_t)1.99861166887218133e+02

We can get that NMC-low in this case is 200.

PT(Gev)

|eta|

NTdata-low

NLdata-low

fake rate

NTMC-low

NLMC-low

fake rate(MC considered)

25->30

0.000->1.000

2857

14457

0.198

221

370

0.187

30->35

0.000->1.000

1625

6870

0.237

244

351

0.212

35->40

0.000->1.000

1088

3758

0.290

267

331

0.240

40->50

0.000->1.000

1356

3346

0.405

403

453

0.329

50->infinity

0.000->1.000

1827

3163

0.578

573

597

0.489

25->30

1.000->1.800

1500

13806

0.109

115

299

0.103

30->35

1.000->1.800

902

6282

0.144

131

262

0.128

35->40

1.000->1.800

683

3436

0.199

144

220

0.168

40->50

1.000->1.800

939

3311

0.284

220

280

0.237

50->infinity

1.000->1.800

1148

2743

0.419

298

333

0.353

25->30

1.800->2.500

2202

14726

0.150

100

269

0.145

30->35

1.800->2.500

1335

6864

0.194

128

228

0.182

35->40

1.800->2.500

858

3778

0.227

136

197

0.202

40->50

1.800->2.500

1121

3508

0.320

200

231

0.281

50->infinity

1.800->2.500

1101

2714

0.406

250

284

0.350

11/05/2011

For the MC part, mc10_7TeV.107681.AlpgenJimmyWenuNp1_pt20.merge.NTUP_TOP.e600_s933_s946_r2302_r2300_p572 (W+1parton) is used.

Luminosity LMC(Run 107681)=Nevts/((FE* x-sec)*Kfactor)=591372/(1304.30*1.20)=377.83 pb-1

Luminosity Ldata(Run 187815)=24.29 pb-1

fake rate(MC considered)=NTdata-low-cTiNTMC-low/NLdata-low-cLiNLMC-low

where NMC-low=NMC evts measured*(Ldata/LMC) and cTiare scale factors for MC(we did not incluce them here)

PT(Gev)

|eta|

NTdata-low

NLdata-low

fake rate

NTMC-low

NLMC-low

fake rate(MC considered)

25->30

0.000->1.000

2857

14457

0.198

103

124

0.192

30->35

0.000->1.000

1625

6870

0.237

125

140

0.223

35->40

0.000->1.000

1088

3758

0.290

137

153

0.264

40->50

0.000->1.000

1356

3346

0.405

293

317

0.351

50->infinity

0.000->1.000

1827

3163

0.578

486

511

0.506

25->30

1.000->1.800

1500

13806

0.109

52

75

0.105

30->35

1.000->1.800

902

6282

0.144

62

84

0.136

35->40

1.000->1.800

683

3436

0.199

68

89

0.184

40->50

1.000->1.800

939

3311

0.284

153

187

0.252

50->infinity

1.000->1.800

1148

2743

0.419

248

283

0.366

25->30

1.800->2.500

2202

14726

0.150

41

58

0.147

30->35

1.800->2.500

1335

6864

0.194

51

66

0.189

35->40

1.800->2.500

858

3778

0.227

60

75

0.216

40->50

1.800->2.500

1121

3508

0.320

135

160

0.295

50->infinity

1.800->2.500

1101

2714

0.406

209

234

0.360

CF_tight_MC.png

CF_loose_MC.png

11/02/2011

A much bigger sample is used to ensure more statistics: data11_7TeV.00187815.physics_Egamma.merge.NTUP_TOPEL.r2713_p705_p694_p722_tid528457_00

PT(Gev)

|eta|

NT(tight electron)

NL(loose electron)

fake rate

25->30

0.000->1.000

2857

14457

0.198

30->35

0.000->1.000

1625

6870

0.237

35->40

0.000->1.000

1088

3758

0.290

40->50

0.000->1.000

1356

3346

0.405

50->infinity

0.000->1.000

1827

3163

0.578

25->30

1.000->1.800

1500

13806

0.109

30->35

1.000->1.800

902

6282

0.144

35->40

1.000->1.800

683

3436

0.199

40->50

1.000->1.800

939

3311

0.284

50->infinity

1.000->1.800

1148

2743

0.419

25->30

1.800->2.500

2202

14726

0.150

30->35

1.800->2.500

1335

6864

0.194

35->40

1.800->2.500

858

3778

0.227

40->50

1.800->2.500

1121

3508

0.320

50->infinity

1.800->2.500

1101

2714

0.406

CF_tight_1.png

CF_loose_1.png

10/28/2011

The fake rates are estimated in the sample "data11_7TeV.00187815.physics_Egamma.merge.NTUP_TOPEL.r2713_p705_p694_p722_tid528457_00/NTUP_TOPE L.528457._000079.root.1" with 1. MET<40GeV(CutC8) 2. at least one jet(CutC9) 3. exactly one lepton(CutC4)

Fake Rate=NT(tight)/NL(loose)

Here is the table of fake rate for electrons:

PT(Gev)

|eta|

NT(tight electron)

NL(loose electron)

fake rate

25->30

0.000->1.000

37

176

0.210

30->35

0.000->1.000

25

88

0.284

35->40

0.000->1.000

7

36

0.194

40->50

0.000->1.000

22

39

0.564

50->infinity

0.000->1.000

22

36

0.611

25->30

1.000->1.800

18

155

0.116

30->35

1.000->1.800

10

63

0.159

35->40

1.000->1.800

7

46

0.152

40->50

1.000->1.800

14

40

0.350

50->infinity

1.000->1.800

19

42

0.452

25->30

1.800->2.500

27

186

0.145

30->35

1.800->2.500

16

68

0.235

35->40

1.800->2.500

11

46

0.239

40->50

1.800->2.500

13

50

0.260

50->infinity

1.800->2.500

10

30

0.333

For tight electrons

CF_tight.png

For loose electrons

CF_loose.png

10/05/2011

For el_trackd0pvunbiased/el_tracksigd0pvunbiased(x-axis) VS el_Etcone20/el_Et(y-axis)

count_iso.png

count_noniso.png

count_all.png

ratio_proj.png iso_ratio_proj.png

all_ratio_proj.png all_iso_ratio_proj.png

original.png

d_all.png d_noniso.png d_iso.png

For el_trackd0pvunbiased VS el_Etcone20/el_Et

d0unbiased_proj.png isolation1_proj.png

10/04/2011

iso1.png

iso.png

d0.png

all.png

isolation_proj.png

distribution.png

d0vsiso.png

09/14/2011

Data: user.AyanaArce.mc10_7TeV.105200.merge.NTUP_BOOST.v2.110624134941/ this is a D3PD with more electrons than the old one

EoverPTRT_proj_EoverP.png

EoverPTRT_proj_TRT.png

EoverPnTRTHighTHits.png

EoverPtrackd0.png

EoverPRhad1.png

trackd0TRTHits.png

EoverP.png

EoverP_Integral.png

Cut

Integral for isolated electrons

Integral for non-isolated electrons

Ratio of isolated electrons over non-isolated electrons

(0.5,0.6)

0.0357961

0.0322341

1.1105

(0.5,0.7)

0.0824742

0.0692639

1.19072

(0.5,0.8)

0.212772

0.211221

1.00735

(0.5,0.9)

0.268328

0.277157

0.968142

(0.5,1)

0.318729

0.351014

0.908023

(0.5,1.1)

0.367698

0.414193

0.887744

(0.5,1.2)

0.367698

0.414193

0.887744

(0.5,1.3)

0.439576

0.508779

0.863982

(0.5,1.4)

0.439576

0.508779

0.863982

(0.5,1.5)

0.490836

0.5752

0.853331

(0.5,1.6)

0.516896

0.60249

0.857933

(0.5,1.7)

0.516896

0.60249

0.857933

(0.5,1.8)

0.559851

0.649916

0.861421

(0.5,1.9)

0.559851

0.649916

0.861421

(0.5,2)

0.602234

0.688531

0.874665

(0.5,2.1)

0.619416

0.704441

0.879301

(0.5,2.2)

0.633734

0.719422

0.880894

(0.5,2.3)

0.633734

0.719422

0.880894

(0.5,2.4)

0.648339

0.732903

0.884618

(0.5,2.5)

0.671535

0.756194

0.888046

(0.5,2.6)

0.684422

0.766691

0.892695

(0.5,2.7)

0.69244

0.776571

0.891663

(0.5,2.8)

0.69244

0.776571

0.891663

(0.5,2.9)

0.699026

0.785233

0.890215

trackd0.png trackd0_integral.png

Cut

Integral for isolated electrons

Integral for non-isolated electrons

Ratio of isolated electrons over non-isolated electrons

(0,0.1)

0.176404

0.0874304

2.01764

(0,0.2)

0.201482

0.11681

1.72487

(0,0.3)

0.201482

0.11681

1.72487

(0,0.4)

0.252209

0.17853

1.41269

(0,0.5)

0.284696

0.211945

1.34326

(0,0.6)

0.284696

0.211945

1.34326

(0,0.7)

0.317184

0.2456

1.29147

(0,0.8)

0.399544

0.324408

1.23161

(0,0.9)

0.445426

0.369414

1.20576

(0,1)

0.507552

0.415329

1.22205

(0,1.1)

0.539755

0.445805

1.21074

(0,1.2)

0.539755

0.445805

1.21074

(0,1.3)

0.548589

0.473248

1.1592

(0,1.4)

0.548589

0.473248

1.1592

(0,1.5)

0.550299

0.486131

1.132

(0,1.6)

0.551724

0.490189

1.12553

(0,1.7)

0.551724

0.490189

1.12553

(0,1.8)

0.552294

0.495374

1.1149

(0,1.9)

0.552294

0.495374

1.1149

(0,2)

0.552579

0.498476

1.10854

09/12/2011

Int_EoverP.png

el_Eoverp.png

Cut

Integral for electrons

Integral for background

Integral for hadrons

Ratio of background over electrons

Ratio of hadrons over electrons

(0.8,0.9)

0.108972

0.261452

0.127442

2.39926

1.1695

(0.8,1)

0.160673

0.269403

0.142603

1.67672

0.887532

(0.8,1.1)

0.2085

0.269403

0.142603

1.2921

0.683946

(0.8,1.2)

0.24664

0.288902

0.21276

1.17135

0.862635

(0.8,1.3)

0.277758

0.299947

0.26651

1.07989

0.959507

(0.8,1.4)

0.277758

0.299947

0.26651

1.07989

0.959507

(0.8,1.5)

0.33188

0.313115

0.32491

0.943456

0.978997

(0.8,1.6)

0.35198

0.34421

0.431764

0.977926

1.22667

(0.8,1.7)

0.3722

0.364393

0.476866

0.979025

1.28121

(0.8,1.8)

0.390604

0.391099

0.515809

1.00127

1.32054

(0.8,1.9)

0.405739

0.420216

0.549887

1.03568

1.35527

(0.8,2)

0.422448

0.420216

0.549887

0.994717

1.30167

(0.8,2.1)

0.435888

0.476639

0.604933

1.09349

1.38782

(0.8,2.2)

0.448117

0.476639

0.604933

1.06365

1.34994

(0.8,2.3)

0.448117

0.527786

0.649069

1.17779

1.44844

(0.8,2.4)

0.472575

0.551962

0.667955

1.16799

1.41344

(0.8,2.5)

0.481777

0.551962

0.667955

1.14568

1.38644

(0.8,2.6)

0.491585

0.595325

0.700799

1.21103

1.42559

(0.8,2.7)

0.499334

0.595325

0.700799

1.19224

1.40347

Int_Rhad1.png

el_Rhad1.png

Cut

Integral for electrons

Integral for background

Integral for hadrons

Ratio of background over electrons

Ratio of hadrons over electrons

(0,0.001)

0.045142

0.052054

0.0402019

1.15312

0.890564

(0,0.002)

0.0705531

0.0737543

0.0594218

1.04537

0.842229

(0,0.003)

0.090583

0.0919434

0.0778736

1.01502

0.859693

(0,0.004)

0.090583

0.0919434

0.0778736

1.01502

0.859693

(0,0.005)

0.111211

0.109212

0.0952304

0.982024

0.856306

(0,0.006)

0.130643

0.124964

0.111736

0.956532

0.855281

(0,0.007)

0.147085

0.140141

0.127745

0.952787

0.868512

(0,0.008)

0.168909

0.154397

0.143303

0.914083

0.848407

(0,0.009)

0.186248

0.166983

0.158061

0.896563

0.848659

(0,0.01)

0.186248

0.166983

0.158061

0.896563

0.848659

(0,0.011)

0.201495

0.179512

0.172313

0.890903

0.855172

(0,0.012)

0.21435

0.190909

0.185888

0.890644

0.867219

(0,0.013)

0.23139

0.201539

0.199565

0.870991

0.86246

(0,0.014)

0.244245

0.212341

0.212579

0.869377

0.870351

(0,0.015)

0.2571

0.222011

0.224747

0.86352

0.87416

(0,0.016)

0.2571

0.222011

0.224747

0.86352

0.87416

(0,0.017)

0.27145

0.232353

0.236675

0.85597

0.871893

(0,0.018)

0.285501

0.242292

0.248673

0.848655

0.871006

(0,0.019)

0.295067

0.250235

0.259907

0.848061

0.880839

(0,0.02)

0.308221

0.259214

0.27132

0.841001

0.880277

/09/09/2011

el_Eoverp.png

I_EoverP.png

Cut

Integral for electrons

Integral for background

Integral for hadrons

Ratio of background over electrons

Ratio of hadrons over electrons

If the left side of the cut is fixed(x1=0 in this case), and the right side of the cut is changing, we could get the plot above.

1

Integral for electrons: 0.500545

Integral for background: 0.391099

Ratio of background over electrons: 0.781347

2

Integral for electrons: 0.76232

Integral for background: 0.635883

Ratio of background over electrons: 0.834142

3

Integral for electrons: 0.861121

Integral for background: 0.777846

Ratio of background over electrons: 0.903294

4

Integral for electrons: 0.914033

Integral for background: 0.857869

Ratio of background over electrons: 0.938554

5

Integral for electrons: 0.942608

Integral for background: 0.906341

Ratio of background over electrons: 0.961525

6

Integral for electrons: 0.962586

Integral for background: 0.937761

Ratio of background over electrons: 0.97421

7

Integral for electrons: 0.976026

Integral for background: 0.961589

Ratio of background over electrons: 0.985208

8

Integral for electrons: 0.984744

Integral for background: 0.977335

Ratio of background over electrons: 0.992476

9

Integral for electrons: 0.994309

Integral for background: 0.990694

Ratio of background over electrons: 0.996364

/08/31/2011

electron_Lei.pdf

electron_Lei.pptx

the defination for considering electrons categories:

Hadron fakes: if they do not match a true elecron, muon or tau and truth types do not equal to 0 or 22."el_truth_type!= 11&&el_truth_type!=13&&el_truth_type!=15&&el_truth_type!=0&&el_truth_type!=22"


Non-isolated electrons: if they match a true electron originating from b-mesons or c-mesons."el_truth_mothertype>410&&el_truth_mothertype<440||el_truth_mothertype>510&&el_truth_mothertype<550||el_truth_mothertype>10410&&el_truth_mothertype<10435||el_truth_mothertype>20410&&el_truth_mothertype<20435"


Background electrons: if they match a true electrons coming from Dalitz decays or from a photon."el_truth_mothertype==111||el_truth_mothertype==22"


truthtype0: "el_truth_type==0"


photon: "el_truth_type==22"


08/25/2011


/!\ Edit conflict - your version:


photon: "el_truth_type==22"


08/25/2011


/!\ End of edit conflict


el_Eoverp.png

Integral for electrons:

0.76232

Integral for background:

0.200017

Integral for Hadrons:

0.268072

Ratio of background over electrons:

0.262379

Ratio of hadrons over electrons:

0.351652

cut:

0.8<E/P<1.5


/!\ Edit conflict - other version:


Integral for electrons:

0.76232

Integral for background: 0.200017 Integral for Hadrons: 0.268072 Ratio of background over electrons: 0.262379 Ratio of hadrons over electrons: 0.351652 cut:0.8<E/P<1.5


/!\ Edit conflict - your version:



/!\ End of edit conflict


el_hasconv2.png

08/24/2011

el_Rhad1.png

el_Rhad1Hadrons.png

el_reta.png

el_retaHadrons.png

08/23/2011

el_weta21.png

el_weta2Hadrons1.png

el_wstot1.png

el_wstotHadrons1.png

el_nTRTHighTHits2.png

08/22/2011

el_deltaPhi2.png

el_deltaEta2.png

el_nTRTHits2.png

el_nTRTHitsEta.png

el_trackd0pvovertracksigd0pv2.png

el_deltaEsoverE2.png

el_deltaEsoverEHadrons.png

el_nPixHits2.png

el_nSiHits2.png

el_trackd02.png

el_TRTHighTHitsRatio2.png

08/21/2011

el_trackd0pvovertracksigd0pv.png

el_weta2.png

el_weta2Hadrons.png

Ratio of the energy difference associated with the largest and second largest energy deposit over the sum of these evergies?

el_deltaEsoverE.png

08/19/2011

Shower width in the first compartment

el_wstot.png

el_wstotHadrons.png

el_nPixHis.png

el_nSiHits.png

08/18/2011

el_trackd0.png

Impact parameter of electron track

el_deltaEta.png

Difference between cluster eta and eta of the extrapolated track

el_deltaPhi.png

Difference between cluster phi and phi of the extrapolated track

el_Eoverp.png

el_TRTHighTHitsRatio.png

el_nTRTHIts.png

el_hasconv.png

Conversion bit is 1 if the electron candidate is considered a conversion

08/17/2011

el_phi_leading_trigger.png

08/16/2011

It turns out to be wrong plots :( see plots after 08/18/2011, which are right

Ratio of the energy difference associated with the largest and second largest energy deposit over the sum of these energies

el_deltaEsoverE.png

Shower width

el_wstot.png

08/15/2011

Number of pixel hits for main backgrounds to isolated electrons(I cannot find any isolated electron..seems it does not have any?):

Number of silicon hits:

el_nSiHits.png

Number of TRT hits:

el_nTRTHits.png

Fraction of high threshold TRT hits

el_TRTHighTHitsRatio.png

ratio of E over p:

08/13/2011

el_nPixHits_comparison.png

el_nSiHits_comparison.png

el_nTRTHits_comparison.png

08/11/2011

el_eta_leading_comparison.png

el_phi_leading_comparison.png

el_phiVs_eta_Good.png

el_phiVs_eta_Bad.png

08/08/2011

el_eta_medium_comparison.png

el_phi_medium_comparison.png

el_phiVs_eta_medium_comparison.png

08/07/2011

Monday_Meeting1.pdf

08/04/2011

plots using files made from "user.leli.mc11_7TeV.105015.NTUP_BOOST.e850_s1273_ath023_r2529.full2_v1/" which comes from "user.AyanaArce.mc11_7TeV.105015.AOD.e850_s1273_ath023_r2529.full2/"

el_eta_comparison.gif

el_phi_comparison.gif

el_phiVS_eta_comparison.gif

jet_eta_comparison.gif

jet_phi_comparison.gif

jet_eta_leading_comparison.gif

jet_phi_leading_comparison.gif

PhiVsEtaGood.gif

PhiVsEtaBad.gif

PhiVsEta_leading_Good.gif

PhiVsEta_leading_Bad.gif

PhiVsEta_leading_comparison.gif

jet_eta_E50_comparison.png

jet_phi_E50_comparison.png

08/03/2011

2-D plots:

For Good:

1.All: Comparison_jet_phiVS_eta_GOOD.gif

2.E>50Gev: Comparison_jet_phiVSjet_eta_E50Gev_GOOD.gif

3 Leading: Comparison_jet_phiVSjet_eta_leading_GOOD.gif

For Bad:

1.All: Comparison_jet_phiVS_eta_BAD.gif

2.E>50Gev: Comparison_jet_phiVSjet_eta_E50Gev_BAD.gif

3.Leading: Comparison_jet_phiVS_eta_leading_BAD.gif

LeiLi (last edited 2014-08-15 16:41:16 by LeiLi)