#include <DTSurvey.h>

Public Member Functions
void	CalculateChambers ()

	DTSurvey (const std::string &, const std::string &, int)

const DTSurveyChamber *	getChamber (int, int) const

int	getId () const

void	ReadChambers (edm::ESHandle< DTGeometry >)

	~DTSurvey ()

Private Member Functions
void	FillWheelInfo ()

Private Attributes
DTSurveyChamber ***	chambers

TMatrixD	delta

int	id

std::string	nameOfChamberInfoFile

std::string	nameOfWheelInfoFile

float	OffsetZ

TMatrixD	Rot

Detailed Description

Implements a set of measurements given by survey, tipically a wheel.
Contains DTSurveyChambers and the displacements and rotations for each are calculated.

Date: 2007/12/06 01:53:28

Revision: 1.3

Author: Pablo Martinez Ruiz del Arbol

Definition at line 25 of file DTSurvey.h.

Constructor & Destructor Documentation

◆ DTSurvey()

DTSurvey::DTSurvey	(	const std::string &	Wheel,
		const std::string &	Chambers,
		int	n
	)

Definition at line 12 of file DTSurvey.cc.

                                                                            : chambers(nullptr) {
   nameOfWheelInfoFile = Wheel;
   nameOfChamberInfoFile = Chambers;
   id = n;
  
   FillWheelInfo();
 }

References FillWheelInfo(), dqmiodumpmetadata::n, nameOfChamberInfoFile, and nameOfWheelInfoFile.

◆ ~DTSurvey()

DTSurvey::~DTSurvey ( )

Definition at line 20 of file DTSurvey.cc.

20 { delete[] chambers; }

References chambers.

Member Function Documentation

◆ CalculateChambers()

void DTSurvey::CalculateChambers ( )

Definition at line 22 of file DTSurvey.cc.

                                  {
   for (int stationCounter = 0; stationCounter < 4; stationCounter++) {
     for (int sectorCounter = 0; sectorCounter < 14; sectorCounter++) {
       if (chambers[stationCounter][sectorCounter]->getNumberPoints() > 2) {
         chambers[stationCounter][sectorCounter]->compute();
       }
     }
   }
 }

References chambers, and DTSurveyChamber::compute().

◆ FillWheelInfo()

void DTSurvey::FillWheelInfo ( )

private

Definition at line 120 of file DTSurvey.cc.

                              {
   std::ifstream wheeltowheel(nameOfWheelInfoFile.c_str());
   float zOffset, deltax, deltay, deltaz, alpha, beta, gamma;
   wheeltowheel >> zOffset >> deltax >> deltay >> deltaz >> alpha >> beta >> gamma;
   wheeltowheel.close();
  
   OffsetZ = zOffset;
  
   //Build displacement vector
   delta.ResizeTo(3, 1);
   delta(0, 0) = deltax / 10.0;
   delta(1, 0) = deltay / 10.0;
   delta(2, 0) = deltaz / 10.0;
  
   //Build rotation matrix
   Rot.ResizeTo(3, 3);
   TMatrixD alpha_m(3, 3);
   TMatrixD beta_m(3, 3);
   TMatrixD gamma_m(3, 3);
   alpha_m.Zero();
   beta_m.Zero();
   gamma_m.Zero();
   for (int k = 0; k < 3; k++) {
     alpha_m(k, k) = 1.0;
     beta_m(k, k) = 1.0;
     gamma_m(k, k) = 1.0;
   }
   alpha /= 1000.0;  //New scale: angles in radians
   beta /= 1000.0;
   gamma /= 1000.0;
   alpha_m(1, 1) = cos(alpha);
   alpha_m(1, 2) = sin(alpha);
   alpha_m(2, 1) = -sin(alpha);
   alpha_m(2, 2) = cos(alpha);
   beta_m(0, 0) = cos(beta);
   beta_m(0, 2) = -sin(beta);
   beta_m(2, 0) = sin(beta);
   beta_m(2, 2) = cos(beta);
   gamma_m(0, 0) = cos(gamma);
   gamma_m(0, 1) = sin(gamma);
   gamma_m(1, 0) = -sin(gamma);
   gamma_m(1, 1) = cos(gamma);
   Rot = alpha_m * beta_m * gamma_m;
 }

References zMuMuMuonUserData::alpha, zMuMuMuonUserData::beta, funct::cos(), delta, CustomPhysics_cfi::gamma, dqmdumpme::k, nameOfWheelInfoFile, OffsetZ, Rot, and funct::sin().

Referenced by DTSurvey().

◆ getChamber()

const DTSurveyChamber * DTSurvey::getChamber	(	int	station,
		int	sector
	)		const

Definition at line 32 of file DTSurvey.cc.

32 { return chambers[station][sector]; }

References chambers, and relativeConstraints::station.

◆ getId()

int DTSurvey::getId ( ) const

inline

Definition at line 35 of file DTSurvey.h.

35 { return id; }

References id.

◆ ReadChambers()

void DTSurvey::ReadChambers ( edm::ESHandle< DTGeometry > pDD )

Definition at line 34 of file DTSurvey.cc.

                                                        {
   //Create the chambers
   chambers = new DTSurveyChamber **[4];
   for (int cont_stat = 0; cont_stat < 4; cont_stat++) {
     chambers[cont_stat] = new DTSurveyChamber *[14];
     for (int cont_sect = 0; cont_sect < 14; cont_sect++) {
       DTChamberId mId(id, cont_stat + 1, cont_sect + 1);
       chambers[cont_stat][cont_sect] = new DTSurveyChamber(id, cont_stat + 1, cont_sect + 1, mId.rawId());
     }
   }
  
   std::cout << nameOfChamberInfoFile << std::endl;
   std::ifstream file(nameOfChamberInfoFile.c_str());
   while (!file.eof()) {
     int code, station, sector;
     double x, y, z, rms, dx, dy, dz;
     file >> code >> x >> y >> z >> rms >> dx >> dy >> dz;
     if (file.eof())
       break;
     x = x / 10.0;
     y = y / 10.0;
     z = z / 10.0;
     dx = dx / 10.0;
     dy = dy / 10.0;
     dz = dz / 10.0;
     rms = rms / 10.0;
     station = code / 10000 - 1;
     sector = (code - (station + 1) * 10000) / 100 - 1;
     //De momento vamos a actuar como si no hubiera otra forma de resolver esto
     TMatrixD r(3, 1);
     r(0, 0) = x;
     r(1, 0) = y;
     r(2, 0) = z + OffsetZ;
     TMatrixD disp(3, 1);
     disp(0, 0) = dx;
     disp(1, 0) = dy;
     disp(2, 0) = dz;
     TMatrixD rp = Rot * r - delta;
     disp = disp - r + rp;
  
     GlobalPoint rg(r(0, 0), r(1, 0), r(2, 0));
     GlobalPoint rt(r(0, 0) - disp(0, 0), r(1, 0) - disp(1, 0), r(2, 0) - disp(2, 0));
     DTChamberId mId(id, station + 1, sector + 1);
     const DTChamber *mChamber = static_cast<const DTChamber *>(pDD->idToDet(mId));
     LocalPoint rl = mChamber->toLocal(rg);
     LocalPoint rtl = mChamber->toLocal(rt);
     TMatrixD rLocal(3, 1);
     rLocal(0, 0) = rl.x();
     rLocal(1, 0) = rl.y();
     rLocal(2, 0) = rl.z();
     TMatrixD rTeo(3, 1);
     rTeo(0, 0) = rtl.x();
     rTeo(1, 0) = rtl.y();
     rTeo(2, 0) = rtl.z();
     TMatrixD diff = rLocal - rTeo;
     TMatrixD errors(3, 1);
     errors(0, 0) = rms;
     errors(1, 0) = rms;
     errors(2, 0) = rms;
     chambers[station][sector]->addPoint(code, rLocal, diff, errors);
   }
   file.close();
 }