CylinderFromSectorMFGrid.cc

Go to the documentation of this file.

#include "CylinderFromSectorMFGrid.h"
#include "MagneticField/VolumeGeometry/interface/MagExceptions.h"
#include <iostream>
 
CylinderFromSectorMFGrid::CylinderFromSectorMFGrid(const GloballyPositioned<float>& vol,
                                                   double phiMin,
                                                   double phiMax,
                                                   MFGrid* sectorGrid)
    : MFGrid(vol),
      thePhiMin(phiMin),
      thePhiMax(phiMax),
      theSectorGrid(sectorGrid)
 
{
  if (thePhiMax < thePhiMin)
    thePhiMax += 2.0 * Geom::pi();
  theDelta = thePhiMax - thePhiMin;
}
 
CylinderFromSectorMFGrid::~CylinderFromSectorMFGrid() { delete theSectorGrid; }
 
MFGrid::LocalVector CylinderFromSectorMFGrid::valueInTesla(const LocalPoint& p) const {
  double phi = p.phi();
  if (phi < thePhiMax && phi > thePhiMin)
    return theSectorGrid->valueInTesla(p);
  else {
    double phiRot = floor((phi - thePhiMin) / theDelta) * theDelta;
    double c = cos(phiRot);
    double s = sin(phiRot);
    double xrot = p.x() * c + p.y() * s;
    double yrot = -p.x() * s + p.y() * c;
 
    // get field in interpolation sector
    MFGrid::LocalVector tmp = theSectorGrid->valueInTesla(LocalPoint(xrot, yrot, p.z()));
 
    // rotate field back to original sector
    return MFGrid::LocalVector(tmp.x() * c - tmp.y() * s, tmp.x() * s + tmp.y() * c, tmp.z());
  }
}
 
void CylinderFromSectorMFGrid::throwUp(const char* message) const {
  std::cout << "Throwing exception " << message << std::endl;
  throw MagGeometryError(message);
}
void CylinderFromSectorMFGrid::toGridFrame(const LocalPoint& p, double& a, double& b, double& c) const {
  throwUp("Not implemented yet");
}
 
MFGrid::LocalPoint CylinderFromSectorMFGrid::fromGridFrame(double a, double b, double c) const {
  throwUp("Not implemented yet");
  return LocalPoint();
}
 
Dimensions CylinderFromSectorMFGrid::dimensions() const { return theSectorGrid->dimensions(); }
 
MFGrid::LocalPoint CylinderFromSectorMFGrid::nodePosition(int i, int j, int k) const {
  throwUp("Not implemented yet");
  return LocalPoint();
}
 
MFGrid::LocalVector CylinderFromSectorMFGrid::nodeValue(int i, int j, int k) const {
  throwUp("Not implemented yet");
  return LocalVector();
}