AlignmentParameterStore.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignmentAlgorithm_AlignmentParameterStore_h
#define Alignment_CommonAlignmentAlgorithm_AlignmentParameterStore_h

#include "Alignment/CommonAlignment/interface/Utilities.h"
#include "Alignment/CommonAlignmentAlgorithm/interface/AlignableData.h"
#include "Alignment/CommonAlignmentAlgorithm/interface/AlignmentCorrelationsStore.h"
#include "Alignment/CommonAlignmentParametrization/interface/CompositeAlignmentParameters.h"


namespace edm { class ParameterSet; }
class AlignmentUserVariables;

class AlignmentParameterStore 
{

public:

  typedef std::vector<AlignmentParameters*> Parameters;
  typedef std::map< std::pair<Alignable*,Alignable*>,AlgebraicMatrix > Correlations;
  typedef std::vector<unsigned int> DetIds;

  AlignmentParameterStore( const align::Alignables &alis, const edm::ParameterSet& config );

  virtual ~AlignmentParameterStore();

  CompositeAlignmentParameters
    selectParameters( const std::vector <AlignableDet*>& alignabledets ) const;
  CompositeAlignmentParameters
    selectParameters( const std::vector <AlignableDetOrUnitPtr>& alignabledets ) const;
  CompositeAlignmentParameters
    selectParameters( const std::vector <Alignable*>& alignables ) const;

  void updateParameters(const CompositeAlignmentParameters& aap, bool updateCorrelations = true);

  const align::Alignables& alignables(void) const { return theAlignables; }

  align::Alignables validAlignables(void) const;

  int numObjects(void) const { return theAlignables.size(); }

  AlignmentCorrelationsStore* correlationsStore( void ) const { return theCorrelationsStore; }

  const unsigned int numCorrelations( void ) const { return theCorrelationsStore->size(); }

/*   Alignable* alignableFromGeomDet( const GeomDet* geomDet ) const; */

  Alignable* alignableFromAlignableDet( AlignableDetOrUnitPtr alignableDet ) const;

/*   Alignable* alignableFromDetId(const unsigned int& detId) const; */

  void applyParameters(void);

  void applyParameters(Alignable* alignable);

  void resetParameters(void);

  void resetParameters(Alignable* ali);

  void cacheTransformations(void);

  void restoreCachedTransformations(void);

  void acquireRelativeParameters(void);

  void applyAlignableAbsolutePositions( const align::Alignables& alis, 
                                        const AlignablePositions& newpos, int& ierr );

  void applyAlignableRelativePositions( const align::Alignables& alivec, 
                                        const AlignableShifts& shifts, int& ierr );

  void attachAlignmentParameters( const align::Alignables& alivec, const Parameters& parvec, int& ierr );

  void attachAlignmentParameters(const Parameters& parvec, int& ierr);

  void attachCorrelations( const align::Alignables& alivec, const Correlations& cormap, 
                           bool overwrite, int& ierr );

  void attachCorrelations( const Correlations& cormap, bool overwrite, int& ierr );

  void attachUserVariables( const align::Alignables& alivec,
                            const std::vector<AlignmentUserVariables*>& uvarvec, int& ierr);

  void setAlignmentPositionError( const align::Alignables& alivec, double valshift, double valrot );

  std::pair<int,int> typeAndLayer( const Alignable* ali ) const;

  typedef std::pair<Alignable*, unsigned int> ParameterId;
  bool hierarchyConstraints(const Alignable *aliMaster, const align::Alignables &aliComps,
                std::vector<std::vector<ParameterId> > &paramIdsVecOut,
                std::vector<std::vector<double> > &factorsVecOut,
                bool all, double epsilon) const;

protected:

  // storage for correlations
  AlignmentCorrelationsStore* theCorrelationsStore;

private:
  enum TypeOfConstraints { NONE, HIERARCHY_CONSTRAINTS, APPROX_AVERAGING_CONSTRAINTS };

  // data members

  align::Alignables theAlignables;

  TypeOfConstraints theTypeOfConstraints;  
};

#endif