BmnMeasureModel.h

Go to the documentation of this file.

/*
    BM@N alignment routine
    BM@N experiment at NICA complex, JINR, 2025
 
    Department: Math & Soft Group of HEP lab
    Author:     Igor Polev, polev@jinr.ru
 
    BmnMeasureModel class header-only implementation
    Parametrized model of linear track
*/
 
#ifndef BMNMEASUREMODEL_H
#define BMNMEASUREMODEL_H
 
#include "BmnAlignDefines.h"
#include "FairHit.h"
 
/*
 * This class is intended to be the only place where the sign selection and rotation
 * sequence scheme of alignment is defined. Namely, the scheme is as follows:
 *
 *      | x_m |   | x |   | s_x |
 *      | y_m | = | y | + | s_y |
 *      | z_0 |   | z |   | s_z |
 *
 * where:
 *      x, y, z       : true values of particle hit coordinates in lab frame of reference
 *      x_m, y_m      : measured values of x and y
 *      z_0           : nominal z coordinate of detector plane
 *      s_x, s_y, s_z : shift placement errors of detector plane
 *      R_zyx         : combined rotation matrix of detector plane in sequence
 *                      over z, over y, over x, composed of rotation placement errors
 */
 
class BmnMeasureModel
{
  public:
    BmnMeasureModel() = default;
    inline Double_t X(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline Double_t Y(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline SVectLC GradXlocal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline SVectLC GradYlocal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline SVectGL GradXglobal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline SVectGL GradYglobal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const;
    inline void AlignHit(FairHit& hit, SVectGL& alignValues) const;
};
 
inline Double_t BmnMeasureModel::X(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return lc(iX0) + lc(iX1) * (z0 - gl(iZ)) + gl(iX);
}
 
inline Double_t BmnMeasureModel::Y(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return lc(iY0) + lc(iY1) * (z0 - gl(iZ)) + gl(iY);
}
 
inline SVectLC BmnMeasureModel::GradXlocal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return SVectLC{1.0, z0 - gl(iZ), 0.0, 0.0};
}
 
inline SVectLC BmnMeasureModel::GradYlocal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return SVectLC{0.0, 0.0, 1.0, z0 - gl(iZ)};
}
 
inline SVectGL BmnMeasureModel::GradXglobal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return SVectGL{1.0, 0.0, -lc(iX1)};
}
 
inline SVectGL BmnMeasureModel::GradYglobal(Double_t z0, const SVectGL& gl, const SVectLC& lc) const
{
    return SVectGL{0.0, 1.0, -lc(iY1)};
}
 
inline void BmnMeasureModel::AlignHit(FairHit& hit, SVectGL& alignValues) const
{
    hit.SetX(hit.GetX() - alignValues(iX));
    hit.SetY(hit.GetY() - alignValues(iY));
    hit.SetZ(hit.GetZ() - alignValues(iZ));
}
 
#endif   // BMNMEASUREMODEL_H