ARICHBtestGeometryPar.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.                  *
 **************************************************************************/
 
#ifndef ARICHBTESTGEOMETRYPAR_H
#define ARICHBTESTGEOMETRYPAR_H
 
#include <utility>
#include <vector>
#include <framework/gearbox/GearDir.h>
#include <Math/Vector3D.h>
#include <Math/Rotation3D.h>
 
namespace Belle2 {
 
  class ARICHTracking {
  public:
    ARICHTracking():
      atdc(0)
    {
      for (int i = 0; i < 4; i++) tdc[i] = 0;
      for (int i = 0; i < 2; i++) {
        slp[i] = 0;
        offset[i] = 0;
        cutll[i] = 0;
        cutul[i] = 0;
        diff[i] = 0;
        sum[i] = 0;
        status[i] = 0;
      }
      for (int i = 0; i < 3; i++) {
        pos[i] = 0;
        reco[i] = 0;
      }
    };
    ~ARICHTracking() {};
    int tdc[4];
    int atdc;
    float slp[2];
    float offset[2];
    int cutll[2];
    int cutul[2];
    float pos[3];
 
    int diff[2];
    int sum[2];
    int status[2];
    double reco[3];
    void Print()
    {
      printf("----------------------------------------------------------------\n");
      for (int i = 0; i < 4; i++) printf("tdc%d= %d\t", i, tdc[i]);
      printf("\n");
      for (int i = 0; i < 1; i++) printf("atdc%d= %d\t", i, atdc);
      printf("\n");
      for (int i = 0; i < 2; i++) printf("slp%d= %3.3f\t", i, slp[i]);
      printf("\n");
      for (int i = 0; i < 2; i++) printf("off%d= %3.3f\t", i, offset[i]);
      printf("\n");
      for (int i = 0; i < 2; i++) printf("cutll%d= %d\t", i, cutll[i]);
      printf("\n");
      for (int i = 0; i < 2; i++) printf("cutul%d= %d\t", i, cutul[i]);
      printf("\n");
      for (int i = 0; i < 3; i++) printf("pos%d= %3.3f\t", i, pos[i]);
      printf("\n");
    };
  };
 
 
  class ARICHBtestGeometryPar {
 
  public:
 
    ARICHBtestGeometryPar();
 
    ARICHBtestGeometryPar(const ARICHBtestGeometryPar& arichPar);
 
    virtual ~ARICHBtestGeometryPar();
 
    ARICHBtestGeometryPar& operator=(const ARICHBtestGeometryPar&) = delete;
 
 
    static ARICHBtestGeometryPar* Instance();
    void Initialize(const GearDir& content);
    void Initialize(const GearDir& modulecontent, const GearDir& mirrorcontent);
 
 
    void clear(void);
 
    void Print(void) const;
 
    void read(const GearDir& content);
 
 
    void setTrackingShift(const ROOT::Math::XYZVector&);
    ROOT::Math::XYZVector getTrackingShift();
    void setRotationCenter(const ROOT::Math::XYZVector&);
    ROOT::Math::XYZVector getRotationCenter();
 
    void setFrameRotation(double);
    ROOT::Math::Rotation3D getFrameRotation();
    void setAverageAgel(bool);
    bool getAverageAgel();
    void setMwpc(ARICHTracking* m_mwpc);
    ARICHTracking* getMwpc();
    std::pair<double, double> GetHapdChannelPosition(int);
    std::pair<int, int> GetHapdElectronicMap(int);
    int AddHapdChannelPositionPair(double, double);
    int AddHapdElectronicMapPair(int, int);
    void setOffset(const ROOT::Math::XYZVector&);
    ROOT::Math::XYZVector getOffset();
 
 
  private:
    ARICHTracking* m_mwpc;
    ROOT::Math::XYZVector   m_offset;
    ROOT::Math::XYZVector   m_trackingShift;
    ROOT::Math::XYZVector   m_rotationCenter;
    ROOT::Math::Rotation3D  m_frameRotation;
    bool           m_averageAgel;
    std::vector < std::pair<double, double> > m_hapdmap;
    std::vector < std::pair<int, int> > m_hapdeid;
 
    int m_nPads;                           
    bool m_init;                          
    bool m_simple;                        
    static ARICHBtestGeometryPar* p_B4ARICHBtestGeometryParDB; 
    // vectors holding information on HAPDs and chips and pads positions.
 
 
  };
 
  //-----------------------------------------------------------------------------
 
} // end of namespace Belle2
 
#endif