SiStripDetSummary.h

Go to the documentation of this file.

#ifndef SiStripDetSummary_h
#define SiStripDetSummary_h

#include "DataFormats/SiStripDetId/interface/TIDDetId.h" 
#include "DataFormats/SiStripDetId/interface/TECDetId.h" 
#include "DataFormats/SiStripDetId/interface/TIBDetId.h" 
#include "DataFormats/SiStripDetId/interface/TOBDetId.h" 
#include "DataFormats/DetId/interface/DetId.h"

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

class SiStripDetSummary
{
public:
  SiStripDetSummary() : computeMean_(true)
  {
    // 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(const DetId & detid, const float & value);
  inline void add(const DetId & detid)
  {
    computeMean_ = false;
    add( detid, 0 );
  }

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

  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_;
};

#endif