storetofTrgWindow.C

// macro to upload tofr5 INL tables to database
//
// based on
//  http://www.star.bnl.gov/STAR/comp/db/StoreDbTable.cc.html
//
// Jing Liu, 02/18/2005 
//


// #include "StDbLib/StDbManager.hh"
// #include "StDbLib/StDbConfigNode.hh"
// #include "StDbLib/StDbTable.h"
// #include "StDbLib/StDbDefs.hh"

#include <iostream>
#include <fstream>
#include <string>
#include <iomanip>
using namespace std;


void storetofTrgWindow() 
{
        // year8
        //  const Int_t NTRAY = 5;
        // year9
        //  const Int_t NTRAY = 94;
        // year10+
        const Int_t NTRAY = 120;
        const Int_t NVPDTRAY = 2;
        const Int_t NMAX = 120;
        const Int_t NMTD = 2; //there is only one MTD, but it has ID=124 (and is only used in Run 10)

        //-- load dBase and Table definition libraries
        gSystem->Load("St_base");
        gSystem->Load("StChain");
        gSystem->Load("StUtilities");
        gSystem->Load("St_Tables.so");
        gSystem->Load("StDbLib.so");
        gSystem->Load("libStDb_Tables.so");

        //-- get the singleton manager
        StDbManager* dbManager = StDbManager::Instance();

        //-- connect to the db & get an empty container
        StDbConfigNode* configNode = dbManager->initConfig("Calibrations_tof");

        // Timestamp - not read from input/timestamp!
        TString ZStoreTime = "2019-01-01 00:00:00";

        //-- add table to the container with descriptor given by Database
        StDbTable* tofTrgWindow = configNode->addDbTable("tofTrgWindow");

        //-- fill structures & store times
        tofTrgWindow_st *trg = new tofTrgWindow_st[NMAX+NVPDTRAY+NMTD];

        double cutlow[122] ={
                1007, 1007, 1007, 1008, 1015, 1015, 1015, 1015, 1014, 1015, 1015, 1015, 1013, 1013, 1015, 1016, 1008, 1008, 1008, 1009, 934, 935, 935, 934, 932, 934, 940, 942, 939, 0, 939, 948, 949, 0, 949, 947, 949, 947, 949, 941, 942, 942, 942, 939, 932, 932, 934, 934, 935, 932, 990, 993, 990, 998, 998, 1000, 1000, 1000, 1007, 1007, 1003, 1003, 994, 995, 995, 999, 997, 999, 997, 999, 999, 1005, 1008, 1004, 1007, 1004, 1012, 1014, 1014, 1013, 1012, 1013, 1013, 1013, 1004, 1004, 1004, 1004, 1004, 997, 997, 998, 998, 998, 998, 1010, 1010, 1010, 1010, 1019, 1019, 0, 1019, 1018, 1019, 1018, 1018, 1019, 1018, 1018, 1018, 1010, 1012, 1010, 1010, 1009, 1010, 1001, 1002, 1003, 958, 1019};

        double cuthi[122] ={
                1122, 1122, 1122, 1123, 1130, 1130, 1130, 1130, 1129, 1130, 1130, 1130, 1128, 1128, 1130, 1131, 1123, 1123, 1123, 1124, 1049, 1050, 1050, 1049, 1047, 1049, 1055, 1057, 1054, 0, 1054, 1063, 1064, 0, 1064, 1062, 1064, 1062, 1064, 1056, 1057, 1057, 1057, 1054, 1047, 1047, 1049, 1049, 1050, 1047, 1105, 1108, 1105, 1113, 1113, 1115, 1115, 1115, 1122, 1122, 1118, 1118, 1109, 1110, 1110, 1114, 1112, 1114, 1112, 1114, 1114, 1120, 1123, 1119, 1122, 1119, 1127, 1129, 1129, 1128, 1127, 1128, 1128, 1128, 1119, 1119, 1119, 1119, 1119, 1112, 1112, 1113, 1113, 1113, 1113, 1125, 1125, 1125, 1125, 1134, 1134, 0, 1134, 1133, 1134, 1133, 1133, 1134, 1133, 1133, 1133, 1125, 1127, 1125, 1125, 1124, 1125, 1116, 1117, 1118, 1073, 1134};

        // ---------------------------------------------
        //- Prepare table
        for(int i=0;i<NMAX+NVPDTRAY;i++) 
        {
                trg[i].trgWindow_Min = cutlow[i];
                trg[i].trgWindow_Max = cuthi[i];
                cout << " tray = " << i+1 << " min = " << trg[i].trgWindow_Min << " max = " << trg[i].trgWindow_Max << endl;
        }

        // Store data in table
        tofTrgWindow->SetTable((char*)trg, NMAX+NVPDTRAY);
        //- set store time
        dbManager->setStoreTime(ZStoreTime.Data());
        // Store table in dBase
        cout<<" here "<<endl;
        dbManager->storeDbTable(tofTrgWindow);
        cout<<"uploaded"<<endl;
}



//----------------------------------------
// History of time stamps (in GMT)
//----------------------------------------
//  TString ZStoreTime = "2007-11-01 00:00:01";
//  TString ZStoreTime = "2009-02-01 00:00:00";  // y9 1st version very loose
//  TString ZStoreTime = "2009-02-01 00:00:01";  // y9 2st version  

//  TString ZStoreTime = "2009-03-26 00:00:00";  // y9 500GeV day 85-92
//  TString ZStoreTime = "2009-04-23 21:00:00";  // 10113066-10146091
//  TString ZStoreTime = "2009-05-27 00:00:00";  // 10147124-10149036
//  TString ZStoreTime = "2009-05-29 14:00:00";  // 10149077-10150023
//  TString ZStoreTime = "2009-05-30 15:30:00";  // 10150024-10152010
//  TString ZStoreTime = "2009-06-03 16:00:00";  // 10154045-...
//
//  TString ZStoreTime = "2009-11-01 00:00:00"; // run 10 preliminary
//  TString ZStoreTime = "2009-11-01 00:01:00"; // updated version, includes MTD
//  TString ZStoreTime = "2009-11-01 00:01:01"; // updated version, (2nd attempt)includes MTD
//  TString ZStoreTime = "2010-12-10 00:00:00"; // run 11 new timing (TOCK)
//  TString ZStoreTime = "2011-12-20 00:00:01"; // run 12 new timing
//  TString ZStoreTime = "2012-03-21 19:50:00"; // run 12 new VPD timing (effective from run 13081013)
//  TString ZStoreTime = "2012-03-21 19:50:01"; // run 12 new VPD timing (corrected, effective from run 13081013)
//  TString ZStoreTime = "2012-12-20 00:00:01"; // run 13 new timing
//  TString ZStoreTime = "2013-12-20 00:00:01"; // run 14 new timing (Au+Au@14.5GeV)
//  TString ZStoreTime = "2014-03-14 00:00:00"; // run 14 new timing (Au+Au@200GeV)
//  TString ZStoreTime = "2015-01-06 00:00:00"; // run 15 initial windows pp200


// year 10
// unsigned short cutlow[124]={
// 2828,2825,2826,2826,2832,2831,2835,2831,2836,2832,2832,2829,2833,2830,2833,
// 2831,2824,2826,2825,2824,2827,2836,2833,2830,2829,2831,2838,2833,2837,2838,
// 2836,2838,2850,2841,2847,2845,2845,2846,2843,2838,2837,2836,2835,2835,2829,
// 2824,2830,2829,2828,2827,2806,2807,2807,2814,2815,2818,2813,2821,2826,2828,
// 2891,2901,2891,2894,2891,2897,2894,2900,2899,2899,2898,2904,2905,2904,2898,
// 2904,2907,2912,2916,2916,2908,2914,2908,2914,2903,2902,2907,2901,2901,2898,
// 2895,2894,2897,2896,2896,2906,2903,2908,2907,2914,2911,2698,2913,2915,2915,
// 2916,2914,2913,2913,2912,2914,2909,2905,2908,2902,2898,2902,2895,2899,2899,
// 2857,2921,0,2900};
// unsigned short cuthi[124]={
// 2903,2900,2901,2901,2907,2906,2910,2906,2911,2907,2907,2904,2908,2905,2908,
// 2906,2899,2901,2900,2899,2902,2911,2908,2905,2904,2906,2913,2908,2912,2913,
// 2911,2913,2925,2916,2922,2920,2920,2921,2918,2913,2912,2911,2910,2910,2904,
// 2899,2905,2904,2903,2902,2881,2882,2882,2889,2890,2893,2888,2896,2901,2903,
// 2966,2976,2966,2969,2966,2972,2969,2975,2974,2974,2973,2979,2980,2979,2973,
// 2979,2982,2987,2991,2991,2983,2989,2983,2989,2978,2977,2982,2976,2976,2973,
// 2970,2969,2972,2971,2971,2981,2978,2983,2982,2989,2986,2773,2988,2990,2990,
// 2991,2989,2988,2988,2987,2989,2984,2980,2983,2977,2973,2977,2970,2974,2974,
// 2932,2996,0,2975};

// year 11
//unsigned short cutlow[124]={
//2816,2817,2813,2817,2820,2822,2820,2823,2823,2813,2820,2824,2820,2817,2825,
//2824,2815,2813,2810,2816,2824,2824,2823,2825,2698,2820,2831,2832,2828,2829,
//2830,2837,2839,2835,2840,2837,2840,2838,2840,2834,2836,2824,2826,2830,2822,
//2808,2822,2823,2824,2823,2799,2800,2799,2798,2809,2811,2810,2793,2803,2816,
//2890,2892,2886,2885,2879,2889,2886,2887,2891,2892,2883,2892,2895,2883,2897,
//2884,2904,2904,2902,2889,2901,2904,2904,2902,2884,2897,2895,2892,2895,2889,
//2885,2881,2887,2890,2886,2898,2900,2901,2902,2908,2910,2698,2907,2909,2910,
//2907,2908,2907,2906,2901,2909,2898,2901,2898,2895,2899,2900,2887,2888,2887,
//2850,2913, 0,0};
//unsigned short cuthi[124]={
//2901,2902,2898,2902,2905,2907,2905,2908,2908,2898,2905,2909,2905,2902,2910,
//2909,2900,2898,2895,2901,2909,2909,2908,2910,2783,2905,2916,2917,2913,2914,
//2915,2922,2924,2920,2925,2922,2925,2923,2925,2919,2921,2909,2911,2915,2907,
//2893,2907,2908,2909,2908,2884,2885,2884,2883,2894,2896,2895,2878,2888,2901,
//2975,2977,2971,2970,2964,2974,2971,2972,2976,2977,2968,2977,2980,2968,2982,
//2969,2989,2989,2987,2974,2986,2989,2989,2987,2969,2982,2980,2977,2980,2974,
//2970,2966,2972,2975,2971,2983,2985,2986,2987,2993,2995,2783,2992,2994,2995,
//2992,2993,2992,2991,2986,2994,2983,2986,2983,2980,2984,2985,2972,2973,2972,
//2935,2998,0,0};

// year 12
// unsigned short cutlow[124]={
// 2881,2873,2873,2873,2893,2876,2882,2878,2888,2881,2885,2878,2884,2876,2887,
// 2886,2877,2876,2880,2865,2804,2802,2799,2802,2798,2800,2810,2811,2807,2813,
// 2817,2812,2796,2806,2815,2812,2808,2809,2815,2807,2805,2811,2805,2803,2799,
// 2790,2801,2801,2804,2807,2860,2862,2862,2868,2871,2865,2870,2867,2872,2871,
// 2872,2864,2866,2864,2863,2868,2865,2871,2871,2862,2861,2867,2875,2871,2873,
// 2871,2878,2880,2888,2883,2883,2872,2884,2882,2874,2873,2874,2871,2867,2860,
// 2861,2866,2863,2863,2698,2876,2880,2879,2881,2878,2877,2698,2891,2890,2887,
// 2880,2891,2892,2881,2887,2885,2874,2877,2878,2874,2877,2888,2876,2874,2865,
// 2826,2889,0,0};
// unsigned short cuthi[124]={
// 2966,2958,2958,2958,2978,2961,2967,2963,2973,2966,2970,2963,2969,2961,2972,
// 2971,2962,2961,2965,2950,2889,2887,2884,2887,2883,2885,2895,2896,2892,2898,
// 2902,2897,2881,2891,2900,2897,2893,2894,2900,2892,2890,2896,2890,2888,2884,
// 2875,2886,2886,2889,2892,2945,2947,2947,2953,2956,2950,2955,2952,2957,2956,
// 2957,2949,2951,2949,2948,2953,2950,2956,2956,2947,2946,2952,2960,2956,2958,
// 2956,2963,2965,2973,2968,2968,2957,2969,2967,2959,2958,2959,2956,2952,2945,
// 2946,2951,2948,2948,2783,2961,2965,2964,2966,2963,2962,2783,2976,2975,2972,
// 2965,2976,2977,2966,2972,2970,2959,2962,2963,2959,2962,2973,2961,2959,2950,
// 2911,2974,0,0};
// 
// // year 12 - updated VPD window to lower TOF deadtime (pp500 -- March 22, 2012)
// //  cutlow[120]= 230 ; cuthi[120]= 290;  // VPD West
// //  cutlow[121]= 245 ; cuthi[121]= 305;  // VPD East
// // year 12 - correct updated VPD window to lower TOF deadtime (pp500 -- March 22, 2012, [5/1'12])
//   cutlow[120]= 2834 ; cuthi[120]= 2919;  // VPD West
//   cutlow[121]= 2896 ; cuthi[121]= 2981;  // VPD East


// year 13
// notes: trays 46, 102, and 106 are interpolated values (trays are dead)
// double cutlow[122]={
// 2815,2814,2815,2815,2822,2821,2824,2823,2821,2824,2822,2821,2821,2820,2823,
// 2822,2814,2816,2815,2814,2740,2741,2742,2740,2739,2741,2746,2749,2746,2746,
// 2745,2754,2755,2754,2756,2754,2756,2753,2755,2750,2750,2748,2749,2745,2738,
// 2738,2739,2739,2741,2738,2796,2801,2795,2804,2806,2807,2807,2807,2814,2814,
// 2808,2808,2801,2803,2802,2806,2803,2806,2804,2806,2806,2811,2814,2810,2814,
// 2810,2819,2820,2822,2820,2820,2819,2818,2820,2810,2812,2810,2811,2810,2804,
// 2804,2804,2804,2804,2805,2817,2816,2817,2817,2825,2826,2825,2825,2825,2825,
// 2825,2826,2826,2825,2824,2825,2816,2818,2817,2816,2818,2816,2807,2808,2808,
// 2764,2827};
// double cuthi[122]={
// 2900,2899,2900,2900,2907,2906,2909,2908,2906,2909,2907,2906,2906,2905,2908,
// 2907,2899,2901,2900,2899,2825,2826,2827,2825,2824,2826,2831,2834,2831,2831,
// 2830,2839,2840,2839,2841,2839,2841,2838,2840,2835,2835,2833,2834,2830,2823,
// 2823,2824,2824,2826,2823,2881,2886,2880,2889,2891,2892,2892,2892,2899,2899,
// 2893,2893,2886,2888,2887,2891,2888,2891,2889,2891,2891,2896,2899,2895,2899,
// 2895,2904,2905,2907,2905,2905,2904,2903,2905,2895,2897,2895,2896,2895,2889,
// 2889,2889,2889,2889,2890,2902,2901,2902,2902,2910,2911,2910,2910,2910,2910,
// 2910,2911,2911,2910,2909,2910,2901,2903,2902,2901,2903,2901,2892,2893,2893,
// 2849,2912};


// year 14 -- AuAu@14.5 initial
// double cutlow[122]={
// 2909,2910,2910,2911,2918,2917,2916,2917,2916,2917,2917,2916,2915,2916,2918,
// 2918,2909,2910,2910,2911,2836,2852,2837,2836,2834,2834,2842,2845,2841,2841,
// 2841,2850,2851,2850,-23,2850,2851,2849,2851,2843,2843,2844,2844,2841,2834,
// 2834,2835,2836,2837,2834,2892,2894,2892,2900,2900,2903,2903,2902,2909,2909,
// 2906,2906,2898,2898,2899,2903,2901,2902,2901,2903,2903,2908,2911,2907,2910,
// 2908,2915,2918,2918,2916,2916,2916,2916,2916,2907,2907,2907,2907,2907,2900,
// 2901,2902,2902,2901,2902,2914,2914,2914,2914,2923,2922,-23,-23,2922,2922,
// 2922,2921,2923,2922,2922,2922,2914,2916,2913,2913,2914,2913,2906,2906,2905,
// 2859,2921};
// double cuthi[122]={
// 2994,2995,2995,2996,3003,3002,3001,3002,3001,3002,3002,3001,3000,3001,3003,
// 3003,2994,2995,2995,2996,2921,2937,2922,2921,2919,2919,2927,2930,2926,2926,
// 2926,2935,2936,2935,62,2935,2936,2934,2936,2928,2928,2929,2929,2926,2919,
// 2919,2920,2921,2922,2919,2977,2979,2977,2985,2985,2988,2988,2987,2994,2994,
// 2991,2991,2983,2983,2984,2988,2986,2987,2986,2988,2988,2993,2996,2992,2995,
// 2993,3000,3003,3003,3001,3001,3001,3001,3001,2992,2992,2992,2992,2992,2985,
// 2986,2987,2987,2986,2987,2999,2999,2999,2999,3008,3007,62,62,3007,3007,
// 3007,3006,3008,3007,3007,3007,2999,3001,2998,2998,2999,2998,2991,2991,2990,
// 2944,3006};


// year 14 -- Au+Au@200GeV
// double cutlow[122]={
// 2927,2927,2927,2927,2933,2934,2933,2934,2933,2934,2933,2934,2932,2932,2935,
// 2934,2927,2927,2927,-23,2853,-23,2854,2853,2852,2852,2859,2861,2858,2858,
// 2858,2866,2868,2866,-23,2866,2868,2866,2868,2861,2861,2861,2861,2858,2851,
// 2852,2853,2853,2854,2852,2908,2912,2908,2916,2917,2919,2919,2919,2927,2927,
// 2922,2922,2914,2914,2914,2918,2916,2918,2916,2918,2918,2924,2927,2923,2926,
// 2923,2931,2933,2933,2931,2931,2931,2931,2931,2923,2924,2923,2923,2923,2916,
// 2916,2916,2917,2916,2917,2930,2930,2930,2930,2938,2938,-23,-23,2938,2937,
// 2937,2937,2938,2938,2937,2937,2929,2931,2929,2929,2929,2929,2921,2921,2922,
// 2877,2937};
// double cuthi[122]={
// 3012,3012,3012,3012,3018,3019,3018,3019,3018,3019,3018,3019,3017,3017,3020,
// 3019,3012,3012,3012,62,2938,62,2939,2938,2937,2937,2944,2946,2943,2943,
// 2943,2951,2953,2951,62,2951,2953,2951,2953,2946,2946,2946,2946,2943,2936,
// 2937,2938,2938,2939,2937,2993,2997,2993,3001,3002,3004,3004,3004,3012,3012,
// 3007,3007,2999,2999,2999,3003,3001,3003,3001,3003,3003,3009,3012,3008,3011,
// 3008,3016,3018,3018,3016,3016,3016,3016,3016,3008,3009,3008,3008,3008,3001,
// 3001,3001,3002,3001,3002,3015,3015,3015,3015,3023,3023,62,62,3023,3022,
// 3022,3022,3023,3023,3022,3022,3014,3016,3014,3014,3014,3014,3006,3006,3007,
// 2962,3022};
// 
// year 15 -- pp 200
// start-time 2015-01-06 00:00:00
// double cutlow[122] ={
// 2819,2819,2819,2819,2826,2827,2827,2827,2826,2826,
// 2826,2826,2825,2826,2828,2828,2819,2819,2820,2821,
// 2746,2747,2747,2746,2745,2745,2752,2755,2751,2751,
// 2751,2759,2762,2760,2760,2759,2761,-25,2761,2753,
// -25,2754,2753,2751,2744,2744,2747,2747,2747,2745,
// 2802,2804,2803,2810,2810,2812,2812,2812,2819,2819,
// 2816,2815,2808,2807,2807,2811,2810,2811,2808,2812,
// 2812,2817,2821,2816,2819,2817,2825,2827,2827,2827,
// 2824,2825,2824,2825,2816,2816,2817,2817,2816,2809,
// 2809,2810,2811,2809,2810,2823,2824,2823,2823,2831,
// 2831,-25,2833,2830,2831,2831,2830,2832,2830,2830,
// 2830,2821,2826,2824,2821,2822,2823,2814,2815,2814,
// 2770,2830 };

// double cuthi[122] ={
// 2929,2929,2929,2929,2936,2937,2937,2937,2936,2936,
// 2936,2936,2935,2936,2938,2938,2929,2929,2930,2931,
// 2856,2857,2857,2856,2855,2855,2862,2865,2861,2861,
// 2861,2869,2872,2870,2870,2869,2871,85,2871,2863,
// 85,2864,2863,2861,2854,2854,2857,2857,2857,2855,
// 2912,2914,2913,2920,2920,2922,2922,2922,2929,2929,
// 2926,2925,2918,2917,2917,2921,2920,2921,2918,2922,
// 2922,2927,2931,2926,2929,2927,2935,2937,2937,2937,
// 2934,2935,2934,2935,2926,2926,2927,2927,2926,2919,
// 2919,2920,2921,2919,2920,2933,2934,2933,2933,2941,
// 2941,85,2943,2940,2941,2941,2940,2942,2940,2940,
// 2940,2931,2936,2934,2931,2932,2933,2924,2925,2924,
// 2880,2940 };
// 
// // Run16 AuAu200 - commissioning
// double cutlow[122] = {2820,2819,2820,2820,2827,2827,2827,2828,2827,2827,2827,2827,2826,2826,2828,2828,2821,2820,2821,2821,2747,2748,2748,2747,2744,2746,2752,2755,2752,2752,2752,2760,2762,2760,2761,2760,2762,2760,2761,2754,2754,2754,2754,2752,2745,2745,2747,2747,2748,2745,2802,2804,2802,2810,2810,2813,2812,2813,2819,2819,2815,2815,2807,2808,2808,2812,2810,2812,2810,2812,2812,2818,2821,2817,2819,2817,2826,2827,2827,2826,2825,2825,2826,2825,2817,2817,2817,2817,2817,2810,2810,2810,2811,2810,2811,2824,2824,2824,2824,2832,2832,-25,2832,2831,-25,2831,2831,2832,2831,2831,2831,2823,2826,2823,2823,2823,2823,2815,2815,2815,2772,2831 };
// double cuthi[122] = {2930,2929,2930,2930,2937,2937,2937,2938,2937,2937,2937,2937,2936,2936,2938,2938,2931,2930,2931,2931,2857,2858,2858,2857,2854,2856,2862,2865,2862,2862,2862,2870,2872,2870,2871,2870,2872,2870,2871,2864,2864,2864,2864,2862,2855,2855,2857,2857,2858,2855,2912,2914,2912,2920,2920,2923,2922,2923,2929,2929,2925,2925,2917,2918,2918,2922,2920,2922,2920,2922,2922,2928,2931,2927,2929,2927,2936,2937,2937,2936,2935,2935,2936,2935,2927,2927,2927,2927,2927,2920,2920,2920,2921,2920,2921,2934,2934,2934,2934,2942,2942,85,2942,2941,85,2941,2941,2942,2941,2941,2941,2933,2936,2933,2933,2933,2933,2925,2925,2925,2882,2941 };
// 

// // Run16 Update due to changes in 17056042
//   double cutlow[122] ={
// 1020,1020,1020,1020,1027,1027,1027,1027,1027,1027,
// 1028,1027,1026,1026,1028,1028,1021,1021,1021,1021,
// 946,947,947,946,944,946,952,954,951,951,
// 951,959,961,960,961,959,961,959,961,954,
// 954,954,954,951,944,944,946,946,947,944,
// 1002,1005,1002,1010,1011,1013,1013,1013,1020,1020,
// 1015,1015,1007,1008,1008,1012,1010,1012,1011,1013,
// 1013,1018,1021,1017,1020,1017,1026,1027,1027,1026,
// 1026,1026,1025,1026,1016,1017,1017,1016,1017,1010,
// 1010,1011,1011,1010,1011,1024,1024,1024,1024,1032,
// 1032,-25,1032,1031,-25,1031,1030,1032,1031,1031,
// 1031,1023,1025,1023,1023,1023,1023,1014,1015,1015,
// 971,1031 };

// double cuthi[122] ={
// 1130,1130,1130,1130,1137,1137,1137,1137,1137,1137,
// 1138,1137,1136,1136,1138,1138,1131,1131,1131,1131,
// 1056,1057,1057,1056,1054,1056,1062,1064,1061,1061,
// 1061,1069,1071,1070,1071,1069,1071,1069,1071,1064,
// 1064,1064,1064,1061,1054,1054,1056,1056,1057,1054,
// 1112,1115,1112,1120,1121,1123,1123,1123,1130,1130,
// 1125,1125,1117,1118,1118,1122,1120,1122,1121,1123,
// 1123,1128,1131,1127,1130,1127,1136,1137,1137,1136,
// 1136,1136,1135,1136,1126,1127,1127,1126,1127,1120,
// 1120,1121,1121,1120,1121,1134,1134,1134,1134,1142,
// 1142,85,1142,1141,85,1141,1140,1142,1141,1141,
// 1141,1133,1135,1133,1133,1133,1133,1124,1125,1125,
// 1081,1141 };
//
//Run17 pp 510
//double cutlow[122]={
// 1018,1018,1018,1019,1027,1026,1027,1027,1025,1026,1026,1027,1026,1025,1028,
// 1027,1018,1020,1019,1019,946,946,946,945,943,945,953,954,951,949,
// 952,960,962,960,960,-23,961,958,960,953,953,952,951,951,-23,
// 943,946,945,947,943,1001,1004,1001,1008,1009,1011,1012,1012,1017,1018,
// 1014,1013,1005,1006,1005,1011,1009,1011,1009,1012,1010,1016,1019,1015,1018,
// 1017,1025,1026,1026,-23,1025,1024,1024,1024,1016,1017,1016,1016,1015,1008,
// 1008,1010,1009,1009,1009,1022,1022,1022,1022,1031,1030,-23,1030,1031,1030,
// 1029,1029,1030,1029,1030,1029,1022,1025,1022,1021,1021,1021,1013,1013,1014,
// 970,1030};
// double cuthi[122]={
// 1103,1103,1103,1104,1112,1111,1112,1112,1110,1111,1111,1112,1111,1110,1113,
// 1112,1103,1105,1104,1104,1031,1031,1031,1030,1028,1030,1038,1039,1036,1034,
// 1037,1045,1047,1045,1045,62,1046,1043,1045,1038,1038,1037,1036,1036,62,
// 1028,1031,1030,1032,1028,1086,1089,1086,1093,1094,1096,1097,1097,1102,1103,
// 1099,1098,1090,1091,1090,1096,1094,1096,1094,1097,1095,1101,1104,1100,1103,
// 1102,1110,1111,1111,62,1110,1109,1109,1109,1101,1102,1101,1101,1100,1093,
// 1093,1095,1094,1094,1094,1107,1107,1107,1107,1116,1115,62,1115,1116,1115,
// 1114,1114,1115,1114,1115,1114,1107,1110,1107,1106,1106,1106,1098,1098,1099,
// 1055,1115};
//
// Run18
//double cutlow[122] ={
//      1014, 1013, 1014, 1014, 1022, 1022, 1022, 0, 1021, 1022, 1022, 1021, 1020, 1020, 1022, 1023, 1015, 1015, 1015, 1015, 941, 941, 941, 941, 938, 940, 946, 949, 945, 946, 946, 954, 955, 954, 955, 954, 955, 954, 955, 948, 948, 948, 0, 945, 938, 938, 941, 940, 941, 938, 996, 999, 996, 1005, 1005, 1007, 1007, 1007, 1015, 1013, 1008, 1009, 999, 1001, 1001, 1006, 1004, 1007, 1004, 1007, 1006, 1011, 1015, 1011, 1014, 1011, 1019, 1021, 1021, 1020, 1020, 1020, 1020, 1020, 1011, 1012, 1011, 1010, 1011, 1003, 1004, 1004, 1004, 1004, 1005, 1017, 1016, 1017, 1017, 1025, 1025, 0, 1025, 1024, 1024, 1024, 1023, 1025, 1024, 1024, 1025, 1015, 1018, 1016, 1015, 0, 1016, 1007, 1007, 1008, 965, 1026};

//double cuthi[122] ={
//      1129, 1128, 1129, 1129, 1137, 1137, 1137, 0, 1136, 1137, 1137, 1136, 1135, 1135, 1137, 1138, 1130, 1130, 1130, 1130, 1056, 1056, 1056, 1056, 1053, 1055, 1061, 1064, 1060, 1061, 1061, 1069, 1070, 1069, 1070, 1069, 1070, 1069, 1070, 1063, 1063, 1063, 0, 1060, 1053, 1053, 1056, 1055, 1056, 1053, 1111, 1114, 1111, 1120, 1120, 1122, 1122, 1122, 1130, 1128, 1123, 1124, 1114, 1116, 1116, 1121, 1119, 1122, 1119, 1122, 1121, 1126, 1130, 1126, 1129, 1126, 1134, 1136, 1136, 1135, 1135, 1135, 1135, 1135, 1126, 1127, 1126, 1125, 1126, 1118, 1119, 1119, 1119, 1119, 1120, 1132, 1131, 1132, 1132, 1140, 1140, 0, 1140, 1139, 1139, 1139, 1138, 1140, 1139, 1139, 1140, 1130, 1133, 1131, 1130, 0, 1131, 1122, 1122, 1123, 1080, 1141};