def load(i90CL=False):
idir = "datacards"
if i90CL:
idir="datacards_90"
#masses = [50,60,75,90,100,110, 125, 135, 150, 165, 180, 200, 250, 300]
masses = [100,110, 125, 135, 150, 165, 180, 200, 250, 300]
KFACTOR = 1.218
results = [];
for i in range(len(masses)):
results.append( getAsymLimits(idir+'/higgsCombineZprime'+str(masses[i])+'.Asymptotic.mH120.root','Zp'+str(masses[i])) );
names = [];
l_obs = [];
l_m2sig = [];
l_m1sig = [];
l_exp = [];
l_p1sig = [];
l_p2sig = [];
ctr = 0
for r in results:
names.append( "Zp"+str(masses[ctr]) );
l_m2sig.append(math.sqrt(r[1]));
l_m1sig.append(math.sqrt(r[2]));
l_exp.append (math.sqrt(r[3]));
l_p1sig.append(math.sqrt(r[4]));
l_p2sig.append(math.sqrt(r[5]));
l_obs.append (math.sqrt(r[0]));
ctr+=1;
gr_mu_exp = makeAGraph( masses, l_exp, 1, 2 );
gr_mu_obs = makeAGraph( masses, l_obs, 1, 1 );
gr_mu_1sigma = makeAFillGraph( masses, l_m1sig, l_p1sig, 0, 3, 1001 );
gr_mu_2sigma = makeAFillGraph( masses, l_m2sig, l_p2sig, 0, 5, 1001 );
return [gr_mu_exp,gr_mu_obs]
def getobs(obslo, obshi, xs):
x = []
y = []
iGDM = 1
iGQ = 0.25
iMDM = 1.
for i in range(obslo.GetN()):
lDMWidth = avtotwidth(2, iGDM, iGQ, obslo.GetX()[i], iMDM)
lWidth = avtotwidth(2, 0., iGQ, obslo.GetX()[i], iMDM)
factor = xs.Eval(obslo.GetX()[i]) / (6. / 5.)
x.append(obslo.GetX()[i])
y.append(
math.sqrt((lWidth / lDMWidth) * obslo.GetY()[i] / factor) / 4.)
for i in range(obshi.GetN()):
lDMWidth = avtotwidth(2, iGDM, iGQ, obshi.GetX()[i], iMDM)
lWidth = avtotwidth(2, 0., iGQ, obshi.GetX()[i], iMDM)
if obshi.GetX()[i] > 3700.: break
factor = xs.Eval(obshi.GetX()[i]) / (6. / 5.)
x.append(obshi.GetX()[i])
y.append(
math.sqrt((lWidth / lDMWidth) * obshi.GetY()[i] / factor) / 4.)
obs_gr = makeAGraph(x, y)
return obs_gr
def getobs(obslo,obshi,xs): x = []; y = []; iGDM=1;iGQ=0.25;iMDM=1.; for i in range( obslo.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,obslo.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,obslo.GetX()[i],iMDM) factor = xs.Eval( obslo.GetX()[i] ) / (6./5.); x.append( obslo.GetX()[i] ); y.append( math.sqrt( (lWidth/lDMWidth)*obslo.GetY()[i] / factor ) / 4. ); for i in range( obshi.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,obshi.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,obshi.GetX()[i],iMDM) if obshi.GetX()[i] > 3700.: break; factor = xs.Eval( obshi.GetX()[i] ) / (6./5.); x.append( obshi.GetX()[i] ); y.append( math.sqrt( (lWidth/lDMWidth)*obshi.GetY()[i] / factor ) / 4. ); obs_gr = makeAGraph( x, y ); return obs_gr
def getCont(iGraph,iVal,color):
lContours = iGraph.GetContourList(iVal);
xnew = array('d', []);
ynew = array('d', []);
for i0 in range(0,lContours.GetSize()):
pCurv = lContours.At(i0)
#if i0 > 0:
#break
for i2 in range(0,pCurv.GetN()):
lY = pCurv.GetY()[i2]
if lY < 10:
lY = -100
xnew.append(pCurv.GetX()[i2])
ynew.append(lY)
xnew.append(xnew[len(xnew)-1])
ynew.append(ynew[len(ynew)-1]+100000)
xnew.append(iGraph.GetX()[0])
ynew.append(ynew[len(ynew)-1])
xnew.append(iGraph.GetX()[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(ynew[0])
lGraph = makeAGraph( xnew, ynew, color, 1 )
lGraph.GetXaxis().SetTitle("m_{med} (GeV)")
lGraph.GetYaxis().SetTitle("m_{dm} (GeV)")
lGraph.SetLineWidth(4)
return lGraph
def getCont(iGraph, iVal, color):
lContours = iGraph.GetContourList(iVal)
xnew = array('d', [])
ynew = array('d', [])
for i0 in range(0, lContours.GetSize()):
pCurv = lContours.At(i0)
#if i0 > 0:
#break
for i2 in range(0, pCurv.GetN()):
lY = pCurv.GetY()[i2]
if lY < 10:
lY = -100
xnew.append(pCurv.GetX()[i2])
ynew.append(lY)
xnew.append(xnew[len(xnew) - 1])
ynew.append(ynew[len(ynew) - 1] + 100000)
xnew.append(iGraph.GetX()[0])
ynew.append(ynew[len(ynew) - 1])
xnew.append(iGraph.GetX()[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(ynew[0])
lGraph = makeAGraph(xnew, ynew, color, 1)
lGraph.GetXaxis().SetTitle("m_{med} (GeV)")
lGraph.GetYaxis().SetTitle("m_{dm} (GeV)")
lGraph.SetLineWidth(4)
return lGraph
def getexp(explo, exphi, xs):
x = []
y = []
yup = []
ydn = []
iGDM = 1
iGQ = 0.25
iMDM = 1.
for i in range(explo.GetN()):
lDMWidth = avtotwidth(2, iGDM, iGQ, explo.GetX()[i], iMDM)
lWidth = avtotwidth(2, 0., iGQ, explo.GetX()[i], iMDM)
if explo.GetX()[i] < 600.0: continue
factor = xs.Eval(explo.GetX()[i]) / (6. / 5.)
cury = explo.GetY()[i]
curyup = explo.GetY()[i] + explo.GetEYhigh()[i]
curydn = explo.GetY()[i] - explo.GetEYlow()[i]
cury /= factor
curyup /= factor
curydn /= factor
x.append(explo.GetX()[i])
y.append(math.sqrt((lWidth / lDMWidth) * cury) / 4.)
yup.append(math.sqrt((lWidth / lDMWidth) * curyup) / 4.)
ydn.append(math.sqrt((lWidth / lDMWidth) * curydn) / 4.)
for i in range(exphi.GetN()):
lDMWidth = avtotwidth(2, iGDM, iGQ, exphi.GetX()[i], iMDM)
lWidth = avtotwidth(2, 0., iGQ, exphi.GetX()[i], iMDM)
print(exphi.GetY()[i])
if exphi.GetX()[i] > 3700.: continue
factor = xs.Eval(exphi.GetX()[i]) / (6. / 5.)
cury = exphi.GetY()[i]
curyup = exphi.GetY()[i] + exphi.GetEYhigh()[i]
curydn = exphi.GetY()[i] - exphi.GetEYlow()[i]
cury /= factor
curyup /= factor
curydn /= factor
x.append(exphi.GetX()[i])
y.append(math.sqrt((lWidth / lDMWidth) * cury) / 4.)
yup.append(math.sqrt((lWidth / lDMWidth) * curyup) / 4.)
ydn.append(math.sqrt((lWidth / lDMWidth) * curydn) / 4.)
exp_gr = makeAGraph(x, y)
exp_gr_band = makeAFillGraph(x, ydn, yup)
return exp_gr, exp_gr_band
def convertgq(iGraph,iGDM,iMDM,iAxial,iConvertGB=True):
xnew = array('d', []);
ynew = array('d', []);
for i0 in range(0,iGraph.GetN()):
xnew.append(iGraph.GetX()[i0])
ynew.append(BRCorrGQ(iGraph.GetY()[i0]/6.,iGDM,iGraph.GetX()[i0],iMDM,iAxial))
if iConvertGB:
ynew[i0]=ynew[i0]
lGraph = makeAGraph( xnew, ynew, 1, 3);
return lGraph
def convertgq(iGraph, iGDM, iMDM, iAxial, iConvertGB=True):
xnew = array('d', [])
ynew = array('d', [])
for i0 in range(0, iGraph.GetN()):
xnew.append(iGraph.GetX()[i0])
ynew.append(
BRCorrGQ(iGraph.GetY()[i0] / 6., iGDM,
iGraph.GetX()[i0], iMDM, iAxial))
if iConvertGB:
ynew[i0] = ynew[i0]
lGraph = makeAGraph(xnew, ynew, 1, 3)
return lGraph
def getexp(explo,exphi,xs): x = []; y = []; yup = []; ydn = []; iGDM=1;iGQ=0.25;iMDM=1.; for i in range( explo.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,explo.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,explo.GetX()[i],iMDM) factor = xs.Eval( explo.GetX()[i] ) / (6./5.); cury = explo.GetY()[i]; curyup = explo.GetY()[i] + explo.GetEYhigh()[i]; curydn = explo.GetY()[i] - explo.GetEYlow()[i]; cury /= factor; curyup /= factor; curydn /= factor; x.append( explo.GetX()[i] ); y.append( math.sqrt((lWidth/lDMWidth)*cury) / 4. ); yup.append( math.sqrt((lWidth/lDMWidth)*curyup) / 4. ); ydn.append( math.sqrt((lWidth/lDMWidth)*curydn) / 4. ); for i in range( exphi.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,exphi.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,exphi.GetX()[i],iMDM) if exphi.GetX()[i] > 3700.: break; factor = xs.Eval( exphi.GetX()[i] ) / (6./5.); cury = exphi.GetY()[i]; curyup = exphi.GetY()[i] + exphi.GetEYhigh()[i]; curydn = exphi.GetY()[i] - exphi.GetEYlow()[i]; cury /= factor; curyup /= factor; curydn /= factor; x.append( exphi.GetX()[i] ); y.append( math.sqrt((lWidth/lDMWidth)*cury) / 4. ); yup.append( math.sqrt((lWidth/lDMWidth)*curyup) / 4. ); ydn.append( math.sqrt((lWidth/lDMWidth)*curydn) / 4. ); exp_gr = makeAGraph( x, y ); exp_gr_band = makeAFillGraph( x, ydn, yup ); return exp_gr, exp_gr_band
def getexp(explo,exphi,xs): x = []; y = []; yup = []; ydn = []; iGDM=1;iGQ=0.25;iMDM=1.; for i in range( explo.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,explo.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,explo.GetX()[i],iMDM) factor = xs.Eval( explo.GetX()[i] ) / (6./5.); cury = explo.GetY()[i]; curyup = explo.GetY()[i] + explo.GetEYhigh()[i]; curydn = explo.GetY()[i] - explo.GetEYlow()[i]; cury /= factor; curyup /= factor; curydn /= factor; x.append( explo.GetX()[i] ); y.append( math.sqrt((lWidth/lDMWidth)*cury) / 4. ); yup.append( math.sqrt((lWidth/lDMWidth)*curyup) / 4. ); ydn.append( math.sqrt((lWidth/lDMWidth)*curydn) / 4. ); for i in range( exphi.GetN() ): lDMWidth = avtotwidth(2,iGDM,iGQ,exphi.GetX()[i],iMDM) lWidth = avtotwidth(2,0. ,iGQ,exphi.GetX()[i],iMDM) if exphi.GetX()[i] > 3700.: break; factor = xs.Eval( exphi.GetX()[i] ) / (6./5.); cury = exphi.GetY()[i]; curyup = exphi.GetY()[i] + exphi.GetEYhigh()[i]; curydn = exphi.GetY()[i] - exphi.GetEYlow()[i]; cury /= factor; curyup /= factor; curydn /= factor; x.append( exphi.GetX()[i] ); y.append( math.sqrt((lWidth/lDMWidth)*cury) / 4. ); yup.append( math.sqrt((lWidth/lDMWidth)*curyup) / 4. ); ydn.append( math.sqrt((lWidth/lDMWidth)*curydn) / 4. ); exp_gr = makeAGraph( x, y ); exp_gr_band = makeAFillGraph( x, ydn, yup ); return exp_gr, exp_gr_band
def load(i90CL=False):
idir = "datacards"
if i90CL:
idir = "datacards_90"
#masses = [50,60,75,90,100,110, 125, 135, 150, 165, 180, 200, 250, 300]
masses = [100, 110, 125, 135, 150, 165, 180, 200, 250, 300]
KFACTOR = 1.218
results = []
for i in range(len(masses)):
results.append(
getAsymLimits(
idir + '/higgsCombineZprime' + str(masses[i]) +
'.Asymptotic.mH120.root', 'Zp' + str(masses[i])))
names = []
l_obs = []
l_m2sig = []
l_m1sig = []
l_exp = []
l_p1sig = []
l_p2sig = []
ctr = 0
for r in results:
names.append("Zp" + str(masses[ctr]))
l_m2sig.append(math.sqrt(r[1]))
l_m1sig.append(math.sqrt(r[2]))
l_exp.append(math.sqrt(r[3]))
l_p1sig.append(math.sqrt(r[4]))
l_p2sig.append(math.sqrt(r[5]))
l_obs.append(math.sqrt(r[0]))
ctr += 1
gr_mu_exp = makeAGraph(masses, l_exp, 1, 2)
gr_mu_obs = makeAGraph(masses, l_obs, 1, 1)
gr_mu_1sigma = makeAFillGraph(masses, l_m1sig, l_p1sig, 0, 3, 1001)
gr_mu_2sigma = makeAFillGraph(masses, l_m2sig, l_p2sig, 0, 5, 1001)
return [gr_mu_exp, gr_mu_obs]
def dijetxs():
## these cross-sections are for gq = 0.25
x = array('d', [])
y = array('d', [])
x.append(500)
y.append(50.11760833)
x.append(1000)
y.append(3.6878968)
x.append(1500)
y.append(0.6309465)
x.append(1800)
y.append(0.25842672)
x.append(1900)
y.append(0.196653793)
x.append(2000)
y.append(0.151528)
x.append(2100)
y.append(0.113179955)
x.append(2200)
y.append(0.088410105)
x.append(2300)
y.append(0.067159167)
x.append(2400)
y.append(0.053936451)
x.append(2500)
y.append(0.041326944)
x.append(2600)
y.append(0.03283815)
x.append(2700)
y.append(0.026234251)
x.append(2800)
y.append(0.020631734)
x.append(3000)
y.append(0.012939516)
x.append(3500)
y.append(0.004114006)
x.append(4000)
y.append(0.001400655)
lGraph = makeAGraph(x, y, 1, 3)
lGraph.SetMarkerColor(1)
return lGraph
def dijetxs():
## these cross-sections are for gq = 0.25
x = array('d', [])
y = array('d', [])
x.append(528.7733914858716)
y.append(57.269239965262344)
x.append(631.2905390292997)
y.append(26.146418729332762)
x.append(795.7485275479014)
y.append(10.943588903119819)
x.append(980.9839334973528)
y.append(4.998039360562294)
x.append(1166.2193394468052)
y.append(2.2826513012206147)
x.append(1331.0361216730032)
y.append(1.3151916426227461)
x.append(1578.3102214269736)
y.append(0.6008150098368491)
x.append(1825.7147916200695)
y.append(0.30829460763931776)
x.append(2073.0541265936026)
y.append(0.1492637583468037)
x.append(2217.3544322672033)
y.append(0.09943936422165509)
x.append(2454.5170729888905)
y.append(0.0556698333169542)
x.append(2701.98687840155)
y.append(0.03027483510044679)
x.append(2928.9075896518298)
y.append(0.018491825891702274)
x.append(3104.1946246179077)
y.append(0.011967992892333293)
x.append(3310.4357712666806)
y.append(0.007309405483957529)
x.append(3475.513494371131)
y.append(0.005313471607526083)
x.append(3568.343211809438)
y.append(0.004337281648532702)
lGraph = makeAGraph(x, y, 1, 3)
lGraph.SetMarkerColor(1)
return lGraph
def dijetxs():
## these cross-sections are for gq = 0.25
x = array('d', [])
y = array('d', [])
x.append(500)
y.append(54.39976667)
x.append(1000)
y.append(3.79482905)
x.append(1500)
y.append(0.64184835)
x.append(1800)
y.append(0.256497348)
x.append(1900)
y.append(0.197961717)
x.append(2000)
y.append(0.15153154)
x.append(2100)
y.append(0.115242345)
x.append(2200)
y.append(0.090241305)
x.append(2300)
y.append(0.06822651)
x.append(2400)
y.append(0.05406678)
x.append(2500)
y.append(0.041757683)
x.append(2600)
y.append(0.033034138)
x.append(2700)
y.append(0.025036183)
x.append(2800)
y.append(0.020569049)
x.append(3000)
y.append(0.012877625)
x.append(3500)
y.append(0.004038073)
x.append(4000)
y.append(0.001372723)
lGraph = makeAGraph( x, y, 1, 3 )
lGraph.SetMarkerColor(1)
return lGraph
def dijetxs():
## these cross-sections are for gq = 0.25
x = array('d', [])
y = array('d', [])
x.append(528.7733914858716)
y.append(57.269239965262344)
x.append(631.2905390292997)
y.append(26.146418729332762)
x.append(795.7485275479014)
y.append(10.943588903119819)
x.append(980.9839334973528)
y.append(4.998039360562294)
x.append(1166.2193394468052)
y.append(2.2826513012206147)
x.append(1331.0361216730032)
y.append(1.3151916426227461)
x.append(1578.3102214269736)
y.append(0.6008150098368491)
x.append(1825.7147916200695)
y.append(0.30829460763931776)
x.append(2073.0541265936026)
y.append(0.1492637583468037)
x.append(2217.3544322672033)
y.append(0.09943936422165509)
x.append(2454.5170729888905)
y.append(0.0556698333169542)
x.append(2701.98687840155)
y.append(0.03027483510044679)
x.append(2928.9075896518298)
y.append(0.018491825891702274)
x.append(3104.1946246179077)
y.append(0.011967992892333293)
x.append(3310.4357712666806)
y.append(0.007309405483957529)
x.append(3475.513494371131)
y.append(0.005313471607526083)
x.append(3568.343211809438)
y.append(0.004337281648532702)
lGraph = makeAGraph( x, y, 1, 3 )
lGraph.SetMarkerColor(1)
return lGraph
def merge(iPGraph,iCGraph,color,iX=-1,iXMax=100000):
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(iCGraph.GetY()[0])
for i0 in range(0,iCGraph.GetN()):
if iCGraph.GetX()[i0] > iX and iX > 0:
break
x.append(iCGraph.GetX()[i0])
y.append(iCGraph.GetY()[i0])
for i0 in range(0,iPGraph.GetN()):
if iPGraph.GetX()[i0] > iXMax:
break
if iX < 0:
x.append(iPGraph.GetX()[i0])
y.append(iPGraph.GetY()[i0])
elif iPGraph.GetX()[i0] > iX:
x.append(iPGraph.GetX()[i0])
y.append(iPGraph.GetY()[i0])
lGraph = makeAGraph( x, y, color, 3 )
lGraph.SetMarkerColor(color)
return lGraph
def merge(iPGraph, iCGraph, color, iX=-1, iXMax=100000):
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(iCGraph.GetY()[0])
for i0 in range(0, iCGraph.GetN()):
if iCGraph.GetX()[i0] > iX and iX > 0:
break
x.append(iCGraph.GetX()[i0])
y.append(iCGraph.GetY()[i0])
for i0 in range(0, iPGraph.GetN()):
if iPGraph.GetX()[i0] > iXMax:
break
if iX < 0:
x.append(iPGraph.GetX()[i0])
y.append(iPGraph.GetY()[i0])
elif iPGraph.GetX()[i0] > iX:
x.append(iPGraph.GetX()[i0])
y.append(iPGraph.GetY()[i0])
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def dijetxs(color):
x = array('d', [])
y = array('d', [])
v = True
if v:
# x.append(500)
# y.append(53.088896)
# x.append(1000)
# y.append(4.36988822)
# x.append(1500)
# y.append(0.79225344)
# x.append(1800)
# y.append(0.328192761)
# x.append(1850)
# y.append(0.291486468)
# x.append(1900)
# y.append(0.255722285)
# x.append(1950)
# y.append(0.221615421)
# x.append(2000)
# y.append(0.197382687)
# x.append(2050)
# y.append(0.172207705)
# x.append(2100)
# y.append(0.151872695)
# x.append(2150)
# y.append(0.132865362)
# x.append(2200)
# y.append(0.119653168)
# x.append(2250)
# y.append(0.10592325)
# x.append(2300)
# y.append(0.093760045)
# x.append(2350)
# y.append(0.082462696)
# x.append(2400)
# y.append(0.072743034)
# x.append(2450)
# y.append(0.066333462)
# x.append(2500)
# y.append(0.058459076)
# x.append(2550)
# y.append(0.052267432)
# x.append(2600)
# y.append(0.0466491248)
# x.append(2650)
# y.append(0.0412177488)
# x.append(2700)
# y.append(0.0376700324)
# x.append(2750)
# y.append(0.0337084692)
# x.append(2800)
# y.append(0.0300574278)
# x.append(3000)
# y.append(0.0197948408)
# x.append(3500)
# y.append(0.0072222648)
# x.append(4000)
# y.append(0.00293975808)
# x.append(4500)
# y.append(0.00125176576)
# x.append(5000)
# y.append(0.00058481674)
# x.append(5500)
# y.append(0.000306748216)
# x.append(6000)
# y.append(0.00017600644)
# x.append(6500)
# y.append(0.0001053256)
# x.append(7000)
# y.append(6.6086196e-05)
x.append(500)
y.append(59.614102)
x.append(510)
y.append(57.1915217478)
x.append(520)
y.append(54.7689414956)
x.append(530)
y.append(52.3463612434)
x.append(540)
y.append(49.9237809912)
x.append(550)
y.append(47.501200739)
x.append(560)
y.append(45.0786204868)
x.append(570)
y.append(42.6560402346)
x.append(580)
y.append(40.2334599824)
x.append(590)
y.append(37.8108797301)
x.append(600)
y.append(35.3882994779)
x.append(610)
y.append(33.9502002263)
x.append(620)
y.append(32.5121009746)
x.append(630)
y.append(31.074001723)
x.append(640)
y.append(29.6359024713)
x.append(650)
y.append(28.1978032197)
x.append(660)
y.append(26.759703968)
x.append(670)
y.append(25.3216047163)
x.append(680)
y.append(23.8835054647)
x.append(690)
y.append(22.445406213)
x.append(700)
y.append(21.0073069614)
x.append(710)
y.append(20.1536182319)
x.append(720)
y.append(19.2999295024)
x.append(730)
y.append(18.4462407729)
x.append(740)
y.append(17.5925520434)
x.append(750)
y.append(16.7388633139)
x.append(760)
y.append(15.8851745844)
x.append(770)
y.append(15.0314858549)
x.append(780)
y.append(14.1777971254)
x.append(790)
y.append(13.3241083959)
x.append(800)
y.append(12.4704196664)
x.append(810)
y.append(11.9636504388)
x.append(820)
y.append(11.4568812113)
x.append(830)
y.append(10.9501119838)
x.append(840)
y.append(10.4433427563)
x.append(850)
y.append(9.93657352872)
x.append(860)
y.append(9.42980430119)
x.append(870)
y.append(8.92303507366)
x.append(880)
y.append(8.41626584613)
x.append(890)
y.append(7.9094966186)
x.append(900)
y.append(7.40272739107)
x.append(910)
y.append(7.10189754396)
x.append(920)
y.append(6.80106769685)
x.append(930)
y.append(6.50023784975)
x.append(940)
y.append(6.19940800264)
x.append(950)
y.append(5.89857815553)
x.append(960)
y.append(5.59774830843)
x.append(970)
y.append(5.29691846132)
x.append(980)
y.append(4.99608861421)
x.append(990)
y.append(4.69525876711)
x.append(1000)
y.append(4.39442892)
x.append(1010)
y.append(4.26341595418)
x.append(1020)
y.append(4.13240298836)
x.append(1030)
y.append(4.00139002254)
x.append(1040)
y.append(3.87037705672)
x.append(1050)
y.append(3.7393640909)
x.append(1060)
y.append(3.60835112507)
x.append(1070)
y.append(3.47733815925)
x.append(1080)
y.append(3.34632519343)
x.append(1090)
y.append(3.21531222761)
x.append(1100)
y.append(3.08429926179)
x.append(1110)
y.append(2.99234574494)
x.append(1120)
y.append(2.90039222808)
x.append(1130)
y.append(2.80843871122)
x.append(1140)
y.append(2.71648519437)
x.append(1150)
y.append(2.62453167751)
x.append(1160)
y.append(2.53257816066)
x.append(1170)
y.append(2.4406246438)
x.append(1180)
y.append(2.34867112695)
x.append(1190)
y.append(2.25671761009)
x.append(1200)
y.append(2.16476409323)
x.append(1210)
y.append(2.10022506682)
x.append(1220)
y.append(2.0356860404)
x.append(1230)
y.append(1.97114701398)
x.append(1240)
y.append(1.90660798757)
x.append(1250)
y.append(1.84206896115)
x.append(1260)
y.append(1.77752993473)
x.append(1270)
y.append(1.71299090831)
x.append(1280)
y.append(1.6484518819)
x.append(1290)
y.append(1.58391285548)
x.append(1300)
y.append(1.51937382906)
x.append(1310)
y.append(1.47407609523)
x.append(1320)
y.append(1.42877836141)
x.append(1330)
y.append(1.38348062758)
x.append(1340)
y.append(1.33818289376)
x.append(1350)
y.append(1.29288515993)
x.append(1360)
y.append(1.2475874261)
x.append(1370)
y.append(1.20228969228)
x.append(1380)
y.append(1.15699195845)
x.append(1390)
y.append(1.11169422463)
x.append(1400)
y.append(1.0663964908)
x.append(1410)
y.append(1.03460356172)
x.append(1420)
y.append(1.00281063264)
x.append(1430)
y.append(0.97101770356)
x.append(1440)
y.append(0.93922477448)
x.append(1450)
y.append(0.9074318454)
x.append(1460)
y.append(0.87563891632)
x.append(1470)
y.append(0.84384598724)
x.append(1480)
y.append(0.81205305816)
x.append(1490)
y.append(0.78026012908)
x.append(1500)
y.append(0.7484672)
x.append(1510)
y.append(0.729957298163)
x.append(1520)
y.append(0.711447396326)
x.append(1530)
y.append(0.692937494489)
x.append(1540)
y.append(0.674427592652)
x.append(1550)
y.append(0.655917690815)
x.append(1560)
y.append(0.637407788978)
x.append(1570)
y.append(0.618897887142)
x.append(1580)
y.append(0.600387985305)
x.append(1590)
y.append(0.581878083468)
x.append(1600)
y.append(0.563368181631)
x.append(1610)
y.append(0.549435854683)
x.append(1620)
y.append(0.535503527735)
x.append(1630)
y.append(0.521571200788)
x.append(1640)
y.append(0.50763887384)
x.append(1650)
y.append(0.493706546892)
x.append(1660)
y.append(0.479774219945)
x.append(1670)
y.append(0.465841892997)
x.append(1680)
y.append(0.451909566049)
x.append(1690)
y.append(0.437977239101)
x.append(1700)
y.append(0.424044912154)
x.append(1710)
y.append(0.413558106989)
x.append(1720)
y.append(0.403071301825)
x.append(1730)
y.append(0.39258449666)
x.append(1740)
y.append(0.382097691496)
x.append(1750)
y.append(0.371610886331)
x.append(1760)
y.append(0.361124081167)
x.append(1770)
y.append(0.350637276002)
x.append(1780)
y.append(0.340150470838)
x.append(1790)
y.append(0.329663665673)
x.append(1800)
y.append(0.319176860509)
x.append(1810)
y.append(0.311283485413)
x.append(1820)
y.append(0.303390110317)
x.append(1830)
y.append(0.295496735221)
x.append(1840)
y.append(0.287603360125)
x.append(1850)
y.append(0.279709985029)
x.append(1860)
y.append(0.271816609933)
x.append(1870)
y.append(0.263923234837)
x.append(1880)
y.append(0.256029859741)
x.append(1890)
y.append(0.248136484645)
x.append(1900)
y.append(0.240243109549)
x.append(1910)
y.append(0.234301798594)
x.append(1920)
y.append(0.228360487639)
x.append(1930)
y.append(0.222419176684)
x.append(1940)
y.append(0.216477865729)
x.append(1950)
y.append(0.210536554775)
x.append(1960)
y.append(0.20459524382)
x.append(1970)
y.append(0.198653932865)
x.append(1980)
y.append(0.19271262191)
x.append(1990)
y.append(0.186771310955)
x.append(2000)
y.append(0.18083)
x.append(2010)
y.append(0.17677271424)
x.append(2020)
y.append(0.17271542848)
x.append(2030)
y.append(0.16865814272)
x.append(2040)
y.append(0.164600856961)
x.append(2050)
y.append(0.160543571201)
x.append(2060)
y.append(0.156486285441)
x.append(2070)
y.append(0.152428999681)
x.append(2080)
y.append(0.148371713921)
x.append(2090)
y.append(0.144314428161)
x.append(2100)
y.append(0.140257142401)
x.append(2110)
y.append(0.137110190532)
x.append(2120)
y.append(0.133963238662)
x.append(2130)
y.append(0.130816286793)
x.append(2140)
y.append(0.127669334923)
x.append(2150)
y.append(0.124522383053)
x.append(2160)
y.append(0.121375431184)
x.append(2170)
y.append(0.118228479314)
x.append(2180)
y.append(0.115081527445)
x.append(2190)
y.append(0.111934575575)
x.append(2200)
y.append(0.108787623705)
x.append(2210)
y.append(0.106346753957)
x.append(2220)
y.append(0.10390588421)
x.append(2230)
y.append(0.101465014462)
x.append(2240)
y.append(0.0990241447138)
x.append(2250)
y.append(0.096583274966)
x.append(2260)
y.append(0.0941424052181)
x.append(2270)
y.append(0.0917015354702)
x.append(2280)
y.append(0.0892606657224)
x.append(2290)
y.append(0.0868197959745)
x.append(2300)
y.append(0.0843789262266)
x.append(2310)
y.append(0.0824857148357)
x.append(2320)
y.append(0.0805925034449)
x.append(2330)
y.append(0.078699292054)
x.append(2340)
y.append(0.0768060806631)
x.append(2350)
y.append(0.0749128692723)
x.append(2360)
y.append(0.0730196578814)
x.append(2370)
y.append(0.0711264464905)
x.append(2380)
y.append(0.0692332350997)
x.append(2390)
y.append(0.0673400237088)
x.append(2400)
y.append(0.0654468123179)
x.append(2410)
y.append(0.0639783810861)
x.append(2420)
y.append(0.0625099498543)
x.append(2430)
y.append(0.0610415186226)
x.append(2440)
y.append(0.0595730873908)
x.append(2450)
y.append(0.058104656159)
x.append(2460)
y.append(0.0566362249272)
x.append(2470)
y.append(0.0551677936954)
x.append(2480)
y.append(0.0536993624636)
x.append(2490)
y.append(0.0522309312318)
x.append(2500)
y.append(0.0507625)
x.append(2510)
y.append(0.0496778788981)
x.append(2520)
y.append(0.0485932577961)
x.append(2530)
y.append(0.0475086366942)
x.append(2540)
y.append(0.0464240155923)
x.append(2550)
y.append(0.0453393944904)
x.append(2560)
y.append(0.0442547733884)
x.append(2570)
y.append(0.0431701522865)
x.append(2580)
y.append(0.0420855311846)
x.append(2590)
y.append(0.0410009100827)
x.append(2600)
y.append(0.0399162889807)
x.append(2610)
y.append(0.0390634143323)
x.append(2620)
y.append(0.0382105396839)
x.append(2630)
y.append(0.0373576650356)
x.append(2640)
y.append(0.0365047903872)
x.append(2650)
y.append(0.0356519157388)
x.append(2660)
y.append(0.0347990410904)
x.append(2670)
y.append(0.033946166442)
x.append(2680)
y.append(0.0330932917936)
x.append(2690)
y.append(0.0322404171452)
x.append(2700)
y.append(0.0313875424968)
x.append(2710)
y.append(0.0307168980067)
x.append(2720)
y.append(0.0300462535166)
x.append(2730)
y.append(0.0293756090264)
x.append(2740)
y.append(0.0287049645363)
x.append(2750)
y.append(0.0280343200462)
x.append(2760)
y.append(0.0273636755561)
x.append(2770)
y.append(0.026693031066)
x.append(2780)
y.append(0.0260223865759)
x.append(2790)
y.append(0.0253517420857)
x.append(2800)
y.append(0.0246810975956)
x.append(2810)
y.append(0.0241537469082)
x.append(2820)
y.append(0.0236263962207)
x.append(2830)
y.append(0.0230990455333)
x.append(2840)
y.append(0.0225716948459)
x.append(2850)
y.append(0.0220443441584)
x.append(2860)
y.append(0.021516993471)
x.append(2870)
y.append(0.0209896427836)
x.append(2880)
y.append(0.0204622920961)
x.append(2890)
y.append(0.0199349414087)
x.append(2900)
y.append(0.0194075907213)
x.append(2910)
y.append(0.0189929168491)
x.append(2920)
y.append(0.018578242977)
x.append(2930)
y.append(0.0181635691049)
x.append(2940)
y.append(0.0177488952328)
x.append(2950)
y.append(0.0173342213606)
x.append(2960)
y.append(0.0169195474885)
x.append(2970)
y.append(0.0165048736164)
x.append(2980)
y.append(0.0160901997443)
x.append(2990)
y.append(0.0156755258721)
x.append(3000)
y.append(0.015260852)
x.append(3010)
y.append(0.0149485495827)
x.append(3020)
y.append(0.0146362471653)
x.append(3030)
y.append(0.014323944748)
x.append(3040)
y.append(0.0140116423306)
x.append(3050)
y.append(0.0136993399133)
x.append(3060)
y.append(0.013387037496)
x.append(3070)
y.append(0.0130747350786)
x.append(3080)
y.append(0.0127624326613)
x.append(3090)
y.append(0.012450130244)
x.append(3100)
y.append(0.0121378278266)
x.append(3110)
y.append(0.0118894358645)
x.append(3120)
y.append(0.0116410439024)
x.append(3130)
y.append(0.0113926519403)
x.append(3140)
y.append(0.0111442599782)
x.append(3150)
y.append(0.0108958680162)
x.append(3160)
y.append(0.0106474760541)
x.append(3170)
y.append(0.010399084092)
x.append(3180)
y.append(0.0101506921299)
x.append(3190)
y.append(0.00990230016778)
x.append(3200)
y.append(0.00965390820569)
x.append(3210)
y.append(0.00945634788144)
x.append(3220)
y.append(0.00925878755719)
x.append(3230)
y.append(0.00906122723295)
x.append(3240)
y.append(0.0088636669087)
x.append(3250)
y.append(0.00866610658445)
x.append(3260)
y.append(0.00846854626021)
x.append(3270)
y.append(0.00827098593596)
x.append(3280)
y.append(0.00807342561171)
x.append(3290)
y.append(0.00787586528747)
x.append(3300)
y.append(0.00767830496322)
x.append(3310)
y.append(0.00752117394582)
x.append(3320)
y.append(0.00736404292843)
x.append(3330)
y.append(0.00720691191103)
x.append(3340)
y.append(0.00704978089364)
x.append(3350)
y.append(0.00689264987624)
x.append(3360)
y.append(0.00673551885885)
x.append(3370)
y.append(0.00657838784145)
x.append(3380)
y.append(0.00642125682406)
x.append(3390)
y.append(0.00626412580666)
x.append(3400)
y.append(0.00610699478927)
x.append(3410)
y.append(0.00598201951034)
x.append(3420)
y.append(0.00585704423142)
x.append(3430)
y.append(0.00573206895249)
x.append(3440)
y.append(0.00560709367356)
x.append(3450)
y.append(0.00548211839463)
x.append(3460)
y.append(0.00535714311571)
x.append(3470)
y.append(0.00523216783678)
x.append(3480)
y.append(0.00510719255785)
x.append(3490)
y.append(0.00498221727893)
x.append(3500)
y.append(0.004857242)
x.append(3510)
y.append(0.00476170076135)
x.append(3520)
y.append(0.00466615952271)
x.append(3530)
y.append(0.00457061828406)
x.append(3540)
y.append(0.00447507704542)
x.append(3550)
y.append(0.00437953580677)
x.append(3560)
y.append(0.00428399456813)
x.append(3570)
y.append(0.00418845332948)
x.append(3580)
y.append(0.00409291209084)
x.append(3590)
y.append(0.00399737085219)
x.append(3600)
y.append(0.00390182961354)
x.append(3610)
y.append(0.00382508119659)
x.append(3620)
y.append(0.00374833277964)
x.append(3630)
y.append(0.00367158436268)
x.append(3640)
y.append(0.00359483594573)
x.append(3650)
y.append(0.00351808752878)
x.append(3660)
y.append(0.00344133911182)
x.append(3670)
y.append(0.00336459069487)
x.append(3680)
y.append(0.00328784227791)
x.append(3690)
y.append(0.00321109386096)
x.append(3700)
y.append(0.00313434544401)
x.append(3710)
y.append(0.00307269332824)
x.append(3720)
y.append(0.00301104121248)
x.append(3730)
y.append(0.00294938909672)
x.append(3740)
y.append(0.00288773698096)
x.append(3750)
y.append(0.0028260848652)
x.append(3760)
y.append(0.00276443274944)
x.append(3770)
y.append(0.00270278063368)
x.append(3780)
y.append(0.00264112851792)
x.append(3790)
y.append(0.00257947640216)
x.append(3800)
y.append(0.0025178242864)
x.append(3810)
y.append(0.0024682990515)
x.append(3820)
y.append(0.0024187738166)
x.append(3830)
y.append(0.0023692485817)
x.append(3840)
y.append(0.0023197233468)
x.append(3850)
y.append(0.0022701981119)
x.append(3860)
y.append(0.002220672877)
x.append(3870)
y.append(0.0021711476421)
x.append(3880)
y.append(0.0021216224072)
x.append(3890)
y.append(0.0020720971723)
x.append(3900)
y.append(0.0020225719374)
x.append(3910)
y.append(0.00198278824366)
x.append(3920)
y.append(0.00194300454992)
x.append(3930)
y.append(0.00190322085618)
x.append(3940)
y.append(0.00186343716244)
x.append(3950)
y.append(0.0018236534687)
x.append(3960)
y.append(0.00178386977496)
x.append(3970)
y.append(0.00174408608122)
x.append(3980)
y.append(0.00170430238748)
x.append(3990)
y.append(0.00166451869374)
x.append(4000)
y.append(0.001624735)
x.append(4010)
y.append(0.00159277670054)
x.append(4020)
y.append(0.00156081840109)
x.append(4030)
y.append(0.00152886010163)
x.append(4040)
y.append(0.00149690180217)
x.append(4050)
y.append(0.00146494350272)
x.append(4060)
y.append(0.00143298520326)
x.append(4070)
y.append(0.0014010269038)
x.append(4080)
y.append(0.00136906860434)
x.append(4090)
y.append(0.00133711030489)
x.append(4100)
y.append(0.00130515200543)
x.append(4110)
y.append(0.00127947985666)
x.append(4120)
y.append(0.0012538077079)
x.append(4130)
y.append(0.00122813555913)
x.append(4140)
y.append(0.00120246341036)
x.append(4150)
y.append(0.0011767912616)
x.append(4160)
y.append(0.00115111911283)
x.append(4170)
y.append(0.00112544696406)
x.append(4180)
y.append(0.0010997748153)
x.append(4190)
y.append(0.00107410266653)
x.append(4200)
y.append(0.00104843051776)
x.append(4210)
y.append(0.00102780804306)
x.append(4220)
y.append(0.00100718556835)
x.append(4230)
y.append(0.000986563093637)
x.append(4240)
y.append(0.000965940618928)
x.append(4250)
y.append(0.000945318144219)
x.append(4260)
y.append(0.00092469566951)
x.append(4270)
y.append(0.000904073194801)
x.append(4280)
y.append(0.000883450720091)
x.append(4290)
y.append(0.000862828245382)
x.append(4300)
y.append(0.000842205770673)
x.append(4310)
y.append(0.000825639706531)
x.append(4320)
y.append(0.00080907364239)
x.append(4330)
y.append(0.000792507578248)
x.append(4340)
y.append(0.000775941514106)
x.append(4350)
y.append(0.000759375449964)
x.append(4360)
y.append(0.000742809385822)
x.append(4370)
y.append(0.00072624332168)
x.append(4380)
y.append(0.000709677257538)
x.append(4390)
y.append(0.000693111193396)
x.append(4400)
y.append(0.000676545129255)
x.append(4410)
y.append(0.000663237585663)
x.append(4420)
y.append(0.000649930042072)
x.append(4430)
y.append(0.000636622498481)
x.append(4440)
y.append(0.000623314954889)
x.append(4450)
y.append(0.000610007411298)
x.append(4460)
y.append(0.000596699867707)
x.append(4470)
y.append(0.000583392324115)
x.append(4480)
y.append(0.000570084780524)
x.append(4490)
y.append(0.000556777236933)
x.append(4500)
y.append(0.000543469693341)
x.append(4510)
y.append(0.000532779724081)
x.append(4520)
y.append(0.000522089754821)
x.append(4530)
y.append(0.000511399785562)
x.append(4540)
y.append(0.000500709816302)
x.append(4550)
y.append(0.000490019847042)
x.append(4560)
y.append(0.000479329877782)
x.append(4570)
y.append(0.000468639908522)
x.append(4580)
y.append(0.000457949939262)
x.append(4590)
y.append(0.000447259970002)
x.append(4600)
y.append(0.000436570000742)
x.append(4610)
y.append(0.000427982732777)
x.append(4620)
y.append(0.000419395464811)
x.append(4630)
y.append(0.000410808196846)
x.append(4640)
y.append(0.000402220928881)
x.append(4650)
y.append(0.000393633660915)
x.append(4660)
y.append(0.00038504639295)
x.append(4670)
y.append(0.000376459124985)
x.append(4680)
y.append(0.000367871857019)
x.append(4690)
y.append(0.000359284589054)
x.append(4700)
y.append(0.000350697321089)
x.append(4710)
y.append(0.000343799156153)
x.append(4720)
y.append(0.000336900991218)
x.append(4730)
y.append(0.000330002826283)
x.append(4740)
y.append(0.000323104661348)
x.append(4750)
y.append(0.000316206496412)
x.append(4760)
y.append(0.000309308331477)
x.append(4770)
y.append(0.000302410166542)
x.append(4780)
y.append(0.000295512001607)
x.append(4790)
y.append(0.000288613836671)
x.append(4800)
y.append(0.000281715671736)
x.append(4810)
y.append(0.000276174365739)
x.append(4820)
y.append(0.000270633059742)
x.append(4830)
y.append(0.000265091753745)
x.append(4840)
y.append(0.000259550447748)
x.append(4850)
y.append(0.000254009141751)
x.append(4860)
y.append(0.000248467835754)
x.append(4870)
y.append(0.000242926529757)
x.append(4880)
y.append(0.00023738522376)
x.append(4890)
y.append(0.000231843917763)
x.append(4900)
y.append(0.000226302611766)
x.append(4910)
y.append(0.000221851272542)
x.append(4920)
y.append(0.000217399933317)
x.append(4930)
y.append(0.000212948594093)
x.append(4940)
y.append(0.000208497254869)
x.append(4950)
y.append(0.000204045915644)
x.append(4960)
y.append(0.00019959457642)
x.append(4970)
y.append(0.000195143237195)
x.append(4980)
y.append(0.000190691897971)
x.append(4990)
y.append(0.000186240558747)
x.append(5000)
y.append(0.000181789219522)
x.append(5010)
y.append(0.000178213452203)
x.append(5020)
y.append(0.000174637684884)
x.append(5030)
y.append(0.000171061917564)
x.append(5040)
y.append(0.000167486150245)
x.append(5050)
y.append(0.000163910382926)
x.append(5060)
y.append(0.000160334615606)
x.append(5070)
y.append(0.000156758848287)
x.append(5080)
y.append(0.000153183080968)
x.append(5090)
y.append(0.000149607313649)
x.append(5100)
y.append(0.000146031546329)
x.append(5110)
y.append(0.000143159127204)
x.append(5120)
y.append(0.000140286708078)
x.append(5130)
y.append(0.000137414288953)
x.append(5140)
y.append(0.000134541869828)
x.append(5150)
y.append(0.000131669450702)
x.append(5160)
y.append(0.000128797031577)
x.append(5170)
y.append(0.000125924612451)
x.append(5180)
y.append(0.000123052193326)
x.append(5190)
y.append(0.0001201797742)
x.append(5200)
y.append(0.000117307355075)
x.append(5210)
y.append(0.000114999936584)
x.append(5220)
y.append(0.000112692518093)
x.append(5230)
y.append(0.000110385099602)
x.append(5240)
y.append(0.000108077681111)
x.append(5250)
y.append(0.00010577026262)
x.append(5260)
y.append(0.000103462844129)
x.append(5270)
y.append(0.000101155425638)
x.append(5280)
y.append(9.88480071469e-05)
x.append(5290)
y.append(9.65405886559e-05)
x.append(5300)
y.append(9.42331701649e-05)
x.append(5310)
y.append(9.23796175112e-05)
x.append(5320)
y.append(9.05260648574e-05)
x.append(5330)
y.append(8.86725122036e-05)
x.append(5340)
y.append(8.68189595499e-05)
x.append(5350)
y.append(8.49654068961e-05)
x.append(5360)
y.append(8.31118542424e-05)
x.append(5370)
y.append(8.12583015886e-05)
x.append(5380)
y.append(7.94047489348e-05)
x.append(5390)
y.append(7.75511962811e-05)
x.append(5400)
y.append(7.56976436273e-05)
x.append(5410)
y.append(7.42086820654e-05)
x.append(5420)
y.append(7.27197205036e-05)
x.append(5430)
y.append(7.12307589417e-05)
x.append(5440)
y.append(6.97417973798e-05)
x.append(5450)
y.append(6.8252835818e-05)
x.append(5460)
y.append(6.67638742561e-05)
x.append(5470)
y.append(6.52749126942e-05)
x.append(5480)
y.append(6.37859511324e-05)
x.append(5490)
y.append(6.22969895705e-05)
else:
x.append(500)
y.append(60.14113)
x.append(1000)
y.append(4.42547616)
x.append(1500)
y.append(0.7571358)
x.append(1800)
y.append(0.310112064)
x.append(1850)
y.append(0.26617037)
x.append(1900)
y.append(0.235984552)
x.append(1950)
y.append(0.20858418)
x.append(2000)
y.append(0.1818336)
x.append(2050)
y.append(0.156303936)
x.append(2100)
y.append(0.135815946)
x.append(2150)
y.append(0.121970058)
x.append(2200)
y.append(0.106092126)
x.append(2250)
y.append(0.093963767)
x.append(2300)
y.append(0.080591)
x.append(2350)
y.append(0.071953414)
x.append(2400)
y.append(0.064723741)
x.append(2450)
y.append(0.056408674)
x.append(2500)
y.append(0.049592333)
x.append(2550)
y.append(0.045233174)
x.append(2600)
y.append(0.039405779)
x.append(2650)
y.append(0.034789025)
x.append(2700)
y.append(0.031481101)
x.append(2750)
y.append(0.028296221)
x.append(2800)
y.append(0.024758081)
x.append(3000)
y.append(0.01552742)
x.append(3500)
y.append(0.004936807)
x.append(4000)
y.append(0.001680786)
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def dijetxs(color):
x = array('d', [])
y = array('d', [])
av = False
if av:
x.append(500)
y.append(50.11760833)
x.append(1000)
y.append(3.6878968)
x.append(1500)
y.append(0.6309465)
x.append(1800)
y.append(0.25842672)
x.append(1850)
y.append(0.221808642)
x.append(1900)
y.append(0.196653793)
x.append(1950)
y.append(0.17382015)
x.append(2000)
y.append(0.151528)
x.append(2050)
y.append(0.13025328)
x.append(2100)
y.append(0.113179955)
x.append(2150)
y.append(0.101641715)
x.append(2200)
y.append(0.088410105)
x.append(2250)
y.append(0.078303139)
x.append(2300)
y.append(0.067159167)
x.append(2350)
y.append(0.059961178)
x.append(2400)
y.append(0.053936451)
x.append(2450)
y.append(0.047007228)
x.append(2500)
y.append(0.041326944)
x.append(2550)
y.append(0.037694312)
x.append(2600)
y.append(0.03283815)
x.append(2650)
y.append(0.028990854)
x.append(2700)
y.append(0.026234251)
x.append(2750)
y.append(0.023580184)
x.append(2800)
y.append(0.020631734)
x.append(3000)
y.append(0.012939516)
x.append(3500)
y.append(0.004114006)
x.append(4000)
y.append(0.001400655)
else:
x.append(500)
y.append(54.39976667)
x.append(1000)
y.append(3.79482905)
x.append(1500)
y.append(0.64184835)
x.append(1800)
y.append(0.256497348)
x.append(1850)
y.append(0.226374873)
x.append(1900)
y.append(0.197961717)
x.append(1950)
y.append(0.169656743)
x.append(2000)
y.append(0.15153154)
x.append(2050)
y.append(0.132613545)
x.append(2100)
y.append(0.115242345)
x.append(2150)
y.append(0.099989991)
x.append(2200)
y.append(0.090241305)
x.append(2250)
y.append(0.077520663)
x.append(2300)
y.append(0.06822651)
x.append(2350)
y.append(0.061410158)
x.append(2400)
y.append(0.05406678)
x.append(2450)
y.append(0.047535109)
x.append(2500)
y.append(0.041757683)
x.append(2550)
y.append(0.036590855)
x.append(2600)
y.append(0.033034138)
x.append(2650)
y.append(0.029252846)
x.append(2700)
y.append(0.025036183)
x.append(2750)
y.append(0.022652202)
x.append(2800)
y.append(0.020569049)
x.append(3000)
y.append(0.012877625)
x.append(3500)
y.append(0.004038073)
x.append(4000)
y.append(0.001372723)
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def getCont(iGraph, iVal, color, label):
lContours = iGraph.GetContourList(iVal)
xnew = array('d', [])
ynew = array('d', [])
for i0 in range(0, lContours.GetSize()):
pCurv = lContours.At(i0)
#if i0 > 0:
#break
for i2 in range(0, pCurv.GetN()):
lY = pCurv.GetY()[i2]
if lY < 10:
lY = -100
# if pCurv.GetX()[i2]<1674.4685364 and pCurv.GetX()[i2]>1674.4685363:
# continue
if pCurv.GetX()[i2] < 1625 and pCurv.GetX(
)[i2] > 1599.9 and lY < 671 and lY > 10:
continue
if pCurv.GetX()[i2] < 2500 and pCurv.GetX(
)[i2] > 1600 and lY > 0 and lY < 3001:
continue
if 1600 == pCurv.GetX(
)[i2] and 759.23182239 < lY and lY < 759.23182240 or (lY < 725 and
lY > 724):
continue
xnew.append(pCurv.GetX()[i2])
ynew.append(lY)
if pCurv.GetX()[i2] > 1638.6097358 and pCurv.GetX(
)[i2] < 1638.6097359 and lY == 3000:
ynew.append(lY)
xnew.append(-10)
xnew.append(-10)
ynew.append(-100)
if pCurv.GetX(
)[i2] == 1600 and lY < 789.53492852 and lY > 789.53492851:
#1600.0 789.534928514
xnew.append(pCurv.GetX()[i2])
ynew.append(-100)
xnew.append(xnew[len(xnew) - 1])
ynew.append(ynew[len(ynew) - 1] + 100000)
xnew.append(iGraph.GetX()[0])
ynew.append(ynew[len(ynew) - 1])
xnew.append(iGraph.GetX()[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(-100)
xnew.append(xnew[0])
ynew.append(ynew[0])
# for i0 in range(0,len(ynew)):
# print(ynew[i0])
# for i in range(len(xnew)):
#print xnew[i]#,':',ynew[i]
# if xnew[i]<2000 and ynew[i]<1000 and xnew[i]>0:
# ynew[i]=-1
content = ''
lGraph = makeAGraph(xnew, ynew, color, 1)
for i in range(len(xnew)):
# if xnew[i]>2000 and ynew[i]<3000:
print str(xnew[i]) + '\t' + str(ynew[i])
content += str(xnew[i]) + '\t' + str(ynew[i]) + '\n'
fout = open(label + '.txt', 'w+')
fout.write(content)
fout.close()
lGraph.GetXaxis().SetTitle("m_{med} (GeV)")
lGraph.GetYaxis().SetTitle("m_{dm} (GeV)")
lGraph.SetLineWidth(4)
return lGraph
def dijetexp(color, i90CL=False):
if i90CL:
lPFile = ROOT.TFile("PF90CL.root")
lPGraph = lPFile.Get("exp")
lCFile = ROOT.TFile("Calo90CL.root")
lCGraph = lCFile.Get("exp")
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
lPFile = ROOT.TFile("limits_freq_qq_calodijet2016_pfdijet2016.root")
lPFile1 = ROOT.TFile(
'limits_freq_gg_qg_qq_pfdijet2016_95CL_moriondDataset_JECV1_35900pb.root'
)
lPFile2 = ROOT.TFile('limits_freq_qq_calodijet2016.root')
lCGraph = lPFile2.Get("exp1sigma_qq_calodijet2016")
lPGraph = lPFile1.Get("exp1sigma_qq_pfdijet2016")
pf = ROOT.TFile(
'limits_freq_gg_qg_qq_pfdijet2016_95CL_moriondDataset_JECV1_35900pb.root'
)
calo = ROOT.TFile('limits_freq_qq_calodijet2016.root')
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(8.619535664753032)
x.append(601.9305019305019)
y.append(calo.Get('exp_qq_calodijet2016').Eval(601.9305019305019))
x.append(650.1930501930501)
y.append(calo.Get('exp_qq_calodijet2016').Eval(650.1930501930501))
x.append(698.4555984555984)
y.append(calo.Get('exp_qq_calodijet2016').Eval(698.4555984555984))
x.append(746.7181467181467)
y.append(calo.Get('exp_qq_calodijet2016').Eval(746.7181467181467))
x.append(812.3552123552123)
y.append(calo.Get('exp_qq_calodijet2016').Eval(812.3552123552123))
x.append(879.9227799227799)
y.append(calo.Get('exp_qq_calodijet2016').Eval(879.9227799227799))
x.append(953.2818532818533)
y.append(calo.Get('exp_qq_calodijet2016').Eval(953.2818532818533))
x.append(1032.4324324324325)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1032.4324324324325))
x.append(1117.3745173745174)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1117.3745173745174))
x.append(1204.2471042471043)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1204.2471042471043))
x.append(1291.1196911196912)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1291.1196911196912))
x.append(1377.992277992278)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1377.992277992278))
x.append(1447.4903474903476)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1447.4903474903476))
x.append(1526.6409266409269)
y.append(calo.Get('exp_qq_calodijet2016').Eval(1526.6409266409269))
x.append(1835.5653097740994)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(1835.5653097740994))
x.append(2093.5706031462005)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(2093.5706031462005))
x.append(2341.333967046895)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(2341.333967046895))
x.append(2578.7249310370526)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(2578.7249310370526))
x.append(2754.175054052038)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(2754.175054052038))
x.append(3043.297547155743)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(3043.297547155743))
x.append(3435.7200104376343)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(3435.7200104376343))
x.append(3797.188362036829)
y.append(pf.Get('exp_qq_pfdijet2016').Eval(3797.188362036829))
#x.append(4138.0097666443)
#y.append(0.010394256012614915)
#x.append(4437.53727726832)
#y.append(0.007783059964888597)
#x.append(4809.37802877805)
#y.append(0.005500486503683233)
#x.append(5201.8657272795035)
#y.append(0.0037763370738787476)
#x.append(5656.294266756131)
#y.append(0.002376824781618676)
#x.append(6162.454335346305)
#y.append(0.0015404057465133833)
#x.append(6658.307239245508)
#y.append(0.0010277152713018512)
#x.append(7081.879519868783)
#y.append(0.0007478854520557367)
#x.append(7515.922053231939)
#y.append(0.0006113493020261009)
#x.append(7826.148139864309)
#y.append(0.0006301833652603328)
#x.append(7971.035562513977)
#y.append(0.0007082739961764732)
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def dijetobs(color,i90CL=False):
if i90CL:
lPFile = ROOT.TFile("PF90CL.root")
lPGraph = lPFile.Get("obs")
lCFile = ROOT.TFile("Calo90CL.root")
lCGraph = lCFile.Get("obs")
lGraph = merge(lPGraph,lCGraph,color,1400,3900)
return lGraph
lPFile = ROOT.TFile("limits_freq_qq_calodijet2016_pfdijet2016.root")
lCGraph = lPFile.Get("obs_qq_calodijet2016")
lPGraph = lPFile.Get("obs_qq_pfdijet2016")
lGraph = merge(lPGraph,lCGraph,color,1400,3900)
return lGraph
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(4.417344703140073)
x.append(598.0694980694981)
y.append(4.417344703140073)
x.append(630.8880308880309)
y.append(4.256300264549507)
x.append(654.0540540540541)
y.append(3.8075460212223757)
x.append(694.5945945945946)
y.append(3.668733193018217)
x.append(729.3436293436293)
y.append(3.8075460212223757)
x.append(764.0926640926641)
y.append(4.101127070551304)
x.append(802.7027027027027)
y.append(4.101127070551304)
x.append(850.965250965251)
y.append(4.101127070551304)
x.append(901.1583011583012)
y.append(3.4061052563594854)
x.append(949.4208494208494)
y.append(2.181271260987366)
x.append(1001.5444015444016)
y.append(1.2040517722974518)
x.append(1053.6679536679537)
y.append(0.742963950759495)
x.append(1103.861003861004)
y.append(0.6404004271197296)
x.append(1152.1235521235521)
y.append(0.6404004271197296)
x.append(1202.3166023166023)
y.append(0.6404004271197296)
x.append(1250.5791505791508)
y.append(0.5318711718664566)
x.append(1300.772200772201)
y.append(0.395161107144985)
x.append(1347.1042471042472)
y.append(0.29359045735093353)
x.append(1406.949806949807)
y.append(0.2530613417967531)
x.append(1451.3513513513512)
y.append(0.21812712609873705)
x.append(1497.6833976833977)
y.append(0.19512934226359663)
x.append(1545.9459459459463)
y.append(0.1880154542858142)
x.append(1601.930501930502)
y.append(0.22638034095214513)
x.append(1773.4939983598003)
y.append(0.19363285556162788)
x.append(1856.1470215462607)
y.append(0.18276433417000734)
x.append(1949.1398270334746)
y.append(0.17251326138790926)
x.append(2104.399649593677)
y.append(0.19961530735946453)
x.append(2228.313949153806)
y.append(0.17271327873762074)
x.append(2279.686684559755)
y.append(0.129191924781803)
x.append(2341.268731827331)
y.append(0.08603798404788117)
x.append(2402.91601431447)
y.append(0.06072699537075284)
x.append(2506.118131663312)
y.append(0.05103381655253201)
x.append(2537.13747856557)
y.append(0.05104038763675874)
x.append(2640.6657720122266)
y.append(0.05735532852726994)
x.append(2764.9388652799516)
y.append(0.0683136621065803)
x.append(2878.676470588234)
y.append(0.06834591974293597)
x.append(2971.6040408558847)
y.append(0.060870506941983885)
x.append(3074.8061582047267)
y.append(0.051154420958441714)
x.append(3167.440169984343)
y.append(0.03507618092268031)
x.append(3229.2179229106086)
y.append(0.0278084825312735)
x.append(3311.773093267725)
y.append(0.024056620832096302)
x.append(3425.412845746663)
y.append(0.022058937089766098)
x.append(3507.9027808842166)
y.append(0.018005500745490594)
x.append(3631.8170804443444)
y.append(0.01557891080700686)
x.append(3797.188362036829)
y.append(0.014709523294944023)
x.append(3880.004473272198)
y.append(0.01605472139402788)
#x.append(4014.291172742861)
#y.append(0.014301172804579109)
#x.append(4096.781107880414)
#y.append(0.011673263156169997)
#x.append(4199.754901960783)
#y.append(0.008004607532669124)
#x.append(4292.519384179526)
#y.append(0.006165135029541548)
#x.append(4395.656266308803)
#y.append(0.004888573359511099)
#x.append(4509.328636397525)
#y.append(0.004614774148116912)
#x.append(4674.797770819354)
#y.append(0.004754085899371195)
#x.append(4778.293446656227)
#y.append(0.005189296395244354)
#x.append(4902.240363826138)
#y.append(0.004622306281508651)
#x.append(5036.527063296801)
#y.append(0.004117443041530253)
#x.append(5139.663945426078)
#y.append(0.003264879400966817)
#x.append(5273.91802728696)
#y.append(0.0028249949084984487)
#x.append(5418.740214717065)
#y.append(0.002995817276430047)
#x.append(5563.56240214717)
#y.append(0.0031769689661235974)
#x.append(5739.240848430627)
#y.append(0.0029139008040259193)
#x.append(5945.840788786996)
#y.append(0.002449838209903935)
#x.append(6121.388764631328)
#y.append(0.0020004411218901654)
#x.append(6266.015246402742)
#y.append(0.0017820229509131172)
#x.append(6266.015246402742)
#y.append(0.0017820229509131172)
#x.append(6431.256057556102)
#y.append(0.0014979640103628745)
#x.append(6431.256057556102)
#y.append(0.0014979640103628745)
#x.append(6596.464251099679)
#y.append(0.0012231257522893248)
#x.append(6782.286774025198)
#y.append(0.0009424133022435063)
#x.append(6906.201073585326)
#y.append(0.0008154048580218117)
#x.append(7092.219302169536)
#y.append(0.000747917549775239)
#x.append(7236.943636770297)
#y.append(0.0007269362219934222)
#x.append(7236.943636770297)
#y.append(0.0007269362219934222)
#x.append(7402.380153582345)
#y.append(0.0007274355606432079)
#x.append(7516.117758890628)
#y.append(0.0007277790549174683)
#x.append(7702.266457913961)
#y.append(0.000749813755546937)
#x.append(7867.800827555356)
#y.append(0.0008186658745549854)
#x.append(8002.446320733617)
#y.append(0.0010038700029024285)
lGraph = makeAGraph( x, y, color, 3 )
lGraph.SetMarkerColor(color)
return lGraph
def dijetexp(color, i90CL=False):
if i90CL:
lPFile = ROOT.TFile("PF90CL.root")
lPGraph = lPFile.Get("exp")
lCFile = ROOT.TFile("Calo90CL.root")
lCGraph = lCFile.Get("exp")
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
lPFile = ROOT.TFile("limits_freq_qq_calodijet2016_pfdijet2016.root")
lCGraph = lPFile.Get("exp1sigma_qq_calodijet2016")
lPGraph = lPFile.Get("exp1sigma_qq_pfdijet2016")
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(8.619535664753032)
x.append(601.9305019305019)
y.append(8.619535664753032)
x.append(650.1930501930501)
y.append(6.170531598835212)
x.append(698.4555984555984)
y.append(3.95161107144985)
x.append(746.7181467181467)
y.append(2.6263635276533353)
x.append(812.3552123552123)
y.append(2.0251239642430052)
x.append(879.9227799227799)
y.append(1.811609194200417)
x.append(953.2818532818533)
y.append(1.5615230060004996)
x.append(1032.4324324324325)
y.append(1.117859177755407)
x.append(1117.3745173745174)
y.append(0.9284145445194764)
x.append(1204.2471042471043)
y.append(0.7158774948942709)
x.append(1291.1196911196912)
y.append(0.5318711718664566)
x.append(1377.992277992278)
y.append(0.4757944314009414)
x.append(1447.4903474903476)
y.append(0.41011270705513125)
x.append(1526.6409266409269)
y.append(0.3807546021222384)
x.append(1835.5653097740994)
y.append(0.19939269222685885)
x.append(2093.5706031462005)
y.append(0.12909216168780502)
x.append(2341.333967046895)
y.append(0.0911857420299053)
x.append(2578.7249310370526)
y.append(0.0625629193564747)
x.append(2754.175054052038)
y.append(0.0468220376016705)
x.append(3043.297547155743)
y.append(0.033078966118145904)
x.append(3435.7200104376343)
y.append(0.021428159256036205)
x.append(3797.188362036829)
y.append(0.014709523294944023)
#x.append(4138.0097666443)
#y.append(0.010394256012614915)
#x.append(4437.53727726832)
#y.append(0.007783059964888597)
#x.append(4809.37802877805)
#y.append(0.005500486503683233)
#x.append(5201.8657272795035)
#y.append(0.0037763370738787476)
#x.append(5656.294266756131)
#y.append(0.002376824781618676)
#x.append(6162.454335346305)
#y.append(0.0015404057465133833)
#x.append(6658.307239245508)
#y.append(0.0010277152713018512)
#x.append(7081.879519868783)
#y.append(0.0007478854520557367)
#x.append(7515.922053231939)
#y.append(0.0006113493020261009)
#x.append(7826.148139864309)
#y.append(0.0006301833652603328)
#x.append(7971.035562513977)
#y.append(0.0007082739961764732)
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def dijetobs(color, i90CL=False):
if i90CL:
lPFile = ROOT.TFile("PF90CL.root")
lPGraph = lPFile.Get("obs")
lCFile = ROOT.TFile("Calo90CL.root")
lCGraph = lCFile.Get("obs")
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
lPFile = ROOT.TFile("limits_freq_qq_calodijet2016_pfdijet2016.root")
lCGraph = lPFile.Get("obs_qq_calodijet2016")
lPGraph = lPFile.Get("obs_qq_pfdijet2016")
lGraph = merge(lPGraph, lCGraph, color, 1400, 3900)
return lGraph
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(4.417344703140073)
x.append(598.0694980694981)
y.append(4.417344703140073)
x.append(630.8880308880309)
y.append(4.256300264549507)
x.append(654.0540540540541)
y.append(3.8075460212223757)
x.append(694.5945945945946)
y.append(3.668733193018217)
x.append(729.3436293436293)
y.append(3.8075460212223757)
x.append(764.0926640926641)
y.append(4.101127070551304)
x.append(802.7027027027027)
y.append(4.101127070551304)
x.append(850.965250965251)
y.append(4.101127070551304)
x.append(901.1583011583012)
y.append(3.4061052563594854)
x.append(949.4208494208494)
y.append(2.181271260987366)
x.append(1001.5444015444016)
y.append(1.2040517722974518)
x.append(1053.6679536679537)
y.append(0.742963950759495)
x.append(1103.861003861004)
y.append(0.6404004271197296)
x.append(1152.1235521235521)
y.append(0.6404004271197296)
x.append(1202.3166023166023)
y.append(0.6404004271197296)
x.append(1250.5791505791508)
y.append(0.5318711718664566)
x.append(1300.772200772201)
y.append(0.395161107144985)
x.append(1347.1042471042472)
y.append(0.29359045735093353)
x.append(1406.949806949807)
y.append(0.2530613417967531)
x.append(1451.3513513513512)
y.append(0.21812712609873705)
x.append(1497.6833976833977)
y.append(0.19512934226359663)
x.append(1545.9459459459463)
y.append(0.1880154542858142)
x.append(1601.930501930502)
y.append(0.22638034095214513)
x.append(1773.4939983598003)
y.append(0.19363285556162788)
x.append(1856.1470215462607)
y.append(0.18276433417000734)
x.append(1949.1398270334746)
y.append(0.17251326138790926)
x.append(2104.399649593677)
y.append(0.19961530735946453)
x.append(2228.313949153806)
y.append(0.17271327873762074)
x.append(2279.686684559755)
y.append(0.129191924781803)
x.append(2341.268731827331)
y.append(0.08603798404788117)
x.append(2402.91601431447)
y.append(0.06072699537075284)
x.append(2506.118131663312)
y.append(0.05103381655253201)
x.append(2537.13747856557)
y.append(0.05104038763675874)
x.append(2640.6657720122266)
y.append(0.05735532852726994)
x.append(2764.9388652799516)
y.append(0.0683136621065803)
x.append(2878.676470588234)
y.append(0.06834591974293597)
x.append(2971.6040408558847)
y.append(0.060870506941983885)
x.append(3074.8061582047267)
y.append(0.051154420958441714)
x.append(3167.440169984343)
y.append(0.03507618092268031)
x.append(3229.2179229106086)
y.append(0.0278084825312735)
x.append(3311.773093267725)
y.append(0.024056620832096302)
x.append(3425.412845746663)
y.append(0.022058937089766098)
x.append(3507.9027808842166)
y.append(0.018005500745490594)
x.append(3631.8170804443444)
y.append(0.01557891080700686)
x.append(3797.188362036829)
y.append(0.014709523294944023)
x.append(3880.004473272198)
y.append(0.01605472139402788)
#x.append(4014.291172742861)
#y.append(0.014301172804579109)
#x.append(4096.781107880414)
#y.append(0.011673263156169997)
#x.append(4199.754901960783)
#y.append(0.008004607532669124)
#x.append(4292.519384179526)
#y.append(0.006165135029541548)
#x.append(4395.656266308803)
#y.append(0.004888573359511099)
#x.append(4509.328636397525)
#y.append(0.004614774148116912)
#x.append(4674.797770819354)
#y.append(0.004754085899371195)
#x.append(4778.293446656227)
#y.append(0.005189296395244354)
#x.append(4902.240363826138)
#y.append(0.004622306281508651)
#x.append(5036.527063296801)
#y.append(0.004117443041530253)
#x.append(5139.663945426078)
#y.append(0.003264879400966817)
#x.append(5273.91802728696)
#y.append(0.0028249949084984487)
#x.append(5418.740214717065)
#y.append(0.002995817276430047)
#x.append(5563.56240214717)
#y.append(0.0031769689661235974)
#x.append(5739.240848430627)
#y.append(0.0029139008040259193)
#x.append(5945.840788786996)
#y.append(0.002449838209903935)
#x.append(6121.388764631328)
#y.append(0.0020004411218901654)
#x.append(6266.015246402742)
#y.append(0.0017820229509131172)
#x.append(6266.015246402742)
#y.append(0.0017820229509131172)
#x.append(6431.256057556102)
#y.append(0.0014979640103628745)
#x.append(6431.256057556102)
#y.append(0.0014979640103628745)
#x.append(6596.464251099679)
#y.append(0.0012231257522893248)
#x.append(6782.286774025198)
#y.append(0.0009424133022435063)
#x.append(6906.201073585326)
#y.append(0.0008154048580218117)
#x.append(7092.219302169536)
#y.append(0.000747917549775239)
#x.append(7236.943636770297)
#y.append(0.0007269362219934222)
#x.append(7236.943636770297)
#y.append(0.0007269362219934222)
#x.append(7402.380153582345)
#y.append(0.0007274355606432079)
#x.append(7516.117758890628)
#y.append(0.0007277790549174683)
#x.append(7702.266457913961)
#y.append(0.000749813755546937)
#x.append(7867.800827555356)
#y.append(0.0008186658745549854)
#x.append(8002.446320733617)
#y.append(0.0010038700029024285)
lGraph = makeAGraph(x, y, color, 3)
lGraph.SetMarkerColor(color)
return lGraph
def dijetexp(color,i90CL=False):
if i90CL:
lPFile = ROOT.TFile("PF90CL.root")
lPGraph = lPFile.Get("exp")
lCFile = ROOT.TFile("Calo90CL.root")
lCGraph = lCFile.Get("exp")
lGraph = merge(lPGraph,lCGraph,color,1400,3900)
return lGraph
lPFile = ROOT.TFile("limits_freq_qq_calodijet2016_pfdijet2016.root")
lCGraph = lPFile.Get("exp1sigma_qq_calodijet2016")
lPGraph = lPFile.Get("exp1sigma_qq_pfdijet2016")
lGraph = merge(lPGraph,lCGraph,color,1400,3900)
return lGraph
x = array('d', [])
y = array('d', [])
#x.append(500)
#y.append(8.619535664753032)
x.append(601.9305019305019)
y.append(8.619535664753032)
x.append(650.1930501930501)
y.append(6.170531598835212)
x.append(698.4555984555984)
y.append(3.95161107144985)
x.append(746.7181467181467)
y.append(2.6263635276533353)
x.append(812.3552123552123)
y.append(2.0251239642430052)
x.append(879.9227799227799)
y.append(1.811609194200417)
x.append(953.2818532818533)
y.append(1.5615230060004996)
x.append(1032.4324324324325)
y.append(1.117859177755407)
x.append(1117.3745173745174)
y.append(0.9284145445194764)
x.append(1204.2471042471043)
y.append(0.7158774948942709)
x.append(1291.1196911196912)
y.append(0.5318711718664566)
x.append(1377.992277992278)
y.append(0.4757944314009414)
x.append(1447.4903474903476)
y.append(0.41011270705513125)
x.append(1526.6409266409269)
y.append(0.3807546021222384)
x.append(1835.5653097740994)
y.append(0.19939269222685885)
x.append(2093.5706031462005)
y.append(0.12909216168780502)
x.append(2341.333967046895)
y.append(0.0911857420299053)
x.append(2578.7249310370526)
y.append(0.0625629193564747)
x.append(2754.175054052038)
y.append(0.0468220376016705)
x.append(3043.297547155743)
y.append(0.033078966118145904)
x.append(3435.7200104376343)
y.append(0.021428159256036205)
x.append(3797.188362036829)
y.append(0.014709523294944023)
#x.append(4138.0097666443)
#y.append(0.010394256012614915)
#x.append(4437.53727726832)
#y.append(0.007783059964888597)
#x.append(4809.37802877805)
#y.append(0.005500486503683233)
#x.append(5201.8657272795035)
#y.append(0.0037763370738787476)
#x.append(5656.294266756131)
#y.append(0.002376824781618676)
#x.append(6162.454335346305)
#y.append(0.0015404057465133833)
#x.append(6658.307239245508)
#y.append(0.0010277152713018512)
#x.append(7081.879519868783)
#y.append(0.0007478854520557367)
#x.append(7515.922053231939)
#y.append(0.0006113493020261009)
#x.append(7826.148139864309)
#y.append(0.0006301833652603328)
#x.append(7971.035562513977)
#y.append(0.0007082739961764732)
lGraph = makeAGraph( x, y, color, 3 )
lGraph.SetMarkerColor(color)
return lGraph
