CERN Twiki Users-Lei Li

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

2012 Fall Semester

ZH2llbb

MC_Tuning

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

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.

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:

W_weight=((FE* x-sec)*K_factor*L_data)/N_evnts

N^MC-low=W_weight*N_measured

After running ZFakes.cxx, ZFakes.MC.ttbar.tight.root and ZFakes.MC.ttbar.loose.root are generated for counting N_measured , 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 N^MC-low in this case is 200.

P_T(Gev)	\|eta\|	N_T^data-low	N_L^data-low	fake rate	N_T^MC-low	N_L^MC-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 L_MC(Run 107681)=N_evts/((FE* x-sec)*K_factor)=591372/(1304.30*1.20)=377.83 pb^-1

Luminosity L_data(Run 187815)=24.29 pb^-1

fake rate(MC considered)=N_T^data-low-c_TⁱN_T^MC-low/N_L^data-low-c_LⁱN_L^MC-low

where N^MC-low=N_{MC evts measured}*(L_data/L_MC) and c_Tⁱare scale factors for MC(we did not incluce them here)

P_T(Gev)	\|eta\|	N_T^data-low	N_L^data-low	fake rate	N_T^MC-low	N_L^MC-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

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

P_T(Gev)	\|eta\|	N_T(tight electron)	N_L(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

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=N_T(tight)/N_L(loose)

Here is the table of fake rate for electrons:

P_T(Gev)	\|eta\|	N_T(tight electron)	N_L(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

For loose electrons

10/05/2011

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

For el_trackd0pvunbiased VS el_Etcone20/el_Et

10/04/2011

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

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

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

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

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

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

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