SiStripDetSummary.h

Go to the documentation of this file.

#ifndef SiStripDetSummary_h
#define SiStripDetSummary_h

#include "DataFormats/DetId/interface/DetId.h"
class TrackerTopology;

#include <sstream>
#include <map>
#include <cmath>
#include <iomanip>
#include <iostream>

class SiStripDetSummary {
public:
  explicit SiStripDetSummary(const TrackerTopology* tTopo) : computeMean_(true), trackerTopo_(tTopo) {
    // Initialize valueMap_ with zeros
    // WARNING: this initialization is strongly connected with how the map is filled in the add method
    // TIB: layers = 4, stereo = the first 2
    // TOB: layers = 6, stereo = the first 2
    // TEC: wheels = 9, stereo = 9
    // TID: wheels = 3, stereo = 3
    unsigned int layers[] = {4, 6, 9, 3};
    unsigned int stereo[] = {2, 2, 9, 3};
    Values initValues;
    for (unsigned int subDet = 0; subDet < 4; ++subDet) {
      // Layers start from 1
      for (unsigned int layer = 1; layer <= layers[subDet]; ++layer) {
        valueMap_[1000 * (subDet + 1) + layer * 10] = initValues;
        if (layer <= stereo[subDet])
          valueMap_[1000 * (subDet + 1) + layer * 10 + 1] = initValues;
      }
    }
  }

  void add(DetId detid, float value);
  inline void add(DetId detid) {
    computeMean_ = false;
    add(detid, 0);
  }

  void print(std::stringstream& ss, const bool mean = true) const;

  inline void clear() { valueMap_.clear(); }

  struct Values {
    Values() : mean(0.), rms(0.), count(0) {}
    double mean;
    double rms;
    unsigned int count;
  };
  std::map<unsigned int, Values> getCounts() { return valueMap_; }

protected:
  // Maps to store the value and the counts
  std::map<unsigned int, Values> valueMap_;
  bool computeMean_;

private:
  const TrackerTopology* trackerTopo_;
};

#endif