TrigInDetTrackMonitoringTool.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022CERN for the benefit of the ATLAS collaboration
*/
 
#include "TMath.h"
#include "TrigInDetTrackMonitoringTool.h"
#include<iostream>
 
//#include "TrigInterfaces/Algo.h"
//
//Constructor-------------------------------------------------------------- 
TrigInDetTrackMonitoringTool::TrigInDetTrackMonitoringTool(const std::string& type,
                                           const std::string& name,
                                           const IInterface* parent):
  AthAlgTool(type,name,parent)
{
  //declareProperty( "MonitoringTool",      m_monTool  );
}
 
//Destructor---------------------------------------------------------------
TrigInDetTrackMonitoringTool::~TrigInDetTrackMonitoringTool(){
  ATH_MSG_DEBUG( "TrigInDetTrackMonitoringTool destructor called");
}
 
//Initialize---------------------------------------------------------------
StatusCode TrigInDetTrackMonitoringTool::initialize(){
  ATH_MSG_DEBUG( "TrigInDetTrackMonitoringTool initialize() called");
 
  //Retrieve monitoring tool
  if (!m_monTool.empty()) ATH_CHECK(m_monTool.retrieve());
   
  return StatusCode::SUCCESS;
 
}
 
StatusCode TrigInDetTrackMonitoringTool::finalize() {
  ATH_MSG_DEBUG("TrigInDetTrackMonitoringTool finalize()");
  return StatusCode::SUCCESS; 
}
 
 
//Monitoring track variables
//void TrigInDetTrackMonitoringTool::monitor_tracks(std::string prefix, std::string suffix, DataVector<xAOD::TrackParticle> tracks ){
//void TrigInDetTrackMonitoringTool::monitor_tracks( const std::string &prefix, const std::string &suffix, const xAOD::TrackParticleContainer &tracks ){
void TrigInDetTrackMonitoringTool::monitor_tracks( const std::string &prefix, const std::string &suffix, const xAOD::TrackParticleContainer &tracks ) const {
   auto  trackPt      = Monitored::Collection(  prefix + "Pt"       + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->pt() * 0.001;});
   auto  trackQOverP  = Monitored::Collection(  prefix + "QOverP"   + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->qOverP() * 1000;});
   auto  trackQPt     = Monitored::Collection(  prefix + "QPt"      + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->pt() * 0.001 * (t->qOverP() > 0 ? 1.0 : -1.0);});
   auto  trackEta     = Monitored::Collection(  prefix + "Eta"      + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->eta();});
   auto  trackPhi     = Monitored::Collection(  prefix + "Phi"      + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->phi();});
   auto  trackTheta   = Monitored::Collection(  prefix + "Theta"    + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->theta();});
   auto  trackZ0      = Monitored::Collection(  prefix + "Z0"       + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->z0();});
   auto  trackD0      = Monitored::Collection(  prefix + "D0"       + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->d0();});
//TODO: append later if needed
//   auto  trackZ0err   = Monitored::Collection(  prefix + "Z0err"    + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->Z0err()   ;});
//   auto  trackD0err   = Monitored::Collection(  prefix + "D0err"    + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->D0err()   ;});
//   auto  trackQual    = Monitored::Collection(  prefix + "Qual"     + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->Qual()    ;});
//
   auto  trackNDF     = Monitored::Collection(  prefix + "NDF"      + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->numberDoF();});
   auto  trackChi2Prob= Monitored::Collection(  prefix + "Chi2Prob" + suffix, tracks, [](const xAOD::TrackParticle *t){ return t->chiSquared();});
 
   //Hits in the layiers
   auto  trackPiHits  = Monitored::Collection(  prefix + "PIXHits"  + suffix, tracks, [](const xAOD::TrackParticle *t){ uint8_t nPixHits = 0; t->summaryValue( nPixHits,  xAOD::numberOfPixelHits); return nPixHits ;});
   auto  trackSCTHits = Monitored::Collection(  prefix + "SCTHits"  + suffix, tracks, [](const xAOD::TrackParticle *t){ uint8_t nSCTHits = 0; t->summaryValue( nSCTHits,  xAOD::numberOfSCTHits); return nSCTHits ;});
 
   auto  trackTRTHits = Monitored::Collection(  prefix + "TRTHits"  + suffix, tracks, [](const xAOD::TrackParticle *t){ uint8_t nTRTHits = 0; t->summaryValue( nTRTHits,  xAOD::numberOfTRTHits); return nTRTHits ;});
 
   // This track counting still doesn't fill histogram
   auto trkCount  = Monitored::Scalar( prefix + "Counting" + suffix, tracks.size());
 
   //If needed add more monitored values from:
   //https://gitlab.cern.ch/atlas/athena/-/blob/master/Tracking/TrkTools/TrkParticleCreator/src/TrackParticleCreatorTool.cxx
 
   auto mon = Monitored::Group(m_monTool, 
              trackPt, trackEta, trackQOverP, trackQPt, trackPhi, trackTheta, trackD0 ,trackZ0, trackNDF, trackChi2Prob, 
              trackPiHits, trackSCTHits, trackTRTHits, //Hits in the det layers
              trkCount
              ); 
}