CDCAlignment.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <calibration/DBObjCalibrationConstMapBase.h>
#include <cdc/dataobjects/WireID.h>
 
namespace Belle2 {
  // N.B. The original was CDCCalibration.h by Tadeas; copied here and removed
  // calibration-related part. This is tentative for code development; to be
  // moved to an appropriate directory later.
  class CDCAlignment: public DBObjCalibrationConstMapBase {
 
  public:
    // Alignment of wires ------------------------------------------------
    static const baseType wireBwdX = 0;
    static const baseType wireBwdY = 1;
    static const baseType wireBwdZ = 2;
    static const baseType wireFwdX = 4;
    static const baseType wireFwdY = 5;
    static const baseType wireFwdZ = 6;
 
    // Tension and additional wire shape parametrization ---------------
    static const baseType wireTension = 21;
 
    // Alignment of layers ------------------------------------------------
    //Numbering taken from Claus ... no need to change, good idea :-)
    //(1: X(bwd), 2: Y(bwd), 6: Phi(bwd), 11: X(fwd-bwd), 12: Y(fwd-bwd), 16: Phi(fwd-bwd)
 
    static const baseType layerX = 1;
    static const baseType layerY = 2;
    static const baseType layerPhi = 6;
 
    static const baseType layerDx = 11;
    static const baseType layerDy = 12;
    static const baseType layerDPhi = 16;
 
    CDCAlignment() : DBObjCalibrationConstMapBase() {}
 
    double getBwdWireShiftX(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireBwdX);
    }
    double getBwdWireShiftY(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireBwdY);
    }
    double getBwdWireShiftZ(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireBwdZ);
    }
 
    double getFwdWireShiftX(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireFwdX);
    }
    double getFwdWireShiftY(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireFwdY);
    }
    double getFwdWireShiftZ(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireFwdZ);
    }
 
    double getWireTension(int ICLayer,  int ILayer) const
    {
      return get(WireID(ICLayer, ILayer), wireTension);
    }
 
    ~CDCAlignment() {}
 
    // ------------- Interface to global Millepede calibration ----------------
    static unsigned short getGlobalUniqueID() {return 27;}
    double getGlobalParam(unsigned short element, unsigned short param) const
    {
      return get(element, param);
    }
    void setGlobalParam(double value, unsigned short element, unsigned short param)
    {
      set(element, param, value);
    }
    std::vector<std::pair<unsigned short, unsigned short>> listGlobalParams()
    {
      return {};
    }
    // ------------------------------------------------------------------------
 
  private:
 
    ClassDef(CDCAlignment, 2); 
  };
 
} // end namespace Belle2