Top Quark Mass Measurement

Top quark analysis within Athena

Here it is how to run ARATopQuarkAnalysis: Release Recipes to run ARATopQuarkAnalysis

Kinematic fitters comparison framework

One can find the guidelines to perform the first kinematic fitters comparison: Giorgio's minutes

Samples

To know more about the mc09 production take a look the to TopMC2009 twiki

D3PD Production

LHC @ 7 TeV (D3PD)
Sample	TopPhys rel.	Framework	Channel	Location
MC (105200)	15.6.13.3.1	ARATopQuarkAnalysis + TopInputs	lepton + jets
Data (Period A-F)	15.6.13.3.1	ARATopQuarkAnalysis + TopInputs	e + jets
Data (Period A-F)	15.6.13.3.1	ARATopQuarkAnalysis + TopInputs	mu + jets

Note the new jet selection requirement: pT(jets) > 25 GeV, >=1 b-jet,... Take a look to this talk.

Documentation

Here one will find all the related notes to the top quark mass reconstruction:

Paper	Date
ATL-COM-PHYS-2009-111. Top Quark Mass Measurements	2009
ATLAS CSC note (T9). Top quark mass measurement with ATLAS	September 22, 2008
ATL-COM-PHYS-2008-117. Kinematic Fit of ttbar Events in the Semi-Leptonic Decay Channel with a Chi2 Method at the ATLAS Experiment	August 22, 2008
ATL-PHYS-INT-2008-021. Top quark mass measurement in the lepton plus jets channel with a kinematic fit on the full final state	April 21, 2008
ATL-COM-PHYS-2009-028. Measurement of the top quark mass using the template method in the lepton plus jets channel with ATLAS	January 19, 2009
ATL-PHYS-INT-2005-002. Top quark mass measurement in the lepton plus jets channel using full simulation	2005

Software

SVN repository

• The ROOT based analysis package can be found here: TopMassAnalysis package

ARA D3PD production

• • follow instructions of ARATopQuarkAnalysis for release 14.5.2. These are basically the following:

source ~/cmthome/setup.sh -tag=AtlasOffline,14.5.2,setup,32,opt
- Check out the following packages:
cmt co -r AthenaROOTAccess-00-05-43  PhysicsAnalysis/AthenaROOTAccess
cmt co -r Navigation-00-06-09-02  Control/Navigation
cmt co -r MissingET-03-01-25-03 Reconstruction/MissingET
cmt co -r RecExCommon-00-10-55 Reconstruction/RecExample/RecExCommon
cd PhysicsAnalysis
cvs co -d ARATopQuarkAnalysis -r ARATopQuarkAnalysis-00-00-76 groups/LPC/ARATopQuarkAnalysis
- Copy to ARATopQuarkAnalysis the following files:
cp /afs/cern.ch/user/m/mjcosta/public/TOP/VERSION2/ARAParticleJet.cxx  ARATopQuarkAnalysis/src/
cp /afs/cern.ch/user/m/mjcosta/public/TOP/VERSION2/ARAParticleJet.h ARATopQuarkAnalysis/ARATopQuarkAnalysis/
- Compile all packages
cd  ARATopQuarkAnalysis/cmt
cmt br cmt config
source setup.sh
cmt br gmake
cd PhysicsAnalysis/AthenaROOTAccess/cmt
cmt br cmt config
source setup.sh
cmt br gmake
cd Reconstruction/MissingET/cmt
cmt config
source setup.sh
gmake
cd Reconstruction/RecExample/RecExCommon/cmt
cmt config
source setup.sh
gmake

• • Copy the following file in ARATopQuarkAnalysis^?/share/: /afs/cern.ch/user/m/mjcosta/public/TOP/VERSION2/MyTopQuarkAnalysis_AODtoD3PD_jobOptions.py which is re-running the missing Et and the b-tagging due to some bugs in the release used to create the AODs, a different selection for the R_* objects stored in D3PD and the following additional lines to store TruthJets and TruthMET:

========================================================
TO ADD TRUTH JETS FOR CONE 0.4
========================================================
thisAlg = ParticleJetHandler("Truth_ParticleJetHandler")
thisAlg.InputParticleJetContainerName    = "Cone4TruthJets"
thisAlg.CutJetPt        =  0*GeV
thisAlg.CutJetEta       =  100.0
thisAlg.OutputParticleJetContainerName   = "Truth_Jets"
topSequence +=[thisAlg]
if m_doD3PDTTree:
    m_D3PDTTreeItems += ["ARAParticleJet#Truth_Jets#truth_jet"]

========================================================
TO ADD TRUTH MISSING ET
=========================================================
m_D3PDTTreeItems += ["ARAMissingET#MET_Truth#truth_missingET"]
 

• • The best is to run over the merged AODs produced in the latest reprocessing (called mc08*merge.AOD.*.r635_t53).
  • Set the corresponding AOD input in TopQuarkAnalysis_AODtoD3PD_jobOptions.py and run the jobOptions: athena TopQuarkAnalysis_AODtoD3PD_jobOptions.py
  • In order to run using Ganga:
    • Copy to ARATopQuarkAnalysis^?/cmt the file: /afs/cern.ch/user/m/mjcosta/public/TOP/VERSION2/ganga.mc08.10TeV.py
    • Choose the location of the jobOptions you want to run:

JobOptionFile = '$HOME/scratch0/MyAthena/14.5.2/PhysicsAnalysis/ARATopQuarkAnalysis/share/MyTopQuarkAnalysis_AODtoD3PD_jobOptions.py'

• • • Choose the dataset you want to process by changing the lines:

Datasets['TTBAR-MCATNLO::IFIC-LCG2_MCDISK'] = []
Datasets['TTBAR-MCATNLO::IFIC-LCG2_MCDISK'] += ['mc08.105200.T1_McAtNlo_Jimmy.merge.AOD.e357_s462_r635_t53/::mc08.105200.T1_McAtNlo_Jimmy.merge.AOD.e357_s462_r635_t53']

ProcessTAG = 'TTBAR-MCATNLO'
Sites = ['IFIC-LCG2_MCDISK']
j.outputdata.datasetname= 'mc08.105200.T1_McAtNlo_Jimmy.merge.DPD.e357_s462_r635_t53.v2'

This last line has to be set by hand because the previous line (j.outputdata.datasetname= outputDS.replace('AOD','DPD') + DatasetTAG) does not work if the name of the dataset is too long.

• • • source /afs/cern.ch/atlas/offline/external/GRID/ddm/DQ2Clients/setup.sh
    • voms-proxy-init -voms atlas
    • source /afs/cern.ch/sw/ganga/install/etc/setup-atlas.sh
    • ganga
    • execfile('ganga.mc08.10TeV.py')

Meetings

• TopPhysics Valencia/CERN Meetings can be found Meetings

-- CarlosEscobar - 18 Oct 2010