FitGeoParams.h

Go to the documentation of this file.

/*
 * FitGeoParams.h
 *
 *  Created on: 20 Feb 2020
 *      Author: Valerio Formato
 */
 
#ifndef HERD_FITGEOPARAMS_H_
#define HERD_FITGEOPARAMS_H_
 
#ifdef HS_USE_ROOT
#include "Rtypes.h"
#endif
 
// C++ headers
#include <optional>
#include <sstream>
#include <vector>
 
// HERD headers
#include "herd/analysis/common/AxesArray.h"
#include "herd/analysis/common/SidesArray.h"
#include "herd/analysis/dataobjects/Parallelepiped.h"
#include "herd/analysis/dataobjects/Point.h"
#include "herd/analysis/dataobjects/TrackingDetectorGeoParams.h"
 
namespace Herd {
 
// clang-format off
  // sipm setup parameters
 
  //  Default Mat and SiPMs
  //  -----------------------------------------------------------
  //  |---------------|    |---------------|    |---------------|
  //  |               |    |               |    |               |
  //  |    SiPM0      |0.09|    SiPM1      |0.09|    SiPM2      |
  //  |   32.54mm     | mm |   32.54mm     | mm |   32.54mm     |
  //  |               |    |               |    |               |
  //  |---------------|    |---------------|    |---------------|
  //  -----------------------------------------------------------
  //  \_________________________97.8mm__________________________/
  //
  //
  // SiPM structure
  //
  //  ==================================================
  //  |    |                |    |                 |    |
  //  |    |                |    |                 |    |
  //  |    |  active die0   |    |  active die0    |    |
  //  |0.16|                |0.22|                 |0.16|
  //  | mm |  0.25mm*64     | mm |  0.25mm*64      | mm |
  //  |    |                |    |                 |    |
  //  |    |                |    |                 |    |
  //  |    |                |    |                 |    |
  //  ==================================================
  //  \__________________32.54mm_______________________/
  //
// clang-format on
 
class FitGeoParams : public TrackingDetectorGeoParams {
public:
  FitGeoParams();
 
  FitGeoParams(RefFrame::Side dside, size_t layerNumber, RefFrame::Direction mainSegm);
 
  ~FitGeoParams() override = default;
 
  const Parallelepiped &BoundingBox() const;
 
  [[nodiscard]] unsigned int NSipmPerMat() const { return _nSipmsPerMat; }
 
  [[nodiscard]] int NChannelsPerSipm() const { return _sipmDieChannels * _sipmNDies; }
 
  [[nodiscard]] int NChannelsPerMat() const { return _nSipmsPerMat * NChannelsPerSipm(); }
 
  [[nodiscard]] const AxesArray<int> &NMats() const { return _nMats; }
 
  [[nodiscard]] unsigned int NTotElements() const {
    return _nMats[RefFrame::Axis::X] * _nMats[RefFrame::Axis::Y] * _nMats[RefFrame::Axis::Z];
  };
 
  [[nodiscard]] const AxesArray<float> &MatSize() const { return _matSize; }
 
  [[nodiscard]] float SipmAirGap() const { return _sipmAirGap; };
 
  [[nodiscard]] Point ElementPosition(unsigned int elementID) const override;
 
  [[nodiscard]] Parallelepiped ElementShape(size_t volID) const;
 
  [[nodiscard]] Point SipmPosition(unsigned int sipmID) const;
 
  [[nodiscard]] Point MatPosition(unsigned int matID) const { return _matPos.at(matID); };
 
  [[nodiscard]] std::optional<unsigned int> MatID(const Point &position) const;
 
  [[nodiscard]] unsigned int MatIDFromChannelID(unsigned int channelID) const;
 
  [[nodiscard]] unsigned int SipmIDFromChannelID(unsigned int channelID) const;
 
  [[nodiscard]] std::optional<unsigned int> ElementID(const Point &position) const override;
 
  [[nodiscard]] std::optional<unsigned int> SipmIDOnMat(float relativePosition) const;
 
  [[nodiscard]] std::optional<unsigned int> ChannelIDOnMat(float relativePosition) const;
 
  [[nodiscard]] float ChannelPositionOnMat(unsigned int channelIDOnMat) const;
 
  void SetSideAndMainSegmDir(RefFrame::Side dside, RefFrame::Direction segmDir) override;
 
  void SetNSipmPerMat(int nSipms) { _nSipmsPerMat = nSipms; };
 
  void SetMatSize(const AxesArray<float> &sizes);
 
  void SetNMats(AxesArray<int> nMats);
 
  void SetMatPosition(size_t matID, Point position);
 
  class LayerZone {
  public:
    LayerZone(float minSegm, float maxSegm) : _min{minSegm}, _max{maxSegm} {}
 
    [[nodiscard]] float MinSegm() const { return _min; }
 
    [[nodiscard]] float MaxSegm() const { return _max; }
 
  private:
    float _min, _max;
  };
 
  using LayerZones = std::vector<LayerZone>;
 
  const LayerZones &ActiveZones();
 
  const LayerZones &DeadZones();
 
private:
  unsigned int _nSipmsPerMat;          
  AxesArray<int> _nMats;               
  AxesArray<float> _matSize;           
  std::vector<Point> _matPos;          
  bool mutable _bbIsValid{false};      
  Parallelepiped mutable _boundingBox; 
 
  LayerZones _activeZones, _deadZones; 
 
  void ComputeBoundingBox() const;
  void ComputeActiveDeadZones();
 
protected:
  float _sipmAirGap = 0.009;       
  float _sipmDieGap = 0.022;       
  float _sipmDeadRegion = 0.016;   
  float _sipmChannelPitch = 0.025; 
  size_t _sipmDieChannels = 64;    
  size_t _sipmNDies = 2;           
  float _sipmWidth =
      2 * _sipmDieChannels * _sipmChannelPitch + 2 * _sipmDeadRegion + _sipmDieGap; 
  std::vector<std::pair<float, float>> _sipmRelativePositions; 
 
#ifdef HS_USE_ROOT
  ClassDefOverride(FitGeoParams, 2)
#endif
};
 
} // namespace Herd
 
#endif /* HERD_FITGEOPARAMS_H_ */