TrkAna
Overview
TrkAna is a track-based ROOT TTree that can be used to help with analysis. Each entry in the TTree corresponds to a single fitted track and contains reconstructed information from the tracker, calorimeter and CRV. There is also the option to write out Monte Carlo truth information and additional information for other track types that might be important to an analysis.
There is a tutorial and a hypernews forum you can sign up to.
How To Run
There is an example fcl configuration in the Offline repository that runs on reco datasets.
mu2e -c TrkDiag/fcl/TrkAnaReco.fcl -s mcs.art --TFileName trkana.root
This runs two instances of the TrackAnalysisReco module. One for negative tracks (TrkAnaNeg) and one for positive tracks (TrkAnaPos). Each tree has the same structure has shown below.
Tree Structure
Here is a very rough description of the tree branches and where to find leaf definitions in the repository:
- evtinfo
- information about the event (TrkDiag/inc/EventInfo.hh)
- hcnt
- count of various hit types (TrkDiag/inc/HitCount.hh)
- tcnt
- count of various track types (TrkDiag/inc/TrkCount.hh)
- de
- global fit information for downstream electron track (TrkDiag/inc/TrkInfo.hh)
- deent
- local fit information for downstream electron track at tracker entrance (TrkDiag/inc/TrkInfo.hh)
- demid
- local fit information for downstream electron track at middle of tracker (TrkDiag/inc/TrkInfo.hh)
- dexit
- local fit information for downstream electron track at tracker exit (TrkDiag/inc/TrkInfo.hh)
- detch
- calorimeter cluster information for cluster used in downstream electron track fit (tch = TrkCaloHit, TrkDiag/inc/TrkCaloHitInfo.hh)
- dequal
- the output values of the TrkQual and TrkPID ANNs
- uetch
- calorimeter cluster information for cluster used in upstream electron track fit (tch = TrkCaloHit, TrkDiag/inc/TrkCaloHitInfo.hh)
- ue
- global fit information for upstream electron track (TrkDiag/inc/TrkInfo.hh)
- dm
- global fit information for downstream muon track (TrkDiag/inc/TrkInfo.hh)
- crvinfo
- information about CRV coincidences (CRVAnalysis/inc/CrvHitInfoReco.hh)
- bestcrv
- element in crvinfo array that is the best
- demc
- MC information about particle that created downstream track (TrkDiag/inc/TrkInfo.hh)
- demcgen
- MC information about particle that started the simulation (TrkDiag/inc/GenInfo.hh)
- demcpri
- MC information about particle that would have ultimately created the GenParticle (TrkDiag/inc/GenInfo.hh)
- demcent
- MC information about step of particle that created downstream track as it enters the tracker (TrkDiag/inc/TrkInfo.hh)
- demcmid
- MC information about step of particle that created downstream track as it passes the middle of the tracker (TrkDiag/inc/TrkInfo.hh)
- demcxit
- MC information about step of particle that created downstream track as it leaves the tracker (TrkDiag/inc/TrkInfo.hh)
- crvinfomc
- MC information about CRV coincidences (CRVAnalysis/inc/CrvHitInfoMC.hh)
- detchmc
- MC information about calorimeter cluster used in the downstream electron fit (TrkDiag/inc/CaloClusterInfoMC.hh)
- uetchmc
- MC information about calorimeter cluster used in the upstream electron fit (TrkDiag/inc/CaloClusterInfoMC.hh)
- detshmc
- MC information about the straw hits used in the downstream electron fit (need diagLevel > 0, TrkDiag/inc/TrkStrawHitInfo.hh)
- detrkqual
- the input variables and output value of the track quality artificial neural network (TrkDiag/inc/TrkQualInfo.hh)
- evtwt
- the values of all EventWeight objects that were in the art even (e.g. proton bunch intensity = PBIWeight)
By setting the diagLevel
to 2, you can get hit level information:
- detsh
- reconstructed straw hit information for downstream electron track fit (TrkDiag/inc/TrkStrawHitInfo.hh)
- detsm
- information about straw materials that the downstream electron track fit goes through (TrkDiag/inc/TrkStrawMatInfo.hh)
Future Developments?
- Support for TTreeReader and RDataFrame?