/* G.Kum 5-9-96 */ /* Detector lookup table for Gammasphere Run GS55 May 1996. Out of 110 possible detector positions 90 where filled. 0 1 2 3 4 5 6 7 8 9 10 ______________________________________________ 0 | 00, 00, 00, 00, 05, 00, 07, 08, 09, 10, 5 1 | 11, 12, 13, 14, 15, 00, 17, 18, 19, 20, 9 2 | 21, 00, 23, 24, 25, 26, 27, 28, 29, 00, 8 3 | 31, 32, 00, 34, 35, 36, 37, 00, 39, 42, 8 4 | 41, 40, 43, 00, 45, 46, 47, 48, 49, 50, 9 5 | 51, 52, 00, 54, 55, 56, 57, 00, 59, 60, 8 6 | 61, 62, 63, 64, 65, 66, 00, 68, 69, 00, 8 7 | 71, 72, 00, 00, 75, 00, 77, 00, 79, 00, 5 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 ___ 90 The electronic channel of det 40 and 42 are interchanged! */ #define UP 4 /* convention for this run */ #define DOWN 2 #define NDET 110 /* full array */ #define NGAM 90 /* actual number */ static int detarray[] = {0, 0, 0, 0, 0, 1, 0, 2, 3, 4, 5, 6, 7, 8, 9,10, 0,11,12,13,14, 15, 0,16,17,18,19,20,21,22, 0, 23,24, 0,25,26,27,28, 0,29,30, 31,32,33, 0,34,35,36,37,38,39, 40,41, 0,42,43,44,45, 0,46,47, 48,49,50,51,52,53, 0,54,55, 0, 56,57, 0, 0,58, 0,59, 0,60, 0, 61,62,63,64,65,66,67,68,69,70, 71,72,73,74,75,76,77,78,79,80, 81,82,83,84,85,86,87,88,89,90 }; 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, 37.38, 37.38, 37.38, 37.38, 50.07, 37.38, 50.07, 50.07, 50.07, 50.07, 50.07, 50.07, 50.07, 50.07, 58.28, 50.07, 58.28, 58.28, 69.82, 58.28, 69.82, 58.28, 69.82, 69.82, 69.82, 69.82, 69.82, 69.82, 79.19, 69.82, 79.19, 69.82, 80.71, 79.19, 80.71, 79.19, 80.71, 79.19, 90.00, 80.71, 90.00, 80.71, 90.00, 90.00, 90.00, 90.00, 90.00, 90.00, 99.29, 90.00, 99.29, 90.00,100.81, 99.29,100.81, 99.29,100.81, 99.29,110.18,100.81, 110.18,100.81,110.18,110.18,110.18,110.18,110.18,110.18,121.72,110.18, 121.72,110.18,121.72,121.72,129.93,121.72,129.93,129.93,129.93,129.93, 129.93,129.93,129.93,129.93,142.62,129.93,142.62,142.62,142.62,142.62, 148.28,148.28,148.28,148.28,162.73,148.28,162.73,162.73,162.73,162.73 }; int ring[17][10] = { { 0 }, /* 17.27 */ { 5, 7, 8, 9,10 }, /* 31.17 */ { 11,12,13,14, }, /* 37.38 */ { 15,17,18,19,20,21 ,23,24,26 }, /* 50.07 */ { 25,27,28, 32 }, /* 58.28 */ { 29,31, 34,35,36,37, 40,42 }, /* 69.82 */ { 39,41, 46,48 }, /* 79.19 */ { 43,45,47,50,52 }, /* 80.71 */ { 49,51, 54,55,56,57, 60,62 }, /* 90.00 */ { 59,61,64,66,68 }, /* 99.29 */ { 63,65, 72 }, /* 100.81 */ { 69,71, 75, 77, 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 */ /* used in output file */ /* here 40 and 42 are interchanged */ #define NSP 90 static int detid[] = { 5, 7, 8, 9,10,11,12,13,14,15,17,18,19,20,21, 23,24,25,26,27,28,29,31,32,34,35,36,37,39,42, 41,40,43,45,46,47,48,49,50,51,52,54,55,56,57, 59,60,61,62,63,64,65,66,68,69,71,72,75,77,79, 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 }; /* 4.18.97 gkum ********************* 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}, /* 42 */ { 69.82,341.51}, /* channels 40 and 42 are interchanged */ /* 41 */ { 79.19,144.00}, /* for this experiment gs55 */ /* 40 */ { 69.82,306.49}, /* 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} }; /* polar coordinates of GS detectors corr_angles[theta][phi] the additional tilt of the magnetic field dir. (72.0 deg) is applied in the calculation, calculated with cvangles.c */ static double corr_angles[111][2] = { /* 0 */ 0.00, 0.00, /* 1 */ 0.00, 17.27, /* 2 */ 16.47, 5.26, /* 3 */ 10.36, 13.90, /* 4 */ 10.36, 13.90, /* 5 */ 19.97, 25.17, /* 6 */ 16.47, 5.26, /* 7 */ 19.97, 25.17, /* 8 */ 0.00, 31.72, /* 9 */ 30.45, 9.35, /* 10 */ 30.45, 9.35, /* 11 */ 0.00, 37.38, /* 12 */ 24.18, 29.42, /* 13 */ 36.00, 10.81, /* 14 */ 24.18, 29.42, /* 15 */ 42.11, 30.10, /* 16 */ 36.00, 10.81, /* 17 */ 24.88, 44.97, /* 18 */ 45.63, 23.38, /* 19 */ 24.88, 44.97, /* 20 */ 15.22, 48.30, /* 21 */ 42.11, 30.10, /* 22 */ 15.22, 48.30, /* 23 */ 49.97, 3.71, /* 24 */ 45.63, 23.38, /* 25 */ 0.00, 58.28, /* 26 */ 49.97, 3.71, /* 27 */ 56.98, 15.24, /* 28 */ 43.56, 43.49, /* 29 */ 65.43, 33.93, /* 30 */ 43.56, 43.49, /* 31 */ 40.79, 62.89, /* 32 */ 56.98, 15.24, /* 33 */ 40.79, 62.89, /* 34 */ 65.70, 33.04, /* 35 */ 65.43, 33.93, /* 36 */ 39.30, 63.52, /* 37 */ 69.82, 0.46, /* 38 */ 39.30, 63.52, /* 39 */ 0.00, 79.19, /* 40 */ 69.82, 0.46, /* 41 */ 78.65, 17.67, /* 42 */ 65.70, 33.04, /* 43 */ 74.45, 52.98, /* 44 */ 72.00, 52.62, /* 45 */ 74.45, 52.98, /* 46 */ 72.00, 52.62, /* 47 */ 80.24, 17.76, /* 48 */ 78.65, 17.67, /* 49 */ 90.00, 36.00, /* 50 */ 0.00, 80.71, /* 51 */ 90.00, 72.00, /* 52 */ 80.24, 17.76, /* 53 */ 90.00, 72.00, /* 54 */ 90.00, 36.00, /* 55 */ 90.00, 36.00, /* 56 */ 90.00, 72.00, /* 57 */ 90.00, 0.00, /* 58 */ 90.00, 72.00, /* 59 */ 180.00, 80.71, /* 60 */ 90.00, 36.00, /* 61 */ 99.76, 17.76, /* 62 */ 90.00, 0.00, /* 63 */ 108.00, 52.62, /* 64 */ 105.55, 52.98, /* 65 */ 108.00, 52.62, /* 66 */ 105.55, 52.98, /* 67 */ 101.35, 17.67, /* 68 */ 99.76, 17.76, /* 69 */ 114.30, 33.04, /* 70 */ 180.00, 79.19, /* 71 */ 140.70, 63.52, /* 72 */ 101.35, 17.67, /* 73 */ 140.70, 63.52, /* 74 */ 114.57, 33.93, /* 75 */ 114.30, 33.04, /* 76 */ 139.21, 62.89, /* 77 */ 110.18, 0.46, /* 78 */ 139.21, 62.89, /* 79 */ 136.44, 43.49, /* 80 */ 114.57, 33.93, /* 81 */ 136.44, 43.49, /* 82 */ 110.18, 0.46, /* 83 */ 123.02, 15.24, /* 84 */ 180.00, 58.28, /* 85 */ 134.37, 23.38, /* 86 */ 123.02, 15.24, /* 87 */ 164.78, 48.30, /* 88 */ 137.89, 30.10, /* 89 */ 164.78, 48.30, /* 90 */ 155.12, 44.97, /* 91 */ 134.37, 23.38, /* 92 */ 155.12, 44.97, /* 93 */ 130.03, 3.71, /* 94 */ 137.89, 30.10, /* 95 */ 155.82, 29.42, /* 96 */ 130.03, 3.71, /* 97 */ 155.82, 29.42, /* 98 */ 180.00, 37.38, /* 99 */ 144.00, 10.81, /* 100 */ 144.00, 10.81, /* 101 */ 180.00, 31.72, /* 102 */ 160.03, 25.17, /* 103 */ 149.55, 9.35, /* 104 */ 160.03, 25.17, /* 105 */ 169.64, 13.90, /* 106 */ 149.55, 9.35, /* 107 */ 169.64, 13.90, /* 108 */ 180.00, 17.27, /* 109 */ 163.53, 5.26, /* 110 */ 163.53, 5.26 }; /********************************************/ /* efficiency parameters for GammaSphere run 55 (May 1996)*/ /* indirect calibration using the target activity only */ /* relative efficiencies for rings 1 through 17 appended * as det 111 to 127 */ /* This parameters were derived with 'automatic' peak integration */ static double par[128][7] = { /* 0 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 1 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 2 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 3 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 4 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 5 chisq = 1.112 */ {5.750213, 13.210927, 0.000000, 6.790561, -1.806645, 0.264147, 0.804853 }, /* 6 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 7 chisq = 1.131 */ {3.314149, 4.415628, 0.000000, 5.088307, -0.655896, -0.172305, 10.510428 }, /* 8 chisq = 1.247 */ {3.635080, 8.268676, 0.000000, 6.800647, -1.836076, 0.727715, 0.936814 }, /* 9 chisq = 1.566 */ {2.079820, 4.987703, 0.000000, 5.121924, -0.664920, -0.273115, 8.060673 }, /* 10 chisq = 1.073 */ {2.418481, 4.729702, 0.000000, 5.083944, -0.650897, -0.166816, 5.149927 }, /* 11 chisq = 0.813 */ {3.632053, 10.373451, 0.000000, 5.239198, -1.436889, 0.020016, 0.946525 }, /* 12 chisq = 1.472 */ {2.581248, 4.548288, 0.000000, 4.916094, -0.501576, -0.011272, 5.268396 }, /* 13 chisq = 0.587 */ {2.332819, 4.872175, -0.367954, 5.159807, -0.601321, -0.147797, 8.182960 }, /* 14 chisq = 3.636 */ {2.346027, 4.634368, 0.000000, 5.035557, -0.527838, -0.079310, 6.818114 }, /* 15 chisq = 0.758 */ {2.382142, 4.656944, -1.166124, 5.144148, -0.550933, -0.198440, 14.865014 }, /* 16 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 17 chisq = 2.764 */ {5.903457, 15.693186, 0.000000, 7.846776, -2.203360, 0.210173, 0.655773 }, /* 18 chisq = 1.076 */ {4.408181, 7.253144, -1.874275, 5.862153, -0.713655, 0.153600, 1.802935 }, /* 19 chisq = 2.403 */ {2.362440, 5.669462, 0.000000, 5.042142, -0.601178, 0.012656, 3.692889 }, /* 20 chisq = 0.696 */ {3.048125, 5.141279, -1.533949, 5.791918, 0.100344, 0.604589, 3.169075 }, /* 21 chisq = 3.213 */ {1.114166, 7.290479, -1.533949, 5.355533, -0.689803, -0.458130, 2.808466 }, /* 22 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 23 chisq = 0.562 */ {2.405696, 6.098998, 0.000000, 5.269316, -0.661988, 0.206124, 2.298205 }, /* 24 chisq = 1.231 */ {3.511189, 8.539850, -1.227277, 5.370478, -0.917837, -0.132290, 1.742103 }, /* 25 chisq = 0.245 */ {1.959560, 6.543041, -1.148486, 4.898231, -0.490627, -0.173397, 7.743435 }, /* 26 chisq = 1.233 */ {2.850857, 4.560819, 0.006141, 4.980702, -0.696470, -0.215576, 11.251479 }, /* 27 chisq = 1.605 */ {1.933759, 5.307440, -0.160849, 5.056849, -0.667606, -0.207250, 6.251122 }, /* 28 chisq = 0.728 */ {2.750837, 4.824002, -1.228662, 4.932102, -0.483000, -0.130142, 13.581328 }, /* 29 chisq = 1.073 */ {3.234175, 9.391302, 0.000000, 6.265216, -1.369051, 0.335062, 1.030776 }, /* 30 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 31 chisq = 2.079 */ {1.831869, 5.546296, -1.603369, 4.621699, -0.350337, -0.026620, 6.865518 }, /* 32 chisq = 1.012 */ {2.972854, 5.939398, 0.000000, 5.063381, -0.580557, 0.030610, 3.458166 }, /* 33 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 34 chisq = 0.613 */ {4.136552, 10.500301, -2.735679, 11.004975, -1.089424, 1.411347, 0.739098 }, /* 35 chisq = 1.183 */ {2.272371, 4.406614, 0.000000, 5.024134, -0.463450, -0.081330, 14.799846 }, /* 36 chisq = 1.719 */ {4.690580, 13.126308, 0.000000, 5.897257, -1.039949, 0.241138, 0.885379 }, /* 37 chisq = 2.344 */ {12.721329, 34.053616, 0.000000, 5.250390, -0.741830, -0.128855, 0.982282 }, /* 38 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 39 chisq = 2.291 */ {0.295099, 5.337376, -0.771608, 4.626912, -0.598705, -0.035524, 1.909027 }, /* 40 chisq = 0.588 */ {8.736829, 23.216066, 0.000000, 5.891793, -1.132144, 0.054811, 0.758395 }, /* 41 chisq = 0.743 */ {2.801573, 5.772323, -1.153520, 5.064821, -0.470740, 0.227218, 3.201044 }, /* 42 chisq = 3.110 */ {1.389064, 6.379728, 0.000000, 4.884761, -0.531152, -0.150485, 9.164878 }, /* 43 chisq = 0.556 */ {1.229283, 11.921393, 0.000000, 5.435559, -0.741080, 0.063042, 0.884616 }, /* 44 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 45 chisq = 2.068 */ {2.133105, 8.096292, 0.000000, 4.857512, -0.492433, -0.031376, 2.203042 }, /* 46 chisq = 2.356 */ {2.254672, 9.771684, 0.000000, 9.137634, -3.176502, 1.388742, 0.618515 }, /* 47 chisq = 2.195 */ {0.871321, 5.145931, 0.000000, 4.888730, -0.652447, -0.423816, 18.927858 }, /* 48 chisq = 7.938 */ {3.836148, 40.221691, 0.000000, 5.199902, -0.876582, -0.118477, 0.849928 }, /* 49 chisq = 1.683 */ {1.145112, 5.532121, 0.000000, 5.720264, -1.401952, 0.354722, 1.010656 }, /* 50 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 51 chisq = 1.191 */ {3.208942, 24.108219, 34.897545, 4.099808, -0.201010, -0.330504, 5.181006 }, /* 52 chisq = 1.043 */ {0.746686, 4.723390, 0.000000, 4.666986, -0.455423, -0.199770, 7.667146 }, /* 53 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 54 chisq = 1.601 */ {3.092335, 3.408154, 0.000000, 4.986257, -0.535889, -0.211274, 18.005024 }, /* 55 chisq = 0.769 */ {2.739344, 3.746556, 0.000000, 5.096318, -0.602019, -0.254327, 17.086578 }, /* 56 chisq = 0.799 */ {1.462646, 5.324323, 0.000000, 5.491094, -1.239940, 0.771551, 1.312135 }, /* 57 chisq = 0.432 */ {2.572002, 7.957870, 0.000000, 5.752131, -1.190189, 0.204412, 1.345068 }, /* 58 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 59 chisq = 0.341 */ {-0.591059, 4.776302, 2.170474, 4.177893, -0.228581, -0.296110, 4.157132 }, /* 60 chisq = 2.198 */ {2.993371, 7.806211, 0.000000, 4.733674, -0.504849, -0.253641, 3.495973 }, /* 61 chisq = 0.728 */ {6.974055, 11.046926, 0.000000, 5.400522, -0.845956, -0.003311, 1.566804 }, /* 62 chisq = 0.624 */ {1.039702, 6.407928, 0.000000, 4.795919, -0.556405, -0.202904, 3.390958 }, /* 63 chisq = 0.159 */ {3.492168, 7.395025, 0.000000, 5.001709, -0.628286, 0.115821, 1.864549 }, /* 64 chisq = 2.837 */ {4.755018, 4.830879, 0.000000, 5.100676, -0.584339, -0.084554, 5.671035 }, /* 65 chisq = 0.908 */ {4.457710, 7.648566, 0.000000, 5.583851, -0.914823, 0.149387, 1.583428 }, /* 66 chisq = 0.760 */ {16.466420, 37.506762, 0.000000, 29.619792, -13.453822, 6.689607, 0.292735 }, /* 67 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 68 chisq = 0.382 */ {5.779137, 8.556554, 0.000000, 5.511747, -0.944686, 0.358148, 1.371945 }, /* 69 chisq = 0.345 */ {4.034569, 3.017766, 0.000000, 5.116028, -0.652882, -0.059847, 14.651407 }, /* 70 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 71 chisq = 0.640 */ {2.951508, 8.678520, 0.000000, 5.047768, -0.728902, 0.083777, 1.748714 }, /* 72 chisq = 2.894 */ {6.777293, 11.909723, 0.000000, 8.988830, -2.876789, 1.112925, 0.629096 }, /* 73 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 74 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 75 chisq = 0.484 */ {3.575914, 4.388129, 0.000000, 5.493224, -0.740422, 0.452501, 2.321947 }, /* 76 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 77 chisq = 1.111 */ {4.219091, 3.148022, 0.000000, 5.213602, -0.645150, -0.093114, 11.489716 }, /* 78 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 79 chisq = 0.690 */ {4.267515, 2.968575, 0.000000, 5.193844, -0.711711, -0.147754, 14.825740 }, /* 80 chisq = 1.0 */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 81 chisq = 0.646 */ {4.068525, 3.760029, 0.000000, 5.210016, -0.605941, 0.047692, 5.233293 }, /* 82 chisq = 0.963 */ {3.737936, 4.034705, 0.000000, 5.344779, -0.710664, 0.101419, 3.749271 }, /* 83 chisq = 0.572 */ {4.769005, 5.441887, 0.000000, 5.225091, -0.764900, 0.152658, 2.499109 }, /* 84 chisq = 0.322 */ {2.100423, 3.939428, 0.000000, 5.049451, -0.576346, -0.283479, 11.473708 }, /* 85 chisq = 1.808 */ {3.118617, 3.616168, 0.000000, 5.282614, -0.607387, 0.046248, 5.545239 }, /* 86 chisq = 2.033 */ {4.099376, 3.114686, 0.000000, 5.180511, -0.644039, -0.079714, 13.794987 }, /* 87 chisq = 1.560 */ {5.654822, 7.941313, 0.000000, 8.396533, -2.416166, 1.554322, 0.780568 }, /* 88 chisq = 1.063 */ {3.272654, 3.422250, 0.000000, 5.144563, -0.643226, -0.147752, 15.670731 }, /* 89 chisq = 1.497 */ {3.857662, 3.680655, 0.000000, 5.276108, -0.667871, -0.105462, 6.934964 }, /* 90 chisq = 0.438 */ {5.123427, 6.270192, 0.000000, 5.044212, -0.723618, -0.007277, 2.547011 }, /* 91 chisq = 0.288 */ {4.787907, 3.746933, 0.000000, 5.299093, -0.591544, 0.051965, 4.950882 }, /* 92 chisq = 0.204 */ {3.369571, 5.241566, 0.000000, 5.334910, -0.696230, 0.215320, 2.109790 }, /* 93 chisq = 0.560 */ {4.071519, 3.212398, 0.000000, 5.117411, -0.641382, -0.083187, 11.416026 }, /* 94 chisq = 1.367 */ {6.371714, 10.609759, 0.000000, 5.324584, -0.641072, 0.091840, 1.595554 }, /* 95 chisq = 1.508 */ {4.199741, 3.972734, 0.000000, 5.330881, -0.600000, 0.146681, 4.014591 }, /* 96 chisq = 0.324 */ {6.766545, 8.847396, 0.000000, 5.181152, -0.634533, 0.025916, 2.455332 }, /* 97 chisq = 0.281 */ {5.629009, 7.823334, 0.000000, 5.574892, -0.858300, 0.196157, 1.624060 }, /* 98 chisq = 1.937 */ {4.722284, 5.339112, 0.000000, 5.149685, -0.634661, -0.043345, 3.741964 }, /* 99 chisq = 0.266 */ {4.122186, 2.707496, 0.000000, 5.191277, -0.722662, -0.146064, 15.567705 }, /* 100 chisq = 1.247 */ {5.239048, 7.126645, 0.000000, 7.623700, -2.480334, 1.118232, 0.867630 }, /* 101 chisq = 1.291 */ {3.968395, 3.166231, 0.000000, 5.221490, -0.589120, -0.016509, 10.289873 }, /* 102 chisq = 0.876 */ {3.297637, 3.603367, 0.000000, 5.156119, -0.666150, -0.152755, 7.251896 }, /* 103 chisq = 0.300 */ {9.236691, 16.579236, 0.000000, 8.642802, -2.343229, 0.871002, 0.611731 }, /* 104 chisq = 0.322 */ {5.904796, 9.754144, 0.000000, 6.350622, -1.304800, 0.410343, 1.035687 }, /* 105 chisq = 1.275 */ {3.787965, 3.376772, 0.000000, 5.183945, -0.640034, -0.102485, 9.230703 }, /* 106 chisq = 0.761 */ {4.026565, 2.946358, 0.000000, 5.289412, -0.607222, -0.074520, 13.254879 }, /* 107 chisq = 1.821 */ {4.066532, 4.581130, 0.000000, 5.231305, -0.687594, -0.035834, 3.840558 }, /* 108 chisq = 0.471 */ {3.874968, 4.221358, 0.000000, 5.214631, -0.675976, 0.027162, 3.836486 }, /* 109 chisq = 1.458 */ {8.002132, 9.468869, 0.000000, 5.942355, -1.313073, 0.011164, 1.078875 }, /* 110 chisq = 0.152 */ {3.584799, 3.603469, 0.000000, 5.179186, -0.667243, -0.186713, 13.007171 }, /* 111 chisq = 1.0 first ring */ {0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000, 0.000000 }, /* 2 chisq = 1.780 */ {3.363495, 4.057447, 0.000000, 5.099786, -0.622330, -0.165568, 5.581331 }, /* 3 chisq = 0.846 */ {3.095003, 2.356619, 0.000000, 4.684074, -0.580081, -0.184067, 20.704745 }, /* 4 chisq = 0.445 */ {4.132085, 7.261742, 0.000000, 5.759001, -0.665104, -0.021448, 2.328372 }, /* 5 chisq = 0.446 */ {3.099310, 8.550535, 0.000000, 4.861103, -0.648275, -0.026817, 2.070629 }, /* 6 chisq = 0.926 */ {3.528680, 3.506597, 0.000000, 5.345194, -0.489192, -0.117327, 8.235985 }, /* 7 chisq = 1.785 */ {1.200624, 9.906573, 0.000000, 4.442448, -0.481633, -0.089079, 2.121392 }, /* 8 chisq = 0.917 */ {1.236546, 6.146697, 0.000000, 4.514243, -0.455956, -0.126097, 3.099657 }, /* 9 chisq = 1.089 */ {3.568001, 10.949971, 0.000000, 5.465161, -0.664510, -0.077517, 1.563468 }, /* 10 chisq = 1.036 */ {4.108838, 2.223499, 0.000000, 4.954450, -0.558626, -0.080420, 11.704817 }, /* 11 chisq = 1.647 */ {2.916001, 3.208322, 0.000000, 4.552745, -0.601580, -0.193408, 10.225078 }, /* 12 chisq = 0.402 */ {4.298961, 4.220965, 0.000000, 5.177173, -0.643644, -0.093431, 4.970709 }, /* 13 chisq = 0.969 */ {4.534408, 3.096070, 0.000000, 5.170505, -0.634157, -0.068246, 6.705700 }, /* 14 chisq = 0.000 */ {5.056166, 2.928681, 0.000000, 5.886193, -0.584274, -0.051499, 8.154615 }, /* 15 chisq = 0.493 */ {4.520018, 2.037066, 0.000000, 5.227524, -0.694607, -0.162062, 16.072930 }, /* 16 chisq = 0.487 */ {4.451283, 1.735243, 0.660288, 5.218744, -0.605471, -0.056662, 9.363573 }, /* 17 chisq = 0.491 */ {4.415205, 6.897258, 0.000000, 5.213504, -0.676006, -0.110160, 3.192231 } }; /* time window for prompt and 2d random peak, not shifted! */ /* 4096 is offsett */ static int twin[NDET+1][4] = { 0, 0, 0, 0, /* 0 */ 0, 0, 0, 0, /* 1 */ 0, 0, 0, 0, /* 2 */ 0, 0, 0, 0, /* 3 */ 0, 0, 0, 0, /* 4 */ 3808, 4123, 3067, 3382, /* 5 */ 3808, 4101, 3067, 3360, /* 7 */ 3811, 4084, 3081, 3354, /* 8 */ 3808, 4101, 3067, 3360, /* 9 */ 3811, 4090, 3067, 3346, /* 10 */ 3811, 4115, 3081, 3385, /* 11 */ 3811, 4082, 3081, 3352, /* 12 */ 3803, 4093, 3067, 3357, /* 13 */ 3797, 4101, 3081, 3385, /* 14 */ 3806, 4093, 3067, 3354, /* 15 */ 0, 0, 0, 0, /* 16 */ 3808, 4132, 3081, 3405, /* 17 */ 3808, 4082, 3081, 3355, /* 18 */ 3808, 4110, 3081, 3383, /* 19 */ 3808, 4110, 3081, 3383, /* 20 */ 3814, 4098, 3081, 3365, /* 21 */ 0, 0, 0, 0, /* 22 */ 3814, 4098, 3081, 3365, /* 23 */ 3803, 4087, 3081, 3365, /* 24 */ 3803, 4101, 3081, 3379, /* 25 */ 3803, 4101, 3081, 3379, /* 26 */ 3803, 4143, 3067, 3407, /* 27 */ 3803, 4121, 3067, 3385, /* 28 */ 3803, 4082, 3067, 3346, /* 29 */ 0, 0, 0, 0, /* 30 */ 3811, 4107, 3067, 3363, /* 31 */ 3811, 4107, 3067, 3363, /* 32 */ 0, 0, 0, 0, /* 33 */ 3811, 4107, 3067, 3363, /* 34 */ 3811, 4107, 3067, 3363, /* 35 */ 3800, 4098, 3067, 3365, /* 36 */ 3811, 4121, 3053, 3363, /* 37 */ 0, 0, 0, 0, /* 38 */ 3811, 4096, 3067, 3352, /* 39 */ 3811, 4084, 3067, 3340, /* 42 */ 3789, 4123, 3073, 3407, /* 41 */ 3814, 4096, 3073, 3355, /* 40 */ 3814, 4096, 3073, 3355, /* 43 */ 0, 0, 0, 0, /* 44 */ 3814, 4118, 3073, 3377, /* 45 */ 3814, 4082, 3073, 3341, /* 46 */ 3797, 4126, 3073, 3402, /* 47 */ 3797, 4154, 3053, 3410, /* 48 */ 3867, 4082, 3415, 3630, /* 49 */ 3820, 4082, 3081, 3343, /* 50 */ 3778, 4110, 3081, 3413, /* 51 */ 3800, 4082, 3081, 3363, /* 52 */ 0, 0, 0, 0, /* 53 */ 3800, 4146, 3081, 3427, /* 54 */ 3814, 4110, 3081, 3377, /* 55 */ 3800, 4132, 3081, 3413, /* 56 */ 3806, 4098, 3081, 3373, /* 57 */ 0, 0, 0, 0, /* 58 */ 3797, 4121, 3081, 3405, /* 59 */ 3803, 4110, 3081, 3388, /* 60 */ 3803, 4110, 3081, 3388, /* 61 */ 3803, 4079, 3081, 3357, /* 62 */ 3808, 4107, 3081, 3380, /* 63 */ 3808, 4132, 3081, 3405, /* 64 */ 3806, 4101, 3081, 3376, /* 65 */ 3806, 4082, 3081, 3357, /* 66 */ 0, 0, 0, 0, /* 67 */ 3806, 4143, 3081, 3418, /* 68 */ 3806, 4132, 3081, 3407, /* 69 */ 0, 0, 0, 0, /* 70 */ 3795, 4123, 3081, 3409, /* 71 */ 3795, 4123, 3059, 3387, /* 72 */ 0, 0, 0, 0, /* 73 */ 0, 0, 0, 0, /* 74 */ 3811, 4093, 3059, 3341, /* 75 */ 0, 0, 0, 0, /* 76 */ 3781, 4115, 3059, 3393, /* 77 */ 0, 0, 0, 0, /* 78 */ 3781, 4132, 3081, 3410, /* 79 */ 0, 0, 0, 0, /* 80 */ 3803, 4118, 3081, 3396, /* 81 */ 3808, 4073, 3081, 3346, /* 82 */ 3795, 4118, 3081, 3404, /* 83 */ 3811, 4098, 3081, 3368, /* 84 */ 3811, 4126, 3081, 3396, /* 85 */ 3800, 4143, 3081, 3424, /* 86 */ 3806, 4123, 3081, 3398, /* 87 */ 3817, 4137, 3081, 3401, /* 88 */ 3803, 4129, 3053, 3379, /* 89 */ 3803, 4132, 3053, 3382, /* 90 */ 3797, 4121, 3081, 3405, /* 91 */ 3811, 4093, 3081, 3363, /* 92 */ 3800, 4123, 3081, 3404, /* 93 */ 3808, 4132, 3081, 3405, /* 94 */ 3808, 4132, 3081, 3405, /* 95 */ 3808, 4132, 3081, 3405, /* 96 */ 3817, 4115, 3081, 3379, /* 97 */ 3814, 4123, 3081, 3390, /* 98 */ 3814, 4123, 3081, 3390, /* 99 */ 3814, 4093, 3081, 3360, /* 100 */ 3800, 4126, 3081, 3407, /* 101 */ 3808, 4096, 3081, 3369, /* 102 */ 3820, 4110, 3081, 3371, /* 103 */ 3808, 4104, 3081, 3377, /* 104 */ 3817, 4112, 3081, 3376, /* 105 */ 3806, 4121, 3081, 3396, /* 106 */ 3806, 4137, 3081, 3412, /* 107 */ 3806, 4115, 3081, 3390, /* 108 */ 3792, 4149, 3081, 3438, /* 109 */ 3817, 4104, 3081, 3368 /* 110 */ };