DTSurveyChamber.h

Go to the documentation of this file.

 
#ifndef Alignment_SurveyAnalysis_DTSurveyChamber_H
#define Alignment_SurveyAnalysis_DTSurveyChamber_H
 
#include <vector>
 
#include "TMath.h"
#include "TMatrixD.h"
 
class DTSurveyChamber {
public:
  DTSurveyChamber(int, int, int, long);
 
  //Add a point to the Chamber
  void addPoint(int, const TMatrixD &, const TMatrixD &, const TMatrixD &);
 
  //Begin the chi2 computation
  int getNumberPoints() const { return pointNumber; }
  void compute();
  void printChamberInfo();
  long getId() const { return rawId; }
  float getDeltaX() const { return Solution(0, 0); }
  float getDeltaY() const { return Solution(1, 0); }
  float getDeltaZ() const { return Solution(2, 0); }
  //My definition of the matrix rotation is not the same as the used in CMSSW
  float getAlpha() const { return Solution(5, 0); }
  float getBeta() const { return -1.0 * Solution(4, 0); }
  float getGamma() const { return Solution(3, 0); }
  float getDeltaXError() const { return TMath::Sqrt(Covariance(0, 0)); }
  float getDeltaYError() const { return TMath::Sqrt(Covariance(1, 1)); }
  float getDeltaZError() const { return TMath::Sqrt(Covariance(2, 2)); }
  float getAlphaError() const { return TMath::Sqrt(Covariance(5, 5)); }
  float getBetaError() const { return TMath::Sqrt(Covariance(4, 4)); }
  float getGammaError() const { return TMath::Sqrt(Covariance(3, 3)); }
 
private:
  //   static const unsigned int MAX_PUNTOS = 16;
 
  TMatrixD &makeMatrix();
  TMatrixD &makeErrors();
  TMatrixD &makeVector();
 
  //Identifiers
  int wheel, station, sector;
 
  long rawId;
 
  //Points data
  std::vector<TMatrixD> points;
  std::vector<TMatrixD> pointsDiff;
  std::vector<TMatrixD> pointsError;
  std::vector<TMatrixD> pointsTheoretical;
 
  TMatrixD Solution;
  TMatrixD Covariance;
 
  //Number of points
  int pointNumber;
};
 
std::ostream &operator<<(std::ostream &, const DTSurveyChamber &);
 
#endif