FastTimeDDDConstants.cc

Go to the documentation of this file.

#include "Geometry/HGCalCommonData/interface/FastTimeDDDConstants.h"

#include "DataFormats/Math/interface/GeantUnits.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"

//#define EDM_ML_DEBUG
using namespace geant_units::operators;

FastTimeDDDConstants::FastTimeDDDConstants(const FastTimeParameters* ft) : ftpar_(ft) {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants::FastTimeDDDConstants "
                                << "( const FastTimeParameters* ft ) constructor";
#endif
  initialize();
}

FastTimeDDDConstants::~FastTimeDDDConstants() {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:destructed!!!";
#endif
}

std::pair<int, int> FastTimeDDDConstants::getZPhi(double z, double phi) const {
  if (phi < 0)
    phi += (2 * geant_units::piRadians);
  int iz = (int)(z / dZBarrel_) + 1;
  if (iz > ftpar_->nZBarrel_)
    iz = ftpar_->nZBarrel_;
  int iphi = (int)(phi / dPhiBarrel_) + 1;
  if (iphi > ftpar_->nPhiBarrel_)
    iphi = 1;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:Barrel z|phi " << z << " " << convertRadToDeg(phi)
                                << " iz|iphi " << iz << " " << iphi;
#endif
  return std::pair<int, int>(iz, iphi);
}

std::pair<int, int> FastTimeDDDConstants::getEtaPhi(double r, double phi) const {
  if (phi < 0)
    phi += (2 * geant_units::piRadians);
  int ir(ftpar_->nEtaEndcap_);
  for (unsigned int k = 1; k < rLimits_.size(); ++k) {
    if (r > rLimits_[k]) {
      ir = k;
      break;
    }
  }
  int iphi = (int)(phi / dPhiEndcap_) + 1;
  if (iphi > ftpar_->nPhiEndcap_)
    iphi = 1;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants:Endcap r|phi " << r << " " << convertRadToDeg(phi)
                                << " ir|iphi " << ir << " " << iphi;
#endif
  return std::pair<int, int>(ir, iphi);
}

GlobalPoint FastTimeDDDConstants::getPosition(int type, int izeta, int iphi, int zside) const {
  double x(0), y(0), z(0);
  if (type == 1) {
    double phi = (iphi - 0.5) * dPhiBarrel_;
    x = ftpar_->geomParBarrel_[2] * cos(phi);
    y = ftpar_->geomParBarrel_[2] * sin(phi);
    z = ftpar_->geomParBarrel_[0] + (izeta - 0.5) * dZBarrel_;
  } else if (type == 2) {
    double phi = (iphi - 0.5) * dPhiEndcap_;
    double r = (izeta <= 0 || izeta >= (int)(rLimits_.size())) ? 0 : 0.5 * (rLimits_[izeta - 1] + rLimits_[izeta]);
    x = (zside < 0) ? -r * cos(phi) : r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParEndcap_[2];
  }
  if (zside < 0)
    z = -z;
  GlobalPoint p(x, y, z);
  return p;
}

std::vector<GlobalPoint> FastTimeDDDConstants::getCorners(int type, int izeta, int iphi, int zside) const {
  double x(0), y(0), z(0), dx(0), dz(0), r(0), phi(0);
  if (type == 1) {
    phi = (iphi - 0.5) * dPhiBarrel_;
    r = ftpar_->geomParBarrel_[2];
    x = r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParBarrel_[0] + (izeta - 0.5) * dZBarrel_;
    dx = 0.5 * ftpar_->geomParBarrel_[3];
    dz = 0.5 * dZBarrel_;
  } else if (type == 2) {
    phi = (iphi - 0.5) * dPhiEndcap_;
    r = (izeta <= 0 || izeta >= (int)(rLimits_.size())) ? 0 : 0.5 * (rLimits_[izeta - 1] + rLimits_[izeta]);
    x = (zside < 0) ? -r * cos(phi) : r * cos(phi);
    y = r * sin(phi);
    z = ftpar_->geomParEndcap_[2];
    dx = 0.5 * r * dPhiEndcap_;
    dz = 0.5 * ftpar_->geomParEndcap_[3];
  }
  if (zside < 0) {
    z = -z;
    dz = -dz;
  }
  static constexpr int signx[8] = {-1, -1, 1, 1, -1, -1, 1, 1};
  static constexpr int signy[8] = {-1, 1, 1, -1, -1, 1, 1, -1};
  static constexpr int signz[8] = {-1, -1, -1, -1, 1, 1, 1, 1};
  std::vector<GlobalPoint> pts;
  for (unsigned int i = 0; i != 8; ++i) {
    GlobalPoint p(x + signx[i] * dx, y + signy[i] * dx, z + signz[i] * dz);
    pts.emplace_back(p);
  }
  return pts;
}

int FastTimeDDDConstants::getCells(int type) const {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nZBarrel_) * (ftpar_->nPhiBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nEtaEndcap_) * (ftpar_->nPhiEndcap_);
  }
  return numb;
}

double FastTimeDDDConstants::getRin(int type) const {
  double value(0);
  if (type == 1) {
    value = (ftpar_->geomParBarrel_[2]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[0]);
  }
  return value;
}

double FastTimeDDDConstants::getRout(int type) const {
  double value(0);
  if (type == 1) {
    value = (ftpar_->geomParBarrel_[2] + ftpar_->geomParBarrel_[3]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[1]);
  }
  return value;
}

double FastTimeDDDConstants::getZHalf(int type) const {
  double value(0);
  if (type == 1) {
    value = 0.5 * (ftpar_->geomParBarrel_[1] - ftpar_->geomParBarrel_[0]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[3]);
  }
  return value;
}

double FastTimeDDDConstants::getZPos(int type) const {
  double value(0);
  if (type == 1) {
    value = 0.5 * (ftpar_->geomParBarrel_[1] + ftpar_->geomParBarrel_[0]);
  } else if (type == 2) {
    value = (ftpar_->geomParEndcap_[2]);
  }
  return value;
}

bool FastTimeDDDConstants::isValidXY(int type, int izeta, int iphi) const {
  bool ok(false);
  if (type == 1) {
    ok = ((izeta > 0) && (izeta <= ftpar_->nZBarrel_) && (iphi > 0) && (iphi <= ftpar_->nPhiBarrel_));
  } else if (type == 2) {
    ok = ((izeta > 0) && (izeta <= ftpar_->nEtaEndcap_) && (iphi > 0) && (iphi <= ftpar_->nPhiEndcap_));
  }
  return ok;
}

int FastTimeDDDConstants::numberEtaZ(int type) const {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nZBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nEtaEndcap_);
  }
  return numb;
}

int FastTimeDDDConstants::numberPhi(int type) const {
  int numb(0);
  if (type == 1) {
    numb = (ftpar_->nPhiBarrel_);
  } else if (type == 2) {
    numb = (ftpar_->nPhiEndcap_);
  }
  return numb;
}

void FastTimeDDDConstants::initialize() {
  double thmin = atan(ftpar_->geomParEndcap_[0] / ftpar_->geomParEndcap_[2]);
  etaMax_ = -log(0.5 * thmin);
  double thmax = atan(ftpar_->geomParEndcap_[1] / ftpar_->geomParEndcap_[2]);
  etaMin_ = -log(0.5 * thmax);
  dEta_ = (etaMax_ - etaMin_) / ftpar_->nEtaEndcap_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Theta range " << convertRadToDeg(thmin) << ":" << convertRadToDeg(thmax)
                                << " Eta range " << etaMin_ << ":" << etaMax_ << ":" << dEta_;
#endif
  for (int k = 0; k <= ftpar_->nEtaEndcap_; ++k) {
    double eta = etaMin_ + k * dEta_;
    double theta = 2.0 * atan(exp(-eta));
    double rval = (ftpar_->geomParEndcap_[2]) * tan(theta);
    rLimits_.emplace_back(rval);
  }
  dZBarrel_ = ftpar_->geomParBarrel_[1] / ftpar_->nZBarrel_;
  dPhiBarrel_ = (2 * geant_units::piRadians) / ftpar_->nPhiBarrel_;
  dPhiEndcap_ = (2 * geant_units::piRadians) / ftpar_->nPhiEndcap_;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "FastTimeDDDConstants initialized with " << ftpar_->nZBarrel_ << ":"
                                << ftpar_->nPhiBarrel_ << ":" << getCells(1) << " cells for barrel; dz|dphi "
                                << dZBarrel_ << "|" << dPhiBarrel_ << " and " << ftpar_->nEtaEndcap_ << ":"
                                << ftpar_->nPhiEndcap_ << ":" << getCells(2) << " cells for endcap; dphi "
                                << dPhiEndcap_ << " The Limits in R are";
  for (unsigned int k = 0; k < rLimits_.size(); ++k)
    edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << rLimits_[k] << " ";
#endif
}

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

TYPELOOKUP_DATA_REG(FastTimeDDDConstants);