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