r12 - 05 Dec 2010 - 10:48:56 - CarlosEscobarYou are here: TWiki >  Atlas Web  >  TopPhysics > TopQuarkMass

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


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) ARATopQuarkAnalysis + TopInputs lepton + jets txt
Data (Period A-F) ARATopQuarkAnalysis + TopInputs e + jets txt
Data (Period A-F) ARATopQuarkAnalysis + TopInputs mu + jets txt

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


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

Paper Date
pdf ATL-COM-PHYS-2009-111. Top Quark Mass Measurements 2009
pdf ATLAS CSC note (T9). Top quark mass measurement with ATLAS September 22, 2008
pdf 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
pdf 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
pdf 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
pdf ATL-PHYS-INT-2005-002. Top quark mass measurement in the lepton plus jets channel using full simulation 2005


SVN repository

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
cd Reconstruction/RecExample/RecExCommon/cmt
cmt config
source setup.sh

    • 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:

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"]

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'] += ['mc08.105200.T1_McAtNlo_Jimmy.merge.AOD.e357_s462_r635_t53/::mc08.105200.T1_McAtNlo_Jimmy.merge.AOD.e357_s462_r635_t53']

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


-- CarlosEscobar - 18 Oct 2010

toggleopenShow attachmentstogglecloseHide attachments
Topic attachments
I Attachment Action Size Date Who Comment
txttxt data_AF_TopPhys1561331_mu+jets.txt manage 2.5 K 18 Oct 2010 - 06:36 CarlosEscobar data_AF_TopPhys1561331_mu+jets.txt
txttxt data_AF_TopPhys1561331_e+jets.txt manage 2.4 K 18 Oct 2010 - 06:39 CarlosEscobar data_AF_TopPhys1561331_e+jets.txt
txttxt old_data_production.txt manage 3.2 K 18 Oct 2010 - 07:24 CarlosEscobar old_data_production.txt
txttxt mc_105200_TopPhys1561331_lepton+jets.txt manage 0.2 K 04 Nov 2010 - 10:21 CarlosEscobar mc_105200_TopPhys1561331_lepton+jets.txt
Edit | WYSIWYG | Attach | PDF | Raw View | Backlinks: Web, All Webs | History: r12 < r11 < r10 < r9 < r8 | More topic actions
Powered by TWiki
This site is powered by the TWiki collaboration platformCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback