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 "DataFormats/CLHEP/interface/AlgebraicObjects.h"
#include "CondFormats/Alignment/interface/Definitions.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