Documentation of the tracking D3PDs for TRT

Introduction

This wiki describes the tracking and vertex ntuples (TrackD3PDs) aimed for the TRT. The content of the D3PD is divided in two parts: default content and additional content. The default content is in all tracking D3PD which are produced by the tracking performance group and in any private production with default settings of the D3PD dumpers (i.e. you do not change any of the TrackD3PDMakerFlags). Additional content is modifying for TRT and contain drift circles information.

Tracking D3PD Variables

VariableName Format Description
ei_RunNumber UInt_t Run number
ei_EventNumber UInt_t Event number
ei_timestamp UInt_t Timestamp, in seconds
ei_timestamp_ns UInt_t Timestamp offset from last whole second, in ns
ei_lbn UInt_t Luminosity block number
ei_bcid UInt_t Bunch crossing ID
ei_detmask0 UInt_t Detector mask 0
ei_detmask1 UInt_t Detector mask 1
ei_actualIntPerXing Float_t Interactions per crossing for the current BCID - for in-time pile-up
ei_averageIntPerXing Float_t Interactions per crossing for all BCIDs - for out-of-time pile-up
trk_n Int_t Number of ntuple rows
trk_pt vector < float >  
trk_eta vector < float >  
trk_d0_wrtPV vector < float >  
trk_z0_wrtPV vector < float >  
trk_phi_wrtPV vector < float >  
trk_theta_wrtPV vector < float >  
trk_qoverp_wrtPV vector < float >  
trk_cov_d0_wrtPV vector < float >  
trk_cov_z0_wrtPV vector < float >  
trk_cov_phi_wrtPV vector < float >  
trk_cov_theta_wrtPV vector < float >  
trk_cov_qoverp_wrtPV vector < float >  
trk_chi2 vector < float >  
trk_ndof vector < int >  
trk_nBLHits vector < int >  
trk_nPixHits vector < int >  
trk_nSCTHits vector < int >  
trk_nTRTHits vector < int >  
trk_nTRTHighTHits vector < int >  
trk_nPixHoles vector < int >  
trk_nSCTHoles vector < int >  
trk_nTRTHoles vector < int >  
trk_nPixelDeadSensors vector < int >  
trk_nSCTDeadSensors vector < int >  
trk_nBLSharedHits vector < int >  
trk_nPixSharedHits vector < int >  
trk_nSCTSharedHits vector < int >  
trk_nBLayerSplitHits vector < int >  
trk_nPixSplitHits vector < int >  
trk_expectBLayerHit vector < int >  
trk_nMDTHits vector < int >  
trk_nCSCEtaHits vector < int >  
trk_nCSCPhiHits vector < int >  
trk_nRPCEtaHits vector < int >  
trk_nRPCPhiHits vector < int >  
trk_nTGCEtaHits vector < int >  
trk_nTGCPhiHits vector < int >  
trk_nHits vector < int >  
trk_nHoles vector < int >  
trk_hitPattern vector < int >  
trk_TRTHighTHitsRatio vector < float >  
trk_TRTHighTOutliersRatio vector < float >  
trk_fitter vector < int >  
trk_patternReco1 vector < int >  
trk_patternReco2 vector < int >  
trk_trackProperties vector < int >  
trk_particleHypothesis vector < int >  
trk_blayerPrediction_x vector < vector < float > > predicted x global position in the B-layer
trk_blayerPrediction_y vector < vector < float > > predicted y global position in the B-layer
trk_blayerPrediction_z vector < vector < float > > predicted z global position in the B-layer
trk_blayerPrediction_locX vector < vector < float > > predicted x local position in the B-layer
trk_blayerPrediction_locY vector < vector < float > > predicted y local position in the B-layer
trk_blayerPrediction_err_locX vector < vector < float > > error on local x prediction in the B-layer
trk_blayerPrediction_err_locY vector < vector < float > > error on local y prediction in the B-layer
trk_blayerPrediction_etaDistToEdge vector < vector < float > > distance of B-layer prediction to module edge along eta
trk_blayerPrediction_phiDistToEdge vector < vector < float > > distance of B-layer prediction to module edge along phi
trk_blayerPrediction_detElementId vector < vector < unsigned int > > detector element identifier of B-layer prediction
trk_blayerPrediction_row vector < vector < int > >  
trk_blayerPrediction_col vector < vector < int > >  
trk_blayerPrediction_type vector < vector < int > >  
trk_TRT_hit_n vector < int > Number of ntuple rows
trk_TRT_hit_id vector < vector < unsigned int > >  
trk_TRT_hit_detElementId vector < vector < unsigned int > >  
trk_TRT_hit_bec vector < vector < int > >  
trk_TRT_hit_layer vector < vector < int > >  
trk_TRT_hit_isCompetingRIO vector < vector < int > >  
trk_TRT_hit_locR vector < vector < float > >  
trk_TRT_hit_err_locR vector < vector < float > >  
trk_TRT_hit_trkLocR vector < vector < float > >  
trk_TRT_hit_trkLocPhi vector < vector < float > >  
trk_TRT_hit_err_trkLocR vector < vector < float > >  
trk_TRT_hit_err_trkLocPhi vector < vector < float > >  
trk_TRT_hit_cov_trkLocRPhi vector < vector < float > >  
trk_TRT_hit_TOT vector < vector < float > >  
trk_TRT_hit_highLevel vector < vector < int > >  
trk_TRT_hit_trailingEdge vector < vector < float > >  
trk_TRT_hit_chi2 vector < vector < float > >  
trk_TRT_hit_ndof vector < vector < int > >  
vx_n Int_t Number of ntuple rows
vx_x vector < float >  
vx_y vector < float >  
vx_z vector < float >  
vx_cov_x vector < float >  
vx_cov_y vector < float >  
vx_cov_z vector < float >  
vx_cov_xy vector < float >  
vx_cov_xz vector < float >  
vx_cov_yz vector < float >  
vx_type vector < int >  
vx_chi2 vector < float >  
vx_ndof vector < int >  
vx_px vector < float >  
vx_py vector < float >  
vx_pz vector < float >  
vx_E vector < float >  
vx_m vector < float >  
vx_nTracks vector < int >  
vx_sumPt vector < float >  
vx_trk_weight vector < vector < float > >  
vx_trk_n vector < int > Number of ntuple rows
vx_trk_index vector < vector < int > > Index in trk

Tracking D3PD Making

Selecting dataset

If there is no existing tracking D3PD, you can create it from ESD or AOD files. As an example, we use 7TeV 2011 data, 8TeV 2012 data or MC data-sets:
data11_7TeV.period*.physics_MinBias.PhysCont.AOD.repro09_v01/ (=D-K)
data11_7TeV.period
.physics_MinBias.PhysCont.AOD.t0pro09_v01/ (=L- )
data11_7TeV.
.physics_MinBias.merge.AOD.r2603_p659/ (B-H)
data11_7TeV.*.physics_MinBias.merge.AOD.r2713_p705/ (I-K)
data12_8TeV.*.physics_MinBias.recon.ESD.f447/ (B10)
data12_8TeV.*.physics_MinBias.merge.AOD.f447_m1153/ (B10)
mc11_7TeV.105001.pythia_minbias.merge.AOD.e815_s1272_s1274_r2727_r2700/
See also COMA Period Documentation Menu about data periods.

Download sub-sample of the data-set
You need to download a sub-sample of the dataset. Here. we assume to download into /data directory.

$ cd /data
$ dq2-get -n 1 data11_7TeV.periodD1.physics_MinBias.PhysCont.AOD.t0pro08_v01/
$ dq2-get -n 1 data11_7TeV.00179725.physics_MinBias.recon.ESD.r2603_tid495504_00
$ dq2-get -n 1 data12_8TeV.00204026.physics_MinBias.merge.AOD.f447_m1153
$ dq2-get -n 1 data12_8TeV.00204240.physics_MinBias.recon.ESD.f447
$ dq2-get -n 1 mc10_7TeV.105001.pythia_minbias.merge.AOD.e574_s932_s946_r2295_r2300/

Configuring jobOptions

Here is listing of TRT-specific flags in the jobOptions for producing D3PDs (which are not presented in ReadInDet _jobOptions.py by default):

InDetFlags.doTrkD3PD = True
InDetFlags.doVtxD3PD = True
from TrackD3PDMaker.TrackD3PDMakerFlags import TrackD3PDFlags
TrackD3PDFlags.storeTRTHitsOnTrack.set_Value_and_Lock(True)
TrackD3PDFlags.storeTRTOutliersOnTrack.set_Value_and_Lock(True)
TrackD3PDFlags.storeTRTHolesOnTrack.set_Value_and_Lock(True)
print TrackD3PDFlags

Producing D3PDs

For now, all you need is a job option file. An example, TRTTrackD3PD_jobOptions.py, is attached at the bottom of this page. Edit this file and run.

$ athena.py TRTTrackD3PD_jobOptions.py
This create InDetTrackD3PD.root

Making & running the Analysis Code

Compilable Analysis Code to run over tracking D3PDs

It contains three C++ files/classes:

InDetTrackTree.h/.C
RealTrackTree.h/.C
TrackingAna.C

and a make and a README file.

To compile it setup your root version and then type make

Input D3PD

TrackingAna.exe require an input D3PD as specified in TrackingAna.C in your analysis code.

Run the TrackingAna.exe

Afrer making analysis code and setting proper input D3PD you can run the TrackingAna.exe ./TrackingAna.exe
TrackingAna.exe produce ROOT file containing tracking performance histograms:

  • Hit residuals
  • Hit occupancies
  • HT fraction

Useful Links


Responsible: AlexeyBoldyrev
Review: Never reviewed

Edit | Attach | Watch | Print version | History: r4 < r3 < r2 < r1 | Backlinks | Raw View | WYSIWYG | More topic actions
Topic revision: r4 - 2012-06-23 - AlexeyBoldyrev
 
    • Cern Search Icon Cern Search
    • TWiki Search Icon TWiki Search
    • Google Search Icon Google Search

    Sandbox All webs login

This site is powered by the TWiki collaboration platform Powered by PerlCopyright & 2008-2020 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