# This script illustrate how to run scripts with HepData (http://hepdata.cedar.ac.uk)
# This script was taken from this entry: 
# http://hepdata.cedar.ac.uk/h8test/view/ins1269454/

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('(1/N)*(1/TAA)*DN/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d1-x1-y1')
p1.add(85.0,56.86,5.0,5.0,5.22,5.22,0.0,0.0,14.04,14.04) 
p1.add(100.0,22.99,10.0,10.0,2.13,2.13,0.0,0.0,5.55,5.55) 
p1.add(120.0,9.076,10.0,10.0,0.863,0.863,0.0,0.0,1.923,1.923) 
p1.add(150.0,3.164,20.0,20.0,0.329,0.329,0.0,0.0,0.673,0.673) 
p1.add(210.0,0.279,40.0,40.0,0.056,0.056,0.0,0.0,0.061,0.061) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d1-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('(1/N)*(1/TAA)*DN/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d2-x1-y1')
p1.add(85.0,49.02,5.0,5.0,8.36,8.36,0.0,0.0,15.05,15.05) 
p1.add(100.0,17.22,10.0,10.0,2.96,2.96,0.0,0.0,5.9,5.9) 
p1.add(120.0,6.992,10.0,10.0,1.275,1.275,0.0,0.0,2.374,2.374) 
p1.add(150.0,1.947,20.0,20.0,0.492,0.492,0.0,0.0,0.803,0.803) 
p1.add(210.0,0.2207,40.0,40.0,0.102,0.102,0.0,0.0,0.1516,0.1516) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d2-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('(1/N)*(1/TAA)*DN/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d3-x1-y1')
p1.add(85.0,58.05,5.0,5.0,6.78,6.78,0.0,0.0,13.0,13.0) 
p1.add(100.0,22.25,10.0,10.0,2.52,2.52,0.0,0.0,4.71,4.71) 
p1.add(120.0,7.689,10.0,10.0,1.012,1.012,0.0,0.0,1.707,1.707) 
p1.add(150.0,3.644,20.0,20.0,0.478,0.478,0.0,0.0,0.786,0.786) 
p1.add(210.0,0.3873,40.0,40.0,0.0858,0.0858,0.0,0.0,0.1434,0.1434) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d3-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('(1/N)*(1/TAA)*DN/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d4-x1-y1')
p1.add(85.0,76.23,5.0,5.0,5.76,5.76,0.0,0.0,19.3,19.3) 
p1.add(100.0,30.4,10.0,10.0,2.4,2.4,0.0,0.0,7.62,7.62) 
p1.add(120.0,11.49,10.0,10.0,1.18,1.18,0.0,0.0,2.96,2.96) 
p1.add(150.0,3.73,20.0,20.0,0.576,0.576,0.0,0.0,0.881,0.881) 
p1.add(210.0,0.6224,40.0,40.0,0.1553,0.1553,0.0,0.0,0.1324,0.1324) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d4-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('(1/N)*(1/TAA)*DN/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d5-x1-y1')
p1.add(85.0,87.18,5.0,5.0,8.73,8.73,0.0,0.0,22.07,22.07) 
p1.add(100.0,38.86,10.0,10.0,4.33,4.33,0.0,0.0,9.74,9.74) 
p1.add(120.0,10.33,10.0,10.0,1.68,1.68,0.0,0.0,2.66,2.66) 
p1.add(150.0,3.736,20.0,20.0,0.945,0.945,0.0,0.0,0.882,0.882) 
p1.add(210.0,0.2226,40.0,40.0,0.1325,0.1325,0.0,0.0,0.0474,0.0474) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d5-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 2')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('DSIG/DPT IN PB/GEV')
p1=P1D('/HepData/10230/d6-x1-y1')
p1.add(85.0,124.8,5.0,5.0,1.8,1.8,0.0,0.0,19.9,19.9) 
p1.add(100.0,52.54,10.0,10.0,0.92,0.92,0.0,0.0,8.42,8.42) 
p1.add(120.0,15.21,10.0,10.0,0.45,0.45,0.0,0.0,2.56,2.56) 
p1.add(150.0,4.827,20.0,20.0,0.24,0.24,0.0,0.0,0.79,0.79) 
p1.add(210.0,0.5237,40.0,40.0,0.062,0.062,0.0,0.0,0.1063,0.1063) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d6-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 3')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d7-x1-y1')
p1.add(85.0,0.4555,5.0,5.0,0.0423,0.0423,0.0,0.0,0.0727,0.0727) 
p1.add(100.0,0.4376,10.0,10.0,0.0412,0.0412,0.0,0.0,0.0803,0.0803) 
p1.add(120.0,0.5968,10.0,10.0,0.0594,0.0594,0.0,0.0,0.1016,0.1016) 
p1.add(150.0,0.6556,20.0,20.0,0.0755,0.0755,0.0,0.0,0.1171,0.1171) 
p1.add(210.0,0.5322,40.0,40.0,0.124,0.124,0.0,0.0,0.1107,0.1107) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d7-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 4')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('CENTRALITY IN PCT')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d8-x1-y1')
p1.add(5.0,0.3927,5.0,5.0,0.0672,0.0672,0.0,0.0,0.09393279512502542,0.09393279512502542) 
p1.add(20.0,0.465,10.0,10.0,0.0547,0.0547,0.0,0.0,0.08989471619622591,0.08989471619622591) 
p1.add(40.0,0.6106,10.0,10.0,0.0469,0.0469,0.0,0.0,0.12674450678431787,0.12674450678431787) 
p1.add(75.0,0.6983,25.0,25.0,0.0706,0.0706,0.0,0.0,0.17828230422562977,0.17828230422562977) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d8-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Figure 4')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('CENTRALITY IN PCT')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d9-x1-y1')
p1.add(5.0,0.3277,5.0,5.0,0.0567,0.0567,0.0,0.0,0.12797909985618747,0.12797909985618747) 
p1.add(20.0,0.4234,10.0,10.0,0.0486,0.0486,0.0,0.0,0.1129998672565592,0.1129998672565592) 
p1.add(40.0,0.5785,10.0,10.0,0.0467,0.0467,0.0,0.0,0.16546782164517668,0.16546782164517668) 
p1.add(75.0,0.7395,25.0,25.0,0.0834,0.0834,0.0,0.0,0.1789160976547387,0.1789160976547387) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d9-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Additional data')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d10-x1-y1')
p1.add(85.0,0.3927,5.0,5.0,0.0672,0.0672,0.0,0.0,0.0917,0.0917) 
p1.add(100.0,0.3277,10.0,10.0,0.0567,0.0567,0.0,0.0,0.0924,0.0924) 
p1.add(120.0,0.4598,10.0,10.0,0.0849,0.0849,0.0,0.0,0.1272,0.1272) 
p1.add(150.0,0.4034,20.0,20.0,0.1038,0.1038,0.0,0.0,0.1529,0.1529) 
p1.add(210.0,0.4214,40.0,40.0,0.201,0.201,0.0,0.0,0.2771,0.2771) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d10-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Additional data')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d11-x1-y1')
p1.add(85.0,0.465,5.0,5.0,0.0547,0.0547,0.0,0.0,0.1062,0.1062) 
p1.add(100.0,0.4234,10.0,10.0,0.0486,0.0486,0.0,0.0,0.0859,0.0859) 
p1.add(120.0,0.5056,10.0,10.0,0.0682,0.0682,0.0,0.0,0.1104,0.1104) 
p1.add(150.0,0.755,20.0,20.0,0.1058,0.1058,0.0,0.0,0.1987,0.1987) 
p1.add(210.0,0.7395,40.0,40.0,0.1858,0.1858,0.0,0.0,0.4168,0.4168) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d11-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Additional data')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d12-x1-y1')
p1.add(85.0,0.6106,5.0,5.0,0.0469,0.0469,0.0,0.0,0.1161,0.1161) 
p1.add(100.0,0.5785,10.0,10.0,0.0467,0.0467,0.0,0.0,0.1129,0.1129) 
p1.add(120.0,0.7557,10.0,10.0,0.0806,0.0806,0.0,0.0,0.1562,0.1562) 
p1.add(150.0,0.7728,20.0,20.0,0.1254,0.1254,0.0,0.0,0.1698,0.1698) 
p1.add(210.0,1.188,40.0,40.0,0.328,0.328,0.0,0.0,0.246,0.246) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d12-x1-y1.eps')

from jhplot import *
c1 = HPlot('HepData')
c1.setGTitle('PRL 113,132301 (2014), Additional data')
c1.visible(1) # can be set to visible(0) to generate figure only without pop-up window
c1.setAutoRange()
c1.setNameX('PT(C=BJET) IN GEV')
c1.setNameY('RAA')
p1=P1D('/HepData/10230/d13-x1-y1')
p1.add(85.0,0.6983,5.0,5.0,0.0706,0.0706,0.0,0.0,0.1328,0.1328) 
p1.add(100.0,0.7395,10.0,10.0,0.0834,0.0834,0.0,0.0,0.1443,0.1443) 
p1.add(120.0,0.6793,10.0,10.0,0.1121,0.1121,0.0,0.0,0.1404,0.1404) 
p1.add(150.0,0.7739,20.0,20.0,0.1995,0.1995,0.0,0.0,0.17,0.17) 
p1.add(210.0,0.425,40.0,40.0,0.258,0.258,0.0,0.0,0.088,0.088) 
# draw and export
p1.setErrAll(1) # show all errors
c1.draw(p1)
c1.export('jhepwork-HepData-10230-d13-x1-y1.eps')