Main Web>TWikiUsers>BingLi>SingleZ (2012-11-12, BingLi)

Single resonance Z->4l analysis

Introduction

The scope of the Z->4l analysis is the measurments of production cross-sections of pp->Z->4l and branching ratios of Z->4l (l=e or mu). This twiki page provides practical information and analysis details for this analysis.

Meetings

The informal weekly Z->4l analysis meetings take place at 5pm on Thursdays (CERN time). Meeting indico agendas can be found at https://indico.cern.ch/categoryDisplay.py?categId=2322

Mailing list

atlas-phys-sm-zsr

Special notes

There are 2 trigger menus for electrons in MC11c D3PD(p956). You may only see trigger variables available for the first trigger menu in some D3PD files and variables only available for the other trigger menu in some other files . In most of the D3PD files, however, trigger variables are available for both trigger menus. This may cause trouble for D3PD skimming and for analysis coding.

2011 versus 2012

There are quite some differences in data taking condition, reconstruction and particle identification between 2011 and 2012.

2011
7 TeV p-p collisions
release 17.0.X
use GSF electrons
use an optimized "loose++" electron ID

2012
8 TeV p-p collisions
more pile up with respect to 2011
release 17.2.X
new trigger menus with respect to 2011
new additional muon chambers (EE) on the Atlas detector with respect to 2012
use a completly new electron recontstruction (GSF electrons by default) (with a lot of optimization on tracking and electron ID)
use an optimized/pile up independant "multilepton" electron ID

2011 and 2012 are treated as two separate analyses and will be combined for Z->4l branching ratio.

Data and MC samples (lists of D3PD files)

The HSG2D3PDs are used for this analysis. Information on filters applied in the D3PD production: here

2012

2012 MC12a: "p1044/"

( MC12a samples list with p1044 )
Some of the MC12a samples with tag 'r3542_r3549_p1044/' are affected by the trigger configuration issue and you should use new digit+reco with tag : 'r3752_r3549_p1044/'
On top of that you should also applied the trigger fix to remove these events for safty : ( trigger fix twiki )

For the MC12a WH samples with tag 'e1217_s1469_s1470_r3752_r3549_p1044' there is configuration issue for the decay that lead to an underestimation of 100 of the cross section and you should use new correct WH samples with tag : 'e1419_s1499_s1504_r3658_r3549_p1044/'

MC12 samples

ttbar
mc12_8TeV.105200.McAtNloJimmy_CT10_ttbar_LeptonFilter.merge.NTUP_HSG2.e1193_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146369.McAtNloJimmy_CT10_ttbar_4LepMass_Mll40GeV12GeV.merge.NTUP_HSG2.e1270_s1469_s1470_r3542_r3549_p1044/
Z+Jets
mc12_8TeV.107650.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp0.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107651.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp1.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107652.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp2.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107653.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp3.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107654.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp4.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107655.AlpgenJimmy_AUET2CTEQ6L1_ZeeNp5.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107660.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp0.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107661.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp1.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107662.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp2.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107663.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp3.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107664.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp4.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107665.AlpgenJimmy_AUET2CTEQ6L1_ZmumuNp5.merge.NTUP_HSG2.e1218_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.107670.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp0.merge.NTUP_HSG2.e1218_s1469_s1470_r3753_r3549_p1044/
mc12_8TeV.107671.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp1.merge.NTUP_HSG2.e1218_s1469_s1470_r3753_r3549_p1044/
mc12_8TeV.107672.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp2.merge.NTUP_HSG2.e1218_s1469_s1470_r3605_r3549_p1044/
mc12_8TeV.107673.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp3.merge.NTUP_HSG2.e1218_s1469_s1470_r3605_r3549_p1044/
mc12_8TeV.107674.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp4.merge.NTUP_HSG2.e1218_s1469_s1470_r3753_r3549_p1044/
mc12_8TeV.107675.AlpgenJimmy_AUET2CTEQ6L1_ZtautauNp5.merge.NTUP_HSG2.e1218_s1469_s1470_r3605_r3549_p1044/
mc12_8TeV.146830.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZeeNp0Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146831.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZeeNp1Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146832.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZeeNp2Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146833.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZeeNp3Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146834.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZeeNp4Excl_Mll10to60.merge.NTUP_HSG2.e1274_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146840.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZmumuNp0Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146841.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZmumuNp1Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146842.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZmumuNp2Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146843.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZmumuNp3Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146844.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZmumuNp4Excl_Mll10to60.merge.NTUP_HSG2.e1274_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146850.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZtautauNp0Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146851.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZtautauNp1Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146852.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZtautauNp2Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146853.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZtautauNp3Excl_Mll10to60.merge.NTUP_HSG2.e1254_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.146854.AlpgenJimmy_Auto_AUET2CTEQ6L1_ZtautauNp4Excl_Mll10to60.merge.NTUP_HSG2.e1348_s1469_s1470_r3542_r3549_p1044/
Zbb
mc12_8TeV.146980.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZeeNp0.merge.NTUP_HSG2.e1315_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146981.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZeeNp1.merge.NTUP_HSG2.e1315_s1499_s1504_r3542_r3549_p1044/
mc12_8TeV.146982.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZbbeeNp2.merge.NTUP_HSG2.e1319_s1499_s1504_r3542_r3549_p1044/
mc12_8TeV.146985.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZmumuNp0.merge.NTUP_HSG2.e1315_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146986.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZmumuNp1.merge.NTUP_HSG2.e1315_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146987.AlpgenJimmy_Auto_AUET2CTEQ6L1_4lFilter_ZbbmumuNp2.merge.NTUP_HSG2.e1319_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146990.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZeeNp0.merge.NTUP_HSG2.e1315_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146991.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZeeNp1.merge.NTUP_HSG2.e1315_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146992.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZbbeeNp2.merge.NTUP_HSG2.e1319_s1499_s1504_r3752_r3549_p1044/
mc12_8TeV.146995.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZmumuNp0.merge.NTUP_HSG2.e1315_s1499_s1504_r3542_r3549_p1044/
mc12_8TeV.146996.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZmumuNp1.merge.NTUP_HSG2.e1315_s1499_s1504_r3542_r3549_p1044/
mc12_8TeV.146997.AlpgenJimmy_Auto_AUET2CTEQ6L1_3lFilter_4lVeto_ZbbmumuNp2.merge.NTUP_HSG2.e1319_s1499_s1504_r3542_r3549_p1044/
WZ
mc12_8TeV.161961.Sherpa_CT10_lllnu_WZ.merge.NTUP_HSG2.e1311_s1469_s1470_r3752_r3549_p1044/
ZZ
mc12_8TeV.116601.gg2ZZJimmy_AUET2CT10_ZZ4e.merge.NTUP_HSG2.e1273_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.116602.gg2ZZJimmy_AUET2CT10_ZZ4mu.merge.NTUP_HSG2.e1273_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.116603.gg2ZZJimmy_AUET2CT10_ZZ2e2mu.merge.NTUP_HSG2.e1273_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.126937.PowhegPythia8_AU2CT10_ZZ_4e_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.126938.PowhegPythia8_AU2CT10_ZZ_2e2mu_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3752_r3549_p1044/
mc12_8TeV.126939.PowhegPythia8_AU2CT10_ZZ_2e2tau_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.126940.PowhegPythia8_AU2CT10_ZZ_4mu_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3752_r3549_p1044/
mc12_8TeV.126941.PowhegPythia8_AU2CT10_ZZ_2mu2tau_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3752_r3549_p1044/
mc12_8TeV.126942.PowhegPythia8_AU2CT10_ZZ_4tau_mll4_2pt5.merge.NTUP_HSG2.e1280_s1469_s1470_r3542_r3549_p1044/
mc12_8TeV.167162.PowhegPythia8_AU2CT10_ZZ_4e_m4l100_150_mll4_4pt3.merge.NTUP_HSG2.e1486_s1499_s1504_r3658_r3549_p1044/
mc12_8TeV.167163.PowhegPythia8_AU2CT10_ZZ_2e2mu_m4l100_150_mll4_4pt3.merge.NTUP_HSG2.e1486_s1499_s1504_r3658_r3549_p1044/
mc12_8TeV.167165.PowhegPythia8_AU2CT10_ZZ_4mu_m4l100_150_mll4_4pt3.merge.NTUP_HSG2.e1486_s1499_s1504_r3658_r3549_p1044/
mc12_8TeV.161988.Sherpa_CT10_llll_ZZ_EW6_noHiggs.merge.NTUP_HSG2.e1346_s1499_s1504_r3658_r3549_p1044/
Other Samles
mc12_8TeV.147806.PowhegPythia8_AU2CT10_Zee.merge.NTUP_HSG2.e1169_s1469_s1470_r3573_r3549_p1044/
mc12_8TeV.147807.PowhegPythia8_AU2CT10_Zmumu.merge.NTUP_HSG2.e1169_s1469_s1470_r3752_r3549_p1044/
mc12_8TeV.161984.Sherpa_CT10_llll_ZZ_EW6.merge.NTUP_HSG2.e1346_s1499_s1504_r3658_r3549_p1044/

For 2012 Data 8 TeV use : "p1044_p1045/"

- Muon Stream : "p1044_p1045/"
- Egamma Stream : "p1044_p1045/"
- Debug Stream : "p1044_p1045/"

[p1084_p1086 is the four lepton filter]

Physics D3PD container for 2012 data: data12_8TeV.current.rtf

Physics D3PD container for 2012 data:

data12_8TeV.periodC8.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC8.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA1.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG2.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD6.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE5.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA7.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG1.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC9.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD2.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC9.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD5.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA5.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA6.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA8.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA8.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodH1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA5.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC2.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA3.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE1.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD3.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodB5.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG4.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD1.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC7.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodH2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD7.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD6.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC3.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE5.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE2.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG5.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC7.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD4.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD8.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC6.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE4.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodH2.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD8.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodH1.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD7.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC2.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE3.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG5.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA6.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD5.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodE.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA4.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG3.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA3.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodD4.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodC6.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodG1.physics_Muons.PhysCont.NTUP_HSG2.grp13_v01_p1044/

data12_8TeV.periodA7.physics_Egamma.PhysCont.NTUP_HSG2.grp13_v01_p1044/

2011

2011 MC11c: "p956/"

( MC11c samples list with p956 )
Some of the MC11c PowHeg samples with tag 'e1180_s1372_s1370_r3108_r3109_p956/' are affected by a bug you should use new one with tag : 'e1271_s1372_s1370_r3108_r3109_p956/'

special D3PD with correct MCTruthClassifier tag : "p1054/"
( MC11c samples list with p1054 )

MC11 samples:

ttbar
mc11_7TeV.105200.T1_McAtNlo_Jimmy.merge.NTUP_HSG2.e835_s1272_s1274_r3043_r2993_p956/
mc11_7TeV.109346.T1_McAtNlo_Jimmy_4LepMass_Mll60GeV12GeV.merge.NTUP_HSG2.e961_s1310_s1300_r3043_r2993_p956/
Z+jets
mc11_7TeV.107650.AlpgenJimmyZeeNp0_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107651.AlpgenJimmyZeeNp1_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107652.AlpgenJimmyZeeNp2_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107653.AlpgenJimmyZeeNp3_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107654.AlpgenJimmyZeeNp4_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107655.AlpgenJimmyZeeNp5_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107660.AlpgenJimmyZmumuNp0_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107661.AlpgenJimmyZmumuNp1_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107662.AlpgenJimmyZmumuNp2_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107663.AlpgenJimmyZmumuNp3_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107664.AlpgenJimmyZmumuNp4_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107665.AlpgenJimmyZmumuNp5_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107670.AlpgenJimmyZtautauNp0_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107671.AlpgenJimmyZtautauNp1_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107672.AlpgenJimmyZtautauNp2_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107673.AlpgenJimmyZtautauNp3_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107674.AlpgenJimmyZtautauNp4_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.107675.AlpgenJimmyZtautauNp5_pt20.merge.NTUP_HSG2.e835_s1299_s1300_r3043_r2993_p956/
mc11_7TeV.116250.AlpgenJimmyZeeNp0_Mll10to40_pt20.merge.NTUP_HSG2.e959_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116251.AlpgenJimmyZeeNp1_Mll10to40_pt20.merge.NTUP_HSG2.e959_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116252.AlpgenJimmyZeeNp2_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116253.AlpgenJimmyZeeNp3_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116254.AlpgenJimmyZeeNp4_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116255.AlpgenJimmyZeeNp5_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116260.AlpgenJimmyZmumuNp0_Mll10to40_pt20.merge.NTUP_HSG2.e959_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116261.AlpgenJimmyZmumuNp1_Mll10to40_pt20.merge.NTUP_HSG2.e959_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116262.AlpgenJimmyZmumuNp2_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116263.AlpgenJimmyZmumuNp3_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116264.AlpgenJimmyZmumuNp4_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116265.AlpgenJimmyZmumuNp5_Mll10to40_pt20.merge.NTUP_HSG2.e944_s1310_s1300_r3043_r2993_p956/
Zbb
mc11_7TeV.116950.AlpgenHWfZeebbNp0_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116951.AlpgenHWfZeebbNp1_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116952.AlpgenHWfZeebbNp2_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116953.AlpgenHWfZeebbNp3_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116955.AlpgenHWfZmumubbNp0_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116956.AlpgenHWfZmumubbNp1_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116957.AlpgenHWfZmumubbNp2_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116958.AlpgenHWfZmumubbNp3_Veto4LepM_Pass3Lep.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116960.AlpgenHWfZeebbNp0_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116961.AlpgenHWfZeebbNp1_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116962.AlpgenHWfZeebbNp2_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116963.AlpgenHWfZeebbNp3_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116965.AlpgenHWfZmumubbNp0_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116966.AlpgenHWfZmumubbNp1_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116967.AlpgenHWfZmumubbNp2_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.116968.AlpgenHWfZmumubbNp3_4LepM.merge.NTUP_HSG2.e835_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.128130.AlpgenJimmyLowMassDYeebbNp0_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128131.AlpgenJimmyLowMassDYeebbNp1_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128132.AlpgenJimmyLowMassDYeebbNp2_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128133.AlpgenJimmyLowMassDYeebbNp3_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128135.AlpgenJimmyLowMassDYmumubbNp0_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128136.AlpgenJimmyLowMassDYmumubbNp1_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128137.AlpgenJimmyLowMassDYmumubbNp2_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128138.AlpgenJimmyLowMassDYmumubbNp3_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128140.AlpgenJimmyLowMassDYtautaubbNp0_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128141.AlpgenJimmyLowMassDYtautaubbNp1_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128142.AlpgenJimmyLowMassDYtautaubbNp2_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128143.AlpgenJimmyLowMassDYtautaubbNp3_nofilter.merge.NTUP_HSG2.e1030_s1372_s1370_r3043_r2993_p956/
WW,WZ
mc11_7TeV.128970.PythiaWW_inclusive.merge.NTUP_HSG2.e988_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.128971.PythiaWZ_inclusive.merge.NTUP_HSG2.e988_s1310_s1300_r3043_r2993_p956/
ZZ
mc11_7TeV.116601.gg2ZZ_JIMMY_ZZ4e.merge.NTUP_HSG2.e1241_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.116602.gg2ZZ_JIMMY_ZZ4mu.merge.NTUP_HSG2.e1241_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.116603.gg2ZZ_JIMMY_ZZ2e2mu.merge.NTUP_HSG2.e1241_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126859.PowHegZZ_4e_trilep5GeV_Pythia.merge.NTUP_HSG2.e1271_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126860.PowHegZZ_4mu_trilep5GeV_Pythia.merge.NTUP_HSG2.e1271_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126861.PowHegZZ_2e2mu_trilep5GeV_Pythia.merge.NTUP_HSG2.e1180_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126862.PowHegZZ_2mu2tau_trilep5GeV_Pythia.merge.NTUP_HSG2.e1180_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126863.PowHegZZ_2e2tau_trilep5GeV_Pythia.merge.NTUP_HSG2.e1180_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.126864.PowHegZZ_4tau_trilep5GeV_Pythia.merge.NTUP_HSG2.e1271_s1372_s1370_r3108_r3109_p956/
mc11_7TeV.128593.PythiaZZall_EF_15_5.merge.NTUP_HSG2.e943_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.161988.Sherpa_CT10_llll_ZZ_EW6_noHiggs.merge.NTUP_HSG2.e1347_s1372_s1370_r3108_r3109_p956/
additional samples
mc11_7TeV.109291.Pythiazz4l.merge.NTUP_HSG2.e972_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.126148.SherpaZZ4l.merge.NTUP_HSG2.e931_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.126148.SherpaZZ4l.merge.NTUP_HSG2.e1039_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.126399.PowHegBoxZZeemm_Pythia_mll025_m4l40.merge.NTUP_HSG2.e995_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.126400.PowHegBoxZZeeee_Pythia_mll025_m4l40.merge.NTUP_HSG2.e995_s1372_s1370_r3043_r2993_p956/
mc11_7TeV.128813.SherpaZZllll.merge.NTUP_HSG2.e950_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.128972.PythiaZZ_inclusive.merge.NTUP_HSG2.e988_s1310_s1300_r3043_r2993_p956/
mc11_7TeV.161984.Sherpa_CT10_llll_ZZ_EW6.merge.NTUP_HSG2.e1347_s1372_s1370_r3108_r3109_p956/

For 2011 Data 7 TeV use : "p997/"

- Muon Stream : "p997/"
- Egamma Stream : "p997/"
- Debug Stream : "p993/"

Physics D3PD container for 2011 data :

data11_7TeV.periodB.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodD.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodE.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodF.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodG.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodH.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodI.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodJ.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodK.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodL.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodM.physics_Egamma.PhysCont.NTUP_HSG2.grp10_v01_p996/

data11_7TeV.periodB.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodD.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodE.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodF.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodG.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodH.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodI.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodJ.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodK.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodL.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/
data11_7TeV.periodM.physics_Muons.PhysCont.NTUP_HSG2.grp10_v01_p996/

MC normalization (cross sections, k-factors, filter efficiencies)

All relevant numbers can be found from here: BkgCrossSection.cxx

GoodRunList

2011

Currently for Data Periods: period B-M + fast reprocessing for few runs with an average total integrated lumonosity of 4583 pb^-1

Channel	Name and Link	Integrated Luminosity pb^-1
All channels	GRL : All_Good.xml	4583

With GRLs "pro10-02" :
You can find the GRLs in /AFS :
/afs/cern.ch/atlas/groups/HSG2/GRL/data11_7TeV.periodAllYear_DetStatus-v36-pro10-02_CoolRunQuery-00-04-08_All_Good.xml

Good Run List twiki page to decide which version of the GRL to use

2012

Currently for Data Periods: period A1-E4 (almost full HCP dataset) with an average total integrated lumonosity of 13022 pb^-1

Channel	Name and Link	Integrated Luminosity pb^-1
All channels	GRL : All_Good.xml	13022

With GRLs "pro13-04" :
You can find the GRLs in /AFS :
/afs/cern.ch/atlas/groups/HSG2/GRL/data12_8TeV.periodAllYear_DetStatus-v53-pro13-04_CoolRunQuery-00-04-08_All_Good.xml

Good Run List twiki page to decide which version of the GRL to use

2012 Data, Luminosity per Period

In this file you can find the luminosities (in μb^-1) for each period contained in the v47-pro13 GRLs (A3 to B14).

Data Period Fraction and PileUpReweighting

Some code for PileUpReweigthing

// Extended Pile up reweighting
Root::TPileupReweighting* PileUpTool = new Root::TPileupReweighting( "PileupReweightingTool" );

2011

string name_PU="";

MC11 c
name_PU="/afs/cern.ch/atlas/groups/HSG2/ExtendedPileUpReweight4l/MC11c.prw.root";

PileUpTool->SetUnrepresentedDataAction(2);
PileUpTool->AddConfigFile(name_PU);
PileUpTool->AddLumiCalcFile("/afs/cern.ch/atlas/groups/HSG2/ExtendedPileUpReweight4l/ilumicalc_2011_AllYear_All_Good.root");
PileUpTool->SetDefaultChannel(109292);
PileUpTool->Initialize();

periodB=PileUpTool->GetIntegratedLumi(177986,178109);
periodD=PileUpTool->GetIntegratedLumi(179710,180481);
periodE=PileUpTool->GetIntegratedLumi(180614,180776);
periodF=PileUpTool->GetIntegratedLumi(182013,182519);
periodG=PileUpTool->GetIntegratedLumi(182726,183462);
periodH=PileUpTool->GetIntegratedLumi(183544,184169);
periodI=PileUpTool->GetIntegratedLumi(185353,186493);
periodJ=PileUpTool->GetIntegratedLumi(186516,186755);
periodK=PileUpTool->GetIntegratedLumi(186873,187815);
periodL=PileUpTool->GetIntegratedLumi(188902,190343);
periodM=PileUpTool->GetIntegratedLumi(190503,191933);

In even loop :
Double_t mc_evt_weight=(*physicsObject->mcevt_weight)[0].at(0);
mu=averageIntPerXing;
double pileupEventWeightPeriofFractionWeight = PileUpTool->GetCombinedWeight(RunNumber,mc_channel_number,mu);

MC11 c
period for H->4l analysis	Run number in MC11c	data period represented	conditions	fraction in MC11c	fraction in data for period B-M(pro10)
B-D	180164	B-D	all FEBs ok	3.2%	3.88 %
E-H	183003	E-H	6 missing FEBs	17.4%	20.63 %
I-K	186169	I-K	4 FEBs restored	25.8%	25.72 %
L-M	189751	L-M	4 FEBs restored	53.5%	49.77 %

2012

MC12a
name_PU="/afs/cern.ch/atlas/groups/HSG2/ExtendedPileUpReweight4l/MC12a.prw.root";

PileUpTool->SetUnrepresentedDataAction(2);
PileUpTool->AddConfigFile(name_PU);
PileUpTool->AddLumiCalcFile("/afs/cern.ch/atlas/groups/HSG2/ExtendedPileUpReweight4l/ilumicalc_2012_AllYear_All_Good.root");
PileUpTool->SetDefaultChannel(160156);
PileUpTool->Initialize();

period2012A=PileUpTool->GetIntegratedLumi(200804,201556);
period2012B=PileUpTool->GetIntegratedLumi(202660,205113);
period2012C=PileUpTool->GetIntegratedLumi(206248,207397);
period2012D=PileUpTool->GetIntegratedLumi(207447,209025);
period2012E=PileUpTool->GetIntegratedLumi(209074,210308); // 2012 period E5

In even loop :
Double_t mc_evt_weight=( physicsObject->mcevt_weight)[0].at(0);
mu=averageIntPerXing;
double pileupEventWeightPeriofFractionWeight = PileUpTool->GetCombinedWeight(RunNumber,mc_channel_number,mu);

MC12 a
period for H->4l analysis	Run number in MC12a	data period represented	conditions
AllYear	195847	AllYear	-

Event Selection

Variables used in invariant mass computation

For the Z mass, use the value from PDG (2010) : 91187.6 MeV

Invariant mass variable recommendations:
2011
Muons:	`(mu__pt)`, (mu__eta), (mu__phi), (mu__E) .
Electrons :	`(el_GSF_Et)`, (el_GSF_tracketa), (el_GSF_trackphi), (el_GSF_cl_E), where `(el_GSF_Et) = (el_GSF_cl_E)/cosh((*el_GSF_tracketa))`
2012
Muons:	`(mu__pt)`, (mu__eta), (mu__phi), (mu__E) .
Electrons :	`(el_Et)`, (el_tracketa), (el_trackphi), (el_cl_E), where `(el_Et) = (el_cl_E)/cosh((*el_tracketa))`

Selection

Event Preselection
GRL	[data only] apply All_Good GRL for ALL final state (4e,4μ,2e2μ). GoodRunLists
LAr error	[data only] remove all events that have LAr EventInfo error flag (larError==2 in D3PD) (data corruption). LAr Event Info twikipage
Rejection of bad / corrupted events	[data only] remove all events that have tileError==2 in D3PD (data corruption). Rejection of bad / corrupted events
Rejection of bad / corrupted events	[data only] remove all events that have coreFlags&0x40000!=0 in D3PD (data corruption).Rejection of bad / corrupted events
PileUp Reweighting	2011-2012 : apply pile up reweighting look table
<strike> Mu Rescaling	<strike>2012 : apply mu rescaling look table
Vertex	at least one reconstructed vertex with 3 associated tracks
z vertex	2012 : apply z vertex reweighting look table
Trigger	lowest p_T single or di lepton unprescaled triggers : 4mu single OR di muon trigger (see the Table) 2e2mu single OR di muon trigger OR single OR di electron trigger (see the Table) 4e single OR di electron trigger (see the Table )
Physics Object Selection
electron smearing and calibraion are applied to all electrons (except for the isLoosePlusPlus/Multilepton macro (see below)
2012 : NO electron crack calibration for 2012 data and MC12a
muon smaring is applied to all muons
Muons (MCP group recommendation link for 2011 )
(MCP group recommendation link for 2012)
μ-ID	Use Staco CB+ST muons (mu_staco_)) ((author==6 or author==7) and isStandAloneMuon==0) Use CaloTag muons (mu_calo_) (author==16) and (`CaloMuonIDTag > 10 or CaloLRLikelihood > 0.9)` Use Staco StandAlone muons (mu_staco_*) (author==6 and isStandAloneMuon==1)
Corrections	2011 MC : smearing 2012 MC : smearing
Kinematics	Staco CB+ST muons p_T>4 GeV, fabs( η)<2.7 CaloTag muons p_T>15 GeV, fabs( η)<0.1 Staco StandAlone muons p_T>4 GeV, fabs( η)>2.5, fabs( η)<2.7
ID/MS hits	Staco CB+ST muons : apply the selection shown here: 2011 Table 2012 Table <nop>CaloTag</nop> muons : apply the selection shown here: 2011 Table 2012 Table Staco <nop>StandAlone</nop> muons : apply the selection shown here: 2011Table 2012 Table
Cosmic Cut	fabs(d0)< 1 mm `(*mu_d0_exPV[i])` (where d₀ is relative to the primary vertex) NOT for StandAloneMuons
Z0 Cut	fabs(z0)<10. mm `(*mu_z0_exPV[i])` NOT for StandAloneMuons
overlap removal	reject CaloMuon (after muon ID cuts) if a Staco (after muon ID cuts) muon is within ΔR < 0.1 (use mu__id_[theta, phi] to compute ΔR) reject StandaloneMuon (after muon ID cuts) if a Staco Segment Tag (after muon ID cuts) muon is within ΔR < 0.2 (use mu__[eta, phi] to compute ΔR)
Electrons (Egamma grouprecommendation link)
(Egamma group recommendation link for 2012)
Algorithm	2011 (author==1 or 3) with GSF and isLoosePlusPlusH4l e-ID (look table) 2012 (author==1 or 3) with GSF (default and only algorithm) and Multilepton e-ID (look table)
Corrections	2011 MC : crack calibration + smearing Data : crack calibration + calibration 2012 MC : smearing Data : calibration
Kinematics	E_T> 7 GeV (EM cluster energy/track direction), fabs(η_Cluster)<2.47
ObjectQuality	(el_GSF_OQ&1446)==0
z₀ Cut	fabs(z0)<10. mm, (*el_GSF_trackz0pvunbiased[i])
Overlap Removal
e-e	Electrons (after electron ID cuts) sharing the same ID-track, keep the electron with highest E_T on the cluster. 2011 : el_*_track[d0, z0, theta, phi, qoverp] 2012 : `el_Unrefittedtrack[d0, z0, theta, phi, qoverp]`
e-e	Electrons (after electron ID cuts) to close each other in eta and phi, keep the electron with highest E_T on the cluster. 2012 : use el_cl_eta,el_cl_phi if fabs(Δeta)< 30.025 && fabs(Δphi) < 50.025 keep the electron with highest E_T on the cluster Δphi = (fabs(phi_1-phi_2) > TMath::Pi()) ? 2*TMath::Pi()-fabs(phi_1-phi_2) : fabs(phi_1-phi_2)
e-μ	Electrons (after electron selection cuts) sharing the same ID-track with a muon candidate (within ΔR<0.02) are removed. 2011 : el__track[eta, phi] and mu__id_[theta, phi]. η = -ln[tan[θ/2]] 2012 : el_Unrefittedtrack[eta, phi] and mu_*_id_[theta, phi]. η = -ln[tan[θ/2]]
e-jets	Jets that overlap with electrons (after electron selection cuts) within ΔR < 0.2 are removed
llll selection
Number of Selected Leptons	4μ : At least 4 selected muons, 4e: At least 4 selected electrons, 2e2mu : At least 2muons and at least 2electrons.
Quadruplet Selection (DILS)	Two same flavour (SF) and opposite sign (OS) lepton pairs. Within a quadruplet, the SFOS pairing that give the mass closest to the Z-boson is considered the primary di-lepton, and the other the secondary di-lepton. The analysis to be split in 4 final states : 2μ2μ, 2e2μ, 2μ2e, 2e2e - where the primary di-lepton is mentioned first. We allow maximum one 1 StandAlone muon OR 1 Calo muon per quadruplet
Kinematics	The leptons in the quaduplet should satisfy the requirements _p_T1>20 GeV, p_T2>15 GeV, p_T3>10 GeV _
Trigger Kinematics and Matching	matching at least one (for single lepton triggers) or two (for dilepton triggers) of the leptons (see table and tool for the configuration) With p956 and p997 some additional lines are needed (see TriggerMatching for p956 and p997) The leptons passed to the trigger matching tool should satisfy the corresponding trigger thresholds. Just one of the chains needs to be matched ,if more than one has fired, for the event to pass. 2011 : use the p_T / E_T theshold listed here detaits 2012 : no p_T / E_T thehsold
Primary Di-Lepton Mass	50 GeV < m₁₂ < 106 GeV
Secondary Di-Lepton Mass	m_{l₃ l₄ >5 (or 0)} GeV
ΔR of leptons in the quadruplet	ΔR(l_i,l_j)>0.10 same flavor ΔR(l_i,l_j)>0.20 opposite flavor
J/ψ veto	not decided now
Selection of one quadruplet	Select only one quadruplet of leptons per event : We select the quadruplets which have the primary di-lepton pair (m_{l₁ l₂}) closest to on-shell Z, we next look for the 'most energetic' off-shell Z (m_{l₃ l₄}), which corresponds to the highest off-shell Z mass
Track Isolation	All leptons of the quadruplet relative track isolation (divided by Et for electron and by pt for muon) 2011 : [cone ΔR<0.20] <15%. 2012 : [cone ΔR<0.20] <15%. Contribution of overlapping leptons from the quadruplet is removed for ΔR<0.20 using el__tracketa,el__trackphi for electrons and mu__eta and mu__phi for muons. For overlapping electron you should substract el__trackpt and for muon 1/fabs(mu__id_qoverp)sin(mu__id_theta) (Since there is no ID track for StandAloneMuon, we don't remove overlapping StandAloneMuon)
Calo Isolation	All leptons of the quadruplet relative calo isolation (divided by Et for electron and by pt for muon) 2011 : [cone ΔR<0.20] < 30% for muons and electrons BUT [cone ΔR<0.20] <15% for StandAlone muon 2012 : [cone ΔR<0.20] < 20% for electron, [cone ΔR<0.20] < 30% for muons BUT [cone ΔR<0.20] <15% for StandAlone muon The contribution of an electron overlapping with another leptons is removed for ΔR<0.18 in all channels with electron : 4e and 2e2mu [see section 5.2.1 of the note] 2011 : For electron : Pile-Up correction (Etcone) : GetNPVCorrectedIsolation (look table) and in case of overlaping electron you need to substract : el_cl_E/cosh(el_tracketa) and for ΔR : use el__tracketa,el__trackphi 2012 : For electron : Pile-Up correction (topo Etcone) : GetPtEDCorrectedTopoIsolation (look table) and in case of overlaping electron you need to substract : el_rawcl_E/cosh(el_tracketa) and for ΔR : use el_etas2,el_phis2 <strike> For muon : Pile-Up correction :- CorrectEtCone </strike>
Impact Parameter Significance	2011 : All leptons of the quadruplet (this cut has no impact on StandAlone muons and is safe) only if m_4l < 190 GeV : 2012 : All leptons of the quadruplet (this cut has no impact on StandAlone muons and is safe) : For muons fabs( d₀/σ(d₀))<3.5 applied to all muons For electrons fabs( d₀/σ(d₀))<6.5 applied to all electrons Recommended variables : `el_trackd0pvunbiased` & `el_tracksigd0pvunbiased` and `mu__trackIPEstimate_sigd0_unbiasedpvunbiased` and`mu__trackIPEstimate_d0_unbiasedpvunbiased` 2012 MC12a :There is also an transverse impact parameter bias to be corrected in MC by subracting 2 μm off the measured value of the d0 : for electrons `el_trackd0pvunbiased-2.e-3` for muons `mu__trackIPEstimate_d0_unbiasedpvunbiased-2.e-3` 2012 MC12a : Also for MC it is recommended that IP is smeared to take into account detector mislaignment effects. Please look at this wiki for details. random number : EventNumber + 100 muon (electron) index eta and pt : from the track execpt for StandAlone Muon (mu_eta, mu_pt)

Trigger details

2011

Trigger details for trigger matching and trigger efficiency can be found in the table below and twiki

2011 MC Period	Fraction with this Trigger	Data Period	Single electron	threshold (matching)	Di electron	threshold (matching)	Single muon	thresold (matching)	Di muon	threshold (matching)
Period B-D	100%	Period B-D	e20_medium	e 21	2e12_medium	2 e 13	EF_mu18_MG	mu 20	EF_2mu10_loose	2 mu 12
Period E-H	100%	Period E-H	e20_medium	e 21	2e12_medium	2 e 13	EF_mu18_MG	mu 20	EF_2mu10_loose	2 mu 12
Period I-K	31.85%	Period I	e20_medium	e 21	2e12_medium	2 e 13	EF_mu18_MG	mu 20	EF_2mu10_loose	2 mu 12
	18.67%	Period J	e20_medium	e 21	2e12_medium	2 e 13	EF_mu18_MG_medium	mu 20	EF_2mu10_loose	2 mu 12
	49.48%	Period K	e22_medium	e 23	2e12T_medium	2 e 13	EF_mu18_MG_medium	mu 20	EF_2mu10_loose	2 mu 12
Period L-M	100%	Period L-M	e22vh_medium1	e 23	2e12Tvh_medium	2 e 13	EF_mu18_MG_medium	mu 20	EF_2mu10_loose	2 mu 12

The fraction for the MC period I-K were we have 2 different trigger for electron and muon is computed using the lumi for period I,J or K /period I+J+K.
period I,J or K luminosty is retrived by the new pile up reweighting tool
Example for handling both electron and muon cases :

 
#include <TRandom3.h>
mutable TRandom3   random3;
//// For MC only
seed=physicsObject->mc_channel_number*physicsObject->EventNumber
random3.SetSeed(seed);
Double_t rd=random3.Uniform();
Double_t frac1=periodI/(periodI+periodJ+periodK);
Double_t frac2=(periodI+periodJ)/(periodI+periodJ+periodK);

if (rd<frac1){
 "trigger_mu18_MG";
}else {
 "trigger_mu18_MG_medium";
}
if (rd<frac2){
 "trigger_e20m"; 
 "trigger_2e12m";
}else {
 "trigger_e22m";
 "trigger_2e12Tm";
}

2011Trigger matcing for D3PD tag p956 and p997

Trigger matching details for trigger matching and trigger efficiency can be found in the table below and twiki
BUT
for D3PD tags p956 and p997
the following need to be set only for muon trigger matching

for "EF_2mu10_loose"

 
triggerNavigationVariables->set_trig_RoI_EF_mu_TrigMuonEFInfoContainerStatus(physicsObject->trig_RoI_EF_mu_TrigMuonEFInfoContainer_MuonEFInfoStatus);
triggerNavigationVariables->set_trig_RoI_EF_mu_TrigMuonEFInfoContainer(physicsObject->trig_RoI_EF_mu_TrigMuonEFInfoContainer_MuonEFInfo);
HERE called the trigger matching function for dimuon trigger

for "EF_mu18_MG" or "EF_mu18_MG_medium"

 
triggerNavigationVariables->set_trig_RoI_EF_mu_TrigMuonEFInfoContainerStatus(physicsObject->trig_RoI_EF_mu_TrigMuonEFInfoContainer_eMuonEFInfoStatus);
triggerNavigationVariables->set_trig_RoI_EF_mu_TrigMuonEFInfoContainer(physicsObject->trig_RoI_EF_mu_TrigMuonEFInfoContainer_eMuonEFInfo);
HERE called the trigger matching function for single muon trigger

2012

Trigger details for trigger matching and trigger efficiency can be found in the table below and twiki

2012 MC Period	Fraction with this Trigger	Data Period	Single electron	threshold (matching) not applied	Di electron	threshold (matching) not applied	Single muon	threshold (matching) not applied	Di muon	threshold (matching) not applied	<strike>electron-muon	threshold (matching) not applied
PeriodAllYear	100%	Period A1-E5	e24vhi_medium1 OR e60_medium1	e 25 OR e 61	2e12Tvh_loose1	2 e 13	mu24i_tight OR mu36_tight	mu 25 OR mu 37	2mu13 OR mu18_tight_mu8_EFFS	2 mu 14 OR mu 19 mu 9	<strike>e12Tvh_medium1_mu8 OR e24vhi_loose1_mu8	e 13 mu 11 OR e 25 mu 9

2012 Trigger matching

For electron and muon trigger matching, several methods such as matchDimuon, matchDimuon, and match are provided in TrigMuonEfficiency (the tag to be used for this package is given below. In 2012, we use both symmetric and asymmetric pT threshold dimuon triggers. In order to see whether a given di-muon object (muon1 and muon2) matches to a dimuon trigger item, one should first call the function

bool matchDimuon(const TLorentzVector& muon1,
                        const TLorentzVector& muon2,
                        const std::string& chain,
                        std::pair<bool, bool>& result1,
                        std::pair<bool, bool>& result2);

and then require

// the first condition is for muon1 as seeded, muon2 as full scan
// the second condition is for muon2 as seeded, muon1 as full scan
bool ismatched = ((result1.first && result2.second) || (result2.first && result1.second));

The higher pT threshold should be applied on the muon that acts the seed. For example adding the pT threshold requirement on EF_mu18_tight_mu8_EFFS trigger item should be done

  
bool ismatched_and_passptthreshold =  (  ( result1.first && result2.second && ( muons1.Pt() > 19000.0 ) && ( muons2.Pt() > 9000.0 ) )
                                                   || ( result2.first && result1.second && ( muons1.Pt() > 9000.0 ) && ( muons2.Pt() > 19000.0 ) ) )

Single electron, single muon or di-electron matching should be done as in 2011.

ID/MS Requirements

2011

ID/MS Hits for 2011

ID/MS Hits for 2011

cuts to apply for Staco/Muid CB+ST muons

!expectBLayerHit || nBLHits > 0
nPixHits + nPixelDeadSensors > 1
nSCTHits + nSCTDeadSensors ≥ 6
nPixHoles + nSCTHoles < 3

n = nTRTHits + nTRTOutliers
if (|η| < 1.9): pass if n > 5 && nTRTOutliers < n*0.9
else if n > 5: pass if nTRTOutliers < n*0.9
else: pass
An example of the TRT + Outliers cut:

   Int_t n = (*mu_staco_nTRTHits)[mu_i] + (*mu_staco_nTRTOutliers)[mu_i];
   Double_t mu_eta = fabs((*mu_staco_eta)[mu_i]);
   Bool_t case1 = ((n > 5) && (*mu_staco_nTRTOutliers)[mu_i] < .9*n) ; //require n>5 && the # of trt outliers < .9n
   Bool_t case2 = (n > 5) ? ((*mu_staco_nTRTOutliers)[mu_i] < .9*n) : true ; //if n>5, we require the # of trt outliers < .9n
   Bool_t trt_ext = ( mu_eta < 1.9) ? case1 : case2 ; // if |eta| < 1.9, then case 1; if |eta|>=1.9, then case 2

cuts to apply for CaloTag muons !expectBLayerHit || nBLHits > 0
nPixHits + nPixelDeadSensors > 1
nSCTHits + nSCTDeadSensors ≥ 6
nPixHoles + nSCTHoles < 3

n = nTRTHits + nTRTOutliers
if (|η| < 0.1): pass if n < 6 || nTRTOutliers < n*0.9

cuts to apply for StandaAlone muons Int_t mu_cscetahits=(*mu_staco_nCSCEtaHits)[mu_i];
Int_t mu_cscphihits=(*mu_staco_nCSCPhiHits)[mu_i];
Int_t mu_emhits=(*mu_staco_nMDTEMHits)[mu_i];
Int_t mu_eohits=(*mu_staco_nMDTEOHits)[mu_i];
mu_cscetahits+mu_cscphihits)>0 && mu_emhits>0 && mu_eohits>0

2012 2012

ID/MS Hits for 2012

ID/MS Hits for 2012

cuts to apply for Staco/Muid CB+ST muons

!expectBLayerHit || nBLHits > 0
nPixHits + nPixelDeadSensors >0 (was 1 in 2011)
nSCTHits + nSCTDeadSensors > 4 (was 5 in 2011)
nPixHoles + nSCTHoles < 3

n = nTRTHits + nTRTOutliers
if (0.1<|η| < 1.9): require if n > 5 && nTRTOutliers < n*0.9
else if (|η| < 0.1 or |η| ≥1.9): if n > 5 require nTRTOutliers < n*0.9
An example of the TRT + Outliers cut:

   Int_t n = (*mu_staco_nTRTHits)[mu_i] + (*mu_staco_nTRTOutliers)[mu_i];
   Double_t mu_eta = fabs((*mu_staco_eta)[mu_i]);
   Bool_t case1 = ((n > 5) && (*mu_staco_nTRTOutliers)[mu_i] < .9*n) ; //require n>5 && the # of trt outliers < .9n
   Bool_t case2 = (n > 5) ? ((*mu_staco_nTRTOutliers)[mu_i] < .9*n) : true ; //if n>5, we require the # of trt outliers < .9n
   Bool_t trt_ext = ( 0.1 < mu_eta < 1.9) ? case1 : case2 ; // if |eta| < 1.9, then case 1; if |eta|>=1.9 or |eta| <0.1 , then case 2

cuts to apply for CaloTag muons !expectBLayerHit || nBLHits > 0
nPixHits + nPixelDeadSensors > 0
nSCTHits + nSCTDeadSensors > 4
nPixHoles + nSCTHoles < 3

n = nTRTHits + nTRTOutliers
if (|η| < 0.1): pass if n < 6 || nTRTOutliers < n*0.9

More information

Package version

Tag	svn Path
`GoodRunsLists-00-01-03`	svn+ssh://svn.cern.ch/reps/atlasoff/DataQuality/GoodRunsLists/tags/
`PileupReweighting-00-02-08`	svn+ssh://svn.cern.ch/reps/atlasoff/PhysicsAnalysis/AnalysisCommon/PileupReweighting/tags/
`MuonMomentumCorrections-00-07-01`	svn+ssh://svn.cern.ch/reps/atlasoff/PhysicsAnalysis/MuonID/MuonIDAnalysis/MuonMomentumCorrections/tags/
`MuonEfficiencyCorrections-02-01-05`	svn+ssh://svn.cern.ch/reps/atlasoff/PhysicsAnalysis/MuonID/MuonIDAnalysis/MuonEfficiencyCorrections/tags/
`MuonIsolationCorrection-00-08`	svn+ssh://svn.cern.ch/reps/penn/schae/MuonIsolationCorrection/tags/
`egammaAnalysisUtils-00-03-61`	svn+ssh://svn.cern.ch/reps/atlasoff/Reconstruction/egamma/egammaAnalysis/egammaAnalysisUtils/tags/
`HiggsZZ4lUtils-00-01-25`	svn+ssh://svn.cern.ch/reps/atlasoff/PhysicsAnalysis/HiggsPhys/HSG2/HiggsZZ4lUtils/tags/
`TrigMuonEfficiency-00-02-17`	svn+ssh://svn.cern.ch/reps/atlasoff/Trigger/TrigAnalysis/TrigMuonEfficiency/tags/ + 'muon_trigger_sf_2012_AtoE.root' is available here /afs/cern.ch/user/t/takashi/public/SF2012/ + 'muon_trigger_sf_mc11c.root' is available /afs/cern.ch/user/t/takashi/public/muon_trigger_sf_mc11c.root

Package usage

Object	goal	when to use it	code	function	package	when to apply it	seed of needed
Muon	*pT resolution (smearing)*	2011	SmearingClass	SmearingClass("Data11","staco" or "muid","q_pT","Rel17","MuonMomentumCorrections/share/")	`MuonMomentumCorrections`	MC	SetSeed(EventNumber, i);
	*pT resolution (smearing)*	2012	SmearingClass	SmearingClass("Data12","staco" or "muid","q_pT","Rel17.2","MuonMomentumCorrections/share/")			SetSeed(EventNumber, i);
			2011-2012	UseScale(1)			i index within mu_ branches
			2011-2012	UseImprovedCombine()
			2011-2012	Event(ptms,ptid,pT,eta) if mu_isCombinedMuon
			2011-2012	Event(pT,eta,"ID") if mu_isSegmentTaggedMuon
			2011-2012	_Event(pT,eta,"ID") for CaloMuons _
			2011-2012	Event(pT,eta,"MS") if isStandAloneMuon
	*Muon efficiencies (all)*	2011	MuonEfficiencyCorrections	config=Analysis::AnalysisMuonConfigurableScaleFactors::AverageOverPeriods; config_sa=Analysis::AnalysisMuonConfigurableScaleFactors::Default; std::string muon_type_staco ("STACO_CB_plus_ST_2011_SF.txt"); std::string muon_type_staco_sa("STACOHighEta.txt"); std::string muon_type_muid ("Muid_CB_plus_ST_2011_SF.txt"); std::string muon_type_muid_sa ("MuidHighEta.txt"); std::string muon_type_calo_mu ("CaloTag_2011_SF.txt"); MuidSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_muid,unit,config); MuidSCF->Initialise(); MuidSASCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_muid_sa,unit,config_sa); MuidSASCF->Initialise(); StacoSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_staco,unit,config); StacoSCF->Initialise(); StacoSASCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_staco_sa,unit,config_sa); StacoSASCF->Initialise(); CaloMuSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_calo_mu,unit,config); CaloMuSCF->Initialise(); Calo/StatoSA/MuidSASCF->scaleFactor(tlv); Staco/MuidSCF->scaleFactor(tlv);	`MuonEfficiencyCorrections`	MC	n.a
	*Muon efficiencies (all)*	2012	MuonEfficiencyCorrections	config=Analysis::AnalysisMuonConfigurableScaleFactors::AverageOverRuns; config_sa=Analysis::AnalysisMuonConfigurableScaleFactors::AverageOverRuns; std::string muon_type_staco ("STACO_CB_plus_ST_2012_SF.txt"); std::string muon_type_staco_sa("STACO_CB_plus_ST_2012_SFms.txt"); std::string muon_type_muid ("Muid_CB_plus_ST_2012_SF.txt"); std::string muon_type_muid_sa ("Muid_CB_plus_ST_2012_SFms.txt"); std::string muon_type_calo_mu ("CaloTag_2012_SF.txt"); MuidSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_muid,unit,config); MuidSCF->Initialise(); MuidSASCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_muid_sa,unit,config_sa); MuidSASCF->Initialise(); StacoSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_staco,unit,config); StacoSCF->Initialise(); StacoSASCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_staco_sa,unit,config_sa); StacoSASCF->Initialise(); CaloMuSCF=new Analysis::AnalysisMuonConfigurableScaleFactors (directory,muon_type_calo_mu,unit,config); CaloMuSCF->Initialise(); Muid/StacoSCF->scaleFactor(charge,tlv); CaloSCF->scaleFactor(tlv);	`MuonEfficiencyCorrections`	MC	n.a
	*Calorimeter Isolation*	not used	CorrectCaloIso	CorrectEtCone(mu_etcone20,nPV(n_tracks>=3),mu_eta,"cone20Comb")	`MuonIsolationCorrection`	Data and MC	n.a
	*Trigger-Offline matching*	2011-2012	TrigMuonEfficiency	muonTriggerMatchTool->match(mu_eta, mu_phi,chain);pt>thresold	`TrigMuonEfficiency`	Data and MC	n.a
				muonTriggerMatchTool->matchDimuon(tlv_mu1, mu2,chain,res1,res2);res1.first && res2.first&&pt1>thresold1 &&pt2>thresold2;
	*muon efficiencies (trigger)*	2011	TrigMuonEfficiency	leptonSF = new LeptonTriggerSF(2011, "TrigMuonEfficiency/share/", "muon_trigger_sf_mc11c.root"); runnumber_sf=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber);p=loosepp;q=loose; pair SF = SF = leptonSF->GetTriggerSF(runnumber_sf,false, muons, q, electrons, p, 0);SF.first	`TrigMuonEfficiency`	MC	n.a
	*muon efficiencies (trigger)*	2012	TrigMuonEfficiency	leptonSF = new LeptonTriggerSF(2012, "TrigMuonEfficiency/share", "muon_trigger_sf_2012_AtoE.root "); /afs/cern.ch/user/t/takashi/public/SF2012/muon_trigger_sf_2012_AtoD.root runnumber_sf=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber);p=ML;q=loose; pair SF = SF = leptonSF->GetTriggerSF(runnumber_sf,false, muons, q, electrons, p, 0);SF.first	`TrigMuonEfficiency`	MC	n.a
Electron	*Energy correction initialisation*	2011	EnergyRescalerUpgrade	Init("egammaAnalysisUtils/share/EnergyRescalerData.root","2011","es2011a");	`egammaAnalysisUtils`	Data/MC	n.a
		2012	EnergyRescalerUpgrade	Init("egammaAnalysisUtils/share/EnergyRescalerData.root","2012","es2012");	`egammaAnalysisUtils`	Data/MC	n.a
	*Energy correction (scale)*	2011-2012	EnergyRescalerUpgrade	applyEnergyCorrection(Eta, E, egRescaler::EnergyRescalerUpgrade::Electron, egRescaler::EnergyRescalerUpgrade::Nominal,1.0,RunNumber);	`egammaAnalysisUtils`	Data	n.a
	*Energy correction for the crack (scale)*	2011	EnergyRescalerUpgrade	applyMCCalibration(Eta, Et, egRescaler::EnergyRescalerUpgrade::Electron)	`egammaAnalysisUtils`	Data AND MC	n.a
	*Energy resolution (smearing)*	2011-2012	EnergyRescalerUpgrade	getSmearingCorrection(Eta, E, egRescaler::EnergyRescalerUpgrade::NOMINAL	`egammaAnalysisUtils`	MC	SetRandomSeed(EventNumber+100*i)
							i index within el_ branches
	*Electron efficiencies (ID)*	2011	egammaSFclass	RunNumber=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber); scaleFactor(el_cl_eta Et, 5, 0, 6, 1,RunNumber)	`egammaAnalysisUtils`	MC	n.a
	*Electron efficiencies (reco)*	2011	egammaSFclass	RunNumber=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber); scaleFactor(el_cl_eta, Et, 4, 0, 6, 1,RunNumber)	`egammaAnalysisUtils`	MC	n.a
	*Electron efficiencies (ID)*	2012	egammaSFclass	RunNumber=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber); scaleFactor(el_cl_eta Et, 30, 0, 8, 1,RunNumber)	`egammaAnalysisUtils`	MC	n.a
	*Electron efficiencies (reco)*	2012	egammaSFclass	RunNumber=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber); scaleFactor(el_cl_eta, Et, 4, 0, 8,1, RunNumber)	`egammaAnalysisUtils`	MC	n.a
	*Calorimeter Isolation*	2011	CaloIsoCorrection	GetNPVCorrectedIsolation(unsigned int nPV, etaS2, radius, is_mc, Etcone_value, CaloIsoCorrection::ELECTRON);	`egammaAnalysisUtils`	Data and MC	n.a
	*Calorimeter Isolation*	2012	CaloIsoCorrection	GetPtEDCorrectedTopoIsolation(argument of this function);	`egammaAnalysisUtils`	Data and MC	n.a
	*Trigger-Offline matching*	2011-2012	TrigMuonEfficiency	electronTriggerMatchTool->match(el_tracketa, el_trackphi,chain);Et>thresold	`TrigMuonEfficiency`	Data and MC	n.a
	*Electron efficiencies (trigger)*	2011	TrigMuonEfficiency	leptonSF = new LeptonTriggerSF(2011, "TrigMuonEfficiency/share/", "muon_trigger_sf_mc11c.root"); runnumber_sf=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber);p=loosepp;q=loose; pair SF = SF = leptonSF->GetTriggerSF(runnumber_sf,false, muons, q, electrons, p, 0);SF.first	`TrigMuonEfficiency`	MC	n.a
	*Electron efficiencies (trigger)*	2012	TrigMuonEfficiency	leptonSF = new LeptonTriggerSF(2012, "TrigMuonEfficiency/share", "muon_trigger_sf_2012_AtoE.root "); /afs/cern.ch/user/t/takashi/public/SF2012/muon_trigger_sf_2012_AtoD.root runnumber_sf=PileUpTool->GetRandomRunNumber(physicsObject->RunNumber);p=ML;q=loose; pair SF = SF = leptonSF->GetTriggerSF(runnumber_sf,false, muons, q, electrons, p, 0);SF.first	`TrigMuonEfficiency`	MC	n.a
	*H4l Loose++ ID macro*	2011	isLoosePlusPlusH4l	isLoosePlusPlusH4l(el_eta_s2,Et_cl_s2_uc,rHad_uc, rHad1_uc,Reta, w2,f1,wstot,DEmaxs1,deltaEta,nSi,nSiOutliers,nPix,nPixOutliers,false,false)	`HiggsZZ4lUtils`	Data and MC	n.a
	*MultiLepton ID macro*	2012	passMultiLepton	passMultiLepton(argument of this function)	`egammaAnalysisUtils`	Data and MC	n.a
Pileup-reweighting	*Pileup-reweighting*	2011-2012	PileupReweighting	TPileupReweighting(see section : "Data Period Fraction and PileUpReweighting")	`PileupReweighting`	MC	n.a
Mu Rescaling	*mu rescaling*	not used	PileupReweighting	PileUpTool->SetDataScaleFactors(1/1.11);	`PileupReweighting`	MC	n.a
				That minbias analysis concludes to scale the MC by a factor of 1.11. Scaling MC can be done with the tool (using SetMCScaleFactors), but, because the mc mu distribution is discrete, unless you are careful about rebinning the mu distribution after rescaling, you will end up with holes in it. This is possible to do with the tool (AddBinning to redefine the binning for the mu distribution, probably using variable bin widths) but in the end it seemed more trouble than it was worth when instead of scaling the MC by 1.11, we could scale the Data by 1/1.11, which is going to be approximately equivalent. So rather than shifting the mc mu distribution up, we shift the data mu distribution down. We do this with: SetDataScaleFactors(1./1.11)
Good Run List	*Selection of good runs and lbs*	2011-2012	GoodRunList	HasRunLumiBlock	`GoodRunLists`	Data	n.a
Z vertex reweighting	*z vertex reweighting*	2012	VertexPositionReweightingTool	ZVertexTool->GetWeight(vxp_z);vxp_z=(physicsObject->mc_vx_z)[2];*	`egammaAnalyisUtils`	MC	n.a

ptms =(1/fabs((*TreeObject->mu_me_qoverp)[i])*sin((*TreeObject->mu_me_theta)[i]));
ptid =(1/fabs((*TreeObject->mu_id_qoverp)[i])*sin((*TreeObject->mu_id_theta)[i]));
int_lum, unit, directory : see MCP twiki for more details
Et (electron): el_cl_E/cosh(el_tracketa)
Et_cl (electron): el_cl_E/cosh(el_cl_eta)
Et_cl_s2_uc (electron): el_cl_E_uc/cosh(el_eta_s2)
el_eta_s2=el_etas2
rHad_uc=((el_Ethad))/Et_cl_s2_uc
rHad1_uc=((el_Ethad1))/Et_cl_s2_uc
Reta=(el_reta)
w2=(el_weta2)
f1=(el_f1)
wstot=(el_wstot)
DEmaxs1=( (el_emaxs1)-(*TreeObject->el_Emax2))/( (el_emaxs1)+(el_Emax2))
deltaEta=(el_deltaeta1)
nSi=(el_nSiHits)
nSiOutliers=(el_nSCTOutliers)+(el_nPixelOutliers)
nPix=(el_nPixHits)
nPixOutliers=(el_nPixelOutliers)
chain="EF_e20_medium", ...
mc_with_constant_term=1 for MC11a,b',b and 0 for MC11c
vector<TLorentzVector> electrons
vector<TLorentzVector> muons
electron_quality p=loosepp
muon_quality q=loose
vxp_z=mc_vx_z)[2]

Cut Flows

Samples not decided

Spreadsheets, divided in three channels (4e, 2e2mu and 4mu, where electrons are GSF and muons are STACO) are organized are follows:

MC12

MC11 c

Background Controls

Systematic Uncertainties

Major updates: -- BingLi - 09-Nov-2012

%REVIEW% Never reviewed

Attachments

Topic attachments
I	Attachment	History	Action	Size	Date	Who	Comment
cxx	BkgCrossSection.cxx	r1	manage	33.2 K	2012-11-10 - 00:01	BingLi	get from HiggsZZ4lUtils

Topic revision: r8 - 2012-11-12 - BingLi

Main

Webs

Public webs

Welcome Guest

- Cern Search
- TWiki Search
- Google Search
Main All webs

Copyright &© 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

Single resonance Z->4l analysis

Introduction

Meetings

Mailing list

Special notes

2011 versus 2012

Data and MC samples (lists of D3PD files)

ttbar

Z+Jets

Zbb

MC normalization (cross sections, k-factors, filter efficiencies)

GoodRunList

2012 Data, Luminosity per Period

Data Period Fraction and PileUpReweighting

Event Selection

Trigger details

2011Trigger matcing for D3PD tag p956 and p997

2012 Trigger matching

ID/MS Requirements

More information

Cut Flows

Background Controls

Systematic Uncertainties