CMSSW Reference Manual

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

Definition at line 10 of file TrapezoidalCylindricalMFGrid.cc.

References a, a1, funct::abs(), b, GenMuonPlsPt100GeV_cfg::cout, Vector3DBase< T, FrameTag >::dot(), lat::endl(), MFGrid::frame(), MFGrid3D::grid_, h, i, j, mapping_, submitDQMOfflineCAF::nLines, PV3DBase< T, PVType, FrameType >::perp(), Geom::pi(), Trapezoid2RectangleMappingX::rectangle(), GloballyPositioned< T >::toLocal(), and PV3DBase< T, PVType, FrameType >::z().

00012   : MFGrid3D(vol)
00013 {
00014   // The parameters read from the data files are given in global coordinates.
00015   // In version 85l, local frame has the same orientation of global frame for the reference
00016   // volume, i.e. the r.f. transformation is only a translation.
00017   // There is therefore no need to convert the field values to local coordinates.
00018   // Check this assumption: 
00019   GlobalVector localXDir(frame().toGlobal(LocalVector(1,0,0)));
00020   GlobalVector localYDir(frame().toGlobal(LocalVector(0,1,0)));
00021 
00022   if (localXDir.dot(GlobalVector(1,0,0)) > 0.999999 &&
00023       localYDir.dot(GlobalVector(0,1,0)) > 0.999999) {
00024     // "null" rotation - requires no conversion...
00025   } else {
00026     cout << "ERROR: TrapezoidalCylindricalMFGrid: unexpected orientation: x: " 
00027          << localXDir << " y: " << localYDir << endl;
00028   }
00029 
00030   int n1, n2, n3;
00031   inFile >> n1 >> n2 >> n3;
00032   double xref, yref, zref;
00033   inFile >> xref >> yref >> zref;
00034   double stepx, stepy, stepz;
00035   inFile >> stepx    >> stepy    >> stepz;
00036 
00037   double BasicDistance1[3][3];  // linear step
00038   double BasicDistance2[3][3];  // linear offset
00039   bool   easya, easyb, easyc;
00040 
00041   inFile >> BasicDistance1[0][0] >> BasicDistance1[1][0] >> BasicDistance1[2][0];
00042   inFile >> BasicDistance1[0][1] >> BasicDistance1[1][1] >> BasicDistance1[2][1];
00043   inFile >> BasicDistance1[0][2] >> BasicDistance1[1][2] >> BasicDistance1[2][2];
00044   inFile >> BasicDistance2[0][0] >> BasicDistance2[1][0] >> BasicDistance2[2][0];
00045   inFile >> BasicDistance2[0][1] >> BasicDistance2[1][1] >> BasicDistance2[2][1];
00046   inFile >> BasicDistance2[0][2] >> BasicDistance2[1][2] >> BasicDistance2[2][2];
00047   inFile >> easya >> easyb >> easyc;
00048 
00049   vector<BVector> fieldValues;
00050   float Bx, By, Bz;
00051   int nLines = n1*n2*n3;
00052   for (int iLine=0; iLine<nLines; ++iLine){
00053     inFile >> Bx >> By >> Bz;
00054     fieldValues.push_back(BVector(Bx,By,Bz));
00055   }
00056   // check completeness
00057   string lastEntry;
00058   inFile >> lastEntry;
00059   if (lastEntry != "complete") {
00060     cout << "ERROR during file reading: file is not complete" << endl;
00061   }
00062 
00063 #ifdef DEBUG_GRID
00064   cout << "easya " << easya << " easyb " << easyb << " easyc " << easyc << endl;
00065 #endif
00066 
00067   if (!easyb || !easyc) {
00068     throw MagGeometryError("TrapezoidalCartesianMFGrid only implemented for first coordinate");
00069   }
00070 
00071 #ifdef DEBUG_GRID
00072   cout << "Grid reference point in grid system: " << xref << "," << yref << "," << zref << endl;
00073   cout << "steps " << stepx << "," <<  stepy << "," << stepz << endl;
00074   cout << "ns " << n1 << "," <<  n2 << "," << n3 << endl;
00075 
00076   for (int i=0; i<3; ++i) for (int j=0; j<3; ++j) {
00077     cout << "BasicDistance1[" << i << "][" << j << "] = " << BasicDistance1[i][j]
00078          << "BasicDistance2[" << i << "][" << j << "] = " << BasicDistance2[i][j] << endl;
00079   }
00080 #endif
00081 
00082   // the "not easy" coordinate is x
00083   double a = stepx * (n1 -1);
00084   double b = a + BasicDistance1[0][1] * (n2-1)*(n1-1) + BasicDistance1[0][2] * (n3-1)*(n1-1);
00085   //  double h = stepy * (n2-1);
00086   double h = stepz * (n3-1);
00087   double delta = -BasicDistance2[0][1] * (n2-1) -BasicDistance2[0][2] * (n3-1);
00088 
00089 #ifdef DEBUG_GRID
00090   cout << "Trapeze size (a,b,h) = " << a << "," << b << "," << h << endl;
00091 #endif
00092 
00093   GlobalPoint grefp( GlobalPoint::Cylindrical( xref, Geom::pi() - yref, zref));
00094   LocalPoint lrefp = frame().toLocal( grefp);
00095 
00096 #ifdef DEBUG_GRID
00097   cout << "Global origin " << grefp << endl;
00098   cout << "Local origin  " << lrefp << endl;
00099 #endif
00100 
00101   double baMinus1 = BasicDistance1[0][2]*(n3-1) / stepx;
00102   if (abs(baMinus1) > 0.000001) {
00103     double b_over_a = 1 + baMinus1;
00104     double a1 = abs(baMinus1) > 0.000001 ? delta / baMinus1 : a/2;
00105 #ifdef DEBUG_GRID
00106    cout << "a1 = " << a1 << endl;
00107 #endif
00108 
00109     // transform reference point to grid frame
00110     double x0 = lrefp.perp() + a1;
00111     double y0 = lrefp.z() + h/2.;
00112     mapping_ = Trapezoid2RectangleMappingX( x0, y0, b_over_a, h);
00113   }
00114   else { // parallelogram
00115     mapping_ = Trapezoid2RectangleMappingX( 0, 0, delta/h);
00116   }
00117   double xrec, yrec;
00118   mapping_.rectangle( lrefp.perp(), lrefp.z(), xrec, yrec);
00119 
00120   Grid1D<double> gridX( xrec, xrec + (a+b)/2., n1);
00121   Grid1D<double> gridY( yref, yref + stepy*(n2-1), n2);
00122   Grid1D<double> gridZ( yrec, yrec + h, n3);
00123   grid_ = GridType( gridX, gridY, gridZ, fieldValues);
00124     
00125   // Activate/deactivate timers
00126 //   static SimpleConfigurable<bool> timerOn(false,"MFGrid:timing");
00127 //   (*TimingReport::current()).switchOn("MagneticFieldProvider::valueInTesla(TrapezoidalCylindricalMFGrid)",timerOn);
00128 }

void TrapezoidalCylindricalMFGrid::dump ( void ) const [virtual]

Reimplemented from MFGrid.

Definition at line 130 of file TrapezoidalCylindricalMFGrid.cc.

00131 {}

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

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

Implements MFGrid.

Definition at line 156 of file TrapezoidalCylindricalMFGrid.cc.

References mapping_, and Trapezoid2RectangleMappingX::trapezoid().

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

void TrapezoidalCylindricalMFGrid::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 146 of file TrapezoidalCylindricalMFGrid.cc.

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

Referenced by uncheckedValueInTesla().

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

MFGrid::LocalVector TrapezoidalCylindricalMFGrid::uncheckedValueInTesla ( const LocalPoint & p ) const [virtual]

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

Implements MFGrid3D.

Definition at line 134 of file TrapezoidalCylindricalMFGrid.cc.

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

00135 {
00136 //   static TimingReport::Item & timer= (*TimingReport::current())["MagneticFieldProvider::valueInTesla(TrapezoidalCylindricalMFGrid)"];
00137 //   TimeMe t(timer,false);
00138 
00139   LinearGridInterpolator3D<GridType::ValueType, GridType::Scalar> interpol( grid_);
00140   double a, b, c;
00141   toGridFrame( p, a, b, c);
00142   GlobalVector gv( interpol.interpolate( a, b, c)); // grid in global frame
00143   return frame().toLocal(gv);           // must return a local vector
00144 }

Trapezoid2RectangleMappingX TrapezoidalCylindricalMFGrid::mapping_ [private]

Definition at line 25 of file TrapezoidalCylindricalMFGrid.h.

Referenced by fromGridFrame(), toGridFrame(), and TrapezoidalCylindricalMFGrid().


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

TrapezoidalCylindricalMFGrid Class Reference

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation