Ast.h

Go to the documentation of this file.

#ifndef AST_H
#define AST_H

// Ast                                                                //
//                                                                    //
// Package: Ast (AstroUtil).                                          //
//                                                                    //
// Concept:  Define data types and static utility methods used by     //
//           AstroUtil package.                                       //
//                                                                    //
// S. Kasahara  04/2003                                               //
//                                                                    //

#include "Conventions/Detector.h"

namespace AstUtil {

// Values for longitude, latitude, and offset between true north and
// local detector z-axis are from Wes Smart minossoft correspondence 
// of 4/7/2003, in which he states:
// "These positions are given in the datum NAD 83, which uses the reference
//  ellipsoid GRS 80.  The semi-major axis of GRS 80 is 6,378,137 m and
//  the semi-minor axis is 6,356,752.3141 m."
// This data can also be found in the NUMI Technical Design Handbook.
const double kFarDetLongitude = -92.24141202;  // degrees  
const double kFarDetLatitude  =  47.82026653;  // degrees

// Rotation matrix is determined from the measured -26.554823 degree offset 
// of the local far detector z-axis from true north measured clockwise.
// The rotation matrix R[9] from local (O-frame) to ideal (O' frame) 
// coordinates is set up such that:
//  x' = R[0]*x + R[1]*y + R[2]*z  
//  y' = R[3]*x + R[4]*y + R[5]*z
//  z' - R[6]*x + R[7]*y + R[8]*z
// where local and ideal coordinates are defined as:
//  Local Det. Coordinate direction cosines (dcosx,dcosy,dcosz) are 
//  defined in the local coordinate system of a specified detector
//  type (Detector::Detector_t):
//    +z: horizontal component of the beam direction.
//    +y: local vertical.
//    +x: chosen to make coordinate system right-handed.
//
//  Ideal Det. Coordinate direction cosines (dcosx_ideal,dcosy_ideal,
//  dcosz_ideal) are defined such that:
//    +z: points North.
//    +y: radius vector originating at center of earth.
//    +x: points West.   
//

const double kFarDetRotMatrixLocalToIdeal[9] 
             = { 0.894507011,0,0.447053919,
                 0,          1,          0,
                -0.447053919,0,0.894507011};
const double kFarDetRotMatrixIdealToLocal[9] =
 { kFarDetRotMatrixLocalToIdeal[0],
   kFarDetRotMatrixLocalToIdeal[3],
   kFarDetRotMatrixLocalToIdeal[6],
   kFarDetRotMatrixLocalToIdeal[1],
   kFarDetRotMatrixLocalToIdeal[4],
   kFarDetRotMatrixLocalToIdeal[7],
   kFarDetRotMatrixLocalToIdeal[2],
   kFarDetRotMatrixLocalToIdeal[5],
   kFarDetRotMatrixLocalToIdeal[8] };
const double kNearDetLongitude = -88.27062086;  // degrees
const double kNearDetLatitude  =  41.84056333;  // degrees
// Rotation matrix is determined from the measured -23.909201 degree offset 
// of the local near detector z-axis from true north measured clockwise.
const double kNearDetRotMatrixLocalToIdeal[9] 
                          = { 0.914188882,0,0.405288399,
                              0,          1,          0,
                             -0.405288399,0,0.914188882};
const double kNearDetRotMatrixIdealToLocal[9] =
{ kNearDetRotMatrixLocalToIdeal[0],
  kNearDetRotMatrixLocalToIdeal[3],
  kNearDetRotMatrixLocalToIdeal[6],
  kNearDetRotMatrixLocalToIdeal[1],
  kNearDetRotMatrixLocalToIdeal[4],
  kNearDetRotMatrixLocalToIdeal[7],
  kNearDetRotMatrixLocalToIdeal[2],
  kNearDetRotMatrixLocalToIdeal[5],
  kNearDetRotMatrixLocalToIdeal[8] };
                                                     
  double GetDetLongitude(Detector::Detector_t detector);
  double GetDetLatitude(Detector::Detector_t detector);
  const double* GetDetRotMatrixLocalToIdeal(Detector::Detector_t detector);
  const double* GetDetRotMatrixIdealToLocal(Detector::Detector_t detector);

}

#endif  // AST_H

---