SurveyAlignment Class Reference

#include <SurveyAlignment.h>

Inheritance diagram for SurveyAlignment:

Public Member Functions
void	iterate (unsigned int nIteration, const std::string &fileName, bool bias=false)
	Run the iteration: find residuals, write to output, shift sensors. More...
	SurveyAlignment (const align::Alignables &sensors, const std::vector< align::StructureType > &levels)
	Constructor to set the sensors and residual levels. More...
virtual	~SurveyAlignment ()

Protected Member Functions
virtual void	findAlignPars (bool bias=false)=0
	Find the alignment parameters for all sensors. More...
virtual void	shiftSensors ()
	Apply the alignment parameters to all sensors. More...

Protected Attributes
const std::vector< align::StructureType > &	theLevels
const align::Alignables &	theSensors

Detailed Description

Alignment using only survey info (no tracks) as a proof of principle.

Date

2007/05/03 20:58:58

Revision

1.4

Author

Chung Khim Lae

Definition at line 16 of file SurveyAlignment.h.

Constructor & Destructor Documentation

SurveyAlignment()

SurveyAlignment::SurveyAlignment	(	const align::Alignables &	sensors,
		const std::vector< align::StructureType > &	levels
	)

Constructor to set the sensors and residual levels.

Definition at line 9 of file SurveyAlignment.cc.

    : theSensors(sensors), theLevels(levels) {}

~SurveyAlignment()

virtual SurveyAlignment::~SurveyAlignment ( )

inlinevirtual

Definition at line 22 of file SurveyAlignment.h.

{}

Member Function Documentation

findAlignPars()

virtual void SurveyAlignment::findAlignPars ( bool  bias = false )

protectedpure virtual

Find the alignment parameters for all sensors.

Implemented in SurveyAlignmentPoints, and SurveyAlignmentSensor.

Referenced by iterate().

iterate()

void SurveyAlignment::iterate	(	unsigned int	nIteration,
		const std::string &	fileName,
		bool	bias = false
	)

Run the iteration: find residuals, write to output, shift sensors.

Definition at line 36 of file SurveyAlignment.cc.

References gather_cfg::cout, AlignmentParameters::covariance(), MillePedeFileConverter_cfg::e, MillePedeFileConverter_cfg::fileName, findAlignPars(), mps_fire::i, dqmiolumiharvest::j, MillePedeFileConverter_cfg::out, AlignmentParameters::parameters(), shiftSensors(), theSensors, and tolerance.

                                                                                           {
  static const double tolerance = 1e-4;  // convergence criteria

  SurveyOutput out(theSensors, fileName);

  out.write(0);

  for (unsigned int i = 1; i <= nIteration; ++i) {
    std::cout << "***** Iteration " << i << " *****\n";
    findAlignPars(bias);
    shiftSensors();
    out.write(i);

    // Check convergence

    double parChi2 = 0.;

    unsigned int nSensor = theSensors.size();

    for (unsigned int j = 0; j < nSensor; ++j) {
      AlignmentParameters* alignPar = theSensors[j]->alignmentParameters();

      const AlgebraicVector& par = alignPar->parameters();
      const AlgebraicSymMatrix& cov = alignPar->covariance();

      int dummy;

      parChi2 += cov.inverse(dummy).similarity(par);
    }

    parChi2 /= static_cast<double>(nSensor);
    std::cout << "chi2 = " << parChi2 << std::endl;
    if (parChi2 < tolerance)
      break;  // converges, so exit loop
  }
}

shiftSensors()

void SurveyAlignment::shiftSensors ( )

protectedvirtual

Apply the alignment parameters to all sensors.

Definition at line 12 of file SurveyAlignment.cc.

References Alignable::alignmentParameters(), mps_fire::i, Alignable::move(), AlignmentParameters::parameters(), align::rectify(), makeMuonMisalignmentScenario::rot, Alignable::rotateInGlobalFrame(), Alignable::surface(), theSensors, AlignableSurface::toGlobal(), and align::toMatrix().

Referenced by iterate().

                                   {
  unsigned int nSensor = theSensors.size();

  for (unsigned int i = 0; i < nSensor; ++i) {
    Alignable* ali = theSensors[i];

    const AlignableSurface& surf = ali->surface();
    const AlgebraicVector& pars = ali->alignmentParameters()->parameters();

    EulerAngles angles(3);

    angles(1) = pars[3];
    angles(2) = pars[4];
    angles(3) = pars[5];

    RotationType rot = surf.toGlobal(toMatrix(angles));

    rectify(rot);  // correct for rounding errors

    ali->move(surf.toGlobal(align::LocalVector(pars[0], pars[1], pars[2])));
    ali->rotateInGlobalFrame(rot);
  }
}

Member Data Documentation

theLevels

const std::vector<align::StructureType>& SurveyAlignment::theLevels

protected

Definition at line 39 of file SurveyAlignment.h.

Referenced by SurveyAlignmentPoints::findAlignPars(), and SurveyAlignmentSensor::findAlignPars().

theSensors

const align::Alignables& SurveyAlignment::theSensors

protected

Definition at line 38 of file SurveyAlignment.h.

Referenced by SurveyAlignmentPoints::findAlignPars(), SurveyAlignmentSensor::findAlignPars(), iterate(), and shiftSensors().