SiHitPlots.h

Go to the documentation of this file.

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

 
#define HIT_PARAMS(X)            \
    X(NInnerMostPixelHits)           \
    X(NInnerMostPixelEndcapHits)     \
    X(NNextToInnerMostPixelHits)     \
    X(NNextToInnerMostPixelEndcapHits)   \
    X(NInnerMostPixelSharedHits)     \
    X(NInnerMostPixelSharedEndcapHits)   \
    X(NPixelHits)            \
    X(NPixelHoles)               \
    X(NPixelSharedHits)              \
    X(NPixeldEdx)            \
    X(NSCTHits)                  \
    X(NSCTHoles)             \
    X(NSCTSharedHits)                \
                         \
    X(NInnerMostPixelOutliers)           \
    X(NInnerMostPixelEndcapOutliers)     \
    X(NInnerMostPixelSplitHits)          \
    X(NInnerMostPixelSplitEndcapHits)    \
    X(NExpectedInnerMostPixelHits)   \
    X(NExpectedNextToInnerMostPixelHits) \
    X(NPixelOutliers)                \
    X(NPixelContribLayers)           \
    X(NPixelSplitHits)               \
    X(NPixelGangedHits)              \
    X(NPixelGangedHitsFlaggedFakes)      \
    X(NPixelDeadSensors)             \
    X(NSCTOutliers)              \
    X(NSCTDoubleHoles)               \
    X(NSCTDeadSensors)
    
 
namespace IDTPM2 {
  
class SiHitPlots : public IPlotBase { 
    
public:
 
  SiHitPlots( const std::string& name, ToolHandle<GenericMonitoringTool>* m );
  
  virtual ~SiHitPlots() { } 
 
  void initialize();
 
  void execute( const TrackAnalysisStore& collections ) const;
  
  void finalize();
 
private:


  
  struct Histograms {


 
    enum HitPar {
#      define X(x) x,
       HIT_PARAMS(X)
#      undef X
       NHitPars
    };
    
    std::vector<std::string> m_hitParNames = {
#     define X(x) #x,
      HIT_PARAMS(X)
#     undef X
    };
 
    size_t                 nHitParNames() const { return m_hitParNames.size(); } 
 
    const std::string&  hitParName(int i) const { return m_hitParNames.at(i); }
    
    void initialize( SiHitPlots* p, const std::string& s="" );
    
    void  fill( const TrackView* t ) const;
 
#if 1
 
    TIDA::Histogram<float> m_hits[NHitPars];

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

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

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