CMS 3D CMS Logo

All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
List of all members | Public Types | Public Member Functions
RectangularCartesianMFGrid Class Reference

#include <RectangularCartesianMFGrid.h>

Inheritance diagram for RectangularCartesianMFGrid:
MFGrid3D MFGrid MagneticFieldProvider< float >

Public Types

using binary_ifstream = magneticfield::interpolation::binary_ifstream
 
- 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, GlobalTagGlobalPointType
 
typedef Vector3DBase< float, GlobalTagGlobalVectorType
 
typedef Point3DBase< float, LocalTagLocalPointType
 
typedef Vector3DBase< float, LocalTagLocalVectorType
 

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...
 
 RectangularCartesianMFGrid (binary_ifstream &istr, const GloballyPositioned< float > &vol)
 
void toGridFrame (const LocalPoint &p, double &a, double &b, double &c) const override
 find grid coordinates for point. For debugging and validation only. More...
 
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 ()
 

Additional Inherited Members

- 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 11 of file RectangularCartesianMFGrid.h.

Member Typedef Documentation

◆ binary_ifstream

Definition at line 13 of file RectangularCartesianMFGrid.h.

Constructor & Destructor Documentation

◆ RectangularCartesianMFGrid()

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

Definition at line 9 of file RectangularCartesianMFGrid.cc.

References gather_cfg::cout, MFGrid::frame(), MFGrid3D::grid_, testHGCalDigi_cfg::inFile, submitDQMOfflineCAF::nLines, and GloballyPositioned< T >::toLocal().

10  : MFGrid3D(vol) {
11  // The parameters read from the data files are given in global coordinates.
12  // In version 85l, local frame has the same orientation of global frame for the reference
13  // volume, i.e. the r.f. transformation is only a translation.
14  // There is therefore no need to convert the field values to local coordinates.
15  // Check this assumption:
16  GlobalVector localXDir(frame().toGlobal(LocalVector(1, 0, 0)));
17  GlobalVector localYDir(frame().toGlobal(LocalVector(0, 1, 0)));
18 
19  if (localXDir.dot(GlobalVector(1, 0, 0)) > 0.999999 && localYDir.dot(GlobalVector(0, 1, 0)) > 0.999999) {
20  // "null" rotation - requires no conversion...
21  } else {
22  cout << "ERROR: RectangularCartesianMFGrid: unexpected orientation: x: " << localXDir << " y: " << localYDir
23  << endl;
24  }
25 
26  int n1, n2, n3;
27  inFile >> n1 >> n2 >> n3;
28  double xref, yref, zref;
29  inFile >> xref >> yref >> zref;
30  double stepx, stepy, stepz;
31  inFile >> stepx >> stepy >> stepz;
32 
33  vector<BVector> fieldValues;
34  float Bx, By, Bz;
35  int nLines = n1 * n2 * n3;
36  fieldValues.reserve(nLines);
37  for (int iLine = 0; iLine < nLines; ++iLine) {
38  inFile >> Bx >> By >> Bz;
39  fieldValues.push_back(BVector(Bx, By, Bz));
40  }
41  // check completeness
42  string lastEntry;
43  inFile >> lastEntry;
44  if (lastEntry != "complete") {
45  cout << "ERROR during file reading: file is not complete" << endl;
46  }
47 
48  GlobalPoint grefp(xref, yref, zref);
49  LocalPoint lrefp = frame().toLocal(grefp);
50 
51  Grid1D gridX(lrefp.x(), lrefp.x() + stepx * (n1 - 1), n1);
52  Grid1D gridY(lrefp.y(), lrefp.y() + stepy * (n2 - 1), n2);
53  Grid1D gridZ(lrefp.z(), lrefp.z() + stepz * (n3 - 1), n3);
54  grid_ = GridType(gridX, gridY, gridZ, fieldValues);
55 
56  // Activate/deactivate timers
57  // static SimpleConfigurable<bool> timerOn(false,"MFGrid:timing");
58  // (*TimingReport::current()).switchOn("MagneticFieldProvider::valueInTesla(RectangularCartesianMFGrid)",timerOn);
59 }
MFGrid3D(const GloballyPositioned< float > &vol)
Definition: MFGrid3D.h:19
T z() const
Definition: PV3DBase.h:61
GloballyPositioned< float >::GlobalVector GlobalVector
Definition: MFGrid.h:30
LocalPoint toLocal(const GlobalPoint &gp) const
T x() const
Definition: PV3DBase.h:59
T y() const
Definition: PV3DBase.h:60
const GloballyPositioned< float > & frame() const
Local reference frame.
Definition: MFGrid.h:60
GridType grid_
Definition: MFGrid3D.h:59
Definition: Grid1D.h:7
Grid3D GridType
Definition: MFGrid3D.h:56
Grid3D::BVector BVector
Definition: MFGrid3D.h:57
GloballyPositioned< float >::LocalVector LocalVector
Definition: MFGrid.h:32

Member Function Documentation

◆ dump()

void RectangularCartesianMFGrid::dump ( void  ) const
overridevirtual

Reimplemented from MFGrid.

Definition at line 61 of file RectangularCartesianMFGrid.cc.

References gather_cfg::cout, Grid3D::data(), MFGrid3D::grid_, Grid3D::grida(), Grid3D::gridb(), Grid3D::gridc(), Grid1D::lower(), Grid1D::nodes(), and Grid1D::step().

61  {
62  cout << endl << "Dump of RectangularCartesianMFGrid" << endl;
63  // cout << "Number of points from file "
64  // << n1 << " " << n2 << " " << n3 << endl;
65  cout << "Number of points from Grid1D " << grid_.grida().nodes() << " " << grid_.gridb().nodes() << " "
66  << grid_.gridc().nodes() << endl;
67 
68  // cout << "Reference Point from file "
69  // << xref << " " << yref << " " << zref << endl;
70  cout << "Reference Point from Grid1D " << grid_.grida().lower() << " " << grid_.gridb().lower() << " "
71  << grid_.gridc().lower() << endl;
72 
73  // cout << "Basic Distance from file "
74  // << stepx << " " << stepy << " " << stepz << endl;
75  cout << "Basic Distance from Grid1D " << grid_.grida().step() << " " << grid_.gridb().step() << " "
76  << grid_.gridc().step() << endl;
77 
78  cout << "Dumping " << grid_.data().size() << " field values " << endl;
79  // grid_.dump();
80 }
const Container & data() const
Definition: Grid3D.h:63
const Grid1D & grida() const
Definition: Grid3D.h:59
Scalar lower() const
Definition: Grid1D.h:19
Scalar step() const
Definition: Grid1D.h:18
int nodes() const
Definition: Grid1D.h:21
GridType grid_
Definition: MFGrid3D.h:59
const Grid1D & gridb() const
Definition: Grid3D.h:60
const Grid1D & gridc() const
Definition: Grid3D.h:61

◆ fromGridFrame()

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

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

Implements MFGrid.

Definition at line 97 of file RectangularCartesianMFGrid.cc.

References a, b, and HltBtagPostValidation_cff::c.

97  {
98  return LocalPoint(a, b, c);
99 }
GloballyPositioned< float >::LocalPoint LocalPoint
Definition: MFGrid.h:31
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121

◆ toGridFrame()

void RectangularCartesianMFGrid::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 91 of file RectangularCartesianMFGrid.cc.

References a, b, HltBtagPostValidation_cff::c, and AlCaHLTBitMon_ParallelJobs::p.

91  {
92  a = p.x();
93  b = p.y();
94  c = p.z();
95 }
double b
Definition: hdecay.h:120
double a
Definition: hdecay.h:121

◆ uncheckedValueInTesla()

MFGrid::LocalVector RectangularCartesianMFGrid::uncheckedValueInTesla ( const LocalPoint p) const
overridevirtual

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

Implements MFGrid3D.

Definition at line 82 of file RectangularCartesianMFGrid.cc.

References MFGrid3D::grid_, LinearGridInterpolator3D::interpolate(), and AlCaHLTBitMon_ParallelJobs::p.

82  {
83  // static TimingReport::Item & timer= (*TimingReport::current())["MagneticFieldProvider::valueInTesla(RectangularCartesianMFGrid)"];
84  // TimeMe t(timer,false);
85 
87  GridType::ReturnType value = interpol.interpolate(p.x(), p.y(), p.z());
88  return LocalVector(value);
89 }
Definition: value.py:1
ValueType ReturnType
Definition: Grid3D.h:33
GridType grid_
Definition: MFGrid3D.h:59
GloballyPositioned< float >::LocalVector LocalVector
Definition: MFGrid.h:32