TRTPlots.cxx

Go to the documentation of this file.

 
#include "TRTPlots.h"
 
IDTPM2::TRTPlots::TRTPlots( const std::string& name, ToolHandle<GenericMonitoringTool>* m )
  : IPlotBase( name, m ) { 
}
 

 
void IDTPM2::TRTPlots::initialize() {
  m_reference.initialize( this, "reference");
  m_test.initialize(      this, "test");
 
}
 
 
 
void IDTPM2::TRTPlots::Histograms::initialize( TRTPlots* p, const std::string& s ) {
 
  static const std::array<bool,NTRTPars> proc = { true, false, false, true };
  
  for ( size_t i=0 ; i<nParNames() ; i++ ) {
 
    if ( !proc[i] ) continue;
    
    m_hits[i] = TIDA::Histogram<float>( p->monTool(), parName(i), s );
  
 
    m_hits_vs_pt[i]  = TIDA::Histogram<float>( p->monTool(), parName(i)+"_vs_pt",  s );
    m_hits_vs_eta[i] = TIDA::Histogram<float>( p->monTool(), parName(i)+"_vs_eta", s );
    m_hits_vs_phi[i] = TIDA::Histogram<float>( p->monTool(), parName(i)+"_vs_phi", s );
    m_hits_vs_d0[i]  = TIDA::Histogram<float>( p->monTool(), parName(i)+"_vs_d0",  s );
 
  }
  
}
 
 
void  IDTPM2::TRTPlots::Histograms::fill( const TrackView* t ) const {
  

  std::array<float,NTRTPars> hitP = {};
 
  hitP[NTRTHits] = t->nTRT();
  //  hitP[NTRTHitsXe]   
  //  hitP[ NTRTHitsAr]  
  hitP[NTRTHighThresholdHits]  = t->nTRTHi(); 
  //  hitP[NTRTHighThresholdHitsXe] 
  //  hitP[NTRTHighThresholdHitsAr] 
  
  //  hitP[NTRTOutliers]
  //  hitP[NTRTHighThresholdOutliers] 
  
    
  float pt  = std::fabs(t->pt())*0.001;
  float eta = t->eta();
  float phi = t->phi();
  float d0  = t->d0();
 
  
  //  std::cout << "TRTs: track: " << *t << std::endl; 
 
  static const std::array<bool,NTRTPars> proc = { true, false, false, true };
  
  for ( int i=0 ; i<NTRTPars ; i++ ) {
 
    if ( !proc[i] ) continue;
    
    //    std::cout << "\tfill par; " << hitParName(i) << "\t: " << hitP[i] << std::endl;
    
    m_hits[i].Fill( hitP[i] );
 
    m_hits_vs_pt[i].Fill(  pt,  hitP[i] );
    m_hits_vs_eta[i].Fill( eta, hitP[i] );
    m_hits_vs_phi[i].Fill( phi, hitP[i] );
    m_hits_vs_d0[i].Fill(  d0,  hitP[i] );
  }
  
}
 
 
 
void IDTPM2::TRTPlots::execute( const TrackAnalysisStore& collections ) const {
  
#if 0
  
  double  mu = 0;
  
  const xAOD::EventInfo* eventinfo = collections.eventInfo();
  
  if (eventinfo) {
    mu  = collections.eventInfo()->averageInteractionsPerCrossing();
  }
#endif
  
  const std::vector<const TrackView*>& ref  = collections.reference_tracks();
  const std::vector<const TrackView*>& test = collections.test_tracks();
 
  for ( size_t i=0 ;  i<ref.size() ; i++ ) {  m_reference.fill(  ref[i] ); }
  for ( size_t i=0 ; i<test.size() ; i++ ) {       m_test.fill( test[i] ); }
  
 
}
  
  
void IDTPM2::TRTPlots::finalize() { }