CalibPPS/AlignmentRelative/interface/AlignmentGeometry.h - master - githubmirror/cmssw

CalibPPS/AlignmentRelative/interface/AlignmentGeometry.h
82 строки3 KB
jan-kaspar
PPS: relative alignment (#31629)
06 ноя 2020, 22:59
06 ноя 2020, 22:595a81488
/****************************************************************************
* Authors: 
*  Jan Kašpar (jan.kaspar@gmail.com) 
****************************************************************************/

#ifndef CalibPPS_AlignmentRelative_AlignmentGeometry_h
#define CalibPPS_AlignmentRelative_AlignmentGeometry_h

#include "FWCore/Utilities/interface/Exception.h"

#include <map>
#include <set>
#include <string>

//----------------------------------------------------------------------------------------------------

/**
 *\brief A structure to hold relevant geometrical information about one detector/sensor.
 **/
struct DetGeometry {
  double z;  ///< z postion at detector centre; mm

  double sx, sy;  ///< detector nominal shift = detector center in global coordinates; in mm

  struct DirectionData {
    double dx, dy, dz;  ///< x, y and z components of the direction unit vector in global coordinates
    double s;           ///< projection of (sx, sy) to (dx, dy)
  };

  std::map<unsigned int, DirectionData> directionData;

  bool isU;  ///< only relevant for strips: true for U detectors, false for V detectors
             ///< global U, V frame is used - that matches with u, v frame of the 1200 detector

  DetGeometry(double _z = 0., double _sx = 0., double _sy = 0., bool _isU = false)
      : z(_z), sx(_sx), sy(_sy), isU(_isU) {}

  void setDirection(unsigned int idx, double dx, double dy, double dz) {
    directionData[idx] = {dx, dy, dz, dx * sx + dy * sy};
  }

  const DirectionData& getDirectionData(unsigned int idx) const {
    auto it = directionData.find(idx);
    if (it == directionData.end())
      throw cms::Exception("PPS") << "direction index " << idx << " not in the mapping.";

    return it->second;
  }
};

//----------------------------------------------------------------------------------------------------

/**
 * A collection of geometrical information.
 **/
class AlignmentGeometry {
protected:
  std::map<unsigned int, DetGeometry> sensorGeometry;

public:
  /// a characteristic z in mm
  double z0;

  /// puts an element to the map
  void insert(unsigned int id, const DetGeometry& g);

  /// retrieves sensor geometry
  const DetGeometry& get(unsigned int id) const;

  const std::map<unsigned int, DetGeometry>& getSensorMap() const { return sensorGeometry; }

  /// returns the number of detectors in the collection
  unsigned int getNumberOfDetectors() const { return sensorGeometry.size(); }

  /// check whether the sensor Id is valid (present in the map)
  bool isValidSensorId(unsigned int id) const { return (sensorGeometry.find(id) != sensorGeometry.end()); }

  /// Prints the geometry.
  void print() const;
};

#endif