#ifndef __CLAS_G12_ECOR_HPP__
#define __CLAS_G12_ECOR_HPP__

#include <fstream>
#include <iostream>
#include <map>
#include <sstream>
#include <stdexcept>
#include <string>
#include <vector>

namespace clas
{
namespace g12 {


double corrected_beam_energy(int run, double uncorrected_beam_energy)
{
double ret = 0;
    if(run == 56363){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997942;
    ret = E_gamma_new;
    }
    else if(run == 56365){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996246;
    ret = E_gamma_new;
    }
    else if(run == 56369){
    double E_scattered = 5.68436 - uncorrected_beam_energy;
    double E_gamma_new = 5.68436 - E_scattered*0.995988;
    ret = E_gamma_new;
    }
    else if(run == 56384){
    double E_scattered = 5.71141 - uncorrected_beam_energy;
    double E_gamma_new = 5.71141 - E_scattered*0.996996;
    ret = E_gamma_new;
    }
    else if(run == 56386){
    double E_scattered = 5.71048 - uncorrected_beam_energy;
    double E_gamma_new = 5.71048 - E_scattered*0.996971;
    ret = E_gamma_new;
    }
    else if(run == 56401){
    double E_scattered = 5.7416 - uncorrected_beam_energy;
    double E_gamma_new = 5.7416 - E_scattered*0.997435;
    ret = E_gamma_new;
    }
    else if(run == 56403){
    double E_scattered = 5.74294 - uncorrected_beam_energy;
    double E_gamma_new = 5.74294 - E_scattered*0.997747;
    ret = E_gamma_new;
    }
    else if(run == 56404){
    double E_scattered = 5.75506 - uncorrected_beam_energy;
    double E_gamma_new = 5.75506 - E_scattered*0.998692;
    ret = E_gamma_new;
    }
    else if(run == 56405){
    double E_scattered = 5.7574 - uncorrected_beam_energy;
    double E_gamma_new = 5.7574 - E_scattered*0.998753;
    ret = E_gamma_new;
    }
    else if(run == 56406){
    double E_scattered = 5.71479 - uncorrected_beam_energy;
    double E_gamma_new = 5.71479 - E_scattered*0.996315;
    ret = E_gamma_new;
    }
    else if(run == 56408){
    double E_scattered = 5.71479 - uncorrected_beam_energy;
    double E_gamma_new = 5.71479 - E_scattered*0.995965;
    ret = E_gamma_new;
    }
    else if(run == 56410){
    double E_scattered = 5.7246 - uncorrected_beam_energy;
    double E_gamma_new = 5.7246 - E_scattered*0.996965;
    ret = E_gamma_new;
    }
    else if(run == 56420){
    double E_scattered = 5.7246 - uncorrected_beam_energy;
    double E_gamma_new = 5.7246 - E_scattered*0.997612;
    ret = E_gamma_new;
    }
    else if(run == 56421){
    double E_scattered = 5.7246 - uncorrected_beam_energy;
    double E_gamma_new = 5.7246 - E_scattered*0.996465;
    ret = E_gamma_new;
    }
    else if(run == 56422){
    double E_scattered = 5.74235 - uncorrected_beam_energy;
    double E_gamma_new = 5.74235 - E_scattered*0.997511;
    ret = E_gamma_new;
    }
    else if(run == 56435){
    double E_scattered = 5.71543 - uncorrected_beam_energy;
    double E_gamma_new = 5.71543 - E_scattered*0.995574;
    ret = E_gamma_new;
    }
    else if(run == 56436){
    double E_scattered = 5.7162 - uncorrected_beam_energy;
    double E_gamma_new = 5.7162 - E_scattered*0.996023;
    ret = E_gamma_new;
    }
    else if(run == 56441){
    double E_scattered = 5.71571 - uncorrected_beam_energy;
    double E_gamma_new = 5.71571 - E_scattered*0.99651;
    ret = E_gamma_new;
    }
    else if(run == 56442){
    double E_scattered = 5.71611 - uncorrected_beam_energy;
    double E_gamma_new = 5.71611 - E_scattered*0.996479;
    ret = E_gamma_new;
    }
    else if(run == 56443){
    double E_scattered = 5.71548 - uncorrected_beam_energy;
    double E_gamma_new = 5.71548 - E_scattered*0.99648;
    ret = E_gamma_new;
    }
    else if(run == 56445){
    double E_scattered = 5.71578 - uncorrected_beam_energy;
    double E_gamma_new = 5.71578 - E_scattered*0.996491;
    ret = E_gamma_new;
    }
    else if(run == 56446){
    double E_scattered = 5.71578 - uncorrected_beam_energy;
    double E_gamma_new = 5.71578 - E_scattered*0.994837;
    ret = E_gamma_new;
    }
    else if(run == 56447){
    double E_scattered = 5.71583 - uncorrected_beam_energy;
    double E_gamma_new = 5.71583 - E_scattered*0.996164;
    ret = E_gamma_new;
    }
    else if(run == 56448){
    double E_scattered = 5.7156 - uncorrected_beam_energy;
    double E_gamma_new = 5.7156 - E_scattered*0.995801;
    ret = E_gamma_new;
    }
    else if(run == 56449){
    double E_scattered = 5.71581 - uncorrected_beam_energy;
    double E_gamma_new = 5.71581 - E_scattered*0.995676;
    ret = E_gamma_new;
    }
    else if(run == 56450){
    double E_scattered = 5.71733 - uncorrected_beam_energy;
    double E_gamma_new = 5.71733 - E_scattered*0.995902;
    ret = E_gamma_new;
    }
    else if(run == 56453){
    double E_scattered = 5.717 - uncorrected_beam_energy;
    double E_gamma_new = 5.717 - E_scattered*0.99608;
    ret = E_gamma_new;
    }
    else if(run == 56454){
    double E_scattered = 5.7173 - uncorrected_beam_energy;
    double E_gamma_new = 5.7173 - E_scattered*0.995104;
    ret = E_gamma_new;
    }
    else if(run == 56455){
    double E_scattered = 5.71709 - uncorrected_beam_energy;
    double E_gamma_new = 5.71709 - E_scattered*0.995819;
    ret = E_gamma_new;
    }
    else if(run == 56456){
    double E_scattered = 5.7174 - uncorrected_beam_energy;
    double E_gamma_new = 5.7174 - E_scattered*0.994909;
    ret = E_gamma_new;
    }
    else if(run == 56457){
    double E_scattered = 5.7174 - uncorrected_beam_energy;
    double E_gamma_new = 5.7174 - E_scattered*0.995918;
    ret = E_gamma_new;
    }
    else if(run == 56458){
    double E_scattered = 5.72061 - uncorrected_beam_energy;
    double E_gamma_new = 5.72061 - E_scattered*0.996474;
    ret = E_gamma_new;
    }
    else if(run == 56459){
    double E_scattered = 5.72226 - uncorrected_beam_energy;
    double E_gamma_new = 5.72226 - E_scattered*0.996352;
    ret = E_gamma_new;
    }
    else if(run == 56460){
    double E_scattered = 5.71832 - uncorrected_beam_energy;
    double E_gamma_new = 5.71832 - E_scattered*0.996672;
    ret = E_gamma_new;
    }
    else if(run == 56461){
    double E_scattered = 5.71731 - uncorrected_beam_energy;
    double E_gamma_new = 5.71731 - E_scattered*0.996044;
    ret = E_gamma_new;
    }
    else if(run == 56462){
    double E_scattered = 5.71653 - uncorrected_beam_energy;
    double E_gamma_new = 5.71653 - E_scattered*0.995326;
    ret = E_gamma_new;
    }
    else if(run == 56465){
    double E_scattered = 5.71616 - uncorrected_beam_energy;
    double E_gamma_new = 5.71616 - E_scattered*0.996238;
    ret = E_gamma_new;
    }
    else if(run == 56467){
    double E_scattered = 5.7174 - uncorrected_beam_energy;
    double E_gamma_new = 5.7174 - E_scattered*0.995519;
    ret = E_gamma_new;
    }
    else if(run == 56468){
    double E_scattered = 5.71741 - uncorrected_beam_energy;
    double E_gamma_new = 5.71741 - E_scattered*0.994413;
    ret = E_gamma_new;
    }
    else if(run == 56469){
    double E_scattered = 5.71713 - uncorrected_beam_energy;
    double E_gamma_new = 5.71713 - E_scattered*0.995584;
    ret = E_gamma_new;
    }
    else if(run == 56470){
    double E_scattered = 5.71696 - uncorrected_beam_energy;
    double E_gamma_new = 5.71696 - E_scattered*0.995532;
    ret = E_gamma_new;
    }
    else if(run == 56471){
    double E_scattered = 5.71668 - uncorrected_beam_energy;
    double E_gamma_new = 5.71668 - E_scattered*0.996027;
    ret = E_gamma_new;
    }
    else if(run == 56472){
    double E_scattered = 5.71694 - uncorrected_beam_energy;
    double E_gamma_new = 5.71694 - E_scattered*0.99621;
    ret = E_gamma_new;
    }
    else if(run == 56476){
    double E_scattered = 5.71303 - uncorrected_beam_energy;
    double E_gamma_new = 5.71303 - E_scattered*0.997191;
    ret = E_gamma_new;
    }
    else if(run == 56478){
    double E_scattered = 5.71227 - uncorrected_beam_energy;
    double E_gamma_new = 5.71227 - E_scattered*0.997087;
    ret = E_gamma_new;
    }
    else if(run == 56479){
    double E_scattered = 5.71298 - uncorrected_beam_energy;
    double E_gamma_new = 5.71298 - E_scattered*0.997148;
    ret = E_gamma_new;
    }
    else if(run == 56480){
    double E_scattered = 5.71452 - uncorrected_beam_energy;
    double E_gamma_new = 5.71452 - E_scattered*0.997387;
    ret = E_gamma_new;
    }
    else if(run == 56481){
    double E_scattered = 5.71452 - uncorrected_beam_energy;
    double E_gamma_new = 5.71452 - E_scattered*0.996949;
    ret = E_gamma_new;
    }
    else if(run == 56482){
    double E_scattered = 5.71209 - uncorrected_beam_energy;
    double E_gamma_new = 5.71209 - E_scattered*0.996923;
    ret = E_gamma_new;
    }
    else if(run == 56483){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.99693;
    ret = E_gamma_new;
    }
    else if(run == 56485){
    double E_scattered = 5.71223 - uncorrected_beam_energy;
    double E_gamma_new = 5.71223 - E_scattered*0.997183;
    ret = E_gamma_new;
    }
    else if(run == 56486){
    double E_scattered = 5.71223 - uncorrected_beam_energy;
    double E_gamma_new = 5.71223 - E_scattered*0.996961;
    ret = E_gamma_new;
    }
    else if(run == 56487){
    double E_scattered = 5.71279 - uncorrected_beam_energy;
    double E_gamma_new = 5.71279 - E_scattered*0.997184;
    ret = E_gamma_new;
    }
    else if(run == 56489){
    double E_scattered = 5.71176 - uncorrected_beam_energy;
    double E_gamma_new = 5.71176 - E_scattered*0.997078;
    ret = E_gamma_new;
    }
    else if(run == 56490){
    double E_scattered = 5.71175 - uncorrected_beam_energy;
    double E_gamma_new = 5.71175 - E_scattered*0.996366;
    ret = E_gamma_new;
    }
    else if(run == 56499){
    double E_scattered = 5.7116 - uncorrected_beam_energy;
    double E_gamma_new = 5.7116 - E_scattered*0.997012;
    ret = E_gamma_new;
    }
    else if(run == 56501){
    double E_scattered = 5.71118 - uncorrected_beam_energy;
    double E_gamma_new = 5.71118 - E_scattered*0.997096;
    ret = E_gamma_new;
    }
    else if(run == 56502){
    double E_scattered = 5.71458 - uncorrected_beam_energy;
    double E_gamma_new = 5.71458 - E_scattered*0.996106;
    ret = E_gamma_new;
    }
    else if(run == 56503){
    double E_scattered = 5.71468 - uncorrected_beam_energy;
    double E_gamma_new = 5.71468 - E_scattered*0.996174;
    ret = E_gamma_new;
    }
    else if(run == 56504){
    double E_scattered = 5.71475 - uncorrected_beam_energy;
    double E_gamma_new = 5.71475 - E_scattered*0.996286;
    ret = E_gamma_new;
    }
    else if(run == 56505){
    double E_scattered = 5.7153 - uncorrected_beam_energy;
    double E_gamma_new = 5.7153 - E_scattered*0.995905;
    ret = E_gamma_new;
    }
    else if(run == 56506){
    double E_scattered = 5.7153 - uncorrected_beam_energy;
    double E_gamma_new = 5.7153 - E_scattered*0.996172;
    ret = E_gamma_new;
    }
    else if(run == 56508){
    double E_scattered = 5.7153 - uncorrected_beam_energy;
    double E_gamma_new = 5.7153 - E_scattered*0.996752;
    ret = E_gamma_new;
    }
    else if(run == 56509){
    double E_scattered = 5.71333 - uncorrected_beam_energy;
    double E_gamma_new = 5.71333 - E_scattered*0.996561;
    ret = E_gamma_new;
    }
    else if(run == 56510){
    double E_scattered = 5.71455 - uncorrected_beam_energy;
    double E_gamma_new = 5.71455 - E_scattered*0.99702;
    ret = E_gamma_new;
    }
    else if(run == 56513){
    double E_scattered = 5.71373 - uncorrected_beam_energy;
    double E_gamma_new = 5.71373 - E_scattered*0.996505;
    ret = E_gamma_new;
    }
    else if(run == 56514){
    double E_scattered = 5.71445 - uncorrected_beam_energy;
    double E_gamma_new = 5.71445 - E_scattered*0.996591;
    ret = E_gamma_new;
    }
    else if(run == 56515){
    double E_scattered = 5.71283 - uncorrected_beam_energy;
    double E_gamma_new = 5.71283 - E_scattered*0.996658;
    ret = E_gamma_new;
    }
    else if(run == 56516){
    double E_scattered = 5.71357 - uncorrected_beam_energy;
    double E_gamma_new = 5.71357 - E_scattered*0.99621;
    ret = E_gamma_new;
    }
    else if(run == 56517){
    double E_scattered = 5.71409 - uncorrected_beam_energy;
    double E_gamma_new = 5.71409 - E_scattered*0.996724;
    ret = E_gamma_new;
    }
    else if(run == 56519){
    double E_scattered = 5.71595 - uncorrected_beam_energy;
    double E_gamma_new = 5.71595 - E_scattered*0.995108;
    ret = E_gamma_new;
    }
    else if(run == 56520){
    double E_scattered = 5.71355 - uncorrected_beam_energy;
    double E_gamma_new = 5.71355 - E_scattered*0.996536;
    ret = E_gamma_new;
    }
    else if(run == 56521){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996239;
    ret = E_gamma_new;
    }
    else if(run == 56522){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996639;
    ret = E_gamma_new;
    }
    else if(run == 56523){
    double E_scattered = 5.71392 - uncorrected_beam_energy;
    double E_gamma_new = 5.71392 - E_scattered*0.996428;
    ret = E_gamma_new;
    }
    else if(run == 56524){
    double E_scattered = 5.71466 - uncorrected_beam_energy;
    double E_gamma_new = 5.71466 - E_scattered*0.996483;
    ret = E_gamma_new;
    }
    else if(run == 56525){
    double E_scattered = 5.7133 - uncorrected_beam_energy;
    double E_gamma_new = 5.7133 - E_scattered*0.997143;
    ret = E_gamma_new;
    }
    else if(run == 56526){
    double E_scattered = 5.71402 - uncorrected_beam_energy;
    double E_gamma_new = 5.71402 - E_scattered*0.996094;
    ret = E_gamma_new;
    }
    else if(run == 56527){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996361;
    ret = E_gamma_new;
    }
    else if(run == 56528){
    double E_scattered = 5.71402 - uncorrected_beam_energy;
    double E_gamma_new = 5.71402 - E_scattered*0.996597;
    ret = E_gamma_new;
    }
    else if(run == 56529){
    double E_scattered = 5.7146 - uncorrected_beam_energy;
    double E_gamma_new = 5.7146 - E_scattered*0.996366;
    ret = E_gamma_new;
    }
    else if(run == 56530){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996524;
    ret = E_gamma_new;
    }
    else if(run == 56531){
    double E_scattered = 5.71404 - uncorrected_beam_energy;
    double E_gamma_new = 5.71404 - E_scattered*0.996477;
    ret = E_gamma_new;
    }
    else if(run == 56532){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997222;
    ret = E_gamma_new;
    }
    else if(run == 56533){
    double E_scattered = 5.7139 - uncorrected_beam_energy;
    double E_gamma_new = 5.7139 - E_scattered*0.996217;
    ret = E_gamma_new;
    }
    else if(run == 56534){
    double E_scattered = 5.71387 - uncorrected_beam_energy;
    double E_gamma_new = 5.71387 - E_scattered*0.996164;
    ret = E_gamma_new;
    }
    else if(run == 56535){
    double E_scattered = 5.71306 - uncorrected_beam_energy;
    double E_gamma_new = 5.71306 - E_scattered*0.995909;
    ret = E_gamma_new;
    }
    else if(run == 56536){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996113;
    ret = E_gamma_new;
    }
    else if(run == 56537){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996566;
    ret = E_gamma_new;
    }
    else if(run == 56538){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996402;
    ret = E_gamma_new;
    }
    else if(run == 56539){
    double E_scattered = 5.7135 - uncorrected_beam_energy;
    double E_gamma_new = 5.7135 - E_scattered*0.996222;
    ret = E_gamma_new;
    }
    else if(run == 56540){
    double E_scattered = 5.71758 - uncorrected_beam_energy;
    double E_gamma_new = 5.71758 - E_scattered*0.997044;
    ret = E_gamma_new;
    }
    else if(run == 56541){
    double E_scattered = 5.71506 - uncorrected_beam_energy;
    double E_gamma_new = 5.71506 - E_scattered*0.996206;
    ret = E_gamma_new;
    }
    else if(run == 56542){
    double E_scattered = 5.71408 - uncorrected_beam_energy;
    double E_gamma_new = 5.71408 - E_scattered*0.996108;
    ret = E_gamma_new;
    }
    else if(run == 56544){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.995668;
    ret = E_gamma_new;
    }
    else if(run == 56545){
    double E_scattered = 5.71267 - uncorrected_beam_energy;
    double E_gamma_new = 5.71267 - E_scattered*0.996235;
    ret = E_gamma_new;
    }
    else if(run == 56546){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.996336;
    ret = E_gamma_new;
    }
    else if(run == 56547){
    double E_scattered = 5.71393 - uncorrected_beam_energy;
    double E_gamma_new = 5.71393 - E_scattered*0.996144;
    ret = E_gamma_new;
    }
    else if(run == 56548){
    double E_scattered = 5.71354 - uncorrected_beam_energy;
    double E_gamma_new = 5.71354 - E_scattered*0.996294;
    ret = E_gamma_new;
    }
    else if(run == 56549){
    double E_scattered = 5.71288 - uncorrected_beam_energy;
    double E_gamma_new = 5.71288 - E_scattered*0.996271;
    ret = E_gamma_new;
    }
    else if(run == 56550){
    double E_scattered = 5.71289 - uncorrected_beam_energy;
    double E_gamma_new = 5.71289 - E_scattered*0.996369;
    ret = E_gamma_new;
    }
    else if(run == 56555){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996216;
    ret = E_gamma_new;
    }
    else if(run == 56556){
    double E_scattered = 5.71383 - uncorrected_beam_energy;
    double E_gamma_new = 5.71383 - E_scattered*0.996713;
    ret = E_gamma_new;
    }
    else if(run == 56559){
    double E_scattered = 5.7155 - uncorrected_beam_energy;
    double E_gamma_new = 5.7155 - E_scattered*0.996437;
    ret = E_gamma_new;
    }
    else if(run == 56561){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.995985;
    ret = E_gamma_new;
    }
    else if(run == 56562){
    double E_scattered = 5.7123 - uncorrected_beam_energy;
    double E_gamma_new = 5.7123 - E_scattered*0.997588;
    ret = E_gamma_new;
    }
    else if(run == 56563){
    double E_scattered = 5.71287 - uncorrected_beam_energy;
    double E_gamma_new = 5.71287 - E_scattered*0.997283;
    ret = E_gamma_new;
    }
    else if(run == 56564){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997457;
    ret = E_gamma_new;
    }
    else if(run == 56573){
    double E_scattered = 5.71415 - uncorrected_beam_energy;
    double E_gamma_new = 5.71415 - E_scattered*0.997052;
    ret = E_gamma_new;
    }
    else if(run == 56574){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997439;
    ret = E_gamma_new;
    }
    else if(run == 56575){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997097;
    ret = E_gamma_new;
    }
    else if(run == 56576){
    double E_scattered = 5.71392 - uncorrected_beam_energy;
    double E_gamma_new = 5.71392 - E_scattered*0.99705;
    ret = E_gamma_new;
    }
    else if(run == 56577){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996992;
    ret = E_gamma_new;
    }
    else if(run == 56578){
    double E_scattered = 5.71381 - uncorrected_beam_energy;
    double E_gamma_new = 5.71381 - E_scattered*0.99709;
    ret = E_gamma_new;
    }
    else if(run == 56579){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997416;
    ret = E_gamma_new;
    }
    else if(run == 56580){
    double E_scattered = 5.712 - uncorrected_beam_energy;
    double E_gamma_new = 5.712 - E_scattered*0.997735;
    ret = E_gamma_new;
    }
    else if(run == 56581){
    double E_scattered = 5.7141 - uncorrected_beam_energy;
    double E_gamma_new = 5.7141 - E_scattered*0.997062;
    ret = E_gamma_new;
    }
    else if(run == 56582){
    double E_scattered = 5.71345 - uncorrected_beam_energy;
    double E_gamma_new = 5.71345 - E_scattered*0.997156;
    ret = E_gamma_new;
    }
    else if(run == 56583){
    double E_scattered = 5.71397 - uncorrected_beam_energy;
    double E_gamma_new = 5.71397 - E_scattered*0.997138;
    ret = E_gamma_new;
    }
    else if(run == 56585){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997118;
    ret = E_gamma_new;
    }
    else if(run == 56586){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997148;
    ret = E_gamma_new;
    }
    else if(run == 56587){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997232;
    ret = E_gamma_new;
    }
    else if(run == 56588){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997114;
    ret = E_gamma_new;
    }
    else if(run == 56589){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997007;
    ret = E_gamma_new;
    }
    else if(run == 56590){
    double E_scattered = 5.71444 - uncorrected_beam_energy;
    double E_gamma_new = 5.71444 - E_scattered*0.997109;
    ret = E_gamma_new;
    }
    else if(run == 56591){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997296;
    ret = E_gamma_new;
    }
    else if(run == 56592){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997054;
    ret = E_gamma_new;
    }
    else if(run == 56593){
    double E_scattered = 5.71403 - uncorrected_beam_energy;
    double E_gamma_new = 5.71403 - E_scattered*0.996945;
    ret = E_gamma_new;
    }
    else if(run == 56605){
    double E_scattered = 5.71399 - uncorrected_beam_energy;
    double E_gamma_new = 5.71399 - E_scattered*0.997163;
    ret = E_gamma_new;
    }
    else if(run == 56608){
    double E_scattered = 5.71423 - uncorrected_beam_energy;
    double E_gamma_new = 5.71423 - E_scattered*0.997385;
    ret = E_gamma_new;
    }
    else if(run == 56609){
    double E_scattered = 5.71397 - uncorrected_beam_energy;
    double E_gamma_new = 5.71397 - E_scattered*0.997136;
    ret = E_gamma_new;
    }
    else if(run == 56610){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997215;
    ret = E_gamma_new;
    }
    else if(run == 56611){
    double E_scattered = 5.71407 - uncorrected_beam_energy;
    double E_gamma_new = 5.71407 - E_scattered*0.997199;
    ret = E_gamma_new;
    }
    else if(run == 56612){
    double E_scattered = 5.71478 - uncorrected_beam_energy;
    double E_gamma_new = 5.71478 - E_scattered*0.997152;
    ret = E_gamma_new;
    }
    else if(run == 56614){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997203;
    ret = E_gamma_new;
    }
    else if(run == 56615){
    double E_scattered = 5.71456 - uncorrected_beam_energy;
    double E_gamma_new = 5.71456 - E_scattered*0.997159;
    ret = E_gamma_new;
    }
    else if(run == 56616){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997353;
    ret = E_gamma_new;
    }
    else if(run == 56617){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997199;
    ret = E_gamma_new;
    }
    else if(run == 56618){
    double E_scattered = 5.71396 - uncorrected_beam_energy;
    double E_gamma_new = 5.71396 - E_scattered*0.9972;
    ret = E_gamma_new;
    }
    else if(run == 56619){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996857;
    ret = E_gamma_new;
    }
    else if(run == 56620){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.996932;
    ret = E_gamma_new;
    }
    else if(run == 56621){
    double E_scattered = 5.71399 - uncorrected_beam_energy;
    double E_gamma_new = 5.71399 - E_scattered*0.997046;
    ret = E_gamma_new;
    }
    else if(run == 56622){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997042;
    ret = E_gamma_new;
    }
    else if(run == 56623){
    double E_scattered = 5.71395 - uncorrected_beam_energy;
    double E_gamma_new = 5.71395 - E_scattered*0.997285;
    ret = E_gamma_new;
    }
    else if(run == 56624){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997294;
    ret = E_gamma_new;
    }
    else if(run == 56625){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997303;
    ret = E_gamma_new;
    }
    else if(run == 56626){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.997183;
    ret = E_gamma_new;
    }
    else if(run == 56627){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997134;
    ret = E_gamma_new;
    }
    else if(run == 56628){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997133;
    ret = E_gamma_new;
    }
    else if(run == 56630){
    double E_scattered = 5.71399 - uncorrected_beam_energy;
    double E_gamma_new = 5.71399 - E_scattered*0.997211;
    ret = E_gamma_new;
    }
    else if(run == 56631){
    double E_scattered = 5.71537 - uncorrected_beam_energy;
    double E_gamma_new = 5.71537 - E_scattered*0.997332;
    ret = E_gamma_new;
    }
    else if(run == 56632){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997242;
    ret = E_gamma_new;
    }
    else if(run == 56633){
    double E_scattered = 5.71483 - uncorrected_beam_energy;
    double E_gamma_new = 5.71483 - E_scattered*0.997443;
    ret = E_gamma_new;
    }
    else if(run == 56634){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997315;
    ret = E_gamma_new;
    }
    else if(run == 56635){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997968;
    ret = E_gamma_new;
    }
    else if(run == 56636){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.997068;
    ret = E_gamma_new;
    }
    else if(run == 56637){
    double E_scattered = 5.7136 - uncorrected_beam_energy;
    double E_gamma_new = 5.7136 - E_scattered*0.998167;
    ret = E_gamma_new;
    }
    else if(run == 56638){
    double E_scattered = 5.71398 - uncorrected_beam_energy;
    double E_gamma_new = 5.71398 - E_scattered*0.99707;
    ret = E_gamma_new;
    }
    else if(run == 56639){
    double E_scattered = 5.714 - uncorrected_beam_energy;
    double E_gamma_new = 5.714 - E_scattered*0.998179;
    ret = E_gamma_new;
    }
    else if(run == 56640){
    double E_scattered = 5.71396 - uncorrected_beam_energy;
    double E_gamma_new = 5.71396 - E_scattered*0.99725;
    ret = E_gamma_new;
    }
    else if(run == 56641){
    double E_scattered = 5.71392 - uncorrected_beam_energy;
    double E_gamma_new = 5.71392 - E_scattered*0.997211;
    ret = E_gamma_new;
    }
    else if(run == 56642){
    double E_scattered = 5.71392 - uncorrected_beam_energy;
    double E_gamma_new = 5.71392 - E_scattered*0.997074;
    ret = E_gamma_new;
    }
    else if(run == 56643){
    double E_scattered = 5.71431 - uncorrected_beam_energy;
    double E_gamma_new = 5.71431 - E_scattered*0.997104;
    ret = E_gamma_new;
    }
    else if(run == 56644){
    double E_scattered = 5.71036 - uncorrected_beam_energy;
    double E_gamma_new = 5.71036 - E_scattered*0.99804;
    ret = E_gamma_new;
    }
    else if(run == 56646){
    double E_scattered = 5.71241 - uncorrected_beam_energy;
    double E_gamma_new = 5.71241 - E_scattered*0.997715;
    ret = E_gamma_new;
    }
    else if(run == 56653){
    double E_scattered = 5.71273 - uncorrected_beam_energy;
    double E_gamma_new = 5.71273 - E_scattered*0.997759;
    ret = E_gamma_new;
    }
    else if(run == 56654){
    double E_scattered = 5.71235 - uncorrected_beam_energy;
    double E_gamma_new = 5.71235 - E_scattered*0.997982;
    ret = E_gamma_new;
    }
    else if(run == 56655){
    double E_scattered = 5.71244 - uncorrected_beam_energy;
    double E_gamma_new = 5.71244 - E_scattered*0.997797;
    ret = E_gamma_new;
    }
    else if(run == 56656){
    double E_scattered = 5.71219 - uncorrected_beam_energy;
    double E_gamma_new = 5.71219 - E_scattered*0.997765;
    ret = E_gamma_new;
    }
    else if(run == 56660){
    double E_scattered = 5.7132 - uncorrected_beam_energy;
    double E_gamma_new = 5.7132 - E_scattered*0.997662;
    ret = E_gamma_new;
    }
    else if(run == 56661){
    double E_scattered = 5.71293 - uncorrected_beam_energy;
    double E_gamma_new = 5.71293 - E_scattered*0.997637;
    ret = E_gamma_new;
    }
    else if(run == 56663){
    double E_scattered = 5.71313 - uncorrected_beam_energy;
    double E_gamma_new = 5.71313 - E_scattered*0.997514;
    ret = E_gamma_new;
    }
    else if(run == 56664){
    double E_scattered = 5.71317 - uncorrected_beam_energy;
    double E_gamma_new = 5.71317 - E_scattered*0.997525;
    ret = E_gamma_new;
    }
    else if(run == 56665){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997574;
    ret = E_gamma_new;
    }
    else if(run == 56666){
    double E_scattered = 5.7134 - uncorrected_beam_energy;
    double E_gamma_new = 5.7134 - E_scattered*0.997686;
    ret = E_gamma_new;
    }
    else if(run == 56667){
    double E_scattered = 5.71317 - uncorrected_beam_energy;
    double E_gamma_new = 5.71317 - E_scattered*0.997717;
    ret = E_gamma_new;
    }
    else if(run == 56668){
    double E_scattered = 5.71211 - uncorrected_beam_energy;
    double E_gamma_new = 5.71211 - E_scattered*0.997563;
    ret = E_gamma_new;
    }
    else if(run == 56669){
    double E_scattered = 5.71282 - uncorrected_beam_energy;
    double E_gamma_new = 5.71282 - E_scattered*0.997323;
    ret = E_gamma_new;
    }
    else if(run == 56670){
    double E_scattered = 5.71276 - uncorrected_beam_energy;
    double E_gamma_new = 5.71276 - E_scattered*0.99773;
    ret = E_gamma_new;
    }
    else if(run == 56673){
    double E_scattered = 5.71268 - uncorrected_beam_energy;
    double E_gamma_new = 5.71268 - E_scattered*0.997682;
    ret = E_gamma_new;
    }
    else if(run == 56674){
    double E_scattered = 5.71262 - uncorrected_beam_energy;
    double E_gamma_new = 5.71262 - E_scattered*0.997702;
    ret = E_gamma_new;
    }
    else if(run == 56675){
    double E_scattered = 5.7126 - uncorrected_beam_energy;
    double E_gamma_new = 5.7126 - E_scattered*0.997773;
    ret = E_gamma_new;
    }
    else if(run == 56679){
    double E_scattered = 5.71272 - uncorrected_beam_energy;
    double E_gamma_new = 5.71272 - E_scattered*0.997965;
    ret = E_gamma_new;
    }
    else if(run == 56680){
    double E_scattered = 5.71341 - uncorrected_beam_energy;
    double E_gamma_new = 5.71341 - E_scattered*0.997688;
    ret = E_gamma_new;
    }
    else if(run == 56681){
    double E_scattered = 5.71283 - uncorrected_beam_energy;
    double E_gamma_new = 5.71283 - E_scattered*0.997597;
    ret = E_gamma_new;
    }
    else if(run == 56683){
    double E_scattered = 5.71248 - uncorrected_beam_energy;
    double E_gamma_new = 5.71248 - E_scattered*0.997575;
    ret = E_gamma_new;
    }
    else if(run == 56685){
    double E_scattered = 5.7124 - uncorrected_beam_energy;
    double E_gamma_new = 5.7124 - E_scattered*0.997317;
    ret = E_gamma_new;
    }
    else if(run == 56686){
    double E_scattered = 5.71355 - uncorrected_beam_energy;
    double E_gamma_new = 5.71355 - E_scattered*0.997274;
    ret = E_gamma_new;
    }
    else if(run == 56687){
    double E_scattered = 5.71395 - uncorrected_beam_energy;
    double E_gamma_new = 5.71395 - E_scattered*0.997319;
    ret = E_gamma_new;
    }
    else if(run == 56688){
    double E_scattered = 5.71392 - uncorrected_beam_energy;
    double E_gamma_new = 5.71392 - E_scattered*0.99726;
    ret = E_gamma_new;
    }
    else if(run == 56689){
    double E_scattered = 5.71391 - uncorrected_beam_energy;
    double E_gamma_new = 5.71391 - E_scattered*0.997332;
    ret = E_gamma_new;
    }
    else if(run == 56690){
    double E_scattered = 5.71367 - uncorrected_beam_energy;
    double E_gamma_new = 5.71367 - E_scattered*0.997223;
    ret = E_gamma_new;
    }
    else if(run == 56691){
    double E_scattered = 5.71313 - uncorrected_beam_energy;
    double E_gamma_new = 5.71313 - E_scattered*0.997265;
    ret = E_gamma_new;
    }
    else if(run == 56692){
    double E_scattered = 5.71362 - uncorrected_beam_energy;
    double E_gamma_new = 5.71362 - E_scattered*0.997414;
    ret = E_gamma_new;
    }
    else if(run == 56693){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997211;
    ret = E_gamma_new;
    }
    else if(run == 56694){
    double E_scattered = 5.71332 - uncorrected_beam_energy;
    double E_gamma_new = 5.71332 - E_scattered*0.997395;
    ret = E_gamma_new;
    }
    else if(run == 56695){
    double E_scattered = 5.71377 - uncorrected_beam_energy;
    double E_gamma_new = 5.71377 - E_scattered*0.997391;
    ret = E_gamma_new;
    }
    else if(run == 56696){
    double E_scattered = 5.71368 - uncorrected_beam_energy;
    double E_gamma_new = 5.71368 - E_scattered*0.997241;
    ret = E_gamma_new;
    }
    else if(run == 56697){
    double E_scattered = 5.7137 - uncorrected_beam_energy;
    double E_gamma_new = 5.7137 - E_scattered*0.997275;
    ret = E_gamma_new;
    }
    else if(run == 56700){
    double E_scattered = 5.71371 - uncorrected_beam_energy;
    double E_gamma_new = 5.71371 - E_scattered*0.99723;
    ret = E_gamma_new;
    }
    else if(run == 56701){
    double E_scattered = 5.71383 - uncorrected_beam_energy;
    double E_gamma_new = 5.71383 - E_scattered*0.997336;
    ret = E_gamma_new;
    }
    else if(run == 56702){
    double E_scattered = 5.71379 - uncorrected_beam_energy;
    double E_gamma_new = 5.71379 - E_scattered*0.997426;
    ret = E_gamma_new;
    }
    else if(run == 56703){
    double E_scattered = 5.71374 - uncorrected_beam_energy;
    double E_gamma_new = 5.71374 - E_scattered*0.99738;
    ret = E_gamma_new;
    }
    else if(run == 56704){
    double E_scattered = 5.71384 - uncorrected_beam_energy;
    double E_gamma_new = 5.71384 - E_scattered*0.997152;
    ret = E_gamma_new;
    }
    else if(run == 56705){
    double E_scattered = 5.71388 - uncorrected_beam_energy;
    double E_gamma_new = 5.71388 - E_scattered*0.997201;
    ret = E_gamma_new;
    }
    else if(run == 56706){
    double E_scattered = 5.71385 - uncorrected_beam_energy;
    double E_gamma_new = 5.71385 - E_scattered*0.997148;
    ret = E_gamma_new;
    }
    else if(run == 56707){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997391;
    ret = E_gamma_new;
    }
    else if(run == 56708){
    double E_scattered = 5.71378 - uncorrected_beam_energy;
    double E_gamma_new = 5.71378 - E_scattered*0.997093;
    ret = E_gamma_new;
    }
    else if(run == 56710){
    double E_scattered = 5.71378 - uncorrected_beam_energy;
    double E_gamma_new = 5.71378 - E_scattered*0.997281;
    ret = E_gamma_new;
    }
    else if(run == 56711){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997427;
    ret = E_gamma_new;
    }
    else if(run == 56712){
    double E_scattered = 5.71369 - uncorrected_beam_energy;
    double E_gamma_new = 5.71369 - E_scattered*0.997253;
    ret = E_gamma_new;
    }
    else if(run == 56713){
    double E_scattered = 5.71383 - uncorrected_beam_energy;
    double E_gamma_new = 5.71383 - E_scattered*0.997404;
    ret = E_gamma_new;
    }
    else if(run == 56714){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997487;
    ret = E_gamma_new;
    }
    else if(run == 56715){
    double E_scattered = 5.71383 - uncorrected_beam_energy;
    double E_gamma_new = 5.71383 - E_scattered*0.997171;
    ret = E_gamma_new;
    }
    else if(run == 56716){
    double E_scattered = 5.71382 - uncorrected_beam_energy;
    double E_gamma_new = 5.71382 - E_scattered*0.997283;
    ret = E_gamma_new;
    }
    else if(run == 56717){
    double E_scattered = 5.71378 - uncorrected_beam_energy;
    double E_gamma_new = 5.71378 - E_scattered*0.997246;
    ret = E_gamma_new;
    }
    else if(run == 56718){
    double E_scattered = 5.71363 - uncorrected_beam_energy;
    double E_gamma_new = 5.71363 - E_scattered*0.997231;
    ret = E_gamma_new;
    }
    else if(run == 56719){
    double E_scattered = 5.71379 - uncorrected_beam_energy;
    double E_gamma_new = 5.71379 - E_scattered*0.997231;
    ret = E_gamma_new;
    }
    else if(run == 56720){
    double E_scattered = 5.7138 - uncorrected_beam_energy;
    double E_gamma_new = 5.7138 - E_scattered*0.997442;
    ret = E_gamma_new;
    }
    else if(run == 56721){
    double E_scattered = 5.71363 - uncorrected_beam_energy;
    double E_gamma_new = 5.71363 - E_scattered*0.997509;
    ret = E_gamma_new;
    }
    else if(run == 56722){
    double E_scattered = 5.7134 - uncorrected_beam_energy;
    double E_gamma_new = 5.7134 - E_scattered*0.997441;
    ret = E_gamma_new;
    }
    else if(run == 56723){
    double E_scattered = 5.71333 - uncorrected_beam_energy;
    double E_gamma_new = 5.71333 - E_scattered*0.997347;
    ret = E_gamma_new;
    }
    else if(run == 56724){
    double E_scattered = 5.71314 - uncorrected_beam_energy;
    double E_gamma_new = 5.71314 - E_scattered*0.997091;
    ret = E_gamma_new;
    }
    else if(run == 56725){
    double E_scattered = 5.7134 - uncorrected_beam_energy;
    double E_gamma_new = 5.7134 - E_scattered*0.997093;
    ret = E_gamma_new;
    }
    else if(run == 56726){
    double E_scattered = 5.71333 - uncorrected_beam_energy;
    double E_gamma_new = 5.71333 - E_scattered*0.997498;
    ret = E_gamma_new;
    }
    else if(run == 56727){
    double E_scattered = 5.71331 - uncorrected_beam_energy;
    double E_gamma_new = 5.71331 - E_scattered*0.9975;
    ret = E_gamma_new;
    }
    else if(run == 56728){
    double E_scattered = 5.71348 - uncorrected_beam_energy;
    double E_gamma_new = 5.71348 - E_scattered*0.997569;
    ret = E_gamma_new;
    }
    else if(run == 56729){
    double E_scattered = 5.71358 - uncorrected_beam_energy;
    double E_gamma_new = 5.71358 - E_scattered*0.997306;
    ret = E_gamma_new;
    }
    else if(run == 56730){
    double E_scattered = 5.7136 - uncorrected_beam_energy;
    double E_gamma_new = 5.7136 - E_scattered*0.997375;
    ret = E_gamma_new;
    }
    else if(run == 56731){
    double E_scattered = 5.7136 - uncorrected_beam_energy;
    double E_gamma_new = 5.7136 - E_scattered*0.997409;
    ret = E_gamma_new;
    }
    else if(run == 56732){
    double E_scattered = 5.7136 - uncorrected_beam_energy;
    double E_gamma_new = 5.7136 - E_scattered*0.997392;
    ret = E_gamma_new;
    }
    else if(run == 56733){
    double E_scattered = 5.71373 - uncorrected_beam_energy;
    double E_gamma_new = 5.71373 - E_scattered*0.997397;
    ret = E_gamma_new;
    }
    else if(run == 56734){
    double E_scattered = 5.71322 - uncorrected_beam_energy;
    double E_gamma_new = 5.71322 - E_scattered*0.997451;
    ret = E_gamma_new;
    }
    else if(run == 56735){
    double E_scattered = 5.71339 - uncorrected_beam_energy;
    double E_gamma_new = 5.71339 - E_scattered*0.997329;
    ret = E_gamma_new;
    }
    else if(run == 56736){
    double E_scattered = 5.71339 - uncorrected_beam_energy;
    double E_gamma_new = 5.71339 - E_scattered*0.997029;
    ret = E_gamma_new;
    }
    else if(run == 56737){
    double E_scattered = 5.71339 - uncorrected_beam_energy;
    double E_gamma_new = 5.71339 - E_scattered*0.99734;
    ret = E_gamma_new;
    }
    else if(run == 56738){
    double E_scattered = 5.71337 - uncorrected_beam_energy;
    double E_gamma_new = 5.71337 - E_scattered*0.997559;
    ret = E_gamma_new;
    }
    else if(run == 56739){
    double E_scattered = 5.71338 - uncorrected_beam_energy;
    double E_gamma_new = 5.71338 - E_scattered*0.996985;
    ret = E_gamma_new;
    }
    else if(run == 56740){
    double E_scattered = 5.71339 - uncorrected_beam_energy;
    double E_gamma_new = 5.71339 - E_scattered*0.997211;
    ret = E_gamma_new;
    }
    else if(run == 56741){
    double E_scattered = 5.71339 - uncorrected_beam_energy;
    double E_gamma_new = 5.71339 - E_scattered*0.997503;
    ret = E_gamma_new;
    }
    else if(run == 56742){
    double E_scattered = 5.71338 - uncorrected_beam_energy;
    double E_gamma_new = 5.71338 - E_scattered*0.996727;
    ret = E_gamma_new;
    }
    else if(run == 56743){
    double E_scattered = 5.71338 - uncorrected_beam_energy;
    double E_gamma_new = 5.71338 - E_scattered*0.997202;
    ret = E_gamma_new;
    }
    else if(run == 56744){
    double E_scattered = 5.71337 - uncorrected_beam_energy;
    double E_gamma_new = 5.71337 - E_scattered*0.997326;
    ret = E_gamma_new;
    }
    else if(run == 56747){
    double E_scattered = 5.71155 - uncorrected_beam_energy;
    double E_gamma_new = 5.71155 - E_scattered*0.996944;
    ret = E_gamma_new;
    }
    else if(run == 56748){
    double E_scattered = 5.71167 - uncorrected_beam_energy;
    double E_gamma_new = 5.71167 - E_scattered*0.997454;
    ret = E_gamma_new;
    }
    else if(run == 56749){
    double E_scattered = 5.71167 - uncorrected_beam_energy;
    double E_gamma_new = 5.71167 - E_scattered*0.997072;
    ret = E_gamma_new;
    }
    else if(run == 56750){
    double E_scattered = 5.71167 - uncorrected_beam_energy;
    double E_gamma_new = 5.71167 - E_scattered*0.997251;
    ret = E_gamma_new;
    }
    else if(run == 56751){
    double E_scattered = 5.71158 - uncorrected_beam_energy;
    double E_gamma_new = 5.71158 - E_scattered*0.997279;
    ret = E_gamma_new;
    }
    else if(run == 56752){
    double E_scattered = 5.71158 - uncorrected_beam_energy;
    double E_gamma_new = 5.71158 - E_scattered*0.997534;
    ret = E_gamma_new;
    }
    else if(run == 56753){
    double E_scattered = 5.71185 - uncorrected_beam_energy;
    double E_gamma_new = 5.71185 - E_scattered*0.997036;
    ret = E_gamma_new;
    }
    else if(run == 56754){
    double E_scattered = 5.71094 - uncorrected_beam_energy;
    double E_gamma_new = 5.71094 - E_scattered*0.996978;
    ret = E_gamma_new;
    }
    else if(run == 56755){
    double E_scattered = 5.71087 - uncorrected_beam_energy;
    double E_gamma_new = 5.71087 - E_scattered*0.997694;
    ret = E_gamma_new;
    }
    else if(run == 56756){
    double E_scattered = 5.71218 - uncorrected_beam_energy;
    double E_gamma_new = 5.71218 - E_scattered*0.997203;
    ret = E_gamma_new;
    }
    else if(run == 56757){
    double E_scattered = 5.71214 - uncorrected_beam_energy;
    double E_gamma_new = 5.71214 - E_scattered*0.997374;
    ret = E_gamma_new;
    }
    else if(run == 56758){
    double E_scattered = 5.71209 - uncorrected_beam_energy;
    double E_gamma_new = 5.71209 - E_scattered*0.997065;
    ret = E_gamma_new;
    }
    else if(run == 56759){
    double E_scattered = 5.71198 - uncorrected_beam_energy;
    double E_gamma_new = 5.71198 - E_scattered*0.996821;
    ret = E_gamma_new;
    }
    else if(run == 56760){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.997315;
    ret = E_gamma_new;
    }
    else if(run == 56761){
    double E_scattered = 5.71204 - uncorrected_beam_energy;
    double E_gamma_new = 5.71204 - E_scattered*0.996864;
    ret = E_gamma_new;
    }
    else if(run == 56762){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.997421;
    ret = E_gamma_new;
    }
    else if(run == 56763){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.996934;
    ret = E_gamma_new;
    }
    else if(run == 56764){
    double E_scattered = 5.71211 - uncorrected_beam_energy;
    double E_gamma_new = 5.71211 - E_scattered*0.997051;
    ret = E_gamma_new;
    }
    else if(run == 56765){
    double E_scattered = 5.71117 - uncorrected_beam_energy;
    double E_gamma_new = 5.71117 - E_scattered*0.996819;
    ret = E_gamma_new;
    }
    else if(run == 56766){
    double E_scattered = 5.71157 - uncorrected_beam_energy;
    double E_gamma_new = 5.71157 - E_scattered*0.99716;
    ret = E_gamma_new;
    }
    else if(run == 56767){
    double E_scattered = 5.71176 - uncorrected_beam_energy;
    double E_gamma_new = 5.71176 - E_scattered*0.996813;
    ret = E_gamma_new;
    }
    else if(run == 56768){
    double E_scattered = 5.71126 - uncorrected_beam_energy;
    double E_gamma_new = 5.71126 - E_scattered*0.997314;
    ret = E_gamma_new;
    }
    else if(run == 56769){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.997391;
    ret = E_gamma_new;
    }
    else if(run == 56770){
    double E_scattered = 5.71117 - uncorrected_beam_energy;
    double E_gamma_new = 5.71117 - E_scattered*0.997525;
    ret = E_gamma_new;
    }
    else if(run == 56771){
    double E_scattered = 5.71119 - uncorrected_beam_energy;
    double E_gamma_new = 5.71119 - E_scattered*0.99726;
    ret = E_gamma_new;
    }
    else if(run == 56772){
    double E_scattered = 5.71121 - uncorrected_beam_energy;
    double E_gamma_new = 5.71121 - E_scattered*0.997177;
    ret = E_gamma_new;
    }
    else if(run == 56774){
    double E_scattered = 5.71119 - uncorrected_beam_energy;
    double E_gamma_new = 5.71119 - E_scattered*0.997396;
    ret = E_gamma_new;
    }
    else if(run == 56775){
    double E_scattered = 5.71115 - uncorrected_beam_energy;
    double E_gamma_new = 5.71115 - E_scattered*0.997286;
    ret = E_gamma_new;
    }
    else if(run == 56776){
    double E_scattered = 5.71113 - uncorrected_beam_energy;
    double E_gamma_new = 5.71113 - E_scattered*0.997384;
    ret = E_gamma_new;
    }
    else if(run == 56777){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.997512;
    ret = E_gamma_new;
    }
    else if(run == 56778){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.997291;
    ret = E_gamma_new;
    }
    else if(run == 56780){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.997004;
    ret = E_gamma_new;
    }
    else if(run == 56781){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.997245;
    ret = E_gamma_new;
    }
    else if(run == 56782){
    double E_scattered = 5.71127 - uncorrected_beam_energy;
    double E_gamma_new = 5.71127 - E_scattered*0.997391;
    ret = E_gamma_new;
    }
    else if(run == 56783){
    double E_scattered = 5.71136 - uncorrected_beam_energy;
    double E_gamma_new = 5.71136 - E_scattered*0.997273;
    ret = E_gamma_new;
    }
    else if(run == 56784){
    double E_scattered = 5.7112 - uncorrected_beam_energy;
    double E_gamma_new = 5.7112 - E_scattered*0.99711;
    ret = E_gamma_new;
    }
    else if(run == 56787){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.99733;
    ret = E_gamma_new;
    }
    else if(run == 56788){
    double E_scattered = 5.71394 - uncorrected_beam_energy;
    double E_gamma_new = 5.71394 - E_scattered*0.997314;
    ret = E_gamma_new;
    }
    else if(run == 56791){
    double E_scattered = 5.71126 - uncorrected_beam_energy;
    double E_gamma_new = 5.71126 - E_scattered*0.996957;
    ret = E_gamma_new;
    }
    else if(run == 56792){
    double E_scattered = 5.71137 - uncorrected_beam_energy;
    double E_gamma_new = 5.71137 - E_scattered*0.997174;
    ret = E_gamma_new;
    }
    else if(run == 56793){
    double E_scattered = 5.71168 - uncorrected_beam_energy;
    double E_gamma_new = 5.71168 - E_scattered*0.997253;
    ret = E_gamma_new;
    }
    else if(run == 56794){
    double E_scattered = 5.7117 - uncorrected_beam_energy;
    double E_gamma_new = 5.7117 - E_scattered*0.997274;
    ret = E_gamma_new;
    }
    else if(run == 56798){
    double E_scattered = 5.71162 - uncorrected_beam_energy;
    double E_gamma_new = 5.71162 - E_scattered*0.99741;
    ret = E_gamma_new;
    }
    else if(run == 56799){
    double E_scattered = 5.71167 - uncorrected_beam_energy;
    double E_gamma_new = 5.71167 - E_scattered*0.99736;
    ret = E_gamma_new;
    }
    else if(run == 56800){
    double E_scattered = 5.7117 - uncorrected_beam_energy;
    double E_gamma_new = 5.7117 - E_scattered*0.997494;
    ret = E_gamma_new;
    }
    else if(run == 56801){
    double E_scattered = 5.71168 - uncorrected_beam_energy;
    double E_gamma_new = 5.71168 - E_scattered*0.997358;
    ret = E_gamma_new;
    }
    else if(run == 56802){
    double E_scattered = 5.71117 - uncorrected_beam_energy;
    double E_gamma_new = 5.71117 - E_scattered*0.997379;
    ret = E_gamma_new;
    }
    else if(run == 56804){
    double E_scattered = 5.71135 - uncorrected_beam_energy;
    double E_gamma_new = 5.71135 - E_scattered*0.997579;
    ret = E_gamma_new;
    }
    else if(run == 56805){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.996911;
    ret = E_gamma_new;
    }
    else if(run == 56806){
    double E_scattered = 5.71215 - uncorrected_beam_energy;
    double E_gamma_new = 5.71215 - E_scattered*0.997212;
    ret = E_gamma_new;
    }
    else if(run == 56807){
    double E_scattered = 5.71236 - uncorrected_beam_energy;
    double E_gamma_new = 5.71236 - E_scattered*0.997428;
    ret = E_gamma_new;
    }
    else if(run == 56808){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.997037;
    ret = E_gamma_new;
    }
    else if(run == 56809){
    double E_scattered = 5.71089 - uncorrected_beam_energy;
    double E_gamma_new = 5.71089 - E_scattered*0.99706;
    ret = E_gamma_new;
    }
    else if(run == 56810){
    double E_scattered = 5.71348 - uncorrected_beam_energy;
    double E_gamma_new = 5.71348 - E_scattered*0.997406;
    ret = E_gamma_new;
    }
    else if(run == 56811){
    double E_scattered = 5.7123 - uncorrected_beam_energy;
    double E_gamma_new = 5.7123 - E_scattered*0.997219;
    ret = E_gamma_new;
    }
    else if(run == 56812){
    double E_scattered = 5.71238 - uncorrected_beam_energy;
    double E_gamma_new = 5.71238 - E_scattered*0.997251;
    ret = E_gamma_new;
    }
    else if(run == 56813){
    double E_scattered = 5.71309 - uncorrected_beam_energy;
    double E_gamma_new = 5.71309 - E_scattered*0.997349;
    ret = E_gamma_new;
    }
    else if(run == 56814){
    double E_scattered = 5.71106 - uncorrected_beam_energy;
    double E_gamma_new = 5.71106 - E_scattered*0.99716;
    ret = E_gamma_new;
    }
    else if(run == 56815){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.997285;
    ret = E_gamma_new;
    }
    else if(run == 56821){
    double E_scattered = 5.71303 - uncorrected_beam_energy;
    double E_gamma_new = 5.71303 - E_scattered*0.997308;
    ret = E_gamma_new;
    }
    else if(run == 56822){
    double E_scattered = 5.7122 - uncorrected_beam_energy;
    double E_gamma_new = 5.7122 - E_scattered*0.997261;
    ret = E_gamma_new;
    }
    else if(run == 56823){
    double E_scattered = 5.71233 - uncorrected_beam_energy;
    double E_gamma_new = 5.71233 - E_scattered*0.996962;
    ret = E_gamma_new;
    }
    else if(run == 56824){
    double E_scattered = 5.7123 - uncorrected_beam_energy;
    double E_gamma_new = 5.7123 - E_scattered*0.996903;
    ret = E_gamma_new;
    }
    else if(run == 56825){
    double E_scattered = 5.71485 - uncorrected_beam_energy;
    double E_gamma_new = 5.71485 - E_scattered*0.997423;
    ret = E_gamma_new;
    }
    else if(run == 56826){
    double E_scattered = 5.7123 - uncorrected_beam_energy;
    double E_gamma_new = 5.7123 - E_scattered*0.997408;
    ret = E_gamma_new;
    }
    else if(run == 56827){
    double E_scattered = 5.71266 - uncorrected_beam_energy;
    double E_gamma_new = 5.71266 - E_scattered*0.997362;
    ret = E_gamma_new;
    }
    else if(run == 56831){
    double E_scattered = 5.7123 - uncorrected_beam_energy;
    double E_gamma_new = 5.7123 - E_scattered*0.997156;
    ret = E_gamma_new;
    }
    else if(run == 56832){
    double E_scattered = 5.71214 - uncorrected_beam_energy;
    double E_gamma_new = 5.71214 - E_scattered*0.997151;
    ret = E_gamma_new;
    }
    else if(run == 56833){
    double E_scattered = 5.71273 - uncorrected_beam_energy;
    double E_gamma_new = 5.71273 - E_scattered*0.99726;
    ret = E_gamma_new;
    }
    else if(run == 56834){
    double E_scattered = 5.71641 - uncorrected_beam_energy;
    double E_gamma_new = 5.71641 - E_scattered*0.997972;
    ret = E_gamma_new;
    }
    else if(run == 56835){
    double E_scattered = 5.71305 - uncorrected_beam_energy;
    double E_gamma_new = 5.71305 - E_scattered*0.997113;
    ret = E_gamma_new;
    }
    else if(run == 56838){
    double E_scattered = 5.71179 - uncorrected_beam_energy;
    double E_gamma_new = 5.71179 - E_scattered*0.997502;
    ret = E_gamma_new;
    }
    else if(run == 56839){
    double E_scattered = 5.71199 - uncorrected_beam_energy;
    double E_gamma_new = 5.71199 - E_scattered*0.997145;
    ret = E_gamma_new;
    }
    else if(run == 56841){
    double E_scattered = 5.7118 - uncorrected_beam_energy;
    double E_gamma_new = 5.7118 - E_scattered*0.997289;
    ret = E_gamma_new;
    }
    else if(run == 56842){
    double E_scattered = 5.71175 - uncorrected_beam_energy;
    double E_gamma_new = 5.71175 - E_scattered*0.997542;
    ret = E_gamma_new;
    }
    else if(run == 56843){
    double E_scattered = 5.71172 - uncorrected_beam_energy;
    double E_gamma_new = 5.71172 - E_scattered*0.996853;
    ret = E_gamma_new;
    }
    else if(run == 56844){
    double E_scattered = 5.71731 - uncorrected_beam_energy;
    double E_gamma_new = 5.71731 - E_scattered*0.99767;
    ret = E_gamma_new;
    }
    else if(run == 56845){
    double E_scattered = 5.71168 - uncorrected_beam_energy;
    double E_gamma_new = 5.71168 - E_scattered*0.99715;
    ret = E_gamma_new;
    }
    else if(run == 56849){
    double E_scattered = 5.71139 - uncorrected_beam_energy;
    double E_gamma_new = 5.71139 - E_scattered*0.997221;
    ret = E_gamma_new;
    }
    else if(run == 56853){
    double E_scattered = 5.71251 - uncorrected_beam_energy;
    double E_gamma_new = 5.71251 - E_scattered*0.99686;
    ret = E_gamma_new;
    }
    else if(run == 56854){
    double E_scattered = 5.71278 - uncorrected_beam_energy;
    double E_gamma_new = 5.71278 - E_scattered*0.997171;
    ret = E_gamma_new;
    }
    else if(run == 56855){
    double E_scattered = 5.71194 - uncorrected_beam_energy;
    double E_gamma_new = 5.71194 - E_scattered*0.997058;
    ret = E_gamma_new;
    }
    else if(run == 56856){
    double E_scattered = 5.71281 - uncorrected_beam_energy;
    double E_gamma_new = 5.71281 - E_scattered*0.997134;
    ret = E_gamma_new;
    }
    else if(run == 56857){
    double E_scattered = 5.71202 - uncorrected_beam_energy;
    double E_gamma_new = 5.71202 - E_scattered*0.996693;
    ret = E_gamma_new;
    }
    else if(run == 56858){
    double E_scattered = 5.71206 - uncorrected_beam_energy;
    double E_gamma_new = 5.71206 - E_scattered*0.997064;
    ret = E_gamma_new;
    }
    else if(run == 56859){
    double E_scattered = 5.7121 - uncorrected_beam_energy;
    double E_gamma_new = 5.7121 - E_scattered*0.997071;
    ret = E_gamma_new;
    }
    else if(run == 56860){
    double E_scattered = 5.71183 - uncorrected_beam_energy;
    double E_gamma_new = 5.71183 - E_scattered*0.996815;
    ret = E_gamma_new;
    }
    else if(run == 56861){
    double E_scattered = 5.71116 - uncorrected_beam_energy;
    double E_gamma_new = 5.71116 - E_scattered*0.997198;
    ret = E_gamma_new;
    }
    else if(run == 56862){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.997174;
    ret = E_gamma_new;
    }
    else if(run == 56864){
    double E_scattered = 5.71083 - uncorrected_beam_energy;
    double E_gamma_new = 5.71083 - E_scattered*0.997437;
    ret = E_gamma_new;
    }
    else if(run == 56865){
    double E_scattered = 5.71025 - uncorrected_beam_energy;
    double E_gamma_new = 5.71025 - E_scattered*0.99712;
    ret = E_gamma_new;
    }
    else if(run == 56866){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997294;
    ret = E_gamma_new;
    }
    else if(run == 56869){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.99745;
    ret = E_gamma_new;
    }
    else if(run == 56870){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997278;
    ret = E_gamma_new;
    }
    else if(run == 56874){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997061;
    ret = E_gamma_new;
    }
    else if(run == 56875){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997491;
    ret = E_gamma_new;
    }
    else if(run == 56877){
    double E_scattered = 5.71065 - uncorrected_beam_energy;
    double E_gamma_new = 5.71065 - E_scattered*0.997502;
    ret = E_gamma_new;
    }
    else if(run == 56879){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997488;
    ret = E_gamma_new;
    }
    else if(run == 56897){
    double E_scattered = 5.71133 - uncorrected_beam_energy;
    double E_gamma_new = 5.71133 - E_scattered*0.997979;
    ret = E_gamma_new;
    }
    else if(run == 56898){
    double E_scattered = 5.7113 - uncorrected_beam_energy;
    double E_gamma_new = 5.7113 - E_scattered*0.997759;
    ret = E_gamma_new;
    }
    else if(run == 56899){
    double E_scattered = 5.71425 - uncorrected_beam_energy;
    double E_gamma_new = 5.71425 - E_scattered*0.997855;
    ret = E_gamma_new;
    }
    else if(run == 56900){
    double E_scattered = 5.7147 - uncorrected_beam_energy;
    double E_gamma_new = 5.7147 - E_scattered*0.998091;
    ret = E_gamma_new;
    }
    else if(run == 56901){
    double E_scattered = 5.71128 - uncorrected_beam_energy;
    double E_gamma_new = 5.71128 - E_scattered*0.997816;
    ret = E_gamma_new;
    }
    else if(run == 56902){
    double E_scattered = 5.71153 - uncorrected_beam_energy;
    double E_gamma_new = 5.71153 - E_scattered*0.997769;
    ret = E_gamma_new;
    }
    else if(run == 56903){
    double E_scattered = 5.71133 - uncorrected_beam_energy;
    double E_gamma_new = 5.71133 - E_scattered*0.997856;
    ret = E_gamma_new;
    }
    else if(run == 56904){
    double E_scattered = 5.71133 - uncorrected_beam_energy;
    double E_gamma_new = 5.71133 - E_scattered*0.997648;
    ret = E_gamma_new;
    }
    else if(run == 56905){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.99765;
    ret = E_gamma_new;
    }
    else if(run == 56906){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.997803;
    ret = E_gamma_new;
    }
    else if(run == 56907){
    double E_scattered = 5.70996 - uncorrected_beam_energy;
    double E_gamma_new = 5.70996 - E_scattered*0.998029;
    ret = E_gamma_new;
    }
    else if(run == 56908){
    double E_scattered = 5.71083 - uncorrected_beam_energy;
    double E_gamma_new = 5.71083 - E_scattered*0.998202;
    ret = E_gamma_new;
    }
    else if(run == 56910){
    double E_scattered = 5.71083 - uncorrected_beam_energy;
    double E_gamma_new = 5.71083 - E_scattered*0.997591;
    ret = E_gamma_new;
    }
    else if(run == 56911){
    double E_scattered = 5.71083 - uncorrected_beam_energy;
    double E_gamma_new = 5.71083 - E_scattered*0.997588;
    ret = E_gamma_new;
    }
    else if(run == 56912){
    double E_scattered = 5.71083 - uncorrected_beam_energy;
    double E_gamma_new = 5.71083 - E_scattered*0.998189;
    ret = E_gamma_new;
    }
    else if(run == 56913){
    double E_scattered = 5.7113 - uncorrected_beam_energy;
    double E_gamma_new = 5.7113 - E_scattered*0.998032;
    ret = E_gamma_new;
    }
    else if(run == 56914){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998049;
    ret = E_gamma_new;
    }
    else if(run == 56915){
    double E_scattered = 5.71173 - uncorrected_beam_energy;
    double E_gamma_new = 5.71173 - E_scattered*0.997785;
    ret = E_gamma_new;
    }
    else if(run == 56916){
    double E_scattered = 5.71109 - uncorrected_beam_energy;
    double E_gamma_new = 5.71109 - E_scattered*0.997806;
    ret = E_gamma_new;
    }
    else if(run == 56917){
    double E_scattered = 5.7106 - uncorrected_beam_energy;
    double E_gamma_new = 5.7106 - E_scattered*0.997779;
    ret = E_gamma_new;
    }
    else if(run == 56918){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.997936;
    ret = E_gamma_new;
    }
    else if(run == 56919){
    double E_scattered = 5.71115 - uncorrected_beam_energy;
    double E_gamma_new = 5.71115 - E_scattered*0.997871;
    ret = E_gamma_new;
    }
    else if(run == 56921){
    double E_scattered = 5.70913 - uncorrected_beam_energy;
    double E_gamma_new = 5.70913 - E_scattered*0.998054;
    ret = E_gamma_new;
    }
    else if(run == 56922){
    double E_scattered = 5.70993 - uncorrected_beam_energy;
    double E_gamma_new = 5.70993 - E_scattered*0.998358;
    ret = E_gamma_new;
    }
    else if(run == 56923){
    double E_scattered = 5.71055 - uncorrected_beam_energy;
    double E_gamma_new = 5.71055 - E_scattered*0.997904;
    ret = E_gamma_new;
    }
    else if(run == 56924){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.997499;
    ret = E_gamma_new;
    }
    else if(run == 56925){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.997496;
    ret = E_gamma_new;
    }
    else if(run == 56926){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.99713;
    ret = E_gamma_new;
    }
    else if(run == 56927){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.997257;
    ret = E_gamma_new;
    }
    else if(run == 56928){
    double E_scattered = 5.71055 - uncorrected_beam_energy;
    double E_gamma_new = 5.71055 - E_scattered*0.997384;
    ret = E_gamma_new;
    }
    else if(run == 56929){
    double E_scattered = 5.7106 - uncorrected_beam_energy;
    double E_gamma_new = 5.7106 - E_scattered*0.997318;
    ret = E_gamma_new;
    }
    else if(run == 56930){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.997361;
    ret = E_gamma_new;
    }
    else if(run == 56932){
    double E_scattered = 5.71067 - uncorrected_beam_energy;
    double E_gamma_new = 5.71067 - E_scattered*1.00114;
    ret = E_gamma_new;
    }
    else if(run == 56933){
    double E_scattered = 5.71057 - uncorrected_beam_energy;
    double E_gamma_new = 5.71057 - E_scattered*0.99838;
    ret = E_gamma_new;
    }
    else if(run == 56934){
    double E_scattered = 5.71055 - uncorrected_beam_energy;
    double E_gamma_new = 5.71055 - E_scattered*0.998314;
    ret = E_gamma_new;
    }
    else if(run == 56935){
    double E_scattered = 5.71195 - uncorrected_beam_energy;
    double E_gamma_new = 5.71195 - E_scattered*0.998502;
    ret = E_gamma_new;
    }
    else if(run == 56936){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.997102;
    ret = E_gamma_new;
    }
    else if(run == 56937){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.998683;
    ret = E_gamma_new;
    }
    else if(run == 56938){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.997271;
    ret = E_gamma_new;
    }
    else if(run == 56939){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.9971;
    ret = E_gamma_new;
    }
    else if(run == 56940){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998248;
    ret = E_gamma_new;
    }
    else if(run == 56948){
    double E_scattered = 5.71093 - uncorrected_beam_energy;
    double E_gamma_new = 5.71093 - E_scattered*0.998451;
    ret = E_gamma_new;
    }
    else if(run == 56949){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.998197;
    ret = E_gamma_new;
    }
    else if(run == 56950){
    double E_scattered = 5.71059 - uncorrected_beam_energy;
    double E_gamma_new = 5.71059 - E_scattered*0.9984;
    ret = E_gamma_new;
    }
    else if(run == 56951){
    double E_scattered = 5.71057 - uncorrected_beam_energy;
    double E_gamma_new = 5.71057 - E_scattered*0.998242;
    ret = E_gamma_new;
    }
    else if(run == 56952){
    double E_scattered = 5.71062 - uncorrected_beam_energy;
    double E_gamma_new = 5.71062 - E_scattered*0.998333;
    ret = E_gamma_new;
    }
    else if(run == 56953){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998353;
    ret = E_gamma_new;
    }
    else if(run == 56954){
    double E_scattered = 5.71127 - uncorrected_beam_energy;
    double E_gamma_new = 5.71127 - E_scattered*0.998387;
    ret = E_gamma_new;
    }
    else if(run == 56955){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998427;
    ret = E_gamma_new;
    }
    else if(run == 56956){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.998327;
    ret = E_gamma_new;
    }
    else if(run == 56958){
    double E_scattered = 5.71028 - uncorrected_beam_energy;
    double E_gamma_new = 5.71028 - E_scattered*0.998552;
    ret = E_gamma_new;
    }
    else if(run == 56960){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.998301;
    ret = E_gamma_new;
    }
    else if(run == 56961){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998371;
    ret = E_gamma_new;
    }
    else if(run == 56962){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998303;
    ret = E_gamma_new;
    }
    else if(run == 56963){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998128;
    ret = E_gamma_new;
    }
    else if(run == 56964){
    double E_scattered = 5.71048 - uncorrected_beam_energy;
    double E_gamma_new = 5.71048 - E_scattered*0.998163;
    ret = E_gamma_new;
    }
    else if(run == 56965){
    double E_scattered = 5.712 - uncorrected_beam_energy;
    double E_gamma_new = 5.712 - E_scattered*0.998259;
    ret = E_gamma_new;
    }
    else if(run == 56966){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.998154;
    ret = E_gamma_new;
    }
    else if(run == 56967){
    double E_scattered = 5.71122 - uncorrected_beam_energy;
    double E_gamma_new = 5.71122 - E_scattered*0.998269;
    ret = E_gamma_new;
    }
    else if(run == 56968){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998219;
    ret = E_gamma_new;
    }
    else if(run == 56969){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998261;
    ret = E_gamma_new;
    }
    else if(run == 56970){
    double E_scattered = 5.71138 - uncorrected_beam_energy;
    double E_gamma_new = 5.71138 - E_scattered*0.998345;
    ret = E_gamma_new;
    }
    else if(run == 56971){
    double E_scattered = 5.7106 - uncorrected_beam_energy;
    double E_gamma_new = 5.7106 - E_scattered*0.99838;
    ret = E_gamma_new;
    }
    else if(run == 56972){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.99836;
    ret = E_gamma_new;
    }
    else if(run == 56973){
    double E_scattered = 5.71065 - uncorrected_beam_energy;
    double E_gamma_new = 5.71065 - E_scattered*0.998189;
    ret = E_gamma_new;
    }
    else if(run == 56974){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998261;
    ret = E_gamma_new;
    }
    else if(run == 56975){
    double E_scattered = 5.71055 - uncorrected_beam_energy;
    double E_gamma_new = 5.71055 - E_scattered*0.998299;
    ret = E_gamma_new;
    }
    else if(run == 56977){
    double E_scattered = 5.71016 - uncorrected_beam_energy;
    double E_gamma_new = 5.71016 - E_scattered*0.998273;
    ret = E_gamma_new;
    }
    else if(run == 56978){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998374;
    ret = E_gamma_new;
    }
    else if(run == 56979){
    double E_scattered = 5.71157 - uncorrected_beam_energy;
    double E_gamma_new = 5.71157 - E_scattered*0.998468;
    ret = E_gamma_new;
    }
    else if(run == 56980){
    double E_scattered = 5.71156 - uncorrected_beam_energy;
    double E_gamma_new = 5.71156 - E_scattered*0.998308;
    ret = E_gamma_new;
    }
    else if(run == 56981){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.99821;
    ret = E_gamma_new;
    }
    else if(run == 56982){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998724;
    ret = E_gamma_new;
    }
    else if(run == 56983){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998164;
    ret = E_gamma_new;
    }
    else if(run == 56985){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998401;
    ret = E_gamma_new;
    }
    else if(run == 56986){
    double E_scattered = 5.71113 - uncorrected_beam_energy;
    double E_gamma_new = 5.71113 - E_scattered*0.998388;
    ret = E_gamma_new;
    }
    else if(run == 56989){
    double E_scattered = 5.71275 - uncorrected_beam_energy;
    double E_gamma_new = 5.71275 - E_scattered*0.998014;
    ret = E_gamma_new;
    }
    else if(run == 56992){
    double E_scattered = 5.71275 - uncorrected_beam_energy;
    double E_gamma_new = 5.71275 - E_scattered*0.998457;
    ret = E_gamma_new;
    }
    else if(run == 56993){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.998376;
    ret = E_gamma_new;
    }
    else if(run == 56994){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998429;
    ret = E_gamma_new;
    }
    else if(run == 56996){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998266;
    ret = E_gamma_new;
    }
    else if(run == 56997){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998404;
    ret = E_gamma_new;
    }
    else if(run == 56998){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998156;
    ret = E_gamma_new;
    }
    else if(run == 56999){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998263;
    ret = E_gamma_new;
    }
    else if(run == 57000){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998441;
    ret = E_gamma_new;
    }
    else if(run == 57001){
    double E_scattered = 5.71049 - uncorrected_beam_energy;
    double E_gamma_new = 5.71049 - E_scattered*0.998239;
    ret = E_gamma_new;
    }
    else if(run == 57002){
    double E_scattered = 5.71015 - uncorrected_beam_energy;
    double E_gamma_new = 5.71015 - E_scattered*0.99832;
    ret = E_gamma_new;
    }
    else if(run == 57003){
    double E_scattered = 5.7119 - uncorrected_beam_energy;
    double E_gamma_new = 5.7119 - E_scattered*0.998351;
    ret = E_gamma_new;
    }
    else if(run == 57004){
    double E_scattered = 5.71052 - uncorrected_beam_energy;
    double E_gamma_new = 5.71052 - E_scattered*0.998205;
    ret = E_gamma_new;
    }
    else if(run == 57005){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.99815;
    ret = E_gamma_new;
    }
    else if(run == 57006){
    double E_scattered = 5.71056 - uncorrected_beam_energy;
    double E_gamma_new = 5.71056 - E_scattered*0.998185;
    ret = E_gamma_new;
    }
    else if(run == 57008){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.998201;
    ret = E_gamma_new;
    }
    else if(run == 57009){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998183;
    ret = E_gamma_new;
    }
    else if(run == 57010){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998114;
    ret = E_gamma_new;
    }
    else if(run == 57011){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998149;
    ret = E_gamma_new;
    }
    else if(run == 57012){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998146;
    ret = E_gamma_new;
    }
    else if(run == 57013){
    double E_scattered = 5.71059 - uncorrected_beam_energy;
    double E_gamma_new = 5.71059 - E_scattered*0.99831;
    ret = E_gamma_new;
    }
    else if(run == 57014){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998295;
    ret = E_gamma_new;
    }
    else if(run == 57015){
    double E_scattered = 5.71059 - uncorrected_beam_energy;
    double E_gamma_new = 5.71059 - E_scattered*0.99822;
    ret = E_gamma_new;
    }
    else if(run == 57016){
    double E_scattered = 5.71057 - uncorrected_beam_energy;
    double E_gamma_new = 5.71057 - E_scattered*0.998146;
    ret = E_gamma_new;
    }
    else if(run == 57017){
    double E_scattered = 5.71059 - uncorrected_beam_energy;
    double E_gamma_new = 5.71059 - E_scattered*0.998122;
    ret = E_gamma_new;
    }
    else if(run == 57021){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.99816;
    ret = E_gamma_new;
    }
    else if(run == 57022){
    double E_scattered = 5.71229 - uncorrected_beam_energy;
    double E_gamma_new = 5.71229 - E_scattered*0.998161;
    ret = E_gamma_new;
    }
    else if(run == 57023){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.99833;
    ret = E_gamma_new;
    }
    else if(run == 57025){
    double E_scattered = 5.71245 - uncorrected_beam_energy;
    double E_gamma_new = 5.71245 - E_scattered*0.998421;
    ret = E_gamma_new;
    }
    else if(run == 57026){
    double E_scattered = 5.71228 - uncorrected_beam_energy;
    double E_gamma_new = 5.71228 - E_scattered*0.99825;
    ret = E_gamma_new;
    }
    else if(run == 57027){
    double E_scattered = 5.71277 - uncorrected_beam_energy;
    double E_gamma_new = 5.71277 - E_scattered*0.998165;
    ret = E_gamma_new;
    }
    else if(run == 57028){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.99797;
    ret = E_gamma_new;
    }
    else if(run == 57030){
    double E_scattered = 5.71069 - uncorrected_beam_energy;
    double E_gamma_new = 5.71069 - E_scattered*0.99808;
    ret = E_gamma_new;
    }
    else if(run == 57031){
    double E_scattered = 5.70865 - uncorrected_beam_energy;
    double E_gamma_new = 5.70865 - E_scattered*0.997975;
    ret = E_gamma_new;
    }
    else if(run == 57032){
    double E_scattered = 5.71063 - uncorrected_beam_energy;
    double E_gamma_new = 5.71063 - E_scattered*0.997968;
    ret = E_gamma_new;
    }
    else if(run == 57036){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.998146;
    ret = E_gamma_new;
    }
    else if(run == 57037){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.997956;
    ret = E_gamma_new;
    }
    else if(run == 57038){
    double E_scattered = 5.71057 - uncorrected_beam_energy;
    double E_gamma_new = 5.71057 - E_scattered*0.997978;
    ret = E_gamma_new;
    }
    else if(run == 57039){
    double E_scattered = 5.71056 - uncorrected_beam_energy;
    double E_gamma_new = 5.71056 - E_scattered*0.998127;
    ret = E_gamma_new;
    }
    else if(run == 57061){
    double E_scattered = 5.71167 - uncorrected_beam_energy;
    double E_gamma_new = 5.71167 - E_scattered*0.998665;
    ret = E_gamma_new;
    }
    else if(run == 57062){
    double E_scattered = 5.71057 - uncorrected_beam_energy;
    double E_gamma_new = 5.71057 - E_scattered*0.998459;
    ret = E_gamma_new;
    }
    else if(run == 57063){
    double E_scattered = 5.71265 - uncorrected_beam_energy;
    double E_gamma_new = 5.71265 - E_scattered*0.998535;
    ret = E_gamma_new;
    }
    else if(run == 57064){
    double E_scattered = 5.71054 - uncorrected_beam_energy;
    double E_gamma_new = 5.71054 - E_scattered*0.998605;
    ret = E_gamma_new;
    }
    else if(run == 57065){
    double E_scattered = 5.71053 - uncorrected_beam_energy;
    double E_gamma_new = 5.71053 - E_scattered*0.99842;
    ret = E_gamma_new;
    }
    else if(run == 57066){
    double E_scattered = 5.71059 - uncorrected_beam_energy;
    double E_gamma_new = 5.71059 - E_scattered*0.998149;
    ret = E_gamma_new;
    }
    else if(run == 57067){
    double E_scattered = 5.71148 - uncorrected_beam_energy;
    double E_gamma_new = 5.71148 - E_scattered*0.998306;
    ret = E_gamma_new;
    }
    else if(run == 57068){
    double E_scattered = 5.71035 - uncorrected_beam_energy;
    double E_gamma_new = 5.71035 - E_scattered*0.998218;
    ret = E_gamma_new;
    }
    else if(run == 57069){
    double E_scattered = 5.71127 - uncorrected_beam_energy;
    double E_gamma_new = 5.71127 - E_scattered*0.99843;
    ret = E_gamma_new;
    }
    else if(run == 57071){
    double E_scattered = 5.71058 - uncorrected_beam_energy;
    double E_gamma_new = 5.71058 - E_scattered*0.998194;
    ret = E_gamma_new;
    }
    else if(run == 57072){
    double E_scattered = 5.7105 - uncorrected_beam_energy;
    double E_gamma_new = 5.7105 - E_scattered*0.998827;
    ret = E_gamma_new;
    }
    else if(run == 57073){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.998149;
    ret = E_gamma_new;
    }
    else if(run == 57075){
    double E_scattered = 5.7104 - uncorrected_beam_energy;
    double E_gamma_new = 5.7104 - E_scattered*0.998494;
    ret = E_gamma_new;
    }
    else if(run == 57076){
    double E_scattered = 5.71055 - uncorrected_beam_energy;
    double E_gamma_new = 5.71055 - E_scattered*0.998955;
    ret = E_gamma_new;
    }
    else if(run == 57077){
    double E_scattered = 5.71071 - uncorrected_beam_energy;
    double E_gamma_new = 5.71071 - E_scattered*0.999218;
    ret = E_gamma_new;
    }
    else if(run == 57078){
    double E_scattered = 5.71061 - uncorrected_beam_energy;
    double E_gamma_new = 5.71061 - E_scattered*0.998815;
    ret = E_gamma_new;
    }
    else if(run == 57079){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998983;
    ret = E_gamma_new;
    }
    else if(run == 57080){
    double E_scattered = 5.71094 - uncorrected_beam_energy;
    double E_gamma_new = 5.71094 - E_scattered*0.998977;
    ret = E_gamma_new;
    }
    else if(run == 57094){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998032;
    ret = E_gamma_new;
    }
    else if(run == 57095){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998245;
    ret = E_gamma_new;
    }
    else if(run == 57096){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998244;
    ret = E_gamma_new;
    }
    else if(run == 57097){
    double E_scattered = 5.7122 - uncorrected_beam_energy;
    double E_gamma_new = 5.7122 - E_scattered*0.998238;
    ret = E_gamma_new;
    }
    else if(run == 57100){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998157;
    ret = E_gamma_new;
    }
    else if(run == 57101){
    double E_scattered = 5.71061 - uncorrected_beam_energy;
    double E_gamma_new = 5.71061 - E_scattered*0.998271;
    ret = E_gamma_new;
    }
    else if(run == 57102){
    double E_scattered = 5.71094 - uncorrected_beam_energy;
    double E_gamma_new = 5.71094 - E_scattered*0.997936;
    ret = E_gamma_new;
    }
    else if(run == 57103){
    double E_scattered = 5.71092 - uncorrected_beam_energy;
    double E_gamma_new = 5.71092 - E_scattered*0.998198;
    ret = E_gamma_new;
    }
    else if(run == 57106){
    double E_scattered = 5.7111 - uncorrected_beam_energy;
    double E_gamma_new = 5.7111 - E_scattered*0.997767;
    ret = E_gamma_new;
    }
    else if(run == 57107){
    double E_scattered = 5.71411 - uncorrected_beam_energy;
    double E_gamma_new = 5.71411 - E_scattered*0.997939;
    ret = E_gamma_new;
    }
    else if(run == 57108){
    double E_scattered = 5.7111 - uncorrected_beam_energy;
    double E_gamma_new = 5.7111 - E_scattered*0.997843;
    ret = E_gamma_new;
    }
    else if(run == 57114){
    double E_scattered = 5.71196 - uncorrected_beam_energy;
    double E_gamma_new = 5.71196 - E_scattered*0.998233;
    ret = E_gamma_new;
    }
    else if(run == 57115){
    double E_scattered = 5.71099 - uncorrected_beam_energy;
    double E_gamma_new = 5.71099 - E_scattered*0.998205;
    ret = E_gamma_new;
    }
    else if(run == 57116){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998178;
    ret = E_gamma_new;
    }
    else if(run == 57117){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998111;
    ret = E_gamma_new;
    }
    else if(run == 57118){
    double E_scattered = 5.7132 - uncorrected_beam_energy;
    double E_gamma_new = 5.7132 - E_scattered*0.998416;
    ret = E_gamma_new;
    }
    else if(run == 57119){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.99815;
    ret = E_gamma_new;
    }
    else if(run == 57120){
    double E_scattered = 5.71106 - uncorrected_beam_energy;
    double E_gamma_new = 5.71106 - E_scattered*0.998182;
    ret = E_gamma_new;
    }
    else if(run == 57121){
    double E_scattered = 5.71116 - uncorrected_beam_energy;
    double E_gamma_new = 5.71116 - E_scattered*0.998114;
    ret = E_gamma_new;
    }
    else if(run == 57122){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.998121;
    ret = E_gamma_new;
    }
    else if(run == 57123){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998499;
    ret = E_gamma_new;
    }
    else if(run == 57124){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.99835;
    ret = E_gamma_new;
    }
    else if(run == 57125){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998143;
    ret = E_gamma_new;
    }
    else if(run == 57126){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.99827;
    ret = E_gamma_new;
    }
    else if(run == 57127){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998435;
    ret = E_gamma_new;
    }
    else if(run == 57128){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998136;
    ret = E_gamma_new;
    }
    else if(run == 57129){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.998298;
    ret = E_gamma_new;
    }
    else if(run == 57130){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.99823;
    ret = E_gamma_new;
    }
    else if(run == 57131){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998124;
    ret = E_gamma_new;
    }
    else if(run == 57132){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998295;
    ret = E_gamma_new;
    }
    else if(run == 57133){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998244;
    ret = E_gamma_new;
    }
    else if(run == 57134){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998184;
    ret = E_gamma_new;
    }
    else if(run == 57135){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998233;
    ret = E_gamma_new;
    }
    else if(run == 57136){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.99837;
    ret = E_gamma_new;
    }
    else if(run == 57137){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.998185;
    ret = E_gamma_new;
    }
    else if(run == 57138){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998379;
    ret = E_gamma_new;
    }
    else if(run == 57139){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998324;
    ret = E_gamma_new;
    }
    else if(run == 57140){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998147;
    ret = E_gamma_new;
    }
    else if(run == 57141){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998173;
    ret = E_gamma_new;
    }
    else if(run == 57142){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.99836;
    ret = E_gamma_new;
    }
    else if(run == 57143){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998883;
    ret = E_gamma_new;
    }
    else if(run == 57144){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.99829;
    ret = E_gamma_new;
    }
    else if(run == 57145){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998228;
    ret = E_gamma_new;
    }
    else if(run == 57146){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998075;
    ret = E_gamma_new;
    }
    else if(run == 57147){
    double E_scattered = 5.71102 - uncorrected_beam_energy;
    double E_gamma_new = 5.71102 - E_scattered*0.998593;
    ret = E_gamma_new;
    }
    else if(run == 57148){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998504;
    ret = E_gamma_new;
    }
    else if(run == 57149){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998156;
    ret = E_gamma_new;
    }
    else if(run == 57150){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998033;
    ret = E_gamma_new;
    }
    else if(run == 57151){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.998127;
    ret = E_gamma_new;
    }
    else if(run == 57152){
    double E_scattered = 5.71101 - uncorrected_beam_energy;
    double E_gamma_new = 5.71101 - E_scattered*0.997795;
    ret = E_gamma_new;
    }
    else if(run == 57155){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.99812;
    ret = E_gamma_new;
    }
    else if(run == 57156){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998228;
    ret = E_gamma_new;
    }
    else if(run == 57159){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998297;
    ret = E_gamma_new;
    }
    else if(run == 57160){
    double E_scattered = 5.7118 - uncorrected_beam_energy;
    double E_gamma_new = 5.7118 - E_scattered*0.998474;
    ret = E_gamma_new;
    }
    else if(run == 57161){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998184;
    ret = E_gamma_new;
    }
    else if(run == 57162){
    double E_scattered = 5.71093 - uncorrected_beam_energy;
    double E_gamma_new = 5.71093 - E_scattered*0.998019;
    ret = E_gamma_new;
    }
    else if(run == 57163){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998216;
    ret = E_gamma_new;
    }
    else if(run == 57164){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998302;
    ret = E_gamma_new;
    }
    else if(run == 57165){
    double E_scattered = 5.71403 - uncorrected_beam_energy;
    double E_gamma_new = 5.71403 - E_scattered*0.99828;
    ret = E_gamma_new;
    }
    else if(run == 57166){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998203;
    ret = E_gamma_new;
    }
    else if(run == 57167){
    double E_scattered = 5.7108 - uncorrected_beam_energy;
    double E_gamma_new = 5.7108 - E_scattered*0.998273;
    ret = E_gamma_new;
    }
    else if(run == 57168){
    double E_scattered = 5.71093 - uncorrected_beam_energy;
    double E_gamma_new = 5.71093 - E_scattered*0.998052;
    ret = E_gamma_new;
    }
    else if(run == 57170){
    double E_scattered = 5.71113 - uncorrected_beam_energy;
    double E_gamma_new = 5.71113 - E_scattered*0.998085;
    ret = E_gamma_new;
    }
    else if(run == 57171){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998137;
    ret = E_gamma_new;
    }
    else if(run == 57172){
    double E_scattered = 5.71307 - uncorrected_beam_energy;
    double E_gamma_new = 5.71307 - E_scattered*0.998227;
    ret = E_gamma_new;
    }
    else if(run == 57173){
    double E_scattered = 5.71307 - uncorrected_beam_energy;
    double E_gamma_new = 5.71307 - E_scattered*0.998232;
    ret = E_gamma_new;
    }
    else if(run == 57174){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998199;
    ret = E_gamma_new;
    }
    else if(run == 57175){
    double E_scattered = 5.7119 - uncorrected_beam_energy;
    double E_gamma_new = 5.7119 - E_scattered*0.998303;
    ret = E_gamma_new;
    }
    else if(run == 57176){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998145;
    ret = E_gamma_new;
    }
    else if(run == 57177){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998085;
    ret = E_gamma_new;
    }
    else if(run == 57178){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998253;
    ret = E_gamma_new;
    }
    else if(run == 57179){
    double E_scattered = 5.71303 - uncorrected_beam_energy;
    double E_gamma_new = 5.71303 - E_scattered*0.998406;
    ret = E_gamma_new;
    }
    else if(run == 57180){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998033;
    ret = E_gamma_new;
    }
    else if(run == 57181){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998187;
    ret = E_gamma_new;
    }
    else if(run == 57182){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.9982;
    ret = E_gamma_new;
    }
    else if(run == 57183){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.997856;
    ret = E_gamma_new;
    }
    else if(run == 57184){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998096;
    ret = E_gamma_new;
    }
    else if(run == 57185){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998069;
    ret = E_gamma_new;
    }
    else if(run == 57189){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998049;
    ret = E_gamma_new;
    }
    else if(run == 57190){
    double E_scattered = 5.71104 - uncorrected_beam_energy;
    double E_gamma_new = 5.71104 - E_scattered*0.998093;
    ret = E_gamma_new;
    }
    else if(run == 57191){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.997871;
    ret = E_gamma_new;
    }
    else if(run == 57192){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998242;
    ret = E_gamma_new;
    }
    else if(run == 57193){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998042;
    ret = E_gamma_new;
    }
    else if(run == 57194){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998208;
    ret = E_gamma_new;
    }
    else if(run == 57195){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998145;
    ret = E_gamma_new;
    }
    else if(run == 57196){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998063;
    ret = E_gamma_new;
    }
    else if(run == 57197){
    double E_scattered = 5.70994 - uncorrected_beam_energy;
    double E_gamma_new = 5.70994 - E_scattered*0.997885;
    ret = E_gamma_new;
    }
    else if(run == 57198){
    double E_scattered = 5.71203 - uncorrected_beam_energy;
    double E_gamma_new = 5.71203 - E_scattered*0.998079;
    ret = E_gamma_new;
    }
    else if(run == 57199){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998087;
    ret = E_gamma_new;
    }
    else if(run == 57200){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998089;
    ret = E_gamma_new;
    }
    else if(run == 57201){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998086;
    ret = E_gamma_new;
    }
    else if(run == 57202){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.998131;
    ret = E_gamma_new;
    }
    else if(run == 57203){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.99812;
    ret = E_gamma_new;
    }
    else if(run == 57204){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998078;
    ret = E_gamma_new;
    }
    else if(run == 57205){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.997864;
    ret = E_gamma_new;
    }
    else if(run == 57206){
    double E_scattered = 5.70935 - uncorrected_beam_energy;
    double E_gamma_new = 5.70935 - E_scattered*0.997868;
    ret = E_gamma_new;
    }
    else if(run == 57207){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.997988;
    ret = E_gamma_new;
    }
    else if(run == 57208){
    double E_scattered = 5.71403 - uncorrected_beam_energy;
    double E_gamma_new = 5.71403 - E_scattered*0.998364;
    ret = E_gamma_new;
    }
    else if(run == 57209){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998012;
    ret = E_gamma_new;
    }
    else if(run == 57210){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.997768;
    ret = E_gamma_new;
    }
    else if(run == 57211){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998084;
    ret = E_gamma_new;
    }
    else if(run == 57212){
    double E_scattered = 5.71078 - uncorrected_beam_energy;
    double E_gamma_new = 5.71078 - E_scattered*0.997945;
    ret = E_gamma_new;
    }
    else if(run == 57213){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.9978;
    ret = E_gamma_new;
    }
    else if(run == 57214){
    double E_scattered = 5.7151 - uncorrected_beam_energy;
    double E_gamma_new = 5.7151 - E_scattered*0.998414;
    ret = E_gamma_new;
    }
    else if(run == 57215){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998085;
    ret = E_gamma_new;
    }
    else if(run == 57216){
    double E_scattered = 5.71533 - uncorrected_beam_energy;
    double E_gamma_new = 5.71533 - E_scattered*0.998309;
    ret = E_gamma_new;
    }
    else if(run == 57217){
    double E_scattered = 5.72218 - uncorrected_beam_energy;
    double E_gamma_new = 5.72218 - E_scattered*0.998731;
    ret = E_gamma_new;
    }
    else if(run == 57218){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.998797;
    ret = E_gamma_new;
    }
    else if(run == 57219){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998827;
    ret = E_gamma_new;
    }
    else if(run == 57220){
    double E_scattered = 5.72218 - uncorrected_beam_energy;
    double E_gamma_new = 5.72218 - E_scattered*0.999001;
    ret = E_gamma_new;
    }
    else if(run == 57221){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.998883;
    ret = E_gamma_new;
    }
    else if(run == 57222){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998858;
    ret = E_gamma_new;
    }
    else if(run == 57223){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998849;
    ret = E_gamma_new;
    }
    else if(run == 57224){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.998893;
    ret = E_gamma_new;
    }
    else if(run == 57225){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998903;
    ret = E_gamma_new;
    }
    else if(run == 57226){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.99883;
    ret = E_gamma_new;
    }
    else if(run == 57227){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.998841;
    ret = E_gamma_new;
    }
    else if(run == 57228){
    double E_scattered = 5.72218 - uncorrected_beam_energy;
    double E_gamma_new = 5.72218 - E_scattered*0.998806;
    ret = E_gamma_new;
    }
    else if(run == 57229){
    double E_scattered = 5.72218 - uncorrected_beam_energy;
    double E_gamma_new = 5.72218 - E_scattered*0.998655;
    ret = E_gamma_new;
    }
    else if(run == 57233){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998793;
    ret = E_gamma_new;
    }
    else if(run == 57234){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998879;
    ret = E_gamma_new;
    }
    else if(run == 57235){
    double E_scattered = 5.72216 - uncorrected_beam_energy;
    double E_gamma_new = 5.72216 - E_scattered*0.998639;
    ret = E_gamma_new;
    }
    else if(run == 57236){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.998816;
    ret = E_gamma_new;
    }
    else if(run == 57237){
    double E_scattered = 5.72215 - uncorrected_beam_energy;
    double E_gamma_new = 5.72215 - E_scattered*0.999086;
    ret = E_gamma_new;
    }
    else if(run == 57238){
    double E_scattered = 5.72216 - uncorrected_beam_energy;
    double E_gamma_new = 5.72216 - E_scattered*0.998827;
    ret = E_gamma_new;
    }
    else if(run == 57239){
    double E_scattered = 5.72216 - uncorrected_beam_energy;
    double E_gamma_new = 5.72216 - E_scattered*0.998186;
    ret = E_gamma_new;
    }
    else if(run == 57241){
    double E_scattered = 5.72217 - uncorrected_beam_energy;
    double E_gamma_new = 5.72217 - E_scattered*0.998256;
    ret = E_gamma_new;
    }
    else if(run == 57249){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.998201;
    ret = E_gamma_new;
    }
    else if(run == 57250){
    double E_scattered = 5.7111 - uncorrected_beam_energy;
    double E_gamma_new = 5.7111 - E_scattered*0.99804;
    ret = E_gamma_new;
    }
    else if(run == 57251){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.998241;
    ret = E_gamma_new;
    }
    else if(run == 57252){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998059;
    ret = E_gamma_new;
    }
    else if(run == 57253){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998105;
    ret = E_gamma_new;
    }
    else if(run == 57255){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998012;
    ret = E_gamma_new;
    }
    else if(run == 57256){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998081;
    ret = E_gamma_new;
    }
    else if(run == 57257){
    double E_scattered = 5.7107 - uncorrected_beam_energy;
    double E_gamma_new = 5.7107 - E_scattered*0.998247;
    ret = E_gamma_new;
    }
    else if(run == 57258){
    double E_scattered = 5.71107 - uncorrected_beam_energy;
    double E_gamma_new = 5.71107 - E_scattered*0.998081;
    ret = E_gamma_new;
    }
    else if(run == 57260){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.997925;
    ret = E_gamma_new;
    }
    else if(run == 57261){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998076;
    ret = E_gamma_new;
    }
    else if(run == 57262){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998401;
    ret = E_gamma_new;
    }
    else if(run == 57263){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998127;
    ret = E_gamma_new;
    }
    else if(run == 57264){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998077;
    ret = E_gamma_new;
    }
    else if(run == 57265){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998082;
    ret = E_gamma_new;
    }
    else if(run == 57266){
    double E_scattered = 5.7111 - uncorrected_beam_energy;
    double E_gamma_new = 5.7111 - E_scattered*0.999738;
    ret = E_gamma_new;
    }
    else if(run == 57267){
    double E_scattered = 5.7111 - uncorrected_beam_energy;
    double E_gamma_new = 5.7111 - E_scattered*0.997869;
    ret = E_gamma_new;
    }
    else if(run == 57268){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998034;
    ret = E_gamma_new;
    }
    else if(run == 57270){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998104;
    ret = E_gamma_new;
    }
    else if(run == 57271){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.99798;
    ret = E_gamma_new;
    }
    else if(run == 57272){
    double E_scattered = 5.71123 - uncorrected_beam_energy;
    double E_gamma_new = 5.71123 - E_scattered*0.997976;
    ret = E_gamma_new;
    }
    else if(run == 57273){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998166;
    ret = E_gamma_new;
    }
    else if(run == 57274){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.997893;
    ret = E_gamma_new;
    }
    else if(run == 57275){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998227;
    ret = E_gamma_new;
    }
    else if(run == 57276){
    double E_scattered = 5.71233 - uncorrected_beam_energy;
    double E_gamma_new = 5.71233 - E_scattered*0.998174;
    ret = E_gamma_new;
    }
    else if(run == 57277){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998186;
    ret = E_gamma_new;
    }
    else if(run == 57278){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998216;
    ret = E_gamma_new;
    }
    else if(run == 57279){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998265;
    ret = E_gamma_new;
    }
    else if(run == 57280){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998103;
    ret = E_gamma_new;
    }
    else if(run == 57281){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998001;
    ret = E_gamma_new;
    }
    else if(run == 57282){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998088;
    ret = E_gamma_new;
    }
    else if(run == 57283){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998216;
    ret = E_gamma_new;
    }
    else if(run == 57284){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998223;
    ret = E_gamma_new;
    }
    else if(run == 57285){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.997984;
    ret = E_gamma_new;
    }
    else if(run == 57286){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.99809;
    ret = E_gamma_new;
    }
    else if(run == 57287){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998023;
    ret = E_gamma_new;
    }
    else if(run == 57288){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998237;
    ret = E_gamma_new;
    }
    else if(run == 57290){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998331;
    ret = E_gamma_new;
    }
    else if(run == 57291){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.99827;
    ret = E_gamma_new;
    }
    else if(run == 57293){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.998115;
    ret = E_gamma_new;
    }
    else if(run == 57294){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*1.00051;
    ret = E_gamma_new;
    }
    else if(run == 57295){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998115;
    ret = E_gamma_new;
    }
    else if(run == 57296){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998091;
    ret = E_gamma_new;
    }
    else if(run == 57297){
    double E_scattered = 5.71088 - uncorrected_beam_energy;
    double E_gamma_new = 5.71088 - E_scattered*0.998152;
    ret = E_gamma_new;
    }
    else if(run == 57298){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998156;
    ret = E_gamma_new;
    }
    else if(run == 57299){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998202;
    ret = E_gamma_new;
    }
    else if(run == 57300){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998142;
    ret = E_gamma_new;
    }
    else if(run == 57301){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998012;
    ret = E_gamma_new;
    }
    else if(run == 57302){
    double E_scattered = 5.71098 - uncorrected_beam_energy;
    double E_gamma_new = 5.71098 - E_scattered*0.998292;
    ret = E_gamma_new;
    }
    else if(run == 57303){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.998118;
    ret = E_gamma_new;
    }
    else if(run == 57304){
    double E_scattered = 5.71095 - uncorrected_beam_energy;
    double E_gamma_new = 5.71095 - E_scattered*0.998038;
    ret = E_gamma_new;
    }
    else if(run == 57305){
    double E_scattered = 5.71105 - uncorrected_beam_energy;
    double E_gamma_new = 5.71105 - E_scattered*0.99812;
    ret = E_gamma_new;
    }
    else if(run == 57306){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.997931;
    ret = E_gamma_new;
    }
    else if(run == 57307){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.997962;
    ret = E_gamma_new;
    }
    else if(run == 57308){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998147;
    ret = E_gamma_new;
    }
    else if(run == 57309){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998134;
    ret = E_gamma_new;
    }
    else if(run == 57310){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.997919;
    ret = E_gamma_new;
    }
    else if(run == 57311){
    double E_scattered = 5.711 - uncorrected_beam_energy;
    double E_gamma_new = 5.711 - E_scattered*0.998107;
    ret = E_gamma_new;
    }
    else if(run == 57312){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998161;
    ret = E_gamma_new;
    }
    else if(run == 57315){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.998106;
    ret = E_gamma_new;
    }
    else if(run == 57316){
    double E_scattered = 5.71097 - uncorrected_beam_energy;
    double E_gamma_new = 5.71097 - E_scattered*0.997952;
    ret = E_gamma_new;
    }
    else if(run == 57317){
    double E_scattered = 5.71103 - uncorrected_beam_energy;
    double E_gamma_new = 5.71103 - E_scattered*0.997881;
    ret = E_gamma_new;
    }
    else {
     throw std::domain_error("Unknown run number");
    }
    return ret;
}


} /* namespace clas::g12 */
} /* namespace clas */






#endif /* __CLAS_G12_ECOR_HPP__ */