/*  G.Kum   5-9-96 */
/*  Detector lookup table for Gammasphere Run GS120 July 1997.
    Out of 110 possible detector positions 103 were filled.

       0   1   2   3   4   5   6   7   8   9   10
     ______________________________________________
  0  |     01, 02, 03, 04, 05, 06, 07, 08, 09, 10,      10
  1  |     11, 12, 13, 14, 15, 16, 17, 18, 19, 20,      10
  2  |     21, 22, 23, 24, 25, 26, 27, 28, 29, 30,      10
  3  |     31, 32, 33, 34, 35, 36, 37, 38, 39, 40,      10
  4  |     41, 42, 43, 44, 45, 00, 47, 48, 49, 00,       8
  5  |     51, 52, 00, 54, 55, 56, 57, 00, 59, 60,       8
  6  |     61, 62, 63, 64, 65, 00, 67, 68, 69, 70,       9
  7  |     71, 72, 00, 74, 75, 76, 77, 00, 79, 80,       8
  8  |     81, 82, 83, 84, 85, 86, 87, 88, 89, 90,      10
  9  |     91, 92, 93, 94, 95, 96, 97, 98, 99,100,      10
 10  |    101,102,103,104,105,106,107,108,109,110       10
                                                       ___
                                                       103
Check timing of det 77
*/
/* convention for this run */
#define UP   32			/* (2) on flip box */
#define DOWN 16			/* (1) on flip box */


#define NDET 110		/* full array */
#define NGAM 103		/* actual number */

static int detarray[] = {0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,
			   11,12,13,14,15,16,17,18,19,20,
			   21,22,23,24,25,26,27,28,29,30,
                           31,32,33,34,35,36,37,38,39,40,
			   41,42,43,44,45, 0,46,47,48, 0,
			   49,50, 0,51,52,53,54, 0,55,56,
			   57,58,59,60,61, 0,62,63,64,65,
			   66,67, 0,68,69,70,71, 0,72,73,
			   74,75,76,77,78,79,80,81,82,83,
			   84,85,86,87,88,89,90,91,92,93,
			   94,95,96,97,98,99,100,101,102,103
                        };

/* used in output file */

#define NSP  103

static int detid[NSP] =  { 1, 2, 3, 4, 5, 6, 7, 8, 9,10,
			  11,12,13,14,15,16,17,18,19,20,
			  21,22,23,24,25,26,27,28,29,30,
			  31,32,33,34,35,36,37,38,39,40,
			  41,42,43,44,45,   47,48,49,  
			  51,52,   54,55,56,57,   59,60,
			  61,62,63,64,65,   67,68,69,70,
			  71,72,   74,75,76,77,   79,80,
			  81,82,83,84,85,86,87,88,89,90,
			  91,92,93,94,95,96,97,98,99,100,
			  101,102,103,104,105,106,107,108,109,110};

static float angle[] =
/*  0     1      2      3      4      5      6      7      8      9     10   */
{  0.0, 17.27, 17.27, 17.27, 17.27, 31.72, 17.27, 31.72, 31.72, 31.72, 31.72,
/* 1*/  37.38, 37.38, 37.38, 37.38, 50.07, 37.38, 50.07, 50.07, 50.07, 50.07,
/* 2*/  50.07, 50.07, 50.07, 50.07, 58.28, 50.07, 58.28, 58.28, 69.82, 58.28,
/* 3*/  69.82, 58.28, 69.82, 69.82, 69.82, 69.82, 69.82, 69.82, 79.19, 69.82,
/* 4*/  79.19, 69.82, 80.71, 79.19, 80.71, 79.19, 80.71, 79.19, 90.00, 80.71,
/* 5*/  90.00, 80.71, 90.00, 90.00, 90.00, 90.00, 90.00, 90.00, 99.29, 90.00,
/* 6*/  99.29, 90.00,100.81, 99.29,100.81, 99.29,100.81, 99.29,110.18,100.81,
/* 7*/ 110.18,100.81,110.18,110.18,110.18,110.18,110.18,110.18,121.72,110.18,
/* 8*/ 121.72,110.18,121.72,121.72,129.93,121.72,129.93,129.93,129.93,129.93,
/* 9*/ 129.93,129.93,129.93,129.93,142.62,129.93,142.62,142.62,142.62,142.62,
/*10*/ 148.28,148.28,148.28,148.28,162.73,148.28,162.73,162.73,162.73,162.73
};


int ring[17][10] = {  {  1, 2, 3, 4, 6 },                 /*  17.27 */
		   {  5, 7, 8, 9,10 },			  /*  31.72 */
		   { 11,12,13,14,16 },			  /*  37.38 */
		   { 15,17,18,19,20,21,22,23,24,26 },	  /*  50.07 */
		   { 25,27,28,30,32 },			  /*  58.28 */
		   { 29,31,33,34,35,36,37,38,40,42 },     /*  69.82 */
		   { 39,41,44,   48 },			  /*  79.19 */
		   { 43,45,47,  52 },			  /*  80.71 */
		   { 49,51,   54,55,56,57,   60,62 },	  /*  90.00 */
		   { 59,61,64,   68 },			  /*  99.29 */
		   { 63,65,67,70,72 },			  /* 100.81 */
		   { 69,71,   74,75,76,77,   80,82 },	  /* 110.18 */
		   { 79,81,83,84,86 },			  /* 121.72 */
		   { 85,87,88,89,90,91,92,93,94,96 },	  /* 129.93 */
		   { 95,97,98,99,100 },			  /* 142.62 */
		   {101,102,103,104,106},		  /* 148.28 */
		   {105,107,108,109,110} };		  /* 162.73 */

/*
   complete angle set  gs(alfa,beta) in Gamma Sphere notation.
              This are not polar coordinates !
   Transformation to polar cordinates as used in angular correlations
   are :
         tan(theta) = tan(alfa) * cos(gamma);
	 sin(phi)   = sin(alfa) * sin(gamma);

	 beta   0 -  90 ==> gamma =  90 - beta;
	       90 - 180 ==>       = beta -  90;
	      180 - 270 ==>       = 270 - beta;
	      270 - 360 ==>       = beta - 270;
   with phi symmetric above and below a horizontal reaction plane at 90 deg.
   to the magnetic field direction.
 */

static double gs_angles[111][2] = {
/*   0 */ {  0.00,  0.00},
/*   1 */ { 17.27, 72.00},
/*   2 */ { 17.27,216.00},
/*   3 */ { 17.27,144.00},
/*   4 */ { 17.27,288.00},
/*   5 */ { 31.72, 36.00},
/*   6 */ { 17.27,360.00},
/*   7 */ { 31.72,108.00},
/*   8 */ { 31.72,252.00},
/*   9 */ { 31.72,180.00},
/*  10 */ { 31.72,324.00},
/*  11 */ { 37.38, 72.00},
/*  12 */ { 37.38,216.00},
/*  13 */ { 37.38,144.00},
/*  14 */ { 37.38,288.00},
/*  15 */ { 50.07, 22.84},
/*  16 */ { 37.38,360.00},
/*  17 */ { 50.07, 49.16},
/*  18 */ { 50.07,193.16},
/*  19 */ { 50.07, 94.84},
/*  20 */ { 50.07,238.84},
/*  21 */ { 50.07,121.16},
/*  22 */ { 50.07,265.16},
/*  23 */ { 50.07,166.84},
/*  24 */ { 50.07,310.84},
/*  25 */ { 58.28, 72.00},
/*  26 */ { 50.07,337.16},
/*  27 */ { 58.28,144.00},
/*  28 */ { 58.28,216.00},
/*  29 */ { 69.82, 18.49},
/*  30 */ { 58.28,288.00},
/*  31 */ { 69.82, 53.51},
/*  32 */ { 58.28,360.00},
/*  33 */ { 69.82, 90.49},
/*  34 */ { 69.82,197.51},
/*  35 */ { 69.82,125.51},
/*  36 */ { 69.82,234.49},
/*  37 */ { 69.82,162.49},
/*  38 */ { 69.82,269.51},
/*  39 */ { 79.19, 72.00},
/*  40 */ { 69.82,306.49},
/*  41 */ { 79.19,144.00},
/*  42 */ { 69.82,341.51},
/*  43 */ { 80.71, 36.00},
/*  44 */ { 79.19,216.00},
/*  45 */ { 80.71,108.00},
/*  46 */ { 79.19,288.00},
/*  47 */ { 80.71,180.00},
/*  48 */ { 79.19,360.00},
/*  49 */ { 90.00, 18.00},
/*  50 */ { 80.71,252.00},
/*  51 */ { 90.00, 54.00},
/*  52 */ { 80.71,324.00},
/*  53 */ { 90.00, 90.00},
/*  54 */ { 90.00,198.00},
/*  55 */ { 90.00,126.00},
/*  56 */ { 90.00,234.00}, 
/*  57 */ { 90.00,162.00},
/*  58 */ { 90.00,270.00},
/*  59 */ { 99.29, 72.00},
/*  60 */ { 90.00,306.00},
/*  61 */ { 99.29,144.00},
/*  62 */ { 90.00,342.00},
/*  63 */ {100.81, 36.00},
/*  64 */ { 99.29,216.00},
/*  65 */ {100.81,108.00},
/*  66 */ { 99.29,288.00},
/*  67 */ {100.81,180.00},
/*  68 */ { 99.29,360.00},
/*  69 */ {110.18, 17.51},
/*  70 */ {100.81,252.00},
/*  71 */ {110.18, 54.49},
/*  72 */ {100.81,324.00},
/*  73 */ {110.18, 89.51},
/*  74 */ {110.18,198.49},
/*  75 */ {110.18,126.49},
/*  76 */ {110.18,233.51},
/*  77 */ {110.18,161.51},
/*  78 */ {110.18,270.49},
/*  79 */ {121.72, 36.00},
/*  80 */ {110.18,305.51},
/*  81 */ {121.72,108.00},
/*  82 */ {110.18,342.49},
/*  83 */ {121.72,180.00},
/*  84 */ {121.72,252.00},
/*  85 */ {129.93, 13.16},
/*  86 */ {121.72,324.00},
/*  87 */ {129.93, 58.84},
/*  88 */ {129.93,202.84},
/*  89 */ {129.93, 85.16},
/*  90 */ {129.93,229.16},
/*  91 */ {129.93,130.84},
/*  92 */ {129.93,274.84},
/*  93 */ {129.93,157.16},
/*  94 */ {129.93,301.16},
/*  95 */ {142.62, 36.00},
/*  96 */ {129.93,346.84},
/*  97 */ {142.62,108.00},
/*  98 */ {142.62,252.00},
/*  99 */ {142.62,180.00},
/* 100 */ {142.62,324.00},
/* 101 */ {148.28, 72.00},
/* 102 */ {148.28,216.00},
/* 103 */ {148.28,144.00},
/* 104 */ {148.28,288.00},
/* 105 */ {162.73, 36.00},
/* 106 */ {148.28,360.00},
/* 107 */ {162.73,108.00},
/* 108 */ {162.73,252.00},
/* 109 */ {162.73,180.00}, 
/* 110 */ {162.73,324.00}
          };




