#include "thdm.h"

#include <stdlib.h>
#include <iostream>

#include <TFile.h>
#include <TNtuple.h>
#include <TH2.h>
#include <TProfile.h>
#include <TCanvas.h>
#include <TFrame.h>
#include <TROOT.h>
#include <TSystem.h>
#include <TRandom.h>
#include <TBenchmark.h>
#include <TCint.h>
#include <RooRandom.h>
#include <TGraph.h>
#include <TGraphErrors.h>
#include <TGraphAsymmErrors.h>
#include <TMultiGraph.h>
#include <TLegend.h>
#include <TPaveText.h>
#include <TString.h>
#include <TStyle.h>
#include <TGraph2D.h>
#include <TH2.h>
#include <TH2Poly.h>
#include "TColor.h"
#include "TExec.h"
#include "TPaletteAxis.h"
#include "TGaxis.h"

using namespace std;

thdm::thdm()
{
  tree = init("/hepstore/gwilliam/thdm_grid_v163_13TeV.root");
  //tree = init("/hepstore/gwilliam/thdm_grid_v165_AZHHZA-mCHmin.root");
  //tree = init("/hepstore/gwilliam/thdm_grid_v165_AZHHZA-mCHmax.root");
  //tree = init("/hepstore/gwilliam/thdm_grid_v164.root");
  //tree = init("/hepstore/gwilliam/thdm_grid_v164_lowmH.root");
}

thdm::thdm(const char* name) {
  tree = init(name);
}

thdm::~thdm() {
  delete label;
}

double thdm::cosba() {
  double cosba = cba*sba/fabs(sba);
  return cosba;
}

bool thdm::next() {
  long n = tree->GetEntries();
  if (m_ievent < n) {
    loadbar(m_ievent+1, n);
    tree->GetEntry(m_ievent++);    
    return m_ievent;
  }

  std::cout << std::endl;
  return 0;
}

TTree* thdm::init(const char* name) {

    TFile *f = TFile::Open(name);
    TTree *fChain = (TTree*) f->Get("thdm");
   
    fChain->SetBranchAddress("mCh", &mCh, &b_mCh);						 	      
    fChain->SetBranchAddress("mA", &mA, &b_mA);						 	      
    fChain->SetBranchAddress("mH", &mH, &b_mH);						 	      
    fChain->SetBranchAddress("mh", &mh, &b_mh);						 	      
    fChain->SetBranchAddress("tanb", &tanb, &b_tanb);					  
    fChain->SetBranchAddress("m12_2", &m12_2, &b_m12_2);
    fChain->SetBranchAddress("sba", &sba, &b_sba);					 	      
    fChain->SetBranchAddress("cba", &cba, &b_cba);					 	      
    fChain->SetBranchAddress("xsec_H_gg", &xsec_H_gg, &b_xsec_H_gg);			 	      
    fChain->SetBranchAddress("xsec_H_bb", &xsec_H_bb, &b_xsec_H_bb);			 	      
    fChain->SetBranchAddress("xsec_H_VBF", &xsec_H_VBF, &b_xsec_H_VBF);			 	      
    fChain->SetBranchAddress("xsec_H_WH", &xsec_H_WH, &b_xsec_H_WH);			 	      
    fChain->SetBranchAddress("xsec_H_ZH", &xsec_H_ZH, &b_xsec_H_ZH);
    fChain->SetBranchAddress("xsec_H_ttH", &xsec_H_ttH, &b_xsec_H_ttH);			 	      
    fChain->SetBranchAddress("br_h_ZZ", &br_h_ZZ, &b_br_h_ZZ);				  
    fChain->SetBranchAddress("br_h_WW", &br_h_WW, &b_br_h_WW);				  
    fChain->SetBranchAddress("br_h_gammagamma", &br_h_gammagamma, &b_br_h_gammagamma);
    fChain->SetBranchAddress("br_h_bb", &br_h_bb, &b_br_h_bb);
    fChain->SetBranchAddress("br_H_ZZ", &br_H_ZZ, &b_br_H_ZZ);				  
    fChain->SetBranchAddress("br_H_WW", &br_H_WW, &b_br_H_WW);
    fChain->SetBranchAddress("br_H_tautau", &br_H_tautau, &b_br_H_tautau);
    fChain->SetBranchAddress("br_H_gammagamma", &br_H_gammagamma, &b_br_H_gammagamma);
    fChain->SetBranchAddress("br_H_hh", &br_H_hh, &b_br_H_hh);
    fChain->SetBranchAddress("br_H_tt", &br_H_tt, &b_br_H_tt);
    fChain->SetBranchAddress("br_H_bb", &br_H_bb, &b_br_H_bb);
    fChain->SetBranchAddress("xsec_h_gg", &xsec_h_gg, &b_xsec_h_gg);			 	      
    fChain->SetBranchAddress("xsec_h_VBF", &xsec_h_VBF, &b_xsec_h_VBF);			 	      
    fChain->SetBranchAddress("width_H", &width_H, &b_width_H);			 	      
    fChain->SetBranchAddress("width_A", &width_A, &b_width_A);			 	      
    fChain->SetBranchAddress("type", &type, &b_type);					  
    fChain->SetBranchAddress("status", &status, &b_status);	                               				  
    fChain->SetBranchAddress("xsec_A_gg", &xsec_A_gg, &b_xsec_A_gg);			 	      
    fChain->SetBranchAddress("xsec_A_VBF", &xsec_A_VBF, &b_xsec_A_VBF);			 	      
    fChain->SetBranchAddress("xsec_A_bb", &xsec_A_bb, &b_xsec_A_bb);			 	      
    fChain->SetBranchAddress("br_A_Zh", &br_A_Zh, &b_br_A_Zh);				  
    fChain->SetBranchAddress("br_A_ZH", &br_A_ZH, &b_br_A_ZH);				  
    fChain->SetBranchAddress("br_H_ZA", &br_H_ZA, &b_br_H_ZA);				  
    fChain->SetBranchAddress("br_A_bb", &br_H_ZA, &b_br_A_bb);				  

    fChain->SetBranchStatus("*",0); 

    fChain->SetBranchStatus("mCh", 1);
    fChain->SetBranchStatus("mA", 1);
    fChain->SetBranchStatus("mH", 1);
    fChain->SetBranchStatus("mh", 1);
    fChain->SetBranchStatus("tanb", 1);
    fChain->SetBranchStatus("m12_2", 1);
    fChain->SetBranchStatus("sba", 1);
    fChain->SetBranchStatus("cba", 1);
    fChain->SetBranchStatus("xsec_H_gg", 1);
    fChain->SetBranchStatus("xsec_H_bb", 1);
    fChain->SetBranchStatus("xsec_H_VBF", 1);
    fChain->SetBranchStatus("xsec_H_WH", 1);
    fChain->SetBranchStatus("xsec_H_ZH", 1);
    fChain->SetBranchStatus("xsec_H_ttH", 1);
    fChain->SetBranchStatus("br_h_ZZ", 1);
    fChain->SetBranchStatus("br_h_WW", 1);
    fChain->SetBranchStatus("br_h_gammagamma", 1);
    fChain->SetBranchStatus("br_h_bb", 1);
    fChain->SetBranchStatus("br_H_ZZ", 1);
    fChain->SetBranchStatus("br_H_WW", 1);
    fChain->SetBranchStatus("br_H_tautau", 1);
    fChain->SetBranchStatus("br_H_gammagamma", 1);
    fChain->SetBranchStatus("br_H_hh", 1);
    fChain->SetBranchStatus("br_H_tt", 1);
    fChain->SetBranchStatus("br_H_bb", 1);
    fChain->SetBranchStatus("xsec_h_gg", 1);
    fChain->SetBranchStatus("xsec_h_VBF", 1);
    fChain->SetBranchStatus("width_H", 1);
    fChain->SetBranchStatus("width_A", 1);
    fChain->SetBranchStatus("type", 1);
    fChain->SetBranchStatus("status", 1);
    fChain->SetBranchStatus("xsec_A_gg", 1);
    fChain->SetBranchStatus("xsec_A_VBF", 1);
    fChain->SetBranchStatus("xsec_A_bb", 1);
    fChain->SetBranchStatus("br_A_Zh", 1);
    fChain->SetBranchStatus("br_A_ZH", 1);
    fChain->SetBranchStatus("br_H_ZA", 1);
    fChain->SetBranchStatus("br_A_bb", 1);

    label = new TPaveText(0.10,0.85,0.90,0.99,"ndc");
    label->SetBorderSize(0);
    label->SetFillStyle(0);
    label->SetFillColor(0);   

    m_ievent = 0;

    return fChain;
}


TGraphErrors* thdm::smwidth(){
  TGraphErrors *sm = new TGraphErrors(231); 
  sm->SetLineWidth(3);  sm->SetFillColor(1); sm->SetFillStyle(3005); sm->SetTitle("; m_{H} [GeV]; #Gamma [GeV]");
  int i=0;
  sm->SetPoint(i, 90.0,2.22E-03); sm->SetPointError(i++,0,2.22E-03 *(5.2)/100.);
  sm->SetPoint(i, 95.0,2.35E-03); sm->SetPointError(i++,0,2.35E-03 *(5.1)/100.);
  sm->SetPoint(i,100.0,2.49E-03); sm->SetPointError(i++,0,2.49E-03 *(5.1)/100.);
  sm->SetPoint(i,105.0,2.65E-03); sm->SetPointError(i++,0,2.65E-03 *(5.0)/100.);
  sm->SetPoint(i,110.0,2.85E-03); sm->SetPointError(i++,0,2.85E-03 *(4.9)/100.);
  sm->SetPoint(i,110.5,2.87E-03); sm->SetPointError(i++,0,2.87E-03 *(4.9)/100.);
  sm->SetPoint(i,111.0,2.90E-03); sm->SetPointError(i++,0,2.90E-03 *(4.8)/100.);
  sm->SetPoint(i,111.5,2.92E-03); sm->SetPointError(i++,0,2.92E-03 *(4.8)/100.);
  sm->SetPoint(i,112.0,2.95E-03); sm->SetPointError(i++,0,2.95E-03 *(4.8)/100.);
  sm->SetPoint(i,112.5,2.98E-03); sm->SetPointError(i++,0,2.98E-03 *(4.8)/100.);
  sm->SetPoint(i,113.0,3.00E-03); sm->SetPointError(i++,0,3.00E-03 *(4.8)/100.);
  sm->SetPoint(i,113.5,3.03E-03); sm->SetPointError(i++,0,3.03E-03 *(4.7)/100.);
  sm->SetPoint(i,114.0,3.06E-03); sm->SetPointError(i++,0,3.06E-03 *(4.7)/100.);
  sm->SetPoint(i,114.5,3.09E-03); sm->SetPointError(i++,0,3.09E-03 *(4.7)/100.);
  sm->SetPoint(i,115.0,3.12E-03); sm->SetPointError(i++,0,3.12E-03 *(4.7)/100.);
  sm->SetPoint(i,115.5,3.16E-03); sm->SetPointError(i++,0,3.16E-03 *(4.6)/100.);
  sm->SetPoint(i,116.0,3.19E-03); sm->SetPointError(i++,0,3.19E-03 *(4.6)/100.);
  sm->SetPoint(i,116.5,3.23E-03); sm->SetPointError(i++,0,3.23E-03 *(4.6)/100.);
  sm->SetPoint(i,117.0,3.26E-03); sm->SetPointError(i++,0,3.26E-03 *(4.5)/100.);
  sm->SetPoint(i,117.5,3.30E-03); sm->SetPointError(i++,0,3.30E-03 *(4.5)/100.);
  sm->SetPoint(i,118.0,3.34E-03); sm->SetPointError(i++,0,3.34E-03 *(4.5)/100.);
  sm->SetPoint(i,118.5,3.37E-03); sm->SetPointError(i++,0,3.37E-03 *(4.4)/100.);
  sm->SetPoint(i,119.0,3.41E-03); sm->SetPointError(i++,0,3.41E-03 *(4.4)/100.);
  sm->SetPoint(i,119.5,3.46E-03); sm->SetPointError(i++,0,3.46E-03 *(4.4)/100.);
  sm->SetPoint(i,120.0,3.50E-03); sm->SetPointError(i++,0,3.50E-03 *(4.3)/100.);
  sm->SetPoint(i,120.5,3.55E-03); sm->SetPointError(i++,0,3.55E-03 *(4.3)/100.);
  sm->SetPoint(i,121.0,3.60E-03); sm->SetPointError(i++,0,3.60E-03 *(4.2)/100.);
  sm->SetPoint(i,121.5,3.65E-03); sm->SetPointError(i++,0,3.65E-03 *(4.2)/100.);
  sm->SetPoint(i,122.0,3.70E-03); sm->SetPointError(i++,0,3.70E-03 *(4.2)/100.);
  sm->SetPoint(i,122.5,3.76E-03); sm->SetPointError(i++,0,3.76E-03 *(4.1)/100.);
  sm->SetPoint(i,123.0,3.82E-03); sm->SetPointError(i++,0,3.82E-03 *(4.1)/100.);
  sm->SetPoint(i,123.5,3.88E-03); sm->SetPointError(i++,0,3.88E-03 *(4.0)/100.);
  sm->SetPoint(i,124.0,3.94E-03); sm->SetPointError(i++,0,3.94E-03 *(4.0)/100.);
  sm->SetPoint(i,124.5,4.00E-03); sm->SetPointError(i++,0,4.00E-03 *(4.0)/100.);
  sm->SetPoint(i,125.0,4.07E-03); sm->SetPointError(i++,0,4.07E-03 *(3.9)/100.);
  sm->SetPoint(i,125.5,4.14E-03); sm->SetPointError(i++,0,4.14E-03 *(3.9)/100.);
  sm->SetPoint(i,126.0,4.21E-03); sm->SetPointError(i++,0,4.21E-03 *(3.8)/100.);
  sm->SetPoint(i,126.5,4.29E-03); sm->SetPointError(i++,0,4.29E-03 *(3.8)/100.);
  sm->SetPoint(i,127.0,4.37E-03); sm->SetPointError(i++,0,4.37E-03 *(3.7)/100.);
  sm->SetPoint(i,127.5,4.45E-03); sm->SetPointError(i++,0,4.45E-03 *(3.7)/100.);
  sm->SetPoint(i,128.0,4.53E-03); sm->SetPointError(i++,0,4.53E-03 *(3.6)/100.);
  sm->SetPoint(i,128.5,4.62E-03); sm->SetPointError(i++,0,4.62E-03 *(3.6)/100.);
  sm->SetPoint(i,129.0,4.71E-03); sm->SetPointError(i++,0,4.71E-03 *(3.5)/100.);
  sm->SetPoint(i,129.5,4.81E-03); sm->SetPointError(i++,0,4.81E-03 *(3.5)/100.);
  sm->SetPoint(i,130.0,4.91E-03); sm->SetPointError(i++,0,4.91E-03 *(3.4)/100.);
  sm->SetPoint(i,130.5,5.01E-03); sm->SetPointError(i++,0,5.01E-03 *(3.4)/100.);
  sm->SetPoint(i,131.0,5.12E-03); sm->SetPointError(i++,0,5.12E-03 *(3.3)/100.);
  sm->SetPoint(i,131.5,5.23E-03); sm->SetPointError(i++,0,5.23E-03 *(3.3)/100.);
  sm->SetPoint(i,132.0,5.35E-03); sm->SetPointError(i++,0,5.35E-03 *(3.2)/100.);
  sm->SetPoint(i,132.5,5.47E-03); sm->SetPointError(i++,0,5.47E-03 *(3.2)/100.);
  sm->SetPoint(i,133.0,5.60E-03); sm->SetPointError(i++,0,5.60E-03 *(3.1)/100.);
  sm->SetPoint(i,133.5,5.74E-03); sm->SetPointError(i++,0,5.74E-03 *(3.1)/100.);
  sm->SetPoint(i,134.0,5.88E-03); sm->SetPointError(i++,0,5.88E-03 *(3.0)/100.);
  sm->SetPoint(i,134.5,6.03E-03); sm->SetPointError(i++,0,6.03E-03 *(3.0)/100.);
  sm->SetPoint(i,135.0,6.18E-03); sm->SetPointError(i++,0,6.18E-03 *(2.9)/100.);
  sm->SetPoint(i,135.5,6.34E-03); sm->SetPointError(i++,0,6.34E-03 *(2.8)/100.);
  sm->SetPoint(i,136.0,6.51E-03); sm->SetPointError(i++,0,6.51E-03 *(2.7)/100.);
  sm->SetPoint(i,136.5,6.69E-03); sm->SetPointError(i++,0,6.69E-03 *(2.6)/100.);
  sm->SetPoint(i,137.0,6.87E-03); sm->SetPointError(i++,0,6.87E-03 *(2.5)/100.);
  sm->SetPoint(i,137.5,7.06E-03); sm->SetPointError(i++,0,7.06E-03 *(2.3)/100.);
  sm->SetPoint(i,138.0,7.27E-03); sm->SetPointError(i++,0,7.27E-03 *(2.2)/100.);
  sm->SetPoint(i,138.5,7.48E-03); sm->SetPointError(i++,0,7.48E-03 *(2.1)/100.);
  sm->SetPoint(i,139.0,7.70E-03); sm->SetPointError(i++,0,7.70E-03 *(2.0)/100.);
  sm->SetPoint(i,139.5,7.93E-03); sm->SetPointError(i++,0,7.93E-03 *(1.9)/100.);
  sm->SetPoint(i,140.0,8.18E-03); sm->SetPointError(i++,0,8.18E-03 *(1.7)/100.);
  sm->SetPoint(i,141.0,8.70E-03); sm->SetPointError(i++,0,8.70E-03 *(1.7)/100.);
  sm->SetPoint(i,142.0,9.28E-03); sm->SetPointError(i++,0,9.28E-03 *(1.6)/100.);
  sm->SetPoint(i,143.0,9.93E-03); sm->SetPointError(i++,0,9.93E-03 *(1.5)/100.);
  sm->SetPoint(i,144.0,1.07E-02); sm->SetPointError(i++,0,1.07E-02 *(1.5)/100.);
  sm->SetPoint(i,145.0,1.15E-02); sm->SetPointError(i++,0,1.15E-02 *(1.4)/100.);
  sm->SetPoint(i,146.0,1.23E-02); sm->SetPointError(i++,0,1.23E-02 *(1.3)/100.);
  sm->SetPoint(i,147.0,1.34E-02); sm->SetPointError(i++,0,1.34E-02 *(1.3)/100.);
  sm->SetPoint(i,148.0,1.45E-02); sm->SetPointError(i++,0,1.45E-02 *(1.2)/100.);
  sm->SetPoint(i,149.0,1.58E-02); sm->SetPointError(i++,0,1.58E-02 *(1.2)/100.);
  sm->SetPoint(i,150.0,1.73E-02); sm->SetPointError(i++,0,1.73E-02 *(1.1)/100.);
  sm->SetPoint(i,151.0,1.91E-02); sm->SetPointError(i++,0,1.91E-02 *(1.1)/100.);
  sm->SetPoint(i,152.0,2.11E-02); sm->SetPointError(i++,0,2.11E-02 *(1.0)/100.);
  sm->SetPoint(i,153.0,2.36E-02); sm->SetPointError(i++,0,2.36E-02 *(1.0)/100.);
  sm->SetPoint(i,154.0,2.66E-02); sm->SetPointError(i++,0,2.66E-02 *(0.9)/100.);
  sm->SetPoint(i,155.0,3.03E-02); sm->SetPointError(i++,0,3.03E-02 *(0.9)/100.);
  sm->SetPoint(i,156.0,3.51E-02); sm->SetPointError(i++,0,3.51E-02 *(0.8)/100.);
  sm->SetPoint(i,157.0,4.14E-02); sm->SetPointError(i++,0,4.14E-02 *(0.8)/100.);
  sm->SetPoint(i,158.0,5.02E-02); sm->SetPointError(i++,0,5.02E-02 *(0.7)/100.);
  sm->SetPoint(i,159.0,6.32E-02); sm->SetPointError(i++,0,6.32E-02 *(0.7)/100.);
  sm->SetPoint(i,160.0,8.31E-02); sm->SetPointError(i++,0,8.31E-02 *(0.6)/100.);
  sm->SetPoint(i,162.0,1.47E-01); sm->SetPointError(i++,0,1.47E-01 *(0.6)/100.);
  sm->SetPoint(i,164.0,2.15E-01); sm->SetPointError(i++,0,2.15E-01 *(0.6)/100.);
  sm->SetPoint(i,165.0,2.46E-01); sm->SetPointError(i++,0,2.46E-01 *(0.6)/100.);
  sm->SetPoint(i,166.0,2.76E-01); sm->SetPointError(i++,0,2.76E-01 *(0.6)/100.);
  sm->SetPoint(i,168.0,3.30E-01); sm->SetPointError(i++,0,3.30E-01 *(0.6)/100.);
  sm->SetPoint(i,170.0,3.80E-01); sm->SetPointError(i++,0,3.80E-01 *(0.6)/100.);
  sm->SetPoint(i,172.0,4.29E-01); sm->SetPointError(i++,0,4.29E-01 *(0.6)/100.);
  sm->SetPoint(i,174.0,4.77E-01); sm->SetPointError(i++,0,4.77E-01 *(0.6)/100.);
  sm->SetPoint(i,175.0,5.01E-01); sm->SetPointError(i++,0,5.01E-01 *(0.6)/100.);
  sm->SetPoint(i,176.0,5.25E-01); sm->SetPointError(i++,0,5.25E-01 *(0.6)/100.);
  sm->SetPoint(i,178.0,5.75E-01); sm->SetPointError(i++,0,5.75E-01 *(0.6)/100.);
  sm->SetPoint(i,180.0,6.31E-01); sm->SetPointError(i++,0,6.31E-01 *(0.6)/100.);
  sm->SetPoint(i,182.0,7.00E-01); sm->SetPointError(i++,0,7.00E-01 *(0.6)/100.);
  sm->SetPoint(i,184.0,7.88E-01); sm->SetPointError(i++,0,7.88E-01 *(0.6)/100.);
  sm->SetPoint(i,185.0,8.32E-01); sm->SetPointError(i++,0,8.32E-01 *(0.6)/100.);
  sm->SetPoint(i,186.0,8.76E-01); sm->SetPointError(i++,0,8.76E-01 *(0.6)/100.);
  sm->SetPoint(i,188.0,9.60E-01); sm->SetPointError(i++,0,9.60E-01 *(0.5)/100.);
  sm->SetPoint(i,190.0,1.04E+00); sm->SetPointError(i++,0,1.04E+00 *(0.5)/100.);
  sm->SetPoint(i,192.0,1.12E+00); sm->SetPointError(i++,0,1.12E+00 *(0.5)/100.);
  sm->SetPoint(i,194.0,1.20E+00); sm->SetPointError(i++,0,1.20E+00 *(0.5)/100.);
  sm->SetPoint(i,195.0,1.24E+00); sm->SetPointError(i++,0,1.24E+00 *(0.5)/100.);
  sm->SetPoint(i,196.0,1.28E+00); sm->SetPointError(i++,0,1.28E+00 *(0.5)/100.);
  sm->SetPoint(i,198.0,1.35E+00); sm->SetPointError(i++,0,1.35E+00 *(0.5)/100.);
  sm->SetPoint(i,200.0,1.43E+00); sm->SetPointError(i++,0,1.43E+00 *(0.5)/100.);
  sm->SetPoint(i,202.0,1.51E+00); sm->SetPointError(i++,0,1.51E+00 *(0.5)/100.);
  sm->SetPoint(i,204.0,1.59E+00); sm->SetPointError(i++,0,1.59E+00 *(0.5)/100.);
  sm->SetPoint(i,206.0,1.68E+00); sm->SetPointError(i++,0,1.68E+00 *(0.5)/100.);
  sm->SetPoint(i,208.0,1.76E+00); sm->SetPointError(i++,0,1.76E+00 *(0.5)/100.);
  sm->SetPoint(i,210.0,1.85E+00); sm->SetPointError(i++,0,1.85E+00 *(0.5)/100.);
  sm->SetPoint(i,212.0,1.93E+00); sm->SetPointError(i++,0,1.93E+00 *(0.5)/100.);
  sm->SetPoint(i,214.0,2.02E+00); sm->SetPointError(i++,0,2.02E+00 *(0.5)/100.);
  sm->SetPoint(i,216.0,2.12E+00); sm->SetPointError(i++,0,2.12E+00 *(0.5)/100.);
  sm->SetPoint(i,218.0,2.21E+00); sm->SetPointError(i++,0,2.21E+00 *(0.5)/100.);
  sm->SetPoint(i,220.0,2.31E+00); sm->SetPointError(i++,0,2.31E+00 *(0.5)/100.);
  sm->SetPoint(i,222.0,2.40E+00); sm->SetPointError(i++,0,2.40E+00 *(0.5)/100.);
  sm->SetPoint(i,224.0,2.50E+00); sm->SetPointError(i++,0,2.50E+00 *(0.5)/100.);
  sm->SetPoint(i,226.0,2.61E+00); sm->SetPointError(i++,0,2.61E+00 *(0.5)/100.);
  sm->SetPoint(i,228.0,2.71E+00); sm->SetPointError(i++,0,2.71E+00 *(0.5)/100.);
  sm->SetPoint(i,230.0,2.82E+00); sm->SetPointError(i++,0,2.82E+00 *(0.5)/100.);
  sm->SetPoint(i,232.0,2.93E+00); sm->SetPointError(i++,0,2.93E+00 *(0.5)/100.);
  sm->SetPoint(i,234.0,3.04E+00); sm->SetPointError(i++,0,3.04E+00 *(0.5)/100.);
  sm->SetPoint(i,236.0,3.16E+00); sm->SetPointError(i++,0,3.16E+00 *(0.5)/100.);
  sm->SetPoint(i,238.0,3.27E+00); sm->SetPointError(i++,0,3.27E+00 *(0.5)/100.);
  sm->SetPoint(i,240.0,3.40E+00); sm->SetPointError(i++,0,3.40E+00 *(0.5)/100.);
  sm->SetPoint(i,242.0,3.52E+00); sm->SetPointError(i++,0,3.52E+00 *(0.5)/100.);
  sm->SetPoint(i,244.0,3.64E+00); sm->SetPointError(i++,0,3.64E+00 *(0.5)/100.);
  sm->SetPoint(i,246.0,3.77E+00); sm->SetPointError(i++,0,3.77E+00 *(0.5)/100.);
  sm->SetPoint(i,248.0,3.91E+00); sm->SetPointError(i++,0,3.91E+00 *(0.5)/100.);
  sm->SetPoint(i,250.0,4.04E+00); sm->SetPointError(i++,0,4.04E+00 *(0.5)/100.);
  sm->SetPoint(i,252.0,4.18E+00); sm->SetPointError(i++,0,4.18E+00 *(0.5)/100.);
  sm->SetPoint(i,254.0,4.32E+00); sm->SetPointError(i++,0,4.32E+00 *(0.5)/100.);
  sm->SetPoint(i,256.0,4.46E+00); sm->SetPointError(i++,0,4.46E+00 *(0.5)/100.);
  sm->SetPoint(i,258.0,4.61E+00); sm->SetPointError(i++,0,4.61E+00 *(0.5)/100.);
  sm->SetPoint(i,260.0,4.76E+00); sm->SetPointError(i++,0,4.76E+00 *(0.5)/100.);
  sm->SetPoint(i,262.0,4.91E+00); sm->SetPointError(i++,0,4.91E+00 *(0.5)/100.);
  sm->SetPoint(i,264.0,5.07E+00); sm->SetPointError(i++,0,5.07E+00 *(0.5)/100.);
  sm->SetPoint(i,266.0,5.23E+00); sm->SetPointError(i++,0,5.23E+00 *(0.5)/100.);
  sm->SetPoint(i,268.0,5.39E+00); sm->SetPointError(i++,0,5.39E+00 *(0.5)/100.);
  sm->SetPoint(i,270.0,5.55E+00); sm->SetPointError(i++,0,5.55E+00 *(0.5)/100.);
  sm->SetPoint(i,272.0,5.72E+00); sm->SetPointError(i++,0,5.72E+00 *(0.5)/100.);
  sm->SetPoint(i,274.0,5.89E+00); sm->SetPointError(i++,0,5.89E+00 *(0.5)/100.);
  sm->SetPoint(i,276.0,6.07E+00); sm->SetPointError(i++,0,6.07E+00 *(0.5)/100.); 
  sm->SetPoint(i,278.0,6.25E+00); sm->SetPointError(i++,0,6.25E+00 *(0.5)/100.); 
  sm->SetPoint(i,280.0,6.43E+00); sm->SetPointError(i++,0,6.43E+00 *(0.5)/100.); 
  sm->SetPoint(i,282.0,6.61E+00); sm->SetPointError(i++,0,6.61E+00 *(0.5)/100.); 
  sm->SetPoint(i,284.0,6.80E+00); sm->SetPointError(i++,0,6.80E+00 *(0.5)/100.); 
  sm->SetPoint(i,286.0,6.99E+00); sm->SetPointError(i++,0,6.99E+00 *(0.5)/100.); 
  sm->SetPoint(i,288.0,7.19E+00); sm->SetPointError(i++,0,7.19E+00 *(0.5)/100.); 
  sm->SetPoint(i,290.0,7.39E+00); sm->SetPointError(i++,0,7.39E+00 *(0.5)/100.); 
  sm->SetPoint(i,295.0,7.90E+00); sm->SetPointError(i++,0,7.90E+00 *(0.5)/100.); 
  sm->SetPoint(i,300.0,8.43E+00); sm->SetPointError(i++,0,8.43E+00 *(0.5)/100.); 
  sm->SetPoint(i,305.0,8.99E+00); sm->SetPointError(i++,0,8.99E+00 *(0.5)/100.); 
  sm->SetPoint(i,310.0,9.57E+00); sm->SetPointError(i++,0,9.57E+00 *(0.5)/100.); 
  sm->SetPoint(i,315.0,1.02E+01); sm->SetPointError(i++,0,1.02E+01 *(0.5)/100.); 
  sm->SetPoint(i,320.0,1.08E+01); sm->SetPointError(i++,0,1.08E+01 *(0.5)/100.); 
  sm->SetPoint(i,325.0,1.14E+01); sm->SetPointError(i++,0,1.14E+01 *(0.6)/100.); 
  sm->SetPoint(i,330.0,1.21E+01); sm->SetPointError(i++,0,1.21E+01 *(0.6)/100.); 
  sm->SetPoint(i,335.0,1.28E+01); sm->SetPointError(i++,0,1.28E+01 *(0.7)/100.); 
  sm->SetPoint(i,340.0,1.35E+01); sm->SetPointError(i++,0,1.35E+01 *(0.7)/100.); 
  sm->SetPoint(i,345.0,1.42E+01); sm->SetPointError(i++,0,1.42E+01 *(1.4)/100.); 
  sm->SetPoint(i,350.0,1.52E+01); sm->SetPointError(i++,0,1.52E+01 *(2.1)/100.); 
  sm->SetPoint(i,360.0,1.76E+01); sm->SetPointError(i++,0,1.76E+01 *(3.4)/100.); 
  sm->SetPoint(i,370.0,2.02E+01); sm->SetPointError(i++,0,2.02E+01 *(3.7)/100.); 
  sm->SetPoint(i,380.0,2.31E+01); sm->SetPointError(i++,0,2.31E+01 *(3.7)/100.); 
  sm->SetPoint(i,390.0,2.61E+01); sm->SetPointError(i++,0,2.61E+01 *(3.7)/100.); 
  sm->SetPoint(i,400.0,2.92E+01); sm->SetPointError(i++,0,2.92E+01 *(3.6)/100.); 
  sm->SetPoint(i,410.0,3.25E+01); sm->SetPointError(i++,0,3.25E+01 *(3.6)/100.); 
  sm->SetPoint(i,420.0,3.59E+01); sm->SetPointError(i++,0,3.59E+01 *(3.5)/100.); 
  sm->SetPoint(i,430.0,3.94E+01); sm->SetPointError(i++,0,3.94E+01 *(3.4)/100.); 
  sm->SetPoint(i,440.0,4.30E+01); sm->SetPointError(i++,0,4.30E+01 *(3.3)/100.); 
  sm->SetPoint(i,450.0,4.68E+01); sm->SetPointError(i++,0,4.68E+01 *(3.2)/100.); 
  sm->SetPoint(i,460.0,5.08E+01); sm->SetPointError(i++,0,5.08E+01 *(3.1)/100.); 
  sm->SetPoint(i,470.0,5.49E+01); sm->SetPointError(i++,0,5.49E+01 *(3.1)/100.); 
  sm->SetPoint(i,480.0,5.91E+01); sm->SetPointError(i++,0,5.91E+01 *(3.0)/100.); 
  sm->SetPoint(i,490.0,6.35E+01); sm->SetPointError(i++,0,6.35E+01 *(2.9)/100.); 
  sm->SetPoint(i,500.0,6.80E+01); sm->SetPointError(i++,0,6.80E+01 *(3.0)/100.); 
  sm->SetPoint(i,510.0,7.27E+01); sm->SetPointError(i++,0,7.27E+01 *(3.0)/100.); 
  sm->SetPoint(i,520.0,7.75E+01); sm->SetPointError(i++,0,7.75E+01 *(3.1)/100.); 
  sm->SetPoint(i,530.0,8.25E+01); sm->SetPointError(i++,0,8.25E+01 *(3.3)/100.); 
  sm->SetPoint(i,540.0,8.77E+01); sm->SetPointError(i++,0,8.77E+01 *(3.4)/100.); 
  sm->SetPoint(i,550.0,9.30E+01); sm->SetPointError(i++,0,9.30E+01 *(3.5)/100.); 
  sm->SetPoint(i,560.0,9.86E+01); sm->SetPointError(i++,0,9.86E+01 *(3.6)/100.); 
  sm->SetPoint(i,570.0,1.04E+02); sm->SetPointError(i++,0,1.04E+02 *(3.7)/100.); 
  sm->SetPoint(i,580.0,1.10E+02); sm->SetPointError(i++,0,1.10E+02 *(3.8)/100.); 
  sm->SetPoint(i,590.0,1.16E+02); sm->SetPointError(i++,0,1.16E+02 *(3.9)/100.); 
  sm->SetPoint(i,600.0,1.23E+02); sm->SetPointError(i++,0,1.23E+02 *(4.0)/100.); 
  sm->SetPoint(i,610.0,1.29E+02); sm->SetPointError(i++,0,1.29E+02 *(4.2)/100.); 
  sm->SetPoint(i,620.0,1.36E+02); sm->SetPointError(i++,0,1.36E+02 *(4.3)/100.); 
  sm->SetPoint(i,630.0,1.43E+02); sm->SetPointError(i++,0,1.43E+02 *(4.5)/100.); 
  sm->SetPoint(i,640.0,1.50E+02); sm->SetPointError(i++,0,1.50E+02 *(4.6)/100.); 
  sm->SetPoint(i,650.0,1.58E+02); sm->SetPointError(i++,0,1.58E+02 *(4.8)/100.); 
  sm->SetPoint(i,660.0,1.65E+02); sm->SetPointError(i++,0,1.65E+02 *(4.9)/100.); 
  sm->SetPoint(i,670.0,1.73E+02); sm->SetPointError(i++,0,1.73E+02 *(5.1)/100.); 
  sm->SetPoint(i,680.0,1.82E+02); sm->SetPointError(i++,0,1.82E+02 *(5.3)/100.); 
  sm->SetPoint(i,690.0,1.90E+02); sm->SetPointError(i++,0,1.90E+02 *(5.5)/100.); 
  sm->SetPoint(i,700.0,1.99E+02); sm->SetPointError(i++,0,1.99E+02 *(5.6)/100.); 
  sm->SetPoint(i,710.0,2.08E+02); sm->SetPointError(i++,0,2.08E+02 *(5.9)/100.); 
  sm->SetPoint(i,720.0,2.17E+02); sm->SetPointError(i++,0,2.17E+02 *(6.1)/100.); 
  sm->SetPoint(i,730.0,2.27E+02); sm->SetPointError(i++,0,2.27E+02 *(6.3)/100.); 
  sm->SetPoint(i,740.0,2.37E+02); sm->SetPointError(i++,0,2.37E+02 *(6.5)/100.); 
  sm->SetPoint(i,750.0,2.47E+02); sm->SetPointError(i++,0,2.47E+02 *(6.8)/100.); 
  sm->SetPoint(i,760.0,2.58E+02); sm->SetPointError(i++,0,2.58E+02 *(7.0)/100.); 
  sm->SetPoint(i,770.0,2.69E+02); sm->SetPointError(i++,0,2.69E+02 *(7.3)/100.); 
  sm->SetPoint(i,780.0,2.80E+02); sm->SetPointError(i++,0,2.80E+02 *(7.6)/100.); 
  sm->SetPoint(i,790.0,2.92E+02); sm->SetPointError(i++,0,2.92E+02 *(7.9)/100.); 
  sm->SetPoint(i,800.0,3.04E+02); sm->SetPointError(i++,0,3.04E+02 *(8.2)/100.); 
  sm->SetPoint(i,810.0,3.17E+02); sm->SetPointError(i++,0,3.17E+02 *(8.6)/100.); 
  sm->SetPoint(i,820.0,3.30E+02); sm->SetPointError(i++,0,3.30E+02 *(8.9)/100.); 
  sm->SetPoint(i,830.0,3.43E+02); sm->SetPointError(i++,0,3.43E+02 *(9.2)/100.); 
  sm->SetPoint(i,840.0,3.57E+02); sm->SetPointError(i++,0,3.57E+02 *(9.6)/100.); 
  sm->SetPoint(i,850.0,3.71E+02); sm->SetPointError(i++,0,3.71E+02 *(10.0)/100.);
  sm->SetPoint(i,860.0,3.86E+02); sm->SetPointError(i++,0,3.86E+02 *(10.4)/100.);
  sm->SetPoint(i,870.0,4.01E+02); sm->SetPointError(i++,0,4.01E+02 *(10.8)/100.);
  sm->SetPoint(i,880.0,4.16E+02); sm->SetPointError(i++,0,4.16E+02 *(11.3)/100.);
  sm->SetPoint(i,890.0,4.32E+02); sm->SetPointError(i++,0,4.32E+02 *(11.7)/100.);
  sm->SetPoint(i,900.0,4.49E+02); sm->SetPointError(i++,0,4.49E+02 *(12.2)/100.);
  sm->SetPoint(i,910.0,4.66E+02); sm->SetPointError(i++,0,4.66E+02 *(12.7)/100.);
  sm->SetPoint(i,920.0,4.84E+02); sm->SetPointError(i++,0,4.84E+02 *(13.2)/100.);
  sm->SetPoint(i,930.0,5.02E+02); sm->SetPointError(i++,0,5.02E+02 *(13.7)/100.);
  sm->SetPoint(i,940.0,5.21E+02); sm->SetPointError(i++,0,5.21E+02 *(14.2)/100.);
  sm->SetPoint(i,950.0,5.40E+02); sm->SetPointError(i++,0,5.40E+02 *(14.8)/100.);
  sm->SetPoint(i,960.0,5.60E+02); sm->SetPointError(i++,0,5.60E+02 *(15.3)/100.);
  sm->SetPoint(i,970.0,5.81E+02); sm->SetPointError(i++,0,5.81E+02 *(15.9)/100.);
  sm->SetPoint(i,980.0,6.02E+02); sm->SetPointError(i++,0,6.02E+02 *(16.6)/100.);
  sm->SetPoint(i,990.0,6.24E+02); sm->SetPointError(i++,0,6.24E+02 *(17.2)/100.);
  sm->SetPoint(i,1000.0,6.47E+02);sm->SetPointError(i++,0,6.47E+02 *(17.8)/100.);
  return sm;
};