TRTPlots.h

Go to the documentation of this file.

/* emacs: this is -*- c++ -*- */
 
 
#ifndef  IDTPM2_TRTPLOTS_H
#define  IDTPM2_TRTPLOTS_H
 
#include "IPlotBase.h"

 
 
  
#define TRT_PARAMS(X)  \
  X(NTRTHits)           \
  X(NTRTHitsXe)         \
  X(NTRTHitsAr)         \
  X(NTRTHighThresholdHits)  \
  X(NTRTHighThresholdHitsXe)    \
  X(NTRTHighThresholdHitsAr)    \
                \
  X(NTRTOutliers)       \
  X(NTRTHighThresholdOutliers)
 
    
namespace IDTPM2 {
  
class TRTPlots : public IPlotBase { 
    
public:
 
  TRTPlots( const std::string& name, ToolHandle<GenericMonitoringTool>* m );
  
  virtual ~TRTPlots() { } 
 
  void initialize();
 
  void execute( const TrackAnalysisStore& collections ) const;
  
  void finalize();
 
private:


  
  struct Histograms {


 
    enum TRTPar {
#      define X(x) x,
       TRT_PARAMS(X)
#      undef X
       NTRTPars
    };
    
    std::vector<std::string> m_parNames = {
#     define X(x) #x,
      TRT_PARAMS(X)
#     undef X
    };
 
    size_t            nParNames() const { return m_parNames.size(); } 
 
    const std::string&  parName(int i) const { return m_parNames.at(i); }
    
    void initialize( TRTPlots* p, const std::string& s="" );
    
    void  fill( const TrackView* t ) const;
 
#if 1
 
    TIDA::Histogram<float> m_hits[NTRTPars];

 
    TIDA::Histogram<float>  m_hits_vs_pt[NTRTPars];
    TIDA::Histogram<float>  m_hits_vs_eta[NTRTPars];
    TIDA::Histogram<float>  m_hits_vs_phi[NTRTPars];
    TIDA::Histogram<float>  m_hits_vs_d0[NTRTPars];
 
#else   

    
#   define X(x) TIDA::Histogram<float> m_##x;
    TRT_PARAMS(X)
#   undef X
 
#   define X(x)  TIDA::Histogram2D<float> m_##x##_vs_pt;
    TRT_PARAMS(X)
#   undef X
 
#   define X(x)  TIDA::Histogram2D<float> m_##x##_vs_eta;
    TRT_PARAMS(X)
#   undef X
 
#   define X(x) TIDA::Histogram2D<float> m_##x##_vs_phi;
    TRT_PARAMS(X)
#   undef X
 
#   define X(x) TIDA::Histogram2D<float> m_##x##_vs_d0;
    TRT_PARAMS(X)
#   undef X
  
#endif
  
  };
 
private:

  
  Histograms m_reference;
  Histograms m_test;
  
};
 
} // namespace IDTPM2
 
 
inline std::ostream& operator<<( std::ostream& s, const IDTPM2::TRTPlots&  ) { 
  return s;
}
 
 
#endif  // IDTPM2_TRTPLOTS_H