CompositeAlignmentParameters.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignmentParametrization_CompositeAlignmentParameters_h
#define Alignment_CommonAlignmentParametrization_CompositeAlignmentParameters_h


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

#include <map>
#include <vector>


class AlignableDet;
class Alignable;

class CompositeAlignmentParameters
{

public:

  typedef std::vector<Alignable*> 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 std::vector<Alignable*>& vec ) const;

  AlgebraicSymMatrix covarianceSubset( const std::vector<Alignable*>& vec ) const;

  AlgebraicMatrix covarianceSubset ( const std::vector<Alignable*>& veci,
                                     const std::vector<Alignable*>& vecj ) const;

protected:
  DataContainer theData;

private:

  bool extractPositionAndLength( const std::vector<Alignable*>& alignables,
                 std::vector<int>& posvec,
                 std::vector<int>& lenvec,
                 int& length ) const;

  std::vector< Alignable* > extractAlignables( const std::vector< Alignable* >& alignables ) const;

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

  Components theComponents;

  AlignableDetToAlignableMap theAlignableDetToAlignableMap;

  Aliposmap theAliposmap;
  Alilenmap theAlilenmap;

};

#endif