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 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