TOPutil.h

//*****************************************************************************
//-----------------------------------------------------------------------------
// Fast simulation and reconstruction package for TOP counter (F77 core)
// M. Staric, March-2009
//-----------------------------------------------------------------------------
//
// TOPutil.h
// C++ interface to F77 functions: data conversion utility (for debugging etc.)
// units: cm, ns
// mostly opsolete now !
//-----------------------------------------------------------------------------
//*****************************************************************************
 
#ifndef _TOPutil_h
#define _TOPutil_h
 
#include <top/reconstruction/TOPf77fun.h>
 
namespace Belle2 {
  namespace TOP {
 
    inline void PointToLocal(int moduleID, double X, double Y, double Z,
                             double& x, double& y, double& z)
    {
      moduleID--;
      float r[3] = {(float) X, (float) Y, (float) Z};
      float dir[3] = {0, 0, 0};
      to_local_(r, dir, &moduleID);
      x = r[0]; y = r[1]; z = r[2];
    }
 
    inline void PointToGlobal(int moduleID, double X, double Y, double Z,
                              double& x, double& y, double& z)
    {
      moduleID--;
      float r[3] = {(float) X, (float) Y, (float) Z};
      float dir[3] = {0, 0, 0};
      to_global_(r, dir, &moduleID);
      x = r[0]; y = r[1]; z = r[2];
    }
 
    inline void MomentumToLocal(int moduleID, double Px, double Py, double Pz,
                                double& px, double& py, double& pz)
    {
      moduleID--;
      float r[3] = {0, 0, 0};
      float dir[3] = {(float) Px, (float) Py, (float) Pz};
      to_local_(r, dir, &moduleID);
      px = dir[0]; py = dir[1]; pz = dir[2];
    }
 
    inline void MomentumToGlobal(int moduleID, double Px, double Py, double Pz,
                                 double& px, double& py, double& pz)
    {
      moduleID--;
      float r[3] = {0, 0, 0};
      float dir[3] = {(float) Px, (float) Py, (float) Pz};
      to_global_(r, dir, &moduleID);
      px = dir[0]; py = dir[1]; pz = dir[2];
    }
 
    inline bool InsideQbar(double X, double Y, double Z, int moduleID)
    {
      moduleID--;
      float x = (float) X; float y = (float) Y; float z = (float) Z;
      return inside_qbar_(&x, &y, &z, &moduleID) != 0;
    }
 
    inline bool InsideExpansionVolume(double X, double Y, double Z, int moduleID)
    {
      moduleID--;
      float x = (float) X; float y = (float) Y; float z = (float) Z;
      return inside_extvol_(&x, &y, &z, &moduleID) != 0;
    }
 
    inline double PhaseIndex(double E)
    {
      float e = (float) E;
      return phase_index_(&e);
    }
 
    inline double PhaseIndex()
    {
      float refind, gruind;
      get_top_par_(&refind, &gruind);
      return refind;
    }
 
    inline double GroupIndex(double E)
    {
      float e = (float) E;
      return group_index_(&e);
    }
 
    inline double GroupIndex()
    {
      float refind, gruind;
      get_top_par_(&refind, &gruind);
      return gruind;
    }
 
    inline double AtenuLength(double E)
    {
      float e = (float) E;
      return abspt_leng_(&e);
    }
 
    inline double MirrorRefl(double E)
    {
      float e = (float) E;
      return mirror_refl_(&e);
    }
 
    inline double QEffi(double E)
    {
      float e = (float) E;
      return qe_pmt_(&e);
    }
 
    inline void QErange(double& Emin, double& Emax)
    {
      float emin, emax;
      qe_range_(&emin, &emax);
      Emin = emin; Emax = emax;
    }
 
    inline double QEmax()
    {
      return qe_max_();
    }
 
    inline double Nzero()
    {
      return qe_nzero_();
    }
 
    inline double CerAngle(double P, double Mass, double RefrIndex)
    {
      float p = (float) P; float mass = (float) Mass;
      float refind = (float) RefrIndex;
      return cher_angle_(&p, &mass, &refind);
    }
 
 
  } // top namespace
} // Belle2 namespace
 
 
#endif