Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F8724495
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
135 KB
Subscribers
None
View Options
diff --git a/PDF/PomeronPDF.cc b/PDF/PomeronPDF.cc
--- a/PDF/PomeronPDF.cc
+++ b/PDF/PomeronPDF.cc
@@ -1,452 +1,451 @@
// -*- C++ -*-
//
// DifractivePDF.cc is a part of Herwig++ - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#include "PomeronPDF.h"
#include <ThePEG/PDT/ParticleData.h>
#include <ThePEG/PDT/EnumParticles.h>
#include <ThePEG/Persistency/PersistentOStream.h>
#include <ThePEG/Persistency/PersistentIStream.h>
#include <ThePEG/Repository/EventGenerator.h>
#include <ThePEG/Interface/ClassDocumentation.h>
#include <ThePEG/Interface/Parameter.h>
#include <ThePEG/Interface/Switch.h>
#include <istream>
#include <iostream>
#include <string>
using namespace std;
using namespace ThePEG;
using namespace Herwig;
const int PomeronPDF::nPDFFlavour_ = 3;
PomeronPDF::PomeronPDF() :
pdfTable_(nPDFFlavour_, vector<vector<double> >() ),
lxGrid_(nPDFFlavour_, vector<double>() ),
lqqGrid_ (nPDFFlavour_, vector<double>() ),
fileName_(3),
nxPoints_(100),
nqPoints_(88),
PDFFit_(0),
boundary_(0)
{}
bool PomeronPDF::canHandleParticle(tcPDPtr particle) const {
return ( abs(particle->id()) == ParticleID::pomeron );
}
cPDVector PomeronPDF::partons(tcPDPtr p) const {
// Return the parton types which are described by these parton
// densities.
cPDVector ret;
if ( canHandleParticle(p) ) {
ret.push_back(getParticleData( ParticleID::g));
if(PDFFit_==0) {
ret.push_back(getParticleData( ParticleID::c));
ret.push_back(getParticleData( ParticleID::cbar));
}
ret.push_back(getParticleData( ParticleID::d));
ret.push_back(getParticleData( ParticleID::dbar));
ret.push_back(getParticleData( ParticleID::u));
ret.push_back(getParticleData( ParticleID::ubar));
ret.push_back(getParticleData( ParticleID::s));
ret.push_back(getParticleData( ParticleID::sbar));
}
return ret;
}
#ifdef NDEBUG
double PomeronPDF::xfx(tcPDPtr, tcPDPtr parton, Energy2 qq,
double x, double, Energy2) const {
#else
double PomeronPDF::xfx(tcPDPtr particle, tcPDPtr parton, Energy2 qq,
double x, double, Energy2) const {
#endif
// assert particle is pomeron
assert(particle->id()==ParticleID::pomeron);
switch(parton->id()) {
case ParticleID::g:
return getPDFValue(gluon , x, qq);
case ParticleID::u: case ParticleID::ubar:
case ParticleID::d: case ParticleID::dbar:
case ParticleID::s: case ParticleID::sbar:
return PDFFit_==0 ? getPDFValue(singlet, x, qq)/6. : getPDFValue(singlet, x, qq);
case ParticleID::c: case ParticleID::cbar:
return PDFFit_==0 ? getPDFValue(charm , x, qq)*9./8. : 0.;
case ParticleID::b: case ParticleID::bbar:
return 0;
default:
- assert(false);
return 0.;
}
}
#ifdef NDEBUG
double PomeronPDF::xfvx(tcPDPtr, tcPDPtr parton, Energy2 qq,
double x, double, Energy2) const {
#else
double PomeronPDF::xfvx(tcPDPtr particle, tcPDPtr parton, Energy2 qq,
double x, double, Energy2) const {
#endif
// assert particle is pomeron
assert(particle->id()==ParticleID::pomeron);
// valence parton is just gluon
if(parton->id() == ParticleID::g){
return getPDFValue(gluon , x, qq);
}
else
return 0.;
}
double PomeronPDF::getPDFValue(PDFFlavour flPDF, double x, Energy2 qq) const {
// compute log(x) and log(q2)
// Fit 2007 contains l= logx. Fit 2006 contains l = logl/log10
double l, lmax, lmin, lqqmin, lqqmax;
if(!PDFFit_) {
l = log(x);
lmin = lxGrid_.at(flPDF).at(0);
lmax = lxGrid_.at(flPDF).at(nxPoints_- 1);
lqqmin = lqqGrid_.at(flPDF).at(0);
lqqmax = lqqGrid_.at(flPDF).at(nqPoints_ -1);
}
else {
l = log(x)/log(10);
lmin = lxGrid_.at(flPDF).at(0);
lmax = lxGrid_.at(flPDF).at(nxPoints_- 1);
lqqmin = lqqGrid_.at(flPDF).at(0);
lqqmax = lqqGrid_.at(flPDF).at(nqPoints_ -1);
}
double lqq = log(qq/GeV2);
// find the x location in the table
if(l <= lmin ) {
if( lqq <= lqqmin) {
double pdf00 = pdfTable_.at(flPDF).at(0).at(0);
double pdf10 = pdfTable_.at(flPDF).at(1).at(0);
if(boundary_==0){
double lqq1 = lqqGrid_.at(flPDF).at(1);
double lqq0 = lqqGrid_.at(flPDF).at(0);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
return pdf00 + qdq*(pdf10 - pdf00);
}
else
return pdf00;
}
else if( lqq >= lqqmax) {
double pdf10 = pdfTable_.at(flPDF).at(nqPoints_-1).at(0);
double pdf00 = pdfTable_.at(flPDF).at(nqPoints_-2).at(0);
if(boundary_==0){
double lqq0 = lqqGrid_.at(flPDF).at(nqPoints_-2);
double lqq1 = lqqGrid_.at(flPDF).at(nqPoints_-1);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq1)/dq;
return pdf10 + qdq*(pdf10 - pdf00);
}
else return pdf10;
} // (x< xmin)
unsigned int qIndx = 0;
for(; qIndx < lqqGrid_[flPDF].size(); ++qIndx )
if(lqq < lqqGrid_[flPDF][qIndx] ) break;
assert(qIndx > 0);
double pdf10 = pdfTable_.at(flPDF).at(qIndx).at(0);
double pdf00 = pdfTable_.at(flPDF).at(qIndx-1).at(0);
double lqq0 = lqqGrid_.at(flPDF).at(qIndx-1);
double lqq1 = lqqGrid_.at(flPDF).at(qIndx);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
return pdf00 + qdq*(pdf10 - pdf00);
}
else if (l >= lmax ) {
if( lqq <= lqqmin){
double pdf01 = pdfTable_.at(flPDF).at(0).at(nxPoints_-1);
double pdf11 = pdfTable_.at(flPDF).at(1).at(nxPoints_-1);
if(boundary_==0){
double lqq1 = lqqGrid_.at(flPDF).at(1);
double lqq0 = lqqGrid_.at(flPDF).at(0);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
return pdf01 + qdq*(pdf11 - pdf01);
}
else return pdf01;
}
else if( lqq >= lqqmax){
double pdf11 = pdfTable_.at(flPDF).at(nqPoints_-1).at(nxPoints_-1);
double pdf01 = pdfTable_.at(flPDF).at(nqPoints_-2).at(nxPoints_-1);
if(boundary_==0){
double lqq1 = lqqGrid_.at(flPDF).at(nqPoints_-1);
double lqq0 = lqqGrid_.at(flPDF).at(nqPoints_-2);
double dq = lqq1 - lqq0;
double qdq = ( lqq - lqq1 )/dq;
return pdf11 + qdq*( pdf11 - pdf01);
} else return pdf11;
}
unsigned int qIndx = 0;
for(; qIndx < lqqGrid_.at(flPDF).size(); ++qIndx )
if(lqq < lqqGrid_.at(flPDF).at(qIndx) ) break;
assert(qIndx > 0);
double pdf01 = pdfTable_.at(flPDF).at(qIndx-1).at(nxPoints_-1);
double pdf11 = pdfTable_.at(flPDF).at(qIndx ).at(nxPoints_-1);
double lqq0 = lqqGrid_.at(flPDF).at(qIndx-1);
double lqq1 = lqqGrid_.at(flPDF).at(qIndx);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
return pdf01 + qdq*(pdf11 - pdf01);
}
else {
unsigned int xIndx = 0;
for(; xIndx < lxGrid_.at(flPDF).size(); ++xIndx )
if (l < lxGrid_.at(flPDF).at(xIndx) ) break;
assert(xIndx > 0);
double lx0 = lxGrid_.at(flPDF).at(xIndx-1);
double lx1 = lxGrid_.at(flPDF).at(xIndx);
double dl = lx1 - lx0;
double ldl = (l - lx0)/dl;
if( lqq <= lqqmin) {
double pdf00 = pdfTable_.at(flPDF).at(0).at(xIndx-1);
double pdf10 = pdfTable_.at(flPDF).at(1).at(xIndx-1);
double pdf01 = pdfTable_.at(flPDF).at(0).at(xIndx);
double pdf11 = pdfTable_.at(flPDF).at(1).at(xIndx);
if(boundary_==0){
double lqq0 = lqqGrid_.at(flPDF).at(0);
double lqq1 = lqqGrid_.at(flPDF).at(1);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
return pdf00 + ldl*(pdf01 - pdf00) + qdq*(pdf10 - pdf00)
+ qdq*ldl*(pdf11 - pdf01 + pdf00 - pdf10);
}
else return pdf00 + ldl*(pdf01 - pdf00);
}
else if(lqq >= lqqmax){
double pdf00 = pdfTable_.at(flPDF).at(nqPoints_-2).at(xIndx-1);
double pdf10 = pdfTable_.at(flPDF).at(nqPoints_-1).at(xIndx-1);
double pdf01 = pdfTable_.at(flPDF).at(nqPoints_-2).at(xIndx);
double pdf11 = pdfTable_.at(flPDF).at(nqPoints_-1).at(xIndx);
if(boundary_==0){
double lqq0 = lqqGrid_.at(flPDF).at(nqPoints_-2);
double lqq1 = lqqGrid_.at(flPDF).at(nqPoints_-1);
double dq = lqq1 - lqq0 ;
double qdq = ( lqq - lqq0)/dq;
return pdf10 + ldl*(pdf11 - pdf10) + qdq*(pdf10 - pdf00)
+ qdq*ldl*(pdf11 - pdf01 + pdf00 - pdf10);
}
else return pdf10 + ldl*(pdf11 - pdf10);
}
else {
unsigned int qIndx = 0;
for(; qIndx < lqqGrid_.at(flPDF).size(); ++qIndx )
if(lqq < lqqGrid_.at(flPDF).at(qIndx) ) break;
double lqq0 = lqqGrid_.at(flPDF).at(qIndx-1);
double lqq1 = lqqGrid_.at(flPDF).at(qIndx);
double dq = lqq1 - lqq0;
double qdq = (lqq - lqq0)/dq;
double pdf00 = pdfTable_.at(flPDF).at(qIndx-1).at(xIndx-1);
double pdf01 = pdfTable_.at(flPDF).at(qIndx-1).at(xIndx);
double pdf10 = pdfTable_.at(flPDF).at(qIndx).at(xIndx-1);
double pdf11 = pdfTable_.at(flPDF).at(qIndx).at(xIndx);
return pdf00 + ldl*(pdf01 - pdf00) + qdq*(pdf10 - pdf00)
+ ldl*qdq*(pdf11 - pdf01 + pdf00 - pdf10);
}
}
throw Exception() << "Problem in PomeronPDF::getPDFValue() "
<< Exception::runerror;
return 0;
}
void PomeronPDF::loadTables() const {
double dtemp; // helper variable for reading
ifstream datafile;
for(int pdfIndx = 0; pdfIndx < nPDFFlavour_; ++pdfIndx ) {
// Set correct size for vectors
pdfTable_.at(pdfIndx).resize(nqPoints_);
lqqGrid_.at(pdfIndx).resize(nqPoints_);
lxGrid_.at(pdfIndx).resize(nxPoints_);
string name = rootName_ + fileName_.at(pdfIndx);
datafile.open(name.c_str(), ios::in );
if(!datafile.is_open())
throw Exception() << "Could not open file '" << name
<< "' in PomeronPDF::loadTables()"
<< Exception::runerror;
// read lx points
for (unsigned int xIndx = 0; xIndx < lxGrid_ .at(pdfIndx).size(); ++xIndx)
datafile >> lxGrid_.at(pdfIndx).at(xIndx);
// read qq points
for (unsigned int qIndx = 0; qIndx < lqqGrid_.at(pdfIndx).size(); ++qIndx) {
datafile >> dtemp;
lqqGrid_.at(pdfIndx).at(qIndx) = log(dtemp);
}
// read xf(lx, qq) table
for (unsigned int qIndx = 0; qIndx < pdfTable_.at(pdfIndx).size(); ++qIndx) {
pdfTable_.at(pdfIndx).at(qIndx).resize(nxPoints_);
for (unsigned int xIndx = 0; xIndx < pdfTable_.at(pdfIndx).at(qIndx).size();
++xIndx) {
datafile >> pdfTable_.at(pdfIndx).at(qIndx).at(xIndx);
if(datafile.eof())
throw Exception() << "Error while reading " << fileName_.at(pdfIndx)
<< " too few data points in file"
<< "in PomeronPDF::loadTables()"
<< Exception::runerror;
}
}
datafile >> dtemp;
if(!datafile.eof())
throw Exception() << "Error reading end of " << fileName_.at(pdfIndx)
<< " too many data points in file"
<< "in PomeronPDF::loadTables()"
<< Exception::runerror;
datafile.close();
}
}
void PomeronPDF::doinit() {
PDFBase::doinit();
switch(PDFFit_){
case 0:
fileName_.at(0) = "2007/h12007jetsdpdf_charm.data";
fileName_.at(1) = "2007/h12007jetsdpdf_gluon.data";
fileName_.at(2) = "2007/h12007jetsdpdf_singlet.data";
break;
case 1:
fileName_.at(0) = "2006/h12006jetspdf_singlet_fitA.data";
fileName_.at(1) = "2006/h12006jetspdf_gluon_fitA.data";
fileName_.at(2) = "2006/h12006jetspdf_singlet_fitA.data";
nxPoints_ = 100;
nqPoints_ = 30;
break;
case 2:
fileName_.at(0) = "2006/h12006jetspdf_singlet_fitB.data";
fileName_.at(1) = "2006/h12006jetspdf_gluon_fitB.data";
fileName_.at(2) = "2006/h12006jetspdf_singlet_fitB.data";
nxPoints_ = 100;
nqPoints_ = 30;
break;
default:
assert(false);
}
loadTables();
}
ClassDescription<PomeronPDF> PomeronPDF::initPomeronPDF;
// Definition of the static class description member.
void PomeronPDF::Init(){
static ClassDocumentation<PomeronPDF> documentation
("Implementation of the Diffractive PDFs");
static Switch<PomeronPDF,int> interfacePDFFit
("PDFFit",
"Switch between different PDF fits.",
&PomeronPDF::PDFFit_, 0, true, false);
static SwitchOption interfacePDFFit2007
(interfacePDFFit,
"2007",
"Hera fit 2007",
0);
static SwitchOption interfacePDFFit2006A
(interfacePDFFit,
"2006A",
"Hera fit A 2006",
1);
static SwitchOption interfacePDFFit2006B
(interfacePDFFit,
"2006B",
"Hera fit B 2006",
2);
static Parameter<PomeronPDF,string> interfaceRootName
("RootName",
"Root name for the input files",
&PomeronPDF::rootName_, "",
false, false);
static Switch<PomeronPDF,int> interfaceBoundary
("Boundary",
"Switch between different beheviours at out of PDF range.",
&PomeronPDF::boundary_, 0, true, false);
static SwitchOption interfaceBounderyExtrapolate
(interfaceBoundary,
"Extrapolate",
"Extrapolate the PDF when it is out of range.",
0);
static SwitchOption interfaceBounderyFreeze
(interfaceBoundary,
"Freeze",
"Freeze the value of PDF at the boundary.",
1);
}
IBPtr PomeronPDF::clone() const {
return new_ptr(*this);
}
IBPtr PomeronPDF::fullclone() const {
return new_ptr(*this);
}
void PomeronPDF::persistentOutput(PersistentOStream & os) const {
os << lxGrid_ << lqqGrid_ << pdfTable_
<< PDFFit_ << nxPoints_ << nqPoints_
<< rootName_ << boundary_;
}
void PomeronPDF::persistentInput(PersistentIStream & is, int) {
is >> lxGrid_ >> lqqGrid_ >> pdfTable_
>> PDFFit_ >> nxPoints_>> nqPoints_
>> rootName_ >> boundary_;
}
diff --git a/Tests/python/merge-DIS b/Tests/python/merge-DIS
--- a/Tests/python/merge-DIS
+++ b/Tests/python/merge-DIS
@@ -1,152 +1,152 @@
#! /usr/bin/env python
import logging
import sys
if sys.version_info[:3] < (2,4,0):
print "rivet scripts require Python version >= 2.4.0... exiting"
sys.exit(1)
import os, yoda
# #############################################
def fillAbove(desthisto, sourcehistosbyq2):
q2high=1e100
q2low=-1e100
for q2,h in sorted(sourcehistosbyq2.iteritems(),reverse=True) :
q2low=q2
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= q2low and
- h.bins[i].lowEdge < q2high ) :
+ if(h.bins[i].xMin >= q2low and
+ h.bins[i].xMin < q2high ) :
desthisto.bins[i] += h.bins[i]
q2high=q2
def merge(hpath):
global inhistos
global outhistos
try:
fillAbove(outhistos[hpath], inhistos[hpath])
except:
pass
def useOne(hpath, q2):
global inhistos
global outhistos
try:
outhistos[hpath] = inhistos[hpath][float(q2)]
except:
pass
if __name__ == "__main__":
import logging
from optparse import OptionParser, OptionGroup
parser = OptionParser(usage="%prog name")
verbgroup = OptionGroup(parser, "Verbosity control")
verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
default=logging.INFO, help="print debug (very verbose) messages")
verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
default=logging.INFO, help="be very quiet")
parser.add_option_group(verbgroup)
(opts, args) = parser.parse_args()
logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
## Check args
if len(args) < 1:
logging.error("Must specify at least the name of the files")
sys.exit(1)
# #######################################
yodafiles=["-e--LowQ2.yoda","-e+-LowQ2.yoda","-e+-HighQ2.yoda"]
## Get histos
inhistos = {}
outhistos={}
weights = {}
for f in yodafiles:
file = args[0]+f
if(file.find("Low")> 0) :
q2=0
elif(file.find("High")>0) :
q2=50
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
if(aopath.find("4129130")>0 ) :
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(q2):
tmpE[q2] = ao
else:
raise Exception("A set with q2 = %s already exists" % ( q2 ))
elif(aopath.find("2919893")>0 ) :
outhistos[aopath] = ao
else :
logging.error("Analysis %s is unknown" % aopath)
sys.exit(1)
## Make empty output histos if needed
for hpath,hsets in inhistos.iteritems():
if( hpath.find("4129130")>0 and (hpath.find("d33")>0 or hpath.find("d34")>0)) :
if(type(hsets.values()[0])==yoda.core.Scatter2D) :
outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].path,hsets.values()[0].title)
elif(type(hsets.values()[0])==yoda.core.Profile1D) :
outhistos[hpath] = yoda.core.Profile1D(hsets.values()[0].path,hsets.values()[0].title)
for i in range(0,hsets.values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].bins[i].lowEdge,hsets.values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].bins[i].xMin,hsets.values()[0].bins[i].xMax)
else :
logging.error("Histogram %s is of unknown type" % hpath)
sys.exit(1)
# H1 transverse energy
useOne("/H1_2000_S4129130/d01-x01-y01","0")
useOne("/H1_2000_S4129130/d02-x01-y01","0")
useOne("/H1_2000_S4129130/d03-x01-y01","0")
useOne("/H1_2000_S4129130/d04-x01-y01","0")
useOne("/H1_2000_S4129130/d05-x01-y01","0")
useOne("/H1_2000_S4129130/d06-x01-y01","0")
useOne("/H1_2000_S4129130/d07-x01-y01","0")
useOne("/H1_2000_S4129130/d08-x01-y01","0")
useOne("/H1_2000_S4129130/d09-x01-y01","0")
useOne("/H1_2000_S4129130/d10-x01-y01","0")
useOne("/H1_2000_S4129130/d11-x01-y01","0")
useOne("/H1_2000_S4129130/d12-x01-y01","0")
useOne("/H1_2000_S4129130/d13-x01-y01","0")
useOne("/H1_2000_S4129130/d14-x01-y01","0")
useOne("/H1_2000_S4129130/d15-x01-y01","0")
useOne("/H1_2000_S4129130/d16-x01-y01","50")
useOne("/H1_2000_S4129130/d17-x01-y01","50")
useOne("/H1_2000_S4129130/d18-x01-y01","50")
useOne("/H1_2000_S4129130/d19-x01-y01","50")
useOne("/H1_2000_S4129130/d20-x01-y01","50")
useOne("/H1_2000_S4129130/d21-x01-y01","50")
useOne("/H1_2000_S4129130/d22-x01-y01","50")
useOne("/H1_2000_S4129130/d23-x01-y01","50")
useOne("/H1_2000_S4129130/d24-x01-y01","50")
useOne("/H1_2000_S4129130/d25-x01-y01","0")
useOne("/H1_2000_S4129130/d26-x01-y01","0")
useOne("/H1_2000_S4129130/d27-x01-y01","0")
useOne("/H1_2000_S4129130/d28-x01-y01","0")
useOne("/H1_2000_S4129130/d29-x01-y01","50")
useOne("/H1_2000_S4129130/d30-x01-y01","50")
useOne("/H1_2000_S4129130/d31-x01-y01","50")
useOne("/H1_2000_S4129130/d32-x01-y01","50")
merge("/H1_2000_S4129130/d33-x01-y01")
merge("/H1_2000_S4129130/d34-x01-y01")
# Choose output file
name = args[0]+".yoda"
# output the yoda file
yoda.writeYODA(outhistos,name)
sys.exit(0)
diff --git a/Tests/python/merge-LHC-Jets b/Tests/python/merge-LHC-Jets
--- a/Tests/python/merge-LHC-Jets
+++ b/Tests/python/merge-LHC-Jets
@@ -1,1550 +1,1550 @@
#! /usr/bin/env python
import logging
import sys
if sys.version_info[:3] < (2,4,0):
print "rivet scripts require Python version >= 2.4.0... exiting"
sys.exit(1)
import os, yoda
# #############################################
def fillAbove(scale,desthisto, sourcehistosbyptmin):
pthigh= 1e100
ptlow =-1e100
for pt, h in sorted(sourcehistosbyptmin.iteritems(),reverse=True):
ptlow=pt
if(type(desthisto)==yoda.core.Scatter2D) :
for i in range(0,h.numPoints) :
- lowEdge = h.points[i].x-h.points[i].xErrs.minus
- if( lowEdge*scale >= ptlow and
- lowEdge*scale < pthigh ) :
+ xMin = h.points[i].x-h.points[i].xErrs.minus
+ if( xMin*scale >= ptlow and
+ xMin*scale < pthigh ) :
desthisto.addPoint(h.points[i])
elif(type(desthisto)==yoda.core.Profile1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge*scale >= ptlow and
- h.bins[i].lowEdge*scale < pthigh ) :
+ if(h.bins[i].xMin*scale >= ptlow and
+ h.bins[i].xMin*scale < pthigh ) :
desthisto.bins[i] += h.bins[i]
elif(type(desthisto)==yoda.core.Histo1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge*scale >= ptlow and
- h.bins[i].lowEdge*scale < pthigh ) :
+ if(h.bins[i].xMin*scale >= ptlow and
+ h.bins[i].xMin*scale < pthigh ) :
desthisto.bins[i] += h.bins[i]
else :
logging.error("Can't merge %s, unknown type" % desthisto.path)
sys.exit(1)
pthigh=pt
def mergeByPt(hpath, sqrts, scale=1.):
global inhistos_pt
global outhistos
try:
fillAbove(scale,outhistos[hpath], inhistos_pt[hpath][float(sqrts)])
except:
pass
def mergeByMass(hpath, sqrts, scale=1.):
global inhistos_mass
global outhistos
try:
fillAbove(scale,outhistos[hpath], inhistos_mass[hpath][float(sqrts)])
except:
pass
def useOnePt(hpath, sqrts, ptmin):
global inhistos_pt
global outhistos
try:
## Find best pT_min match
ptmins = inhistos_pt[hpath][float(sqrts)].keys()
closest_ptmin = None
for ptm in ptmins:
if closest_ptmin is None or \
abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
closest_ptmin = ptm
if closest_ptmin != float(ptmin):
logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
"using pTmin=%e instead" % closest_ptmin)
outhistos[hpath] = inhistos_pt[hpath][float(sqrts)][closest_ptmin]
except:
pass
def useOneMass(hpath, sqrts, ptmin):
global inhistos_pt
global outhistos
try:
## Find best pT_min match
ptmins = inhistos_mass[hpath][float(sqrts)].keys()
closest_ptmin = None
for ptm in ptmins:
if closest_ptmin is None or \
abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
closest_ptmin = ptm
if closest_ptmin != float(ptmin):
logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
"using pTmin=%e instead" % closest_ptmin)
outhistos[hpath] = inhistos_mass[hpath][float(sqrts)][closest_ptmin]
except:
pass
# #######################################
if __name__ == "__main__":
import logging
from optparse import OptionParser, OptionGroup
parser = OptionParser(usage="%prog name")
verbgroup = OptionGroup(parser, "Verbosity control")
verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
default=logging.INFO, help="print debug (very verbose) messages")
verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
default=logging.INFO, help="be very quiet")
parser.add_option_group(verbgroup)
(opts, args) = parser.parse_args()
logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
## Prefix used in dat file headers
headerprefix = "# "
yodafiles=["-900-UE.yoda" ,"-2360-UE.yoda" ,"-2760-UE.yoda" ,
"-7-UE.yoda" ,"-900-UE-Long.yoda",
"-7-UE-Long.yoda","-7-Jets-0.yoda",
"-7-Jets-All-1.yoda","-7-Jets-All-2.yoda",
"-7-Jets-All-3.yoda","-7-Jets-All-4.yoda","-7-Jets-All-5.yoda",
"-7-Jets-1.yoda","-7-Jets-2.yoda",
"-7-Jets-3.yoda","-7-Jets-4.yoda","-7-Jets-5.yoda",
"-7-Jets-6.yoda","-7-Jets-7.yoda","-7-Jets-8.yoda",
"-7-Heavy-0.yoda","-7-Heavy-1.yoda","-7-Heavy-2.yoda",
"-7-Heavy-3.yoda","-7-Heavy-4.yoda","-7-Heavy-5.yoda",
"-7-Top-SL.yoda","-7-Top-L.yoda"]
## Get histos
inhistos_pt = {}
inhistos_mass = {}
outhistos={}
weights = {}
for f in yodafiles:
file="LHC"+f
sqrts=7000
if(file.find("-900-UE")>0) :
sqrts=900
ptmin=0.
elif(file.find("-2360-UE")>0) :
sqrts=2360
ptmin=0.
elif(file.find("-2760-UE")>0) :
sqrts=2760
ptmin=0.
elif(file.find("-7-UE")>0) :
ptmin=0.
elif(file.find("-7-Jets-0")>0) :
ptmin=4.
elif(file.find("-7-Jets-1")>0) :
ptmin=30.
elif(file.find("-7-Jets-2")>0) :
ptmin=40.
elif(file.find("-7-Jets-3")>0) :
ptmin=80.
elif(file.find("-7-Jets-4")>0) :
ptmin=110.
elif(file.find("-7-Jets-5")>0) :
ptmin=210.
elif(file.find("-7-Jets-6")>0) :
ptmin=260.
elif(file.find("-7-Jets-7")>0) :
ptmin=400.
elif(file.find("-7-Jets-8")>0) :
ptmin=600.
elif(file.find("-7-Heavy-0.yoda")>0) :
ptmin=0.
elif(file.find("-7-Heavy-1.yoda")>0) :
ptmin=10.
elif(file.find("-7-Heavy-2.yoda")>0) :
ptmin=30.
elif(file.find("-7-Heavy-3.yoda")>0) :
ptmin=70.
elif(file.find("-7-Heavy-4.yoda")>0) :
ptmin=100.
elif(file.find("-7-Heavy-5.yoda")>0) :
ptmin=130.
elif(file.find("-7-Top-SL.yoda")>0 or file.find("-7-Top-L.yoda")>0):
ptmin=0.
elif(file.find("-7-Jets-All-1")>0) :
ptmin=25.
elif(file.find("-7-Jets-All-2")>0) :
ptmin=50.
elif(file.find("-7-Jets-All-3")>0) :
ptmin=100.
elif(file.find("-7-Jets-All-4")>0) :
ptmin=170.
elif(file.find("-7-Jets-All-5")>0) :
ptmin=320.
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
if(aopath.find("S8924791")>0 or
aopath.find("S8971293")>0 or
aopath.find("S8817804")>0 or
aopath.find("I1082936")>0 or
aopath.find("S8994773")>0 or
aopath.find("S8918562")>0 or
aopath.find("S8624100")>0 or
aopath.find("S8625980")>0 or
aopath.find("S8894728")>0 or
aopath.find("S8957746")>0 or
aopath.find("S9126244")>0 or
aopath.find("S9120041")>0 or
aopath.find("S8950903")>0 or
aopath.find("S9086218")>0 or
aopath.find("S9088458")>0 or
aopath.find("I919017" )>0 or
aopath.find("I926145" )>0 or
aopath.find("S8941262")>0 or
aopath.find("S8973270")>0 or
aopath.find("I1118269")>0 or
aopath.find("I1188891")>0 or
aopath.find("I1082009")>0 or
aopath.find("I1087342")>0 or
aopath.find("S9035664")>0 or
aopath.find("I1125575")>0 or
aopath.find("I1094564")>0 or
aopath.find("I930220")>0 or
aopath.find("I1224539")>0 or
aopath.find("I1273574")>0 or
aopath.find("I1261026")>0 or
aopath.find("CMS_2012_PAS_QCD_11_010")>0) :
if not inhistos_pt.has_key(aopath):
inhistos_pt[aopath] = {}
tmpE = inhistos_pt[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
if not tmpE[sqrts].has_key(ptmin):
tmpE[sqrts][ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
else :
if(aopath.find("I1243871")>0) :
if(aopath.find("x01")>0 and file.find("-7-Top-L.yoda")>0 ) :
outhistos[aopath] = ao
elif(aopath.find("x02")>0 and file.find("-7-Top-SL.yoda")>0 ) :
outhistos[aopath] = ao
else :
outhistos[aopath] = ao
yodafiles=["-7-Jets-1.yoda" ,"-7-Jets-9.yoda" ,"-7-Jets-10.yoda",
"-7-Jets-11.yoda","-7-Jets-12.yoda","-7-Jets-13.yoda",
"-7-Heavy-6.yoda","-7-Heavy-7.yoda","-7-Heavy-8.yoda"]
for f in yodafiles:
file="LHC"+f
if(file.find("-7-Jets-1.yoda")>0) :
sqrts=7000
mass=0
if(file.find("-7-Jets-9.yoda")>0) :
sqrts=7000
mass=100
elif(file.find("-7-Jets-10.yoda")>0) :
sqrts=7000
mass=300
elif(file.find("-7-Jets-11.yoda")>0) :
sqrts=7000
mass=800
elif(file.find("-7-Jets-12.yoda")>0) :
sqrts=7000
mass=1200
elif(file.find("-7-Jets-13.yoda")>0) :
sqrts=7000
mass=2000
elif(file.find("-7-Heavy-6.yoda")>0) :
sqrts=7000
mass=110
elif(file.find("-7-Heavy-7.yoda")>0) :
sqrts=7000
mass=370
elif(file.find("-7-Heavy-8.yoda")>0) :
sqrts=7000
mass=550
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
if(aopath.find("8817804")>0 or
aopath.find("1082936")>0 or
aopath.find("I930220")>0 or
aopath.find("1261026")>0 or
aopath.find("1090423")>0) :
if not inhistos_mass.has_key(aopath):
inhistos_mass[aopath] = {}
tmpE = inhistos_mass[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(mass):
tmpP[mass] = ao
else:
raise Exception("A set with mass = %s already exists" % ( mass))
## Make empty output histos if needed
for hpath,hsets in inhistos_pt.iteritems():
if( hpath.find("8924791")>0 or
hpath.find("8971293")>0 or
hpath.find("8817804")>0 or
hpath.find("8968497")>0 or
hpath.find("9120041")>0 or
hpath.find("9126244")>0 or
hpath.find("926145") >0 or
hpath.find("9086218")>0 or
hpath.find("1082936")>0 or
hpath.find("8941262")>0 or
hpath.find("1118269")>0 or
hpath.find("1087342")>0 or
hpath.find("1188891")>0 or
hpath.find("919017")>0 or
hpath.find("9035664")>0 or
hpath.find("1125575")>0 or
hpath.find("1094564")>0 or
hpath.find("I930220")>0 or
hpath.find("S9088458")>0 or
hpath.find("I1273574")>0 or
hpath.find("I1261026")>0 or
hpath.find("I1090423")>0 or
hpath.find("CMS_2012_PAS_QCD_11_010")>0) :
if(type(hsets.values()[0].values()[0])==yoda.core.Scatter2D) :
outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
elif(type(hsets.values()[0].values()[0])==yoda.core.Profile1D) :
outhistos[hpath] = yoda.core.Profile1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
elif(type(hsets.values()[0].values()[0])==yoda.core.Histo1D) :
outhistos[hpath] = yoda.core.Histo1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
else :
logging.error("Histogram %s is of unknown type" % hpath)
print hpath,type(hsets.values()[0].values()[0])
sys.exit(1)
logging.info("Processing CMS_2011_S8957746")
useOnePt("/CMS_2011_S8957746/d01-x01-y01", "7000", "80" )
useOnePt("/CMS_2011_S8957746/d02-x01-y01", "7000", "80" )
useOnePt("/CMS_2011_S8957746/d03-x01-y01", "7000", "110" )
useOnePt("/CMS_2011_S8957746/d04-x01-y01", "7000", "110" )
useOnePt("/CMS_2011_S8957746/d05-x01-y01", "7000", "210" )
useOnePt("/CMS_2011_S8957746/d06-x01-y01", "7000", "210" )
logging.info("Processing ATLAS_2010_S8894728")
useOnePt("/ATLAS_2010_S8894728/d01-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d01-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d01-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d02-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d02-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d02-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d03-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d03-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d03-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d04-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d04-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d04-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d05-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d06-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d07-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d08-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d09-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d09-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d09-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d10-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d10-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d10-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d11-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d11-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d11-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d12-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d12-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d12-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d13-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d13-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d13-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d13-x01-y04", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d14-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d14-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d14-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d14-x01-y04", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d15-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d15-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d15-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d15-x01-y04", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d16-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d16-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d16-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d16-x01-y04", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d17-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d17-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d17-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d18-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d18-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d18-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d19-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d19-x01-y02", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d19-x01-y03", "900", "0" )
useOnePt("/ATLAS_2010_S8894728/d20-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d20-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d20-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d21-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8894728/d22-x01-y01", "7000", "0" )
logging.info("Processing ATLAS_2011_S8994773")
useOnePt("/ATLAS_2011_S8994773/d01-x01-y01", "900", "0" )
useOnePt("/ATLAS_2011_S8994773/d02-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2011_S8994773/d03-x01-y01", "900", "0" )
useOnePt("/ATLAS_2011_S8994773/d04-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2011_S8994773/d13-x01-y01", "900", "0" )
useOnePt("/ATLAS_2011_S8994773/d13-x01-y02", "900", "0" )
useOnePt("/ATLAS_2011_S8994773/d13-x01-y03", "900", "0" )
useOnePt("/ATLAS_2011_S8994773/d14-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2011_S8994773/d14-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2011_S8994773/d14-x01-y03", "7000", "0" )
logging.info("Processing ALICE_2010_S8624100")
useOnePt("/ALICE_2010_S8624100/d11-x01-y01", "900", "0" )
useOnePt("/ALICE_2010_S8624100/d12-x01-y01", "900", "0" )
useOnePt("/ALICE_2010_S8624100/d13-x01-y01", "900", "0" )
useOnePt("/ALICE_2010_S8624100/d17-x01-y01","2360", "0" )
useOnePt("/ALICE_2010_S8624100/d18-x01-y01","2360", "0" )
useOnePt("/ALICE_2010_S8624100/d19-x01-y01","2360", "0" )
logging.info("Processing ALICE_2010_S8625980")
useOnePt("/ALICE_2010_S8625980/d03-x01-y01", "7000", "0" )
useOnePt("/ALICE_2010_S8625980/d04-x01-y01", "900", "0" )
useOnePt("/ALICE_2010_S8625980/d05-x01-y01", "2360", "0" )
useOnePt("/ALICE_2010_S8625980/d06-x01-y01", "7000", "0" )
logging.info("Processing ATLAS_2010_S8918562")
useOnePt("/ATLAS_2010_S8918562/d01-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d02-x01-y01", "2360", "0" )
useOnePt("/ATLAS_2010_S8918562/d03-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d04-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d05-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d06-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d07-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d08-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d09-x01-y01", "2360", "0" )
useOnePt("/ATLAS_2010_S8918562/d10-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d11-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d12-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d13-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d14-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d15-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d16-x01-y01", "2360", "0" )
useOnePt("/ATLAS_2010_S8918562/d17-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d18-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d19-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d20-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d21-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d22-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d23-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d24-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d25-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d26-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d27-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d28-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d29-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d30-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d31-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d32-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d33-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d34-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d35-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d36-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d37-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2010_S8918562/d38-x01-y01", "900", "0" )
useOnePt("/ATLAS_2010_S8918562/d39-x01-y01", "7000", "0" )
logging.info("Processing ATLAS_2011_S8971293")
useOnePt("/ATLAS_2011_S8971293/d01-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y03", "7000", "210" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y04", "7000", "260" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y05", "7000", "260" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y06", "7000", "400" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y07", "7000", "400" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y08", "7000", "600" )
useOnePt("/ATLAS_2011_S8971293/d01-x01-y09", "7000", "600" )
logging.info("Processing ATLAS_2011_S8924791")
useOnePt("/ATLAS_2011_S8924791/d01-x01-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x01-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x02-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x02-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x03-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x03-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x04-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x04-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x05-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x05-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x06-y01", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d01-x06-y02", "7000", "30" )
useOnePt("/ATLAS_2011_S8924791/d02-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x02-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x02-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x03-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x03-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x04-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x04-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x05-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x05-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x06-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d02-x06-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x02-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x02-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x03-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x03-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x04-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x04-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x05-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x05-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x06-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d03-x06-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S8924791/d04-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x01-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x02-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x02-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x03-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x03-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x04-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x04-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x05-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x05-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x06-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d04-x06-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S8924791/d05-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x02-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x02-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x03-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x03-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x04-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x04-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x05-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x05-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x06-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d05-x06-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x02-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x02-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x03-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x03-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x04-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x04-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x05-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x05-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x06-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d06-x06-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S8924791/d07-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x02-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x02-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x03-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x03-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x04-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x04-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x05-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x05-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x06-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d07-x06-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S8924791/d08-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x01-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x02-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x02-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x03-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x03-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x04-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x04-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x05-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x05-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x06-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d08-x06-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x01-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x02-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x02-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x03-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x03-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x04-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x04-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x05-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x05-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x06-y01", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d09-x06-y02", "7000", "260" )
useOnePt("/ATLAS_2011_S8924791/d10-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x01-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x02-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x02-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x03-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x03-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x04-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x04-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x05-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x05-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x06-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d10-x06-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x01-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x02-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x02-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x03-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x03-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x04-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x04-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x05-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x05-y02", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x06-y01", "7000", "400" )
useOnePt("/ATLAS_2011_S8924791/d11-x06-y02", "7000", "400" )
logging.info("Processing ATLAS_2010_S8817804")
mergeByPt("/ATLAS_2010_S8817804/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d03-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d04-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d05-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d06-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d07-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d08-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d09-x01-y01", "7000")
mergeByPt("/ATLAS_2010_S8817804/d10-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d11-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d12-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d13-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d14-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d15-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d16-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d17-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d18-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d19-x01-y01", "7000")
mergeByMass("/ATLAS_2010_S8817804/d20-x01-y01", "7000")
useOneMass("/ATLAS_2010_S8817804/d21-x01-y01", "7000", "300" )
useOneMass("/ATLAS_2010_S8817804/d22-x01-y01", "7000", "300" )
useOneMass("/ATLAS_2010_S8817804/d23-x01-y01", "7000", "800" )
useOneMass("/ATLAS_2010_S8817804/d24-x01-y01", "7000", "300" )
useOneMass("/ATLAS_2010_S8817804/d25-x01-y01", "7000", "300" )
useOneMass("/ATLAS_2010_S8817804/d26-x01-y01", "7000", "800" )
logging.info("Processing ATLAS_2011_I930220")
mergeByPt("/ATLAS_2011_I930220/d01-x01-y01", "7000" )
mergeByPt("/ATLAS_2011_I930220/d02-x01-y01", "7000" )
mergeByPt("/ATLAS_2011_I930220/d03-x01-y01", "7000" )
mergeByPt("/ATLAS_2011_I930220/d04-x01-y01", "7000" )
mergeByPt("/ATLAS_2011_I930220/d05-x01-y01", "7000" )
mergeByPt("/ATLAS_2011_I930220/d06-x01-y01", "7000" )
mergeByMass("/ATLAS_2011_I930220/d07-x01-y01", "7000")
useOneMass("/ATLAS_2011_I930220/d08-x01-y01", "7000", "110" )
useOneMass("/ATLAS_2011_I930220/d09-x01-y01", "7000", "110" )
useOneMass("/ATLAS_2011_I930220/d10-x01-y01", "7000", "370" )
logging.info("Processing ATLAS_2012_I1082936")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y02", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y03", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y04", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y05", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y06", "7000")
mergeByPt("/ATLAS_2012_I1082936/d01-x01-y07", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y02", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y03", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y04", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y05", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y06", "7000")
mergeByPt("/ATLAS_2012_I1082936/d02-x01-y07", "7000")
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y01", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y02", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y03", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y04", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y05", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y06", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y07", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y08", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d03-x01-y09", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y01", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y02", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y03", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y04", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y05", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y06", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y07", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y08", "7000", 1000.)
mergeByMass("/ATLAS_2012_I1082936/d04-x01-y09", "7000", 1000.)
logging.info("Processing CMS_2011_S8968497")
useOneMass("/CMS_2011_S8968497/d01-x01-y01", "7000", "2000" )
useOneMass("/CMS_2011_S8968497/d02-x01-y01", "7000", "1200" )
useOneMass("/CMS_2011_S8968497/d03-x01-y01", "7000", "1200" )
useOneMass("/CMS_2011_S8968497/d04-x01-y01", "7000", "800" )
useOneMass("/CMS_2011_S8968497/d05-x01-y01", "7000", "800" )
useOneMass("/CMS_2011_S8968497/d06-x01-y01", "7000", "300" )
useOneMass("/CMS_2011_S8968497/d07-x01-y01", "7000", "300" )
useOneMass("/CMS_2011_S8968497/d08-x01-y01", "7000", "300" )
useOneMass("/CMS_2011_S8968497/d09-x01-y01", "7000", "100" )
logging.info("Processing ATLAS_2011_S9126244")
mergeByPt("/ATLAS_2011_S9126244/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d01-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d02-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d03-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d03-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d04-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d04-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d05-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d05-x01-y02", "7000")
useOnePt("/ATLAS_2011_S9126244/d06-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d06-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d07-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d07-x01-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d08-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d08-x01-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d09-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d09-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d10-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d10-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d11-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d11-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d12-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d12-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d13-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d13-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d14-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d14-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d15-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d15-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d16-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d16-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d17-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d17-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d18-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d18-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d19-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d20-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d21-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d22-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d23-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d24-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d25-x01-y01", "7000", "210" )
mergeByPt("/ATLAS_2011_S9126244/d26-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d26-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d27-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d27-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d28-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d28-x01-y02", "7000")
mergeByPt("/ATLAS_2011_S9126244/d29-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9126244/d29-x01-y02", "7000")
useOnePt("/ATLAS_2011_S9126244/d30-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d31-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d32-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d33-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d34-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d35-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d36-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d37-x01-y01", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d37-x01-y02", "7000", "40" )
useOnePt("/ATLAS_2011_S9126244/d38-x01-y01", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d38-x01-y02", "7000", "80" )
useOnePt("/ATLAS_2011_S9126244/d39-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d39-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d40-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d40-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d41-x01-y01", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d41-x01-y02", "7000", "110" )
useOnePt("/ATLAS_2011_S9126244/d42-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d42-x01-y02", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d43-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2011_S9126244/d43-x01-y02", "7000", "210" )
# CMS_2011_S9120041 UE analysis
logging.info("Processing CMS_2011_S9120041")
mergeByPt("/CMS_2011_S9120041/d01-x01-y01", "7000")
mergeByPt("/CMS_2011_S9120041/d02-x01-y01", "7000")
useOnePt("/CMS_2011_S9120041/d03-x01-y01", "900", "0" )
useOnePt("/CMS_2011_S9120041/d04-x01-y01", "900", "0" )
useOnePt("/CMS_2011_S9120041/d05-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d06-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d07-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d08-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d09-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d10-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d11-x01-y01", "900", "0" )
useOnePt("/CMS_2011_S9120041/d12-x01-y01", "900", "0" )
useOnePt("/CMS_2011_S9120041/d13-x01-y01", "900", "0" )
# CMS dijet decorrelation
logging.info("Processing CMS_2011_S8950903")
useOnePt("/CMS_2011_S8950903/d01-x01-y01", "7000", "80" )
useOnePt("/CMS_2011_S8950903/d02-x01-y01", "7000", "110" )
useOnePt("/CMS_2011_S8950903/d03-x01-y01", "7000", "110" )
useOnePt("/CMS_2011_S8950903/d04-x01-y01", "7000", "210" )
useOnePt("/CMS_2011_S8950903/d05-x01-y01", "7000", "260" )
# CMS jet cross section
logging.info("Processing CMS_2011_S9086218")
mergeByPt("/CMS_2011_S9086218/d01-x01-y01", "7000")
mergeByPt("/CMS_2011_S9086218/d02-x01-y01", "7000")
mergeByPt("/CMS_2011_S9086218/d03-x01-y01", "7000")
mergeByPt("/CMS_2011_S9086218/d04-x01-y01", "7000")
mergeByPt("/CMS_2011_S9086218/d05-x01-y01", "7000")
mergeByPt("/CMS_2011_S9086218/d06-x01-y01", "7000")
# CMS 2/3 jet cross section ratio
logging.info("Processing CMS_2011_S9086218")
mergeByPt("/CMS_2011_S9088458/d01-x01-y01", "7000",500.)
# ATLAS track jet
logging.info("Processing ATLAS_2011_I919017")
mergeByPt("/ATLAS_2011_I919017/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2011_I919017/d01-x01-y02", "7000")
mergeByPt("/ATLAS_2011_I919017/d01-x01-y03", "7000")
mergeByPt("/ATLAS_2011_I919017/d01-x01-y04", "7000")
useOnePt("/ATLAS_2011_I919017/d01-x02-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x02-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x03-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x04-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d01-x05-y25", "7000", "30" )
mergeByPt("/ATLAS_2011_I919017/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2011_I919017/d02-x01-y02", "7000")
mergeByPt("/ATLAS_2011_I919017/d02-x01-y03", "7000")
mergeByPt("/ATLAS_2011_I919017/d02-x01-y04", "7000")
useOnePt("/ATLAS_2011_I919017/d02-x02-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x02-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x03-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x04-y25", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y01", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y02", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y03", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y04", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y05", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y06", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y07", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y08", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y09", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y10", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y11", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y12", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y13", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y14", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y15", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y16", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y17", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y18", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y19", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y20", "7000", "30" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y21", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y22", "7000", "0" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y23", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y24", "7000", "4" )
useOnePt("/ATLAS_2011_I919017/d02-x05-y25", "7000", "30" )
logging.info("Processing ATLAS_2011_I926145")
mergeByPt("/ATLAS_2011_I926145/d01-x01-y01", "7000",1.5)
mergeByPt("/ATLAS_2011_I926145/d02-x01-y01", "7000",1.5)
mergeByPt("/ATLAS_2011_I926145/d03-x01-y01", "7000",1.5)
logging.info("Processing CMS_2011_S8941262")
useOnePt("/CMS_2011_S8941262/d01-x01-y01", "7000", "10" )
useOnePt("/CMS_2011_S8941262/d03-x01-y01", "7000", "10" )
mergeByPt("/CMS_2011_S8941262/d02-x01-y01", "7000",1.5)
logging.info("Processing CMS_2011_S8973270")
useOnePt("/CMS_2011_S8973270/d01-x01-y01", "7000", "70" )
useOnePt("/CMS_2011_S8973270/d02-x01-y01", "7000", "100" )
useOnePt("/CMS_2011_S8973270/d03-x01-y01", "7000", "130" )
useOnePt("/CMS_2011_S8973270/d04-x01-y01", "7000", "70" )
useOnePt("/CMS_2011_S8973270/d05-x01-y01", "7000", "100" )
useOnePt("/CMS_2011_S8973270/d06-x01-y01", "7000", "130" )
logging.info("Processing ATLAS_2012_I1082009")
useOnePt("/ATLAS_2012_I1082009/d08-x01-y01", "7000", "25" )
useOnePt("/ATLAS_2012_I1082009/d09-x01-y01", "7000", "25" )
useOnePt("/ATLAS_2012_I1082009/d10-x01-y01", "7000", "25" )
useOnePt("/ATLAS_2012_I1082009/d11-x01-y01", "7000", "50" )
useOnePt("/ATLAS_2012_I1082009/d12-x01-y01", "7000", "50" )
useOnePt("/ATLAS_2012_I1082009/d13-x01-y01", "7000", "25" )
logging.info("Processing ATLAS_2012_I1118269")
mergeByPt("/ATLAS_2012_I1118269/d01-x01-y01", "7000")
useOnePt("/ATLAS_2012_I1118269/d02-x01-y01", "7000", "10" )
logging.info("Processing ATLAS_2012_I1188891")
mergeByPt("/ATLAS_2012_I1188891/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1188891/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1188891/d03-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1188891/d04-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1188891/d05-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1188891/d06-x01-y01", "7000")
logging.info("Processing CMS_2012_I1087342")
mergeByPt("/CMS_2012_I1087342/d01-x01-y01", "7000")
mergeByPt("/CMS_2012_I1087342/d02-x01-y01", "7000")
mergeByPt("/CMS_2012_I1087342/d03-x01-y01", "7000")
logging.info("Processing CMS_2012_PAS_QCD_11_010")
mergeByPt("/CMS_2012_PAS_QCD_11_010/d01-x01-y01", "7000")
mergeByPt("/CMS_2012_PAS_QCD_11_010/d02-x01-y01", "7000")
mergeByPt("/CMS_2012_PAS_QCD_11_010/d03-x01-y01", "7000")
mergeByPt("/CMS_2012_PAS_QCD_11_010/d04-x01-y01", "7000")
logging.info("Processing ATLAS_2011_S9035664")
mergeByPt("/ATLAS_2011_S9035664/d11-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d12-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d13-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d14-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d15-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d16-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d17-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d18-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d20-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d21-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d22-x01-y01", "7000")
mergeByPt("/ATLAS_2011_S9035664/d23-x01-y01", "7000")
logging.info("Processing ATLAS_2012_I1125575")
mergeByPt("/ATLAS_2012_I1125575/d01-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x01-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x02-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x02-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x03-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x03-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x04-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x04-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x05-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d01-x05-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x01-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x02-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x02-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x03-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x03-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x04-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x04-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x05-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d02-x05-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x01-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x01-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x02-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x02-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x03-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x03-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x04-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x04-y02", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x05-y01", "7000")
mergeByPt("/ATLAS_2012_I1125575/d03-x05-y02", "7000")
useOnePt("/ATLAS_2012_I1125575/d04-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x01-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x02-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x03-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x04-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d04-x05-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x01-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x02-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x03-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x04-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d05-x05-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x01-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x02-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x03-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x04-y20", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y01", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y02", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y03", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y04", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y05", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y06", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y07", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y08", "7000", "0" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y09", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y10", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y11", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y12", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y13", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y14", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y15", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y16", "7000", "4" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y17", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y18", "7000", "30" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y19", "7000", "40" )
useOnePt("/ATLAS_2012_I1125575/d06-x05-y20", "7000", "40" )
# ATLAS_2012_I1094564
useOnePt("/ATLAS_2012_I1094564/d01-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d02-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d03-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d04-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d05-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d06-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d07-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d08-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d09-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d10-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d11-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d12-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d13-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d14-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d15-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d16-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d17-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d18-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d19-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d20-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d21-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d22-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d23-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d24-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d25-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d26-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d27-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d28-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d29-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d30-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d31-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d32-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d33-x01-y01", "7000", "210" )
useOnePt("/ATLAS_2012_I1094564/d34-x01-y01", "7000", "260" )
useOnePt("/ATLAS_2012_I1094564/d35-x01-y01", "7000", "400" )
useOnePt("/ATLAS_2012_I1094564/d36-x01-y01", "7000", "400" )
logging.info("Processing CMS_2013_I1224539_DIJET")
useOnePt("/CMS_2013_I1224539_DIJET/d01-x01-y01", "7000", "210" )
useOnePt("/CMS_2013_I1224539_DIJET/d02-x01-y01", "7000", "260" )
useOnePt("/CMS_2013_I1224539_DIJET/d03-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d04-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d05-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d06-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d07-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d08-x01-y01", "7000", "210" )
useOnePt("/CMS_2013_I1224539_DIJET/d09-x01-y01", "7000", "260" )
useOnePt("/CMS_2013_I1224539_DIJET/d10-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d11-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d12-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d13-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d14-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d15-x01-y01", "7000", "210" )
useOnePt("/CMS_2013_I1224539_DIJET/d16-x01-y01", "7000", "260" )
useOnePt("/CMS_2013_I1224539_DIJET/d17-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d18-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d19-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d20-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d21-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d22-x01-y01", "7000", "210" )
useOnePt("/CMS_2013_I1224539_DIJET/d23-x01-y01", "7000", "260" )
useOnePt("/CMS_2013_I1224539_DIJET/d24-x01-y01", "7000", "400" )
useOnePt("/CMS_2013_I1224539_DIJET/d25-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d26-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d27-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1224539_DIJET/d28-x01-y01", "7000", "600" )
useOnePt("/CMS_2013_I1273574/d01-x01-y01", "7000", "80" )
mergeByPt("/CMS_2013_I1273574/d02-x01-y01", "7000",1.)
useOnePt("/CMS_2013_I1273574/d03-x01-y01", "7000", "80" )
useOnePt("/CMS_2013_I1273574/d04-x01-y01", "7000", "80" )
useOnePt("/CMS_2013_I1273574/d05-x01-y01", "7000", "80" )
useOnePt("/CMS_2013_I1273574/d06-x01-y01", "7000", "80" )
mergeByPt("/CMS_2013_I1273574/d07-x01-y01", "7000",1.)
useOnePt("/CMS_2013_I1273574/d08-x01-y01", "7000", "80" )
mergeByPt("/CMS_2013_I1273574/d09-x01-y01", "7000",1.)
useOnePt("/CMS_2013_I1273574/d10-x01-y01", "7000", "80" )
mergeByPt("/CMS_2013_I1273574/d11-x01-y01", "7000",1.)
useOnePt("/CMS_2013_I1261026/d01-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d02-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d03-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d04-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d05-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d06-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d07-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d08-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d09-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d10-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d11-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d12-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d13-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d14-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d15-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d16-x01-y01", "7000", "0" )
useOnePt("/CMS_2013_I1261026/d17-x01-y01", "7000", "0" )
useOneMass("/CMS_2012_I1090423/d01-x01-y01", "7000", "2000" )
useOneMass("/CMS_2012_I1090423/d02-x01-y01", "7000", "2000" )
useOneMass("/CMS_2012_I1090423/d03-x01-y01", "7000", "1200" )
useOneMass("/CMS_2012_I1090423/d04-x01-y01", "7000", "1200" )
useOneMass("/CMS_2012_I1090423/d05-x01-y01", "7000", "1200" )
useOneMass("/CMS_2012_I1090423/d06-x01-y01", "7000", "800" )
useOneMass("/CMS_2012_I1090423/d07-x01-y01", "7000", "800" )
useOneMass("/CMS_2012_I1090423/d08-x01-y01", "7000", "300" )
useOneMass("/CMS_2012_I1090423/d09-x01-y01", "7000", "300" )
# Choose output file
name = "LHC"+"-Jets.yoda"
yoda.writeYODA(outhistos,name)
sys.exit(0)
diff --git a/Tests/python/merge-Star b/Tests/python/merge-Star
--- a/Tests/python/merge-Star
+++ b/Tests/python/merge-Star
@@ -1,167 +1,167 @@
#! /usr/bin/env python
import logging
import sys
if sys.version_info[:3] < (2,4,0):
print "rivet scripts require Python version >= 2.4.0... exiting"
sys.exit(1)
import os, yoda
# #############################################
def fillAbove(desthisto, sourcehistosbyptmin):
pthigh= 1e100
ptlow =-1e100
for pt, h in sorted(sourcehistosbyptmin.iteritems(),reverse=True):
ptlow=pt
if(type(desthisto)==yoda.core.Scatter2D) :
for i in range(0,h.numPoints) :
- lowEdge = h.points[i].x-h.points[i].xErrs.minus
- if( lowEdge >= ptlow and
- lowEdge < pthigh ) :
+ xMin = h.points[i].x-h.points[i].xErrs.minus
+ if( xMin >= ptlow and
+ xMin < pthigh ) :
desthisto.addPoint(h.points[i])
elif(type(desthisto)==yoda.core.Profile1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= ptlow and
- h.bins[i].lowEdge < pthigh ) :
+ if(h.bins[i].xMin >= ptlow and
+ h.bins[i].xMin < pthigh ) :
desthisto.bins[i] += h.bins[i]
elif(type(desthisto)==yoda.core.Histo1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= ptlow and
- h.bins[i].lowEdge < pthigh ) :
+ if(h.bins[i].xMin >= ptlow and
+ h.bins[i].xMin < pthigh ) :
desthisto.bins[i] += h.bins[i]
else :
logging.error("Can't merge %s, unknown type" % desthisto.path)
sys.exit(1)
pthigh=pt
def mergeByPt(hpath, sqrts):
global inhistos
global outhistos
try:
fillAbove(outhistos[hpath], inhistos[hpath][float(sqrts)])
except:
pass
def useOnePt(hpath, sqrts, ptmin):
global inhistos
global outhistos
try:
## Find best pT_min match
ptmins = inhistos[hpath][float(sqrts)].keys()
closest_ptmin = None
for ptm in ptmins:
if closest_ptmin is None or \
abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
closest_ptmin = ptm
if closest_ptmin != float(ptmin):
logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
"using pTmin=%e instead" % closest_ptmin)
outhistos[hpath] = inhistos[hpath][float(sqrts)][closest_ptmin]
except:
pass
# #######################################
if __name__ == "__main__":
import logging
from optparse import OptionParser, OptionGroup
parser = OptionParser(usage="%prog name")
verbgroup = OptionGroup(parser, "Verbosity control")
verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
default=logging.INFO, help="print debug (very verbose) messages")
verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
default=logging.INFO, help="be very quiet")
parser.add_option_group(verbgroup)
(opts, args) = parser.parse_args()
logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
## Check args
if len(args) < 1:
logging.error("Must specify at least the name of the files")
sys.exit(1)
yodafiles=["-UE.yoda","-Jets-1.yoda","-Jets-2.yoda","-Jets-3.yoda","-Jets-4.yoda"]
## Get histos
inhistos = {}
outhistos={}
weights = {}
for f in yodafiles:
file=args[0]+f
if(file.find("-UE")>0) :
sqrts=200
ptmin=0.
elif(file.find("-Jets-1")>0) :
sqrts=200
ptmin=2.
elif(file.find("-Jets-2")>0) :
sqrts=200
ptmin=5.
elif(file.find("-Jets-3")>0) :
sqrts=200
ptmin=26.
elif(file.find("-Jets-4")>0) :
sqrts=200
ptmin=15.
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
if(aopath.find("6870392")>0 or aopath.find("HELEN")>0) :
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(ptmin):
tmpP[ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
else :
outhistos[aopath] = ao
## Make empty output histos if needed
for hpath,hsets in inhistos.iteritems():
if( (hpath.find("6870392")>0 and hpath.find("d02-x01-y01" )>0 ) or
hpath.find("HELEN")>0 ) :
if(type(hsets.values()[0].values()[0])==yoda.core.Scatter2D) :
outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
elif(type(hsets.values()[0].values()[0])==yoda.core.Profile1D) :
outhistos[hpath] = yoda.core.Profile1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
elif(type(hsets.values()[0].values()[0])==yoda.core.Histo1D) :
outhistos[hpath] = yoda.core.Histo1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
else :
logging.error("Histogram %s is of unknown type" % hpath)
print hpath,type(hsets.values()[0].values()[0])
sys.exit(1)
useOnePt("/STAR_2006_S6870392/d01-x01-y01","200","2")
mergeByPt("/STAR_2006_S6870392/d02-x01-y01","200")
mergeByPt("/STAR_2009_UE_HELEN/d01-x01-y01","200")
mergeByPt("/STAR_2009_UE_HELEN/d02-x01-y01","200")
mergeByPt("/STAR_2009_UE_HELEN/d03-x01-y01","200")
# Choose output file
name = args[0]+".yoda"
print "write to ",name
yoda.writeYODA(outhistos,name)
sys.exit(0)
diff --git a/Tests/python/merge-TVT-Energy b/Tests/python/merge-TVT-Energy
--- a/Tests/python/merge-TVT-Energy
+++ b/Tests/python/merge-TVT-Energy
@@ -1,144 +1,144 @@
#! /usr/bin/env python
import logging
import sys
if sys.version_info[:3] < (2,4,0):
print "rivet scripts require Python version >= 2.4.0... exiting"
sys.exit(1)
import os, yoda
# #############################################
def fillAbove(desthisto, sourcehistosbyptmin):
pthigh= 1e100
ptlow =-1e100
for pt, h in sorted(sourcehistosbyptmin.iteritems(),reverse=True):
ptlow=pt
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= ptlow and
- h.bins[i].lowEdge < pthigh ) :
+ if(h.bins[i].xMin >= ptlow and
+ h.bins[i].xMin < pthigh ) :
desthisto.bins[i] += h.bins[i]
pthigh=pt
def mergeByPt(hpath, sqrts):
global inhistos
global outhistos
try:
fillAbove(outhistos[hpath], inhistos[hpath][float(sqrts)])
except:
pass
# #######################################
if __name__ == "__main__":
import logging
from optparse import OptionParser, OptionGroup
parser = OptionParser(usage="%prog name")
verbgroup = OptionGroup(parser, "Verbosity control")
verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
default=logging.INFO, help="print debug (very verbose) messages")
verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
default=logging.INFO, help="be very quiet")
parser.add_option_group(verbgroup)
(opts, args) = parser.parse_args()
logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
## Check args
if len(args) < 1:
logging.error("Must specify at least the name of the files")
sys.exit(1)
yodafiles=["-Run-II-UE.yoda","-Run-II-Jets-0.yoda",
"-300-UE.yoda" ,"-300-Jets-1.yoda" ,
"-900-UE.yoda" ,"-900-Jets-1.yoda" ]
## Get histos
inhistos = {}
outhistos={}
weights = {}
for f in yodafiles:
file=args[0]+f
if(file.find("Run-II-UE")>0) :
sqrts=1960
ptmin=0.
elif(file.find("Run-II-Jets-0")>0) :
sqrts=1960
ptmin=15
elif(file.find("300-UE")>0) :
sqrts=300
ptmin=0.
elif(file.find("300-Jets-1")>0) :
sqrts=300
ptmin=6.
elif(file.find("900-UE")>0) :
sqrts=900
ptmin=0.
elif(file.find("900-Jets-1")>0) :
sqrts=900
ptmin=10.
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
if(aopath.find("NOTE10874")<0) : continue
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(ptmin):
tmpP[ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
## Make empty output histos if needed
for hpath,hsets in inhistos.iteritems():
if(type(hsets.values()[0].values()[0])==yoda.core.Scatter2D) :
outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
elif(type(hsets.values()[0].values()[0])==yoda.core.Profile1D) :
outhistos[hpath] = yoda.core.Profile1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
else :
logging.error("Histogram %s is of unknown type" % hpath)
sys.exit(1)
## Field analysis
logging.info("Processing CDF_2012_NOTE10874")
## Angular distributions in different pT bins
mergeByPt("/CDF_2012_NOTE10874/d01-x01-y01", "300")
mergeByPt("/CDF_2012_NOTE10874/d01-x01-y02", "900")
mergeByPt("/CDF_2012_NOTE10874/d01-x01-y03","1960")
mergeByPt("/CDF_2012_NOTE10874/d02-x01-y01", "300")
mergeByPt("/CDF_2012_NOTE10874/d02-x01-y02", "900")
mergeByPt("/CDF_2012_NOTE10874/d02-x01-y03","1960")
mergeByPt("/CDF_2012_NOTE10874/d03-x01-y01", "300")
mergeByPt("/CDF_2012_NOTE10874/d03-x01-y02", "900")
mergeByPt("/CDF_2012_NOTE10874/d03-x01-y03","1960")
# Choose output file
for i in range(1,4) :
if(i==1) :
name = args[0]+"-300-Energy.yoda"
elif(i==2) :
name = args[0]+"-900-Energy.yoda"
elif(i==3) :
name = args[0]+"-1960-Energy.yoda"
temphistos={}
for hpath, h in sorted(outhistos.iteritems()):
if(hpath.find("y0"+str(i))>0) :
temphistos[hpath] = h
yoda.writeYODA(temphistos,name)
sys.exit(0)
diff --git a/Tests/python/merge-TVT-Jets b/Tests/python/merge-TVT-Jets
--- a/Tests/python/merge-TVT-Jets
+++ b/Tests/python/merge-TVT-Jets
@@ -1,505 +1,505 @@
#! /usr/bin/env python
import logging
import sys
if sys.version_info[:3] < (2,4,0):
print "rivet scripts require Python version >= 2.4.0... exiting"
sys.exit(1)
import os, yoda
# #############################################
def fillAbove(desthisto, sourcehistosbyptmin):
pthigh= 1e100
ptlow =-1e100
for pt, h in sorted(sourcehistosbyptmin.iteritems(),reverse=True):
ptlow=pt
if(type(desthisto)==yoda.core.Scatter2D) :
for i in range(0,h.numPoints) :
- lowEdge = h.points[i].x-h.points[i].xErrs.minus
- if( lowEdge >= ptlow and
- lowEdge < pthigh ) :
+ xMin = h.points[i].x-h.points[i].xErrs.minus
+ if( xMin >= ptlow and
+ xMin < pthigh ) :
desthisto.addPoint(h.points[i])
elif(type(desthisto)==yoda.core.Profile1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= ptlow and
- h.bins[i].lowEdge < pthigh ) :
+ if(h.bins[i].xMin >= ptlow and
+ h.bins[i].xMin < pthigh ) :
desthisto.bins[i] += h.bins[i]
elif(type(desthisto)==yoda.core.Histo1D) :
for i in range(0,h.numBins) :
- if(h.bins[i].lowEdge >= ptlow and
- h.bins[i].lowEdge < pthigh ) :
+ if(h.bins[i].xMin >= ptlow and
+ h.bins[i].xMin < pthigh ) :
desthisto.bins[i] += h.bins[i]
else :
logging.error("Can't merge %s, unknown type" % desthisto.path)
sys.exit(1)
pthigh=pt
def mergeByPt(hpath, sqrts):
global inhistos
global outhistos
try:
fillAbove(outhistos[hpath], inhistos[hpath][float(sqrts)])
except:
pass
def useOnePt(hpath, sqrts, ptmin):
global inhistos
global outhistos
try:
## Find best pT_min match
ptmins = inhistos[hpath][float(sqrts)].keys()
closest_ptmin = None
for ptm in ptmins:
if closest_ptmin is None or \
abs(ptm-float(ptmin)) < abs(closest_ptmin-float(ptmin)):
closest_ptmin = ptm
if closest_ptmin != float(ptmin):
logging.warning("Inexact match for requested pTmin=%s: " % ptmin + \
"using pTmin=%e instead" % closest_ptmin)
outhistos[hpath] = inhistos[hpath][float(sqrts)][closest_ptmin]
except:
pass
# #######################################
if __name__ == "__main__":
import logging
from optparse import OptionParser, OptionGroup
parser = OptionParser(usage="%prog name")
verbgroup = OptionGroup(parser, "Verbosity control")
verbgroup.add_option("-v", "--verbose", action="store_const", const=logging.DEBUG, dest="LOGLEVEL",
default=logging.INFO, help="print debug (very verbose) messages")
verbgroup.add_option("-q", "--quiet", action="store_const", const=logging.WARNING, dest="LOGLEVEL",
default=logging.INFO, help="be very quiet")
parser.add_option_group(verbgroup)
(opts, args) = parser.parse_args()
logging.basicConfig(level=opts.LOGLEVEL, format="%(message)s")
## Check args
if len(args) < 1:
logging.error("Must specify at least the name of the files")
sys.exit(1)
yodafiles=["-Run-II-UE.yoda" ,"-Run-II-Jets-1.yoda","-Run-II-Jets-2.yoda",\
"-Run-II-Jets-3.yoda","-Run-II-Jets-4.yoda","-Run-II-Jets-5.yoda",\
"-Run-II-Jets-6.yoda","-Run-II-Jets-7.yoda",\
"-Run-I-UE.yoda" ,"-Run-I-Jets-1.yoda","-Run-I-Jets-2.yoda",\
"-Run-I-Jets-3.yoda","-Run-I-Jets-4.yoda","-Run-I-Jets-5.yoda",\
"-630-UE.yoda" ,"-630-Jets-1.yoda" ,"-630-Jets-2.yoda" ,\
"-630-Jets-3.yoda", "-300-Energy.yoda", "-900-Energy.yoda", "-1960-Energy.yoda",
# "-RatioPlots.yoda"
]
## Get histos
inhistos = {}
outhistos={}
for f in yodafiles:
file=args[0]+f
if(file.find("Run-II-UE")>0) :
sqrts=1960
ptmin=0.
elif(file.find("Run-II-Jets-1")>0) :
sqrts=1960
ptmin=36.
elif(file.find("Run-II-Jets-2")>0) :
sqrts=1960
ptmin=55.
elif(file.find("Run-II-Jets-3")>0) :
sqrts=1960
ptmin=75.
elif(file.find("Run-II-Jets-4")>0) :
sqrts=1960
ptmin=100.
elif(file.find("Run-II-Jets-5")>0) :
sqrts=1960
ptmin=125.
elif(file.find("Run-II-Jets-6")>0) :
ptmin=175.
sqrts=1960
elif(file.find("Run-II-Jets-7")>0) :
sqrts=1960
ptmin=265.
elif(file.find("630-UE")>0) :
sqrts=630
ptmin=0.
elif(file.find("630-Jets-1")>0) :
sqrts=630
ptmin=30.
elif(file.find("630-Jets-2")>0) :
sqrts=630
ptmin=55.
elif(file.find("630-Jets-3")>0) :
sqrts=630
ptmin=90.
elif(file.find("Run-I-UE")>0) :
sqrts=1800
ptmin=0.
elif(file.find("Run-I-Jets-1")>0) :
sqrts=1800
ptmin=30.
elif(file.find("Run-I-Jets-2")>0) :
sqrts=1800
ptmin=55.
elif(file.find("Run-I-Jets-3")>0) :
sqrts=1800
ptmin=80.
elif(file.find("Run-I-Jets-4")>0) :
sqrts=1800
ptmin=105.
elif(file.find("Run-I-Jets-5")>0) :
sqrts=1800
ptmin=175.
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this YODA file
for aopath, ao in aos.iteritems() :
# di-jet decorrelations
# jet shapes
if(aopath.find("5992206")>0 or aopath.find("6217184")>0 or
aopath.find("LEADINGJETS")>0 or aopath.find("7662670")>0 or
aopath.find("7057202")>0 or aopath.find("6450792")>0 or
aopath.find("7828950")>0 or aopath.find("4751469")>0 or
aopath.find("5839831")>0 or aopath.find("4563131")>0 or
aopath.find("4517016")>0 or aopath.find("3618439")>0 or
aopath.find("8591881")>0) :
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(ptmin):
tmpP[ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
elif(aopath.find("8233977")>0 or aopath.find("NOTE_9936")>0 or
aopath.find("3905616")>0 or aopath.find("3324664")>0 or
aopath.find("4796047")>0 or aopath.find("1865951")>0 or
aopath.find("2089246")>0 or aopath.find("3108457")>0 or
aopath.find("3349578")>0 or aopath.find("3541940")>0 or
aopath.find("3214044")>0 or aopath.find("2952106")>0 or
aopath.find("NOTE10874")>0 or aopath.find("895662")>0 ) :
outhistos[aopath] = ao
else :
print aopath
quit()
yodafiles=["-Run-II-Jets-8.yoda","-Run-II-Jets-9.yoda","-Run-II-Jets-10.yoda","-Run-II-Jets-11.yoda",\
"-Run-I-Jets-6.yoda","-Run-I-Jets-7.yoda","-Run-I-Jets-8.yoda"]
for f in yodafiles:
file=args[0]+f
if(file.find("Run-II-Jets-8")>0) :
sqrts=1960
ptmin=0.150
elif(file.find("Run-II-Jets-9")>0) :
sqrts=1960
ptmin=0.400
elif(file.find("Run-II-Jets-10")>0) :
sqrts=1960
ptmin=0.600
elif(file.find("Run-II-Jets-11")>0) :
sqrts=1960
ptmin=1.000
elif(file.find("Run-I-Jets-6")>0) :
sqrts=1800
ptmin=0.150
elif(file.find("Run-I-Jets-7")>0) :
sqrts=1800
ptmin=0.5
elif(file.find("Run-I-Jets-8")>0) :
sqrts=1800
ptmin=0.8
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
aos = yoda.read(file)
except:
logging.error("%s can not be parsed as YODA" % file)
break
## Get histos from this AIDA file
for aopath, ao in aos.iteritems() :
if(aopath.find("8566488")>0 or aopath.find("8320160")>0) :
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(ptmin):
tmpP[ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
elif(aopath.find("8093652")>0 or aopath.find("3418421")>0 or
aopath.find("4266730")>0) :
if not inhistos.has_key(aopath):
inhistos[aopath] = {}
tmpE = inhistos[aopath]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(1000.*ptmin):
tmpP[1000.*ptmin] = ao
else:
raise Exception("A set with ptmin = %s already exists" % ( 1000.*ptmin))
## Make empty output histos if needed
for hpath,hsets in inhistos.iteritems():
if( (hpath.find("6217184")>0 and hpath.find("d13-x01-y01")>0 ) or
hpath.find("LEADINGJETS")>0 or hpath.find("7662670")>0 or
hpath.find("7057202")>0 or hpath.find("6450792")>0 or
hpath.find("7828950")>0 or hpath.find("8566488")>0 or
hpath.find("8320160")>0 or hpath.find("8093652")>0 or
hpath.find("4751469")>0 or hpath.find("5839831")>0 or
hpath.find("4563131")>0 or hpath.find("4517016")>0 or
hpath.find("3618439")>0 or hpath.find("4266730")>0 or
hpath.find("3418421")>0 or hpath.find("8591881")>0) :
if(type(hsets.values()[0].values()[0])==yoda.core.Scatter2D) :
outhistos[hpath] = yoda.core.Scatter2D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
elif(type(hsets.values()[0].values()[0])==yoda.core.Profile1D) :
outhistos[hpath] = yoda.core.Profile1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
elif(type(hsets.values()[0].values()[0])==yoda.core.Histo1D) :
outhistos[hpath] = yoda.core.Histo1D(hsets.values()[0].values()[0].path,
hsets.values()[0].values()[0].title)
for i in range(0,hsets.values()[0].values()[0].numBins) :
- outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].lowEdge,
- hsets.values()[0].values()[0].bins[i].highEdge)
+ outhistos[hpath].addBin(hsets.values()[0].values()[0].bins[i].xMin,
+ hsets.values()[0].values()[0].bins[i].xMax)
else :
logging.error("Histogram %s is of unknown type" % hpath)
print hpath,type(hsets.values()[0].values()[0])
sys.exit(1)
## Field analysis
logging.info("Processing CDF_2001_S4751469")
## Angular distributions in different pT bins
useOnePt("/CDF_2001_S4751469/d01-x01-y01", "1800", "0")
useOnePt("/CDF_2001_S4751469/d01-x01-y02", "1800", "0")
useOnePt("/CDF_2001_S4751469/d01-x01-y03", "1800", "30")
useOnePt("/CDF_2001_S4751469/d02-x01-y01", "1800", "0")
useOnePt("/CDF_2001_S4751469/d02-x01-y02", "1800", "0")
useOnePt("/CDF_2001_S4751469/d02-x01-y03", "1800", "30")
## Number, profile in pT_lead (True?)
useOnePt("/CDF_2001_S4751469/d03-x01-y01", "1800", "0")
useOnePt("/CDF_2001_S4751469/d03-x01-y02", "1800", "0")
useOnePt("/CDF_2001_S4751469/d03-x01-y03", "1800", "0")
#mergeByPt("/CDF_2001_S4751469/d04-x01-y01", "1800")
#mergeByPt("/CDF_2001_S4751469/d04-x01-y02", "1800")
#mergeByPt("/CDF_2001_S4751469/d04-x01-y03", "1800")
useOnePt("/CDF_2001_S4751469/d04-x01-y01", "1800", "30")
useOnePt("/CDF_2001_S4751469/d04-x01-y02", "1800", "30")
useOnePt("/CDF_2001_S4751469/d04-x01-y03", "1800", "30")
## pT sums, profile in pT_lead (True?)
useOnePt("/CDF_2001_S4751469/d05-x01-y01", "1800", "0")
useOnePt("/CDF_2001_S4751469/d05-x01-y02", "1800", "0")
useOnePt("/CDF_2001_S4751469/d05-x01-y03", "1800", "0")
# mergeByPt("/CDF_2001_S4751469/d06-x01-y01", "1800")
# mergeByPt("/CDF_2001_S4751469/d06-x01-y02", "1800")
# mergeByPt("/CDF_2001_S4751469/d06-x01-y03", "1800")
useOnePt("/CDF_2001_S4751469/d06-x01-y01", "1800", "30")
useOnePt("/CDF_2001_S4751469/d06-x01-y02", "1800", "30")
useOnePt("/CDF_2001_S4751469/d06-x01-y03", "1800", "30")
## pT distributions (use a specific pT cut run?)
useOnePt("/CDF_2001_S4751469/d07-x01-y01", "1800", "0")
useOnePt("/CDF_2001_S4751469/d07-x01-y02", "1800", "0")
useOnePt("/CDF_2001_S4751469/d07-x01-y03", "1800", "30")
## Acosta analysis
logging.info("Processing CDF_2004_S5839831")
## Mean pT, profile in ET_lead
mergeByPt("/CDF_2004_S5839831/d01-x01-y01", "1800")
mergeByPt("/CDF_2004_S5839831/d01-x01-y02", "1800")
## pT_max,min, profiles in ET_lead
mergeByPt("/CDF_2004_S5839831/d02-x01-y01", "1800")
mergeByPt("/CDF_2004_S5839831/d02-x01-y02", "1800")
mergeByPt("/CDF_2004_S5839831/d02-x01-y03", "1800")
## pT distributions (want to use a specific pT cut run)
useOnePt("/CDF_2004_S5839831/d03-x01-y01", "1800", "30")
useOnePt("/CDF_2004_S5839831/d03-x01-y02", "1800", "80")
useOnePt("/CDF_2004_S5839831/d03-x01-y03", "1800", "105")
useOnePt("/CDF_2004_S5839831/d03-x01-y04", "1800", "105")
useOnePt("/CDF_2004_S5839831/d03-x01-y05", "1800", "175")
## N_max,min, profiles in ET_lead
mergeByPt("/CDF_2004_S5839831/d04-x01-y01", "1800")
mergeByPt("/CDF_2004_S5839831/d04-x01-y02", "1800")
## Min bias dbs (want to use min bias pT cut)
useOnePt("/CDF_2004_S5839831/d05-x01-y01", "1800", "0")
useOnePt("/CDF_2004_S5839831/d06-x01-y01", "1800", "0")
## Swiss Cheese, profile in ET_lead
mergeByPt("/CDF_2004_S5839831/d07-x01-y01", "1800")
mergeByPt("/CDF_2004_S5839831/d07-x01-y02", "1800")
## pT_max,min, profiles in ET_lead
mergeByPt("/CDF_2004_S5839831/d08-x01-y01", "630")
mergeByPt("/CDF_2004_S5839831/d08-x01-y02", "630")
mergeByPt("/CDF_2004_S5839831/d08-x01-y03", "630")
## Swiss Cheese, profile in ET_lead
mergeByPt("/CDF_2004_S5839831/d09-x01-y01", "630")
mergeByPt("/CDF_2004_S5839831/d09-x01-y02", "630")
## Min bias dbs (want to use min bias pT cut)
useOnePt("/CDF_2004_S5839831/d10-x01-y01", "630", "0")
useOnePt("/CDF_2004_S5839831/d11-x01-y01", "630", "0")
## CDF jet shape analysis
logging.info("Processing CDF_2005_S6217184")
useOnePt("/CDF_2005_S6217184/d01-x01-y01", "1960", "36" )
useOnePt("/CDF_2005_S6217184/d01-x01-y02", "1960", "36" )
useOnePt("/CDF_2005_S6217184/d01-x01-y03", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d02-x01-y01", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d02-x01-y02", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d02-x01-y03", "1960", "75" )
useOnePt("/CDF_2005_S6217184/d03-x01-y01", "1960", "75" )
useOnePt("/CDF_2005_S6217184/d03-x01-y02", "1960", "100")
useOnePt("/CDF_2005_S6217184/d03-x01-y03", "1960", "100")
useOnePt("/CDF_2005_S6217184/d04-x01-y01", "1960", "125")
useOnePt("/CDF_2005_S6217184/d04-x01-y02", "1960", "125")
useOnePt("/CDF_2005_S6217184/d04-x01-y03", "1960", "175")
useOnePt("/CDF_2005_S6217184/d05-x01-y01", "1960", "175")
useOnePt("/CDF_2005_S6217184/d05-x01-y02", "1960", "175")
useOnePt("/CDF_2005_S6217184/d05-x01-y03", "1960", "175")
useOnePt("/CDF_2005_S6217184/d06-x01-y01", "1960", "265")
useOnePt("/CDF_2005_S6217184/d06-x01-y02", "1960", "265")
useOnePt("/CDF_2005_S6217184/d06-x01-y03", "1960", "265")
useOnePt("/CDF_2005_S6217184/d07-x01-y01", "1960", "36" )
useOnePt("/CDF_2005_S6217184/d07-x01-y02", "1960", "36" )
useOnePt("/CDF_2005_S6217184/d07-x01-y03", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d08-x01-y01", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d08-x01-y02", "1960", "55" )
useOnePt("/CDF_2005_S6217184/d08-x01-y03", "1960", "75" )
useOnePt("/CDF_2005_S6217184/d09-x01-y01", "1960", "75" )
useOnePt("/CDF_2005_S6217184/d09-x01-y02", "1960", "100")
useOnePt("/CDF_2005_S6217184/d09-x01-y03", "1960", "100")
useOnePt("/CDF_2005_S6217184/d10-x01-y01", "1960", "125")
useOnePt("/CDF_2005_S6217184/d10-x01-y02", "1960", "125")
useOnePt("/CDF_2005_S6217184/d10-x01-y03", "1960", "175")
useOnePt("/CDF_2005_S6217184/d11-x01-y01", "1960", "175")
useOnePt("/CDF_2005_S6217184/d11-x01-y02", "1960", "175")
useOnePt("/CDF_2005_S6217184/d11-x01-y03", "1960", "175")
useOnePt("/CDF_2005_S6217184/d12-x01-y01", "1960", "265")
useOnePt("/CDF_2005_S6217184/d12-x01-y02", "1960", "265")
useOnePt("/CDF_2005_S6217184/d12-x01-y03", "1960", "265")
mergeByPt("/CDF_2005_S6217184/d13-x01-y01", "1960")
# ## CDF dijet mass spectrum
mergeByPt("/CDF_2008_S8093652/d01-x01-y01", "1960")
# ## Rick Field Run-II Leading Jets analysis
# logging.info("Processing CDF_2008_LEADINGJETS")
# ## charged particle density
# mergeByPt("/CDF_2008_LEADINGJETS/d01-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d02-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d03-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d04-x01-y01", "1960")
# ## pT sum density
# mergeByPt("/CDF_2008_LEADINGJETS/d05-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d06-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d07-x01-y01", "1960")
# mergeByPt("/CDF_2008_LEADINGJETS/d08-x01-y01", "1960")
# ## mean pT
# mergeByPt("/CDF_2008_LEADINGJETS/d09-x01-y01", "1960")
## newer version
logging.info("Processing CDF_2010_S8591881_QCD")
mergeByPt("/CDF_2010_S8591881_QCD/d10-x01-y01", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d10-x01-y02", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d10-x01-y03", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d11-x01-y01", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d11-x01-y02", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d11-x01-y03", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d12-x01-y01", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d12-x01-y02", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d12-x01-y03", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d13-x01-y01", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d13-x01-y02", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d13-x01-y03", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d14-x01-y01", "1960")
mergeByPt("/CDF_2010_S8591881_QCD/d15-x01-y01", "1960")
## D0 dijet correlation analysis
logging.info("Processing D0_2004_S5992206")
useOnePt("/D0_2004_S5992206/d01-x02-y01", "1960", "75")
useOnePt("/D0_2004_S5992206/d02-x02-y01", "1960", "100")
useOnePt("/D0_2004_S5992206/d03-x02-y01", "1960", "125")
useOnePt("/D0_2004_S5992206/d04-x02-y01", "1960", "175")
## D0 incl jet cross-section analysis
logging.info("Processing D0_2008_S7662670")
mergeByPt("/D0_2008_S7662670/d01-x01-y01", "1960")
mergeByPt("/D0_2008_S7662670/d02-x01-y01", "1960")
mergeByPt("/D0_2008_S7662670/d03-x01-y01", "1960")
mergeByPt("/D0_2008_S7662670/d04-x01-y01", "1960")
mergeByPt("/D0_2008_S7662670/d05-x01-y01", "1960")
mergeByPt("/D0_2008_S7662670/d06-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d01-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d02-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d03-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d04-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d05-x01-y01", "1960")
mergeByPt("/D0_2010_S8566488/d06-x01-y01", "1960")
# CDF jet cross section
mergeByPt("/CDF_2001_S4563131/d01-x01-y01", "1800")
mergeByPt("/CDF_2001_S4517016/d01-x01-y01", "1800")
mergeByPt("/CDF_2001_S4517016/d02-x01-y01", "1800")
mergeByPt("/CDF_2001_S4517016/d03-x01-y01", "1800")
mergeByPt("/CDF_2001_S4517016/d04-x01-y01", "1800")
useOnePt("/CDF_1998_S3618439/d01-x01-y01", "1800","105")
useOnePt("/CDF_1998_S3618439/d01-x01-y02", "1800","105")
mergeByPt("/CDF_2008_S7828950/d01-x01-y01", "1960")
mergeByPt("/CDF_2008_S7828950/d02-x01-y01", "1960")
mergeByPt("/CDF_2008_S7828950/d03-x01-y01", "1960")
mergeByPt("/CDF_2008_S7828950/d04-x01-y01", "1960")
mergeByPt("/CDF_2008_S7828950/d05-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d01-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d02-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d03-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d04-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d05-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d06-x01-y01", "1960")
mergeByPt("/CDF_2007_S7057202/d07-x01-y01", "1960")
mergeByPt("/CDF_2006_S6450792/d01-x01-y01", "1960")
mergeByPt("/CDF_2000_S4266730/d01-x01-y01", "1800")
useOnePt("/CDF_1996_S3418421/d01-x01-y01","1800","150")
useOnePt("/CDF_1996_S3418421/d01-x01-y02","1800","150")
useOnePt("/CDF_1996_S3418421/d01-x01-y03","1800","150")
useOnePt("/CDF_1996_S3418421/d01-x01-y04","1800","500")
useOnePt("/CDF_1996_S3418421/d01-x01-y05","1800","500")
mergeByPt("/CDF_1996_S3418421/d02-x01-y01","1800")
useOnePt("/D0_2009_S8320160/d01-x01-y01", "1960", "0.15" )
useOnePt("/D0_2009_S8320160/d02-x01-y01", "1960", "0.15" )
useOnePt("/D0_2009_S8320160/d03-x01-y01", "1960", "0.4" )
useOnePt("/D0_2009_S8320160/d04-x01-y01", "1960", "0.4" )
useOnePt("/D0_2009_S8320160/d05-x01-y01", "1960", "0.6" )
useOnePt("/D0_2009_S8320160/d06-x01-y01", "1960", "0.6" )
useOnePt("/D0_2009_S8320160/d07-x01-y01", "1960", "0.6" )
useOnePt("/D0_2009_S8320160/d08-x01-y01", "1960", "0.6" )
useOnePt("/D0_2009_S8320160/d09-x01-y01", "1960", "1.0" )
useOnePt("/D0_2009_S8320160/d10-x01-y01", "1960", "1.0" )
# Choose output file
name = args[0]+"-Jets.yoda"
yoda.writeYODA(outhistos,name)
sys.exit(0)
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Mon, Jan 20, 11:33 PM (1 d, 9 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4220262
Default Alt Text
(135 KB)
Attached To
rHERWIGHG herwighg
Event Timeline
Log In to Comment