babymaker  e95a6a9342d4604277fe7cc6149b6b5b24447d89
cross_sections.cc
Go to the documentation of this file.
1 // Lists the MC cross sections
2 
3 #include <iostream>
4 #include "babymaker/bmaker/interface/cross_sections.hh"
5 
6 using namespace std;
7 
8 namespace xsec{
9 
10  float crossSection(const TString &file){
11  float xsec(1.), Htobb(0.5824);
12 
13  if(file.Contains("Run201")) xsec = 1.;
14 
15 
16  //NOTE: When you add more topologies, make sure to also add them to list of things to parse
17  //for model parameters (search "cit" in bmaker_basic.cc)
18  if(file.Contains("SMS")) xsec = 1.0; // Filler to avoid crashes in the SMS scans
19 
20  // From https://twiki.cern.ch/twiki/bin/view/LHCPhysics/SUSYCrossSections13TeVgluglu
21  if(file.Contains("T1tttt")) xsec = 0.0856418; //just to avoid crash on scans
22  if(file.Contains("T1tttt") && file.Contains("825_")) xsec = 1.2167;
23  if(file.Contains("T1tttt") && file.Contains("1025_")) xsec = 0.272778;
24  if(file.Contains("T1tttt") && file.Contains("1150_")) xsec = 0.117687;
25  if(file.Contains("T1tttt") && file.Contains("1200_")) xsec = 0.0856418;
26  if(file.Contains("T1tttt") && file.Contains("1500_")) xsec = 0.0141903;
27 
28  if(file.Contains("T2tt") && file.Contains("650_")) xsec = 0.107045;
29  if(file.Contains("T2tt") && file.Contains("850_")) xsec = 0.0189612;
30 
31  if(file.Contains("TChiHH") && file.Contains("-400")) xsec = 0.0887325*Htobb*Htobb;
32 
33 
34  if(file.Contains("SMS-T2tt_2J_mStop-425_mLSP-325")) xsec = 1.31169;
35  if(file.Contains("SMS-T2tt_2J_mStop-500_mLSP-325")) xsec = 0.51848;
36  if(file.Contains("SMS-T1bbbb_2J_mGl-1500_mLSP-100")) xsec = 0.0141903;
37  if(file.Contains("SMS-T1bbbb_2J_mGl-1000_mLSP-900")) xsec = 0.325388;
38  if(file.Contains("SMS-T1qqqq_2J_mGl-1400_mLSP-100")) xsec = 0.0252977;
39  if(file.Contains("SMS-T1qqqq_2J_mGl-1000_mLSP-800")) xsec = 0.325388;
40  if(file.Contains("SMS-T2bb_2J_mStop-600_mLSP-580")) xsec = 0.174599;
41  if(file.Contains("SMS-T2bb_2J_mStop-900_mLSP-100")) xsec = 0.0128895;
42 
43  if(file.Contains("RPV") && file.Contains("1000")) xsec = 0.325388;
44  if(file.Contains("RPV") && file.Contains("1100")) xsec = 0.163491;
45  if(file.Contains("RPV") && file.Contains("1200")) xsec = 0.0856418;
46  if(file.Contains("RPV") && file.Contains("1300")) xsec = 0.0460525;
47  if(file.Contains("RPV") && file.Contains("1400")) xsec = 0.0252977;
48 
49 
50  // Cross-section taken from https://twiki.cern.ch/twiki/bin/view/LHCPhysics/TtbarNNLO
51  // Alternative option: https://twiki.cern.ch/twiki/bin/view/Sandbox/FullNNLOcrossSections#Top_cross_section_for_13_TeV
52  if(file.Contains("TTJets_Tune") || file.Contains("TT_")) xsec = 815.96;
53  if(file.Contains("TTJets_HT")){//LO cross sections with k-factor of 1.625 already applied
54  if(file.Contains("2500toInf")) xsec = 0.0023234211;
55  if(file.Contains("1200to2500")) xsec = 0.194972521;
56  if(file.Contains("800to1200")) xsec = 1.07722318;
57  if(file.Contains("600to800")) xsec = 2.61537118;
58 
59  }
60  // The efficiency of the mtt>1000 cut is taken from sigma(mtt>1000)/sigma(inclusive) from mcm
61  double mtt_1000_eff=(11.16/670.3);
62  if(file.Contains("TTJets_Mtt-1000toInf")) xsec = 815.96*mtt_1000_eff;
63 
64  if(file.Contains("TTJets_DiLept") || file.Contains("TTTo2L2Nu")) xsec = 85.66; // (3*0.108)^2*815.96
65  if(file.Contains("TTJets_SingleLept")) xsec = 178.7; //(1- ((1-3*0.108)^2+(3*0.108)^2))*815.96*0.5 per half
66  if(file.Contains("TTToSemiLeptonic")) xsec = 357.4;
67 
68  if(file.Contains("TTJets_DiLept_genMET-150")) xsec = 0.06392*85.66; // filter_eff*(3*0.108)^2*815.96
69  if(file.Contains("TTJets_SingleLept") && file.Contains("genMET-150")) xsec = 0.05217*178.7; //filter_eff*(1- ((1-3*0.108)^2+(3*0.108)^2))*815.96*0.5 per half
70 
71  // cross sections from mcm
72  if(file.Contains("TTG")) xsec = 3.697;
73  if(file.Contains("TTTT_TuneCUETP8M1_13TeV-amcatnlo")) xsec = 0.009103;
74  // mcm cross section times the same kfactors used for leptonic samples
75  if(file.Contains("WJetsToQQ_HT-600ToInf")) xsec = 95.14*1.21;
76  if(file.Contains("ZJetsToQQ_HT600toInf")) xsec = 5.67*1.23;
77 
78  // From https://cms-pdmv.cern.ch/mcm
79  // k-factors from https://mangano.web.cern.ch/mangano/public/MECCA/samples_50ns_miniaod.txt
80  // k-factors are ratio of https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat13TeV
81  // NLO/NNLO cross-sections to that of an inclusive sample in mcm at lower order (LO/NLO)
82 
83  if(file.Contains("WJetsToLNu_TuneCUET")) xsec=61526.7; //NNLO from Lesya's summary table
84 
85  //cross-section per slice changed due to change in genHT definition
86  if(file.Contains("WJetsToLNu_HT-70To100")) xsec = 1372.*1.21;
87  if(file.Contains("WJetsToLNu_HT-100To200")) xsec = 1347.*1.21;
88  if(file.Contains("WJetsToLNu_HT-200To400")) xsec = 360.*1.21;
89  if(file.Contains("WJetsToLNu_HT-400To600")) xsec = 48.98*1.21;
90  if(file.Contains("WJetsToLNu_HT-600ToInf")) xsec = 18.77*1.21;
91  if(file.Contains("WJetsToLNu_HT-600To800")) xsec = 12.05*1.21;
92  if(file.Contains("WJetsToLNu_HT-800To1200")) xsec = 5.501*1.21;
93  if(file.Contains("WJetsToLNu_HT-1200To2500")) xsec = 1.329*1.21;
94  if(file.Contains("WJetsToLNu_HT-2500ToInf")) xsec = 0.03216*1.21;
95  if(file.Contains("WToENu")) xsec = 16000.0;
96  if(file.Contains("WToMuNu")) xsec = 16100.0;
97 
98  if(file.Contains("QCD_HT-100To250_13TeV-madgraph")) xsec = 28730000.;
99  if(file.Contains("QCD_HT_250To500_13TeV-madgraph")) xsec = 670500.0;
100  if(file.Contains("QCD_HT-500To1000_13TeV-madgraph")) xsec = 26740.0;
101  if(file.Contains("QCD_HT_1000ToInf_13TeV-madgraph")) xsec = 769.7;
102 
103  if(file.Contains("QCD_HT50to100_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1; // not in MCM, need to derive
104  if(file.Contains("QCD_HT100to200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 27540000;
105  if(file.Contains("QCD_HT200to300_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1735000;
106  if(file.Contains("QCD_HT300to500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 366800;
107  if(file.Contains("QCD_HT500to700_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 29370;
108  if(file.Contains("QCD_HT700to1000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 6524;
109  if(file.Contains("QCD_HT1000to1500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1064;
110  if(file.Contains("QCD_HT1500to2000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 121.5;
111  if(file.Contains("QCD_HT2000toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 25.42;
112 
113  if(file.Contains("QCD_bEnriched_HT100to200_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1321000;
114  if(file.Contains("QCD_bEnriched_HT200to300_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 87950;
115  if(file.Contains("QCD_bEnriched_HT300to500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 17570;
116  if(file.Contains("QCD_bEnriched_HT500to700_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1512;
117  if(file.Contains("QCD_bEnriched_HT700to1000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 320;
118  if(file.Contains("QCD_bEnriched_HT1000to1500_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 50.26;
119  if(file.Contains("QCD_bEnriched_HT1500to2000_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 4.41;
120  if(file.Contains("QCD_bEnriched_HT2000toInf_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 0.78;
121 
122  if(file.Contains("QCD_HT100to200_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1899000;
123  if(file.Contains("QCD_HT200to300_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 156500;
124  if(file.Contains("QCD_HT300to500_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 38970;
125  if(file.Contains("QCD_HT500to700_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 4150;
126  if(file.Contains("QCD_HT700to1000_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 1000;
127  if(file.Contains("QCD_HT1000to1500_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 184.4;
128  if(file.Contains("QCD_HT1500to2000_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 21.31;
129  if(file.Contains("QCD_HT2000toInf_BGenFilter_TuneCUETP8M1_13TeV-madgraphMLM-pythia8")) xsec = 4.16;
130 
131  if(file.Contains("QCD_Pt-1800_")) xsec = 0.1091;
132 
133  // only thing that changed for Spring 15 is "QCD_Pt-" --> "QCD_Pt_"
134  if (file.Contains("QCD_Pt")){
135  if(file.Contains("5to10")) xsec = 80710000000;
136  if(file.Contains("10to15")) xsec = 7528000000;
137  if(file.Contains("15to30")) xsec = 2237000000;
138  if(file.Contains("30to50")) xsec = 161500000;
139  if(file.Contains("50to80")) xsec = 22110000;
140  if(file.Contains("80to120")) xsec = 2762530;// xsec = 3000114.3;
141  if(file.Contains("120to170")) xsec = 471100;//493200;
142  if(file.Contains("170to300")) xsec = 117276;//120300;
143  if(file.Contains("300to470")) xsec = 7823;//7475;
144  if(file.Contains("470to600")) xsec = 648.2;//587.1;
145  if(file.Contains("600to800")) xsec = 186.9;//167;
146  if(file.Contains("800to1000")) xsec = 32.293;//28.25;
147  if(file.Contains("1000to1400")) xsec = 9.4183;//8.195;
148  if(file.Contains("1400to1800")) xsec = 0.84265;// 0.7346;
149  if(file.Contains("1800to2400")) xsec = 0.114943;//0.102;
150  if(file.Contains("2400to3200")) xsec = 0.00682981;//0.00644;
151  if(file.Contains("3200toInf")) xsec = 0.000165445;// 0.000163;
152  }
153 
154  // Cross sections from https://twiki.cern.ch/twiki/bin/view/LHCPhysics/SingleTopRefXsec
155  // multiplied by BF(W->mu,e,tau) = 0.324
156  if (file.Contains("ST_s-channel_4f_leptonDecays_13TeV-amcatnlo-pythia8")) xsec = 3.34;
157  if (file.Contains("ST_t-channel_antitop_4f_leptonDecays_13TeV-powheg-pythia8")) xsec = 26.23;
158  if (file.Contains("ST_t-channel_top_4f_leptonDecays_13TeV-powheg-pythia8")) xsec = 44.07;
159  if (file.Contains("ST_t-channel_antitop_4f_inclusiveDecays_TuneCUETP8M2T4_13TeV-powhegV2")) xsec = 80.95;
160  if (file.Contains("ST_t-channel_top_4f_inclusiveDecays_TuneCUETP8M2T4_13TeV-powhegV2")) xsec = 136.02;
161  if (file.Contains("ST_tW_antitop_5f_inclusiveDecays_13TeV-powheg-pythia8")) xsec = 35.85;
162  if (file.Contains("ST_tW_top_5f_inclusiveDecays_13TeV-powheg-pythia8")) xsec = 35.85;
163  if (file.Contains("ST_tW_antitop_5f_NoFullyHadronicDecays_13TeV-powheg")) xsec = 35.85*0.543;
164  if (file.Contains("ST_tW_top_5f_NoFullyHadronicDecays_13TeV-powheg")) xsec = 35.85*0.543;
165 
166  if(file.Contains("DYJetsToLL_M-10to50_Tune")) xsec = 18610*1.23;
167  if(file.Contains("DYJetsToLL_M-50_TuneCUETP8M1_13TeV")) xsec = 4895*1.23;
168 
169  if(file.Contains("DYJetsToLL_M-50_HT-70to100")) xsec = 175.3*1.23;
170  if(file.Contains("DYJetsToLL_M-50_HT-100to200")) xsec = 139.4*1.23;
171  if(file.Contains("DYJetsToLL_M-50_HT-200to400")) xsec = 42.75*1.23;
172  if(file.Contains("DYJetsToLL_M-50_HT-400to600")) xsec = 5.497*1.23;
173  if(file.Contains("DYJetsToLL_M-50_HT-600to800")) xsec = 1.363*1.23;
174  if(file.Contains("DYJetsToLL_M-50_HT-800to1200")) xsec = 0.6759*1.23;
175  if(file.Contains("DYJetsToLL_M-50_HT-1200to2500")) xsec = 0.116*1.23;
176  if(file.Contains("DYJetsToLL_M-50_HT-2500toInf")) xsec = 0.002592*1.23;
177  if(file.Contains("DYJetsToLL_M-50_HT-600toInf")) xsec = 2.21*1.23;
178 
179  if(file.Contains("DYJetsToLL_Pt-50To100")) xsec = 354.3;
180  if(file.Contains("DYJetsToLL_Pt-100To250")) xsec = 83.12;
181  if(file.Contains("DYJetsToLL_Pt-250To400")) xsec = 3.047;
182  if(file.Contains("DYJetsToLL_Pt-400To650")) xsec = 0.3921;
183  if(file.Contains("DYJetsToLL_Pt-650ToInf")) xsec = 0.03636;
184 
185  if(file.Contains("ZJetsToNuNu_HT-100To200_")) xsec = 280.35*1.27;
186  if(file.Contains("ZJetsToNuNu_HT-200To400_")) xsec = 77.67*1.27;
187  if(file.Contains("ZJetsToNuNu_HT-400To600_")) xsec = 10.73*1.27;
188  if(file.Contains("ZJetsToNuNu_HT-600To800_")) xsec = 2.536*1.27;
189  if(file.Contains("ZJetsToNuNu_HT-800To1200_")) xsec = 1.161*1.27;
190  if(file.Contains("ZJetsToNuNu_HT-1200To2500_")) xsec = 0.2824*1.27;
191  if(file.Contains("ZJetsToNuNu_HT-2500ToInf_")) xsec = 0.006459*1.27;
192  if(file.Contains("ZJetsToNuNu_HT-600ToInf_")) xsec = 3.986*1.27;
193 
194  if(file.Contains("TTZJets_Tune4C_13TeV-madgraph-tauola")) xsec = 0.7598;
195  if(file.Contains("TTWJets_Tune4C_13TeV-madgraph-tauola")) xsec = 0.5662;
196  if(file.Contains("TTZToQQ")) xsec = 0.5297;
197  if(file.Contains("TTZToLLNuNu_M-10")) xsec = 0.2529;
198  if(file.Contains("TTWJetsToQQ")) xsec = 0.4062;
199  if(file.Contains("TTWJetsToLNu")) xsec = 0.2043;
200 
201  //https://twiki.cern.ch/twiki/bin/viewauth/CMS/SummaryTable1G25ns#Diboson
202  if(file.Contains("WWTo2L2Nu_13TeV-powheg")) xsec = 12.178; //NNLO
203  if(file.Contains("WWToLNuQQ_13TeV-powheg")) xsec = 49.997; //NNLO
204  if(file.Contains("ttHJetTobb_M125_13TeV_amcatnloFXFX")) xsec = 0.2934;
205 
206  if(file.Contains("WZTo1L3Nu_13TeV_amcatnloFXFX")) xsec = 3.05;
207  if(file.Contains("WZTo1L1Nu2Q_13TeV_amcatnloFXFX")) xsec = 10.96;
208  if(file.Contains("WZTo2L2Q_13TeV_amcatnloFXFX")) xsec = 5.595;
209  if(file.Contains("WZTo3LNu_TuneCUETP8M1_13TeV-powheg-")) xsec = 4.42965;
210  if(file.Contains("VVTo2L2Nu_13TeV_amcatnloFXFX")) xsec = 11.95;
211  if(file.Contains("ZZ_TuneCUETP8M1")) xsec = 16.523;
212 
213  // Calculated at 13 TeV in
214  // https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageAt1314TeV
215  // Higgs branching ratios from
216  // https://twiki.cern.ch/twiki/bin/view/LHCPhysics/CERNYellowReportPageBR
217  if(file.Contains("ZH_HToBB_ZToLL_M-125_13TeV")) xsec = 0.883*Htobb*0.033658;
218  if(file.Contains("ZH_HToBB_ZToNuNu_M-125_13TeV")) xsec = 0.883*Htobb*0.2;
219  if(file.Contains("WH_HToBB_WToLNu_M-125_13TeV")) xsec = 1.373*Htobb*(0.1071+0.1063+0.1138);
220  if(file.Contains("ZH_HToBB_ZToNuNu_M125_13TeV")) xsec = 0.883*Htobb*0.2;
221  if(file.Contains("WH_HToBB_WToLNu_M125_13TeV")) xsec = 1.373*Htobb*(0.1071+0.1063+0.1138);
222 
223  if(xsec<=0) std::cout<<"BABYMAKER: Cross section not found for "<<file<<std::endl;
224 
225  return xsec;
226  }
227 
228  float fractionNegWeights(const TString &file){
229  float fneg(0.);
230 
231  // ttH, ttZ, ttW, ttGamma
232  if(file.Contains("ttHJetTobb_M125_13TeV_amcatnloFXFX")) fneg = 0.3515;
233  if(file.Contains("TTZToQQ")) fneg = 0.2657;
234  if(file.Contains("TTZToLLNuNu_M-10")) fneg = 0.2676;
235  if(file.Contains("TTWJetsToQQ")) fneg = 0.2412;
236  if(file.Contains("TTWJetsToLNu")) fneg = 0.2433;
237  if(file.Contains("TTG")) fneg = 0.342; // from MCM
238 
239  if(file.Contains("TTTT_TuneCUETP8M1_13TeV-amcatnlo")) fneg = 0.41; // from MCM
240  if(file.Contains("VVTo2L2Nu_13TeV_amcatnloFXFX")) fneg = 0.20; // from MCM
241  if(file.Contains("TTJets_Mtt-1000toInf")) fneg = 0.376996;
242 
243  // Single top
244  if (file.Contains("ST_s-channel_4f_leptonDecays_13TeV-amcatnlo-pythia8")) fneg = 0.1884;
245 
246  return fneg;
247  }
248 
249 
250  void signalCrossSection(int glu_mass, double &xsec, double &xsec_unc){
251  if (glu_mass == 595) { xsec = 0.; xsec_unc = 0.; return; } // we shouldn't have these points
252  else if (glu_mass == 600) { xsec = 9.20353; xsec_unc = 0.137185; return; }
253  else if (glu_mass == 605) { xsec = 8.74315; xsec_unc = 0.137502; return; }
254  else if (glu_mass == 610) { xsec = 8.30988; xsec_unc = 0.136818; return; }
255  else if (glu_mass == 615) { xsec = 7.9012; xsec_unc = 0.137122; return; }
256  else if (glu_mass == 620) { xsec = 7.51811; xsec_unc = 0.136123; return; }
257  else if (glu_mass == 625) { xsec = 7.15194; xsec_unc = 0.136874; return; }
258  else if (glu_mass == 630) { xsec = 6.80558; xsec_unc = 0.136655; return; }
259  else if (glu_mass == 635) { xsec = 6.47541; xsec_unc = 0.136459; return; }
260  else if (glu_mass == 640) { xsec = 6.17196; xsec_unc = 0.136449; return; }
261  else if (glu_mass == 645) { xsec = 5.87366; xsec_unc = 0.136392; return; }
262  else if (glu_mass == 650) { xsec = 5.60048; xsec_unc = 0.136262; return; }
263  else if (glu_mass == 655) { xsec = 5.33799; xsec_unc = 0.137031; return; }
264  else if (glu_mass == 660) { xsec = 5.09822; xsec_unc = 0.137329; return; }
265  else if (glu_mass == 665) { xsec = 4.86409; xsec_unc = 0.137702; return; }
266  else if (glu_mass == 670) { xsec = 4.64349; xsec_unc = 0.138022; return; }
267  else if (glu_mass == 675) { xsec = 4.43132; xsec_unc = 0.138749; return; }
268  else if (glu_mass == 680) { xsec = 4.23046; xsec_unc = 0.139166; return; }
269  else if (glu_mass == 685) { xsec = 4.03841; xsec_unc = 0.139934; return; }
270  else if (glu_mass == 690) { xsec = 3.85666; xsec_unc = 0.139917; return; }
271  else if (glu_mass == 695) { xsec = 3.68567; xsec_unc = 0.140759; return; }
272  else if (glu_mass == 700) { xsec = 3.5251; xsec_unc = 0.141034; return; }
273  else if (glu_mass == 705) { xsec = 3.3737; xsec_unc = 0.141609; return; }
274  else if (glu_mass == 710) { xsec = 3.22336; xsec_unc = 0.141972; return; }
275  else if (glu_mass == 715) { xsec = 3.0811; xsec_unc = 0.142311; return; }
276  else if (glu_mass == 720) { xsec = 2.9509; xsec_unc = 0.142518; return; }
277  else if (glu_mass == 725) { xsec = 2.81957; xsec_unc = 0.143333; return; }
278  else if (glu_mass == 730) { xsec = 2.7; xsec_unc = 0.143772; return; }
279  else if (glu_mass == 735) { xsec = 2.57737; xsec_unc = 0.144452; return; }
280  else if (glu_mass == 740) { xsec = 2.47729; xsec_unc = 0.144485; return; }
281  else if (glu_mass == 745) { xsec = 2.3661; xsec_unc = 0.145381; return; }
282  else if (glu_mass == 750) { xsec = 2.26585; xsec_unc = 0.145653; return; }
283  else if (glu_mass == 755) { xsec = 2.17436; xsec_unc = 0.145861; return; }
284  else if (glu_mass == 760) { xsec = 2.08446; xsec_unc = 0.146279; return; }
285  else if (glu_mass == 765) { xsec = 1.99341; xsec_unc = 0.147278; return; }
286  else if (glu_mass == 770) { xsec = 1.91352; xsec_unc = 0.147424; return; }
287  else if (glu_mass == 775) { xsec = 1.83188; xsec_unc = 0.147835; return; }
288  else if (glu_mass == 780) { xsec = 1.76145; xsec_unc = 0.148078; return; }
289  else if (glu_mass == 785) { xsec = 1.68078; xsec_unc = 0.148956; return; }
290  else if (glu_mass == 790) { xsec = 1.62071; xsec_unc = 0.149017; return; }
291  else if (glu_mass == 795) { xsec = 1.54896; xsec_unc = 0.149976; return; }
292  else if (glu_mass == 800) { xsec = 1.4891; xsec_unc = 0.150167; return; }
293  else if (glu_mass == 805) { xsec = 1.42888; xsec_unc = 0.150599; return; }
294  else if (glu_mass == 810) { xsec = 1.36759; xsec_unc = 0.151122; return; }
295  else if (glu_mass == 815) { xsec = 1.31749; xsec_unc = 0.151184; return; }
296  else if (glu_mass == 820) { xsec = 1.26659; xsec_unc = 0.151928; return; }
297  else if (glu_mass == 825) { xsec = 1.2167; xsec_unc = 0.152141; return; }
298  else if (glu_mass == 830) { xsec = 1.16617; xsec_unc = 0.152437; return; }
299  else if (glu_mass == 835) { xsec = 1.12555; xsec_unc = 0.153009; return; }
300  else if (glu_mass == 840) { xsec = 1.07523; xsec_unc = 0.15367; return; }
301  else if (glu_mass == 845) { xsec = 1.03426; xsec_unc = 0.154018; return; }
302  else if (glu_mass == 850) { xsec = 0.996137; xsec_unc = 0.154252; return; }
303  else if (glu_mass == 855) { xsec = 0.957975; xsec_unc = 0.154597; return; }
304  else if (glu_mass == 860) { xsec = 0.921447; xsec_unc = 0.155362; return; }
305  else if (glu_mass == 865) { xsec = 0.885917; xsec_unc = 0.155643; return; }
306  else if (glu_mass == 870) { xsec = 0.852433; xsec_unc = 0.156368; return; }
307  else if (glu_mass == 875) { xsec = 0.820259; xsec_unc = 0.156742; return; }
308  else if (glu_mass == 880) { xsec = 0.788789; xsec_unc = 0.156746; return; }
309  else if (glu_mass == 885) { xsec = 0.759346; xsec_unc = 0.157507; return; }
310  else if (glu_mass == 890) { xsec = 0.731213; xsec_unc = 0.157879; return; }
311  else if (glu_mass == 895) { xsec = 0.703532; xsec_unc = 0.158276; return; }
312  else if (glu_mass == 900) { xsec = 0.677478; xsec_unc = 0.158762; return; }
313  else if (glu_mass == 905) { xsec = 0.652317; xsec_unc = 0.15914; return; }
314  else if (glu_mass == 910) { xsec = 0.627695; xsec_unc = 0.159569; return; }
315  else if (glu_mass == 915) { xsec = 0.605596; xsec_unc = 0.159838; return; }
316  else if (glu_mass == 920) { xsec = 0.58302; xsec_unc = 0.16029; return; }
317  else if (glu_mass == 925) { xsec = 0.561889; xsec_unc = 0.160626; return; }
318  else if (glu_mass == 930) { xsec = 0.540533; xsec_unc = 0.161499; return; }
319  else if (glu_mass == 935) { xsec = 0.521159; xsec_unc = 0.161607; return; }
320  else if (glu_mass == 940) { xsec = 0.501865; xsec_unc = 0.16245; return; }
321  else if (glu_mass == 945) { xsec = 0.483546; xsec_unc = 0.162492; return; }
322  else if (glu_mass == 950) { xsec = 0.466352; xsec_unc = 0.163378; return; }
323  else if (glu_mass == 955) { xsec = 0.45012; xsec_unc = 0.163303; return; }
324  else if (glu_mass == 960) { xsec = 0.433842; xsec_unc = 0.164161; return; }
325  else if (glu_mass == 965) { xsec = 0.418744; xsec_unc = 0.164473; return; }
326  else if (glu_mass == 970) { xsec = 0.403514; xsec_unc = 0.164538; return; }
327  else if (glu_mass == 975) { xsec = 0.389266; xsec_unc = 0.165308; return; }
328  else if (glu_mass == 980) { xsec = 0.375053; xsec_unc = 0.165398; return; }
329  else if (glu_mass == 985) { xsec = 0.36182; xsec_unc = 0.16619; return; }
330  else if (glu_mass == 990) { xsec = 0.349764; xsec_unc = 0.166462; return; }
331  else if (glu_mass == 995) { xsec = 0.337454; xsec_unc = 0.166888; return; }
332  else if (glu_mass == 1000) { xsec = 0.325388; xsec_unc = 0.16758; return; }
333  else if (glu_mass == 1005) { xsec = 0.314329; xsec_unc = 0.167865; return; }
334  else if (glu_mass == 1010) { xsec = 0.30314; xsec_unc = 0.168766; return; }
335  else if (glu_mass == 1015) { xsec = 0.292987; xsec_unc = 0.168793; return; }
336  else if (glu_mass == 1020) { xsec = 0.282927; xsec_unc = 0.169098; return; }
337  else if (glu_mass == 1025) { xsec = 0.272778; xsec_unc = 0.169917; return; }
338  else if (glu_mass == 1030) { xsec = 0.263724; xsec_unc = 0.170244; return; }
339  else if (glu_mass == 1035) { xsec = 0.254721; xsec_unc = 0.170758; return; }
340  else if (glu_mass == 1040) { xsec = 0.245426; xsec_unc = 0.171325; return; }
341  else if (glu_mass == 1045) { xsec = 0.237403; xsec_unc = 0.171542; return; }
342  else if (glu_mass == 1050) { xsec = 0.229367; xsec_unc = 0.171975; return; }
343  else if (glu_mass == 1055) { xsec = 0.221273; xsec_unc = 0.172482; return; }
344  else if (glu_mass == 1060) { xsec = 0.214167; xsec_unc = 0.173167; return; }
345  else if (glu_mass == 1065) { xsec = 0.207025; xsec_unc = 0.173211; return; }
346  else if (glu_mass == 1070) { xsec = 0.199967; xsec_unc = 0.173603; return; }
347  else if (glu_mass == 1075) { xsec = 0.193881; xsec_unc = 0.174329; return; }
348  else if (glu_mass == 1080) { xsec = 0.186836; xsec_unc = 0.174816; return; }
349  else if (glu_mass == 1085) { xsec = 0.180783; xsec_unc = 0.175245; return; }
350  else if (glu_mass == 1090) { xsec = 0.174652; xsec_unc = 0.175336; return; }
351  else if (glu_mass == 1095) { xsec = 0.168526; xsec_unc = 0.176231; return; }
352  else if (glu_mass == 1100) { xsec = 0.163491; xsec_unc = 0.176402; return; }
353  else if (glu_mass == 1105) { xsec = 0.158451; xsec_unc = 0.176564; return; }
354  else if (glu_mass == 1110) { xsec = 0.153298; xsec_unc = 0.177266; return; }
355  else if (glu_mass == 1115) { xsec = 0.148246; xsec_unc = 0.177755; return; }
356  else if (glu_mass == 1120) { xsec = 0.143169; xsec_unc = 0.17813; return; }
357  else if (glu_mass == 1125) { xsec = 0.139009; xsec_unc = 0.178569; return; }
358  else if (glu_mass == 1130) { xsec = 0.133972; xsec_unc = 0.179205; return; }
359  else if (glu_mass == 1135) { xsec = 0.129938; xsec_unc = 0.17938; return; }
360  else if (glu_mass == 1140) { xsec = 0.125799; xsec_unc = 0.179658; return; }
361  else if (glu_mass == 1145) { xsec = 0.121755; xsec_unc = 0.180222; return; }
362  else if (glu_mass == 1150) { xsec = 0.117687; xsec_unc = 0.180655; return; }
363  else if (glu_mass == 1155) { xsec = 0.11358; xsec_unc = 0.181327; return; }
364  else if (glu_mass == 1160) { xsec = 0.110557; xsec_unc = 0.181465; return; }
365  else if (glu_mass == 1165) { xsec = 0.107532; xsec_unc = 0.181655; return; }
366  else if (glu_mass == 1170) { xsec = 0.10339; xsec_unc = 0.182421; return; }
367  else if (glu_mass == 1175) { xsec = 0.10036; xsec_unc = 0.182686; return; }
368  else if (glu_mass == 1180) { xsec = 0.0971485; xsec_unc = 0.183142; return; }
369  else if (glu_mass == 1185) { xsec = 0.0942072; xsec_unc = 0.183623; return; }
370  else if (glu_mass == 1190) { xsec = 0.0912756; xsec_unc = 0.183957; return; }
371  else if (glu_mass == 1195) { xsec = 0.0883712; xsec_unc = 0.184467; return; }
372  else if (glu_mass == 1200) { xsec = 0.0856418; xsec_unc = 0.184814; return; }
373  else if (glu_mass == 1205) { xsec = 0.0830236; xsec_unc = 0.185276; return; }
374  else if (glu_mass == 1210) { xsec = 0.0804313; xsec_unc = 0.185714; return; }
375  else if (glu_mass == 1215) { xsec = 0.0779039; xsec_unc = 0.186096; return; }
376  else if (glu_mass == 1220) { xsec = 0.0755801; xsec_unc = 0.186429; return; }
377  else if (glu_mass == 1225) { xsec = 0.0732255; xsec_unc = 0.187227; return; }
378  else if (glu_mass == 1230) { xsec = 0.0709683; xsec_unc = 0.187266; return; }
379  else if (glu_mass == 1235) { xsec = 0.0688462; xsec_unc = 0.187544; return; }
380  else if (glu_mass == 1240) { xsec = 0.0666928; xsec_unc = 0.188404; return; }
381  else if (glu_mass == 1245) { xsec = 0.0646423; xsec_unc = 0.188414; return; }
382  else if (glu_mass == 1250) { xsec = 0.0627027; xsec_unc = 0.189328; return; }
383  else if (glu_mass == 1255) { xsec = 0.0607803; xsec_unc = 0.189693; return; }
384  else if (glu_mass == 1260) { xsec = 0.0589319; xsec_unc = 0.189695; return; }
385  else if (glu_mass == 1265) { xsec = 0.0571859; xsec_unc = 0.190561; return; }
386  else if (glu_mass == 1270) { xsec = 0.0554225; xsec_unc = 0.191806; return; }
387  else if (glu_mass == 1275) { xsec = 0.0536906; xsec_unc = 0.192452; return; }
388  else if (glu_mass == 1280) { xsec = 0.052051; xsec_unc = 0.192396; return; }
389  else if (glu_mass == 1285) { xsec = 0.0504982; xsec_unc = 0.193577; return; }
390  else if (glu_mass == 1290) { xsec = 0.0489404; xsec_unc = 0.194903; return; }
391  else if (glu_mass == 1295) { xsec = 0.047474; xsec_unc = 0.195871; return; }
392  else if (glu_mass == 1300) { xsec = 0.0460525; xsec_unc = 0.1964; return; }
393  else if (glu_mass == 1305) { xsec = 0.0447038; xsec_unc = 0.197627; return; }
394  else if (glu_mass == 1310) { xsec = 0.0433373; xsec_unc = 0.198601; return; }
395  else if (glu_mass == 1315) { xsec = 0.0420362; xsec_unc = 0.198634; return; }
396  else if (glu_mass == 1320) { xsec = 0.0407723; xsec_unc = 0.199586; return; }
397  else if (glu_mass == 1325) { xsec = 0.0395728; xsec_unc = 0.19951; return; }
398  else if (glu_mass == 1330) { xsec = 0.0383587; xsec_unc = 0.19993; return; }
399  else if (glu_mass == 1335) { xsec = 0.0372043; xsec_unc = 0.201012; return; }
400  else if (glu_mass == 1340) { xsec = 0.0361694; xsec_unc = 0.202191; return; }
401  else if (glu_mass == 1345) { xsec = 0.0350586; xsec_unc = 0.201714; return; }
402  else if (glu_mass == 1350) { xsec = 0.0340187; xsec_unc = 0.203088; return; }
403  else if (glu_mass == 1355) { xsec = 0.0330251; xsec_unc = 0.202807; return; }
404  else if (glu_mass == 1360) { xsec = 0.0320787; xsec_unc = 0.203682; return; }
405  else if (glu_mass == 1365) { xsec = 0.0311325; xsec_unc = 0.205466; return; }
406  else if (glu_mass == 1370) { xsec = 0.0302294; xsec_unc = 0.204724; return; }
407  else if (glu_mass == 1375) { xsec = 0.0292919; xsec_unc = 0.206217; return; }
408  else if (glu_mass == 1380) { xsec = 0.0284627; xsec_unc = 0.207773; return; }
409  else if (glu_mass == 1385) { xsec = 0.0276679; xsec_unc = 0.206729; return; }
410  else if (glu_mass == 1390) { xsec = 0.0268339; xsec_unc = 0.208251; return; }
411  else if (glu_mass == 1395) { xsec = 0.0260313; xsec_unc = 0.207488; return; }
412  else if (glu_mass == 1400) { xsec = 0.0252977; xsec_unc = 0.209163; return; }
413  else if (glu_mass == 1405) { xsec = 0.0245679; xsec_unc = 0.210704; return; }
414  else if (glu_mass == 1410) { xsec = 0.0238741; xsec_unc = 0.209586; return; }
415  else if (glu_mass == 1415) { xsec = 0.0231433; xsec_unc = 0.211204; return; }
416  else if (glu_mass == 1420) { xsec = 0.0225194; xsec_unc = 0.212481; return; }
417  else if (glu_mass == 1425) { xsec = 0.0218959; xsec_unc = 0.214183; return; }
418  else if (glu_mass == 1430) { xsec = 0.0211928; xsec_unc = 0.21365; return; }
419  else if (glu_mass == 1435) { xsec = 0.0206244; xsec_unc = 0.217574; return; }
420  else if (glu_mass == 1440) { xsec = 0.0200458; xsec_unc = 0.216629; return; }
421  else if (glu_mass == 1445) { xsec = 0.0194648; xsec_unc = 0.215531; return; }
422  else if (glu_mass == 1450) { xsec = 0.0188887; xsec_unc = 0.219548; return; }
423  else if (glu_mass == 1455) { xsec = 0.018364; xsec_unc = 0.221266; return; }
424  else if (glu_mass == 1460) { xsec = 0.0178858; xsec_unc = 0.220054; return; }
425  else if (glu_mass == 1465) { xsec = 0.0173622; xsec_unc = 0.221916; return; }
426  else if (glu_mass == 1470) { xsec = 0.0168403; xsec_unc = 0.223972; return; }
427  else if (glu_mass == 1475) { xsec = 0.0163556; xsec_unc = 0.222173; return; }
428  else if (glu_mass == 1480) { xsec = 0.0159386; xsec_unc = 0.223581; return; }
429  else if (glu_mass == 1485) { xsec = 0.0154568; xsec_unc = 0.222281; return; }
430  else if (glu_mass == 1490) { xsec = 0.0150345; xsec_unc = 0.224111; return; }
431  else if (glu_mass == 1495) { xsec = 0.0146102; xsec_unc = 0.225293; return; }
432  else if (glu_mass == 1500) { xsec = 0.0141903; xsec_unc = 0.227296; return; }
433  else if (glu_mass == 1505) { xsec = 0.01377; xsec_unc = 0.229402; return; }
434  else if (glu_mass == 1510) { xsec = 0.0133923; xsec_unc = 0.226528; return; }
435  else if (glu_mass == 1515) { xsec = 0.0130286; xsec_unc = 0.232697; return; }
436  else if (glu_mass == 1520) { xsec = 0.012649; xsec_unc = 0.230194; return; }
437  else if (glu_mass == 1525) { xsec = 0.0123374; xsec_unc = 0.231801; return; }
438  else if (glu_mass == 1530) { xsec = 0.0119628; xsec_unc = 0.229449; return; }
439  else if (glu_mass == 1535) { xsec = 0.0116378; xsec_unc = 0.231293; return; }
440  else if (glu_mass == 1540) { xsec = 0.0113183; xsec_unc = 0.233535; return; }
441  else if (glu_mass == 1545) { xsec = 0.0110039; xsec_unc = 0.23456; return; }
442  else if (glu_mass == 1550) { xsec = 0.0107027; xsec_unc = 0.234971; return; }
443  else if (glu_mass == 1555) { xsec = 0.0103967; xsec_unc = 0.23505; return; }
444  else if (glu_mass == 1560) { xsec = 0.0101149; xsec_unc = 0.236723; return; }
445  else if (glu_mass == 1565) { xsec = 0.00984079; xsec_unc = 0.237486; return; }
446  else if (glu_mass == 1570) { xsec = 0.00956216; xsec_unc = 0.238011; return; }
447  else if (glu_mass == 1575) { xsec = 0.00930893; xsec_unc = 0.238712; return; }
448  else if (glu_mass == 1580) { xsec = 0.00905112; xsec_unc = 0.239145; return; }
449  else if (glu_mass == 1585) { xsec = 0.00880102; xsec_unc = 0.24088; return; }
450  else if (glu_mass == 1590) { xsec = 0.00856388; xsec_unc = 0.241033; return; }
451  else if (glu_mass == 1595) { xsec = 0.00832287; xsec_unc = 0.242052; return; }
452  else if (glu_mass == 1600) { xsec = 0.00810078; xsec_unc = 0.242679; return; }
453  else if (glu_mass == 1605) { xsec = 0.0078785; xsec_unc = 0.243322; return; }
454  else if (glu_mass == 1610) { xsec = 0.00767087; xsec_unc = 0.244839; return; }
455  else if (glu_mass == 1615) { xsec = 0.00745579; xsec_unc = 0.245137; return; }
456  else if (glu_mass == 1620) { xsec = 0.00725443; xsec_unc = 0.24569; return; }
457  else if (glu_mass == 1625) { xsec = 0.00705942; xsec_unc = 0.246853; return; }
458  else if (glu_mass == 1630) { xsec = 0.00687457; xsec_unc = 0.24804; return; }
459  else if (glu_mass == 1635) { xsec = 0.00668418; xsec_unc = 0.248672; return; }
460  else if (glu_mass == 1640) { xsec = 0.00651001; xsec_unc = 0.249776; return; }
461  else if (glu_mass == 1645) { xsec = 0.00633268; xsec_unc = 0.250679; return; }
462  else if (glu_mass == 1650) { xsec = 0.00616072; xsec_unc = 0.25138; return; }
463  else if (glu_mass == 1655) { xsec = 0.00599673; xsec_unc = 0.252591; return; }
464  else if (glu_mass == 1660) { xsec = 0.00583243; xsec_unc = 0.253829; return; }
465  else if (glu_mass == 1665) { xsec = 0.00567868; xsec_unc = 0.255006; return; }
466  else if (glu_mass == 1670) { xsec = 0.00553066; xsec_unc = 0.255203; return; }
467  else if (glu_mass == 1675) { xsec = 0.00538094; xsec_unc = 0.255439; return; }
468  else if (glu_mass == 1680) { xsec = 0.00523764; xsec_unc = 0.256602; return; }
469  else if (glu_mass == 1685) { xsec = 0.00509647; xsec_unc = 0.258745; return; }
470  else if (glu_mass == 1690) { xsec = 0.0049577; xsec_unc = 0.258847; return; }
471  else if (glu_mass == 1695) { xsec = 0.00483094; xsec_unc = 0.260944; return; }
472  else if (glu_mass == 1700) { xsec = 0.00470323; xsec_unc = 0.261021; return; }
473  else if (glu_mass == 1705) { xsec = 0.0045807; xsec_unc = 0.262095; return; }
474  else if (glu_mass == 1710) { xsec = 0.00445824; xsec_unc = 0.263238; return; }
475  else if (glu_mass == 1715) { xsec = 0.0043369; xsec_unc = 0.263092; return; }
476  else if (glu_mass == 1720) { xsec = 0.00422488; xsec_unc = 0.264093; return; }
477  else if (glu_mass == 1725) { xsec = 0.00411276; xsec_unc = 0.26513; return; }
478  else if (glu_mass == 1730) { xsec = 0.00400698; xsec_unc = 0.267386; return; }
479  else if (glu_mass == 1735) { xsec = 0.00389655; xsec_unc = 0.267109; return; }
480  else if (glu_mass == 1740) { xsec = 0.00379497; xsec_unc = 0.268072; return; }
481  else if (glu_mass == 1745) { xsec = 0.00370003; xsec_unc = 0.2704; return; }
482  else if (glu_mass == 1750) { xsec = 0.00359842; xsec_unc = 0.271502; return; }
483  else if (glu_mass == 1755) { xsec = 0.00350486; xsec_unc = 0.27229; return; }
484  else if (glu_mass == 1760) { xsec = 0.00341375; xsec_unc = 0.273209; return; }
485  else if (glu_mass == 1765) { xsec = 0.00332255; xsec_unc = 0.27416; return; }
486  else if (glu_mass == 1770) { xsec = 0.00323809; xsec_unc = 0.276458; return; }
487  else if (glu_mass == 1775) { xsec = 0.00314866; xsec_unc = 0.275834; return; }
488  else if (glu_mass == 1780) { xsec = 0.00306841; xsec_unc = 0.276481; return; }
489  else if (glu_mass == 1785) { xsec = 0.00298808; xsec_unc = 0.277145; return; }
490  else if (glu_mass == 1790) { xsec = 0.00291365; xsec_unc = 0.279548; return; }
491  else if (glu_mass == 1795) { xsec = 0.0028312; xsec_unc = 0.280642; return; }
492  else if (glu_mass == 1800) { xsec = 0.00276133; xsec_unc = 0.28108; return; }
493  else if (glu_mass == 1805) { xsec = 0.00269156; xsec_unc = 0.281566; return; }
494  else if (glu_mass == 1810) { xsec = 0.00262156; xsec_unc = 0.282017; return; }
495  else if (glu_mass == 1815) { xsec = 0.00254938; xsec_unc = 0.282755; return; }
496  else if (glu_mass == 1820) { xsec = 0.00248581; xsec_unc = 0.285102; return; }
497  else if (glu_mass == 1825) { xsec = 0.00241549; xsec_unc = 0.285869; return; }
498  else if (glu_mass == 1830) { xsec = 0.00235625; xsec_unc = 0.286103; return; }
499  else if (glu_mass == 1835) { xsec = 0.00229576; xsec_unc = 0.28596; return; }
500  else if (glu_mass == 1840) { xsec = 0.00223603; xsec_unc = 0.286654; return; }
501  else if (glu_mass == 1845) { xsec = 0.00218302; xsec_unc = 0.288949; return; }
502  else if (glu_mass == 1850) { xsec = 0.00212345; xsec_unc = 0.289167; return; }
503  else if (glu_mass == 1855) { xsec = 0.00207; xsec_unc = 0.291835; return; }
504  else if (glu_mass == 1860) { xsec = 0.00200972; xsec_unc = 0.291901; return; }
505  else if (glu_mass == 1865) { xsec = 0.00196025; xsec_unc = 0.292103; return; }
506  else if (glu_mass == 1870) { xsec = 0.00191132; xsec_unc = 0.291893; return; }
507  else if (glu_mass == 1875) { xsec = 0.00185789; xsec_unc = 0.294928; return; }
508  else if (glu_mass == 1880) { xsec = 0.00181527; xsec_unc = 0.29723; return; }
509  else if (glu_mass == 1885) { xsec = 0.00176658; xsec_unc = 0.297236; return; }
510  else if (glu_mass == 1890) { xsec = 0.00172274; xsec_unc = 0.299813; return; }
511  else if (glu_mass == 1895) { xsec = 0.00167806; xsec_unc = 0.296455; return; }
512  else if (glu_mass == 1900) { xsec = 0.00163547; xsec_unc = 0.299045; return; }
513  else if (glu_mass == 1905) { xsec = 0.0015925; xsec_unc = 0.302039; return; }
514  else if (glu_mass == 1910) { xsec = 0.00155445; xsec_unc = 0.301015; return; }
515  else if (glu_mass == 1915) { xsec = 0.00151503; xsec_unc = 0.300356; return; }
516  else if (glu_mass == 1920) { xsec = 0.00147199; xsec_unc = 0.303575; return; }
517  else if (glu_mass == 1925) { xsec = 0.0014401; xsec_unc = 0.305951; return; }
518  else if (glu_mass == 1930) { xsec = 0.0014016; xsec_unc = 0.305171; return; }
519  else if (glu_mass == 1935) { xsec = 0.00136297; xsec_unc = 0.304873; return; }
520  else if (glu_mass == 1940) { xsec = 0.001331; xsec_unc = 0.307414; return; }
521  else if (glu_mass == 1945) { xsec = 0.001299; xsec_unc = 0.310066; return; }
522  else if (glu_mass == 1950) { xsec = 0.0012642; xsec_unc = 0.304581; return; }
523  else if (glu_mass == 1955) { xsec = 0.00123087; xsec_unc = 0.308644; return; }
524  else if (glu_mass == 1960) { xsec = 0.00120048; xsec_unc = 0.309669; return; }
525  else if (glu_mass == 1965) { xsec = 0.00117053; xsec_unc = 0.310216; return; }
526  else if (glu_mass == 1970) { xsec = 0.00114051; xsec_unc = 0.310814; return; }
527  else if (glu_mass == 1975) { xsec = 0.00111722; xsec_unc = 0.315357; return; }
528  else if (glu_mass == 1980) { xsec = 0.00108758; xsec_unc = 0.315568; return; }
529  else if (glu_mass == 1985) { xsec = 0.00105813; xsec_unc = 0.315103; return; }
530  else if (glu_mass == 1990) { xsec = 0.00102936; xsec_unc = 0.314167; return; }
531  else if (glu_mass == 1995) { xsec = 0.00100614; xsec_unc = 0.317628; return; }
532  else if (glu_mass == 2000) { xsec = 0.000981077; xsec_unc = 0.318422; return; }
533  else if (glu_mass == 2005) { xsec = 0.000956613; xsec_unc = 0.319751; return; }
534  else if (glu_mass == 2010) { xsec = 0.000932325; xsec_unc = 0.320882; return; }
535  else if (glu_mass == 2015) { xsec = 0.000909408; xsec_unc = 0.321466; return; }
536  else if (glu_mass == 2020) { xsec = 0.000886417; xsec_unc = 0.322175; return; }
537  else if (glu_mass == 2025) { xsec = 0.000863806; xsec_unc = 0.323142; return; }
538  else if (glu_mass == 2030) { xsec = 0.00084231; xsec_unc = 0.324635; return; }
539  else if (glu_mass == 2035) { xsec = 0.000821016; xsec_unc = 0.325832; return; }
540  else if (glu_mass == 2040) { xsec = 0.00080083; xsec_unc = 0.326879; return; }
541  else if (glu_mass == 2045) { xsec = 0.000781098; xsec_unc = 0.327197; return; }
542  else if (glu_mass == 2050) { xsec = 0.000761286; xsec_unc = 0.329341; return; }
543  else if (glu_mass == 2055) { xsec = 0.000742602; xsec_unc = 0.329529; return; }
544  else if (glu_mass == 2060) { xsec = 0.000724162; xsec_unc = 0.331309; return; }
545  else if (glu_mass == 2065) { xsec = 0.000706078; xsec_unc = 0.332221; return; }
546  else if (glu_mass == 2070) { xsec = 0.000688773; xsec_unc = 0.332606; return; }
547  else if (glu_mass == 2075) { xsec = 0.000671106; xsec_unc = 0.3347; return; }
548  else if (glu_mass == 2080) { xsec = 0.000654555; xsec_unc = 0.334624; return; }
549  else if (glu_mass == 2085) { xsec = 0.000638154; xsec_unc = 0.336239; return; }
550  else if (glu_mass == 2090) { xsec = 0.000622299; xsec_unc = 0.336944; return; }
551  else if (glu_mass == 2095) { xsec = 0.000606772; xsec_unc = 0.338846; return; }
552  else if (glu_mass == 2100) { xsec = 0.000591918; xsec_unc = 0.339326; return; }
553  else if (glu_mass == 2105) { xsec = 0.000577016; xsec_unc = 0.339891; return; }
554  else if (glu_mass == 2110) { xsec = 0.000562775; xsec_unc = 0.34126; return; }
555  else if (glu_mass == 2115) { xsec = 0.000548803; xsec_unc = 0.341942; return; }
556  else if (glu_mass == 2120) { xsec = 0.000535625; xsec_unc = 0.343143; return; }
557  else if (glu_mass == 2125) { xsec = 0.000522481; xsec_unc = 0.344322; return; }
558  else if (glu_mass == 2130) { xsec = 0.000509317; xsec_unc = 0.345795; return; }
559  else if (glu_mass == 2135) { xsec = 0.000496197; xsec_unc = 0.346013; return; }
560  else if (glu_mass == 2140) { xsec = 0.000484078; xsec_unc = 0.34712; return; }
561  else if (glu_mass == 2145) { xsec = 0.000472404; xsec_unc = 0.346967; return; }
562  else if (glu_mass == 2150) { xsec = 0.000460941; xsec_unc = 0.349082; return; }
563  else if (glu_mass == 2155) { xsec = 0.00044976; xsec_unc = 0.350413; return; }
564  else if (glu_mass == 2160) { xsec = 0.000438717; xsec_unc = 0.351403; return; }
565  else if (glu_mass == 2165) { xsec = 0.000427663; xsec_unc = 0.352429; return; }
566  else if (glu_mass == 2170) { xsec = 0.000417087; xsec_unc = 0.352067; return; }
567  else if (glu_mass == 2175) { xsec = 0.000406541; xsec_unc = 0.354714; return; }
568  else if (glu_mass == 2180) { xsec = 0.000396549; xsec_unc = 0.355491; return; }
569  else if (glu_mass == 2185) { xsec = 0.000387227; xsec_unc = 0.357425; return; }
570  else if (glu_mass == 2190) { xsec = 0.000377685; xsec_unc = 0.356672; return; }
571  else if (glu_mass == 2195) { xsec = 0.00036824; xsec_unc = 0.359098; return; }
572  else if (glu_mass == 2200) { xsec = 0.000359318; xsec_unc = 0.359623; return; }
573  else if (glu_mass == 2205) { xsec = 0.000350402; xsec_unc = 0.360185; return; }
574  else if (glu_mass == 2210) { xsec = 0.000342071; xsec_unc = 0.36245; return; }
575  else if (glu_mass == 2215) { xsec = 0.000333108; xsec_unc = 0.363203; return; }
576  else if (glu_mass == 2220) { xsec = 0.000325306; xsec_unc = 0.363282; return; }
577  else if (glu_mass == 2225) { xsec = 0.000316997; xsec_unc = 0.365413; return; }
578  else if (glu_mass == 2230) { xsec = 0.000309731; xsec_unc = 0.367432; return; }
579  else if (glu_mass == 2235) { xsec = 0.000301864; xsec_unc = 0.367687; return; }
580  else if (glu_mass == 2240) { xsec = 0.000294675; xsec_unc = 0.369412; return; }
581  else if (glu_mass == 2245) { xsec = 0.000287398; xsec_unc = 0.371824; return; }
582  else if (glu_mass == 2250) { xsec = 0.00028065; xsec_unc = 0.371485; return; }
583  else if (glu_mass == 2255) { xsec = 0.000273425; xsec_unc = 0.373577; return; }
584  else if (glu_mass == 2260) { xsec = 0.000266633; xsec_unc = 0.373399; return; }
585  else if (glu_mass == 2265) { xsec = 0.000260448; xsec_unc = 0.375246; return; }
586  else if (glu_mass == 2270) { xsec = 0.000254453; xsec_unc = 0.377047; return; }
587  else if (glu_mass == 2275) { xsec = 0.000248358; xsec_unc = 0.378684; return; }
588  else if (glu_mass == 2280) { xsec = 0.00024216; xsec_unc = 0.38071; return; }
589  else if (glu_mass == 2285) { xsec = 0.000235996; xsec_unc = 0.3826; return; }
590  else if (glu_mass == 2290) { xsec = 0.000230291; xsec_unc = 0.381896; return; }
591  else if (glu_mass == 2295) { xsec = 0.000224084; xsec_unc = 0.384083; return; }
592  else if (glu_mass == 2300) { xsec = 0.000219049; xsec_unc = 0.385249; return; }
593  else if (glu_mass == 2305) { xsec = 0.000214032; xsec_unc = 0.386523; return; }
594  else if (glu_mass == 2310) { xsec = 0.000208248; xsec_unc = 0.386168; return; }
595  else if (glu_mass == 2315) { xsec = 0.000203199; xsec_unc = 0.38742; return; }
596  else if (glu_mass == 2320) { xsec = 0.00019809; xsec_unc = 0.389115; return; }
597  else if (glu_mass == 2325) { xsec = 0.000193715; xsec_unc = 0.392614; return; }
598  else if (glu_mass == 2330) { xsec = 0.000189162; xsec_unc = 0.391278; return; }
599  else if (glu_mass == 2335) { xsec = 0.000184075; xsec_unc = 0.393222; return; }
600  else if (glu_mass == 2340) { xsec = 0.000180106; xsec_unc = 0.393788; return; }
601  else if (glu_mass == 2345) { xsec = 0.000175678; xsec_unc = 0.398124; return; }
602  else if (glu_mass == 2350) { xsec = 0.000171031; xsec_unc = 0.396435; return; }
603  else if (glu_mass == 2355) { xsec = 0.000167051; xsec_unc = 0.39707; return; }
604  else if (glu_mass == 2360) { xsec = 0.000163052; xsec_unc = 0.398228; return; }
605  else if (glu_mass == 2365) { xsec = 0.000159079; xsec_unc = 0.398948; return; }
606  else if (glu_mass == 2370) { xsec = 0.000155058; xsec_unc = 0.400085; return; }
607  else if (glu_mass == 2375) { xsec = 0.000151091; xsec_unc = 0.400903; return; }
608  else if (glu_mass == 2380) { xsec = 0.000147865; xsec_unc = 0.404762; return; }
609  else if (glu_mass == 2385) { xsec = 0.000143836; xsec_unc = 0.406052; return; }
610  else if (glu_mass == 2390) { xsec = 0.00014101; xsec_unc = 0.40553; return; }
611  else if (glu_mass == 2395) { xsec = 0.000137035; xsec_unc = 0.406244; return; }
612  else if (glu_mass == 2400) { xsec = 0.000133965; xsec_unc = 0.407945; return; }
613  else if (glu_mass == 2405) { xsec = 0.000130879; xsec_unc = 0.409847; return; }
614  else if (glu_mass == 2410) { xsec = 0.000127783; xsec_unc = 0.411703; return; }
615  else if (glu_mass == 2415) { xsec = 0.0001241; xsec_unc = 0.409566; return; }
616  else if (glu_mass == 2420) { xsec = 0.000121089; xsec_unc = 0.410638; return; }
617  else if (glu_mass == 2425) { xsec = 0.000118777; xsec_unc = 0.41508; return; }
618  else if (glu_mass == 2430) { xsec = 0.000115804; xsec_unc = 0.415799; return; }
619  else if (glu_mass == 2435) { xsec = 0.00011285; xsec_unc = 0.416479; return; }
620  else if (glu_mass == 2440) { xsec = 0.000109965; xsec_unc = 0.416973; return; }
621  else if (glu_mass == 2445) { xsec = 0.000107077; xsec_unc = 0.416702; return; }
622  else if (glu_mass == 2450) { xsec = 0.000104886; xsec_unc = 0.419997; return; }
623  else if (glu_mass == 2455) { xsec = 0.000102035; xsec_unc = 0.419184; return; }
624  else if (glu_mass == 2460) { xsec = 9.97969e-05; xsec_unc = 0.421891; return; }
625  else if (glu_mass == 2465) { xsec = 9.73046e-05; xsec_unc = 0.422923; return; }
626  else if (glu_mass == 2470) { xsec = 9.50093e-05; xsec_unc = 0.423683; return; }
627  else if (glu_mass == 2475) { xsec = 9.26878e-05; xsec_unc = 0.425345; return; }
628  else if (glu_mass == 2480) { xsec = 9.04143e-05; xsec_unc = 0.425981; return; }
629  else if (glu_mass == 2485) { xsec = 8.82972e-05; xsec_unc = 0.42816; return; }
630  else if (glu_mass == 2490) { xsec = 8.61566e-05; xsec_unc = 0.428228; return; }
631  else if (glu_mass == 2495) { xsec = 8.40135e-05; xsec_unc = 0.429963; return; }
632  else if (glu_mass == 2500) { xsec = 8.20068e-05; xsec_unc = 0.431071; return; }
633  else {xsec = 0.; xsec_unc = 0.;}
634  }
635 
636  void higgsinoCrossSection(int hig_mass, float &xsec, float &xsec_unc) {
637  if(hig_mass == 127) { xsec = .5824*.5824*7602.2/1000; xsec_unc = 0.0393921; return;}
638  else if(hig_mass == 150) { xsec = .5824*.5824*3832.31/1000; xsec_unc = 0.0413612; return;}
639  else if(hig_mass == 175) { xsec = .5824*.5824*2267.94/1000; xsec_unc = 0.044299; return;}
640  else if(hig_mass == 200) { xsec = .5824*.5824*1335.62/1000; xsec_unc = 0.0474362; return;}
641  else if(hig_mass == 225) { xsec = .5824*.5824*860.597/1000; xsec_unc = 0.0504217; return;}
642  else if(hig_mass == 250) { xsec = .5824*.5824*577.314/1000; xsec_unc = 0.0532731; return;}
643  else if(hig_mass == 275) { xsec = .5824*.5824*400.107/1000; xsec_unc = 0.0560232; return;}
644  else if(hig_mass == 300) { xsec = .5824*.5824*284.855/1000; xsec_unc = 0.0586867; return;}
645  else if(hig_mass == 325) { xsec = .5824*.5824* 207.36/1000; xsec_unc = 0.0613554; return;}
646  else if(hig_mass == 350) { xsec = .5824*.5824*153.841/1000; xsec_unc = 0.0640598; return;}
647  else if(hig_mass == 375) { xsec = .5824*.5824*116.006/1000; xsec_unc = 0.066892; return;}
648  else if(hig_mass == 400) { xsec = .5824*.5824*88.7325/1000; xsec_unc = 0.0697517; return;}
649  else if(hig_mass == 425) { xsec = .5824*.5824*68.6963/1000; xsec_unc = 0.0723531; return;}
650  else if(hig_mass == 450) { xsec = .5824*.5824*53.7702/1000; xsec_unc = 0.0748325; return;}
651  else if(hig_mass == 475) { xsec = .5824*.5824*42.4699/1000; xsec_unc = 0.0775146; return;}
652  else if(hig_mass == 500) { xsec = .5824*.5824*33.8387/1000; xsec_unc = 0.0802572; return;}
653  else if(hig_mass == 525) { xsec = .5824*.5824*27.1867/1000; xsec_unc = 0.0825803; return;}
654  else if(hig_mass == 550) { xsec = .5824*.5824*21.9868/1000; xsec_unc = 0.0849278; return;}
655  else if(hig_mass == 575) { xsec = .5824*.5824*17.9062/1000; xsec_unc = 0.087561; return;}
656  else if(hig_mass == 600) { xsec = .5824*.5824*14.6677/1000; xsec_unc = 0.0900693; return;}
657  else if(hig_mass == 625) { xsec = .5824*.5824* 12.062/1000; xsec_unc = 0.091959; return;}
658  else if(hig_mass == 650) { xsec = .5824*.5824*9.96406/1000; xsec_unc = 0.094065; return;}
659  else if(hig_mass == 675) { xsec = .5824*.5824*8.28246/1000; xsec_unc = 0.0957436; return;}
660  else if(hig_mass == 700) { xsec = .5824*.5824*6.89981/1000; xsec_unc = 0.0982894; return;}
661  else if(hig_mass == 725) { xsec = .5824*.5824*5.78355/1000; xsec_unc = 0.0999915; return;}
662  else if(hig_mass == 750) { xsec = .5824*.5824* 4.8731/1000; xsec_unc = 0.101211; return;}
663  else if(hig_mass == 775) { xsec = .5824*.5824*4.09781/1000; xsec_unc = 0.104646; return;}
664  else if(hig_mass == 800) { xsec = .5824*.5824*3.46143/1000; xsec_unc = 0.107618; return;}
665  else if(hig_mass == 825) { xsec = .5824*.5824* 2.9337/1000; xsec_unc = 0.108353; return;}
666  else if(hig_mass == 850) { xsec = .5824*.5824* 2.4923/1000; xsec_unc = 0.110016; return;}
667  else if(hig_mass == 875) { xsec = .5824*.5824*2.13679/1000; xsec_unc = 0.112636; return;}
668  else if(hig_mass == 900) { xsec = .5824*.5824*1.80616/1000; xsec_unc = 0.1134; return;}
669  else if(hig_mass == 925) { xsec = .5824*.5824*1.55453/1000; xsec_unc = 0.116949; return;}
670  else if(hig_mass == 950) { xsec = .5824*.5824*1.32692/1000; xsec_unc = 0.117027; return;}
671  else if(hig_mass == 975) { xsec = .5824*.5824*1.12975/1000; xsec_unc = 0.121244; return;}
672  else if(hig_mass ==1000) { xsec = .5824*.5824*0.968853/1000; xsec_unc = 0.126209; return;}
673  else{ xsec = 0; xsec_unc = 0;}
674  }
675 
676  void stopCrossSection(int stop_mass, double &xsec, double &xsec_unc){
677 
678  if(stop_mass == 100) { xsec = 1521.11; xsec_unc = 0.154038; return; }
679  else if(stop_mass == 105) { xsec = 1233.18; xsec_unc = 0.154059; return; }
680  else if(stop_mass == 110) { xsec = 1013.76; xsec_unc = 0.154088; return; }
681  else if(stop_mass == 115) { xsec = 832.656; xsec_unc = 0.151503; return; }
682  else if(stop_mass == 120) { xsec = 689.799; xsec_unc = 0.15044; return; }
683  else if(stop_mass == 125) { xsec = 574.981; xsec_unc = 0.149895; return; }
684  else if(stop_mass == 130) { xsec = 481.397; xsec_unc = 0.148906; return; }
685  else if(stop_mass == 135) { xsec = 405.159; xsec_unc = 0.148952; return; }
686  else if(stop_mass == 140) { xsec = 342.865; xsec_unc = 0.149119; return; }
687  else if(stop_mass == 145) { xsec = 291.752; xsec_unc = 0.148022; return; }
688  else if(stop_mass == 150) { xsec = 249.409; xsec_unc = 0.147477; return; }
689  else if(stop_mass == 155) { xsec = 214.221; xsec_unc = 0.145928; return; }
690  else if(stop_mass == 160) { xsec = 184.623; xsec_unc = 0.145821; return; }
691  else if(stop_mass == 165) { xsec = 159.614; xsec_unc = 0.147859; return; }
692  else if(stop_mass == 170) { xsec = 139.252; xsec_unc = 0.14547; return; }
693  else if(stop_mass == 175) { xsec = 121.416; xsec_unc = 0.146341; return; }
694  else if(stop_mass == 180) { xsec = 106.194; xsec_unc = 0.142033; return; }
695  else if(stop_mass == 185) { xsec = 93.3347; xsec_unc = 0.144893; return; }
696  else if(stop_mass == 190) { xsec = 82.2541; xsec_unc = 0.144677; return; }
697  else if(stop_mass == 195) { xsec = 72.7397; xsec_unc = 0.144452; return; }
698  else if(stop_mass == 200) { xsec = 64.5085; xsec_unc = 0.144098; return; }
699  else if(stop_mass == 205) { xsec = 57.2279; xsec_unc = 0.144191; return; }
700  else if(stop_mass == 210) { xsec = 50.9226; xsec_unc = 0.142457; return; }
701  else if(stop_mass == 215) { xsec = 45.3761; xsec_unc = 0.14344; return; }
702  else if(stop_mass == 220) { xsec = 40.5941; xsec_unc = 0.142634; return; }
703  else if(stop_mass == 225) { xsec = 36.3818; xsec_unc = 0.142189; return; }
704  else if(stop_mass == 230) { xsec = 32.6679; xsec_unc = 0.141592; return; }
705  else if(stop_mass == 235) { xsec = 29.3155; xsec_unc = 0.142233; return; }
706  else if(stop_mass == 240) { xsec = 26.4761; xsec_unc = 0.141723; return; }
707  else if(stop_mass == 245) { xsec = 23.8853; xsec_unc = 0.139482; return; }
708  else if(stop_mass == 250) { xsec = 21.5949; xsec_unc = 0.140595; return; }
709  else if(stop_mass == 255) { xsec = 19.5614; xsec_unc = 0.138755; return; }
710  else if(stop_mass == 260) { xsec = 17.6836; xsec_unc = 0.139505; return; }
711  else if(stop_mass == 265) { xsec = 16.112; xsec_unc = 0.139531; return; }
712  else if(stop_mass == 270) { xsec = 14.6459; xsec_unc = 0.139278; return; }
713  else if(stop_mass == 275) { xsec = 13.3231; xsec_unc = 0.142549; return; }
714  else if(stop_mass == 280) { xsec = 12.1575; xsec_unc = 0.141584; return; }
715  else if(stop_mass == 285) { xsec = 11.0925; xsec_unc = 0.140904; return; }
716  else if(stop_mass == 290) { xsec = 10.1363; xsec_unc = 0.138967; return; }
717  else if(stop_mass == 295) { xsec = 9.29002; xsec_unc = 0.139107; return; }
718  else if(stop_mass == 300) { xsec = 8.51615; xsec_unc = 0.139223; return; }
719  else if(stop_mass == 305) { xsec = 7.81428; xsec_unc = 0.138996; return; }
720  else if(stop_mass == 310) { xsec = 7.17876; xsec_unc = 0.139357; return; }
721  else if(stop_mass == 315) { xsec = 6.60266; xsec_unc = 0.139256; return; }
722  else if(stop_mass == 320) { xsec = 6.08444; xsec_unc = 0.137957; return; }
723  else if(stop_mass == 325) { xsec = 5.60471; xsec_unc = 0.138144; return; }
724  else if(stop_mass == 330) { xsec = 5.17188; xsec_unc = 0.136954; return; }
725  else if(stop_mass == 335) { xsec = 4.77871; xsec_unc = 0.137554; return; }
726  else if(stop_mass == 340) { xsec = 4.41629; xsec_unc = 0.137945; return; }
727  else if(stop_mass == 345) { xsec = 4.08881; xsec_unc = 0.137075; return; }
728  else if(stop_mass == 350) { xsec = 3.78661; xsec_unc = 0.136877; return; }
729  else if(stop_mass == 355) { xsec = 3.50911; xsec_unc = 0.138089; return; }
730  else if(stop_mass == 360) { xsec = 3.25619; xsec_unc = 0.138002; return; }
731  else if(stop_mass == 365) { xsec = 3.02472; xsec_unc = 0.137093; return; }
732  else if(stop_mass == 370) { xsec = 2.8077; xsec_unc = 0.138064; return; }
733  else if(stop_mass == 375) { xsec = 2.61162; xsec_unc = 0.138477; return; }
734  else if(stop_mass == 380) { xsec = 2.43031; xsec_unc = 0.136999; return; }
735  else if(stop_mass == 385) { xsec = 2.26365; xsec_unc = 0.13728; return; }
736  else if(stop_mass == 390) { xsec = 2.10786; xsec_unc = 0.13732; return; }
737  else if(stop_mass == 395) { xsec = 1.9665; xsec_unc = 0.134737; return; }
738  else if(stop_mass == 400) { xsec = 1.83537; xsec_unc = 0.136985; return; }
739  else if(stop_mass == 405) { xsec = 1.70927; xsec_unc = 0.137114; return; }
740  else if(stop_mass == 410) { xsec = 1.60378; xsec_unc = 0.135468; return; }
741  else if(stop_mass == 415) { xsec = 1.49798; xsec_unc = 0.134453; return; }
742  else if(stop_mass == 420) { xsec = 1.39688; xsec_unc = 0.136719; return; }
743  else if(stop_mass == 425) { xsec = 1.31169; xsec_unc = 0.135013; return; }
744  else if(stop_mass == 430) { xsec = 1.22589; xsec_unc = 0.133237; return; }
745  else if(stop_mass == 435) { xsec = 1.14553; xsec_unc = 0.135478; return; }
746  else if(stop_mass == 440) { xsec = 1.07484; xsec_unc = 0.137238; return; }
747  else if(stop_mass == 445) { xsec = 1.01019; xsec_unc = 0.134187; return; }
748  else if(stop_mass == 450) { xsec = 0.948333; xsec_unc = 0.134559; return; }
749  else if(stop_mass == 455) { xsec = 0.890847; xsec_unc = 0.134587; return; }
750  else if(stop_mass == 460) { xsec = 0.836762; xsec_unc = 0.134468; return; }
751  else if(stop_mass == 465) { xsec = 0.787221; xsec_unc = 0.134149; return; }
752  else if(stop_mass == 470) { xsec = 0.740549; xsec_unc = 0.134127; return; }
753  else if(stop_mass == 475) { xsec = 0.697075; xsec_unc = 0.133926; return; }
754  else if(stop_mass == 480) { xsec = 0.655954; xsec_unc = 0.134392; return; }
755  else if(stop_mass == 485) { xsec = 0.618562; xsec_unc = 0.133705; return; }
756  else if(stop_mass == 490) { xsec = 0.582467; xsec_unc = 0.133914; return; }
757  else if(stop_mass == 495) { xsec = 0.549524; xsec_unc = 0.133691; return; }
758  else if(stop_mass == 500) { xsec = 0.51848; xsec_unc = 0.133797; return; }
759  else if(stop_mass == 505) { xsec = 0.489324; xsec_unc = 0.133608; return; }
760  else if(stop_mass == 510) { xsec = 0.462439; xsec_unc = 0.133046; return; }
761  else if(stop_mass == 515) { xsec = 0.436832; xsec_unc = 0.133703; return; }
762  else if(stop_mass == 520) { xsec = 0.412828; xsec_unc = 0.13272; return; }
763  else if(stop_mass == 525) { xsec = 0.390303; xsec_unc = 0.133443; return; }
764  else if(stop_mass == 530) { xsec = 0.368755; xsec_unc = 0.133769; return; }
765  else if(stop_mass == 535) { xsec = 0.348705; xsec_unc = 0.132706; return; }
766  else if(stop_mass == 540) { xsec = 0.330157; xsec_unc = 0.132981; return; }
767  else if(stop_mass == 545) { xsec = 0.312672; xsec_unc = 0.13277; return; }
768  else if(stop_mass == 550) { xsec = 0.296128; xsec_unc = 0.132687; return; }
769  else if(stop_mass == 555) { xsec = 0.280734; xsec_unc = 0.132363; return; }
770  else if(stop_mass == 560) { xsec = 0.266138; xsec_unc = 0.13193; return; }
771  else if(stop_mass == 565) { xsec = 0.251557; xsec_unc = 0.131731; return; }
772  else if(stop_mass == 570) { xsec = 0.238537; xsec_unc = 0.133409; return; }
773  else if(stop_mass == 575) { xsec = 0.226118; xsec_unc = 0.132741; return; }
774  else if(stop_mass == 580) { xsec = 0.214557; xsec_unc = 0.131697; return; }
775  else if(stop_mass == 585) { xsec = 0.203566; xsec_unc = 0.133257; return; }
776  else if(stop_mass == 590) { xsec = 0.193079; xsec_unc = 0.132037; return; }
777  else if(stop_mass == 595) { xsec = 0.183604; xsec_unc = 0.130973; return; }
778  else if(stop_mass == 600) { xsec = 0.174599; xsec_unc = 0.132074; return; }
779  else if(stop_mass == 605) { xsec = 0.166131; xsec_unc = 0.130154; return; }
780  else if(stop_mass == 610) { xsec = 0.158242; xsec_unc = 0.13142; return; }
781  else if(stop_mass == 615) { xsec = 0.150275; xsec_unc = 0.13285; return; }
782  else if(stop_mass == 620) { xsec = 0.142787; xsec_unc = 0.130642; return; }
783  else if(stop_mass == 625) { xsec = 0.136372; xsec_unc = 0.127962; return; }
784  else if(stop_mass == 630) { xsec = 0.129886; xsec_unc = 0.132957; return; }
785  else if(stop_mass == 635) { xsec = 0.123402; xsec_unc = 0.13016; return; }
786  else if(stop_mass == 640) { xsec = 0.11795; xsec_unc = 0.127132; return; }
787  else if(stop_mass == 645) { xsec = 0.112008; xsec_unc = 0.12808; return; }
788  else if(stop_mass == 650) { xsec = 0.107045; xsec_unc = 0.129232; return; }
789  else if(stop_mass == 655) { xsec = 0.102081; xsec_unc = 0.130012; return; }
790  else if(stop_mass == 660) { xsec = 0.09725; xsec_unc = 0.129038; return; }
791  else if(stop_mass == 665) { xsec = 0.0927515; xsec_unc = 0.129548; return; }
792  else if(stop_mass == 670) { xsec = 0.0885084; xsec_unc = 0.130218; return; }
793  else if(stop_mass == 675) { xsec = 0.0844877; xsec_unc = 0.130703; return; }
794  else if(stop_mass == 680) { xsec = 0.0806192; xsec_unc = 0.131131; return; }
795  else if(stop_mass == 685) { xsec = 0.0769099; xsec_unc = 0.131517; return; }
796  else if(stop_mass == 690) { xsec = 0.0734901; xsec_unc = 0.132344; return; }
797  else if(stop_mass == 695) { xsec = 0.0701805; xsec_unc = 0.132716; return; }
798  else if(stop_mass == 700) { xsec = 0.0670476; xsec_unc = 0.133429; return; }
799  else if(stop_mass == 705) { xsec = 0.0641426; xsec_unc = 0.13363; return; }
800  else if(stop_mass == 710) { xsec = 0.0612942; xsec_unc = 0.133941; return; }
801  else if(stop_mass == 715) { xsec = 0.0585678; xsec_unc = 0.134663; return; }
802  else if(stop_mass == 720) { xsec = 0.0560753; xsec_unc = 0.134984; return; }
803  else if(stop_mass == 725) { xsec = 0.0536438; xsec_unc = 0.135804; return; }
804  else if(stop_mass == 730) { xsec = 0.0513219; xsec_unc = 0.135682; return; }
805  else if(stop_mass == 735) { xsec = 0.0491001; xsec_unc = 0.136268; return; }
806  else if(stop_mass == 740) { xsec = 0.0470801; xsec_unc = 0.136895; return; }
807  else if(stop_mass == 745) { xsec = 0.045061; xsec_unc = 0.136816; return; }
808  else if(stop_mass == 750) { xsec = 0.0431418; xsec_unc = 0.137455; return; }
809  else if(stop_mass == 755) { xsec = 0.0413447; xsec_unc = 0.137833; return; }
810  else if(stop_mass == 760) { xsec = 0.0396264; xsec_unc = 0.138518; return; }
811  else if(stop_mass == 765) { xsec = 0.0379036; xsec_unc = 0.138537; return; }
812  else if(stop_mass == 770) { xsec = 0.0363856; xsec_unc = 0.139334; return; }
813  else if(stop_mass == 775) { xsec = 0.0348796; xsec_unc = 0.139597; return; }
814  else if(stop_mass == 780) { xsec = 0.0334669; xsec_unc = 0.140267; return; }
815  else if(stop_mass == 785) { xsec = 0.0320548; xsec_unc = 0.140406; return; }
816  else if(stop_mass == 790) { xsec = 0.0307373; xsec_unc = 0.14115; return; }
817  else if(stop_mass == 795) { xsec = 0.0295348; xsec_unc = 0.141397; return; }
818  else if(stop_mass == 800) { xsec = 0.0283338; xsec_unc = 0.14171; return; }
819  else if(stop_mass == 805) { xsec = 0.0272206; xsec_unc = 0.14241; return; }
820  else if(stop_mass == 810) { xsec = 0.0261233; xsec_unc = 0.142891; return; }
821  else if(stop_mass == 815) { xsec = 0.0251107; xsec_unc = 0.143632; return; }
822  else if(stop_mass == 820) { xsec = 0.0241099; xsec_unc = 0.143805; return; }
823  else if(stop_mass == 825) { xsec = 0.0230866; xsec_unc = 0.144428; return; }
824  else if(stop_mass == 830) { xsec = 0.0221834; xsec_unc = 0.144791; return; }
825  else if(stop_mass == 835) { xsec = 0.0213766; xsec_unc = 0.145511; return; }
826  else if(stop_mass == 840) { xsec = 0.0204715; xsec_unc = 0.146131; return; }
827  else if(stop_mass == 845) { xsec = 0.0197653; xsec_unc = 0.146602; return; }
828  else if(stop_mass == 850) { xsec = 0.0189612; xsec_unc = 0.14702; return; }
829  else if(stop_mass == 855) { xsec = 0.0182516; xsec_unc = 0.147648; return; }
830  else if(stop_mass == 860) { xsec = 0.0175509; xsec_unc = 0.147944; return; }
831  else if(stop_mass == 865) { xsec = 0.0168336; xsec_unc = 0.148528; return; }
832  else if(stop_mass == 870) { xsec = 0.0162314; xsec_unc = 0.148772; return; }
833  else if(stop_mass == 875) { xsec = 0.015625; xsec_unc = 0.149567; return; }
834  else if(stop_mass == 880) { xsec = 0.0150143; xsec_unc = 0.150389; return; }
835  else if(stop_mass == 885) { xsec = 0.0144112; xsec_unc = 0.150614; return; }
836  else if(stop_mass == 890) { xsec = 0.0138979; xsec_unc = 0.151; return; }
837  else if(stop_mass == 895) { xsec = 0.0133962; xsec_unc = 0.151325; return; }
838  else if(stop_mass == 900) { xsec = 0.0128895; xsec_unc = 0.152026; return; }
839  else if(stop_mass == 905) { xsec = 0.0123843; xsec_unc = 0.152968; return; }
840  else if(stop_mass == 910) { xsec = 0.0119837; xsec_unc = 0.153089; return; }
841  else if(stop_mass == 915) { xsec = 0.0114713; xsec_unc = 0.153678; return; }
842  else if(stop_mass == 920) { xsec = 0.0110688; xsec_unc = 0.154082; return; }
843  else if(stop_mass == 925) { xsec = 0.0106631; xsec_unc = 0.154806; return; }
844  else if(stop_mass == 930) { xsec = 0.0102629; xsec_unc = 0.155313; return; }
845  else if(stop_mass == 935) { xsec = 0.0098874; xsec_unc = 0.156066; return; }
846  else if(stop_mass == 940) { xsec = 0.00952142; xsec_unc = 0.156055; return; }
847  else if(stop_mass == 945) { xsec = 0.00916636; xsec_unc = 0.156849; return; }
848  else if(stop_mass == 950) { xsec = 0.00883465; xsec_unc = 0.157177; return; }
849  else if(stop_mass == 955) { xsec = 0.00851073; xsec_unc = 0.158094; return; }
850  else if(stop_mass == 960) { xsec = 0.00820884; xsec_unc = 0.15844; return; }
851  else if(stop_mass == 965) { xsec = 0.00791403; xsec_unc = 0.159216; return; }
852  else if(stop_mass == 970) { xsec = 0.00763112; xsec_unc = 0.159742; return; }
853  else if(stop_mass == 975) { xsec = 0.00735655; xsec_unc = 0.160548; return; }
854  else if(stop_mass == 980) { xsec = 0.00710317; xsec_unc = 0.160626; return; }
855  else if(stop_mass == 985) { xsec = 0.00684867; xsec_unc = 0.16144; return; }
856  else if(stop_mass == 990) { xsec = 0.00660695; xsec_unc = 0.161813; return; }
857  else if(stop_mass == 995) { xsec = 0.00637546; xsec_unc = 0.162158; return; }
858  else if(stop_mass == 1000) { xsec = 0.00615134; xsec_unc = 0.162953; return; }
859  else if(stop_mass == 1005) { xsec = 0.00593765; xsec_unc = 0.163716; return; }
860  else if(stop_mass == 1010) { xsec = 0.00572452; xsec_unc = 0.163857; return; }
861  else if(stop_mass == 1015) { xsec = 0.00553094; xsec_unc = 0.164628; return; }
862  else if(stop_mass == 1020) { xsec = 0.00533968; xsec_unc = 0.164963; return; }
863  else if(stop_mass == 1025) { xsec = 0.00514619; xsec_unc = 0.165762; return; }
864  else if(stop_mass == 1030) { xsec = 0.00497235; xsec_unc = 0.165838; return; }
865  else if(stop_mass == 1035) { xsec = 0.00479906; xsec_unc = 0.166646; return; }
866  else if(stop_mass == 1040) { xsec = 0.00463806; xsec_unc = 0.166947; return; }
867  else if(stop_mass == 1045) { xsec = 0.00447537; xsec_unc = 0.167071; return; }
868  else if(stop_mass == 1050) { xsec = 0.00432261; xsec_unc = 0.167859; return; }
869  else if(stop_mass == 1055) { xsec = 0.00417983; xsec_unc = 0.168637; return; }
870  else if(stop_mass == 1060) { xsec = 0.00403886; xsec_unc = 0.168981; return; }
871  else if(stop_mass == 1065) { xsec = 0.0038962; xsec_unc = 0.169794; return; }
872  else if(stop_mass == 1070) { xsec = 0.00376343; xsec_unc = 0.169764; return; }
873  else if(stop_mass == 1075) { xsec = 0.00364174; xsec_unc = 0.170634; return; }
874  else if(stop_mass == 1080) { xsec = 0.00352093; xsec_unc = 0.170908; return; }
875  else if(stop_mass == 1085) { xsec = 0.00339813; xsec_unc = 0.171929; return; }
876  else if(stop_mass == 1090) { xsec = 0.00328695; xsec_unc = 0.172274; return; }
877  else if(stop_mass == 1095) { xsec = 0.00317628; xsec_unc = 0.172617; return; }
878  else if(stop_mass == 1100) { xsec = 0.00307413; xsec_unc = 0.173377; return; }
879  else if(stop_mass == 1105) { xsec = 0.00297377; xsec_unc = 0.173822; return; }
880  else if(stop_mass == 1110) { xsec = 0.00287148; xsec_unc = 0.174725; return; }
881  else if(stop_mass == 1115) { xsec = 0.00278078; xsec_unc = 0.175091; return; }
882  else if(stop_mass == 1120) { xsec = 0.00268873; xsec_unc = 0.175883; return; }
883  else if(stop_mass == 1125) { xsec = 0.00260821; xsec_unc = 0.176126; return; }
884  else if(stop_mass == 1130) { xsec = 0.00251529; xsec_unc = 0.176836; return; }
885  else if(stop_mass == 1135) { xsec = 0.00243484; xsec_unc = 0.177128; return; }
886  else if(stop_mass == 1140) { xsec = 0.00236295; xsec_unc = 0.177977; return; }
887  else if(stop_mass == 1145) { xsec = 0.00228192; xsec_unc = 0.178507; return; }
888  else if(stop_mass == 1150) { xsec = 0.00221047; xsec_unc = 0.179259; return; }
889  else if(stop_mass == 1155) { xsec = 0.00213907; xsec_unc = 0.180255; return; }
890  else if(stop_mass == 1160) { xsec = 0.00206845; xsec_unc = 0.180518; return; }
891  else if(stop_mass == 1165) { xsec = 0.0020063; xsec_unc = 0.180954; return; }
892  else if(stop_mass == 1170) { xsec = 0.00194569; xsec_unc = 0.181194; return; }
893  else if(stop_mass == 1175) { xsec = 0.0018741; xsec_unc = 0.182145; return; }
894  else if(stop_mass == 1180) { xsec = 0.00182266; xsec_unc = 0.183074; return; }
895  else if(stop_mass == 1185) { xsec = 0.00176211; xsec_unc = 0.183375; return; }
896  else if(stop_mass == 1190) { xsec = 0.00170006; xsec_unc = 0.184075; return; }
897  else if(stop_mass == 1195) { xsec = 0.00164968; xsec_unc = 0.184438; return; }
898  else if(stop_mass == 1200) { xsec = 0.00159844; xsec_unc = 0.185209; return; }
899  else if(stop_mass == 1205) { xsec = 0.0015472; xsec_unc = 0.185977; return; }
900  else if(stop_mass == 1210) { xsec = 0.00149657; xsec_unc = 0.186485; return; }
901  else if(stop_mass == 1215) { xsec = 0.00145544; xsec_unc = 0.187347; return; }
902  else if(stop_mass == 1220) { xsec = 0.00140288; xsec_unc = 0.188774; return; }
903  else if(stop_mass == 1225) { xsec = 0.00136155; xsec_unc = 0.18989; return; }
904  else if(stop_mass == 1230) { xsec = 0.00131271; xsec_unc = 0.188763; return; }
905  else if(stop_mass == 1235) { xsec = 0.0012717; xsec_unc = 0.189588; return; }
906  else if(stop_mass == 1240) { xsec = 0.00123066; xsec_unc = 0.19049; return; }
907  else if(stop_mass == 1245) { xsec = 0.00119994; xsec_unc = 0.191442; return; }
908  else if(stop_mass == 1250) { xsec = 0.0011583; xsec_unc = 0.193006; return; }
909  else if(stop_mass == 1255) { xsec = 0.00112694; xsec_unc = 0.194441; return; }
910  else if(stop_mass == 1260) { xsec = 0.00108716; xsec_unc = 0.194141; return; }
911  else if(stop_mass == 1265) { xsec = 0.00105517; xsec_unc = 0.196361; return; }
912  else if(stop_mass == 1270) { xsec = 0.00102241; xsec_unc = 0.196297; return; }
913  else if(stop_mass == 1275) { xsec = 0.000991293; xsec_unc = 0.19762; return; }
914  else if(stop_mass == 1280) { xsec = 0.000961012; xsec_unc = 0.197926; return; }
915  else if(stop_mass == 1285) { xsec = 0.000932394; xsec_unc = 0.198682; return; }
916  else if(stop_mass == 1290) { xsec = 0.000903404; xsec_unc = 0.199924; return; }
917  else if(stop_mass == 1295) { xsec = 0.000876957; xsec_unc = 0.200777; return; }
918  else if(stop_mass == 1300) { xsec = 0.000850345; xsec_unc = 0.201604; return; }
919  else if(stop_mass == 1305) { xsec = 0.00082443; xsec_unc = 0.202883; return; }
920  else if(stop_mass == 1310) { xsec = 0.00079983; xsec_unc = 0.20373; return; }
921  else if(stop_mass == 1315) { xsec = 0.000775222; xsec_unc = 0.204622; return; }
922  else if(stop_mass == 1320) { xsec = 0.000751372; xsec_unc = 0.205919; return; }
923  else if(stop_mass == 1325) { xsec = 0.000728912; xsec_unc = 0.206884; return; }
924  else if(stop_mass == 1330) { xsec = 0.000706867; xsec_unc = 0.207763; return; }
925  else if(stop_mass == 1335) { xsec = 0.000685372; xsec_unc = 0.208587; return; }
926  else if(stop_mass == 1340) { xsec = 0.000664649; xsec_unc = 0.209879; return; }
927  else if(stop_mass == 1345) { xsec = 0.000644804; xsec_unc = 0.211487; return; }
928  else if(stop_mass == 1350) { xsec = 0.000625155; xsec_unc = 0.212761; return; }
929  else if(stop_mass == 1355) { xsec = 0.000606802; xsec_unc = 0.213529; return; }
930  else if(stop_mass == 1360) { xsec = 0.000588512; xsec_unc = 0.214428; return; }
931  else if(stop_mass == 1365) { xsec = 0.000570506; xsec_unc = 0.216584; return; }
932  else if(stop_mass == 1370) { xsec = 0.000553379; xsec_unc = 0.216036; return; }
933  else if(stop_mass == 1375) { xsec = 0.000536646; xsec_unc = 0.21775; return; }
934  else if(stop_mass == 1380) { xsec = 0.000521404; xsec_unc = 0.218383; return; }
935  else if(stop_mass == 1385) { xsec = 0.000505008; xsec_unc = 0.219675; return; }
936  else if(stop_mass == 1390) { xsec = 0.000490353; xsec_unc = 0.221444; return; }
937  else if(stop_mass == 1395) { xsec = 0.000476164; xsec_unc = 0.222016; return; }
938  else if(stop_mass == 1400) { xsec = 0.000461944; xsec_unc = 0.222704; return; }
939  else if(stop_mass == 1405) { xsec = 0.000448172; xsec_unc = 0.224911; return; }
940  else if(stop_mass == 1410) { xsec = 0.000435082; xsec_unc = 0.225606; return; }
941  else if(stop_mass == 1415) { xsec = 0.000422967; xsec_unc = 0.226095; return; }
942  else if(stop_mass == 1420) { xsec = 0.000410381; xsec_unc = 0.22797; return; }
943  else if(stop_mass == 1425) { xsec = 0.000398106; xsec_unc = 0.228949; return; }
944  else if(stop_mass == 1430) { xsec = 0.000386792; xsec_unc = 0.231319; return; }
945  else if(stop_mass == 1435) { xsec = 0.000375724; xsec_unc = 0.231724; return; }
946  else if(stop_mass == 1440) { xsec = 0.000364616; xsec_unc = 0.232234; return; }
947  else if(stop_mass == 1445) { xsec = 0.000353965; xsec_unc = 0.234637; return; }
948  else if(stop_mass == 1450) { xsec = 0.000343923; xsec_unc = 0.234948; return; }
949  else if(stop_mass == 1455) { xsec = 0.000333885; xsec_unc = 0.235468; return; }
950  else if(stop_mass == 1460) { xsec = 0.000324344; xsec_unc = 0.23771; return; }
951  else if(stop_mass == 1465) { xsec = 0.0003153; xsec_unc = 0.238004; return; }
952  else if(stop_mass == 1470) { xsec = 0.00030583; xsec_unc = 0.240064; return; }
953  else if(stop_mass == 1475) { xsec = 0.000296811; xsec_unc = 0.240314; return; }
954  else if(stop_mass == 1480) { xsec = 0.000288149; xsec_unc = 0.239248; return; }
955  else if(stop_mass == 1485) { xsec = 0.000279711; xsec_unc = 0.241257; return; }
956  else if(stop_mass == 1490) { xsec = 0.000271724; xsec_unc = 0.241274; return; }
957  else if(stop_mass == 1495) { xsec = 0.000264275; xsec_unc = 0.243545; return; }
958  else if(stop_mass == 1500) { xsec = 0.000256248; xsec_unc = 0.24372; return; }
959  else if(stop_mass == 1505) { xsec = 0.000248853; xsec_unc = 0.245827; return; }
960  else if(stop_mass == 1510) { xsec = 0.000241844; xsec_unc = 0.246187; return; }
961  else if(stop_mass == 1515) { xsec = 0.000234438; xsec_unc = 0.248442; return; }
962  else if(stop_mass == 1520) { xsec = 0.000227374; xsec_unc = 0.248909; return; }
963  else if(stop_mass == 1525) { xsec = 0.000221045; xsec_unc = 0.250895; return; }
964  else if(stop_mass == 1530) { xsec = 0.000214431; xsec_unc = 0.248728; return; }
965  else if(stop_mass == 1535) { xsec = 0.000208092; xsec_unc = 0.251043; return; }
966  else if(stop_mass == 1540) { xsec = 0.000201748; xsec_unc = 0.253207; return; }
967  else if(stop_mass == 1545) { xsec = 0.000196399; xsec_unc = 0.255641; return; }
968  else if(stop_mass == 1550) { xsec = 0.000190474; xsec_unc = 0.255213; return; }
969  else if(stop_mass == 1555) { xsec = 0.000185188; xsec_unc = 0.257329; return; }
970  else if(stop_mass == 1560) { xsec = 0.000179263; xsec_unc = 0.256931; return; }
971  else if(stop_mass == 1565) { xsec = 0.000174021; xsec_unc = 0.259111; return; }
972  else if(stop_mass == 1570) { xsec = 0.000169176; xsec_unc = 0.258106; return; }
973  else if(stop_mass == 1575) { xsec = 0.000163861; xsec_unc = 0.260597; return; }
974  else if(stop_mass == 1580) { xsec = 0.000159583; xsec_unc = 0.262958; return; }
975  else if(stop_mass == 1585) { xsec = 0.000154719; xsec_unc = 0.26195; return; }
976  else if(stop_mass == 1590) { xsec = 0.000150506; xsec_unc = 0.264111; return; }
977  else if(stop_mass == 1595) { xsec = 0.000145626; xsec_unc = 0.263077; return; }
978  else if(stop_mass == 1600) { xsec = 0.000141382; xsec_unc = 0.265291; return; }
979  else if(stop_mass == 1605) { xsec = 0.000137131; xsec_unc = 0.267424; return; }
980  else if(stop_mass == 1610) { xsec = 0.000132187; xsec_unc = 0.26668; return; }
981  else if(stop_mass == 1615) { xsec = 0.000127929; xsec_unc = 0.269117; return; }
982  else if(stop_mass == 1620) { xsec = 0.000124086; xsec_unc = 0.267738; return; }
983  else if(stop_mass == 1625) { xsec = 0.00011982; xsec_unc = 0.270483; return; }
984  else if(stop_mass == 1630) { xsec = 0.000116042; xsec_unc = 0.268071; return; }
985  else if(stop_mass == 1635) { xsec = 0.000112767; xsec_unc = 0.27127; return; }
986  else if(stop_mass == 1640) { xsec = 0.000108936; xsec_unc = 0.269351; return; }
987  else if(stop_mass == 1645) { xsec = 0.000105746; xsec_unc = 0.271783; return; }
988  else if(stop_mass == 1650) { xsec = 0.000102693; xsec_unc = 0.27292; return; }
989  else if(stop_mass == 1655) { xsec = 0.000100112; xsec_unc = 0.274445; return; }
990  else if(stop_mass == 1660) { xsec = 9.75763e-05; xsec_unc = 0.275431; return; }
991  else if(stop_mass == 1665) { xsec = 9.52062e-05; xsec_unc = 0.276946; return; }
992  else if(stop_mass == 1670) { xsec = 9.29857e-05; xsec_unc = 0.277869; return; }
993  else if(stop_mass == 1675) { xsec = 9.08285e-05; xsec_unc = 0.279347; return; }
994  else if(stop_mass == 1680) { xsec = 8.87433e-05; xsec_unc = 0.281539; return; }
995  else if(stop_mass == 1685) { xsec = 8.66618e-05; xsec_unc = 0.283509; return; }
996  else if(stop_mass == 1690) { xsec = 8.46535e-05; xsec_unc = 0.284432; return; }
997  else if(stop_mass == 1695) { xsec = 8.27102e-05; xsec_unc = 0.28591; return; }
998  else if(stop_mass == 1700) { xsec = 8.07774e-05; xsec_unc = 0.287497; return; }
999  else if(stop_mass == 1705) { xsec = 7.8666e-05; xsec_unc = 0.288194; return; }
1000  else if(stop_mass == 1710) { xsec = 7.6572e-05; xsec_unc = 0.290265; return; }
1001  else if(stop_mass == 1715) { xsec = 7.45994e-05; xsec_unc = 0.291193; return; }
1002  else if(stop_mass == 1720) { xsec = 7.25199e-05; xsec_unc = 0.293013; return; }
1003  else if(stop_mass == 1725) { xsec = 7.05189e-05; xsec_unc = 0.293697; return; }
1004  else if(stop_mass == 1730) { xsec = 6.85712e-05; xsec_unc = 0.294972; return; }
1005  else if(stop_mass == 1735) { xsec = 6.67296e-05; xsec_unc = 0.296167; return; }
1006  else if(stop_mass == 1740) { xsec = 6.49184e-05; xsec_unc = 0.297686; return; }
1007  else if(stop_mass == 1745) { xsec = 6.30949e-05; xsec_unc = 0.298524; return; }
1008  else if(stop_mass == 1750) { xsec = 6.13637e-05; xsec_unc = 0.299789; return; }
1009  else if(stop_mass == 1755) { xsec = 5.97301e-05; xsec_unc = 0.300928; return; }
1010  else if(stop_mass == 1760) { xsec = 5.80751e-05; xsec_unc = 0.302585; return; }
1011  else if(stop_mass == 1765) { xsec = 5.65479e-05; xsec_unc = 0.30366; return; }
1012  else if(stop_mass == 1770) { xsec = 5.49998e-05; xsec_unc = 0.305241; return; }
1013  else if(stop_mass == 1775) { xsec = 5.35686e-05; xsec_unc = 0.306718; return; }
1014  else if(stop_mass == 1780) { xsec = 5.20828e-05; xsec_unc = 0.306799; return; }
1015  else if(stop_mass == 1785) { xsec = 5.07079e-05; xsec_unc = 0.309201; return; }
1016  else if(stop_mass == 1790) { xsec = 4.93948e-05; xsec_unc = 0.310043; return; }
1017  else if(stop_mass == 1795) { xsec = 4.80635e-05; xsec_unc = 0.31138; return; }
1018  else if(stop_mass == 1800) { xsec = 4.67492e-05; xsec_unc = 0.312291; return; }
1019  else if(stop_mass == 1805) { xsec = 4.55055e-05; xsec_unc = 0.314321; return; }
1020  else if(stop_mass == 1810) { xsec = 4.42835e-05; xsec_unc = 0.315499; return; }
1021  else if(stop_mass == 1815) { xsec = 4.30744e-05; xsec_unc = 0.316302; return; }
1022  else if(stop_mass == 1820) { xsec = 4.19954e-05; xsec_unc = 0.317151; return; }
1023  else if(stop_mass == 1825) { xsec = 4.08527e-05; xsec_unc = 0.319048; return; }
1024  else if(stop_mass == 1830) { xsec = 3.97561e-05; xsec_unc = 0.319718; return; }
1025  else if(stop_mass == 1835) { xsec = 3.87041e-05; xsec_unc = 0.322028; return; }
1026  else if(stop_mass == 1840) { xsec = 3.76008e-05; xsec_unc = 0.32268; return; }
1027  else if(stop_mass == 1845) { xsec = 3.66914e-05; xsec_unc = 0.324529; return; }
1028  else if(stop_mass == 1850) { xsec = 3.56995e-05; xsec_unc = 0.325039; return; }
1029  else if(stop_mass == 1855) { xsec = 3.47689e-05; xsec_unc = 0.326767; return; }
1030  else if(stop_mass == 1860) { xsec = 3.38528e-05; xsec_unc = 0.328878; return; }
1031  else if(stop_mass == 1865) { xsec = 3.29644e-05; xsec_unc = 0.328975; return; }
1032  else if(stop_mass == 1870) { xsec = 3.20679e-05; xsec_unc = 0.329608; return; }
1033  else if(stop_mass == 1875) { xsec = 3.12583e-05; xsec_unc = 0.331541; return; }
1034  else if(stop_mass == 1880) { xsec = 3.04342e-05; xsec_unc = 0.333117; return; }
1035  else if(stop_mass == 1885) { xsec = 2.96516e-05; xsec_unc = 0.332866; return; }
1036  else if(stop_mass == 1890) { xsec = 2.88952e-05; xsec_unc = 0.336279; return; }
1037  else if(stop_mass == 1895) { xsec = 2.81145e-05; xsec_unc = 0.336845; return; }
1038  else if(stop_mass == 1900) { xsec = 2.73974e-05; xsec_unc = 0.338247; return; }
1039  else if(stop_mass == 1905) { xsec = 2.66796e-05; xsec_unc = 0.339708; return; }
1040  else if(stop_mass == 1910) { xsec = 2.59941e-05; xsec_unc = 0.339526; return; }
1041  else if(stop_mass == 1915) { xsec = 2.52784e-05; xsec_unc = 0.341137; return; }
1042  else if(stop_mass == 1920) { xsec = 2.46598e-05; xsec_unc = 0.342714; return; }
1043  else if(stop_mass == 1925) { xsec = 2.39932e-05; xsec_unc = 0.342328; return; }
1044  else if(stop_mass == 1930) { xsec = 2.33737e-05; xsec_unc = 0.34394; return; }
1045  else if(stop_mass == 1935) { xsec = 2.27623e-05; xsec_unc = 0.345138; return; }
1046  else if(stop_mass == 1940) { xsec = 2.21454e-05; xsec_unc = 0.346933; return; }
1047  else if(stop_mass == 1945) { xsec = 2.15924e-05; xsec_unc = 0.350815; return; }
1048  else if(stop_mass == 1950) { xsec = 2.10232e-05; xsec_unc = 0.349444; return; }
1049  else if(stop_mass == 1955) { xsec = 2.05211e-05; xsec_unc = 0.350155; return; }
1050  else if(stop_mass == 1960) { xsec = 1.98996e-05; xsec_unc = 0.352135; return; }
1051  else if(stop_mass == 1965) { xsec = 1.9408e-05; xsec_unc = 0.353328; return; }
1052  else if(stop_mass == 1970) { xsec = 1.88974e-05; xsec_unc = 0.354643; return; }
1053  else if(stop_mass == 1975) { xsec = 1.84612e-05; xsec_unc = 0.357904; return; }
1054  else if(stop_mass == 1980) { xsec = 1.79562e-05; xsec_unc = 0.358898; return; }
1055  else if(stop_mass == 1985) { xsec = 1.75673e-05; xsec_unc = 0.35989; return; }
1056  else if(stop_mass == 1990) { xsec = 1.70612e-05; xsec_unc = 0.360953; return; }
1057  else if(stop_mass == 1995) { xsec = 1.66228e-05; xsec_unc = 0.364709; return; }
1058  else if(stop_mass == 2000) { xsec = 1.62355e-05; xsec_unc = 0.365277; return; }
1059  else { xsec = 0.; xsec_unc = 0.; }
1060  }
1061 
1062 }
xsec::fractionNegWeights
float fractionNegWeights(const TString &file)
Definition: cross_sections.cc:158
xsec::crossSection
float crossSection(const TString &file)
Definition: cross_sections.cc:10
std
STL namespace.
xsec::higgsinoCrossSection
void higgsinoCrossSection(int hig_mass, float &xsec, float &xsec_unc)
Definition: cross_sections.cc:178
make_dataset_table.xsec
int xsec
Definition: make_dataset_table.py:54
xsec::stopCrossSection
void stopCrossSection(int stop_mass, double &xsec, double &xsec_unc)
Definition: cross_sections.cc:676
cross_sections.hh
xsec::signalCrossSection
void signalCrossSection(int glu_mass, float &xsec, float &xsec_unc)
Definition: cross_sections.cc:218
Generated by   doxygen 1.8.11