void outopen() {
  gROOT->Reset();
  gStyle->SetOptStat(0);

  TString DataSetName = "SANE2_new_CH"; // SANE1(1.05,19), SANE2(2.01,11), SANE3(2.04,9)
  Double_t Wmin = 2.01;//1.070;//1.050;//1.155;//2.025;//1.05;//2.01;//1.05;//2.0;//0.90;//1.05; //To apply dilution factor in [0.0,3.0] with 200 bins, Wmin should be divided by 0.015 GeV.
  Int_t Wbin = 22;//19;//28;//19;//34;//18;//40;//22;//33;//75; // For 15MeV resoultion deltaW=1.05 for Wbin=70, 30MeV for Wbin=35
  Double_t Wbinsize = 0.015; //0.010; // 0.015 for 15MeV, 0.030 for 30MeV
  Double_t Wmax = Wmin + Wbin*Wbinsize;
  TString DataQ2Name = "1.31";
  TString CLASQ2Name = "1.2";//"0.844";//"1.71";//"1.2";//"1.71";
  Double_t CLAS_Q2 = 1.2;//0.844;//1.2;//1.71;


  TObjArray Hlist(0);

  ifstream outfile;
  outfile.open("rss_proton_b1.dat");
  
  TString line;
  Int_t i = 0;
  Int_t j = 0;
  const Int_t n = 50;
  Double_t W[n],W_mid[n],x_Bj[n],Q2[n];
  Double_t Apar_data[n],eApar_data[n],Apar_Born[n],eApar_Born[n];
  Double_t Aper_data[n],eAper_data[n],Aper_Born[n],eAper_Born[n];
  Double_t A1_data[n],eA1_data[n],A1_Born[n],eA1_Born[n];
  Double_t A2_data[n],eA2_data[n],A2_Born[n],eA2_Born[n];
  Double_t g1_Born[n],eg1_Born[n];
  Double_t g2_Born[n],eg2_Born[n];
  Double_t g2WW_Born[n];
  Double_t sigTT_Born[n],esigTT_Born[n];
  Double_t sigLT_Born[n],esigLT_Born[n];
  Double_t smA1[n],smA2[n];

  while(!outfile.eof()) {
    
    line.ReadLine(outfile);

    if(!line.Contains("#") && !line.IsNull()) {

      if(line_elem(line,0)==0) {
	Double_t W_kin = line_elem(line,49);
	Double_t Apar = line_elem(line,8);
	Double_t eApar = line_elem(line,9);
	Double_t Aper = line_elem(line,10);
	Double_t eAper = line_elem(line,11);
	Double_t A1 = line_elem(line,20);
	Double_t eA1 = line_elem(line,21);
	Double_t A2 = line_elem(line,22);
	Double_t eA2 = line_elem(line,23);
	Double_t g1 = line_elem(line,16);
	Double_t eg1 = line_elem(line,17);
	Double_t g2 = line_elem(line,18);
	Double_t eg2 = line_elem(line,19);
	Double_t gT = line_elem(line,43);
	Double_t egT = line_elem(line,44);

        W[i] = line_elem(line,2)/1.e3;
        W_mid[i] = line_elem(line,49)/1.e3;
        x_Bj[i] = line_elem(line,4);
        Q2[i] = line_elem(line,3);
        Apar_data[i] = line_elem(line,8);
        eApar_data[i] = line_elem(line,9);
        Aper_data[i] = line_elem(line,10);
        eAper_data[i] = line_elem(line,11);
        A1_data[i] = line_elem(line,20);
        eA1_data[i] = line_elem(line,21);
        A2_data[i] = line_elem(line,22);
        eA2_data[i] = line_elem(line,23);

        i++;
      }

      if(line_elem(line,0)==2) {
	Double_t W_kin = line_elem(line,49);
	Double_t Apar = line_elem(line,8);
	Double_t eApar = line_elem(line,9);
	Double_t Aper = line_elem(line,10);
	Double_t eAper = line_elem(line,11);
	Double_t A1 = line_elem(line,20);
	Double_t eA1 = line_elem(line,21);
	Double_t A2 = line_elem(line,22);
	Double_t eA2 = line_elem(line,23);
	Double_t g1 = line_elem(line,16);
	Double_t eg1 = line_elem(line,17);
	Double_t g2 = line_elem(line,18);
	Double_t eg2 = line_elem(line,19);
	Double_t gT = line_elem(line,43);
	Double_t egT = line_elem(line,44);
	Double_t sigTT = line_elem(line,26);
	Double_t esigTT = line_elem(line,27);
	Double_t sigLT = line_elem(line,28);
	Double_t esigLT = line_elem(line,29);
        smA1[j] = line_elem(line,59); // Systematic error of model dependence
        smA2[j] = line_elem(line,60);

        Apar_Born[j] = line_elem(line,8);
        eApar_Born[j] = line_elem(line,9);
        Aper_Born[j] = line_elem(line,10);
        eAper_Born[j] = line_elem(line,11);
        A1_Born[j] = line_elem(line,20);
        eA1_Born[j] = line_elem(line,21);
        A2_Born[j] = line_elem(line,22);
        eA2_Born[j] = line_elem(line,23);
        g1_Born[j] = line_elem(line,16);
        eg1_Born[j] = line_elem(line,17);
        g2_Born[j] = line_elem(line,18);
        eg2_Born[j] = line_elem(line,19);
        g2WW_Born[j] = line_elem(line,42);
        sigTT_Born[j] = line_elem(line,26);
        esigTT_Born[j] = line_elem(line,27);
        sigLT_Born[j] = line_elem(line,28);
        esigLT_Born[j] = line_elem(line,29);

        j++;
      }
    }
  }
  outfile.close();

// Systematics
  ifstream SEfile;
  SEfile.open("sane_systematics.txt");
  
  TString line;
  Int_t k = 0;
  const Int_t n = 50;
  Double_t seA1[n],seA2[n],seg1[n],seg2[n];
  Double_t A1pos[n],A2pos[n],g1pos[n],g2pos[n];

  while(!SEfile.eof()) {

    line.ReadLine(SEfile);

    if(!line.Contains("#") && !line.IsNull()) {
      seA1[k] = line_elem(line,8);
      seA2[k] = line_elem(line,10);
      seg1[k] = line_elem(line,12);
      seg2[k] = line_elem(line,14);
      //cout << seA1[k] << seg1[k] << endl;
      A1pos[k] = -0.6; A2pos[k] = -1.0; g1pos[k] = -0.1; g2pos[k] = -0.15;
      k++;      
    }
  }

  if(i!=k||j!=k) cout << "systematic file error" << endl;


  TGraphErrors *Apar_noRC = new TGraphErrors(i,W,Apar_data,NULL,NULL);
  TGraphErrors *Aper_noRC = new TGraphErrors(i,W,Aper_data,NULL,NULL);
  TGraphErrors *Apar_RC = new TGraphErrors(j,W,Apar_Born,NULL,eApar_Born);
  TGraphErrors *Aper_RC = new TGraphErrors(j,W,Aper_Born,NULL,eAper_Born);
  TGraphErrors *A1_noRC = new TGraphErrors(i,W,A1_data,NULL,NULL);
  TGraphErrors *A2_noRC = new TGraphErrors(i,W,A2_data,NULL,NULL);
  TGraphErrors *A1x_noRC = new TGraphErrors(i,x_Bj,A1_data,NULL,NULL);
  TGraphErrors *A2x_noRC = new TGraphErrors(i,x_Bj,A2_data,NULL,NULL);
  TGraphErrors *A1_RC = new TGraphErrors(j,W,A1_Born,NULL,eA1_Born);
  TGraphErrors *A2_RC = new TGraphErrors(j,W,A2_Born,NULL,eA2_Born);
  TGraph2DErrors *A1_WQ2 = new TGraph2DErrors(j,W,Q2,A1_Born,NULL,NULL,eA1_Born);
  TGraph2DErrors *A2_WQ2 = new TGraph2DErrors(j,W,Q2,A2_Born,NULL,NULL,eA2_Born);
  TGraphErrors *A1x_RC = new TGraphErrors(j,x_Bj,A1_Born,NULL,eA1_Born);
  TGraphErrors *A2x_RC = new TGraphErrors(j,x_Bj,A2_Born,NULL,eA2_Born);
  TGraphErrors *g1_RC = new TGraphErrors(j,x_Bj,g1_Born,NULL,eg1_Born);
  TGraphErrors *g2_RC = new TGraphErrors(j,x_Bj,g2_Born,NULL,eg2_Born);
  TGraphErrors *g1W_RC = new TGraphErrors(j,W,g1_Born,NULL,eg1_Born);
  TGraphErrors *g2W_RC = new TGraphErrors(j,W,g2_Born,NULL,eg2_Born);
  TGraphErrors *g2WW_RC = new TGraphErrors(j,x_Bj,g2WW_Born,NULL,NULL);
  TGraphErrors *sigTT_RC = new TGraphErrors(j,W,sigTT_Born,NULL,esigTT_Born);
  TGraphErrors *sigLT_RC = new TGraphErrors(j,W,sigLT_Born,NULL,esigLT_Born);

  TGraphAsymmErrors *A1sysb_RC = new TGraphAsymmErrors(j,W,A1pos,NULL,NULL,NULL,seA1);
  TGraphErrors *A1sys_RC = new TGraphErrors(j,W,A1_Born,NULL,seA1);
  TGraphErrors *A2sys_RC = new TGraphErrors(j,W,A2_Born,NULL,seA2);
  TGraphErrors *A1modsys_RC = new TGraphErrors(j,W,A1_Born,NULL,smA1);
  TGraphErrors *A2modsys_RC = new TGraphErrors(j,W,A2_Born,NULL,smA2);
  TGraphErrors *A1xsys_RC = new TGraphErrors(j,x_Bj,A1_Born,NULL,seA1);
  TGraphErrors *A2xsys_RC = new TGraphErrors(j,x_Bj,A2_Born,NULL,seA2);
  TGraphErrors *g1sys_RC = new TGraphErrors(j,x_Bj,g1_Born,NULL,seg1);
  TGraphErrors *g2sys_RC = new TGraphErrors(j,x_Bj,g2_Born,NULL,seg2);

  TF1 *hzl=new TF1("hzl","0.",0.,10.);
  hzl->SetLineStyle(3); hzl->SetLineColor(kBlack);

  TCanvas *W_Apar_c = new TCanvas("W_Apar_c","W_Apar");
  W_Apar_c->SetGrid();
  W_Apar_c->cd();
  Apar_noRC->SetTitle("A_{180} along W");
  Apar_noRC->GetXaxis()->SetTitle("W [GeV]");
  Apar_noRC->SetMaximum(0.2); Apar_noRC->SetMinimum(-0.2);
  Apar_noRC->Draw("ac");
  Apar_RC->SetMarkerStyle(20); Apar_RC->SetMarkerColor(kBlue); Apar_RC->SetLineColor(kBlue);
  Apar_RC->Draw("p");
  legApar = new TLegend(0.5,0.15,0.8,0.25);
  legApar->AddEntry(Apar_noRC,"This work without r.c","l");
  legApar->AddEntry(Apar_RC,"This work with r.c","lep");
  legApar->Draw();
  W_Apar_c->Print(DataSetName+"_A180.png");

  TCanvas *W_Aper_c = new TCanvas("W_Aper_c","W_Aper");
  W_Aper_c->SetGrid();
  W_Aper_c->cd();
  Aper_noRC->SetTitle("A_{80} along W");
  Aper_noRC->GetXaxis()->SetTitle("W [GeV]");
  Aper_noRC->SetMaximum(0.2); Aper_noRC->SetMinimum(-0.2);
  Aper_noRC->Draw("ac");
  Aper_RC->SetMarkerStyle(20); Aper_RC->SetMarkerColor(kBlue); Aper_RC->SetLineColor(kBlue); 
  Aper_RC->Draw("p");
  legAper = new TLegend(0.5,0.6,0.8,0.7);
  legAper->AddEntry(Aper_noRC,"This work without r.c","l");
  legAper->AddEntry(Aper_RC,"This work with r.c","lep");
  legAper->Draw();
  W_Aper_c->Print(DataSetName+"_A80.png");

  TCanvas *W_A180_A80_c = new TCanvas("W_A180_A80_c","W_A180_A80");
  //W_A180_A80_c->SetGridy();
  W_A180_A80_c->cd();
  Apar_RC->SetTitle(NULL);
  Apar_RC->GetXaxis()->SetTitle("W [GeV]");
  Apar_RC->SetMaximum(1.0); Apar_RC->SetMinimum(-1.0);
  Apar_RC->SetMarkerStyle(20); Apar_RC->SetMarkerColor(kBlue); Apar_RC->SetLineColor(kBlue);
  Apar_RC->Draw("ap");
  Apar_noRC->Draw("c");
  Aper_RC->SetMarkerStyle(21); Aper_RC->SetMarkerColor(kRed); Aper_RC->SetLineColor(kRed);
  Aper_RC->Draw("p");
  Aper_noRC->Draw("c");
  hzl->Draw("same");
  legA180_A80 = new TLegend(0.5,0.15,0.6,0.25);
  legA180_A80->AddEntry(Apar_RC,"A_{#parallel}","lep");
  legA180_A80->AddEntry(Aper_RC,"A_{#perp}","lep");
  legA180_A80->Draw();
  W_A180_A80_c->Print(DataSetName+"_A180_A80.png");

  //TCanvas *W_g1_c = new TCanvas("W_g1_c","W_g1");
  //W_g1_c->SetGrid();
  //W_g1_c->cd();
  //W_g1_RC->SetMaximum(0.2); W_g1_RC->SetMinimum(-0.05);
  //W_g1_noRC->SetMarkerStyle(20); W_g1_noRC->Draw("C");
  //W_g1_RC->SetMarkerStyle(20); W_g1_RC->SetMarkerColor(kBlue); W_g1_RC->SetLineColor(kBlue);
  //W_g1_RC->Draw();
  //W_g1_c->Print(DataSetName+"_g1.png");

  //TCanvas *W_g2_c = new TCanvas("W_g2_c","W_g2");
  //W_g2_c->SetGrid();
  //W_g2_c->cd();
  //W_g2_RC->SetMaximum(0.03); W_g2_RC->SetMinimum(-0.15);
  //W_g2_noRC->SetMarkerStyle(20); W_g2_noRC->Draw("C");
  //W_g2_RC->SetMarkerStyle(20); W_g2_RC->SetMarkerColor(kBlue); W_g2_RC->SetLineColor(kBlue);
  //W_g2_RC->Draw();
  //W_g2_c->Print(DataSetName+"_g2.png");

  // ***********
  // CLAS result
  // ***********

  ifstream CLASfile1;
  CLASfile1.open("clasdb_E93M7.txt");
  
  TString line;
  Int_t i = 0;
  Int_t j = 0;
  Int_t k = 0;
  Int_t l = 0;
  const Int_t n = 50;
  Double_t W_CLAS[n],x_Bj_CLAS[n],A1_CLAS[n],eA1_CLAS[n];
  Double_t W_CLASl[n],x_Bj_CLASl[n],A1_CLASl[n],eA1_CLASl[n];
  Double_t W_CLASh[n],x_Bj_CLASh[n],A1_CLASh[n],eA1_CLASh[n];
  Double_t W_CLASll[n],x_Bj_CLASll[n],A1_CLASll[n],eA1_CLASll[n];
  
  while(!CLASfile1.eof()) {

    line.ReadLine(CLASfile1);

    if(!line.Contains("#") && !line.IsNull()) {
      if(line_elem(line,2)==1.2) { //CLAS_Q2
      W_CLAS[i] = line_elem(line,0);
      x_Bj_CLAS[i] = line_elem(line,1);
      A1_CLAS[i] = line_elem(line,3);
      eA1_CLAS[i] = line_elem(line,4);
      //cout << W_CLAS[i] << A1_CLAS[i] << endl;
      i++;
      }
      if(line_elem(line,2)==1.44) {
      W_CLASl[j] = line_elem(line,0);
      x_Bj_CLASl[j] = line_elem(line,1);
      A1_CLASl[j] = line_elem(line,3);
      eA1_CLASl[j] = line_elem(line,4);
      j++;
      }
      if(line_elem(line,2)==2.05) {
      W_CLASh[k] = line_elem(line,0);
      x_Bj_CLASh[k] = line_elem(line,1);
      A1_CLASh[k] = line_elem(line,3);
      eA1_CLASh[k] = line_elem(line,4);
      k++;
      }
      if(line_elem(line,2)==1.2) {
      W_CLASll[l] = line_elem(line,0);
      x_Bj_CLASll[l] = line_elem(line,1);
      A1_CLASll[l] = line_elem(line,3);
      eA1_CLASll[l] = line_elem(line,4);
      l++;
      }
    }
  }

  TGraphErrors *A1_CLAS_gr = new TGraphErrors(i,W_CLAS,A1_CLAS,NULL,eA1_CLAS);
  TGraphErrors *A1_CLASl_gr = new TGraphErrors(j,W_CLASl,A1_CLASl,NULL,eA1_CLASl);
  TGraphErrors *A1_CLASh_gr = new TGraphErrors(k,W_CLASh,A1_CLASh,NULL,eA1_CLASh);
  TGraphErrors *A1_CLASll_gr = new TGraphErrors(l,W_CLASll,A1_CLASll,NULL,eA1_CLASll);

  TGraphErrors *A1x_CLAS_gr = new TGraphErrors(i,x_Bj_CLAS,A1_CLAS,NULL,eA1_CLAS);
  TGraphErrors *A1x_CLASl_gr = new TGraphErrors(j,x_Bj_CLASl,A1_CLASl,NULL,eA1_CLASl);
  TGraphErrors *A1x_CLASh_gr = new TGraphErrors(k,x_Bj_CLASh,A1_CLASh,NULL,eA1_CLASh);
  TGraphErrors *A1x_CLASll_gr = new TGraphErrors(l,x_Bj_CLASll,A1_CLASll,NULL,eA1_CLASll);

  ifstream CLASfile2;
  CLASfile2.open("clasdb_E93M9.txt");

  TString line;
  Int_t i = 0;
  Int_t j = 0;
  Int_t k = 0;
  const Int_t n = 50;
  Double_t W_CLAS[n],x_Bj_CLAS[n],g1_CLAS[n],eg1_CLAS[n];
  Double_t x_Bj_CLASl[n],g1_CLASl[n],eg1_CLASl[n];
  Double_t x_Bj_CLASh[n],g1_CLASh[n],eg1_CLASh[n];

  while(!CLASfile2.eof()) {

    line.ReadLine(CLASfile2);

    if(!line.Contains("#") && !line.IsNull()) {
      if(line_elem(line,2)==1.2) { //CLAS_Q2
      W_CLAS[i] = line_elem(line,0);
      x_Bj_CLAS[i] = line_elem(line,1);
      g1_CLAS[i] = line_elem(line,3);
      eg1_CLAS[i] = line_elem(line,4);
      i++;
      }
      if(line_elem(line,2)==1.44) {
      x_Bj_CLASl[j] = line_elem(line,1);
      g1_CLASl[j] = line_elem(line,3);
      eg1_CLASl[j] = line_elem(line,4);
      j++;
      }
      if(line_elem(line,2)==2.05) {
      x_Bj_CLASh[k] = line_elem(line,1);
      g1_CLASh[k] = line_elem(line,3);
      eg1_CLASh[k] = line_elem(line,4);
      k++;
      }
    }
  }

  TGraphErrors *g1_CLAS_gr = new TGraphErrors(i,x_Bj_CLAS,g1_CLAS,NULL,eg1_CLAS);
  TGraphErrors *g1_CLASl_gr = new TGraphErrors(j,x_Bj_CLASl,g1_CLASl,NULL,eg1_CLASl);
  TGraphErrors *g1_CLASh_gr = new TGraphErrors(k,x_Bj_CLASh,g1_CLASh,NULL,eg1_CLASh);

  TGraphErrors *g1W_CLAS_gr = new TGraphErrors(i,W_CLAS,g1_CLAS,NULL,eg1_CLAS);

  // **********
  // RSS result
  // **********

  Double_t Wmin = 1.075;
  Int_t Wbin = 28;
  Double_t Wbinsize = 0.030;
  Double_t Wmax = Wmin + Wbin*Wbinsize;

  ifstream RSSfile;
  RSSfile.open("rss_p_v6c_Born.txt");
  
  TString line;
  Int_t i = 0;
  const Int_t n = 50;
  Double_t W_RSS[n],x_Bj_RSS[n],A1_RSS[n],eA1_RSS[n],A2_RSS[n],eA2_RSS[n],g1_RSS[n],eg1_RSS[n],g2_RSS[n],eg2_RSS[n],sigTT_RSS[n],esigTT_RSS[n],sigLT_RSS[n],esigLT_RSS[n];

  while(!RSSfile.eof()) {
    
    line.ReadLine(RSSfile);

    if(!line.Contains("#") && !line.IsNull()) {

      if(line_elem(line,0)==2) {
	Double_t W_kin = line_elem(line,49);
	Double_t A1 = line_elem(line,20);
	Double_t eA1 = line_elem(line,21);
	Double_t A2 = line_elem(line,22);
	Double_t eA2 = line_elem(line,23);
	Double_t Apar = line_elem(line,8);
	Double_t eApar = line_elem(line,9);
	Double_t Aper = line_elem(line,10);
	Double_t eAper = line_elem(line,11);
	Double_t g1 = line_elem(line,16);
	Double_t eg1 = line_elem(line,17);
	Double_t g2 = line_elem(line,18);
	Double_t eg2 = line_elem(line,19);
	Double_t sigTT = line_elem(line,26);
	Double_t esigTT = line_elem(line,27);
	Double_t sigLT = line_elem(line,28);
	Double_t esigLT = line_elem(line,29);

        W_RSS[i] = line_elem(line,2)/1.e3;
        x_Bj_RSS[i] = line_elem(line,4);
        A1_RSS[i] = line_elem(line,20);
        eA1_RSS[i] = line_elem(line,21);
        A2_RSS[i] = line_elem(line,22);
        eA2_RSS[i] = line_elem(line,23);   
        g1_RSS[i] = line_elem(line,16);
        eg1_RSS[i] = line_elem(line,17);
        g2_RSS[i] = line_elem(line,18);
        eg2_RSS[i] = line_elem(line,19);
        sigTT_RSS[i] = line_elem(line,26);
	esigTT_RSS[i] = line_elem(line,27);
	sigLT_RSS[i] = line_elem(line,28);
	esigLT_RSS[i] = line_elem(line,29);

        i++;
      }
    }
  }

  TGraphErrors *A1_RSS_gr = new TGraphErrors(i,W_RSS,A1_RSS,NULL,eA1_RSS);
  TGraphErrors *A2_RSS_gr = new TGraphErrors(i,W_RSS,A2_RSS,NULL,eA2_RSS);
  TGraphErrors *A1x_RSS_gr = new TGraphErrors(i,x_Bj_RSS,A1_RSS,NULL,eA1_RSS);
  TGraphErrors *A2x_RSS_gr = new TGraphErrors(i,x_Bj_RSS,A2_RSS,NULL,eA2_RSS);
  TGraphErrors *g1_RSS_gr = new TGraphErrors(i,x_Bj_RSS,g1_RSS,NULL,eg1_RSS);
  TGraphErrors *g2_RSS_gr = new TGraphErrors(i,x_Bj_RSS,g2_RSS,NULL,eg2_RSS);
  TGraphErrors *g1W_RSS_gr = new TGraphErrors(i,W_RSS,g1_RSS,NULL,eg1_RSS);
  TGraphErrors *g2W_RSS_gr = new TGraphErrors(i,W_RSS,g2_RSS,NULL,eg2_RSS);
  TGraphErrors *sigTT_RSS_gr = new TGraphErrors(i,W_RSS,sigTT_RSS,NULL,esigTT_RSS);
  TGraphErrors *sigLT_RSS_gr = new TGraphErrors(i,W_RSS,sigLT_RSS,NULL,esigLT_RSS);


  A1_CLAS_gr->SetMarkerStyle(22);
  A1_CLASh_gr->SetMarkerStyle(23);
  A1x_CLAS_gr->SetMarkerStyle(22);
  g1_CLAS_gr->SetMarkerStyle(22);
  g1W_CLAS_gr->SetMarkerStyle(22);
  g1_CLASh_gr->SetMarkerStyle(23);
  A1_RSS_gr->SetMarkerStyle(21);
  A2_RSS_gr->SetMarkerStyle(21);
  A1x_RSS_gr->SetMarkerStyle(21);
  A2x_RSS_gr->SetMarkerStyle(21);
  g1_RSS_gr->SetMarkerStyle(21);
  g2_RSS_gr->SetMarkerStyle(21);
  g1W_RSS_gr->SetMarkerStyle(21);
  g2W_RSS_gr->SetMarkerStyle(21);
  sigTT_RSS_gr->SetMarkerStyle(21);
  sigLT_RSS_gr->SetMarkerStyle(21);

  A1_CLAS_gr->SetLineColor(kOrange);
  A1_CLASh_gr->SetLineColor(kCyan);
  A1x_CLAS_gr->SetLineColor(kOrange);
  g1_CLAS_gr->SetLineColor(kOrange);
  g1W_CLAS_gr->SetLineColor(kOrange);
  g1_CLASh_gr->SetLineColor(kCyan);
  A1_RSS_gr->SetLineColor(kRed);
  A2_RSS_gr->SetLineColor(kRed);
  A1x_RSS_gr->SetLineColor(kRed);
  A2x_RSS_gr->SetLineColor(kRed);
  g1_RSS_gr->SetLineColor(kRed);
  g2_RSS_gr->SetLineColor(kRed);
  g1W_RSS_gr->SetLineColor(kRed);
  g2W_RSS_gr->SetLineColor(kRed);
  sigTT_RSS_gr->SetLineColor(kRed);
  sigLT_RSS_gr->SetLineColor(kRed);

  A1_CLAS_gr->SetMarkerColor(kOrange);
  A1_CLASh_gr->SetMarkerColor(kCyan);
  A1x_CLAS_gr->SetMarkerColor(kOrange);
  g1_CLAS_gr->SetMarkerColor(kOrange);
  g1W_CLAS_gr->SetMarkerColor(kOrange);
  g1_CLASh_gr->SetMarkerColor(kCyan);
  A1_RSS_gr->SetMarkerColor(kRed);
  A2_RSS_gr->SetMarkerColor(kRed);
  A1x_RSS_gr->SetMarkerColor(kRed);
  A2x_RSS_gr->SetMarkerColor(kRed);
  g1_RSS_gr->SetMarkerColor(kRed);
  g2_RSS_gr->SetMarkerColor(kRed);
  g1W_RSS_gr->SetMarkerColor(kRed);
  g2W_RSS_gr->SetMarkerColor(kRed);
  sigTT_RSS_gr->SetMarkerColor(kRed);
  sigLT_RSS_gr->SetMarkerColor(kRed);

  A1_CLASl_gr->SetMarkerStyle(23); A1_CLASl_gr->SetMarkerColor(kGreen); A1_CLASl_gr->SetLineColor(kGreen);
  A1x_CLASl_gr->SetMarkerStyle(23); A1x_CLASl_gr->SetMarkerColor(kGreen); A1x_CLASl_gr->SetLineColor(kGreen);
  g1_CLASl_gr->SetMarkerStyle(23); g1_CLASl_gr->SetMarkerColor(kGreen); g1_CLASl_gr->SetLineColor(kGreen);
  
  // MAID 2003

  ifstream MAIDfile;
  MAIDfile.open("maid_output/SANE1_MAID.txt");
  
  TString line;
  Int_t i = 0;
  const Int_t n = 19;
  Double_t W_MAID[n],sigTT_MAID[n],sigLT_MAID[n],A1_MAID[n],A2_MAID[n];

  while(!MAIDfile.eof()) {
    
    line.ReadLine(MAIDfile);

    if(!line.Contains("#") && !line.IsNull()) {
      W_MAID[i] = line_elem(line,0)/1.e3;
      sigTT_MAID[i] = line_elem(line,2);
      sigLT_MAID[i] = line_elem(line,3);
      A1_MAID[i] = line_elem(line,4);
      A2_MAID[i] = line_elem(line,5);
      i++;
    }
  }
  
  TGraph *sigTT_MAID_gr = new TGraph(n,W_MAID,sigTT_MAID);
  TGraph *sigLT_MAID_gr = new TGraph(n,W_MAID,sigLT_MAID);
  TGraph *A1_MAID_gr = new TGraph(n,W_MAID,A1_MAID);
  TGraph *A2_MAID_gr = new TGraph(n,W_MAID,A2_MAID);

  sigTT_MAID_gr->SetLineStyle(2);
  sigLT_MAID_gr->SetLineStyle(2);
  A1_MAID_gr->SetLineStyle(2);
  A2_MAID_gr->SetLineStyle(2);
  sigTT_MAID_gr->SetLineColor(kPink);
  sigLT_MAID_gr->SetLineColor(kPink);
  A1_MAID_gr->SetLineColor(kPink);
  A2_MAID_gr->SetLineColor(kPink);

  // MAID 2007

  ifstream MAID7file;
  MAID7file.open("maid_output/SANE1_MAID2007.txt");
  
  TString line;
  Int_t i = 0;
  const Int_t n = 19;
  Double_t W_MAID7[n],sigTT_MAID7[n],sigLT_MAID7[n],A1_MAID7[n],A2_MAID7[n];

  while(!MAID7file.eof()) {
    
    line.ReadLine(MAID7file);

    if(!line.Contains("#") && !line.IsNull()) {
      W_MAID7[i] = line_elem(line,0)/1.e3;
      sigTT_MAID7[i] = line_elem(line,2);
      sigLT_MAID7[i] = line_elem(line,3);
      A1_MAID7[i] = line_elem(line,4);
      A2_MAID7[i] = line_elem(line,5);
      i++;
    }
  }
  
  TGraph *sigTT_MAID7_gr = new TGraph(n,W_MAID7,sigTT_MAID7);
  TGraph *sigLT_MAID7_gr = new TGraph(n,W_MAID7,sigLT_MAID7);
  TGraph *A1_MAID7_gr = new TGraph(n,W_MAID7,A1_MAID7);
  TGraph *A2_MAID7_gr = new TGraph(n,W_MAID7,A2_MAID7);

  sigTT_MAID7_gr->SetLineStyle(9);
  sigLT_MAID7_gr->SetLineStyle(9);
  A1_MAID7_gr->SetLineStyle(9);
  A2_MAID7_gr->SetLineStyle(9);
  sigTT_MAID7_gr->SetLineColor(kGreen);
  sigLT_MAID7_gr->SetLineColor(kGreen);
  A1_MAID7_gr->SetLineColor(kGreen);
  A2_MAID7_gr->SetLineColor(kGreen);

  TCanvas *W_sigTT_c = new TCanvas("W_sigTT_c","W_sigTT");
  W_sigTT_c->SetGrid();
  W_sigTT_c->cd(1);sigTT_RSS_gr->Draw("ap");
  sigTT_RC->SetMarkerStyle(20); sigTT_RC->SetMarkerColor(kBlue); sigTT_RC->SetLineColor(kBlue);
  sigTT_RC->Draw("p");
  sigTT_MAID_gr->Draw();
  sigTT_MAID7_gr->Draw();
  legsigTT = new TLegend(0.6,0.15,0.9,0.35);
  legsigTT->AddEntry(sigTT_RC,"This work with r.c","lep");
  legsigTT->AddEntry(sigTT_RSS_gr,"RSS","lep");
  legsigTT->AddEntry(sigTT_MAID_gr,"MAID2003","l");
  legsigTT->AddEntry(sigTT_MAID7_gr,"MAID2007","l");
  legsigTT->Draw();
  W_sigTT_c->Print(DataSetName+"_sigTT.png");

  TCanvas *W_sigLT_c = new TCanvas("W_sigLT_c","W_sigLT");
  W_sigLT_c->SetGrid();
  W_sigLT_c->cd(1);sigLT_RSS_gr->Draw("ap");
  sigLT_RC->SetMarkerStyle(20); sigLT_RC->SetMarkerColor(kBlue); sigLT_RC->SetLineColor(kBlue);
  sigLT_RC->Draw("p");
  sigLT_MAID_gr->Draw();
  sigLT_MAID7_gr->Draw();
  legsigLT = new TLegend(0.6,0.15,0.9,0.35);
  legsigLT->AddEntry(sigLT_RC,"This work with r.c","lep");
  legsigLT->AddEntry(sigLT_RSS_gr,"RSS","lep");
  legsigLT->AddEntry(sigLT_MAID_gr,"MAID2003","l");
  legsigLT->AddEntry(sigLT_MAID7_gr,"MAID2007","l");
  legsigLT->Draw();
  W_sigLT_c->Print(DataSetName+"_sigLT.png");

  TCanvas *W_A1_c = new TCanvas("W_A1_c","W_A1");
  //W_A1_c->SetGrid();
  W_A1_c->cd();
  A1_RC->SetTitle(NULL);//("A_{1} along W");
  A1_RC->GetXaxis()->SetTitle("W [GeV]");
  A1_RC->GetXaxis()->SetLimits(1.0,2.4);
  A1_RC->SetMaximum(1.0); A1_RC->SetMinimum(-1.0);
  A1_RC->Draw("ac");
  //A1_RC->SetTitle(NULL);
  A1_RC->SetMarkerStyle(20); A1_RC->SetMarkerColor(kBlue); A1_RC->SetLineColor(kBlue);
  A1_RC->Draw("ap");//("ap");
  A1sys_RC->Draw("||");//("p3");
  A1_noRC->Draw("c");
  //A1_MAID_gr->Draw();
  A1_MAID7_gr->Draw();
  A1_RSS_gr->Draw("p");
  A1_CLAS_gr->Draw("p");
  //A1_CLASh_gr->Draw("p");
  //A1_CLASl_gr->Draw("p");
  hzl->Draw("same");
  legA1 = new TLegend(0.5,0.15,0.9,0.35);
  legA1->AddEntry(A1_noRC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2} without r.c.","l");
  legA1->AddEntry(A1_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");//1.86,1.31
  legA1->AddEntry(A1_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2}","lep");
  legA1->AddEntry(A1_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  //legA1->AddEntry(A1_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  //legA1->AddEntry(A1_CLASh_gr,"CLAS eg1b <Q^{2}> = 2.05 GeV^{2}","lep");
  //legA1->AddEntry(A1_MAID_gr,"MAID2003","l");
  legA1->AddEntry(A1_MAID7_gr,"MAID2007","l");
  legA1->Draw();
  W_A1_c->Print(DataSetName+"_A1.png");

  TCanvas *W_A2_c = new TCanvas("W_A2_c","W_A2");
  //W_A2_c->SetGrid();
  W_A2_c->cd();
  A2_RC->GetXaxis()->SetLimits(1.0,2.4);
  A2_RC->SetTitle(NULL);//("A_{2} along W");
  A2_RC->GetXaxis()->SetTitle("W [GeV]");
  A2_RC->SetMaximum(1.0); A2_RC->SetMinimum(-1.0);
  //A2_noRC->SetTitle(NULL);
  A2_RC->SetMarkerStyle(20); A2_RC->SetMarkerColor(kBlue); A2_RC->SetLineColor(kBlue);
  A2_RC->Draw("ap");//("ap");
  A2sys_RC->Draw("||");
  A2_noRC->Draw("c");
  //A2_noRC->Draw("c");
  //A2_MAID_gr->Draw();
  A2_MAID7_gr->Draw();
  A2_RSS_gr->Draw("p");
  hzl->Draw("same");
  legA2 = new TLegend(0.5,0.7,0.8,0.9);
  legA2->AddEntry(A2_noRC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2} without r.c.","l");
  legA2->AddEntry(A2_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");//1.86,1.31
  legA2->AddEntry(A2_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2}","lep");
  //legA2->AddEntry(A2_MAID_gr,"MAID2003","l");
  legA2->AddEntry(A2_MAID7_gr,"MAID2007","l");
  legA2->Draw();
  W_A2_c->Print(DataSetName+"_A2.png");

  TCanvas *x_A1_c = new TCanvas("x_A1_c","x_A1");
  //x_A1_c->SetGrid();
  x_A1_c->cd();
  A1x_RC->GetXaxis()->SetLimits(0.,1.);//(0.15,0.85);
  A1x_RC->SetTitle("A_{1} along x");
  A1x_RC->GetXaxis()->SetTitle("x");
  A1x_RC->GetYaxis()->SetTitle("A_{1}");
  A1x_RC->SetMaximum(1.); A1x_RC->SetMinimum(-1.);
  //A1x_RC->SetMarkerStyle(20); A1x_RC->SetMarkerColor(kBlue); A1x_RC->SetLineColor(kBlue);
  A1x_RC->SetMarkerStyle(24); A1x_RC->SetMarkerColor(kBlack); A1x_RC->SetLineColor(kBlack);
  A1x_RC->Draw("ap");
  //A1xsys_RC->Draw("||");
  A1x_RSS_gr->Draw("p");
  A1x_CLAS_gr->Draw("p");
  //A1x_CLASl_gr->Draw("p");
  legA1x = new TLegend(0.6,0.15,0.9,0.25);//(0.1,0.2,0.4,0.4);
  //legA1x->AddEntry(A1x_noRC,"This work without r.c","l");
  legA1x->AddEntry(A1x_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");
  legA1x->AddEntry(A1x_RSS_gr,"RSS <Q^{2}> = 1.31 GeV^{2}","lep");
  legA1x->AddEntry(A1x_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  //legA1x->AddEntry(A1x_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  legA1x->Draw();
  x_A1_c->Print(DataSetName+"_A1x.png");

  TCanvas *x_A2_c = new TCanvas("x_A2_c","x_A2");
  //x_A2_c->SetGrid();
  x_A2_c->cd();
  A2x_RC->GetXaxis()->SetLimits(0.,1.);//(0.15,0.85);
  A2x_RC->SetTitle("A_{2} along x");
  A2x_RC->GetXaxis()->SetTitle("x");
  A2x_RC->GetYaxis()->SetTitle("A_{2}");
  A2x_RC->SetMaximum(1.); A2x_RC->SetMinimum(-1.);
  A2x_RC->Draw("ap");
  //A2x_RC->SetMarkerStyle(20); A2x_RC->SetMarkerColor(kBlue); A2x_RC->SetLineColor(kBlue);
  A2x_RC->SetMarkerStyle(24); A2x_RC->SetMarkerColor(kBlack); A2x_RC->SetLineColor(kBlack);
  A2x_RC->Draw("p");
  //A2xsys_RC->Draw("||");
  hzl->Draw("same");
  A2x_RSS_gr->Draw("p");
  legA2x = new TLegend(0.2,0.15,0.5,0.25);//(0.5,0.7,0.8,0.9);
  //legA2x->AddEntry(A2x_noRC,"This work without r.c","l");
  legA2x->AddEntry(A2x_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");
  legA2x->AddEntry(A2x_RSS_gr,"RSS <Q^{2}> = 1.31 GeV^{2}","lep");
  legA2x->Draw();
  x_A2_c->Print(DataSetName+"_A2x.png");


  TCanvas *x_g1_c = new TCanvas("x_g1_c","x_g1");
  //x_g1_c->SetGrid();
  x_g1_c->cd();
  g1_RC->SetTitle(NULL);//("g_{1} along x");
  g1_RC->GetXaxis()->SetTitle("x");
  g1_RC->GetXaxis()->SetLimits(0.1,0.8);
  g1_RC->GetYaxis()->SetTitle("g_{1}");
  g1_RC->SetMaximum(0.4); g1_RC->SetMinimum(-0.1);
  g1_RC->SetMarkerStyle(24); g1_RC->SetMarkerColor(kBlack); g1_RC->SetLineColor(kBlack);
  g1_RC->SetMarkerStyle(20); g1_RC->SetMarkerColor(kBlue); g1_RC->SetLineColor(kBlue);
  g1_RC->Draw("ap");
  g1sys_RC->Draw("||");
  g1_CLAS_gr->Draw("p");
  g1_CLASl_gr->Draw("p");
  g1_RSS_gr->Draw("p");
  hzl->Draw("same");
  legg1 = new TLegend(0.6,0.7,0.9,0.85);//(0.5,0.6,0.85,0.8);
  legg1->AddEntry(g1_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");//1.86,1.31
  legg1->AddEntry(g1_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2} ","lep");
  //legg1->AddEntry(g1_CLASh_gr,"CLAS eg1b <Q^{2}> = 2.05 GeV^{2}","lep");
  legg1->AddEntry(g1_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  legg1->AddEntry(g1_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2} ","lep");
  legg1->Draw();
  x_g1_c->Print(DataSetName+"_g1.png");

  TCanvas *x_g2_c = new TCanvas("x_g2_c","x_g2");
  //x_g2_c->SetGrid();
  x_g2_c->cd();
  g2_RC->SetTitle(NULL);//("g_{2} along x");
  g2_RC->GetYaxis()->SetTitle("g_{2}");
  g2_RC->GetXaxis()->SetTitle("x");
  g2_RC->GetXaxis()->SetLimits(0.1,0.8);//(0.15,0.85);
  g2_RC->SetMaximum(0.1); g2_RC->SetMinimum(-0.5);//(-0.15);
  g2_RC->SetMarkerStyle(20); g2_RC->SetMarkerColor(kBlue); g2_RC->SetLineColor(kBlue);
  g2_RC->Draw("ap");
  g2sys_RC->Draw("||");
  g2WW_RC->SetLineStyle(3); g2WW_RC->SetLineColor(kMagenta);
  g2WW_RC->Draw("c");
  g2_RSS_gr->Draw("p");
  hzl->Draw("same");
  legg2 = new TLegend(0.6,0.15,0.9,0.25);//(0.5,0.15,0.85,0.35);
  legg2->AddEntry(g2_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");//1.86,1.31
  legg2->AddEntry(g2WW_RC,"Approximation g_{2}^{WW} of SANE-HMS","l");
  legg2->AddEntry(g2_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2} ","lep");
  legg2->Draw();
  x_g2_c->Print(DataSetName+"_g2.png");

  //Additional investigation
  TCanvas *W_A1_c2 = new TCanvas("W_A1_c2","W_A1");
  W_A1_c2->SetGrid();
  W_A1_c2->cd();
  A1_noRC->SetTitle("A_{1} along W");
  A1_noRC->GetXaxis()->SetTitle("W [GeV]");
  A1_noRC->SetMaximum(1.0); A1_noRC->SetMinimum(-1.0);
  //A1_noRC->Draw("ac");
  A1_RC->SetMarkerStyle(20); A1_RC->SetMarkerColor(kBlue); A1_RC->SetLineColor(kBlue);
  A1_RC->Draw("p");
  A1_RSS_gr->Draw("p");
  A1_CLAS_gr->SetMarkerStyle(21); A1_CLAS_gr->SetMarkerColor(kYellow); A1_CLAS_gr->SetLineColor(kYellow);
  A1_CLAS_gr->Draw("p");
  A1_CLASh_gr->SetMarkerStyle(22); A1_CLASh_gr->SetMarkerColor(kMagenta); A1_CLASh_gr->SetLineColor(kMagenta);
  A1_CLASh_gr->Draw("p");
  A1_CLASl_gr->SetMarkerStyle(23); A1_CLASl_gr->SetMarkerColor(kGreen); A1_CLASl_gr->SetLineColor(kGreen);
  A1_CLASl_gr->Draw("p");
  A1_CLASll_gr->SetMarkerStyle(24); A1_CLASll_gr->SetMarkerColor(kCyan); A1_CLASll_gr->SetLineColor(kCyan);
  A1_CLASll_gr->Draw("p");
  legA1 = new TLegend(0.5,0.15,0.9,0.35);
  legA1->AddEntry(A1_noRC,"This work without r.c","l");
  legA1->AddEntry(A1_RC,"This work with r.c","lep");
  legA1->AddEntry(A1_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2}","lep");
  legA1->AddEntry(A1_CLASll_gr,"CLAS eg1b <Q^{2}> = 1.20 GeV^{2}","lep");
  legA1->AddEntry(A1_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  legA1->AddEntry(A1_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  legA1->AddEntry(A1_CLASh_gr,"CLAS eg1b <Q^{2}> = 2.05 GeV^{2}","lep");
  legA1->Draw();
  //W_A1_c2->Print(DataSetName+"_A1_CLAS.png");

  TCanvas *x_A1_c2 = new TCanvas("x_A1_c2","x_A1");
  x_A1_c2->SetGrid();
  x_A1_c2->cd();
  A1x_noRC->SetTitle("A_{1} along x");
  A1x_noRC->GetXaxis()->SetTitle("x");
  A1x_noRC->GetXaxis()->SetLimits(0.2,0.9);
  A1x_noRC->SetMaximum(1.0); A1x_noRC->SetMinimum(-1.0);
  A1x_noRC->Draw("ac");
  A1x_RC->SetMarkerStyle(20); A1x_RC->SetMarkerColor(kBlue); A1x_RC->SetLineColor(kBlue);
  A1x_RC->Draw("p");
  A1x_RSS_gr->Draw("p");
  A1x_CLAS_gr->SetMarkerStyle(21); A1x_CLAS_gr->SetMarkerColor(kYellow); A1x_CLAS_gr->SetLineColor(kYellow);
  A1x_CLAS_gr->Draw("p");
  A1x_CLASh_gr->SetMarkerStyle(22); A1x_CLASh_gr->SetMarkerColor(kMagenta); A1x_CLASh_gr->SetLineColor(kMagenta);
  A1x_CLASh_gr->Draw("p");
  A1x_CLASl_gr->SetMarkerStyle(23); A1x_CLASl_gr->SetMarkerColor(kGreen); A1x_CLASl_gr->SetLineColor(kGreen);
  A1x_CLASl_gr->Draw("p");
  A1x_CLASll_gr->SetMarkerStyle(24); A1x_CLASll_gr->SetMarkerColor(kCyan); A1x_CLASll_gr->SetLineColor(kCyan);
  A1x_CLASll_gr->Draw("p");
  legA1x = new TLegend(0.1,0.2,0.4,0.4);
  legA1x->AddEntry(A1x_noRC,"This work without r.c","l");
  legA1x->AddEntry(A1x_RC,"This work with r.c","lep");
  legA1x->AddEntry(A1x_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2}","lep");
  legA1x->AddEntry(A1x_CLASll_gr,"CLAS eg1b <Q^{2}> = 1.20 GeV^{2}","lep");
  legA1x->AddEntry(A1x_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  legA1x->AddEntry(A1x_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  legA1x->AddEntry(A1x_CLASh_gr,"CLAS eg1b <Q^{2}> = 2.05 GeV^{2}","lep");
  legA1x->Draw();
  //x_A1_c2->Print(DataSetName+"_A1x_CLAS.png");

  TCanvas *x_g1_c2 = new TCanvas("x_g1_c2","x_g1");
  x_g1_c2->SetGrid();
  x_g1_c2->cd();
  g1_RC->SetTitle("g_{1} along x");
  g1_RC->GetXaxis()->SetTitle("x");
  g1_RC->SetMaximum(0.2); g1_RC->SetMinimum(-0.05);
  g1_RC->SetMarkerStyle(20); g1_RC->SetMarkerColor(kBlue); g1_RC->SetLineColor(kBlue);
  g1_RC->Draw("ap");
  g1_CLAS_gr->SetMarkerStyle(21); g1_CLAS_gr->SetMarkerColor(kBlack); g1_CLAS_gr->SetLineColor(kBlack);
  g1_CLAS_gr->Draw("p");
  g1_CLASh_gr->SetMarkerStyle(22); g1_CLASh_gr->SetMarkerColor(kMagenta); g1_CLASh_gr->SetLineColor(kMagenta);
  //g1_CLASh_gr->Draw("p");
  g1_CLASl_gr->SetMarkerStyle(23); g1_CLASl_gr->SetMarkerColor(kGreen); g1_CLASl_gr->SetLineColor(kGreen);
  //g1_CLASl_gr->Draw("p");
  legg1 = new TLegend(0.5,0.6,0.85,0.85);
  legg1->AddEntry(g1_RC,"This work with r.c","lep");
  legg1->AddEntry(g1_CLASl_gr,"CLAS eg1b <Q^{2}> = 1.44 GeV^{2}","lep");
  legg1->AddEntry(g1_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  legg1->AddEntry(g1_CLASh_gr,"CLAS eg1b <Q^{2}> = 2.05 GeV^{2}","lep");
  legg1->Draw();
//  x_g1_c2->Print(DataSetName+"_g1_CLAS.png");

  TCanvas *W_g1_c = new TCanvas("W_g1_c","W_g1");
//  W_g1_c->SetGrid();
  W_g1_c->cd();
  g1W_RC->SetTitle("g_{1} along W");
  g1W_RC->GetXaxis()->SetTitle("W (GeV)");
  g1W_RC->GetXaxis()->SetLimits(1.0,2.4);
  g1W_RC->SetMaximum(0.4); g1W_RC->SetMinimum(-0.05);
  g1W_RC->SetMarkerStyle(20); g1W_RC->SetMarkerColor(kBlue); g1W_RC->SetLineColor(kBlue);
  g1W_RC->Draw("ap");
  g1W_RSS_gr->Draw("p");
  g1W_CLAS_gr->Draw("p");
  hzl->Draw("same");
  legg1W = new TLegend(0.1,0.75,0.5,0.85);
  legg1W->AddEntry(g1W_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");
  legg1W->AddEntry(g1W_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2} ","lep");
  legg1W->AddEntry(g1W_CLAS_gr,"CLAS eg1b <Q^{2}> = "+CLASQ2Name+" GeV^{2}","lep");
  legg1W->Draw();
  W_g1_c->Print(DataSetName+"_g1W.png");
  
  TCanvas *W_g2_c = new TCanvas("W_g2_c2","W_g2");
//  W_g2_c->SetGrid();
  W_g2_c->cd();
  g2W_RC->SetTitle("g_{2} along W");
  g2W_RC->GetXaxis()->SetTitle("W (GeV)");
  g2W_RC->GetXaxis()->SetLimits(1.0,2.4);
  g2W_RC->SetMaximum(0.05); g2W_RC->SetMinimum(-0.4);
  g2W_RC->SetMarkerStyle(20); g2W_RC->SetMarkerColor(kBlue); g2W_RC->SetLineColor(kBlue);
  g2W_RC->Draw("ap");
  g2W_RSS_gr->Draw("p");
  hzl->Draw("same");
  legg2W = new TLegend(0.1,0.15,0.5,0.25);
  legg2W->AddEntry(g2W_RC,"SANE-HMS <Q^{2}> = "+DataQ2Name+" GeV^{2}","lep");
  legg2W->AddEntry(g2W_RSS_gr,"RSS <Q^{2}> = 1.3 GeV^{2} ","lep");
  legg2W->Draw();
  W_g2_c->Print(DataSetName+"_g2W.png");

  gStyle->SetOptFit(kTRUE);
  TCanvas *A1fit_c = new TCanvas("A1fit_c","A1fit");
  A1fit_c->cd();
  //A1fit_c->SetTheta(0);
  //A1fit_c->SetPhi(0);
  //A1_WQ2->Draw("err p0");
  A1_WQ2->SetTitle("");
  A1_WQ2->GetXaxis()->SetTitle("W [GeV]");
  A1_WQ2->GetYaxis()->SetTitle("Q^2 [GeV]");
  TF2 *f_A1 = new TF2("f_A1","[0]+[1]/x+[2]*y");
  A1_WQ2->Fit(f_A1);
  TF2 *fit_A1 = (TF2*)A1_WQ2->FindObject("f_A1");
  A1_WQ2->Draw("err p0");
  fit_A1->Draw("same surf");
  A1fit_c->Print(DataSetName+"_A1fit.png");

  TCanvas *A2fit_c = new TCanvas("A2fit_c","A2fit");
  A2fit_c->cd();
  //A2fit_c->SetTheta(0);
  //A2fit_c->SetPhi(0);
  //A1_WQ2->Draw("err p0");
  A2_WQ2->SetTitle("");
  A2_WQ2->GetXaxis()->SetTitle("W [GeV]");
  A2_WQ2->GetYaxis()->SetTitle("Q^2 [GeV]");
  TF2 *f_A2 = new TF2("f_A2","[0]+[1]/x+[2]/sqrt(y)");
  A2_WQ2->Fit(f_A2);
  TF2 *fit_A2 = (TF2*)A2_WQ2->FindObject("f_A2");
  A2_WQ2->Draw("err p0");
  fit_A2->Draw("same surf");
  A2fit_c->Print(DataSetName+"_A2fit.png");

  //TFile *result = new TFile("histgrams.root","RECREATE");
  //Hlist.Write();
  //result->Close();
}

Double_t line_elem(TString line, Int_t index) {
  TObjArray *sublist = line.Tokenize("\t: ");
  TObjString *sub = sublist->At(index);
  return atof(sub->GetString().Data());
}