/* The alpha energy in the CsI(Tl) for a given channel is given in terms */ /* of the proton energy in MeV for that chan */ /* E_alpha = A[i]*E_prot + B[i] * ln(1.0 + C[i]*E_prot), in MeV */ /* or if you want to use a single A,B,C for all detectors */ /* a global fit to the E_alpha vs E_proton data for all detectors gives the */ /* parameters A = 0.92986, B = 1.2883, C = 2.5926 */ static float A[96] = { 0.934200, 0.942299, 0.936925, 0.959343, 0.937798, 0.928362, 0.921240, 0.912639, 0.901005, 0.908404, 0.940089, 0.912461, 0.912350, 0.901195, 0.914769, 0.892475, 0.896736, 0.899863, 0.902169, 0.933028, 0.922739, 0.928672, 0.899594, 0.939483, 0.937672, 0.920185, 0.906493, 0.892060, 0.904419, 0.935633, 0.911003, 0.923204, 0.942717, 0.953565, 0.952435, 0.951577, 0.942402, 0.949788, 0.954599, 0.952521, 0.915878, 0.928887, 0.912831, 0.922268, 0.898698, 0.893863, 0.904683, 0.911713, 0.904549, 0.905733, 0.912470, 0.909260, 0.919230, 0.919775, 0.996623, 0.994609, 0.993291, 1.010100, 0.998824, 0.989882, 1.003585, 0.000000, 0.993215, 0.995257, 0.996364, 1.003246, 1.012534, 0.990252, 1.004867, 0.999765, 0.994611, 1.017524, 0.993539, 1.011617, 1.003621, 1.030644, 1.012893, 1.027525, 1.012459, 1.011196, 1.012816, 1.027115, 1.045016, 1.016040, 1.042942, 1.019954, 1.063221, 1.025518, 1.036987, 1.015121, 1.055128, 1.034737, 1.042734, 1.025288, 1.027369, 1.040136}; static float B[96] = { 1.263646, 1.244334, 1.344514, 1.261847, 1.256271, 1.422589, 1.409140, 1.273158, 1.454960, 1.343008, 1.357531, 1.428617, 1.367857, 1.354336, 1.356435, 1.444668, 1.230925, 1.343830, 1.344840, 1.248203, 1.256978, 1.337096, 1.258376, 1.239182, 1.218393, 1.235717, 1.250234, 1.365658, 1.246750, 1.302025, 1.301854, 1.288867, 1.289211, 1.268643, 1.313126, 1.270537, 1.318379, 1.286654, 1.200452, 1.305638, 1.279761, 1.271001, 1.305433, 1.289665, 1.416753, 1.367293, 1.311324, 1.288310, 1.265963, 1.248262, 1.246121, 1.262277, 1.244915, 1.313498, 1.005219, 1.017346, 1.013901, 0.995649, 1.003055, 1.040711, 0.959026, 0.000000, 1.011970, 1.013911, 0.980561, 0.978653, 0.925989, 1.027979, 1.109114, 1.045096, 1.104403, 1.031438, 1.062865, 1.006176, 1.082263, 0.997566, 1.040285, 1.047964, 1.091244, 1.042422, 1.087310, 1.022978, 0.997175, 1.050182, 1.000465, 1.077645, 0.985055, 1.045362, 1.047004, 1.088751, 0.958184, 1.002613, 0.977186, 0.999757, 1.002293, 0.993558}; static float C[96] = { 2.772990, 2.296231, 3.044439, 2.484341, 2.753946, 2.933591, 2.860706, 3.162925, 3.581383, 3.121153, 3.217770, 2.774136, 3.222666, 2.964709, 3.263576, 3.262366, 3.026444, 2.924540, 2.921840, 3.218950, 3.136670, 2.795932, 3.215227, 2.588296, 2.273101, 2.383447, 2.860040, 3.042692, 2.549426, 2.431569, 2.457300, 2.478588, 2.464625, 2.441392, 2.305521, 2.442723, 2.310726, 2.333163, 2.612827, 2.273497, 2.507488, 2.616465, 2.798524, 2.725375, 2.491995, 2.868137, 2.618502, 2.764037, 2.694419, 2.660241, 2.715499, 2.762967, 2.728182, 2.445356, 2.303395, 2.414088, 2.395042, 2.318093, 2.392956, 2.581927, 2.307777, 0.000000, 2.529146, 2.419894, 2.345280, 2.234665, 2.186570, 2.376382, 2.630177, 2.481780, 2.725228, 2.471283, 2.533020, 2.361137, 2.603937, 2.386167, 2.509250, 2.602179, 2.813197, 2.492308, 2.619832, 2.494498, 2.449003, 2.517344, 2.452461, 2.587436, 2.441491, 2.536896, 2.528284, 2.622791, 2.402187, 2.450238, 2.412056, 2.410418, 2.383656, 2.435109};