ZhuTopAnalysis < Main

Main Web>TWikiUsers>ChengguangZhu>ZhuTopAnalysis (2011-08-14, ChengguangZhu) (raw view)
%TOC%

---++ relavant twiki

ZhuEtmiss 

ZhuTopPhysicsAimedPlotsAndTables

---++ W helicity measurement plan, aiming at 2010 and 2011 data

---+++ Man power and responsibilities 
   * Liming Chen: 
      * Lepton selection: implement the lepton selection in code, and update it.
      * KLFitter reco method,  
      * Systematics from PDF, top mass, 
      * Data skiming  
   * Bo Li: 
      * DD Wjets background estimation
      * Systematics from b-tagging efficiency
      * Electrion isEM efficiencies estimation from data
   * Lianyou Shan: 
      * Jets selection: implement the lepton selection in code, and update it.
      * Final parameters extraction 
      * Correction function for each helicity state
      * Systematics from pileup, ISR/FSR,  JES,
   * Chengguang Zhu: 
      * Etmiss selection: implement the lepton selection in code, and update it.
      * Etmiss performance: connection with Etmiss liaison gruop
      * Background shape estimations 
      * Systematics from shape of background met scale 
   * Wang Jin:
      * DD QCD estimation, 

---+++ time scale
   * conference ahead:
      * June 6:  Physics at LHC,
      * July 22:  EPS 2011 
      * Aug 22:  lepton-Photon 2011
      * Sep 25:  Top2011
   * Luminosity schedule:
      * March 18 -- March 25:  11pb-1
      * April  17 -- May 4:  100pb-1
      * May 14 -- June 15:  600pb-1
   * Note aproval:
      * INT note posting to TOPWG approval 3 days
      * TOPWG approval to CIRCULATION 10 days
         *Includes writing of CONF note, EB approval of conf note, PhysCoord signoff
      * ATLAS circulation takes 7 days
      * CIRCULATION to APPROVAL 7 days
         * Includes EB final signoff and final reading

   * Proposed EPS schedule – working backwards
      * Thu July 21- Conference starts
      * Mon July 18 – Final reader approval
      * Tue July 12 – ATLAS approval meeting
      * Mon July 4 – CONF note circulation
      * Tue June 28 – top-wg approval on full data
      * Mon June 27 – INT/CONF on full data posted
      * Wed June 22 – full dataset available for analysis
      * Wed June 15– End of LHC run
      * Tue June 14 – top-wg analysis pre-approval
      * Fri June 10 - INT note posted
      * Fri June 3 – Cut dataset for preapproval
      * Thu May 19– Signoff on TopPhys production cache
      * Sat May 14– Start of LHC run
      * Thu May 12 – Finalization of object definitions

   * Aiming at top2011, inverse time direction:
      * Sep 25: Conference starts
      * Sep 22: – Final reader approval
      * Sep 16: – ATLAS approval meeting
      * Sep 8 : – CONF note circulation
      * Sep 2: – top-wg approval on full data
      * Sep 1:  – INT/CONF on full data posted
      * Aug 25 : – full dataset available for analysis
      * Aug 17:  – top-wg analysis pre-approval
      * Aug 13:  – INT note posted
      * Aug 6:  – Cut dataset for preapproval


---+++ current task
   * include KLFitter into D3PD processing code
   * data driven method to estimate QCD events number and cospsi shape
   * data driven method to estimate Wjets events number and cospsi shape
   * Maximum Likelyhood method to draw the FL-FR 2D fit
   * data slimming and skimming
   * check all needed data is there, otherwise to produce.


---+++ analysis procedure
   * Objects selection: [[https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopCommonObjects][topinputs selection]], any methods to improve the selection to get better S/B
   * Reconstruction method: KLFitter and CSC method.
   * Background estimation: data driven for QCD/Wjets and MC
   * Bacground to be subtrcted from the all observables
   * Correction function from MC to correct the cospsi distribtuion. should verify the correction function for FL, FR and F0 part are same or not
   * Histogram fit or LH fit
   * Systematics Estimation

---+++ Full plan
---++++ 2010 data [[https://twiki.cern.ch/twiki/bin/view/AtlasProtected/AtlasProductionGroup][ref page]]
   * 2010 data: top group central produced D3PD 
      * real data: period E5 - I, 35.3 pb-1 (+/-3.2%)
      * MC: 
         * signal: [[http://arxiv.org/PS_cache/arxiv/pdf/1012/1012.1792v2.pdf][ref]] MC@NLO v3.41 [7], assumed top-quark mass of 172.5 GeV, NLO PDF set CTEQ66
         * backgroup: 
            * Alpgen v2.13 implementing the exact LO matrix elements and LO PDF set CTEQ6L1, &lsquo;MLM&rsquo; matching. 
               * W+jets up to 5 partons, W+c+jets, Wbb, Wcc
               * Z+jets up to 5 partons and dilepton invariant mass mll &gt; 40 GeV, Zbb, Zcc
               * QCD multi-jet up to 6 partons, bb, cc
               * diboson WW+jets, WZ+jets and ZZ+jets events.
            * Pythia: Z boson covering 10 GeV &lt; m\x{2113}\x{2113} &lt; 40 GeV.
            * MC@NLO: single-top, removing overlaps with tt events.

---++++ 2011 data
   * Collision data:
      * Collected in period B and D, please use the GRL: [[%ATTACHURL%/Top_GRL_D.xml][Top_GRL_D.xml]] 
      * The list of the datasets is  [[%ATTACHURL%/datasets_2011_realdata][datasets_2011_realdata]]
         * These data is slimed to 2 datasets, if it satisfied  [[%ATTACHURL%/SlimCuts_leptonMoreThan0_jetsMoreThan3.txt][SlimCuts_leptonMoreThan0_jetsMoreThan3]]:
            * user.zhuchengguang.Slim_physics_Muons_201106080525/,  has transfered to CERN-PROD_SCRATCHDISK,
            * user.zhuchengguang.Slim_physics_Egamma_201106080525/, except one sub dataset, has been transfered to CERN-PROD_SCRATCHDISK

   * MC data. 
      * MC10a: the list of the MC10a data used is [[%ATTACHURL%/datasets_2011_MC10a][datasets_2011_MC10a]].
         * these data has been slimed to  user.zhuchengguang.Slim_MC_201106080525/, most of them copied to  CERN-PROD_SCRATCHDISK, some of them still in the site where it was slimed.
      * MC10b: we can now switch to these data: [[%ATTACHURL%/datasets_2011_MC10b][datasets_2011_MC10b]]
         * MC10 vs MC10a vs MC10b: use MC10 for 2010 data analysis, use MC10a,b for 2011 data analysis.
         * it will be slimed soon. 
         * If you want to use more data than you can taken from [[https://twiki.cern.ch/twiki/pub/AtlasProtected/TopMC10For2011Data/MC10b_D3PD_status_05Jun.txt][this list]]
         * 50ns bunch spacing in 3 trains of 36 bunches with 225ns train separation

       * MC11
            *  planned for release 17, used for beyond summer conference.
 
---++++ Background estimation
   * ä»Ždataé€‰æ‹©ç‰¹å®šäº‹ä¾‹çš„æ–¹æ³•
      * TTbaré€‰æ‹©æ–¹æ¡ˆä¸ºtopç»„æ ‡å‡†æ–¹æ¡ˆ
      * Wjetsé€‰æ‹©æ–¹æ¡ˆå’ŒTTbaré€‰æ‹©æ–¹æ¡ˆé™¤äº†Bjetsé€‰æ‹©å¤–ï¼Œå…¶ä»–ç›¸åŒã€‚Wjetsé€‰æ‹©æ–¹æ¡ˆè¦æ±‚æ²¡æœ‰Bjetsã€‚
         * å¯¹Wjetsäº‹ä¾‹çš„é€‰æ‹©æ•ˆçŽ‡ï¼šå‡è®¾light jetçš„fake rateä¸ºfï¼Œä¸€èˆ¬fæ˜¯ä¸ªå¾ˆå¤§çš„æ•°ï¼Œå…¶ä»–cuté€‰æ‹©äº†4ä¸ªjetæ—¶ï¼Œå†è¦æ±‚å…¶ä¸æ²¡æœ‰bjetï¼Œé€šè¿‡çŽ‡ä¸º~1-4*1/fï¼Œ 5ä¸ªjetæ—¶çš„é€šè¿‡çŽ‡ï¼š1-5*1/fï¼Œè¿œå¤§äºŽTTBaré€‰æ‹©æ–¹æ¡ˆå¯¹Wjetsäº‹ä¾‹çš„é€‰æ‹©ã€‚
         * TTbaräº‹ä¾‹çš„æ··å…¥ï¼šæ ¹æ®Bjet tagæ•ˆçŽ‡50%ï¼Œå¯çŸ¥25%çš„äº‹ä¾‹ä¸å°†æ²¡æœ‰Bjetsï¼Œå› æ¤å¯ä»¥æ··å…¥25%çš„TTbaräº‹ä¾‹ï¼ˆå¿½ç•¥è½»jetçš„fakeï¼‰ï¼Œå’ŒTTbaré€‰æ‹©æ–¹æ¡ˆé€‰æ‹©TTbaräº‹ä¾‹çš„æ•ˆçŽ‡ç›¸å½“ã€‚
         * ä¸ºäº†æ¯”è¾ƒwjetsäº‹ä¾‹çš„é‡å»ºæ•ˆæžœï¼Œä»Žlight jetä¸ä»¥50%çš„æ¦‚çŽ‡ï¼Œä»»æ„æŒ‡å®šä¸ºbjetsï¼ŒéšåŽå†é€šè¿‡top pairé‡å»ºã€‚
      * QCDé€‰æ‹©æ–¹æ¡ˆï¼šè¦æ±‚äº‹ä¾‹ä¸æ²¡æœ‰è½»åï¼ˆä»»ä½•å¯ç–‘è½»åï¼‰ï¼Œæ²¡æœ‰b tagçš„jetsï¼Œå…¶ä»–å’ŒTTbaré€‰æ‹©æ–¹æ¡ˆç›¸åŒã€‚
         * é€šè¿‡é€‰æ‹©çš„äº‹ä¾‹å¾ˆå°‘ï¼Œæ˜¯å¦åº”è¯¥æ”¾æ¾ä¸€ä¸‹ï¼ˆæ²¡æœ‰å¯ç–‘è½»åï¼‰çš„é€‰æ‹©ã€‚
         * æ²¡æœ‰b tagçš„jetçš„é€‰æ‹©æ˜¯å¦å¿…è¦ï¼Ÿ
         * ä¸ºäº†æ¯”è¾ƒQCDäº‹ä¾‹çš„é‡å»ºæ•ˆæžœï¼Œä»Žjetä¸ä»»æ„æŒ‡å®šä¸€ä¸ªjetä¸ºè½»åï¼Œä»Žå…¶ä½™çš„light jetä¸ä»¥50%çš„æ¦‚çŽ‡ï¼Œä»»æ„æŒ‡å®šä¸ºbjetsï¼ŒéšåŽå†é€šè¿‡top pairé‡å»ºã€‚
   * Background estiomation 
      * Data driven estimate QCD background. 
         * method 1 (obsolete): Take control sample: small met & leading pt electron or muon without quality cut as the lepton (remove jets which overlap with the lepton). these should be dominated by QCD events. Pass them through ttbar reconstruction to determine the cospsi shape of QCD. Scale the shape to the estimated number of QCD events. This should be validated whether the shape of cospsi of different met bin are same: it is proved NOT, because the met is used in the W reconstruction and cospsi calculation. 
            * method1 patch (obsolete): on the "top" hadronic side, omitting Pz of one light jet (randomly) to mimic met, the other light jet as lepton. Then reconstruct the W as if W decay leptonically and then "leptonic top". If the cospsi distributution of both "lep and had" top side are same, then we can use low met to take pure QCD events and use the 3 jets to mimic lepton and met to calculate the cospsi. the met pt distribution and the jet pt distribution are not same, which makes the cospsi shape NOT same.
         * <span style="color: red;" class="WYSIWYG_COLOR"><span style="color: blue;" class="WYSIWYG_COLOR"> </span><span style="color: red;" class="WYSIWYG_COLOR"><span style="color: blue;" class="WYSIWYG_COLOR"> </span><strong><span style="color: lime;" class="WYSIWYG_COLOR"><span style="color: blue;" class="WYSIWYG_COLOR"><span style="color: blue;" class="WYSIWYG_COLOR"><span style="color: red;" class="WYSIWYG_COLOR">method 2:</span></span></span></span></strong> </span></span> 
            * MC or control data: cuts (e.g. no medium electron above ) select un-biased dominant QCD events, other backround and signal should be supressed
            * electron in QCD comes mainly from jet-fake, so take random jet as electron. Cut "electron" and met &gt;20GeV.
            * Should prove that electron + met =W and jet+met=W are same in 4m?
            * take one bjet or random jet as bjet (if there is). Shoud test here W+bjet=top has same cospsi shape for jet-electron and electron.
            * Should test extra jets change the cospsi distribution or not?, if yes, need 3 extra jets in the event selection.
            * need the top mass is around 172.5
            * calculate the cospsi, as the background shape of QCD.
            * muon channel, muon in QCD comes mainly from jet decay, hard muon takes most energy of the closest jet? muon is little softer than electron and jets?
            * check 4m of muon, take jets as muon, scale event weight if 4m not same
            * if no much statistical, use looser cuts to test this distribution is same as the standard reco method.
      * Data driven estimate WJets background. 
         * <strong><span style="color: blue;" class="WYSIWYG_COLOR"><span style="color: red;" class="WYSIWYG_COLOR">method 1:</span></span></strong> 
            * MC or data: top group [[https://twiki.cern.ch/twiki/bin/view/AtlasProtected/Winter2011CutFlow][standard cut for version 16]] .
            * ä½¿ç”¨W/Z ratioæ–¹æ³•ä¼°è®¡Wjetsäº‹ä¾‹çš„æ•°é‡
            * ä½¿ç”¨Wjetsé€‰æ‹©æ–¹æ¡ˆé€‰æ‹©Wjetsäº‹ä¾‹ï¼Œä¼°è®¡äº‹ä¾‹çš„cospsiå½¢çŠ¶
               * æ¯”è¾ƒWjetsé€‰æ‹©æ–¹æ¡ˆå’ŒTTbaré€‰æ‹©æ–¹æ¡ˆé€‰æ‹©çš„Wjetsäº‹ä¾‹çš„æ‰€æœ‰å˜é‡å…·æœ‰ç›¸å½“çš„è¿åŠ¨å¦ç‰¹å¾
            * 
               * æ¯”è¾ƒWjetsé€‰æ‹©æ–¹æ¡ˆå’ŒTTbaré€‰æ‹©æ–¹æ¡ˆé€‰æ‹©çš„Wjetsäº‹ä¾‹çš„cospsi shape is same.
      * MC estimate of other non ttbar background
      * MC estimate conbinatorial background
   * Extract measurement 
      * Reconstruct ttbar events with KLFitter and CSC method when statistics increase(?)
      * Remove background of above estimate
      * Correction function take from signal MC method to correct the signal distortion
      * Function fitting (LH fit for the 2D FL-FR)

---++++ Systematics estimation
   * docu: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopSystematicUncertainties
   * relative x-sections of background to signals 
      * in the fit, changing the relative background ratio to evaluate the errors.
   * background shape uncertainties, with psudo-experiment to estimate the errors 
      * data-driven estimated Wjets shape:
      * data-driven estimated QCD shape:
      * other small background
   * ISR/FSR, PDF, top mass, JES, The b-tagging efficiency and mistag fraction, met scale makes the correction function different, which makes the fit results different, as the systematic error 
      * ISR/FRS: more samples varying the parameters controlling ISR and FSR in a range consistent with experimental data. Conventionally using ISR/FSR on/off samples.
      * top mass: difference top mass sample.
      * PDF: PDF reweight method
      * JES: varies in the range 6&ndash;10% as a function of jet pT and Eta;, the varies are used to reproduce the correction function. 
         * 2011.02.17: a recommended preliminary jet energy scale uncertainty for b-tagged jets (using SV0 tagger) of 2.5% for pt&gt;25 GeV.
      * The b-tagging efficiency and mistag fraction: the varies are used to reproduce the correction function.
      * met scale: he varies are used to reproduce the correction function.
   * pileup

---++++ Ensemble test
   * compare the final objects' kinematics of MC ttbar events with data (ttbar control sample) 
      * to make sure whether the combinetorial bagd is same for data and MC

---++++ Code
   * Version_1_0: cut flow OK for version 16 
      * svn co svn+ssh://svn.cern.ch/reps/topana/topD3PD/tags/version_1_0 cutflow
   * Version_2_0: add Wjets and QCD selection procedure for version 16 
      * svn co svn+ssh://svn.cern.ch/reps/topana/topD3PD/tags/version_2_0 topAna
   * develeping version ( which is updated regularly 
      * svn co svn+ssh://svn.cern.ch/reps/topana/topD3PD/trunk MyD3PD_Ana
   * 2011-08-12, Whel note proval meeting, corresponding to our int note   [[%ATTACHURL%/main.pdf][version 4]] at http://cdsweb.cern.ch/record/1361020. The code to get this note is 
      * D3PD process code:  svn co svn+ssh://svn.cern.ch/reps/topana/topD3PD/tags/version_5_0
      * Analysis code:  svn co -r 233 svn+ssh://svn.cern.ch/reps/topana/QCDestimate/trunk 
      * The note tex: svn co -r 232 svn+ssh://svn.cern.ch/reps/topana/WhelIntNote/trunk

---++ Test run
---+++ run small example
   * data file\x{ff1a}user.jsearcy.group.phys-top.D2PD_data10_7TeV.periodI.physics_Muons.NTUP_TOP.TOPMU.repro05_v02_pp_16_0_3_3.002949.EXT0._00495.NTUP_TOP.root

   * run with parameter MC\x{ff1a} 
      * total events\x{ff1a}4995: 1874 el triggered, 1694 mu triggered, 1760 not triggered , 333 both triggered.
      * selected 334 = e 143 + mu 191\x{ff0c} 32 with e+mu trigger

   * run as muon stream: 
      * 1694 mu triggered 3300 not mu triggered
      * selected 191, 178 Topl recoed.

   * run as electron stream: 
      * 1874 e triggered 3020 not mu triggered
      * 143 e triggerd, 89 recoed


---++ needed material:

---+++ reference

https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopReferences

---+++ trigger

trigger search: https://atlas-trigconf.cern.ch/

atlas trigger menu: https://twiki.cern.ch/twiki/bin/viewauth/Atlas/TriggerPhysicsMenu

---+++ GRL

ATLAS Official GRL: http://atlasdqm.web.cern.ch/atlasdqm/grlgen/

GRL: https://atlas-top-grls.web.cern.ch/atlas-top-grls/ and how to use: http://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=93780

object selection: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopCommonObjects

---+++ Data

Data period: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/DataPeriods

top group D3PD variable list https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopInputsD3PDMakerBranches

top crosssection area: https://espace.cern.ch/topxsec/default.aspx

top group dAOD https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopD2PD2010

Top MC and data https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopD3PDProduction

https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopMC2009

Topinput D3PD maker https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopInputsD3PDMaker

lowest un-prescaled trigger chain: https://twiki.cern.ch/twiki/bin/viewauth/Atlas/LowestUnprescaled

remove double event: https://espace.cern.ch/topxsec/twiki/Double%20Events.aspx

---++ Analysis

e_mu overlap removal
   * https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopInputsD3PDMaker#Why_does_el_overlap_mu_not_cut_a

lumi  calculation
   * https://atlas-datasummary.cern.ch/lumicalc/

cut flow:
   * https://twiki.cern.ch/twiki/bin/view/AtlasProtected/Winter2011CutFlow
   * https://twiki.cern.ch/twiki/bin/view/AtlasProtected/Winter2011TopCouplings

KLFitter
   * https://svnweb.cern.ch/trac/atlasgrp/browser/Institutes/Goettingen

pileup reweighting
   * https://twiki.cern.ch/twiki/bin/view/AtlasProtected/PileupReweighting

btag&#65306;
   * https://twiki.cern.ch/twiki/bin/view/AtlasProtected/BTaggingForEarlyData#SV0_Tagger


---+++ regulations:

plot style: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopPlot

---++ relative result

top cross section plot: https://twiki.cern.ch/twiki/bin/view/AtlasProtected/TopCrossSectionPreliminary

---++ Skim and Slim

|  No.  | select condition  |  Data type  |  Data name  |  !DaTRI Site  |
|  1.1  |  topEl_n+ topMu_n &gt; 0  |  data10 Egamma(E4-I)  |  user.limingchen.slimD3PD_electron_201102202133/  |  CERN-PROD_SCRATCHDISK  |
|  1.2  |  topEl_n+ topMu_n &gt; 0  | data10 Muons(E4-I)  | user.limingchen.slimD3PD_muon_201102202133/  |  CERN-PROD_SCRATCHDISK  |
|  1.3  |  topEl_n+ topMu_n &gt; 0  |  MC  |  user.limingchen.slimD3PD_MC_201102202133/  |  CERN-PROD_SCRATCHDISK  |

DataIsHere
Attachments
Topic attachments
I	Attachment	History	Action	Size	Date	Who
txt	SlimCuts_leptonMoreThan0_jetsMoreThan3.txt	r1	manage	2.6 K	2011-06-11 - 18:39	ChengguangZhu
xml	Top_GRL_D.xml	r1	manage	34.8 K	2011-06-11 - 18:29	ChengguangZhu
ext	datasets_2011_MC10a	r1	manage	7.2 K	2011-06-11 - 18:31	ChengguangZhu
ext	datasets_2011_MC10b	r1	manage	7.1 K	2011-06-11 - 18:30	ChengguangZhu
ext	datasets_2011_realdata	r1	manage	8.5 K	2011-06-11 - 18:30	ChengguangZhu
pdf	main.pdf	r1	manage	544.0 K	2011-08-14 - 00:32	ChengguangZhu
Topic revision: r46 - 2011-08-14 - ChengguangZhu
Main
Webs
Public webs
Welcome Guest
- Cern Search
- TWiki Search
- Google Search
Main All webs
Copyright &© 2008-2024 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