TwoBowedSurfacesAlignmentParameters.h

Go to the documentation of this file.

#ifndef Alignment_CommonAlignment_TwoBowedSurfacesAlignmentParameters_h
#define Alignment_CommonAlignment_TwoBowedSurfacesAlignmentParameters_h

#include "Alignment/CommonAlignment/interface/AlignmentParameters.h"
#include "CondFormats/Alignment/interface/Definitions.h"
#include "DataFormats/CLHEP/interface/AlgebraicObjects.h"

#include "Alignment/CommonAlignmentParametrization/interface/BowedSurfaceAlignmentDerivatives.h"


class Alignable;
class AlignableDetOrUnitPtr;
class TrajectoryStateOnSurface;

class TwoBowedSurfacesAlignmentParameters : public AlignmentParameters {
public:
  typedef BowedSurfaceAlignmentDerivatives BowedDerivs;
  enum AlignmentParameterName {
    // 1st surface
    dx1 = BowedDerivs::dx,
    dy1 = BowedDerivs::dy,
    dz1 = BowedDerivs::dz,
    dslopeX1 = BowedDerivs::dslopeX,  // NOTE: slope(u) -> halfWidth*tan(beta),
    dslopeY1 = BowedDerivs::dslopeY,  //       slope(v) -> halfLength*tan(alpha)
    drotZ1 = BowedDerivs::drotZ,      //       rot(w)   -> g-scale*gamma
    dsagittaX1 = BowedDerivs::dsagittaX,
    dsagittaXY1 = BowedDerivs::dsagittaXY,
    dsagittaY1 = BowedDerivs::dsagittaY,
    // 2nd surface
    dx2 = BowedDerivs::dx + BowedDerivs::N_PARAM,
    dy2 = BowedDerivs::dy + BowedDerivs::N_PARAM,
    dz2 = BowedDerivs::dz + BowedDerivs::N_PARAM,
    dslopeX2 = BowedDerivs::dslopeX + BowedDerivs::N_PARAM,  // NOTE: slope(u) -> k*tan(beta),
    dslopeY2 = BowedDerivs::dslopeY + BowedDerivs::N_PARAM,  //       slope(v) -> k*tan(alpha)
    drotZ2 = BowedDerivs::drotZ + BowedDerivs::N_PARAM,      //       rot(w)   -> m*gamma
    dsagittaX2 = BowedDerivs::dsagittaX + BowedDerivs::N_PARAM,
    dsagittaXY2 = BowedDerivs::dsagittaXY + BowedDerivs::N_PARAM,
    dsagittaY2 = BowedDerivs::dsagittaY + BowedDerivs::N_PARAM,
    // number of parameters
    N_PARAM = BowedDerivs::N_PARAM + BowedDerivs::N_PARAM
  };

  TwoBowedSurfacesAlignmentParameters(Alignable *alignable);

  TwoBowedSurfacesAlignmentParameters(Alignable *alignable,
                                      const AlgebraicVector &parameters,
                                      const AlgebraicSymMatrix &covMatrix);

  TwoBowedSurfacesAlignmentParameters(Alignable *alignable,
                                      const AlgebraicVector &parameters,
                                      const AlgebraicSymMatrix &covMatrix,
                                      const std::vector<bool> &selection);

  ~TwoBowedSurfacesAlignmentParameters() override{};
  void apply() override;
  int type() const override;

  TwoBowedSurfacesAlignmentParameters *clone(const AlgebraicVector &parameters,
                                             const AlgebraicSymMatrix &covMatrix) const override;

  TwoBowedSurfacesAlignmentParameters *cloneFromSelected(const AlgebraicVector &parameters,
                                                         const AlgebraicSymMatrix &covMatrix) const override;

  AlgebraicMatrix derivatives(const TrajectoryStateOnSurface &tsos, const AlignableDetOrUnitPtr &aliDet) const override;

  virtual void print() const;

  double ySplit() const { return ySplit_; }

private:
  double ySplitFromAlignable(const Alignable *ali) const;

  double ySplit_;
};

#endif