Index: trunk/examples/heraCompare.cpp
===================================================================
--- trunk/examples/heraCompare.cpp (revision 332)
+++ trunk/examples/heraCompare.cpp (revision 333)
@@ -1,1694 +1,1702 @@
//==============================================================================
// heraCompare.cpp
//
-// Copyright (C) 2010-2013 Tobias Toll and Thomas Ullrich
+// Copyright (C) 2010-2018 Tobias Toll and Thomas Ullrich
//
// This file is part of Sartre.
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation.
// This program is distributed in the hope that it will be useful,
// but without any warranty; without even the implied warranty of
// merchantability or fitness for a particular purpose. See the
// GNU General Public License for more details.
// You should have received a copy of the GNU General Public License
// along with this program. If not, see .
//
// Author: Thomas Ullrich and Tobias Toll
// Last update:
// $Date$
// $Author$
//==============================================================================
//
// Main program to compare Sartre with HERA data (J/psi and DVCS)
// Usage:
// heraCompare rootfile jpsi|phi|dvcs
//
// Output are graphs with HERA and Sartre data to compare.
//==============================================================================
#include
#include
#include
#include
#include "Sartre.h"
#include "TFile.h"
#include "TGraphAsymmErrors.h"
#include "TGraphErrors.h"
#include "Constants.h"
using namespace std;
#define PR(x) cout << #x << " = " << (x) << endl;
double betaMin=0.7;
double xpomMax=1e-2;
+double maxLam=0.35;
+DipoleModelType theDipoleModelType = bNonSat;
+DipoleModelParameterSet theDipoleModelParameterSet = HMPZ;
+TableSetType theTableSetType = total_and_coherent; // or coherent_and_incoherent;
-double maxLam=0.8;
struct HeraDataPoint {
double Q2;
double W;
double t;
double tmin;
double tmax;
double Q2min;
double Q2max;
double Wmin;
double Wmax;
double xsection;
double xsection_stat;
double xsection_sys;
double xsection_sys_up;
double xsection_sys_lo;
HeraDataPoint();
};
HeraDataPoint::HeraDataPoint()
{
t=-1;
tmin=-1; tmax=-1;
Q2=-1;
Q2min=-1; Q2max=-1;
W=-1;
Wmin=-1; Wmax=-1;
xsection=-1;
xsection_stat=-1;
xsection_sys=-1;
xsection_sys_up=-1; xsection_sys_lo=-1;
}
void loadJpsiData(vector&);
void loadPhiData(vector&);
void loadDVCSData(vector&, vector&);
void compareJpsi(string&);
void compareDVCS(string&);
void comparePhi(string&);
double chi2_global = 0;
int ndof_global = 0;
int main(int argc, char** argv)
{
//
// Check command line arguments
//
if (argc != 3) {
cout << "Usage: " << argv[0] << " rootfile jpsi|phi|dvcs|all" << endl;
return 2;
}
//
// Run the requested comparison(s)
//
string rootfile = argv[1];
string what = argv[2];
if (what.find("jpsi") != string::npos)
compareJpsi(rootfile);
else if (what.find("dvcs") != string::npos)
compareDVCS(rootfile);
else if (what.find("phi") != string::npos)
comparePhi(rootfile);
else if (what.find("all") != string::npos){
compareDVCS(rootfile);
comparePhi(rootfile);
compareJpsi(rootfile);
cout << "-------------------------------------------------" << endl;
cout << "Global Comparison HERA-Sartre:" << endl;
cout << "chi2 = " << chi2_global << endl;
cout << "ndof = " << ndof_global << endl;
cout << "chi2/ndof = " << chi2_global/ndof_global << endl;
cout << "-------------------------------------------------" << endl;
}
else {
cout << "Error: '" << what.c_str() << "' no such probe to compare." << endl;
return 2;
}
return 0;
}
void comparePhi(string& rootfile)
{
//
// PHI HERA/ZEUS data
//
double heraProtonEnergy = 920;
double heraElectronEnergy = 27.6;
vector zeusDataPhi_table1;
loadPhiData(zeusDataPhi_table1);
vector graphVec;
vector sgraphVec;
vector sgraphName;
vector graphName;
// Table 1
TGraphErrors *graph1a = new TGraphErrors(7);
graph1a->SetTitle("Q2=2.4");
graphVec.push_back(graph1a);
graphName.push_back("graph1a");
TGraphErrors *graph1b = new TGraphErrors(6);
graph1b->SetTitle("Q2=3.8");
graphVec.push_back(graph1b);
graphName.push_back("graph1b");
TGraphErrors *graph1c = new TGraphErrors(6);
graph1c->SetTitle("Q2=6.5");
graphVec.push_back(graph1c);
graphName.push_back("graph1c");
TGraphErrors *graph1d = new TGraphErrors(6);
graph1d->SetTitle("Q2=13.0");
graphVec.push_back(graph1d);
graphName.push_back("graph1d");
TGraphErrors *sgraph1a = new TGraphErrors(7);
sgraph1a->SetTitle("Q2=2.4");
sgraphVec.push_back(sgraph1a);
sgraphName.push_back("sgraph1a");
TGraphErrors *sgraph1b = new TGraphErrors(6);
sgraph1b->SetTitle("Q2=3.8");
sgraphVec.push_back(sgraph1b);
sgraphName.push_back("sgraph1b");
TGraphErrors *sgraph1c = new TGraphErrors(6);
sgraph1c->SetTitle("Q2=6.5");
sgraphVec.push_back(sgraph1c);
sgraphName.push_back("sgraph1c");
TGraphErrors *sgraph1d = new TGraphErrors(6);
sgraph1d->SetTitle("Q2=13.0");
sgraphVec.push_back(sgraph1d);
sgraphName.push_back("sgraph1d");
// Table 3
TGraphErrors *graph3 = new TGraphErrors(9);
graph3->SetTitle("Q2");
graphVec.push_back(graph3);
graphName.push_back("graph3");
TGraphErrors *sgraph3 = new TGraphErrors(9);
sgraph3->SetTitle("Q2");
sgraphVec.push_back(sgraph3);
sgraphName.push_back("sgraph3");
//
// Loop over data points and for each initialize Sartre and
// calculate the cross-section. We are using no runcard but
// set thigs for each data point
//
Sartre sartre;
EventGeneratorSettings* settings = sartre.runSettings();
//
// Table 1
//
double diff_table1, chi2_table1 = 0;
double diff_table3, chi2_table3 = 0;
int ndof_table1 = 0;
int ndof_table3 = 0;
int graphCounter = 0;
int pointCounter = 0;
for (unsigned int i=0; iSetPoint(pointCounter, zeusDataPhi_table1[i].W, zeusDataPhi_table1[i].xsection);
if(graphCounter==4)
graphVec[graphCounter]->SetPoint(pointCounter, zeusDataPhi_table1[i].Q2, zeusDataPhi_table1[i].xsection);
double error = sqrt(zeusDataPhi_table1[i].xsection_stat*zeusDataPhi_table1[i].xsection_stat + zeusDataPhi_table1[i].xsection_sys*zeusDataPhi_table1[i].xsection_sys);
graphVec[graphCounter]->SetPointError(pointCounter, 0., error);
// set up sartre
- settings->setVerbose(false);
- settings->setVerboseLevel(0);
+ settings->setVerbose(true);
+ settings->setVerboseLevel(2);
settings->setNumberOfEvents(0);
settings->setTimesToShow(0);
settings->setVectorMesonId(333);
settings->setApplyPhotonFlux(false);
settings->setElectronBeamEnergy(heraElectronEnergy);
settings->setHadronBeamEnergy(heraProtonEnergy);
- settings->setDipoleModelType(bSat);
- settings->setA(1);
- settings->setCorrectForRealAmplitude(false);
- settings->setCorrectSkewedness(false);
- settings->setEnableNuclearBreakup(false);
+ settings->setDipoleModelType(theDipoleModelType);
+ settings->setDipoleModelParameterSet(theDipoleModelParameterSet);
+ settings->setA(1);
+ settings->setCorrectForRealAmplitude(true);
+ settings->setCorrectSkewedness(true);
+ settings->setMaxLambdaUsedInCorrections(maxLam);
+ settings->setEnableNuclearBreakup(false);
settings->setQ2min(zeusDataPhi_table1[i].Q2min); double Q2min=settings->Q2min();
settings->setQ2max(zeusDataPhi_table1[i].Q2max); double Q2max=settings->Q2max();
settings->setWmin(zeusDataPhi_table1[i].Wmin); double Wmin=settings->Wmin();
settings->setWmax(zeusDataPhi_table1[i].Wmax); double Wmax=settings->Wmax();
- settings->setMaxLambdaUsedInCorrections(maxLam);
if(Wmax > 142 || Wmin > 142) continue;
//calculate in the t=0 limits, mp=0:
TParticlePDG *vectorMesonPDG = settings->lookupPDG(settings->vectorMesonId());
double vmMass = vectorMesonPDG->Mass();
double beta=Q2min/(vmMass*vmMass + Q2min);
double xpom=(vmMass*vmMass + Q2max) / (Wmin*Wmin + Q2max - protonMass2);
if(xpom > xpomMax) continue;
if(beta < betaMin) continue;
bool ok = sartre.init();
if (!ok) {
cout << "Initialization of sartre failed. Going to next data point." << endl;
continue;
}
- // settings->list();
+ settings->list();
//
// Sartre cross-section
//
vector > limits = sartre.kinematicLimits(); // t, Q2, W
double lower[3], upper[3];
lower[0] = -zeusDataPhi_table1[i].tmax; // t
upper[0] = -zeusDataPhi_table1[i].tmin;
lower[1] = limits[1].first; // Q2
upper[1] = limits[1].second;
lower[2] = limits[2].first; // W
upper[2] = limits[2].second;
double CS = sartre.totalCrossSection(lower, upper);
CS /= upper[0]-lower[0];
CS /= limits[1].second - limits[1].first;
CS /= limits[2].second*limits[2].second - limits[2].first*limits[2].first;
if(graphCounter<4){
sgraphVec[graphCounter]->SetPoint(pointCounter, zeusDataPhi_table1[i].W, CS);
sgraphVec[graphCounter]->SetPointError(pointCounter, 0., 0.);
}
if(graphCounter==4){
sgraphVec[graphCounter]->SetPoint(pointCounter, zeusDataPhi_table1[i].Q2, CS);
sgraphVec[graphCounter]->SetPointError(pointCounter, 0., 0.);
}
//
// chi2 calculation
//
if(graphCounter<4){
diff_table1 = CS - zeusDataPhi_table1[i].xsection;
chi2_table1 += (diff_table1*diff_table1) / (error*error);
ndof_table1++;
chi2_global+=(diff_table1*diff_table1) / (error*error);
ndof_global++;
}
if(graphCounter==4){
diff_table3 = CS - zeusDataPhi_table1[i].xsection;
chi2_table3 += (diff_table3*diff_table3) / (error*error);
ndof_table3++;
chi2_global+=(diff_table3*diff_table3) / (error*error);
ndof_global++;
}
//
// Print out
//
cout << endl;
cout << "dsig/dt at |t|=" << fabs(zeusDataPhi_table1[i].t)
<< ", Q2=" << zeusDataPhi_table1[i].Q2
<< ", W=" << zeusDataPhi_table1[i].W << endl;
cout << "\t\tSartre: " << CS << endl; // in nb/GeV2
cout << "\t\tZeus: " << zeusDataPhi_table1[i].xsection << "+/-" << zeusDataPhi_table1[i].xsection_stat << "+/-" <Write(graphName[i].c_str());
sgraphVec[i]->Write(sgraphName[i].c_str());
}
// Table 1
// graph1a->Write("graph1a");
// sgraph1a->Write("sgraph1a");
// for(int i=0; i<3; i++){
// graphVec.at(i)->Write("graph1a");
// sgraphVec.at(i)->Write("sgraph1a");
// }
hfile->Close();
}
void compareDVCS(string& rootfile)
{
//
// DVCS HERA/H1 data
//
double heraProtonEnergy = 920;
double heraElectronEnergy = 27.6;
vector h1DataDVCS_table4;
vector h1DataDVCS_table1;
loadDVCSData(h1DataDVCS_table1, h1DataDVCS_table4);
vector graphVec;
vector sgraphVec;
vector sgraphName;
vector graphName;
// Table 4
TGraphErrors *graph4a = new TGraphErrors(4); graph4a->SetTitle("Q2=8, W=40"); graphVec.push_back(graph4a); graphName.push_back("graph4a");
TGraphErrors *graph4b = new TGraphErrors(4); graph4b->SetTitle("Q2=8, W=70"); graphVec.push_back(graph4b); graphName.push_back("graph4b");
TGraphErrors *graph4c = new TGraphErrors(4); graph4c->SetTitle("Q2=8, W=100"); graphVec.push_back(graph4c); graphName.push_back("graph4c");
TGraphErrors *graph4d = new TGraphErrors(4); graph4d->SetTitle("Q2=20, W=40"); graphVec.push_back(graph4d); graphName.push_back("graph4d");
TGraphErrors *graph4e = new TGraphErrors(4); graph4e->SetTitle("Q2=20, W=70"); graphVec.push_back(graph4e); graphName.push_back("graph4e");
TGraphErrors *graph4f = new TGraphErrors(4); graph4f->SetTitle("Q2=20, W=100"); graphVec.push_back(graph4f); graphName.push_back("graph4f");
TGraphErrors *sgraph4a = new TGraphErrors(4); sgraph4a->SetTitle("Q2=8, W=40"); sgraphVec.push_back(sgraph4a); sgraphName.push_back("sgraph4a");
TGraphErrors *sgraph4b = new TGraphErrors(4); sgraph4b->SetTitle("Q2=8, W=70"); sgraphVec.push_back(sgraph4b); sgraphName.push_back("sgraph4b");
TGraphErrors *sgraph4c = new TGraphErrors(4); sgraph4c->SetTitle("Q2=8, W=100"); sgraphVec.push_back(sgraph4c); sgraphName.push_back("sgraph4c");
TGraphErrors *sgraph4d = new TGraphErrors(4); sgraph4d->SetTitle("Q2=20, W=40"); sgraphVec.push_back(sgraph4d); sgraphName.push_back("sgraph4d");
TGraphErrors *sgraph4e = new TGraphErrors(4); sgraph4e->SetTitle("Q2=20, W=70"); sgraphVec.push_back(sgraph4e); sgraphName.push_back("sgraph4e");
TGraphErrors *sgraph4f = new TGraphErrors(4); sgraph4f->SetTitle("Q2=20, W=100"); sgraphVec.push_back(sgraph4f); sgraphName.push_back("sgraph4f");
// Table 1
TGraphErrors *graph1a = new TGraphErrors(4); graph1a->SetTitle("30 < W < 140, |t|<1");
TGraphErrors *sgraph1a = new TGraphErrors(4); sgraph1a->SetTitle("30 < W < 140, |t|<1");
TGraphErrors *graph1b = new TGraphErrors(5); graph1b->SetTitle("6.5 < Q2 < 80, |t|<1");
TGraphErrors *sgraph1b = new TGraphErrors(5); sgraph1b->SetTitle("6.5 < Q2 < 80, |t|<1");
//
// Loop over data points and for each initialize Sartre and
// calculate the cross-section. We are using no runcard but
// set thigs for each data point
//
Sartre sartre;
EventGeneratorSettings* settings = sartre.runSettings();
//
// Table 4
//
double diff_table4, chi2_table4 = 0;
int ndof_table4 = 0;
int graphCounter, pointCounter;
for (unsigned int i=0; iSetPoint(pointCounter, fabs(h1DataDVCS_table4[i].t), h1DataDVCS_table4[i].xsection);
double error = sqrt(h1DataDVCS_table4[i].xsection_stat*h1DataDVCS_table4[i].xsection_stat + h1DataDVCS_table4[i].xsection_sys*h1DataDVCS_table4[i].xsection_sys);
graphVec[graphCounter]->SetPointError(pointCounter, 0., error);
// set up sartre
settings->setVerbose(false);
settings->setVerboseLevel(0);
settings->setNumberOfEvents(0);
settings->setTimesToShow(0);
settings->setVectorMesonId(22);
settings->setApplyPhotonFlux(false);
settings->setElectronBeamEnergy(heraElectronEnergy);
settings->setHadronBeamEnergy(heraProtonEnergy);
- settings->setDipoleModelType(bNonSat);
+ settings->setDipoleModelType(theDipoleModelType);
+ settings->setDipoleModelParameterSet(theDipoleModelParameterSet);
settings->setA(1);
settings->setCorrectForRealAmplitude(false);
settings->setCorrectSkewedness(false);
settings->setEnableNuclearBreakup(false);
settings->setQ2min(h1DataDVCS_table4[i].Q2-.1);
settings->setQ2max(h1DataDVCS_table4[i].Q2+.1);
settings->setWmin(h1DataDVCS_table4[i].W-.1);
settings->setWmax(h1DataDVCS_table4[i].W+.1);
if (settings->Wmax() > 142 || settings->Wmin()>142) continue;
settings->setMaxLambdaUsedInCorrections(maxLam);
bool ok = sartre.init();
if (!ok) {
cout << "Initialization of sartre failed. Going to next data point." << endl;
continue;
}
//
// Sartre cross-section
//
vector > limits = sartre.kinematicLimits(); // t, Q2, W
double lower[3], upper[3];
lower[0] = -h1DataDVCS_table4[i].t-0.05; // t
upper[0] = -h1DataDVCS_table4[i].t+0.05;
lower[1] = limits[1].first; // Q2
upper[1] = limits[1].second;
lower[2] = limits[2].first; // W
upper[2] = limits[2].second;
double CS = sartre.totalCrossSection(lower, upper);
CS /= upper[0]-lower[0];
CS /= limits[1].second - limits[1].first;
CS /= limits[2].second*limits[2].second - limits[2].first*limits[2].first;
sgraphVec[graphCounter]->SetPoint(pointCounter, fabs(h1DataDVCS_table4[i].t), CS);
sgraphVec[graphCounter]->SetPointError(pointCounter, 0., 0.);
//
// chi2 calculation
//
diff_table4 = CS - h1DataDVCS_table4[i].xsection;
chi2_table4 += (diff_table4*diff_table4) / (error*error);
ndof_table4++;
//
// Print out
//
cout << endl;
cout << "dsig/dt at |t|=" << fabs(h1DataDVCS_table4[i].t)
<< ", Q2=" << h1DataDVCS_table4[i].Q2
<< ", W=" << h1DataDVCS_table4[i].W << endl;
cout << "\t\tSartre: " << CS << endl; // in nb/GeV2
cout << "\t\tH1: " << h1DataDVCS_table4[i].xsection << "+/-" << h1DataDVCS_table4[i].xsection_stat << "+/-" <SetPoint(i, fabs(h1DataDVCS_table1[i].Q2), h1DataDVCS_table1[i].xsection);
graph1a->SetPointError(i , 0., error);
}
else {
graph1b->SetPoint(i-4, fabs(h1DataDVCS_table1[i].W), h1DataDVCS_table1[i].xsection);
graph1b->SetPointError(i-4 , 0., error);
}
// set up sartre
settings->setVerbose(true);
settings->setVerboseLevel(0);
settings->setNumberOfEvents(0);
settings->setTimesToShow(0);
settings->setVectorMesonId(22);
settings->setApplyPhotonFlux(false);
settings->setElectronBeamEnergy(heraElectronEnergy);
settings->setHadronBeamEnergy(heraProtonEnergy);
- settings->setDipoleModelType(bNonSat);
+ settings->setDipoleModelType(theDipoleModelType);
+ settings->setDipoleModelParameterSet(theDipoleModelParameterSet);
settings->setA(1);
settings->setCorrectForRealAmplitude(false);
settings->setCorrectSkewedness(false);
settings->setEnableNuclearBreakup(false);
settings->setQ2min(h1DataDVCS_table1[i].Q2-.1);
settings->setQ2max(h1DataDVCS_table1[i].Q2+.1);
settings->setWmin(h1DataDVCS_table1[i].W-.1);
settings->setWmax(h1DataDVCS_table1[i].W+.1);
settings->setMaxLambdaUsedInCorrections(maxLam);
if(settings->Wmax() > 142 || settings->Wmin()>142)
continue;
bool ok = sartre.init();
if (!ok) {
cout << "Initialization of sartre failed. Going to next data point." << endl;
continue;
}
//
// Sartre cross-section
//
vector > limits = sartre.kinematicLimits(); // t, Q2, W
double lower[3], upper[3];
lower[0] = -1; // t
upper[0] = limits[0].second;
lower[1] = limits[1].first; // Q2
upper[1] = limits[1].second;
lower[2] = limits[2].first; // W
upper[2] = limits[2].second;
double CS = sartre.totalCrossSection(lower, upper);
// H1 cross-section sum for crosschecks
double h1cs = h1DataDVCS_table1[i].xsection;
h1cs *= h1DataDVCS_table1[i].Q2max - h1DataDVCS_table1[i].Q2min; // Q2
h1cs *= h1DataDVCS_table1[i].Wmax - h1DataDVCS_table1[i].Wmin; // W
if (i < 4) {
sigma_a += CS;
h1_sigma_a += h1cs;
}
else {
sigma_b += CS;
h1_sigma_b += h1cs;
}
limits = sartre.kinematicLimits();
CS /= limits[1].second - limits[1].first; // Q2
CS /= limits[2].second*limits[2].second - limits[2].first*limits[2].first; // W2
// fill graphs
if (i < 4) {
sgraph1a->SetPoint(i, fabs(h1DataDVCS_table1[i].Q2), CS);
sgraph1a->SetPointError(i , 0., 0.);
}
else {
sgraph1b->SetPoint(i-4, fabs(h1DataDVCS_table1[i].W), CS);
sgraph1b->SetPointError(i-4 , 0., 0.);
}
//
// chi2 calculation
//
diff_table1 = CS - h1DataDVCS_table1[i].xsection;
chi2_table1 += (diff_table1*diff_table1) / (error*error);
chi2_global += (diff_table1*diff_table1) / (error*error);
ndof_table1++;
ndof_global++;
//
// Print out
//
cout << endl;
cout << "\t\tSartre: " << CS << endl;
cout << "\t\tH1: " << h1DataDVCS_table1[i].xsection << "+/-" << h1DataDVCS_table1[i].xsection_stat << "+/-" <Write(graphName[i].c_str());
sgraphVec[i]->Write(sgraphName[i].c_str());
}
// Table 1
graph1a->Write("graph1a");
graph1b->Write("graph1b");
sgraph1a->Write("sgraph1a");
sgraph1b->Write("sgraph1b");
hfile->Close();
}
void compareJpsi(string& rootfile)
{
//
// J/Psi HERA/ZEUS data
//
double heraProtonEnergy = 920;
double heraElectronEnergy = 27.5;
vector zeusDataJpsi;
loadJpsiData(zeusDataJpsi);
double chi2=0;
double ndof=0;
// HERA data - table 1
TGraphAsymmErrors *graph1H = new TGraphAsymmErrors(3); graph1H->SetTitle("0.15 < Q2 < 0.8");
TGraphAsymmErrors *graph2H = new TGraphAsymmErrors(6); graph2H->SetTitle("2 < Q2 < 5");
TGraphAsymmErrors *graph3H = new TGraphAsymmErrors(6); graph3H->SetTitle("5 < Q2 < 10");
TGraphAsymmErrors *graph4H = new TGraphAsymmErrors(6); graph4H->SetTitle("10 < Q2 < 100");
// Sartre results - table 1
TGraphAsymmErrors *graph1S = new TGraphAsymmErrors(3); graph1S->SetTitle("0.15 < Q2 < 0.8");
TGraphAsymmErrors *graph2S = new TGraphAsymmErrors(6); graph2S->SetTitle("2 < Q2 < 5");
TGraphAsymmErrors *graph3S = new TGraphAsymmErrors(6); graph3S->SetTitle("5 < Q2 < 10");
TGraphAsymmErrors *graph4S = new TGraphAsymmErrors(6); graph4S->SetTitle("10 < Q2 < 100");
// HERA data - table 2
TGraphAsymmErrors *graph5H = new TGraphAsymmErrors(8); graph5H->SetTitle("30 < W < 220");
TGraphAsymmErrors *graph6H = new TGraphAsymmErrors(7); graph6H->SetTitle("45 < W < 160");
// Sartre results - table 2
TGraphAsymmErrors *graph5S = new TGraphAsymmErrors(8); graph5S->SetTitle("30 < W < 220");
TGraphAsymmErrors *graph6S = new TGraphAsymmErrors(7); graph6S->SetTitle("45 < W < 160");
//
// Loop over data points and for each initialize Sartre and
// calculate the cross-section. We are using no runcard but
// set things for each data point
//
Sartre sartre;
EventGeneratorSettings* settings = sartre.runSettings();
double diff;
for (unsigned int i=0; isetVerbose(false);
- settings->setVerboseLevel(1);
+ settings->setVerbose(true);
+ settings->setVerboseLevel(2);
settings->setNumberOfEvents(0);
settings->setTimesToShow(0);
settings->setQ2min(zeusDataJpsi[i].Q2min);
settings->setQ2max(zeusDataJpsi[i].Q2max);
settings->setWmin(zeusDataJpsi[i].Wmin);
settings->setWmax(zeusDataJpsi[i].Wmax);
settings->setVectorMesonId(443);
settings->setApplyPhotonFlux(true);
settings->setElectronBeamEnergy(heraElectronEnergy);
settings->setHadronBeamEnergy(heraProtonEnergy);
- settings->setDipoleModelType(bNonSat);
+ settings->setDipoleModelType(theDipoleModelType);
+ settings->setDipoleModelParameterSet(theDipoleModelParameterSet);
settings->setA(1);
- settings->setCorrectForRealAmplitude(true);
- settings->setCorrectSkewedness(true);
+ settings->setCorrectForRealAmplitude(true);
+ settings->setCorrectSkewedness(true);
settings->setEnableNuclearBreakup(false);
settings->setMaxLambdaUsedInCorrections(maxLam);
+ settings->setTableSetType(theTableSetType);
- // settings->list();
bool ok = sartre.init();
if (!ok) {
cout << "Initialization of sartre failed. Going to next data point." << endl;
continue;
- }
-
+ }
+ settings->list();
+
cout << "Total cross-sections for Q2=[" << zeusDataJpsi[i].Q2min << ", " << zeusDataJpsi[i].Q2max << "], ";
cout << "W=[" << zeusDataJpsi[i].Wmin << ", " << zeusDataJpsi[i].Wmax<< "], |t| < 1" << endl;
cout << "\t\tSartre: " << sartre.totalCrossSection()*1000 << endl; // in pb
cout << "\t\tZEUS: " << zeusDataJpsi[i].xsection << "+/-" << zeusDataJpsi[i].xsection_stat << "+" <= 3 && i <= 8) {
graphS = graph2S;
graphH = graph2H;
offset = 3;
}
else if (i >= 9 && i <= 14) {
graphS = graph3S;
graphH = graph3H;
offset = 9;
}
else if (i >= 15 && i <= 20) {
graphS = graph4S;
graphH = graph4H;
offset = 15;
}
else if (i >= 21 && i <= 28) {
graphSS = graph5S;
graphHH = graph5H;
offset = 21;
}
else if (i >= 29 && i <= 35) {
graphSS = graph6S;
graphHH = graph6H;
offset = 29;
}
double meanW, meanQ2, exl, exu, eyl, eyu;
eyl = sqrt(zeusDataJpsi[i].xsection_stat*zeusDataJpsi[i].xsection_stat + zeusDataJpsi[i].xsection_sys_lo*zeusDataJpsi[i].xsection_sys_lo);
eyu = sqrt(zeusDataJpsi[i].xsection_stat*zeusDataJpsi[i].xsection_stat + zeusDataJpsi[i].xsection_sys_up*zeusDataJpsi[i].xsection_sys_up);
if (graphS && graphH) {
meanW = (zeusDataJpsi[i].Wmax+zeusDataJpsi[i].Wmin)/2;
graphS->SetPoint(i-offset, meanW, 1000*sartre.totalCrossSection());
graphS->SetPointError(i-offset, 0, 0, 0, 0);
exl = meanW-zeusDataJpsi[i].Wmin;
exu = zeusDataJpsi[i].Wmax-meanW;
graphH->SetPoint(i-offset, meanW, zeusDataJpsi[i].xsection);
graphH->SetPointError(i-offset, exl, exu, eyl, eyu);
}
else if (graphSS && graphHH) {
meanQ2 = (zeusDataJpsi[i].Q2max+zeusDataJpsi[i].Q2min)/2;
graphSS->SetPoint(i-offset, meanQ2, 1000*sartre.totalCrossSection());
graphSS->SetPointError(i-offset, 0, 0, 0, 0);
exl = meanQ2-zeusDataJpsi[i].Q2min;
exu = zeusDataJpsi[i].Q2max-meanQ2;
graphHH->SetPoint(i-offset, meanQ2, zeusDataJpsi[i].xsection);
graphHH->SetPointError(i-offset, exl, exu, eyl, eyu);
}
diff = 1000*sartre.totalCrossSection() - zeusDataJpsi[i].xsection;
if ( zeusDataJpsi[i].xsection > 1000*sartre.totalCrossSection() ) {
chi2 += (diff*diff) / (eyl*eyl);
chi2_global+= (diff*diff) / (eyl*eyl);
}
else {
chi2 += (diff*diff) / (eyu*eyu);
chi2_global += (diff*diff) / (eyu*eyu);
}
ndof++;
}
cout << "-------------------------------------------------" << endl;
cout << "Comparison HERA-Sartre:" << endl;
cout << "chi2 = " << chi2 << endl;
cout << "ndof = " << ndof << endl;
cout << "chi2/ndof = " << chi2/ndof << endl;
cout << "-------------------------------------------------" << endl;
//
// ROOT file
//
TFile *hfile = 0;
hfile = new TFile(rootfile.c_str(),"RECREATE");
cout << "ROOT file is '" << rootfile.c_str() << "'." << endl;
graph1H->Write("graph1H");
graph1S->Write("graph1S");
graph2H->Write("graph2H");
graph2S->Write("graph2S");
graph3H->Write("graph3H");
graph3S->Write("graph3S");
graph4H->Write("graph4H");
graph4S->Write("graph4S");
graph5H->Write("graph5H");
graph5S->Write("graph5S");
graph6H->Write("graph6H");
graph6S->Write("graph6S");
hfile->Close();
return;
}
void loadJpsiData(vector& vec)
{
//
// J/Psi HERA/ZEUS data
//
// Data for e p -> e' p' J/psi
// Journal reference: Nucl.Phys.B695:3-37,2004
// arXiv:hep-ex/0404008v1
//
// All data are for |t|<1 GeV2
// Cross section in pb
//
vector zeusDataJpsi(21+15);
unsigned int i = 0;
// table 1 (21 entries)
zeusDataJpsi[i].Q2min=0.15; zeusDataJpsi[i].Q2max=0.8; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=65; zeusDataJpsi[i].xsection=217; zeusDataJpsi[i].xsection_stat=53; zeusDataJpsi[i].xsection_sys_up=12; zeusDataJpsi[i].xsection_sys_lo=19; i++;
zeusDataJpsi[i].Q2min=0.15; zeusDataJpsi[i].Q2max=0.8; zeusDataJpsi[i].Wmin=65; zeusDataJpsi[i].Wmax=105; zeusDataJpsi[i].xsection=257; zeusDataJpsi[i].xsection_stat=46; zeusDataJpsi[i].xsection_sys_up=18; zeusDataJpsi[i].xsection_sys_lo=17; i++;
zeusDataJpsi[i].Q2min=0.15; zeusDataJpsi[i].Q2max=0.8; zeusDataJpsi[i].Wmin=105; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=498; zeusDataJpsi[i].xsection_stat=89; zeusDataJpsi[i].xsection_sys_up=37; zeusDataJpsi[i].xsection_sys_lo=38; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=45; zeusDataJpsi[i].xsection=41.5; zeusDataJpsi[i].xsection_stat=8.4; zeusDataJpsi[i].xsection_sys_up=5.6; zeusDataJpsi[i].xsection_sys_lo=6.6; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=70; zeusDataJpsi[i].xsection=48.8; zeusDataJpsi[i].xsection_stat=5.2; zeusDataJpsi[i].xsection_sys_up=3.1; zeusDataJpsi[i].xsection_sys_lo=3.9; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=70; zeusDataJpsi[i].Wmax=90; zeusDataJpsi[i].xsection=36.4; zeusDataJpsi[i].xsection_stat=4.1; zeusDataJpsi[i].xsection_sys_up=10.5; zeusDataJpsi[i].xsection_sys_lo=3.0; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=90; zeusDataJpsi[i].Wmax=112; zeusDataJpsi[i].xsection=35.4; zeusDataJpsi[i].xsection_stat=4.0; zeusDataJpsi[i].xsection_sys_up=3.0; zeusDataJpsi[i].xsection_sys_lo=4.5; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=112; zeusDataJpsi[i].Wmax=145; zeusDataJpsi[i].xsection=44.7; zeusDataJpsi[i].xsection_stat=5.0; zeusDataJpsi[i].xsection_sys_up=9.0; zeusDataJpsi[i].xsection_sys_lo=4.3; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=145; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=76.5; zeusDataJpsi[i].xsection_stat=10.3; zeusDataJpsi[i].xsection_sys_up=11.5; zeusDataJpsi[i].xsection_sys_lo=5.1; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=50; zeusDataJpsi[i].xsection=19.6; zeusDataJpsi[i].xsection_stat=4.1; zeusDataJpsi[i].xsection_sys_up=3.9; zeusDataJpsi[i].xsection_sys_lo=1.9; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=50; zeusDataJpsi[i].Wmax=74; zeusDataJpsi[i].xsection=19.3; zeusDataJpsi[i].xsection_stat=2.2; zeusDataJpsi[i].xsection_sys_up=2.9; zeusDataJpsi[i].xsection_sys_lo=1.3; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=74; zeusDataJpsi[i].Wmax=96; zeusDataJpsi[i].xsection=15.6; zeusDataJpsi[i].xsection_stat=1.8; zeusDataJpsi[i].xsection_sys_up=1.6; zeusDataJpsi[i].xsection_sys_lo=1.4; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=96; zeusDataJpsi[i].Wmax=120; zeusDataJpsi[i].xsection=13.5; zeusDataJpsi[i].xsection_stat=1.7; zeusDataJpsi[i].xsection_sys_up=1.1; zeusDataJpsi[i].xsection_sys_lo=0.7; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=120; zeusDataJpsi[i].Wmax=150; zeusDataJpsi[i].xsection=14.9; zeusDataJpsi[i].xsection_stat=1.9; zeusDataJpsi[i].xsection_sys_up=1.1; zeusDataJpsi[i].xsection_sys_lo=1.3; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=150; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=27.9; zeusDataJpsi[i].xsection_stat=4.1; zeusDataJpsi[i].xsection_sys_up=4.5; zeusDataJpsi[i].xsection_sys_lo=1.4; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=55; zeusDataJpsi[i].xsection=10.9; zeusDataJpsi[i].xsection_stat=3.1; zeusDataJpsi[i].xsection_sys_up=0.8; zeusDataJpsi[i].xsection_sys_lo=1.0; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=55; zeusDataJpsi[i].Wmax=78; zeusDataJpsi[i].xsection=8.4; zeusDataJpsi[i].xsection_stat=1.2; zeusDataJpsi[i].xsection_sys_up=1.4; zeusDataJpsi[i].xsection_sys_lo=0.4; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=78; zeusDataJpsi[i].Wmax=100; zeusDataJpsi[i].xsection=8.6; zeusDataJpsi[i].xsection_stat=1.1; zeusDataJpsi[i].xsection_sys_up=0.9; zeusDataJpsi[i].xsection_sys_lo=1.4; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=100; zeusDataJpsi[i].Wmax=124; zeusDataJpsi[i].xsection=8.4; zeusDataJpsi[i].xsection_stat=1.1; zeusDataJpsi[i].xsection_sys_up=0.4; zeusDataJpsi[i].xsection_sys_lo=1.2; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=124; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=10.8; zeusDataJpsi[i].xsection_stat=1.4; zeusDataJpsi[i].xsection_sys_up=2.1; zeusDataJpsi[i].xsection_sys_lo=0.8; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=160; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=25.1; zeusDataJpsi[i].xsection_stat=3.8; zeusDataJpsi[i].xsection_sys_up=1.7; zeusDataJpsi[i].xsection_sys_lo=1.2; i++;
// table 2 (15 entries)
zeusDataJpsi[i].Q2min=0.15; zeusDataJpsi[i].Q2max=0.8; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=954; zeusDataJpsi[i].xsection_stat=108; zeusDataJpsi[i].xsection_sys_up=63; zeusDataJpsi[i].xsection_sys_lo=74; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=3.2; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=150; zeusDataJpsi[i].xsection_stat=14; zeusDataJpsi[i].xsection_sys_up=53; zeusDataJpsi[i].xsection_sys_lo=8; i++;
zeusDataJpsi[i].Q2min=3.2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=132; zeusDataJpsi[i].xsection_stat=12; zeusDataJpsi[i].xsection_sys_up=8; zeusDataJpsi[i].xsection_sys_lo=17; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=7; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=59.9; zeusDataJpsi[i].xsection_stat=6.1; zeusDataJpsi[i].xsection_sys_up=5.5; zeusDataJpsi[i].xsection_sys_lo=3.6; i++;
zeusDataJpsi[i].Q2min=7; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=42.6; zeusDataJpsi[i].xsection_stat=4.3; zeusDataJpsi[i].xsection_sys_up=4.7; zeusDataJpsi[i].xsection_sys_lo=5.0; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=15; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=36.7; zeusDataJpsi[i].xsection_stat=4.0; zeusDataJpsi[i].xsection_sys_up=1.4; zeusDataJpsi[i].xsection_sys_lo=2.9; i++;
zeusDataJpsi[i].Q2min=15; zeusDataJpsi[i].Q2max=40; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=29.3; zeusDataJpsi[i].xsection_stat=3.7; zeusDataJpsi[i].xsection_sys_up=2.0; zeusDataJpsi[i].xsection_sys_lo=4.7; i++;
zeusDataJpsi[i].Q2min=40; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=30; zeusDataJpsi[i].Wmax=220; zeusDataJpsi[i].xsection=4.5; zeusDataJpsi[i].xsection_stat=1.5; zeusDataJpsi[i].xsection_sys_up=0.5; zeusDataJpsi[i].xsection_sys_lo=1.1; i++;
zeusDataJpsi[i].Q2min=2; zeusDataJpsi[i].Q2max=3.2; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=96; zeusDataJpsi[i].xsection_stat=11; zeusDataJpsi[i].xsection_sys_up=5; zeusDataJpsi[i].xsection_sys_lo=14; i++;
zeusDataJpsi[i].Q2min=3.2; zeusDataJpsi[i].Q2max=5; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=91.6; zeusDataJpsi[i].xsection_stat=8.9; zeusDataJpsi[i].xsection_sys_up=12.2; zeusDataJpsi[i].xsection_sys_lo=6.6; i++;
zeusDataJpsi[i].Q2min=5; zeusDataJpsi[i].Q2max=7; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=48.7; zeusDataJpsi[i].xsection_stat=5.2; zeusDataJpsi[i].xsection_sys_up=1.2; zeusDataJpsi[i].xsection_sys_lo=2.5; i++;
zeusDataJpsi[i].Q2min=7; zeusDataJpsi[i].Q2max=10; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=32.5; zeusDataJpsi[i].xsection_stat=3.4; zeusDataJpsi[i].xsection_sys_up=4.4; zeusDataJpsi[i].xsection_sys_lo=2.4; i++;
zeusDataJpsi[i].Q2min=10; zeusDataJpsi[i].Q2max=15; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=24.1; zeusDataJpsi[i].xsection_stat=2.8; zeusDataJpsi[i].xsection_sys_up=1.2; zeusDataJpsi[i].xsection_sys_lo=1.6; i++;
zeusDataJpsi[i].Q2min=15; zeusDataJpsi[i].Q2max=40; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=18.4; zeusDataJpsi[i].xsection_stat=2.4; zeusDataJpsi[i].xsection_sys_up=0.9; zeusDataJpsi[i].xsection_sys_lo=1.4; i++;
zeusDataJpsi[i].Q2min=40; zeusDataJpsi[i].Q2max=100; zeusDataJpsi[i].Wmin=45; zeusDataJpsi[i].Wmax=160; zeusDataJpsi[i].xsection=2.2; zeusDataJpsi[i].xsection_stat=0.9; zeusDataJpsi[i].xsection_sys_up=0.4; zeusDataJpsi[i].xsection_sys_lo=0.6;
vec = zeusDataJpsi;
}
void loadDVCSData(vector& vec1, vector& vec4)
{
//
// DVCS HERA/ZEUS data
//
// Origin:
// Journal reference: Phys.Lett.B659:796-806,2008
// Report number: DESY-07-142
// Cite as: arXiv:0709.4114v1 [hep-ex]
vector h1DataDVCS_table1(4+5);
vector h1DataDVCS_table4(12+12);
unsigned int i = 0;
// Table 1a (4 points versus Q2)
h1DataDVCS_table1[i].Wmin=30; h1DataDVCS_table1[i].Wmax=140; h1DataDVCS_table1[i].W=82; // W is ref only
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=11; h1DataDVCS_table1[i].Q2=8.75; // Q2 is ref value only
h1DataDVCS_table1[i].xsection=3.59; h1DataDVCS_table1[i].xsection_stat=.21; h1DataDVCS_table1[i].xsection_sys=.41;
i++;
h1DataDVCS_table1[i].Wmin=30; h1DataDVCS_table1[i].Wmax=140;; h1DataDVCS_table1[i].W=82; // W is ref only
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=11; h1DataDVCS_table1[i].Q2max=20; h1DataDVCS_table1[i].Q2=15.5;
h1DataDVCS_table1[i].xsection=1.38; h1DataDVCS_table1[i].xsection_stat=.1; h1DataDVCS_table1[i].xsection_sys=.21;
i++;
h1DataDVCS_table1[i].Wmin=30; h1DataDVCS_table1[i].Wmax=140;; h1DataDVCS_table1[i].W=82; // W is ref only
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=20; h1DataDVCS_table1[i].Q2max=30; h1DataDVCS_table1[i].Q2=25;
h1DataDVCS_table1[i].xsection=.58; h1DataDVCS_table1[i].xsection_stat=.09; h1DataDVCS_table1[i].xsection_sys=.09;
i++;
h1DataDVCS_table1[i].Wmin=30; h1DataDVCS_table1[i].Wmax=140;; h1DataDVCS_table1[i].W=82; // W is ref only
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=30; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=55;
h1DataDVCS_table1[i].xsection=.13; h1DataDVCS_table1[i].xsection_stat=.03; h1DataDVCS_table1[i].xsection_sys=.04;
i++;
// Table 1b (5 points versus W)
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=8; // Q2 is ref only
h1DataDVCS_table1[i].Wmin=30; h1DataDVCS_table1[i].Wmax=60; h1DataDVCS_table1[i].W=45;
h1DataDVCS_table1[i].xsection=2.91; h1DataDVCS_table1[i].xsection_stat=.2; h1DataDVCS_table1[i].xsection_sys=.25;
i++;
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=8; // Q2 is ref only
h1DataDVCS_table1[i].Wmin=60; h1DataDVCS_table1[i].Wmax=80; h1DataDVCS_table1[i].W=70;// W is ref value only
h1DataDVCS_table1[i].xsection=3.96; h1DataDVCS_table1[i].xsection_stat=.32; h1DataDVCS_table1[i].xsection_sys=.37;
i++;
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=8; // Q2 is ref only
h1DataDVCS_table1[i].Wmin=80; h1DataDVCS_table1[i].Wmax=100; h1DataDVCS_table1[i].W=90;
h1DataDVCS_table1[i].xsection= 4.78; h1DataDVCS_table1[i].xsection_stat=.41; h1DataDVCS_table1[i].xsection_sys=.57;
i++;
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=8; // Q2 is ref only
h1DataDVCS_table1[i].Wmin=100; h1DataDVCS_table1[i].Wmax=120; h1DataDVCS_table1[i].W=110;
h1DataDVCS_table1[i].xsection= 5.55; h1DataDVCS_table1[i].xsection_stat=.57; h1DataDVCS_table1[i].xsection_sys=.88;
i++;
h1DataDVCS_table1[i].tmin=0; h1DataDVCS_table1[i].tmax=1;
h1DataDVCS_table1[i].Q2min=6.5; h1DataDVCS_table1[i].Q2max=80; h1DataDVCS_table1[i].Q2=8;
h1DataDVCS_table1[i].Wmin=120; h1DataDVCS_table1[i].Wmax=140; h1DataDVCS_table1[i].W=130;
h1DataDVCS_table1[i].xsection= 6.56; h1DataDVCS_table1[i].xsection_stat=1.17; h1DataDVCS_table1[i].xsection_sys=1.77;
i++;
vec1 = h1DataDVCS_table1;
/*
// Table 3a
// PROBLEM: H1 gives W = 82 w/o specifying the bin range
// and the same for Q=10. The call this reference value
// which is pretty meaningless. The bins here are a guess
// and are apparently not working.
// We leave those out for the comparison.
//
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=11; h1DataDVCS[i].Q2=8.;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=13.1; h1DataDVCS[i].xsection_stat=1.10; h1DataDVCS[i].xsection_sys=1.85;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=11; h1DataDVCS[i].Q2=8.;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=4.69; h1DataDVCS[i].xsection_stat=0.45; h1DataDVCS[i].xsection_sys=0.55;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=11; h1DataDVCS[i].Q2=8.;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=1.37; h1DataDVCS[i].xsection_stat=0.21; h1DataDVCS[i].xsection_sys=0.23;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=11; h1DataDVCS[i].Q2=8.;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.19; h1DataDVCS[i].xsection_stat=0.04; h1DataDVCS[i].xsection_sys=0.06;
i++;
//-
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=11; h1DataDVCS[i].Q2max=20; h1DataDVCS[i].Q2=15.5;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=4.37; h1DataDVCS[i].xsection_stat=0.47; h1DataDVCS[i].xsection_sys=0.86;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=11; h1DataDVCS[i].Q2max=20; h1DataDVCS[i].Q2=15.5;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=1.02; h1DataDVCS[i].xsection_stat=0.16; h1DataDVCS[i].xsection_sys=0.18;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=11; h1DataDVCS[i].Q2max=20; h1DataDVCS[i].Q2=15.5;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=0.49; h1DataDVCS[i].xsection_stat=0.08; h1DataDVCS[i].xsection_sys=0.08;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=11; h1DataDVCS[i].Q2max=20; h1DataDVCS[i].Q2=15.5;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.12; h1DataDVCS[i].xsection_stat=0.02; h1DataDVCS[i].xsection_sys=0.02;
i++;
//-
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=20; h1DataDVCS[i].Q2max=80; h1DataDVCS[i].Q2=25;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=1.41; h1DataDVCS[i].xsection_stat=0.40; h1DataDVCS[i].xsection_sys=0.43;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=20; h1DataDVCS[i].Q2max=80; h1DataDVCS[i].Q2=25;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=0.71; h1DataDVCS[i].xsection_stat=0.16; h1DataDVCS[i].xsection_sys=0.08;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=20; h1DataDVCS[i].Q2max=80; h1DataDVCS[i].Q2=25;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=0.28; h1DataDVCS[i].xsection_stat=0.07; h1DataDVCS[i].xsection_sys=0.04;
i++;
h1DataDVCS[i].W=82; h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=140;
h1DataDVCS[i].Q2min=20; h1DataDVCS[i].Q2max=80; h1DataDVCS[i].Q2=25;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.04; h1DataDVCS[i].xsection_stat=0.01; h1DataDVCS[i].xsection_sys=0.02;
i++;
// Table 3b
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=50; h1DataDVCS[i].W=40;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=4.99; h1DataDVCS[i].xsection_stat=0.66; h1DataDVCS[i].xsection_sys=0.54;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=50; h1DataDVCS[i].W=40;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=1.45; h1DataDVCS[i].xsection_stat=0.29; h1DataDVCS[i].xsection_sys=0.18;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=50; h1DataDVCS[i].W=40;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=0.49; h1DataDVCS[i].xsection_stat=0.14; h1DataDVCS[i].xsection_sys=0.08;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=30; h1DataDVCS[i].Wmax=50; h1DataDVCS[i].W=40;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.12; h1DataDVCS[i].xsection_stat=0.03; h1DataDVCS[i].xsection_sys=0.03;
i++;
//-
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=50; h1DataDVCS[i].Wmax=85; h1DataDVCS[i].W=70;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=7.78; h1DataDVCS[i].xsection_stat=0.69; h1DataDVCS[i].xsection_sys=0.87;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=50; h1DataDVCS[i].Wmax=85; h1DataDVCS[i].W=70;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=2.74; h1DataDVCS[i].xsection_stat=0.31; h1DataDVCS[i].xsection_sys=0.3;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=50; h1DataDVCS[i].Wmax=85; h1DataDVCS[i].W=70;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=0.81; h1DataDVCS[i].xsection_stat=0.14; h1DataDVCS[i].xsection_sys=0.11;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=50; h1DataDVCS[i].Wmax=85; h1DataDVCS[i].W=70;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.19; h1DataDVCS[i].xsection_stat=0.03; h1DataDVCS[i].xsection_sys=0.03;
i++;
//-
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=85; h1DataDVCS[i].Wmax=140; h1DataDVCS[i].W=100;
h1DataDVCS[i].tmin=0; h1DataDVCS[i].tmax=0.2; h1DataDVCS[i].t=0.1;
h1DataDVCS[i].xsection=10.9; h1DataDVCS[i].xsection_stat=1.14; h1DataDVCS[i].xsection_sys=2.36;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=85; h1DataDVCS[i].Wmax=140; h1DataDVCS[i].W=100;
h1DataDVCS[i].tmin=0.2; h1DataDVCS[i].tmax=0.4; h1DataDVCS[i].t=0.3;
h1DataDVCS[i].xsection=3.47; h1DataDVCS[i].xsection_stat=0.42; h1DataDVCS[i].xsection_sys=0.53;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=85; h1DataDVCS[i].Wmax=140; h1DataDVCS[i].W=100;
h1DataDVCS[i].tmin=0.4; h1DataDVCS[i].tmax=0.6; h1DataDVCS[i].t=0.5;
h1DataDVCS[i].xsection=1.49; h1DataDVCS[i].xsection_stat=0.21; h1DataDVCS[i].xsection_sys=0.24;
i++;
h1DataDVCS[i].Q2=10; h1DataDVCS[i].Q2min=6.5; h1DataDVCS[i].Q2max=80;
h1DataDVCS[i].Wmin=85; h1DataDVCS[i].Wmax=140; h1DataDVCS[i].W=100;
h1DataDVCS[i].tmin=0.6; h1DataDVCS[i].tmax=1; h1DataDVCS[i].t=0.8;
h1DataDVCS[i].xsection=0.19; h1DataDVCS[i].xsection_stat=0.04; h1DataDVCS[i].xsection_sys=0.06;
i++;
*/
// Table 4a
i=0;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=8.10; h1DataDVCS_table4[i].xsection_stat=1.22; h1DataDVCS_table4[i].xsection_sys=0.82;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=2.30; h1DataDVCS_table4[i].xsection_stat=0.54; h1DataDVCS_table4[i].xsection_sys=0.28;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=0.45; h1DataDVCS_table4[i].xsection_stat=0.22; h1DataDVCS_table4[i].xsection_sys=0.10;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.16; h1DataDVCS_table4[i].xsection_stat=0.06; h1DataDVCS_table4[i].xsection_sys=0.03;
i++;
//-
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=10.0; h1DataDVCS_table4[i].xsection_stat=1.30; h1DataDVCS_table4[i].xsection_sys=1.27;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=4.35; h1DataDVCS_table4[i].xsection_stat=0.63; h1DataDVCS_table4[i].xsection_sys=0.46;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=1.08; h1DataDVCS_table4[i].xsection_stat=0.27; h1DataDVCS_table4[i].xsection_sys=0.17;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.13; h1DataDVCS_table4[i].xsection_stat=0.06; h1DataDVCS_table4[i].xsection_sys=0.04;
i++;
//-
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=16.0; h1DataDVCS_table4[i].xsection_stat=2.11; h1DataDVCS_table4[i].xsection_sys=2.74;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=5.45; h1DataDVCS_table4[i].xsection_stat=0.80; h1DataDVCS_table4[i].xsection_sys=0.73;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=1.96; h1DataDVCS_table4[i].xsection_stat=0.41; h1DataDVCS_table4[i].xsection_sys=0.35;
i++;
h1DataDVCS_table4[i].Q2min=6.5; h1DataDVCS_table4[i].Q2max=11; h1DataDVCS_table4[i].Q2=8.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.21; h1DataDVCS_table4[i].xsection_stat=0.09; h1DataDVCS_table4[i].xsection_sys=0.08;
i++;
// Table 4b
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=1.06; h1DataDVCS_table4[i].xsection_stat=0.28; h1DataDVCS_table4[i].xsection_sys=0.28;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=0.33; h1DataDVCS_table4[i].xsection_stat=0.07; h1DataDVCS_table4[i].xsection_sys=0.07;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=0.22; h1DataDVCS_table4[i].xsection_stat=0.06; h1DataDVCS_table4[i].xsection_sys=0.06;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=30; h1DataDVCS_table4[i].Wmax=50; h1DataDVCS_table4[i].W=40;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.04; h1DataDVCS_table4[i].xsection_stat=0.01; h1DataDVCS_table4[i].xsection_sys=0.01;
i++;
//-
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=2.38; h1DataDVCS_table4[i].xsection_stat=0.29; h1DataDVCS_table4[i].xsection_sys=0.26;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=0.67; h1DataDVCS_table4[i].xsection_stat=0.12; h1DataDVCS_table4[i].xsection_sys=0.07;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=0.24; h1DataDVCS_table4[i].xsection_stat=0.05; h1DataDVCS_table4[i].xsection_sys=0.03;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=50; h1DataDVCS_table4[i].Wmax=85; h1DataDVCS_table4[i].W=70;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.07; h1DataDVCS_table4[i].xsection_stat=0.01; h1DataDVCS_table4[i].xsection_sys=0.02;
i++;
//-
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0; h1DataDVCS_table4[i].tmax=0.2; h1DataDVCS_table4[i].t=0.1;
h1DataDVCS_table4[i].xsection=2.98; h1DataDVCS_table4[i].xsection_stat=0.49; h1DataDVCS_table4[i].xsection_sys=0.85;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.2; h1DataDVCS_table4[i].tmax=0.4; h1DataDVCS_table4[i].t=0.3;
h1DataDVCS_table4[i].xsection=0.89; h1DataDVCS_table4[i].xsection_stat=0.17; h1DataDVCS_table4[i].xsection_sys=0.17;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.4; h1DataDVCS_table4[i].tmax=0.6; h1DataDVCS_table4[i].t=0.5;
h1DataDVCS_table4[i].xsection=0.44; h1DataDVCS_table4[i].xsection_stat=0.08; h1DataDVCS_table4[i].xsection_sys=0.08;
i++;
h1DataDVCS_table4[i].Q2min=11; h1DataDVCS_table4[i].Q2max=80; h1DataDVCS_table4[i].Q2=20.;
h1DataDVCS_table4[i].Wmin=85; h1DataDVCS_table4[i].Wmax=140; h1DataDVCS_table4[i].W=100;
h1DataDVCS_table4[i].tmin=0.6; h1DataDVCS_table4[i].tmax=1; h1DataDVCS_table4[i].t=0.8;
h1DataDVCS_table4[i].xsection=0.06; h1DataDVCS_table4[i].xsection_stat=0.02; h1DataDVCS_table4[i].xsection_sys=0.02;
i++;
vec4 = h1DataDVCS_table4;
}
void loadPhiData(vector& vec1)
{
vector ZEUSDataPhi_table1(7+6+6+6+9);
unsigned int i = 0;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=45;
ZEUSDataPhi_table1[i].W=40;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=76.4;
ZEUSDataPhi_table1[i].xsection_stat=6.5;
ZEUSDataPhi_table1[i].xsection_sys=5.9;
i++;
ZEUSDataPhi_table1[i].Wmin=45;
ZEUSDataPhi_table1[i].Wmax=55;
ZEUSDataPhi_table1[i].W=50;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=101.2;
ZEUSDataPhi_table1[i].xsection_stat=7.7;
ZEUSDataPhi_table1[i].xsection_sys=8.1;
i++;
ZEUSDataPhi_table1[i].Wmin=55;
ZEUSDataPhi_table1[i].Wmax=65;
ZEUSDataPhi_table1[i].W=60;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=101.9;
ZEUSDataPhi_table1[i].xsection_stat=8.3;
ZEUSDataPhi_table1[i].xsection_sys=7.4;
i++;
ZEUSDataPhi_table1[i].Wmin=65;
ZEUSDataPhi_table1[i].Wmax=75;
ZEUSDataPhi_table1[i].W=70;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=112.8;
ZEUSDataPhi_table1[i].xsection_stat=9.4;
ZEUSDataPhi_table1[i].xsection_sys=7.6;
i++;
ZEUSDataPhi_table1[i].Wmin=75;
ZEUSDataPhi_table1[i].Wmax=85;
ZEUSDataPhi_table1[i].W=80;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=107.;
ZEUSDataPhi_table1[i].xsection_stat=11.;
ZEUSDataPhi_table1[i].xsection_sys=10.;
i++;
ZEUSDataPhi_table1[i].Wmin=85;
ZEUSDataPhi_table1[i].Wmax=95;
ZEUSDataPhi_table1[i].W=90;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=122;
ZEUSDataPhi_table1[i].xsection_stat=11.;
ZEUSDataPhi_table1[i].xsection_sys=10.;
i++;
ZEUSDataPhi_table1[i].Wmin=95;
ZEUSDataPhi_table1[i].Wmax=105;
ZEUSDataPhi_table1[i].W=100;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=110.;
ZEUSDataPhi_table1[i].xsection_stat=11.;
ZEUSDataPhi_table1[i].xsection_sys=8.;
i++;
//***** Q2 = 3.8 ******************
ZEUSDataPhi_table1[i].Wmin=40;
ZEUSDataPhi_table1[i].Wmax=50;
ZEUSDataPhi_table1[i].W=45;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=47.;
ZEUSDataPhi_table1[i].xsection_stat=4.1;
ZEUSDataPhi_table1[i].xsection_sys=2.5;
i++;
ZEUSDataPhi_table1[i].Wmin=50;
ZEUSDataPhi_table1[i].Wmax=60;
ZEUSDataPhi_table1[i].W=55;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=44.3;
ZEUSDataPhi_table1[i].xsection_stat=4.5;
ZEUSDataPhi_table1[i].xsection_sys=3.6;
i++;
ZEUSDataPhi_table1[i].Wmin=60;
ZEUSDataPhi_table1[i].Wmax=70;
ZEUSDataPhi_table1[i].W=65;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=56.7;
ZEUSDataPhi_table1[i].xsection_stat=5.1;
ZEUSDataPhi_table1[i].xsection_sys=4.;
i++;
ZEUSDataPhi_table1[i].Wmin=70;
ZEUSDataPhi_table1[i].Wmax=85;
ZEUSDataPhi_table1[i].W=77.5;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=62.3;
ZEUSDataPhi_table1[i].xsection_stat=5.;
ZEUSDataPhi_table1[i].xsection_sys=4.3;
i++;
ZEUSDataPhi_table1[i].Wmin=85;
ZEUSDataPhi_table1[i].Wmax=100;
ZEUSDataPhi_table1[i].W=92.5;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=57.4;
ZEUSDataPhi_table1[i].xsection_stat=5.4;
ZEUSDataPhi_table1[i].xsection_sys=4.3;
i++;
ZEUSDataPhi_table1[i].Wmin=100;
ZEUSDataPhi_table1[i].Wmax=115;
ZEUSDataPhi_table1[i].W=107.5;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=5;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=59.0;
ZEUSDataPhi_table1[i].xsection_stat=6.2;
ZEUSDataPhi_table1[i].xsection_sys=3.;
i++;
//*******Q2=6.5******************
ZEUSDataPhi_table1[i].Wmin=45;
ZEUSDataPhi_table1[i].Wmax=55;
ZEUSDataPhi_table1[i].W=50;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=16.4;
ZEUSDataPhi_table1[i].xsection_stat=1.8;
ZEUSDataPhi_table1[i].xsection_sys=1.1;
i++;
ZEUSDataPhi_table1[i].Wmin=55;
ZEUSDataPhi_table1[i].Wmax=70;
ZEUSDataPhi_table1[i].W=62.5;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=19.1;
ZEUSDataPhi_table1[i].xsection_stat=1.7;
ZEUSDataPhi_table1[i].xsection_sys=1.5;
i++;
ZEUSDataPhi_table1[i].Wmin=70;
ZEUSDataPhi_table1[i].Wmax=85;
ZEUSDataPhi_table1[i].W=77.5;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=19.6;
ZEUSDataPhi_table1[i].xsection_stat=1.9;
ZEUSDataPhi_table1[i].xsection_sys=1.6;
i++;
ZEUSDataPhi_table1[i].Wmin=85;
ZEUSDataPhi_table1[i].Wmax=100;
ZEUSDataPhi_table1[i].W=92.5;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=21.6;
ZEUSDataPhi_table1[i].xsection_stat=2.3;
ZEUSDataPhi_table1[i].xsection_sys=1.1;
i++;
ZEUSDataPhi_table1[i].Wmin=100;
ZEUSDataPhi_table1[i].Wmax=115;
ZEUSDataPhi_table1[i].W=107.5;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=23.1;
ZEUSDataPhi_table1[i].xsection_stat=2.5;
ZEUSDataPhi_table1[i].xsection_sys=1.3;
i++;
ZEUSDataPhi_table1[i].Wmin=115;
ZEUSDataPhi_table1[i].Wmax=135;
ZEUSDataPhi_table1[i].W=125;
ZEUSDataPhi_table1[i].Q2min=5;
ZEUSDataPhi_table1[i].Q2max=9;
ZEUSDataPhi_table1[i].Q2=6.5;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=25.3;
ZEUSDataPhi_table1[i].xsection_stat=3.5;
ZEUSDataPhi_table1[i].xsection_sys=3.9;
i++;
//*******Q2=13.0******************
ZEUSDataPhi_table1[i].Wmin=50;
ZEUSDataPhi_table1[i].Wmax=60;
ZEUSDataPhi_table1[i].W=55;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=5.05;
ZEUSDataPhi_table1[i].xsection_stat=.73;
ZEUSDataPhi_table1[i].xsection_sys=.37;
i++;
ZEUSDataPhi_table1[i].Wmin=60;
ZEUSDataPhi_table1[i].Wmax=75;
ZEUSDataPhi_table1[i].W=67.5;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=4.96;
ZEUSDataPhi_table1[i].xsection_stat=.59;
ZEUSDataPhi_table1[i].xsection_sys=.27;
i++;
ZEUSDataPhi_table1[i].Wmin=75;
ZEUSDataPhi_table1[i].Wmax=90;
ZEUSDataPhi_table1[i].W=82.5;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=6.12;
ZEUSDataPhi_table1[i].xsection_stat=.81;
ZEUSDataPhi_table1[i].xsection_sys=.4;
i++;
ZEUSDataPhi_table1[i].Wmin=90;
ZEUSDataPhi_table1[i].Wmax=105;
ZEUSDataPhi_table1[i].W=97.5;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=6.82;
ZEUSDataPhi_table1[i].xsection_stat=.85;
ZEUSDataPhi_table1[i].xsection_sys=.4;
i++;
ZEUSDataPhi_table1[i].Wmin=105;
ZEUSDataPhi_table1[i].Wmax=125;
ZEUSDataPhi_table1[i].W=115;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=6.33;
ZEUSDataPhi_table1[i].xsection_stat=.78;
ZEUSDataPhi_table1[i].xsection_sys=.5;
i++;
ZEUSDataPhi_table1[i].Wmin=125;
ZEUSDataPhi_table1[i].Wmax=145;
ZEUSDataPhi_table1[i].W=135;
ZEUSDataPhi_table1[i].Q2min=9;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=13.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=6.65;
ZEUSDataPhi_table1[i].xsection_stat=.97;
ZEUSDataPhi_table1[i].xsection_sys=.52;
i++;
//*****Q2 cross-section**********
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=2;
ZEUSDataPhi_table1[i].Q2max=3;
ZEUSDataPhi_table1[i].Q2=2.4;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=105.5;
ZEUSDataPhi_table1[i].xsection_stat=3.4;
ZEUSDataPhi_table1[i].xsection_sys=6.;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=3;
ZEUSDataPhi_table1[i].Q2max=4.5;
ZEUSDataPhi_table1[i].Q2=3.6;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=57.6;
ZEUSDataPhi_table1[i].xsection_stat=2.4;
ZEUSDataPhi_table1[i].xsection_sys=3.5;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=4.5;
ZEUSDataPhi_table1[i].Q2max=6;
ZEUSDataPhi_table1[i].Q2=5.2;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=31.1;
ZEUSDataPhi_table1[i].xsection_stat=1.8;
ZEUSDataPhi_table1[i].xsection_sys=1.8;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=6;
ZEUSDataPhi_table1[i].Q2max=8;
ZEUSDataPhi_table1[i].Q2=6.9;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=17.9;
ZEUSDataPhi_table1[i].xsection_stat=1.1;
ZEUSDataPhi_table1[i].xsection_sys=1.;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=8;
ZEUSDataPhi_table1[i].Q2max=11;
ZEUSDataPhi_table1[i].Q2=9.2;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=11.06;
ZEUSDataPhi_table1[i].xsection_stat=0.73;
ZEUSDataPhi_table1[i].xsection_sys=0.56;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=11;
ZEUSDataPhi_table1[i].Q2max=15;
ZEUSDataPhi_table1[i].Q2=12.6;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=6.42;
ZEUSDataPhi_table1[i].xsection_stat=0.52;
ZEUSDataPhi_table1[i].xsection_sys=0.24;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=15;
ZEUSDataPhi_table1[i].Q2max=20;
ZEUSDataPhi_table1[i].Q2=17.1;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=2.5;
ZEUSDataPhi_table1[i].xsection_stat=0.37;
ZEUSDataPhi_table1[i].xsection_sys=0.22;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=20;
ZEUSDataPhi_table1[i].Q2max=30;
ZEUSDataPhi_table1[i].Q2=24.;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=0.98;
ZEUSDataPhi_table1[i].xsection_stat=0.19;
ZEUSDataPhi_table1[i].xsection_sys=0.05;
i++;
ZEUSDataPhi_table1[i].Wmin=35;
ZEUSDataPhi_table1[i].Wmax=142;
ZEUSDataPhi_table1[i].W=75;
ZEUSDataPhi_table1[i].Q2min=30;
ZEUSDataPhi_table1[i].Q2max=70;
ZEUSDataPhi_table1[i].Q2=38.8;
ZEUSDataPhi_table1[i].tmin=0;
ZEUSDataPhi_table1[i].tmax=1;
ZEUSDataPhi_table1[i].xsection=0.37;
ZEUSDataPhi_table1[i].xsection_stat=0.13;
ZEUSDataPhi_table1[i].xsection_sys=0.04;
i++;
vec1 = ZEUSDataPhi_table1;
}