Trapezoid2RectangleMappingX.h

Go to the documentation of this file.

#ifndef Trapezoid2RectangleMappingX_h
#define Trapezoid2RectangleMappingX_h

//#define DEBUG_GRID_TRM

#ifdef DEBUG_GRID_TRM
#include <iostream>
#endif

#include "FWCore/Utilities/interface/Visibility.h"

class dso_internal Trapezoid2RectangleMappingX {
public:
  Trapezoid2RectangleMappingX() {}

  Trapezoid2RectangleMappingX(double x0, double y0, double bovera, double h) : x0_(x0), y0_(y0), parallel_(false) {
    k_ = 2 / h * (bovera - 1.) / (bovera + 1.);

#ifdef DEBUG_GRID_TRM
    std::cout << "Trapezoid2RectangleMappingX constructed with x0,y0 " << x0 << " " << y0 << " b/a= " << bovera
              << " h= " << h << std::endl;
#endif
  }

  Trapezoid2RectangleMappingX(double x0, double y0, double k) : x0_(x0), y0_(y0), k_(k), parallel_(true) {
#ifdef DEBUG_GRID_TRM
    std::cout << "Trapezoid2RectangleMappingX constructed with x0,y0 " << x0 << " " << y0 << " k= " << k << std::endl;
#endif
  }

  void rectangle(double xtrap, double ytrap, double& xrec, double& yrec) const {
    yrec = ytrap - y0_;
    if (!parallel_)
      xrec = (xtrap - x0_) / (1. + yrec * k_);
    else
      xrec = xtrap - x0_ + k_ * yrec;

#ifdef DEBUG_GRID
    std::cout << xtrap << " " << ytrap << " transformed to rectangle " << xrec << " " << yrec << std::endl;
#endif
  }

  void trapezoid(double xrec, double yrec, double& xtrap, double& ytrap) const {
    if (!parallel_)
      xtrap = x0_ + xrec * (1. + yrec * k_);
    else
      xtrap = x0_ + xrec - k_ * yrec;
    ytrap = y0_ + yrec;

#ifdef DEBUG_GRID
    std::cout << xrec << " " << yrec << " transformed to trapezoid " << xtrap << " " << ytrap << std::endl;
#endif
  }

private:
  double x0_;
  double y0_;
  double k_;
  bool parallel_;
};

#endif