---

RectangularCylindricalMFGrid Class Reference

#include <MagneticField/Interpolation/src/RectangularCylindricalMFGrid.h>

Inheritance diagram for RectangularCylindricalMFGrid:

List of all members.

Public Member Functions
virtual void	dump () const
virtual LocalPoint	fromGridFrame (double a, double b, double c) const
	find grid coordinates for point. For debugging and validation only.
	RectangularCylindricalMFGrid (binary_ifstream &istr, const GloballyPositioned< float > &vol)
virtual void	toGridFrame (const LocalPoint &p, double &a, double &b, double &c) const
	find grid coordinates for point. For debugging and validation only.
virtual LocalVector	uncheckedValueInTesla (const LocalPoint &p) const
	Interpolated field value at given point; does not check for exceptions.

---

Detailed Description

Definition at line 8 of file RectangularCylindricalMFGrid.h.

---

Constructor & Destructor Documentation

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

Definition at line 8 of file RectangularCylindricalMFGrid.cc.

References GenMuonPlsPt100GeV_cfg::cout, Vector3DBase< T, FrameTag >::dot(), lat::endl(), MFGrid::frame(), MFGrid3D::grid_, submitDQMOfflineCAF::nLines, PV3DBase< T, PVType, FrameType >::perp(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::z().

  : 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: RectangularCylindricalMFGrid: 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;

  vector<BVector> fieldValues;
  float Bx, By, Bz;
  int nLines = n1*n2*n3;
  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;
  }

  GlobalPoint grefp( GlobalPoint::Cylindrical( xref, yref, zref));
  LocalPoint lrefp = frame().toLocal( grefp);

#ifdef DEBUG_GRID
  cout << "Grid reference point in grid system: " << xref << "," << yref << "," << zref << endl;
  cout << "Grid reference point in global x,y,z: " << grefp << endl;
  cout << "Grid reference point in local x,y,z: " << lrefp << endl;
  cout << "steps " << stepx << "," <<  stepy << "," << stepz << endl;
#endif

  Grid1D<double> gridX( lrefp.perp(), lrefp.perp() + stepx*(n1-1), n1);
  //Grid1D<double> gridY( lrefp.phi(), lrefp.phi() + stepy*(n2-1), n2); // wrong: gives zero
  Grid1D<double> gridY( yref, yref + stepy*(n2-1), n2);
  Grid1D<double> gridZ( lrefp.z(), lrefp.z() + stepz*(n3-1), n3);

  grid_ = GridType( gridX, gridY, gridZ, fieldValues);
  
}

---

Member Function Documentation

void RectangularCylindricalMFGrid::dump

(

void

)

const [virtual]

Reimplemented from MFGrid.

Definition at line 68 of file RectangularCylindricalMFGrid.cc.

References GenMuonPlsPt100GeV_cfg::cout, Grid3D< Value, T >::data(), lat::endl(), MFGrid3D::grid_, Grid3D< Value, T >::grida(), Grid3D< Value, T >::gridb(), and Grid3D< Value, T >::gridc().

{
  cout << endl << "Dump of RectangularCylindricalMFGrid" << endl;
  cout << "Number of points from Grid1D " 
       << grid_.grida().nodes() << " " << grid_.gridb().nodes() << " " << grid_.gridc().nodes() << endl;

  cout << "Reference Point from Grid1D " 
       << grid_.grida().lower() << " " << grid_.gridb().lower() << " " << grid_.gridc().lower() << endl;

  cout << "Basic Distance from Grid1D "
       << grid_.grida().step() << " " << grid_.gridb().step() << " " << grid_.gridc().step() << endl;


  cout << "Dumping " << grid_.data().size() << " field values " << endl;
  // grid_.dump();
}

MFGrid::LocalPoint RectangularCylindricalMFGrid::fromGridFrame	(	double	a,
		double	b,
		double	c
	)			const [virtual]

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

Implements MFGrid.

Definition at line 115 of file RectangularCylindricalMFGrid.cc.

{

  // FIXME: "OLD" convention of phi.
  //  return LocalPoint( LocalPoint::Cylindrical(a, Geom::pi() - b, c));
  return LocalPoint( LocalPoint::Cylindrical(a, b, c));
}

void RectangularCylindricalMFGrid::toGridFrame	(	const LocalPoint &	p,
		double &	a,
		double &	b,
		double &	c
	)			const [virtual]

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

Implements MFGrid.

Definition at line 105 of file RectangularCylindricalMFGrid.cc.

References PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), and PV3DBase< T, PVType, FrameType >::z().

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

MFGrid::LocalVector RectangularCylindricalMFGrid::uncheckedValueInTesla

(

const LocalPoint &

p

)

const [virtual]

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

Implements MFGrid3D.

Definition at line 85 of file RectangularCylindricalMFGrid.cc.

References e, MFGrid3D::grid_, Grid3D< Value, T >::grida(), Grid3D< Value, T >::gridc(), LinearGridInterpolator3D< Value, T >::interpolate(), k, PV3DBase< T, PVType, FrameType >::perp(), PV3DBase< T, PVType, FrameType >::phi(), dttmaxenums::R, HLT_VtxMuL3::result, value, and PV3DBase< T, PVType, FrameType >::z().

{
  const float minimalSignificantR = 1e-6; // [cm], points below this radius are treated as zero radius
  float R = p.perp();
  if (R < minimalSignificantR) {
    if (grid_.grida().lower() < minimalSignificantR) {
      int k = grid_.gridc().index(p.z());
      double u = (p.z() - grid_.gridc().node(k)) / grid_.gridc().step();
      LocalVector result((1-u)*grid_(0,  0,  k) + u*grid_(0,  0,  k+1)); 
      return result;
    }
  }
  
  LinearGridInterpolator3D<GridType::ValueType, GridType::Scalar> interpol( grid_);
  // FIXME: "OLD" convention of phi.
  // GridType::ValueType value = interpol( R, Geom::pi() - p.phi(), p.z());
  GridType::ValueType value = interpol.interpolate( R, p.phi(), p.z());
  return LocalVector(value);
}

---

The documentation for this class was generated from the following files:

• MagneticField/Interpolation/src/RectangularCylindricalMFGrid.h
• MagneticField/Interpolation/src/RectangularCylindricalMFGrid.cc

---