210 likes | 1.31k Views
NIPHAD meeting, Februari 9 th :. Fully Hadronic Top Anti-Top Decay (Using TopView). Ido Mussche. Outline. 1___ Introduction to Research Subject 2___ Analysis: TopView 3___ Some results 4 ___ Outlook. t t 6 jets. Introduction. Estimated cross section ttbar: 800 pb
E N D
NIPHAD meeting, Februari 9th : Fully Hadronic Top Anti-Top Decay (Using TopView) Ido Mussche
Outline 1___ Introduction to Research Subject 2___ Analysis: TopView 3___ Some results 4 ___ Outlook
t t 6 jets Introduction Estimated cross section ttbar: 800 pb Luminosity: 1fb-1 tt 6jets events: 35000 (in first year of running)
t t 6 jets Introduction • Research subject: • Optimize fully hadronic signal • Investigate and map background • Compare with semi-leptonic channel • Offline in cooperation with online trigger (Menelaos)
RDO Algorithm Configuration ESD Add your own stuff AOD TopView Analysis Results AANtuple Histos Or Small Ntuple AANtuple Analyze with Ntuple data Histos Histos Analysis
Analysis TopView: Collection of common tools for top group built around EventView framework Starting point for implementing analysis or use as NTuple dumper tool Includes C++ tools and Python modules C++ : Variable calculator (TruthLabeller, Total Invariant mass) Python: Inserters (Reco, Truth), Analysis (W reconstruction, Top reconstruction), Matching (Truth to Reco and Reco to Truth) https://twiki.cern.ch/twiki/bin/view/Atlas/TopView
Analysis • Levels of customization in TopView: • Use default top JobOptions file (change/include ntuple variables, include analyses) • Make your own top jobOption using existing python modules and C++ tools, creating your own analysis • Write your own python modules • Write your own C++ tools
AOD Electron Muon ParticleJet (b-tagged) Truth Analysis Step 1: Insertion and kinematics calculation TopView Insertion and overlap removal Insert objects in order, remove particles in same region Final State Objects Electron ParticleJet (b-tagged) Muon Truth Inferred Objects Kinematics calculation Loop over Final State Objects, calculate kinematics, add prefix and dump in User Data User Data El_p_T Mu_E PJet_eta EtMiss
Analysis Step 2: Combine jets to W and calculate variables TopView HadronicWReconstruction_module.py: EVSimpleCombo/HadronicWMaker( Labels[“ParticleJet”] OutputLabel[“HadronicW”] NDaughters[2] LowMass[60] HighMass[100] ) Final State Objects Electron ParticleJet (b-tagged) Muon Truth Inferred Objects Hadronic W User Data El_p_T Mu_E PJet_eta EtMiss HadW_p_T HadW_m
Analysis Step 3: Find closest b-jet to W and combine to top Loop over Ws and find+label closest B-Jet TopView TopReconstruction_module.py: EVGeneralCombo/TopMaker( Labels[“W”, “ChosenBJet”] OutputLabel[“TopLabel”] NDaughters[2] LowMass[0] HighMass[10000] ) Final State Objects Electron ParticleJet (b-tagged) Muon Truth Inferred Objects Hadronic W Top User Data El_p_T Mu_E PJet_eta EtMiss HadW_p_T HadW_m
Analysis Step 4: Calculate top variables, match with truth, and create AANTuple TopView Final State Objects Electron ParticleJet (b-tagged) Muon Truth Inferred Objects Hadronic W Top User Data El_p_T Mu_E HadTop_N PJet_eta EtMiss HadTop_m HadW_p_T HadW_m HadTop_Tru_m
Results All results are extracted from: Csc11.005204.TTbar_FullHad_McAtNlo_Jimmy.recon.AOD.v11004205 Reconstructed top mass fitted with simple Gaussian for 10000 events Mean should be 175 GeV
3269 161 GeV Results Reconstructed mass split up in matched (fitted with Gauss) and unmatched (fitted with Landau) From 20000 tops only 3269 matching are found.
Results Total reconstructed fitted with convoluted Gauss + Landau
Results Reco –True (eta) True (eta)
Results Reco –True (pT) True (pT)
Extra plots Reco –True (pT) True (eta)
Extra plots Reco –True (eta) True (pT)
Future goals • Optimize hadronic top reconstruction: • Map background • Check jet energy scale calibration • Compare with semi-leptonic channel • Kinematic fit