CompositeAlignmentParameters.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignmentParametrization_CompositeAlignmentParameters_h
#define Alignment_CommonAlignmentParametrization_CompositeAlignmentParameters_h

#include "Alignment/CommonAlignment/interface/AlignableDetOrUnitPtr.h"
#include "Alignment/CommonAlignment/interface/AlignmentParametersData.h"
#include "Alignment/CommonAlignment/interface/Utilities.h"
#include "TrackingTools/TrajectoryState/interface/TrajectoryStateOnSurface.h"

#include <map>
#include <vector>


class AlignableDet;
class Alignable;

class CompositeAlignmentParameters {
public:
  typedef align::Alignables Components;

  typedef std::map<AlignableDetOrUnitPtr, Alignable *> AlignableDetToAlignableMap;
  typedef std::map<Alignable *, int> Aliposmap;
  typedef std::map<Alignable *, int> Alilenmap;

  typedef AlignmentParametersData::DataContainer DataContainer;


  CompositeAlignmentParameters(const AlgebraicVector &par, const AlgebraicSymMatrix &cov, const Components &comp);

  CompositeAlignmentParameters(const AlgebraicVector &par,
                               const AlgebraicSymMatrix &cov,
                               const Components &comp,
                               const AlignableDetToAlignableMap &alimap,
                               const Aliposmap &aliposmap,
                               const Alilenmap &alilenmap);

  CompositeAlignmentParameters(const DataContainer &data,
                               const Components &comp,
                               const AlignableDetToAlignableMap &alimap,
                               const Aliposmap &aliposmap,
                               const Alilenmap &alilenmap);

  virtual ~CompositeAlignmentParameters();

  const AlgebraicVector &parameters() const { return theData->parameters(); }

  const AlgebraicSymMatrix &covariance() const { return theData->covariance(); }

  CompositeAlignmentParameters *clone(const AlgebraicVector &par, const AlgebraicSymMatrix &cov) const;

  CompositeAlignmentParameters *clone(const AlgebraicVector &par,
                                      const AlgebraicSymMatrix &cov,
                                      const AlignableDetToAlignableMap &alimap,
                                      const Aliposmap &aliposmap,
                                      const Alilenmap &alilenmap) const;

  Components components() const;

  AlgebraicMatrix derivatives(const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &alidet) const;
  AlgebraicMatrix selectedDerivatives(const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &alidet) const;
  AlgebraicMatrix derivatives(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                              const std::vector<AlignableDet *> &alidetvec) const;
  AlgebraicMatrix derivatives(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                              const std::vector<AlignableDetOrUnitPtr> &alidetvec) const;
  AlgebraicMatrix selectedDerivatives(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                      const std::vector<AlignableDet *> &alidetvec) const;
  AlgebraicMatrix selectedDerivatives(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                      const std::vector<AlignableDetOrUnitPtr> &alidetvec) const;

  AlgebraicVector correctionTerm(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                 const std::vector<AlignableDet *> &alidetvec) const;
  AlgebraicVector correctionTerm(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                 const std::vector<AlignableDetOrUnitPtr> &alidetvec) const;
  AlgebraicMatrix derivativesLegacy(const TrajectoryStateOnSurface &tsos, AlignableDet *alidet) const;
  AlgebraicMatrix selectedDerivativesLegacy(const TrajectoryStateOnSurface &tsos, AlignableDet *alidet) const;
  AlgebraicMatrix derivativesLegacy(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                    const std::vector<AlignableDet *> &alidetvec) const;
  AlgebraicMatrix selectedDerivativesLegacy(const std::vector<TrajectoryStateOnSurface> &tsosvec,
                                            const std::vector<AlignableDet *> &alidetvec) const;

  Alignable *alignableFromAlignableDet(const AlignableDetOrUnitPtr &adet) const;

  AlgebraicVector parameterSubset(const align::Alignables &) const;

  AlgebraicSymMatrix covarianceSubset(const align::Alignables &) const;

  AlgebraicMatrix covarianceSubset(const align::Alignables &, const align::Alignables &) const;

protected:
  DataContainer theData;

private:
  bool extractPositionAndLength(const align::Alignables &alignables,
                                std::vector<int> &posvec,
                                std::vector<int> &lenvec,
                                int &length) const;

  align::Alignables extractAlignables(const align::Alignables &) const;

  void convert(const std::vector<AlignableDet *> &input, std::vector<AlignableDetOrUnitPtr> &output) const;

  Components theComponents;

  AlignableDetToAlignableMap theAlignableDetToAlignableMap;

  Aliposmap theAliposmap;
  Alilenmap theAlilenmap;
};

#endif