SiStripFakeAPVParameters.h

Go to the documentation of this file.

#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
#include "DataFormats/TrackerCommon/interface/TrackerTopology.h"

class SiStripFakeAPVParameters {
public:
  using index = std::pair<int,int>;

  SiStripFakeAPVParameters() {}

  SiStripFakeAPVParameters(const edm::ParameterSet& pset, const std::string& parameterName)
  {
    const int layersTIB = 4;
    const int ringsTID = 3;
    const int layersTOB = 6;
    const int ringsTEC = 7;

    fillSubDetParameter(pset.getParameter<std::vector<double>>(parameterName+"TIB"), int(StripSubdetector::TIB), layersTIB );
    fillSubDetParameter(pset.getParameter<std::vector<double>>(parameterName+"TID"), int(StripSubdetector::TID), ringsTID );
    fillSubDetParameter(pset.getParameter<std::vector<double>>(parameterName+"TOB"), int(StripSubdetector::TOB), layersTOB );
    fillSubDetParameter(pset.getParameter<std::vector<double>>(parameterName+"TEC"), int(StripSubdetector::TEC), ringsTEC );
  }

  inline double get(const index& idx) const
  {
    return m_data.at(idx.first)[idx.second];
  }

  static index getIndex(const TrackerTopology* tTopo, DetId id)
  {
    int layerId{0};
    const int subId =  StripSubdetector(id).subdetId();
    switch(subId) {
      case int(StripSubdetector::TIB):
        layerId = tTopo->tibLayer(id) - 1;
        break;
      case int(StripSubdetector::TOB):
        layerId = tTopo->tobLayer(id) - 1;
        break;
      case int(StripSubdetector::TID):
        layerId = tTopo->tidRing(id) - 1;
        break;
      case int(StripSubdetector::TEC):
        layerId = tTopo->tecRing(id) - 1;
        break;
      default:
        break;
    }
    return std::make_pair(subId, layerId);
  }
private:
  using LayerParameters = std::vector<double>;
  using SubdetParameters = std::map<int,LayerParameters>;
  SubdetParameters m_data;

  void fillSubDetParameter(const std::vector<double>& v, const int subDet, const unsigned short layers)
  {
    if( v.size() == layers ) {
      m_data.insert(std::make_pair( subDet, v ));
    }
    else if( v.size() == 1 ) {
      LayerParameters parV(layers, v[0]);
      m_data.insert(std::make_pair( subDet, parV ));
    }
    else {
      throw cms::Exception("Configuration") << "ERROR: number of parameters for subDet " << subDet << " are " << v.size() << ". They must be either 1 or " << layers << std::endl;
    }
  }
};