#include <TrapezoidalCylindricalMFGrid.h>

Inheritance diagram for TrapezoidalCylindricalMFGrid:

Public Member Functions
void	dump () const override

LocalPoint	fromGridFrame (double a, double b, double c) const override
	find grid coordinates for point. For debugging and validation only. More...

void	toGridFrame (const LocalPoint &p, double &a, double &b, double &c) const override
	find grid coordinates for point. For debugging and validation only. More...

	TrapezoidalCylindricalMFGrid (binary_ifstream &istr, const GloballyPositioned< float > &vol)

LocalVector	uncheckedValueInTesla (const LocalPoint &p) const override
	Interpolated field value at given point; does not check for exceptions. More...

Public Member Functions inherited from MFGrid3D
Dimensions	dimensions (void) const override

Indexes	index (const LocalPoint &p) const override

	MFGrid3D (const GloballyPositioned< float > &vol)

LocalPoint	nodePosition (int i, int j, int k) const override
	Position of node in local frame. More...

LocalVector	nodeValue (int i, int j, int k) const override
	Field value at node. More...

LocalVector	valueInTesla (const LocalPoint &p) const override
	Interpolated field value at given point. More...

Public Member Functions inherited from MFGrid
const GloballyPositioned< float > &	frame () const
	Local reference frame. More...

	MFGrid (const GloballyPositioned< float > &vol)

	~MFGrid () override

Public Member Functions inherited from MagneticFieldProvider< float >
virtual LocalVectorType	derivativeInTeslaPerMeter (const LocalPointType &p, int N) const

virtual int	hasDerivatives () const

virtual LocalVectorType	valueInTesla (const LocalPointType &p) const=0

virtual	~MagneticFieldProvider ()

Private Attributes
Trapezoid2RectangleMappingX	mapping_

Additional Inherited Members
Public Types inherited from MFGrid
typedef GloballyPositioned< float >::GlobalPoint	GlobalPoint

typedef GloballyPositioned< float >::GlobalVector	GlobalVector

typedef GloballyPositioned< float >::LocalPoint	LocalPoint

typedef GloballyPositioned< float >::LocalVector	LocalVector

Public Types inherited from MagneticFieldProvider< float >
typedef Point3DBase< float, GlobalTag >	GlobalPointType

typedef Vector3DBase< float, GlobalTag >	GlobalVectorType

typedef Point3DBase< float, LocalTag >	LocalPointType

typedef Vector3DBase< float, LocalTag >	LocalVectorType

Protected Types inherited from MFGrid3D
using	BVector = Grid3D::BVector

using	GridType = Grid3D

Protected Member Functions inherited from MFGrid3D
void	setGrid (const GridType &grid)

Protected Attributes inherited from MFGrid3D
GridType	grid_

Detailed Description

Definition at line 10 of file TrapezoidalCylindricalMFGrid.h.

Constructor & Destructor Documentation

◆ TrapezoidalCylindricalMFGrid()

TrapezoidalCylindricalMFGrid::TrapezoidalCylindricalMFGrid	(	binary_ifstream &	istr,
		const GloballyPositioned< float > &	vol
	)

Definition at line 10 of file TrapezoidalCylindricalMFGrid.cc.

References a, funct::abs(), b, gather_cfg::cout, dumpMFGeometry_cfg::delta, MFGrid::frame(), MFGrid3D::grid_, h, mps_fire::i, dqmiolumiharvest::j, mapping_, submitDQMOfflineCAF::nLines, Geom::pi(), Trapezoid2RectangleMappingX::rectangle(), and GloballyPositioned< T >::toLocal().

     : MFGrid3D(vol) {
   // The parameters read from the data files are given in global coordinates.
   // In version 85l, local frame has the same orientation of global frame for the reference
   // volume, i.e. the r.f. transformation is only a translation.
   // There is therefore no need to convert the field values to local coordinates.
   // Check this assumption:
   GlobalVector localXDir(frame().toGlobal(LocalVector(1, 0, 0)));
   GlobalVector localYDir(frame().toGlobal(LocalVector(0, 1, 0)));
 
   if (localXDir.dot(GlobalVector(1, 0, 0)) > 0.999999 && localYDir.dot(GlobalVector(0, 1, 0)) > 0.999999) {
     // "null" rotation - requires no conversion...
   } else {
     cout << "ERROR: TrapezoidalCylindricalMFGrid: unexpected orientation: x: " << localXDir << " y: " << localYDir
          << endl;
   }
 
   int n1, n2, n3;
   inFile >> n1 >> n2 >> n3;
   double xref, yref, zref;
   inFile >> xref >> yref >> zref;
   double stepx, stepy, stepz;
   inFile >> stepx >> stepy >> stepz;
 
   double BasicDistance1[3][3];  // linear step
   double BasicDistance2[3][3];  // linear offset
   bool easya, easyb, easyc;
 
   inFile >> BasicDistance1[0][0] >> BasicDistance1[1][0] >> BasicDistance1[2][0];
   inFile >> BasicDistance1[0][1] >> BasicDistance1[1][1] >> BasicDistance1[2][1];
   inFile >> BasicDistance1[0][2] >> BasicDistance1[1][2] >> BasicDistance1[2][2];
   inFile >> BasicDistance2[0][0] >> BasicDistance2[1][0] >> BasicDistance2[2][0];
   inFile >> BasicDistance2[0][1] >> BasicDistance2[1][1] >> BasicDistance2[2][1];
   inFile >> BasicDistance2[0][2] >> BasicDistance2[1][2] >> BasicDistance2[2][2];
   inFile >> easya >> easyb >> easyc;
 
   vector<BVector> fieldValues;
   float Bx, By, Bz;
   int nLines = n1 * n2 * n3;
   fieldValues.reserve(nLines);
   for (int iLine = 0; iLine < nLines; ++iLine) {
     inFile >> Bx >> By >> Bz;
     fieldValues.push_back(BVector(Bx, By, Bz));
   }
   // check completeness
   string lastEntry;
   inFile >> lastEntry;
   if (lastEntry != "complete") {
     cout << "ERROR during file reading: file is not complete" << endl;
   }
 
 #ifdef DEBUG_GRID
   cout << "easya " << easya << " easyb " << easyb << " easyc " << easyc << endl;
 #endif
 
   if (!easyb || !easyc) {
     throw MagGeometryError("TrapezoidalCartesianMFGrid only implemented for first coordinate");
   }
 
 #ifdef DEBUG_GRID
   cout << "Grid reference point in grid system: " << xref << "," << yref << "," << zref << endl;
   cout << "steps " << stepx << "," << stepy << "," << stepz << endl;
   cout << "ns " << n1 << "," << n2 << "," << n3 << endl;
 
   for (int i = 0; i < 3; ++i)
     for (int j = 0; j < 3; ++j) {
       cout << "BasicDistance1[" << i << "][" << j << "] = " << BasicDistance1[i][j] << "BasicDistance2[" << i << "]["
            << j << "] = " << BasicDistance2[i][j] << endl;
     }
 #endif
 
   // the "not easy" coordinate is x
   double a = stepx * (n1 - 1);
   double b = a + BasicDistance1[0][1] * (n2 - 1) * (n1 - 1) + BasicDistance1[0][2] * (n3 - 1) * (n1 - 1);
   //  double h = stepy * (n2-1);
   double h = stepz * (n3 - 1);
   double delta = -BasicDistance2[0][1] * (n2 - 1) - BasicDistance2[0][2] * (n3 - 1);
 
 #ifdef DEBUG_GRID
   cout << "Trapeze size (a,b,h) = " << a << "," << b << "," << h << endl;
 #endif
 
   GlobalPoint grefp(GlobalPoint::Cylindrical(xref, Geom::pi() - yref, zref));
   LocalPoint lrefp = frame().toLocal(grefp);
 
 #ifdef DEBUG_GRID
   cout << "Global origin " << grefp << endl;
   cout << "Local origin  " << lrefp << endl;
 #endif
 
   double baMinus1 = BasicDistance1[0][2] * (n3 - 1) / stepx;
   if (std::abs(baMinus1) > 0.000001) {
     double b_over_a = 1 + baMinus1;
     double a1 = std::abs(baMinus1) > 0.000001 ? delta / baMinus1 : a / 2;
 #ifdef DEBUG_GRID
     cout << "a1 = " << a1 << endl;
 #endif
 
     // transform reference point to grid frame
     double x0 = lrefp.perp() + a1;
     double y0 = lrefp.z() + h / 2.;
     mapping_ = Trapezoid2RectangleMappingX(x0, y0, b_over_a, h);
   } else {  // parallelogram
     mapping_ = Trapezoid2RectangleMappingX(0, 0, delta / h);
   }
   double xrec, yrec;
   mapping_.rectangle(lrefp.perp(), lrefp.z(), xrec, yrec);
 
   Grid1D gridX(xrec, xrec + (a + b) / 2., n1);
   Grid1D gridY(yref, yref + stepy * (n2 - 1), n2);
   Grid1D gridZ(yrec, yrec + h, n3);
   grid_ = GridType(gridX, gridY, gridZ, fieldValues);
 
   // Activate/deactivate timers
   //   static SimpleConfigurable<bool> timerOn(false,"MFGrid:timing");
   //   (*TimingReport::current()).switchOn("MagneticFieldProvider::valueInTesla(TrapezoidalCylindricalMFGrid)",timerOn);
 }

Member Function Documentation

◆ dump()

void TrapezoidalCylindricalMFGrid::dump ( void ) const

overridevirtual

Reimplemented from MFGrid.

Definition at line 129 of file TrapezoidalCylindricalMFGrid.cc.

129 {}

◆ fromGridFrame()

MFGrid::LocalPoint TrapezoidalCylindricalMFGrid::fromGridFrame	(	double	a,
		double	b,
		double	c
	)		const

overridevirtual

find grid coordinates for point. For debugging and validation only.

Implements MFGrid.

Definition at line 149 of file TrapezoidalCylindricalMFGrid.cc.

References a, b, c, mapping_, and Trapezoid2RectangleMappingX::trapezoid().

                                                                                                {
   double rtrap, ztrap;
   mapping_.trapezoid(a, c, rtrap, ztrap);
   // FIXME: "OLD" convention of phi.
   //  return LocalPoint(LocalPoint::Cylindrical(rtrap, Geom::pi() - b, ztrap));
   return LocalPoint(LocalPoint::Cylindrical(rtrap, b, ztrap));
 }

◆ toGridFrame()

void TrapezoidalCylindricalMFGrid::toGridFrame	(	const LocalPoint &	p,
		double &	a,
		double &	b,
		double &	c
	)		const

overridevirtual

find grid coordinates for point. For debugging and validation only.

Implements MFGrid.

Definition at line 142 of file TrapezoidalCylindricalMFGrid.cc.

References a, b, c, mapping_, AlCaHLTBitMon_ParallelJobs::p, and Trapezoid2RectangleMappingX::rectangle().

Referenced by uncheckedValueInTesla().

                                                                                                          {
   mapping_.rectangle(p.perp(), p.z(), a, c);
   // FIXME: "OLD" convention of phi.
   //  b = Geom::pi() - p.phi();
   b = p.phi();
 }

◆ uncheckedValueInTesla()

MFGrid::LocalVector TrapezoidalCylindricalMFGrid::uncheckedValueInTesla ( const LocalPoint & p ) const

overridevirtual

Interpolated field value at given point; does not check for exceptions.

Implements MFGrid3D.

Definition at line 131 of file TrapezoidalCylindricalMFGrid.cc.

References a, b, c, MFGrid::frame(), MFGrid3D::grid_, LinearGridInterpolator3D::interpolate(), AlCaHLTBitMon_ParallelJobs::p, toGridFrame(), and GloballyPositioned< T >::toLocal().

                                                                                                {
   //   static TimingReport::Item & timer= (*TimingReport::current())["MagneticFieldProvider::valueInTesla(TrapezoidalCylindricalMFGrid)"];
   //   TimeMe t(timer,false);
 
   LinearGridInterpolator3D interpol(grid_);
   double a, b, c;
   toGridFrame(p, a, b, c);
   GlobalVector gv(interpol.interpolate(a, b, c));  // grid in global frame
   return frame().toLocal(gv);                      // must return a local vector
 }