#include "Geometry/GEMGeometry/interface/ME0EtaPartitionSpecs.h"
#include "Geometry/CommonTopologies/interface/RectangularStripTopology.h"
#include "Geometry/CommonTopologies/interface/TrapezoidalStripTopology.h"

using namespace GeomDetEnumerators;

ME0EtaPartitionSpecs::ME0EtaPartitionSpecs(SubDetector rss, const std::string& name, const ME0Specs& pars)
    : GeomDetType(name, rss), _p(pars), _n(name) {
  if (rss == ME0) {
    float b = _p[0];
    float B = _p[1];
    float h = _p[2];
    float r0 = h * (B + b) / (B - b);
    float striplength = h * 2;
    float strips = _p[3];
    float pitch = (b + B) / strips;
    int nstrip = static_cast<int>(strips);
    _top = new TrapezoidalStripTopology(nstrip, pitch, striplength, r0);

    float pads = _p[4];
    float pad_pitch = (b + B) / pads;
    int npad = static_cast<int>(pads);
    _top_pad = new TrapezoidalStripTopology(npad, pad_pitch, striplength, r0);
  } else {
    _top = nullptr;
    _top_pad = nullptr;
  }
}

ME0EtaPartitionSpecs::~ME0EtaPartitionSpecs() {
  if (_top)
    delete _top;
  if (_top_pad)
    delete _top_pad;
}

const Topology& ME0EtaPartitionSpecs::topology() const { return *_top; }

const StripTopology& ME0EtaPartitionSpecs::specificTopology() const { return *_top; }

const Topology& ME0EtaPartitionSpecs::padTopology() const { return *_top_pad; }

const StripTopology& ME0EtaPartitionSpecs::specificPadTopology() const { return *_top_pad; }

const std::string& ME0EtaPartitionSpecs::detName() const { return _n; }

const std::vector<float>& ME0EtaPartitionSpecs::parameters() const { return _p; }