file /home/anarendran/Documents/temp/rivet/pyext/build/lib.linux-x86_64-3.10/rivet/hepdatapatches/DELPHI_2003_I620250.py
/home/anarendran/Documents/temp/rivet/pyext/build/lib.linux-x86_64-3.10/rivet/hepdatapatches/DELPHI_2003_I620250.py
Namespaces
Name |
---|
rivet |
rivet::hepdatapatches |
rivet::hepdatapatches::DELPHI_2003_I620250 |
Source code
def patch(path, ao):
# fix bin widths
if "DELPHI_2003_I620250" in path :
if "d30" in path or "d31" in path:
step=0.25
for p in ao.points():
p.setXErrs(step)
ao.points()[-1].setXErrs(0.5)
elif "d32" in path or "d33" in path:
step=0.2
for p in ao.points() : p.setXErrs(step)
elif "d34" in path or "d35" in path:
for p in ao.points() :
if(p.x()>.5 and p.x()<1.) :
p.setXErrs(0.125)
elif(p.x()>.9 and p.x()<1.1) :
p.setXErrs(0.1)
elif(p.x()>1.1 and p.x()<1.3) :
p.setXErrs(0.15)
elif(p.x()>.2 and p.x()<.3) :
p.setXErrs(0.15)
elif(p.x()<.2) :
p.setXErrs((p.x(),.1-p.x()))
elif(p.x()>2 and p.x()<4) :
p.setXErrs(0.5)
elif(p.x()>4 and p.x()<8) :
p.setXErrs(1)
elif(p.x()>1.6 and p.x()<1.7) :
p.setXErrs((p.x()-1.4,2.-p.x()))
else :
p.setXErrs(2)
elif "d36" in path or "d37" in path:
for p in ao.points() :
if(p.x()<0.6) :
p.setXErrs(0.1)
elif(p.x()<0.8) :
p.setXErrs(0.125)
elif(p.x()<1.1) :
p.setXErrs(0.175)
elif(p.x()<1.9) :
p.setXErrs(0.2)
else :
p.setXErrs(0.5)
elif "d38" in path or "d39" in path:
for p in ao.points() :
if(p.x()<0.1) :
p.setXErrs(0.005)
elif (p.x()<0.2) :
p.setXErrs(0.01)
else :
p.setXErrs(0.02)
elif ("d40" in path or "d41" in path or
"d42" in path or "d43" in path):
for p in ao.points() :
if(p.x()<0.04) :
p.setXErrs(0.01)
elif (p.x()<0.08) :
p.setXErrs(0.005)
elif (p.x()<0.17) :
p.setXErrs(0.01)
else :
p.setXErrs(0.02)
elif "d44" in path or "d45" in path:
for p in ao.points() :
if(p.x()<0.2) :
p.setXErrs(0.01)
else :
p.setXErrs(0.02)
elif "d46" in path or "d47" in path:
for p in ao.points() :
if(p.x()<0.08) :
p.setXErrs(0.005)
elif (p.x()<0.15) :
p.setXErrs(0.01)
elif (p.x()<0.2) :
p.setXErrs(0.015)
else :
p.setXErrs(0.02)
elif "d48" in path or "d49" in path:
for p in ao.points() :
if(p.x()<0.11) :
p.setXErrs(0.005)
elif(p.x()<0.21) :
p.setXErrs(0.01)
else :
p.setXErrs(0.015)
elif "d50" in path or "d51" in path:
for p in ao.points() :
if(p.x()<0.10) :
p.setXErrs(0.005)
elif(p.x()<0.21) :
p.setXErrs(0.01)
else :
p.setXErrs(0.02)
elif "d52" in path or "d53" in path:
for p in ao.points() :
p.setXErrs(0.02)
elif "d54" in path or "d55" in path:
for p in ao.points() :
if(p.x()<0.16) :
p.setXErrs(0.01)
elif(p.x()<0.5) :
p.setXErrs(0.02)
else :
p.setXErrs(0.03)
elif ("d56" in path or "d57" in path or "d58" in path or
"d59" in path or "d60" in path or "d61" in path or
"d62" in path or "d63" in path ):
for p in ao.points() :
if(p.x()<0.06) :
p.setXErrs(0.005)
elif(p.x()<0.16) :
p.setXErrs(0.01)
else :
p.setXErrs(0.02)
elif "d64" in path or "d65" in path:
for p in ao.points() :
if(p.x()<0.04) :
p.setXErrs(0.005)
elif(p.x()<0.08) :
p.setXErrs(0.01)
elif(p.x()<0.2) :
p.setXErrs(0.02)
else :
p.setXErrs(0.025)
elif ("d66" in path or "d67" in path or
"d68" in path or "d69" in path ):
for p in ao.points() :
if(p.x()<0.06) :
p.setXErrs(0.005)
elif(p.x()<0.12) :
p.setXErrs(0.01)
elif(p.x()<0.2) :
p.setXErrs(0.02)
elif(p.x()<0.41) :
p.setXErrs(0.025)
else :
p.setXErrs(0.05)
elif "d70" in path or "d71" in path:
for p in ao.points() :
if(p.x()>.04 and p.x()<.12) :
p.setXErrs(0.01)
elif(p.x()>.12) :
p.setXErrs(0.02)
elif(p.x()>0.02) :
p.setXErrs(0.005)
elif(p.x()>0.017) :
p.setXErrs(0.002)
elif(p.x()>0.006) :
p.setXErrs(0.003)
else :
p.setXErrs(0.002)
return ao
Updated on 2022-08-07 at 20:46:08 +0100