Page Menu
Home
HEPForge
Search
Configure Global Search
Log In
Files
F9514737
No One
Temporary
Actions
View File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Flag For Later
Size
130 KB
Subscribers
None
View Options
diff --git a/MatrixElement/MEfftoVH.cc b/MatrixElement/MEfftoVH.cc
--- a/MatrixElement/MEfftoVH.cc
+++ b/MatrixElement/MEfftoVH.cc
@@ -1,422 +1,422 @@
// -*- C++ -*-
//
// This is the implementation of the non-inlined, non-templated member
// functions of the MEfftoVH class.
//
#include "MEfftoVH.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Interface/Parameter.h"
#include "ThePEG/Interface/Switch.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/PDT/EnumParticles.h"
#include "ThePEG/MatrixElement/Tree2toNDiagram.h"
#include "ThePEG/Utilities/SimplePhaseSpace.h"
#include "ThePEG/Handlers/StandardXComb.h"
#include "ThePEG/Cuts/Cuts.h"
#include "Herwig++/Models/StandardModel/StandardModel.h"
#include "Herwig++/MatrixElement/HardVertex.h"
#include "Herwig++/Utilities/Kinematics.h"
#include "ThePEG/PDF/PolarizedBeamParticleData.h"
using namespace Herwig;
void MEfftoVH::persistentOutput(PersistentOStream & os) const {
os << _shapeopt << _wplus << _wminus << _z0 << _higgs
<< _vertexFFW << _vertexFFZ << _vertexWWH << _maxflavour
<< ounit(_mh,GeV) << ounit(_wh,GeV) << _hmass;
}
void MEfftoVH::persistentInput(PersistentIStream & is, int) {
is >> _shapeopt >> _wplus >> _wminus >> _z0 >> _higgs
>> _vertexFFW >> _vertexFFZ >> _vertexWWH >> _maxflavour
>> iunit(_mh,GeV) >> iunit(_wh,GeV) >> _hmass;
}
AbstractClassDescription<MEfftoVH> MEfftoVH::initMEfftoVH;
// Definition of the static class description member.
void MEfftoVH::Init() {
static ClassDocumentation<MEfftoVH> documentation
("The MEfftoVH class is the base class for the Bjirken process f fbar -> V H");
static Switch<MEfftoVH,unsigned int> interfaceShapeOption
("ShapeScheme",
"Option for the treatment of the Higgs resonance shape",
&MEfftoVH::_shapeopt, 2, false, false);
static SwitchOption interfaceStandardShapeFixed
(interfaceShapeOption,
"FixedBreitWigner",
"Breit-Wigner s-channel resonanse",
1);
static SwitchOption interfaceStandardShapeRunning
(interfaceShapeOption,
"MassGenerator",
"Use the mass generator to give the shape",
2);
static SwitchOption interfaceStandardShapeOnShell
(interfaceShapeOption,
"OnShell",
"Produce an on-shell Higgs boson",
0);
static Parameter<MEfftoVH,unsigned int> interfaceMaxFlavour
( "MaxFlavour",
"The heaviest incoming quark flavour this matrix element is allowed to handle "
"(if applicable).",
&MEfftoVH::_maxflavour, 5, 1, 5, false, false, true);
}
unsigned int MEfftoVH::orderInAlphaS() const {
return 0;
}
unsigned int MEfftoVH::orderInAlphaEW() const {
return 3;
}
Energy2 MEfftoVH::scale() const {
return sHat();
}
int MEfftoVH::nDim() const {
return 4 + int(_shapeopt>0);
}
void MEfftoVH::setKinematics() {
DrellYanBase::setKinematics();
}
Selector<MEBase::DiagramIndex>
MEfftoVH::diagrams(const DiagramVector & diags) const {
Selector<DiagramIndex> sel;
for ( DiagramIndex i = 0; i < diags.size(); ++i )
sel.insert(1.0, i);
return sel;
}
Selector<const ColourLines *>
MEfftoVH::colourGeometries(tcDiagPtr ) const {
static ColourLines c1("");
static ColourLines c2("6 -7");
static ColourLines c3("1 -2");
static ColourLines c4("1 -2, 6 -7");
Selector<const ColourLines *> sel;
if(mePartonData()[0]->coloured()) {
if(mePartonData()[4]->coloured()) sel.insert(1.0, &c4);
else sel.insert(1.0, &c3);
}
else {
if(mePartonData()[4]->coloured()) sel.insert(1.0, &c2);
else sel.insert(1.0, &c1);
}
return sel;
}
void MEfftoVH::doinit() {
DrellYanBase::doinit();
// get the vedrtex pointers from the SM object
tcHwSMPtr hwsm= dynamic_ptr_cast<tcHwSMPtr>(standardModel());
// do the initialisation
if(hwsm) {
_vertexFFW = hwsm->vertexFFW();
_vertexFFZ = hwsm->vertexFFZ();
}
else throw InitException() << "Wrong type of StandardModel object in "
<< "MEfftoVH::doinit() the Herwig++"
<< " version must be used"
<< Exception::runerror;
// get the particle data objects for the intermediates
_wplus = getParticleData(ParticleID::Wplus );
_wminus = getParticleData(ParticleID::Wminus);
_z0 = getParticleData(ParticleID::Z0);
// higgs stuff
_mh = _higgs->mass();
_wh = _higgs->width();
if(_higgs->massGenerator()) {
_hmass=dynamic_ptr_cast<GenericMassGeneratorPtr>(_higgs->massGenerator());
}
if(_shapeopt==2&&!_hmass)
throw InitException()
<< "If using the mass generator for the line shape in MEfftoVH::doinit()"
<< "the mass generator must be an instance of the GenericMassGenerator class"
<< Exception::runerror;
}
double MEfftoVH::me2() const {
vector<SpinorWaveFunction> fin,aout;
vector<SpinorBarWaveFunction> ain,fout;
SpinorWaveFunction q(meMomenta()[0],mePartonData()[0],incoming);
SpinorBarWaveFunction qbar(meMomenta()[1],mePartonData()[1],incoming);
SpinorBarWaveFunction f(meMomenta()[3],mePartonData()[3],outgoing);
SpinorWaveFunction fbar(meMomenta()[4],mePartonData()[4],outgoing);
for(unsigned int ix=0;ix<2;++ix) {
q.reset(ix) ; fin.push_back(q);
qbar.reset(ix); ain.push_back(qbar);
f.reset(ix) ;fout.push_back(f);
fbar.reset(ix);aout.push_back(fbar);
}
return helicityME(fin,ain,fout,aout,false)*sHat()*UnitRemoval::InvE2;
}
double MEfftoVH::helicityME(vector<SpinorWaveFunction> & fin ,
vector<SpinorBarWaveFunction> & ain ,
vector<SpinorBarWaveFunction> & fout,
vector<SpinorWaveFunction> & aout,
bool calc) const {
// scale
Energy2 mb2(scale());
// matrix element to be stored
ProductionMatrixElement menew(PDT::Spin1Half,PDT::Spin1Half,PDT::Spin0,
PDT::Spin1Half,PDT::Spin1Half);
// work out the id of the vector boson
int incharge = mePartonData()[0]->iCharge()+mePartonData()[1]->iCharge();
tcPDPtr vec;
if(incharge==0) vec = _z0;
else if(incharge>0) vec = _wplus;
else vec = _wminus;
// vertex for vector boson
AbstractFFVVertexPtr vertex = vec==_z0 ? _vertexFFZ : _vertexFFW;
// wavefunction for the scalar
ScalarWaveFunction higgs(meMomenta()[2],mePartonData()[2],1.,outgoing);
// calculate the matrix element
VectorWaveFunction inter[2];
unsigned int ihel1,ihel2,ohel1,ohel2;
Complex diag;
double me(0.);
for(ihel1=0;ihel1<2;++ihel1) {
for(ihel2=0;ihel2<2;++ihel2) {
// wavefunction for the intermediate 1st vector
inter[0] = vertex->evaluate(mb2,1,vec,fin[ihel1],ain[ihel2]);
// after the emission of the higgs
inter[1] = _vertexWWH->evaluate(mb2,1,vec,inter[0],higgs);
// boson decay piece
for(ohel1=0;ohel1<2;++ohel1) {
for(ohel2=0;ohel2<2;++ohel2) {
diag = vertex->evaluate(sqr(inter[1].particle()->mass()),
aout[ohel2],fout[ohel1],inter[1]);
me += norm(diag);
menew(ihel1,ihel2,0,ohel1,ohel2) = diag;
}
}
}
}
// spin factor
me *=0.25;
tcPolarizedBeamPDPtr beam[2] =
{dynamic_ptr_cast<tcPolarizedBeamPDPtr>(mePartonData()[0]),
dynamic_ptr_cast<tcPolarizedBeamPDPtr>(mePartonData()[1])};
if( beam[0] || beam[1] ) {
RhoDMatrix rho[2] =
{beam[0] ? beam[0]->rhoMatrix() : RhoDMatrix(mePartonData()[0]->iSpin()),
beam[1] ? beam[1]->rhoMatrix() : RhoDMatrix(mePartonData()[1]->iSpin())};
me = menew.average(rho[0],rho[1]);
}
// incoming colour factor
if(mePartonData()[0]->coloured()) me /= 3.;
// outgoing colour factor
if(mePartonData()[3]->coloured()) me *= 3.;
if(calc) _me.reset(menew);
return me;
}
void MEfftoVH::constructVertex(tSubProPtr sub) {
// extract the particles in the hard process
ParticleVector hard;
hard.push_back(sub->incoming().first);
hard.push_back(sub->incoming().second);
hard.push_back(sub->outgoing()[0]);
hard.push_back(sub->outgoing()[1]);
hard.push_back(sub->outgoing()[2]);
// ensure right order
if(hard[0]->id()<0) swap(hard[0],hard[1]);
if(hard[3]->dataPtr()->iSpin()==PDT::Spin0) swap(hard[2],hard[3]);
if(hard[4]->dataPtr()->iSpin()==PDT::Spin0) swap(hard[2],hard[4]);
if(hard[3]->id()<0) swap(hard[3],hard[4]);
vector<SpinorWaveFunction> fin,aout;
vector<SpinorBarWaveFunction> ain,fout;
SpinorWaveFunction( fin ,hard[0],incoming,false,true);
SpinorBarWaveFunction(ain ,hard[1],incoming,false,true);
ScalarWaveFunction( hard[2],outgoing,true);
SpinorBarWaveFunction(fout,hard[3],outgoing,true ,true);
SpinorWaveFunction( aout,hard[4],outgoing,true ,true);
helicityME(fin,ain,fout,aout,true);
// construct the vertex
HardVertexPtr hardvertex=new_ptr(HardVertex());
// set the matrix element for the vertex
hardvertex->ME(_me);
// set the pointers and to and from the vertex
for(unsigned int ix=0;ix<5;++ix) {
tcSpinPtr spin = hard[ix]->spinInfo();
if(ix<2) {
tcPolarizedBeamPDPtr beam =
dynamic_ptr_cast<tcPolarizedBeamPDPtr>(hard[ix]->dataPtr());
if(beam) spin->rhoMatrix() = beam->rhoMatrix();
}
spin->productionVertex(hardvertex);
}
}
bool MEfftoVH::generateKinematics(const double * r) {
using Constants::pi;
// workout the ID of the vector boson
tcPDPtr vec = mePartonData()[0]->iCharge()+mePartonData()[1]->iCharge()!=0
? _wplus : _z0;
// order determined randomly
Energy e = sqrt(sHat())/2.0;
Energy mh(_mh),mv;
double jac(1.);
if(UseRandom::rndbool()) {
double rhomax,rhomin;
// generate the mass of the Higgs
if(_shapeopt!=0) {
Energy mhmax = min(2.*e-vec->massMin(),mePartonData()[2]->massMax());
Energy mhmin = max(ZERO ,mePartonData()[2]->massMin());
if(mhmax<=mhmin) return false;
rhomin = atan2((sqr(mhmin)-sqr(_mh)), _mh*_wh);
rhomax = atan2((sqr(mhmax)-sqr(_mh)), _mh*_wh);
mh = sqrt(_mh*_wh*tan(rhomin+r[3]*(rhomax-rhomin))+sqr(_mh));
jac *= rhomax-rhomin;
}
// generate the mass of the vector boson
Energy2 mvmax2 = sqr(min(2.*e-mh,vec->massMax()));
Energy2 mvmin2 = sqr(vec->massMin());
if(mvmax2<=mvmin2) return false;
rhomin = atan2((mvmin2-sqr(vec->mass())), vec->mass()*vec->width());
rhomax = atan2((mvmax2-sqr(vec->mass())), vec->mass()*vec->width());
mv = sqrt(vec->mass()*vec->width()*tan(rhomin+r[1]*(rhomax-rhomin))
+sqr(vec->mass()));
jac *= rhomax-rhomin;
}
else {
// generate the mass of the vector boson
Energy2 mvmax2 = sqr(min(2.*e,vec->massMax()));
Energy2 mvmin2 = sqr(vec->massMin());
if(mvmax2<=mvmin2) return false;
double rhomin = atan2((mvmin2-sqr(vec->mass())), vec->mass()*vec->width());
double rhomax = atan2((mvmax2-sqr(vec->mass())), vec->mass()*vec->width());
mv = sqrt(vec->mass()*vec->width()*tan(rhomin+r[1]*(rhomax-rhomin))
+sqr(vec->mass()));
jac *= rhomax-rhomin;
// generate the mass of the Higgs
if(_shapeopt!=0) {
Energy mhmax = min(2.*e-mv,mePartonData()[2]->massMax());
Energy mhmin = max(ZERO ,mePartonData()[2]->massMin());
if(mhmax<=mhmin) return false;
rhomin = atan2((sqr(mhmin)-sqr(_mh)), _mh*_wh);
rhomax = atan2((sqr(mhmax)-sqr(_mh)), _mh*_wh);
mh = sqrt(_mh*_wh*tan(rhomin+r[3]*(rhomax-rhomin))+sqr(_mh));
jac *= rhomax-rhomin;
}
}
if(mh+mv>2.*e) return false;
// assign masses
meMomenta()[2].setMass(mh);
for(unsigned int ix=3;ix<5;++ix)
meMomenta()[ix] = Lorentz5Momentum(mePartonData()[ix]->generateMass());
Energy q = ZERO;
Lorentz5Momentum pvec(mv);
try {
q = SimplePhaseSpace::
getMagnitude(sHat(), meMomenta()[2].mass(), mv);
}
catch ( ImpossibleKinematics ) {
return false;
}
Energy ptmin = max(lastCuts().minKT(mePartonData()[2]),
lastCuts().minKT(vec));
Energy2 m22 = meMomenta()[2].mass2();
Energy2 m32 = pvec .mass2();
Energy2 e0e2 = 2.0*e*sqrt(sqr(q) + m22);
Energy2 e1e2 = 2.0*e*sqrt(sqr(q) + m22);
Energy2 e0e3 = 2.0*e*sqrt(sqr(q) + m32);
Energy2 e1e3 = 2.0*e*sqrt(sqr(q) + m32);
Energy2 pq = 2.0*e*q;
double ctmin = -1.0,ctmax = 1.0;
Energy2 thmin = lastCuts().minTij(mePartonData()[0], mePartonData()[2]);
if ( thmin > ZERO ) ctmax = min(ctmax, (e0e2 - m22 - thmin)/pq);
thmin = lastCuts().minTij(mePartonData()[1], mePartonData()[2]);
if ( thmin > ZERO ) ctmin = max(ctmin, (thmin + m22 - e1e2)/pq);
thmin = lastCuts().minTij(mePartonData()[1], mePartonData()[3]);
if ( thmin > ZERO ) ctmax = min(ctmax, (e1e3 - m32 - thmin)/pq);
thmin = lastCuts().minTij(mePartonData()[0], mePartonData()[3]);
if ( thmin > ZERO ) ctmin = max(ctmin, (thmin + m32 - e0e3)/pq);
if ( ptmin > ZERO ) {
double ctm = 1.0 - sqr(ptmin/q);
if ( ctm <= 0.0 ) return false;
ctmin = max(ctmin, -sqrt(ctm));
ctmax = min(ctmax, sqrt(ctm));
}
if ( ctmin >= ctmax ) return false;
-
+ jacobian(1.);
double cth = getCosTheta(ctmin, ctmax, r[0]);
Energy pt = q*sqrt(1.0-sqr(cth));
double phi = UseRandom::rnd()*2.0*pi;
meMomenta()[2].setX(pt*sin(phi));
meMomenta()[2].setY(pt*cos(phi));
meMomenta()[2].setZ(q*cth);
pvec .setX(-pt*sin(phi));
pvec .setY(-pt*cos(phi));
pvec .setZ(-q*cth);
meMomenta()[2].rescaleEnergy();
pvec .rescaleEnergy();
// decay of the vector boson
bool test=Kinematics::twoBodyDecay(pvec,meMomenta()[3].mass(),
meMomenta()[4].mass(),
-1.+2*r[2],r[3]*2.*pi,
meMomenta()[3],meMomenta()[4]);
if(!test) return false;
// check cuts
vector<LorentzMomentum> out;
tcPDVector tout;
for(unsigned int ix=2;ix<5;++ix) {
out .push_back(meMomenta()[ix]);
tout.push_back(mePartonData()[ix]);
}
if ( !lastCuts().passCuts(tout, out, mePartonData()[0], mePartonData()[1]) )
return false;
// jacobian factors
// main piece
jacobian((pq/sHat())*pi*jacobian());
// mass piece
jacobian(jac*jacobian());
// decay piece
Energy p2 = Kinematics::pstarTwoBodyDecay(mv,meMomenta()[3].mass(),
meMomenta()[4].mass());
jacobian(p2/mv/8./sqr(pi)*jacobian());
// jacobian factor for the gauge boson
jacobian((sqr(sqr(mv)-sqr(vec->mass()))+sqr(vec->mass()*vec->width()))/
(vec->mass()*vec->width())*jacobian()/sHat());
return true;
}
CrossSection MEfftoVH::dSigHatDR() const {
using Constants::pi;
// jacobian factor for the higgs
InvEnergy2 bwfact(ZERO);
Energy moff =meMomenta()[2].mass();
if(_shapeopt==1) {
tcPDPtr h0 = mePartonData()[2]->iSpin()==PDT::Spin0 ?
mePartonData()[2] : mePartonData()[3];
bwfact = h0->generateWidth(moff)*moff/pi/
(sqr(sqr(moff)-sqr(_mh))+sqr(_mh*_wh));
}
else if(_shapeopt==2) {
bwfact = _hmass->BreitWignerWeight(moff);
}
double jac1 = _shapeopt!=0 ?
double(bwfact*(sqr(sqr(moff)-sqr(_mh))+sqr(_mh*_wh))/(_mh*_wh)) : 1.;
// answer
return jac1*me2()*jacobian()/(16.0*sqr(pi)*sHat())*sqr(hbarc);
}
diff --git a/Tests/Makefile.am b/Tests/Makefile.am
--- a/Tests/Makefile.am
+++ b/Tests/Makefile.am
@@ -1,307 +1,310 @@
AUTOMAKE_OPTIONS = -Wno-portability
AM_LDFLAGS += -module -avoid-version -rpath /dummy/path/not/used
EXTRA_DIST = Inputs python Rivet
dist-hook:
rm -rf $(distdir)/Inputs/.svn
rm -rf $(distdir)/python/.svn
rm -rf $(distdir)/Rivet/.svn
EXTRA_LTLIBRARIES = LeptonTest.la GammaTest.la HadronTest.la DISTest.la
if WANT_LIBFASTJET
EXTRA_LTLIBRARIES += HadronJetTest.la LeptonJetTest.la
HadronJetTest_la_SOURCES = \
Hadron/VHTest.h Hadron/VHTest.cc\
Hadron/VTest.h Hadron/VTest.cc\
Hadron/HTest.h Hadron/HTest.cc
HadronJetTest_la_CPPFLAGS = $(AM_CPPFLAGS) $(FASTJETINCLUDE) \
-I$(FASTJETPATH)
HadronJetTest_la_LIBADD = $(FASTJETLIBS)
LeptonJetTest_la_SOURCES = \
Lepton/TopDecay.h Lepton/TopDecay.cc
LeptonJetTest_la_CPPFLAGS = $(AM_CPPFLAGS) $(FASTJETINCLUDE) \
-I$(FASTJETPATH)
LeptonJetTest_la_LIBADD = $(FASTJETLIBS)
endif
LeptonTest_la_SOURCES = \
Lepton/VVTest.h Lepton/VVTest.cc \
Lepton/VBFTest.h Lepton/VBFTest.cc \
Lepton/VHTest.h Lepton/VHTest.cc \
Lepton/FermionTest.h Lepton/FermionTest.cc
GammaTest_la_SOURCES = \
Gamma/GammaMETest.h Gamma/GammaMETest.cc \
Gamma/GammaPMETest.h Gamma/GammaPMETest.cc
DISTest_la_SOURCES = \
DIS/DISTest.h DIS/DISTest.cc
HadronTest_la_SOURCES = \
Hadron/HadronVVTest.h Hadron/HadronVVTest.cc\
Hadron/HadronVBFTest.h Hadron/HadronVBFTest.cc\
Hadron/WHTest.h Hadron/WHTest.cc\
Hadron/ZHTest.h Hadron/ZHTest.cc\
Hadron/VGammaTest.h Hadron/VGammaTest.cc\
Hadron/ZJetTest.h Hadron/ZJetTest.cc\
Hadron/WJetTest.h Hadron/WJetTest.cc\
Hadron/QQHTest.h Hadron/QQHTest.cc
REPO = $(top_builddir)/src/HerwigDefaults.rpo
HERWIG = $(top_builddir)/src/Herwig++
HWREAD = $(HERWIG) read --repo $(REPO) -L $(builddir)/.libs
HWRUN = $(HERWIG) run
tests : tests-LEP tests-DIS tests-LHC tests-Gamma
if WANT_LIBFASTJET
tests-LEP : test-LEP-VV test-LEP-VH test-LEP-VBF test-LEP-BB test-LEP-Quarks test-LEP-Leptons \
test-LEP-default test-LEP-Powheg test-LEP-TopDecay
else
tests-LEP : test-LEP-VV test-LEP-VH test-LEP-VBF test-LEP-BB test-LEP-Quarks test-LEP-Leptons
endif
tests-DIS : test-DIS-Charged test-DIS-Neutral
if WANT_LIBFASTJET
tests-LHC : test-LHC-WW test-LHC-WZ test-LHC-ZZ test-LHC-ZGamma test-LHC-WGamma \
test-LHC-ZH test-LHC-WH test-LHC-ZJet test-LHC-WJet test-LHC-Z test-LHC-W test-LHC-ZZVBF test-LHC-VBF \
test-LHC-WWVBF test-LHC-bbH test-LHC-ttH test-LHC-GammaGamma test-LHC-GammaJet test-LHC-Higgs \
test-LHC-HiggsJet test-LHC-QCDFast test-LHC-QCD test-LHC-Top test-LHC-Bottom \
test-LHC-WHJet test-LHC-ZHJet test-LHC-HJet test-LHC-ZShower test-LHC-WShower\
test-LHC-WHJet-Powheg test-LHC-ZHJet-Powheg test-LHC-HJet-Powheg \
test-LHC-ZShower-Powheg test-LHC-WShower-Powheg
else
tests-LHC : test-LHC-WW test-LHC-WZ test-LHC-ZZ test-LHC-ZGamma test-LHC-WGamma \
test-LHC-ZH test-LHC-WH test-LHC-ZJet test-LHC-WJet test-LHC-Z test-LHC-W test-LHC-ZZVBF test-LHC-VBF \
test-LHC-WWVBF test-LHC-bbH test-LHC-ttH test-LHC-GammaGamma test-LHC-GammaJet test-LHC-Higgs \
test-LHC-HiggsJet test-LHC-QCDFast test-LHC-QCD test-LHC-Top
endif
tests-Gamma : test-Gamma-FF test-Gamma-WW test-Gamma-P
if WANT_LIBFASTJET
test-LEP-% : Inputs/LEP-%.in LeptonTest.la LeptonJetTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
else
test-LEP-% : Inputs/LEP-%.in LeptonTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
endif
Rivet-LEP-% : Rivet/LEP-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-BFactory-% : Rivet/BFactory-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-TVT-% : Rivet/TVT-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-DIS-% : Rivet/DIS-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-LHC-% : Rivet/LHC-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-Star-% : Rivet/Star-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-SppS-% : Rivet/SppS-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-ISR-% : Rivet/ISR-%.in
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
Rivet-LEP: Rivet-LEP-22 Rivet-LEP-35 Rivet-LEP-44 Rivet-LEP-91 \
Rivet-LEP-130 Rivet-LEP-133 Rivet-LEP-136 Rivet-LEP-161 \
Rivet-LEP-172 Rivet-LEP-177 Rivet-LEP-183 Rivet-LEP-189 \
Rivet-LEP-192 Rivet-LEP-196 Rivet-LEP-197 Rivet-LEP-200 \
Rivet-LEP-202 Rivet-LEP-206 Rivet-LEP-14 Rivet-LEP-10\
Rivet-LEP-Powheg-14 Rivet-LEP-Powheg-22 \
Rivet-LEP-Powheg-35 Rivet-LEP-Powheg-44 \
Rivet-LEP-Powheg-91 Rivet-LEP-Powheg-130 \
Rivet-LEP-Powheg-133 Rivet-LEP-Powheg-136 \
Rivet-LEP-Powheg-161 Rivet-LEP-Powheg-172 \
Rivet-LEP-Powheg-177 Rivet-LEP-Powheg-183 \
Rivet-LEP-Powheg-189 Rivet-LEP-Powheg-192 \
Rivet-LEP-Powheg-196 Rivet-LEP-Powheg-197 \
Rivet-LEP-Powheg-200 Rivet-LEP-Powheg-202 \
Rivet-LEP-Powheg-206 Rivet-LEP-Powheg-10
for i in LEP-*.aida; do rivet-rmgaps $$i; done;
rm -rf Rivet-LEP
python/merge-LEP LEP
python/merge-LEP LEP-Powheg
rivet-mkhtml -o Rivet-LEP LEP.aida:Hw++ LEP-Powheg.aida:Hw++-Powheg
Rivet-BFactory: Rivet-BFactory-10.52 Rivet-BFactory-Powheg-10.52 \
Rivet-BFactory-10.52-sym Rivet-BFactory-Powheg-10.52-sym \
Rivet-BFactory-10.54 Rivet-BFactory-Powheg-10.54 \
Rivet-BFactory-10.45 Rivet-BFactory-Powheg-10.45 \
Rivet-BFactory-10.58 Rivet-BFactory-Powheg-10.58 \
Rivet-BFactory-Upsilon Rivet-BFactory-Upsilon2 Rivet-BFactory-Upsilon4 \
Rivet-BFactory-Tau
for i in BFactory-*.aida; do rivet-rmgaps $$i; done;
rm -rf Rivet-BFactory
python/merge-BFactory BFactory
python/merge-BFactory BFactory-Powheg
rivet-mkhtml -o Rivet-BFactory BFactory.aida:Hw++ BFactory-Powheg.aida:Hw++-Powheg
Rivet-DIS: Rivet-DIS-e--LowQ2 \
Rivet-DIS-e+-LowQ2 Rivet-DIS-e+-HighQ2\
Rivet-DIS-Powheg-e--LowQ2 \
Rivet-DIS-Powheg-e+-LowQ2 Rivet-DIS-Powheg-e+-HighQ2\
Rivet-DIS-NoME-e--LowQ2 \
Rivet-DIS-NoME-e+-LowQ2 Rivet-DIS-NoME-e+-HighQ2
rivet-rmgaps DIS-e+-LowQ2.aida
rivet-rmgaps DIS-e--LowQ2.aida
rivet-rmgaps DIS-e+-HighQ2.aida
rivet-rmgaps DIS-Powheg-e+-LowQ2.aida
rivet-rmgaps DIS-Powheg-e--LowQ2.aida
rivet-rmgaps DIS-Powheg-e+-HighQ2.aida
rivet-rmgaps DIS-NoME-e+-LowQ2.aida
rivet-rmgaps DIS-NoME-e--LowQ2.aida
rivet-rmgaps DIS-NoME-e+-HighQ2.aida
python/merge-DIS DIS
python/merge-DIS DIS-Powheg
python/merge-DIS DIS-NoME
rivet-mkhtml -o Rivet-DIS DIS.aida:Hw++ DIS-Powheg.aida:Hw++-Powheg DIS-NoME.aida:Hw++-NoME
Rivet-TVT-WZ: Rivet-TVT-Run-I-Z Rivet-TVT-Powheg-Run-I-Z \
Rivet-TVT-Run-I-W Rivet-TVT-Powheg-Run-I-W \
Rivet-TVT-Run-I-WZ Rivet-TVT-Powheg-Run-I-WZ\
Rivet-TVT-Run-II-Z-e Rivet-TVT-Powheg-Run-II-Z-e \
Rivet-TVT-Run-II-Z-mu Rivet-TVT-Powheg-Run-II-Z-mu \
Rivet-TVT-Run-II-W Rivet-TVT-Powheg-Run-II-W
rivet-rmgaps TVT-Run-II-Z-e.aida;
rivet-rmgaps TVT-Powheg-Run-II-Z-e.aida;
rm -rf Rivet-TVT-WZ
python/merge-TVT-EW TVT-Run-II-W.aida TVT-Run-II-Z-{e,mu}.aida\
TVT-Run-I-{W,Z,WZ}.aida -o TVT-WZ.aida
python/merge-TVT-EW TVT-Powheg-Run-II-W.aida TVT-Powheg-Run-II-Z-{e,mu}.aida\
TVT-Powheg-Run-I-{W,Z,WZ}.aida -o TVT-Powheg-WZ.aida
rivet-mkhtml -o Rivet-TVT-WZ TVT-WZ.aida:Hw++ TVT-Powheg-WZ.aida:Hw++-Powheg
Rivet-TVT-Photon: Rivet-TVT-Run-II-DiPhoton Rivet-TVT-Run-II-PromptPhoton
# Rivet-TVT-Run-I-PromptPhoton
rm -rf Rivet-TVT-Photon
python/merge-aida TVT-Run-II-DiPhoton.aida TVT-Run-II-PromptPhoton.aida\
-o TVT-Photon.aida
rivet-mkhtml -o Rivet-TVT-Photon TVT-Photon.aida:Hw++
Rivet-TVT-Jets: Rivet-TVT-Run-II-Jets-1 Rivet-TVT-Run-II-Jets-2 \
Rivet-TVT-Run-II-Jets-3 Rivet-TVT-Run-II-Jets-4 \
Rivet-TVT-Run-II-Jets-5 Rivet-TVT-Run-II-Jets-6 \
Rivet-TVT-Run-II-Jets-7 Rivet-TVT-Run-II-Jets-8 \
Rivet-TVT-Run-II-Jets-9 Rivet-TVT-Run-II-Jets-10\
Rivet-TVT-Run-II-Jets-11 Rivet-TVT-Run-II-UE \
Rivet-TVT-Run-I-Jets-1 Rivet-TVT-Run-I-Jets-2 \
Rivet-TVT-Run-I-Jets-3 Rivet-TVT-Run-I-Jets-4 \
Rivet-TVT-Run-I-Jets-5 Rivet-TVT-Run-I-Jets-6 \
Rivet-TVT-Run-I-Jets-7 Rivet-TVT-Run-I-Jets-8\
Rivet-TVT-Run-I-UE\
Rivet-TVT-630-UE Rivet-TVT-630-Jets-1 \
Rivet-TVT-630-Jets-2 Rivet-TVT-630-Jets-3
rivet-rmgaps TVT-Run-I-Jets-4.aida
rm -rf Rivet-TVT-Jets
python/merge-TVT-Jets TVT
rivet-mkhtml -o Rivet-TVT-Jets TVT-Jets.aida:Hw++
Rivet-LHC-Jets: Rivet-LHC-7-Jets-0 Rivet-LHC-7-Jets-1 Rivet-LHC-7-Jets-2 \
Rivet-LHC-7-Jets-3 Rivet-LHC-7-Jets-4 \
Rivet-LHC-7-Jets-5 Rivet-LHC-7-Jets-6 \
Rivet-LHC-7-Jets-7 Rivet-LHC-7-Jets-8 \
Rivet-LHC-7-Jets-9 Rivet-LHC-7-Jets-10 \
Rivet-LHC-7-Jets-11 Rivet-LHC-7-Jets-12 \
Rivet-LHC-7-Jets-13 Rivet-LHC-7-UE \
- Rivet-LHC-2360-UE Rivet-LHC-900-UE \
+ Rivet-LHC-2360-UE Rivet-LHC-900-UE Rivet-LHC-2760-UE \
+ Rivet-LHC-7-Heavy-0 \
Rivet-LHC-7-Heavy-1 Rivet-LHC-7-Heavy-2 \
Rivet-LHC-7-Heavy-3 Rivet-LHC-7-Heavy-4 \
Rivet-LHC-7-Heavy-5 Rivet-LHC-7-Top-L \
Rivet-LHC-7-Top-SL \
Rivet-LHC-7-UE-Long Rivet-LHC-900-UE-Long \
- Rivet-LHC-7-Jets-All-1 Rivet-LHC-7-Jets-All-2
+ Rivet-LHC-7-Jets-All-1 Rivet-LHC-7-Jets-All-2 \
+ Rivet-LHC-7-Jets-All-3 Rivet-LHC-7-Jets-All-4 \
+ Rivet-LHC-7-Jets-All-5
rm -rf Rivet-LHC-Jets
python/merge-LHC-Jets
rivet-mkhtml -o Rivet-LHC-Jets LHC-Jets.aida:Hw++
Rivet-Star: Rivet-Star-UE Rivet-Star-Jets-1 \
Rivet-Star-Jets-2 Rivet-Star-Jets-3 \
Rivet-Star-Jets-4
rm -rf Rivet-Star
rivet-rmgaps Star-UE.aida
python/merge-Star Star
rivet-mkhtml -o Rivet-Star Star.aida
Rivet-SppS: Rivet-ISR-44-UE Rivet-ISR-63-UE \
Rivet-SppS-63-UE Rivet-ISR-30-UE \
Rivet-ISR-53-UE Rivet-SppS-200-UE \
Rivet-SppS-500-UE Rivet-SppS-900-UE \
Rivet-SppS-546-UE Rivet-SppS-53-UE
rm -rf Rivet-SppS
python/merge-aida ISR-44-UE.aida ISR-63-UE.aida SppS-63-UE.aida ISR-30-UE.aida \
ISR-53-UE.aida SppS-200-UE.aida SppS-500-UE.aida SppS-900-UE.aida \
SppS-546-UE.aida SppS-53-UE.aida -o SppS.aida
rivet-mkhtml -o Rivet-SppS SppS.aida
Rivet-LHC-EW: Rivet-LHC-W-e Rivet-LHC-Powheg-W-e \
Rivet-LHC-W-mu Rivet-LHC-Powheg-W-mu \
Rivet-LHC-Z-e Rivet-LHC-Powheg-Z-e \
Rivet-LHC-Z-mu Rivet-LHC-Powheg-Z-mu \
Rivet-LHC-WZ Rivet-LHC-Powheg-WZ\
Rivet-LHC-WW Rivet-LHC-Powheg-WW\
Rivet-LHC-ZZ Rivet-LHC-Powheg-ZZ
rm -rf Rivet-LHC-EW;
python/merge-LHC-EW LHC-{W-e,W-mu,Z-e,Z-mu,WW,WZ,ZZ}.aida -o LHC-EW.aida;
python/merge-LHC-EW LHC-Powheg-{W-e,W-mu,Z-e,Z-mu,WW,WZ,ZZ}.aida -o LHC-Powheg-EW.aida;
rivet-mkhtml -o Rivet-LHC-EW LHC-EW.aida:Hw++ LHC-Powheg-EW.aida:Hw++-Powheg;
Rivet-LHC-Photon: Rivet-LHC-7-PromptPhoton Rivet-LHC-7-DiPhoton
rm -rf Rivet-LHC-Photon
python/merge-aida LHC-7-PromptPhoton.aida LHC-7-DiPhoton.aida -o LHC-Photon.aida
rivet-mkhtml -o Rivet-LHC-Photon LHC-Photon.aida:Hw++
Rivet-LHC-Higgs: Rivet-LHC-Powheg-ggH Rivet-LHC-ggH Rivet-LHC-ggHJet \
Rivet-LHC-Powheg-VBF Rivet-LHC-VBF Rivet-LHC-WH Rivet-LHC-ZH \
Rivet-LHC-Powheg-WH Rivet-LHC-Powheg-ZH
rm -rf Rivet-LHC-Higgs
rivet-mkhtml -o Rivet-LHC-Higgs LHC-Powheg-ggH.aida:gg-Powheg LHC-ggH.aida:gg LHC-ggHJet.aida:HJet \
LHC-Powheg-VBF.aida:VBF-Powheg LHC-VBF.aida:VBF LHC-WH.aida:WH LHC-ZH.aida:ZH \
LHC-Powheg-WH.aida:WH-Powheg LHC-Powheg-ZH.aida:ZH-Powheg
tests-Rivet : Rivet-LEP Rivet-BFactory Rivet-DIS Rivet-TVT-WZ Rivet-TVT-Photon Rivet-TVT-Jets Rivet-LHC-Jets Rivet-Star Rivet-SppS Rivet-LHC-EW Rivet-LHC-Photon
test-Gamma-% : Inputs/Gamma-%.in GammaTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
test-DIS-% : Inputs/DIS-%.in DISTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
if WANT_LIBFASTJET
test-LHC-% : Inputs/LHC-%.in HadronTest.la GammaTest.la HadronJetTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
else
test-LHC-% : Inputs/LHC-%.in HadronTest.la GammaTest.la
$(HWREAD) $<
$(HWRUN) $(notdir $(subst .in,.run,$<)) -N $${NUMEVENTS:-10000}
endif
clean-local:
rm -f *.out *.log *.tex *.top *.run *.dump *.mult *.Bmult *.aida
diff --git a/Tests/Rivet/LHC-2360-UE.in b/Tests/Rivet/LHC-2760-UE.in
copy from Tests/Rivet/LHC-2360-UE.in
copy to Tests/Rivet/LHC-2760-UE.in
--- a/Tests/Rivet/LHC-2360-UE.in
+++ b/Tests/Rivet/LHC-2760-UE.in
@@ -1,50 +1,42 @@
##################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
# override the beam energy
cd /Herwig/Generators
-set LHCGenerator:EventHandler:LuminosityFunction:Energy 2360.0
+set LHCGenerator:EventHandler:LuminosityFunction:Energy 2760.0
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
# min bias
insert SimpleQCD:MatrixElements[0] MEMinBias
# same as ue process
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
# min ibas cuts
set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
set /Herwig/Cuts/QCDCuts:X1Min 0.01
set /Herwig/Cuts/QCDCuts:X2Min 0.01
##################################################
# select the analyses
##################################################
-# ATLAS UE
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562
-# ALICE charged particles
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8625980
-# ALICE charged particles
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8624100
-# CMS charged particle multiplicity
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8884919
-# CMS charged particle pT and rapidity
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2010_S8547297
+# underlying event forward rapidities
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_FWD_11_003
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
# get the cross section right
create Herwig::MPIXSecReweighter MPIXSecReweighter
insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
-saverun LHC-2360-UE LHCGenerator
+saverun LHC-2760-UE LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-1.in b/Tests/Rivet/LHC-7-Heavy-1.in
--- a/Tests/Rivet/LHC-7-Heavy-1.in
+++ b/Tests/Rivet/LHC-7-Heavy-1.in
@@ -1,35 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
cp MEHeavyQuark MEBottom
set MEBottom:QuarkType Bottom
cp MEHeavyQuark MECharm
set MECharm:QuarkType Charm
insert SimpleQCD:MatrixElements[0] MEBottom
insert SimpleQCD:MatrixElements[0] MECharm
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 5.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I926145
# CMS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8941262
+# b-hadron cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1118269
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Heavy-1 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-2.in b/Tests/Rivet/LHC-7-Heavy-2.in
--- a/Tests/Rivet/LHC-7-Heavy-2.in
+++ b/Tests/Rivet/LHC-7-Heavy-2.in
@@ -1,35 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
cp MEHeavyQuark MEBottom
set MEBottom:QuarkType Bottom
cp MEHeavyQuark MECharm
set MECharm:QuarkType Charm
insert SimpleQCD:MatrixElements[0] MEBottom
insert SimpleQCD:MatrixElements[0] MECharm
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I926145
# CMS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8941262
+# b-hadron cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1118269
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Heavy-2 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-3.in b/Tests/Rivet/LHC-7-Heavy-3.in
--- a/Tests/Rivet/LHC-7-Heavy-3.in
+++ b/Tests/Rivet/LHC-7-Heavy-3.in
@@ -1,36 +1,38 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
cp MEHeavyQuark MEBottom
set MEBottom:QuarkType Bottom
cp MEHeavyQuark MECharm
set MECharm:QuarkType Charm
insert SimpleQCD:MatrixElements[0] MEBottom
insert SimpleQCD:MatrixElements[0] MECharm
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 50.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I926145
# CMS inclusive lepton
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8941262
# CMS displaced vertex
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8973270
+# b-hadron cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1118269
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Heavy-3 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-4.in b/Tests/Rivet/LHC-7-Heavy-4.in
--- a/Tests/Rivet/LHC-7-Heavy-4.in
+++ b/Tests/Rivet/LHC-7-Heavy-4.in
@@ -1,32 +1,34 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
cp MEHeavyQuark MEBottom
set MEBottom:QuarkType Bottom
cp MEHeavyQuark MECharm
set MECharm:QuarkType Charm
insert SimpleQCD:MatrixElements[0] MEBottom
insert SimpleQCD:MatrixElements[0] MECharm
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 80.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# CMS displaced vertex
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8973270
+# b-hadron cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1118269
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Heavy-4 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Heavy-5.in b/Tests/Rivet/LHC-7-Heavy-5.in
--- a/Tests/Rivet/LHC-7-Heavy-5.in
+++ b/Tests/Rivet/LHC-7-Heavy-5.in
@@ -1,32 +1,34 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
cp MEHeavyQuark MEBottom
set MEBottom:QuarkType Bottom
cp MEHeavyQuark MECharm
set MECharm:QuarkType Charm
insert SimpleQCD:MatrixElements[0] MEBottom
insert SimpleQCD:MatrixElements[0] MECharm
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 110.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# CMS displaced vertex
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8973270
+# b-hadron cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1118269
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Heavy-5 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-0.in b/Tests/Rivet/LHC-7-Jets-0.in
--- a/Tests/Rivet/LHC-7-Jets-0.in
+++ b/Tests/Rivet/LHC-7-Jets-0.in
@@ -1,30 +1,34 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 5.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-0 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-1.in b/Tests/Rivet/LHC-7-Jets-1.in
--- a/Tests/Rivet/LHC-7-Jets-1.in
+++ b/Tests/Rivet/LHC-7-Jets-1.in
@@ -1,41 +1,47 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 10.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
+# strange particles in underlying event
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_QCD_11_010
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-1 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-10.in b/Tests/Rivet/LHC-7-Jets-10.in
--- a/Tests/Rivet/LHC-7-Jets-10.in
+++ b/Tests/Rivet/LHC-7-Jets-10.in
@@ -1,33 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
set /Herwig/Cuts/QCDCuts:MHatMin 200.*GeV
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8968497
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-10 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-11.in b/Tests/Rivet/LHC-7-Jets-11.in
--- a/Tests/Rivet/LHC-7-Jets-11.in
+++ b/Tests/Rivet/LHC-7-Jets-11.in
@@ -1,33 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
set /Herwig/Cuts/QCDCuts:MHatMin 600.*GeV
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8968497
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-11 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-12.in b/Tests/Rivet/LHC-7-Jets-12.in
--- a/Tests/Rivet/LHC-7-Jets-12.in
+++ b/Tests/Rivet/LHC-7-Jets-12.in
@@ -1,33 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
set /Herwig/Cuts/QCDCuts:MHatMin 1000.*GeV
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8968497
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-12 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-13.in b/Tests/Rivet/LHC-7-Jets-13.in
--- a/Tests/Rivet/LHC-7-Jets-13.in
+++ b/Tests/Rivet/LHC-7-Jets-13.in
@@ -1,33 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
set /Herwig/Cuts/QCDCuts:MHatMin 1800.*GeV
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8968497
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-13 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-2.in b/Tests/Rivet/LHC-7-Jets-2.in
--- a/Tests/Rivet/LHC-7-Jets-2.in
+++ b/Tests/Rivet/LHC-7-Jets-2.in
@@ -1,40 +1,49 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijets with large rapidity intervals
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1102908
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
+# strange particles in underlying event
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_QCD_11_010
+
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-2 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-3.in b/Tests/Rivet/LHC-7-Jets-3.in
--- a/Tests/Rivet/LHC-7-Jets-3.in
+++ b/Tests/Rivet/LHC-7-Jets-3.in
@@ -1,45 +1,51 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 40.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS event shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8957746
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS dijet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8950903
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
+# strange particles in underlying event
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_QCD_11_010
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-3 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-4.in b/Tests/Rivet/LHC-7-Jets-4.in
--- a/Tests/Rivet/LHC-7-Jets-4.in
+++ b/Tests/Rivet/LHC-7-Jets-4.in
@@ -1,49 +1,53 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 70.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS di-jet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8971293
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS event shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8957746
# ATLAS multijet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9128077
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS dijet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8950903
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-4 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-6.in b/Tests/Rivet/LHC-7-Jets-6.in
--- a/Tests/Rivet/LHC-7-Jets-6.in
+++ b/Tests/Rivet/LHC-7-Jets-6.in
@@ -1,43 +1,47 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 200.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS di-jet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8971293
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS dijet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8950903
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-6 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-7.in b/Tests/Rivet/LHC-7-Jets-7.in
--- a/Tests/Rivet/LHC-7-Jets-7.in
+++ b/Tests/Rivet/LHC-7-Jets-7.in
@@ -1,41 +1,45 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 300.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS di-jet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8971293
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-7 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-8.in b/Tests/Rivet/LHC-7-Jets-8.in
--- a/Tests/Rivet/LHC-7-Jets-8.in
+++ b/Tests/Rivet/LHC-7-Jets-8.in
@@ -1,39 +1,43 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 500.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS di-jet decorrelation
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8971293
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS dijet with veto
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9126244
# CMS jet cross section
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9086218
# CMS 3/2 jet cross section ratio
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9088458
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-8 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-9.in b/Tests/Rivet/LHC-7-Jets-9.in
--- a/Tests/Rivet/LHC-7-Jets-9.in
+++ b/Tests/Rivet/LHC-7-Jets-9.in
@@ -1,33 +1,37 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
set /Herwig/Cuts/QCDCuts:MHatMin 90.*GeV
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# CMS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8968497
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-9 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-All-1.in b/Tests/Rivet/LHC-7-Jets-All-1.in
--- a/Tests/Rivet/LHC-7-Jets-All-1.in
+++ b/Tests/Rivet/LHC-7-Jets-All-1.in
@@ -1,29 +1,31 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set MEQCD2to2:MaximumFlavour 5
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 20.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082009
+# flavour composition of jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1188891
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-All-1 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-Jets-All-2.in b/Tests/Rivet/LHC-7-Jets-All-2.in
--- a/Tests/Rivet/LHC-7-Jets-All-2.in
+++ b/Tests/Rivet/LHC-7-Jets-All-2.in
@@ -1,29 +1,31 @@
###################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
insert SimpleQCD:MatrixElements[0] MEQCD2to2
set MEQCD2to2:MaximumFlavour 5
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
set /Herwig/Cuts/JetKtCut:MinKT 40.
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082009
+# flavour composition of jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1188891
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-7-Jets-All-2 LHCGenerator
diff --git a/Tests/Rivet/LHC-7-UE.in b/Tests/Rivet/LHC-7-UE.in
--- a/Tests/Rivet/LHC-7-UE.in
+++ b/Tests/Rivet/LHC-7-UE.in
@@ -1,74 +1,91 @@
##################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
# min bias
insert SimpleQCD:MatrixElements[0] MEMinBias
# same as ue process
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
# min ibas cuts
set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
set /Herwig/Cuts/QCDCuts:X1Min 0.01
set /Herwig/Cuts/QCDCuts:X2Min 0.01
##################################################
# select the analyses
##################################################
# ATLAS jet shapes
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8924791
# ATLAS jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8817804
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1082936
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8994773
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8894728
# ATLAS fragmentation function
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_CONF_2010_049
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562
# ALICE charged particles
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8625980
# CMS particle spectra
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8978280
# CMS charged particle multiplicity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8884919
# CMS charged particle pT and rapidity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2010_S8656010
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
# ATLAS track jet
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I919017
# LHCB phi production
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2011_I919315
+# LHCB promt hadron productiom
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2012_I1119400
# ATLAS rap gap
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1084540
# CMS forward energy flow
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9215166
# CMS charged particle rapidity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_QCD_10_024
# LHC K0s/Lambda
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2011_I917009
# TOTEM at large rapidity
-insert /Herwig/Analysis/RivetAnalysis:Analyses 0 TOTEM_2012_DNDETA
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 TOTEM_2012_I1115294
# ATLAS Azimuthal ordering of charged hadrons
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1091481
-
+# ALICE single/double diffractive and inelastic sigma
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2012_I1181770
+# ATLAS inelastic cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I894867
+# CMS inelastic cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1193338
+# total transverse energy
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1183818
+# underlying event forward rapidities
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_FWD_11_003
+# strange particles in underlying event
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_QCD_11_010
+# CMS central and forward jets
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1087342
+# CMS dijet
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1184941
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
# get the cross section right
create Herwig::MPIXSecReweighter MPIXSecReweighter
insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
saverun LHC-7-UE LHCGenerator
diff --git a/Tests/Rivet/LHC-900-UE.in b/Tests/Rivet/LHC-900-UE.in
--- a/Tests/Rivet/LHC-900-UE.in
+++ b/Tests/Rivet/LHC-900-UE.in
@@ -1,75 +1,83 @@
##################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
# override the beam energy
cd /Herwig/Generators
set LHCGenerator:EventHandler:LuminosityFunction:Energy 900.0
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
# min bias
insert SimpleQCD:MatrixElements[0] MEMinBias
# same as ue process
set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0
# No QED radiation
erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
# min ibas cuts
set /Herwig/Cuts/JetKtCut:MinKT 0.0*GeV
set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV
set /Herwig/Cuts/QCDCuts:X1Min 0.01
set /Herwig/Cuts/QCDCuts:X2Min 0.01
##################################################
# select the analyses
##################################################
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S8994773
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8918562
# ATLAS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8894728
# ATLAS charged particles
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2010_S8591806
# ALICE charged particles
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8706239
# ALICE charged particles
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8625980
# ALICE charged particles
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2010_S8624100
# ALICE particle spectra
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2011_S8945144
# ALICE strange particle production
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2011_S8909580
# CMS particle spectra
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8978280
# CMS charged particle multiplicity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S8884919
# CMS charged particle pT and rapidity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2010_S8547297
# CMS UE
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9120041
# LHCB k_s0
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2010_S8758301
+# LHCB promt hadron production
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2012_I1119400
# CMS forward energy flow
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2011_S9215166
# CMS charged particle rapidity
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_QCD_10_024
# LHC K0s/Lambda
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 LHCB_2011_I917009
# ATLAS Azimuthal ordering of charged hadrons
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2012_I1091481
+# ALICE single/double diffractive and inelastic sigma
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ALICE_2012_I1181770
+# ATLAS inelastic cross section
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I894867
+# underlying event forward rapidities
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_PAS_FWD_11_003
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
# get the cross section right
create Herwig::MPIXSecReweighter MPIXSecReweighter
insert LHCGenerator:EventHandler:PostSubProcessHandlers 0 MPIXSecReweighter
saverun LHC-900-UE LHCGenerator
diff --git a/Tests/Rivet/LHC-Powheg-Z-mu.in b/Tests/Rivet/LHC-Powheg-Z-mu.in
--- a/Tests/Rivet/LHC-Powheg-Z-mu.in
+++ b/Tests/Rivet/LHC-Powheg-Z-mu.in
@@ -1,34 +1,36 @@
##################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase-Powheg.in
##################################################
# LEP physics parameters (override defaults)
##################################################
cd /Herwig/Generators
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
# Drell-Yan Z/gamma
insert SimpleQCD:MatrixElements[0] PowhegMEqq2gZ2ff
set PowhegMEqq2gZ2ff:Process Muon
# No QED radiation
#erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# ATLAS pT
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9131140
# ATLAS Z+jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I945498
+# CMS underlying event in Z events
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1107658
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-Powheg-Z-mu LHCGenerator
diff --git a/Tests/Rivet/LHC-Z-mu.in b/Tests/Rivet/LHC-Z-mu.in
--- a/Tests/Rivet/LHC-Z-mu.in
+++ b/Tests/Rivet/LHC-Z-mu.in
@@ -1,35 +1,37 @@
##################################################
# Rivet analyses at the Z pole
##################################################
read LHCBase.in
##################################################
# LEP physics parameters (override defaults)
##################################################
cd /Herwig/Generators
##################################################
# Use the q qbar matrix element
##################################################
cd /Herwig/MatrixElements
# Drell-Yan Z/gamma
insert SimpleQCD:MatrixElements[0] MEqq2gZ2ff
set MEqq2gZ2ff:Process Muon
# No QED radiation
#erase /Herwig/EventHandlers/LHCHandler:PostSubProcessHandlers 0
##################################################
# select the analyses
##################################################
# General analysis
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_S9131140
# ATLAS Z+jets
insert /Herwig/Analysis/RivetAnalysis:Analyses 0 ATLAS_2011_I945498
+# CMS underlying event in Z events
+insert /Herwig/Analysis/RivetAnalysis:Analyses 0 CMS_2012_I1107658
##################################################
# Save run for later usage with 'Herwig++ run'
##################################################
cd /Herwig/Generators
saverun LHC-Z-mu LHCGenerator
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,1032 +1,1074 @@
#! /usr/bin/env python
import lighthisto
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, copy, re
from math import sqrt
## Try to load faster but non-standard cElementTree module
try:
import xml.etree.cElementTree as ET
except ImportError:
try:
import cElementTree as ET
except ImportError:
try:
import xml.etree.ElementTree as ET
except:
sys.stderr.write("Can't load the ElementTree XML parser: please install it!\n")
sys.exit(1)
# #############################################
def fillAbove(desthisto, sourcehistosbyptmin):
for i,b in enumerate(desthisto.getBins()):
## Fill bins with pT-ordered histos (so that 'highest always wins')
for ptmin, h in sorted(sourcehistosbyptmin.iteritems()):
newb = h.getBin(i)
if newb.xlow >= ptmin:
b.val = newb.val
b.errplus = newb.errplus
b.errminus = newb.errminus
b._focus= newb._focus
def fillAboveTeV(desthisto, sourcehistosbyptmin):
for i,b in enumerate(desthisto.getBins()):
## Fill bins with pT-ordered histos (so that 'highest always wins')
for ptmin, h in sorted(sourcehistosbyptmin.iteritems()):
newb = h.getBin(i)
if newb.xlow*500. >= ptmin:
b.val = newb.val
b.errplus = newb.errplus
b.errminus = newb.errminus
b._focus= newb._focus
def mergeByPt(hpath, sqrts):
global inhistos_pt
global outhistos
try:
fillAbove(outhistos[hpath], inhistos_pt[hpath][float(sqrts)])
except:
pass
def mergeByPtTeV(hpath, sqrts):
global inhistos_pt
global outhistos
try:
fillAboveTeV(outhistos[hpath], inhistos_pt[hpath][float(sqrts)])
except:
pass
def mergeByMass(hpath, sqrts):
global inhistos_mass
global outhistos
try:
fillAbove(outhistos[hpath], inhistos_mass[hpath][float(sqrts)])
except:
pass
def mergeByMassTeV(hpath, sqrts):
global inhistos_mass
global outhistos
try:
fillAboveTeV(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 = "# "
-aidafiles=["-900-UE.aida" ,"-2360-UE.aida" ,
+aidafiles=["-900-UE.aida" ,"-2360-UE.aida" ,"-2760-UE.aida" ,
"-7-UE.aida" ,"-900-UE-Long.aida",
"-7-UE-Long.aida","-7-Jets-0.aida",
"-7-Jets-All-1.aida","-7-Jets-All-2.aida",
+ "-7-Jets-All-3.aida","-7-Jets-All-4.aida","-7-Jets-All-5.aida",
"-7-Jets-1.aida","-7-Jets-2.aida",
"-7-Jets-3.aida","-7-Jets-4.aida","-7-Jets-5.aida",
"-7-Jets-6.aida","-7-Jets-7.aida","-7-Jets-8.aida",
- "-7-Heavy-1.aida","-7-Heavy-2.aida","-7-Heavy-3.aida",
- "-7-Heavy-4.aida","-7-Heavy-5.aida"]
-
+ "-7-Heavy-0.aida","-7-Heavy-1.aida","-7-Heavy-2.aida",
+ "-7-Heavy-3.aida","-7-Heavy-4.aida","-7-Heavy-5.aida"]
## Get histos
inhistos_pt = {}
inhistos_mass = {}
outhistos={}
weights = {}
for f in aidafiles:
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.aida")>0) :
+ ptmin=0.
elif(file.find("-7-Heavy-1.aida")>0) :
- ptmin=0.
+ ptmin=5.
elif(file.find("-7-Heavy-2.aida")>0) :
ptmin=20.
elif(file.find("-7-Heavy-3.aida")>0) :
ptmin=50.
elif(file.find("-7-Heavy-4.aida")>0) :
ptmin=80.
elif(file.find("-7-Heavy-5.aida")>0) :
ptmin=110.
elif(file.find("-7-Top-SL.aida")>0 or file.find("-7-Top-L.aida")>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:
tree = ET.parse(file)
except:
logging.error("%s can not be parsed as XML" % file)
break
tree = ET.parse(file)
## Get histos from this AIDA file
for dps in tree.findall("dataPointSet"):
h = lighthisto.Histo.fromDPS(dps)
- # di-jet decorrelations
- # jet shapes
if(h.fullPath().find("S8924791")>0 or
h.fullPath().find("S8971293")>0 or
h.fullPath().find("S8817804")>0 or
h.fullPath().find("I1082936")>0 or
h.fullPath().find("S8994773")>0 or
h.fullPath().find("S8918562")>0 or
h.fullPath().find("S8624100")>0 or
h.fullPath().find("S8625980")>0 or
h.fullPath().find("S8894728")>0 or
h.fullPath().find("S8957746")>0 or
h.fullPath().find("S9126244")>0 or
h.fullPath().find("S9120041")>0 or
h.fullPath().find("S8950903")>0 or
h.fullPath().find("S9086218")>0 or
h.fullPath().find("S9088458")>0 or
h.fullPath().find("I919017" )>0 or
h.fullPath().find("I926145" )>0 or
h.fullPath().find("S8941262")>0 or
h.fullPath().find("S8973270")>0 or
- h.fullPath().find("I1082009")>0) :
+ h.fullPath().find("I1118269")>0 or
+ h.fullPath().find("I1188891")>0 or
+ h.fullPath().find("I1082009")>0 or
+ h.fullPath().find("I1087342")>0 or
+ h.fullPath().find("CMS_2012_PAS_QCD_11_010")>0) :
if not inhistos_pt.has_key(h.fullPath()):
inhistos_pt[h.fullPath()] = {}
tmpE = inhistos_pt[h.fullPath()]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
- tmpP = tmpE[sqrts]
- if not tmpP.has_key(ptmin):
- tmpP[ptmin] = h
+ if not tmpE[sqrts].has_key(ptmin):
+ tmpE[sqrts][ptmin] = h
else:
raise Exception("A set with ptmin = %s already exists" % ( ptmin))
else :
outhistos[h.fullPath()] = h
aidafiles=["-7-Jets-9.aida" ,"-7-Jets-10.aida",
"-7-Jets-11.aida","-7-Jets-12.aida","-7-Jets-13.aida"]
for f in aidafiles:
file="LHC"+f
if(file.find("-7-Jets-9.aida")>0) :
sqrts=7000
mass=110
elif(file.find("-7-Jets-10.aida")>0) :
sqrts=7000
mass=300
elif(file.find("-7-Jets-11.aida")>0) :
sqrts=7000
mass=800
elif(file.find("-7-Jets-12.aida")>0) :
sqrts=7000
mass=1200
elif(file.find("-7-Jets-13.aida")>0) :
sqrts=7000
mass=2000
if not os.access(file, os.R_OK):
logging.error("%s can not be read" % file)
break
try:
tree = ET.parse(file)
except:
logging.error("%s can not be parsed as XML" % file)
break
tree = ET.parse(file)
## Get histos from this AIDA file
for dps in tree.findall("dataPointSet"):
h = lighthisto.Histo.fromDPS(dps)
if(h.fullPath().find("8817804")>0 or
h.fullPath().find("1082936")>0 or
h.fullPath().find("8968497")>0) :
if not inhistos_mass.has_key(h.fullPath()):
inhistos_mass[h.fullPath()] = {}
tmpE = inhistos_mass[h.fullPath()]
if not tmpE.has_key(sqrts):
tmpE[sqrts] = {}
tmpP = tmpE[sqrts]
if not tmpP.has_key(mass):
tmpP[mass] = h
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
- h.fullPath().find("8971293")>0 or
- h.fullPath().find("8817804")>0 or
- h.fullPath().find("8968497")>0 or
- h.fullPath().find("9126244")>0 or
- h.fullPath().find("1082936")>0 ) :
+ hpath.find("8971293")>0 or
+ hpath.find("8817804")>0 or
+ hpath.find("8968497")>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("CMS_2012_PAS_QCD_11_010")>0) :
workhisto = copy.deepcopy(hsets.values()[0].values()[0])
outhistos[hpath] = workhisto
## Empty the bin set for histos which we're going to merge
for b in outhistos[hpath]:
b.val = 0
b.errplus = 0
b.errminus = 0
b._focus= None
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_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")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y01", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y02", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y03", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y04", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y05", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y06", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y07", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y08", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d03-x01-y09", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y01", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y02", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y03", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y04", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y05", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y06", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y07", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y08", "7000")
mergeByMassTeV("/ATLAS_2012_I1082936/d04-x01-y09", "7000")
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", "110" )
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")
useOnePt("/CMS_2011_S9120041/d01-x01-y01", "7000", "0" )
useOnePt("/CMS_2011_S9120041/d02-x01-y01", "7000", "0" )
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")
mergeByPtTeV("/CMS_2011_S9088458/d01-x01-y01", "7000")
# 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")
mergeByPt("/ATLAS_2011_I926145/d02-x01-y01", "7000")
mergeByPt("/ATLAS_2011_I926145/d03-x01-y01", "7000")
logging.info("Processing CMS_2011_S8941262")
-useOnePt("/CMS_2011_S8941262/d01-x01-y01", "7000", "0" )
-useOnePt("/CMS_2011_S8941262/d03-x01-y01", "7000", "0" )
+useOnePt("/CMS_2011_S8941262/d01-x01-y01", "7000", "5" )
+useOnePt("/CMS_2011_S8941262/d03-x01-y01", "7000", "5" )
mergeByPt("/CMS_2011_S8941262/d02-x01-y01", "7000")
logging.info("Processing CMS_2011_S8973270")
useOnePt("/CMS_2011_S8973270/d01-x01-y01", "7000", "50" )
useOnePt("/CMS_2011_S8973270/d02-x01-y01", "7000", "80" )
useOnePt("/CMS_2011_S8973270/d03-x01-y01", "7000", "110" )
useOnePt("/CMS_2011_S8973270/d04-x01-y01", "7000", "50" )
useOnePt("/CMS_2011_S8973270/d05-x01-y01", "7000", "80" )
useOnePt("/CMS_2011_S8973270/d06-x01-y01", "7000", "110" )
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", "5" )
+
+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")
+
# Choose output file
name = "LHC"+"-Jets.aida"
out = open(name, "w")
## Write out merged histos
out.write('<?xml version="1.0" encoding="ISO-8859-1" ?>\n')
out.write('<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">\n')
out.write('<aida version="3.3">\n')
out.write(' <implementation version="1.1" package="FreeHEP"/>\n')
for hpath, h in sorted(outhistos.iteritems()):
logging.debug("hpath = %s" % hpath)
out.write(h.asAIDA() + "\n\n")
out.write('</aida>\n')
sys.exit(0)
diff --git a/UnderlyingEvent/MPIXSecReweighter.h b/UnderlyingEvent/MPIXSecReweighter.h
--- a/UnderlyingEvent/MPIXSecReweighter.h
+++ b/UnderlyingEvent/MPIXSecReweighter.h
@@ -1,138 +1,133 @@
// -*- C++ -*-
//
// MPIXSecReweighter.h 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.
//
#ifndef Herwig_MPIXSecReweighter_H
#define Herwig_MPIXSecReweighter_H
//
// This is the declaration of the MPIXSecReweighter class.
//
#include "ThePEG/Handlers/StepHandler.h"
namespace Herwig {
using namespace ThePEG;
/**
* \ingroup UnderlyingEvent
*
* \brief MPIXSecReweighter sets up the proper minimum bias cross
* section.
*
*/
class MPIXSecReweighter: public StepHandler {
public:
/** @name Standard constructors and destructors. */
//@{
/**
* The default constructor.
*/
MPIXSecReweighter();
/**
* The destructor.
*/
virtual ~MPIXSecReweighter();
//@}
public:
/** @name Virtual functions required by the StepHandler class. */
//@{
/**
* The main function called by the EventHandler class to
* perform a step. Given the current state of an Event, this function
* performs the event generation step and includes the result in a new
* Step object int the Event record.
* @param eh the EventHandler in charge of the Event generation.
* @param tagged if not empty these are the only particles which should
* be considered by the StepHandler.
* @param hint a Hint object with possible information from previously
* performed steps.
* @throws Veto if the StepHandler requires the current step to be discarded.
* @throws Stop if the generation of the current Event should be stopped
* after this call.
* @throws Exception if something goes wrong.
*/
virtual void handle(EventHandler & eh, const tPVector & tagged,
const Hint & hint);
//@}
public:
/** @name Functions used by the persistent I/O system. */
//@{
/**
* Function used to write out object persistently.
* @param os the persistent output stream written to.
*/
void persistentOutput(PersistentOStream & os) const;
/**
* Function used to read in object persistently.
* @param is the persistent input stream read from.
* @param version the version number of the object when written.
*/
void persistentInput(PersistentIStream & is, int version);
//@}
/**
* The standard Init function used to initialize the interfaces.
* Called exactly once for each class by the class description system
* before the main function starts or
* when this class is dynamically loaded.
*/
static void Init();
protected:
/** @name Clone Methods. */
//@{
/**
* Make a simple clone of this object.
* @return a pointer to the new object.
*/
virtual IBPtr clone() const;
/** Make a clone of this object, possibly modifying the cloned object
* to make it sane.
* @return a pointer to the new object.
*/
virtual IBPtr fullclone() const;
//@}
-
-// If needed, insert declarations of virtual function defined in the
-// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
-
-
private:
/**
* The assignment operator is private and must never be called.
* In fact, it should not even be implemented.
*/
MPIXSecReweighter & operator=(const MPIXSecReweighter &);
/**
* The sum of weights currently accumulated.
*/
double sumWeights;
/**
* The integrated (ME) cross section currently accumulated.
*/
CrossSection xSec;
};
}
#endif /* Herwig_MPIXSecReweighter_H */
File Metadata
Details
Attached
Mime Type
text/x-diff
Expires
Mon, Feb 24, 6:39 AM (1 d, 14 h)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4485003
Default Alt Text
(130 KB)
Attached To
rHERWIGHG herwighg
Event Timeline
Log In to Comment