[Rivet-svn] r2583 - in trunk: data/plotinfo src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Tue Jul 13 11:54:00 BST 2010


Author: richardn
Date: Tue Jul 13 11:54:36 2010
New Revision: 2583

Log:
fixed non-LEP1 plots and added missing distributions. Validation plots available here http://www.ippp.dur.ac.uk/~richardn/Rivet-LEP/ALEPH_2004_S5765862/index.html

Modified:
   trunk/data/plotinfo/ALEPH_2004_S5765862.plot
   trunk/src/Analyses/ALEPH_2004_S5765862.cc

Modified: trunk/data/plotinfo/ALEPH_2004_S5765862.plot
==============================================================================
--- trunk/data/plotinfo/ALEPH_2004_S5765862.plot	Mon Jul 12 23:28:32 2010	(r2582)
+++ trunk/data/plotinfo/ALEPH_2004_S5765862.plot	Tue Jul 13 11:54:36 2010	(r2583)
@@ -1,3 +1,395 @@
+# BEGIN PLOT /ALEPH_2004_S5765862/d01-x01-y01
+Title=Charged multiplicity at a function of energy
+XLabel=$E_\mathrm{CMS}/GeV$
+YLabel=$N_\text{ch}$
+LegendXPos=0.65
+LegendYPos=0.85
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d02-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d03-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d04-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d05-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d06-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d07-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d08-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d09-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$x_p$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_p}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d11-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d12-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d13-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d14-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d15-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d16-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d17-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d18-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\xi_p$
+YLabel=$1/\sigma \, \text{d}\sigma/\text{d}\xi_p$
+FullRange=1
+LegendYPos=0.65
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d19-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d20-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d21-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d22-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d23-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d24-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d25-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d26-x01-y01
+Title=Charged particle spectrum ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$x_E$
+YLabel=$1/\sigma \, \text{d}{\sigma}/\text{d}{x_E}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d27-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d28-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d29-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d30-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d31-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d32-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d33-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d34-x01-y01
+Title=In-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$p_\perp^\text{in}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{in}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d35-x01-y01
+Title=Out-of-plane $p_\perp$ in GeV w.r.t. thrust axes ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$p_\perp^\text{out}$ / GeV
+YLabel=$N \, \text{d}{\sigma}/\text{d}{p_\perp^\text{out}}$
+FullRange=1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d36-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d37-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d38-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d39-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d40-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d41-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d42-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d43-x01-y01
+Title=Rapidity w.r.t. thrust axes, $y_T$ ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$y_T$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_T}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d44-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d45-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d46-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d47-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d48-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d49-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=196$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d50-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d51-x01-y01
+Title=Rapidity w.r.t. sphericity axes, $y_S$ ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$y_S$
+YLabel=$N \, \text{d}{\sigma}/\text{d}{y_S}$
+FullRange=1
+LegendXPos=0.5
+LegendYPos=0.2
+# END PLOT
+
 # BEGIN PLOT /ALEPH_2004_S5765862/d54-x01-y01
 Title=Thrust ($E_\mathrm{CMS}=91.2$ GeV)
 XLabel=$T$
@@ -5,139 +397,1187 @@
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d62-x01-y01
-Title=Heavy jet mass ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$\rho$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+# BEGIN PLOT /ALEPH_2004_S5765862/d55-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d56-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d57-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d58-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d59-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d60-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d61-x01-y01
+Title=Thrust ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$1-T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}(1-T)$
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d62-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=0.7
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d63-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=1.1
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d64-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d65-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d66-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d67-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d68-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d69-x01-y01
+Title=Heavy jet mass ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\rho$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\rho$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d70-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d71-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d72-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d73-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d74-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d75-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d76-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d77-x01-y01
+Title=Total jet broadening ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$B_T$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d78-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d79-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d80-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d81-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d82-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d83-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d84-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d85-x01-y01
+Title=Wide jet broadening ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$B_W$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d86-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d87-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d88-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d89-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d90-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d91-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d92-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d93-x01-y01
+Title=C-Parameter ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$C$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
+LegendXPos=0.6
+LegendYPos=0.3
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d94-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d95-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d96-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d97-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d98-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d99-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d100-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d101-x01-y01
+Title=Thrust major ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$T_\mathrm{major}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d102-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$\ln(T_\mathrm{minor})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(T_\mathrm{minor})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d103-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d104-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d105-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d106-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d107-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d108-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d109-x01-y01
+Title=Thrust minor ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$T_\mathrm{minor}$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{minor}$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d110-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d111-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d112-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d113-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=1722$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d114-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d115-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d116-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d117-x01-y01
+Title=Jet mass difference ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$M_D$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d118-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d119-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d120-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d121-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d122-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d123-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d124-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d125-x01-y01
+Title=Aplanarity ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$A$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d126-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d127-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d128-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d129-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d130-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d131-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d132-x01-y01
+Title=Planarity ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$P$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}P$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d133-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d134-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d135-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d136-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d137-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d138-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d139-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d140-x01-y01
+Title=Oblateness ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$O$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d141-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d142-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d143-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d144-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d145-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d146-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d147-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d148-x01-y01
+Title=Sphericity ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$S$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d149-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d150-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d151-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d152-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=173$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d153-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d154-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d155-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d156-x01-y01
+Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$-\ln(y_{12})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
+LegendXPos=1.1
+LegendYPos=0.8
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d157-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d158-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d159-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d160-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d161-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d162-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d163-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d164-x01-y01
+Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$-\ln(y_{23})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
+LegendXPos=0.6
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d165-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d166-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d167-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d168-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d169-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d170-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d171-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d172-x01-y01
+Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$-\ln(y_{34})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d173-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d174-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d175-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d176-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d177-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d178-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d179-x01-y01
+Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$-\ln(y_{45})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d180-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d181-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d182-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d183-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d184-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d185-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d186-x01-y01
+Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$-\ln(y_{56})$
+YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+LegendXPos=0.7
+LegendYPos=0.4
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d187-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$\ln(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d188-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
 LegendXPos=0.5
-LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d70-x01-y01
-Title=Total jet broadening ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$B_T$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_T$
-LegendXPos=0.6
-LegendYPos=0.6
+# BEGIN PLOT /ALEPH_2004_S5765862/d189-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d78-x01-y01
-Title=Wide jet broadening ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$B_W$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}B_W$
-LegendXPos=0.6
+# BEGIN PLOT /ALEPH_2004_S5765862/d190-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d191-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d192-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d193-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d194-x01-y01
+Title=1-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.8
+LegendXPos=0.5
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d195-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$\ln(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d86-x01-y01
-Title=C-Parameter ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$C$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}C$
-LegendXPos=0.6
+# BEGIN PLOT /ALEPH_2004_S5765862/d196-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d94-x01-y01
-Title=Thrust major ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$T_\mathrm{major}$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}T_\mathrm{major}$
-LegendXPos=0.6
+# BEGIN PLOT /ALEPH_2004_S5765862/d197-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d102-x01-y01
-Title=Thrust minor ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$\ln(T_\mathrm{minor})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(T_\mathrm{minor})$
-LegendXPos=0.8
+# BEGIN PLOT /ALEPH_2004_S5765862/d198-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d110-x01-y01
-Title=Jet mass difference ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$M_D$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}M_D$
-LegendXPos=0.5
+# BEGIN PLOT /ALEPH_2004_S5765862/d199-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d118-x01-y01
-Title=Aplanarity ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$A$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}A$
+# BEGIN PLOT /ALEPH_2004_S5765862/d200-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d133-x01-y01
-Title=Oblateness ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$O$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}O$
-LegendXPos=0.5
+# BEGIN PLOT /ALEPH_2004_S5765862/d201-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d141-x01-y01
-Title=Sphericity ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$S$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}S$
+# BEGIN PLOT /ALEPH_2004_S5765862/d202-x01-y01
+Title=2-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
+LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d149-x01-y01
-Title=Durham jet resolution $2 \to 1$ ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$-\ln(y_{12})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{12})$
-LegendXPos=0.5
+# BEGIN PLOT /ALEPH_2004_S5765862/d203-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
+XLabel=$\ln(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d157-x01-y01
-Title=Durham jet resolution $3 \to 2$ ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$-\ln(y_{23})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{23})$
-LegendXPos=0.6
+# BEGIN PLOT /ALEPH_2004_S5765862/d204-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d165-x01-y01
-Title=Durham jet resolution $4 \to 3$ ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$-\ln(y_{34})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{34})$
+# BEGIN PLOT /ALEPH_2004_S5765862/d205-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
 LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d173-x01-y01
-Title=Durham jet resolution $5 \to 4$ ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$-\ln(y_{45})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{45})$
+# BEGIN PLOT /ALEPH_2004_S5765862/d206-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
 LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d180-x01-y01
-Title=Durham jet resolution $6 \to 5$ ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$-\ln(y_{56})$
-YLabel=$1/\sigma \; \mathrm{d}\sigma/\mathrm{d}\ln(y_{56})$
+# BEGIN PLOT /ALEPH_2004_S5765862/d207-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
 LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d187-x01-y01
-Title=1-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$\ln(y_\mathrm{cut})$
-YLabel=$\sigma(\text{1 jet})/\sigma(\text{inclusive})$
-LegendYPos=0.5
+# BEGIN PLOT /ALEPH_2004_S5765862/d208-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.7
+LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d195-x01-y01
-Title=2-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$\ln(y_\mathrm{cut})$
-YLabel=$\sigma(\text{2 jet})/\sigma(\text{inclusive})$
+# BEGIN PLOT /ALEPH_2004_S5765862/d209-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.7
 LegendYPos=0.6
 # END PLOT
 
-# BEGIN PLOT /ALEPH_2004_S5765862/d203-x01-y01
-Title=3-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
-XLabel=$\ln(y_\mathrm{cut})$
+# BEGIN PLOT /ALEPH_2004_S5765862/d210-x01-y01
+Title=3-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
 YLabel=$\sigma(\text{3 jet})/\sigma(\text{inclusive})$
 LegendXPos=0.7
 LegendYPos=0.6
@@ -151,6 +1591,62 @@
 LegendYPos=0.6
 # END PLOT
 
+# BEGIN PLOT /ALEPH_2004_S5765862/d212-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d213-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d214-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d215-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d216-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d217-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d218-x01-y01
+Title=4-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{4 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
 # BEGIN PLOT /ALEPH_2004_S5765862/d219-x01-y01
 Title=5-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
 XLabel=$\ln(y_\mathrm{cut})$
@@ -159,6 +1655,62 @@
 LegendYPos=0.6
 # END PLOT
 
+# BEGIN PLOT /ALEPH_2004_S5765862/d220-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d221-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d222-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d223-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d224-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d225-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d226-x01-y01
+Title=5-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma(\text{5 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.8
+LegendYPos=0.6
+# END PLOT
+
 # BEGIN PLOT /ALEPH_2004_S5765862/d227-x01-y01
 Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=91.2$ GeV)
 XLabel=$\ln(y_\mathrm{cut})$
@@ -166,3 +1718,59 @@
 LegendXPos=0.5
 LegendYPos=0.6
 # END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d228-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=133$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d229-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=161$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d230-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=172$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d231-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=183$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d232-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=189$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d233-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=200$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT
+
+# BEGIN PLOT /ALEPH_2004_S5765862/d234-x01-y01
+Title=$\geq$6-jet fraction ($E_\mathrm{CMS}=206$ GeV)
+XLabel=$\log(y_\mathrm{cut})$
+YLabel=$\sigma($\geq$\text{6 jet})/\sigma(\text{inclusive})$
+LegendXPos=0.5
+LegendYPos=0.6
+# END PLOT

Modified: trunk/src/Analyses/ALEPH_2004_S5765862.cc
==============================================================================
--- trunk/src/Analyses/ALEPH_2004_S5765862.cc	Mon Jul 12 23:28:32 2010	(r2582)
+++ trunk/src/Analyses/ALEPH_2004_S5765862.cc	Tue Jul 13 11:54:36 2010	(r2583)
@@ -9,6 +9,7 @@
 #include "Rivet/Projections/Sphericity.hh"
 #include "Rivet/Projections/ParisiTensor.hh"
 #include "Rivet/Projections/Hemispheres.hh"
+#include "Rivet/Projections/Beam.hh"
 
 namespace Rivet {
 
@@ -18,7 +19,8 @@
   public:
 
     ALEPH_2004_S5765862()
-      : Analysis("ALEPH_2004_S5765862") , _initialised(false)
+      : Analysis("ALEPH_2004_S5765862") , _initialisedJets(false),
+	_initialisedSpectra(false), _weightedTotalChargedPartNum(0)
     {
       setBeams(ELECTRON, POSITRON);
     }
@@ -27,8 +29,8 @@
   public:
 
     void init() {
-      _initialised=true;
-
+      _initialisedJets    = true;
+      _initialisedSpectra = true;
       const FinalState fs;
       addProjection(fs, "FS");
       addProjection(FastJets(fs, FastJets::DURHAM, 0.7), "DurhamJets");
@@ -39,7 +41,14 @@
       addProjection(ParisiTensor(fs), "Parisi");
       addProjection(Hemispheres(thrust), "Hemispheres");
 
+      const ChargedFinalState cfs;
+      addProjection(Beam(), "Beams");
+      addProjection(cfs, "CFS");
+      addProjection(Thrust(cfs), "ChargedThrust");
+      addProjection(Sphericity(cfs), "ChargedSphericity");
+
       // Histos
+      // offset for the event shapes and jets
       int offset = 0;
       switch (int(sqrtS()/GeV + 0.5)) {
       case 91: offset = 0; break;
@@ -51,170 +60,243 @@
       case 200: offset = 6; break;
       case 206: offset = 7; break;
       default:
+	_initialisedJets = false;
+      }
+      // event shapes
+      if(_initialisedJets) {
+	_h_thrust = bookHistogram1D(offset+54, 1, 1);
+	_h_heavyjetmass = bookHistogram1D(offset+62, 1, 1);
+	_h_totaljetbroadening = bookHistogram1D(offset+70, 1, 1);
+	_h_widejetbroadening = bookHistogram1D(offset+78, 1, 1);
+	_h_cparameter = bookHistogram1D(offset+86, 1, 1);
+	_h_thrustmajor = bookHistogram1D(offset+94, 1, 1);
+	_h_thrustminor = bookHistogram1D(offset+102, 1, 1);
+	_h_jetmassdifference = bookHistogram1D(offset+110, 1, 1);
+	_h_aplanarity = bookHistogram1D(offset+118, 1, 1);
+	_h_planarity  = offset==0 ? NULL : bookHistogram1D(offset+125, 1, 1);
+	_h_oblateness = bookHistogram1D(offset+133, 1, 1);
+	_h_sphericity = bookHistogram1D(offset+141, 1, 1);
+	
+	// Durham n->m jet resolutions
+	_h_y_Durham[0] = bookHistogram1D(offset+149, 1, 1);   // y12 d149 ... d156
+	_h_y_Durham[1] = bookHistogram1D(offset+157, 1, 1);   // y23 d157 ... d164
+	if (offset<6) { // there is no y34, y45 and y56 for 200 gev
+	  _h_y_Durham[2] = bookHistogram1D(offset+165, 1, 1); // y34 d165 ... d172, but not 171
+	  _h_y_Durham[3] = bookHistogram1D(offset+173, 1, 1); // y45 d173 ... d179
+	  _h_y_Durham[4] = bookHistogram1D(offset+180, 1, 1); // y56 d180 ... d186
+	}
+	else if (offset==6) {
+	  _h_y_Durham[2] = NULL;
+	  _h_y_Durham[3] = NULL;
+	  _h_y_Durham[4] = NULL;
+	}
+	else if (offset==7) {
+	  _h_y_Durham[2] = bookHistogram1D(172, 1, 1);
+	  _h_y_Durham[3] = bookHistogram1D(179, 1, 1);
+	  _h_y_Durham[4] = bookHistogram1D(186, 1, 1);
+	}
+	
+	// Durham n-jet fractions
+	_h_R_Durham[0] = bookDataPointSet(offset+187, 1, 1); // R1 d187 ... d194
+	_h_R_Durham[1] = bookDataPointSet(offset+195, 1, 1); // R2 d195 ... d202
+	_h_R_Durham[2] = bookDataPointSet(offset+203, 1, 1); // R3 d203 ... d210
+	_h_R_Durham[3] = bookDataPointSet(offset+211, 1, 1); // R4 d211 ... d218
+	_h_R_Durham[4] = bookDataPointSet(offset+219, 1, 1); // R5 d219 ... d226
+	_h_R_Durham[5] = bookDataPointSet(offset+227, 1, 1); // R>=6 d227 ... d234
+      }
+      // offset for the charged particle distributions
+      offset = 0;
+      switch (int(sqrtS()/GeV + 0.5)) {
+      case 133: offset = 0; break;
+      case 161: offset = 1; break;
+      case 172: offset = 2; break;
+      case 183: offset = 3; break;
+      case 189: offset = 4; break;
+      case 196: offset = 5; break;
+      case 200: offset = 6; break;
+      case 206: offset = 7; break;
+      default:
+	_initialisedSpectra=false;
+      }
+      if(_initialisedSpectra) {
+	_h_xp = bookHistogram1D( 2+offset, 1, 1);
+	_h_xi = bookHistogram1D(11+offset, 1, 1);
+	_h_xe = bookHistogram1D(19+offset, 1, 1);
+	_h_pTin  = bookHistogram1D(27+offset, 1, 1);
+	_h_pTout = offset!=7 ? NULL : bookHistogram1D(35, 1, 1);
+	_h_rapidityT = bookHistogram1D(36+offset, 1, 1);
+	_h_rapidityS = bookHistogram1D(44+offset, 1, 1);
+      }
+
+      if(!_initialisedSpectra && !_initialisedJets) {
         getLog() << Log::WARNING
                  << "CMS energy of events sqrt(s) = " << sqrtS()/GeV
-                 <<" doesn't match any available analysis energy." << endl;
-        _initialised=false;
-        return;
+                 <<" doesn't match any available analysis energy ." << endl;
       }
-
-      // event shapes
-      _h_thrust = bookHistogram1D(offset+54, 1, 1);
-      _h_heavyjetmass = bookHistogram1D(offset+62, 1, 1);
-      _h_totaljetbroadening = bookHistogram1D(offset+70, 1, 1);
-      _h_widejetbroadening = bookHistogram1D(offset+78, 1, 1);
-      _h_cparameter = bookHistogram1D(offset+86, 1, 1);
-      _h_thrustmajor = bookHistogram1D(offset+94, 1, 1);
-      _h_thrustminor = bookHistogram1D(offset+102, 1, 1);
-      _h_jetmassdifference = bookHistogram1D(offset+110, 1, 1);
-      _h_aplanarity = bookHistogram1D(offset+118, 1, 1);
-      // planarity is missing the 91 gev data, so left out
-      _h_oblateness = bookHistogram1D(offset+133, 1, 1);
-      _h_sphericity = bookHistogram1D(offset+141, 1, 1);
-
-
-      // Durham n->m jet resolutions
-      _h_y_Durham[0] = bookHistogram1D(offset+149, 1, 1);   // y12 d149 ... d156
-      _h_y_Durham[1] = bookHistogram1D(offset+157, 1, 1);   // y23 d157 ... d164
-      if (offset<6) { // there is no y34, y45 and y56 for 200 gev
-        _h_y_Durham[2] = bookHistogram1D(offset+165, 1, 1); // y34 d165 ... d172, but not 171
-        _h_y_Durham[3] = bookHistogram1D(offset+173, 1, 1); // y45 d173 ... d179
-        _h_y_Durham[4] = bookHistogram1D(offset+180, 1, 1); // y56 d180 ... d186
-      }
-      else if (offset==6) {
-        _h_y_Durham[2] = NULL;
-        _h_y_Durham[3] = NULL;
-        _h_y_Durham[4] = NULL;
-      }
-      else if (offset==7) {
-        _h_y_Durham[2] = bookHistogram1D(172, 1, 1);
-        _h_y_Durham[3] = bookHistogram1D(179, 1, 1);
-        _h_y_Durham[4] = bookHistogram1D(186, 1, 1);
-      }
-
-      // Durham n-jet fractions
-      _h_R_Durham[0] = bookDataPointSet(offset+187, 1, 1); // R1 d187 ... d194
-      _h_R_Durham[1] = bookDataPointSet(offset+195, 1, 1); // R2 d195 ... d202
-      _h_R_Durham[2] = bookDataPointSet(offset+203, 1, 1); // R3 d203 ... d210
-      _h_R_Durham[3] = bookDataPointSet(offset+211, 1, 1); // R4 d211 ... d218
-      _h_R_Durham[4] = bookDataPointSet(offset+219, 1, 1); // R5 d219 ... d226
-      _h_R_Durham[5] = bookDataPointSet(offset+227, 1, 1); // R>=6 d227 ... d234
     }
 
 
     void analyze(const Event& e) {
-      if (!_initialised) return;
       const double weight = e.weight();
+      if(_initialisedJets) {
+	bool LEP1 = fuzzyEquals(sqrtS(),91.2*GeV,0.01);
+	// event shapes
+	const Thrust& thrust = applyProjection<Thrust>(e, "Thrust");
+	double thr = LEP1 ? thrust.thrust() : 1.0 - thrust.thrust();
+	_h_thrust->fill(thr,weight);
+	_h_thrustmajor->fill(thrust.thrustMajor(),weight);
+	if(LEP1)
+	  _h_thrustminor->fill(log(thrust.thrustMinor()),weight);
+	else
+	  _h_thrustminor->fill(thrust.thrustMinor(),weight);
+	_h_oblateness->fill(thrust.oblateness(),weight);
+	
+	const Hemispheres& hemi = applyProjection<Hemispheres>(e, "Hemispheres");
+	_h_heavyjetmass->fill(hemi.scaledM2high(),weight);
+	_h_jetmassdifference->fill(hemi.scaledM2diff(),weight);
+	_h_totaljetbroadening->fill(hemi.Bsum(),weight);
+	_h_widejetbroadening->fill(hemi.Bmax(),weight);
+	
+	const ParisiTensor& parisi = applyProjection<ParisiTensor>(e, "Parisi");
+	_h_cparameter->fill(parisi.C(),weight);
+	
+	const Sphericity& sphericity = applyProjection<Sphericity>(e, "Sphericity");
+	_h_aplanarity->fill(sphericity.aplanarity(),weight);
+	if(_h_planarity)
+	  _h_planarity->fill(sphericity.planarity(),weight);
+	_h_sphericity->fill(sphericity.sphericity(),weight);
+	
+	// jet rates
+	const FastJets& durjet = applyProjection<FastJets>(e, "DurhamJets");
+	double log10e = log10(exp(1.));
+	if (durjet.clusterSeq()) {
+	  double logynm1=0.;
+	  double logyn;
+	  for (size_t i=0; i<5; ++i) {
+	    logyn = -log(durjet.clusterSeq()->exclusive_ymerge_max(i+1));
+	    if (_h_y_Durham[i]) {
+	      _h_y_Durham[i]->fill(logyn, weight);
+	    }
+	    if(!LEP1) logyn *= log10e;
+	    for (size_t j = 0; j < _h_R_Durham[i]->size(); ++j) {
+	      IDataPoint* dp = _h_R_Durham[i]->point(j);
+	      double val = -dp->coordinate(0)->value()+dp->coordinate(0)->errorMinus();
+	      if(val<=logynm1) break;
+	      if(val<logyn) {
+		dp->coordinate(1)->setValue(dp->coordinate(1)->value()+weight);
+	      }
+	    }
+	    logynm1 = logyn;
+	  }
+	  for (size_t j = 0; j < _h_R_Durham[5]->size(); ++j) {
+	    IDataPoint* dp = _h_R_Durham[5]->point(j);
+	    double val = -dp->coordinate(0)->value()+dp->coordinate(0)->errorMinus();
+	    if(val<=logynm1) break;
+	    dp->coordinate(1)->setValue(dp->coordinate(1)->value()+weight);
+	  }
+	}
+	if( !_initialisedSpectra) {
+	  const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(e, "CFS");
+	  const size_t numParticles = cfs.particles().size();
+	  _weightedTotalChargedPartNum += numParticles * weight;
+	}
+      }
 
-      // event shapes
-      const Thrust& thrust = applyProjection<Thrust>(e, "Thrust");
-      double thr = (fuzzyEquals(sqrtS(),91.2*GeV,0.5)?thrust.thrust():1.0-thrust.thrust());
-      _h_thrust->fill(thr,weight);
-      _h_thrustmajor->fill(thrust.thrustMajor(),weight);
-      _h_thrustminor->fill(log(thrust.thrustMinor()),weight);
-      _h_oblateness->fill(thrust.oblateness(),weight);
-
-      const Hemispheres& hemi = applyProjection<Hemispheres>(e, "Hemispheres");
-      _h_heavyjetmass->fill(hemi.scaledM2high(),weight);
-      _h_jetmassdifference->fill(hemi.scaledM2diff(),weight);
-      _h_totaljetbroadening->fill(hemi.Bsum(),weight);
-      _h_widejetbroadening->fill(hemi.Bmax(),weight);
-
-      const ParisiTensor& parisi = applyProjection<ParisiTensor>(e, "Parisi");
-      _h_cparameter->fill(parisi.C(),weight);
-
-      const Sphericity& sphericity = applyProjection<Sphericity>(e, "Sphericity");
-      _h_aplanarity->fill(sphericity.aplanarity(),weight);
-      _h_sphericity->fill(sphericity.sphericity(),weight);
-
-      // jet rates
-      const FastJets& durjet = applyProjection<FastJets>(e, "DurhamJets");
-      if (durjet.clusterSeq()) {
-        for (size_t i=0; i<5; ++i) {
-          if (_h_y_Durham[i]) {
-            _h_y_Durham[i]->fill(-log(durjet.clusterSeq()->exclusive_ymerge_max(i+1)), weight);
-          }
-        }
+      // charged particle distributions
+      if(_initialisedSpectra) {
+	const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(e, "CFS");
+	const size_t numParticles = cfs.particles().size();
+	_weightedTotalChargedPartNum += numParticles * weight;
+	const ParticlePair& beams = applyProjection<Beam>(e, "Beams").beams();
+	const double meanBeamMom = ( beams.first.momentum().vector3().mod() +
+				     beams.second.momentum().vector3().mod() ) / 2.0;
+	const Thrust& cthrust = applyProjection<Thrust>(e, "ChargedThrust");
+	const Sphericity& csphere = applyProjection<Sphericity>(e, "ChargedSphericity");
+	foreach (const Particle& p, cfs.particles()) {
+	  const double xp = p.momentum().vector3().mod()/meanBeamMom;
+	  _h_xp->fill(xp   , weight);
+	  const double logxp = -std::log(xp);
+	  _h_xi->fill(logxp, weight);
+	  const double xe = p.momentum().E()/meanBeamMom;
+	  _h_xe->fill(xe   , weight);
+	  const double pTinT  = dot(p.momentum().vector3(), cthrust.thrustMajorAxis());
+	  const double pToutT = dot(p.momentum().vector3(), cthrust.thrustMinorAxis());
+	  _h_pTin->fill(fabs(pTinT/GeV), weight);
+	  if(_h_pTout) _h_pTout->fill(fabs(pToutT/GeV), weight);
+	  const double momT = dot(cthrust.thrustAxis()        ,p.momentum().vector3());
+	  const double rapidityT = 0.5 * std::log((p.momentum().E() + momT) / 
+						  (p.momentum().E() - momT));
+	  _h_rapidityT->fill(rapidityT, weight);
+	  const double momS = dot(csphere.sphericityAxis(),p.momentum().vector3());
+	  const double rapidityS = 0.5 * std::log((p.momentum().E() + momS) / 
+						  (p.momentum().E() - momS));
+	  _h_rapidityS->fill(rapidityS, weight);
+	}
       }
     }
 
-
     void finalize() {
-      if (!_initialised) return;
+      if(!_initialisedJets && !_initialisedSpectra) return;
 
-      normalize(_h_thrust);
-      normalize(_h_heavyjetmass);
-      normalize(_h_totaljetbroadening);
-      normalize(_h_widejetbroadening);
-      normalize(_h_cparameter);
-      normalize(_h_thrustmajor);
-      normalize(_h_thrustminor);
-      normalize(_h_jetmassdifference);
-      normalize(_h_aplanarity);
-      normalize(_h_oblateness);
-      normalize(_h_sphericity);
-
-      for (int N=1; N<7; ++N) {
-        // calculate the N jet fraction from the jet resolution histograms
-
-        for (size_t i = 0; i < _h_R_Durham[N-1]->size(); ++i) {
-          IDataPoint* dp = _h_R_Durham[N-1]->point(i);
-          // get ycut at which the njet-fraction is to be calculated
-          /// @todo HepData has binwidths here, which doesn't make sense at all
-          /// I assume the low edge to be the ycut
-          double ycut = dp->coordinate(0)->value()-dp->coordinate(0)->errorMinus();
-
-          // sum all >=N jet events
-          double sigmaNinclusive = 0.0;
-          if (N>1) {
-            if (_h_y_Durham[N-2]) {
-              AIDA::IHistogram1D* y_Nminus1_N = _h_y_Durham[N-2];
-              // watch out, y_NM is negatively binned
-              int cutbin=y_Nminus1_N->coordToIndex(-ycut);
-              if (cutbin==AIDA::IAxis::UNDERFLOW_BIN) cutbin=0;
-              if (cutbin==AIDA::IAxis::OVERFLOW_BIN) cutbin=y_Nminus1_N->axis().bins()-1;
-              for (int ibin=0; ibin<cutbin; ++ibin) {
-                sigmaNinclusive += y_Nminus1_N->binHeight(ibin);
-              }
-            }
-          }
-          else sigmaNinclusive = sumOfWeights();
-
-          // sum all >=N+1 jet events
-          double sigmaNplus1inclusive = 0.0;
-          if (N<6) {
-            if (_h_y_Durham[N-1]) {
-              AIDA::IHistogram1D* y_N_Nplus1 = _h_y_Durham[N-1];
-              // watch out, y_NM is negatively binned
-              int cutbin=y_N_Nplus1->coordToIndex(-ycut);
-              if (cutbin==AIDA::IAxis::UNDERFLOW_BIN) cutbin=0;
-              if (cutbin==AIDA::IAxis::OVERFLOW_BIN) cutbin=y_N_Nplus1->axis().bins();
-              for (int ibin=0; ibin<cutbin; ++ibin) {
-                sigmaNplus1inclusive += y_N_Nplus1->binHeight(ibin);
-              }
-            }
-          }
-
-          // njetfraction = (sigma(>=N jet) - sigma(>=N+1 jet)) / sigma(tot)
-          double njetfraction = (sigmaNinclusive-sigmaNplus1inclusive)/sumOfWeights();
-          dp->coordinate(1)->setValue(njetfraction);
-        }
-      }
-
-
-      for (size_t n = 0; n < 5; ++n) {
-        if (_h_y_Durham[n]) {
-          scale(_h_y_Durham[n], 1.0/sumOfWeights());
-        }
+      if (_initialisedJets) {
+	normalize(_h_thrust);
+	normalize(_h_heavyjetmass);
+	normalize(_h_totaljetbroadening);
+	normalize(_h_widejetbroadening);
+	normalize(_h_cparameter);
+	normalize(_h_thrustmajor);
+	normalize(_h_thrustminor);
+	normalize(_h_jetmassdifference);
+	normalize(_h_aplanarity);
+	if(_h_planarity) normalize(_h_planarity);
+	normalize(_h_oblateness);
+	normalize(_h_sphericity);
+	
+	for (int N=1; N<7; ++N) {
+	  for (size_t i = 0; i < _h_R_Durham[N-1]->size(); ++i) {
+	    _h_R_Durham[N-1]->point(i)->coordinate(1)->
+	      setValue(_h_R_Durham[N-1]->point(i)->coordinate(1)->value()/sumOfWeights());
+	  }
+	}
+	
+	for (size_t n = 0; n < 5; ++n) {
+	  if (_h_y_Durham[n]) {
+	    scale(_h_y_Durham[n], 1.0/sumOfWeights());
+	  }
+	}
       }
 
-    }
+      const double avgNumParts = _weightedTotalChargedPartNum / sumOfWeights();
+      AIDA::IDataPointSet * mult = bookDataPointSet(1, 1, 1);
+      for(size_t i = 0; i < mult->size(); ++i) {
+	if(fuzzyEquals(sqrtS(),mult->point(i)->coordinate(0)->value(),0.01)) {
+	  mult->point(i)->coordinate(1)->setValue(avgNumParts);
+	}
+      }
 
+      if(_initialisedSpectra) {
+	normalize(_h_xp, avgNumParts);
+	normalize(_h_xi, avgNumParts);
+	normalize(_h_xe, avgNumParts);
+	normalize(_h_pTin , avgNumParts);
+	if(_h_pTout) normalize(_h_pTout, avgNumParts);
+	normalize(_h_rapidityT, avgNumParts);
+	normalize(_h_rapidityS, avgNumParts);
+      }
+    }
 
   private:
 
-    bool _initialised;
+    bool _initialisedJets;
+    bool _initialisedSpectra;
 
+    AIDA::IHistogram1D *_h_xp;
+    AIDA::IHistogram1D *_h_xi;
+    AIDA::IHistogram1D *_h_xe;
+    AIDA::IHistogram1D *_h_pTin;
+    AIDA::IHistogram1D *_h_pTout;
+    AIDA::IHistogram1D *_h_rapidityT;
+    AIDA::IHistogram1D *_h_rapidityS;
     AIDA::IHistogram1D *_h_thrust;
     AIDA::IHistogram1D *_h_heavyjetmass;
     AIDA::IHistogram1D *_h_totaljetbroadening;
@@ -224,12 +306,15 @@
     AIDA::IHistogram1D *_h_thrustminor;
     AIDA::IHistogram1D *_h_jetmassdifference;
     AIDA::IHistogram1D *_h_aplanarity;
+    AIDA::IHistogram1D *_h_planarity;
     AIDA::IHistogram1D *_h_oblateness;
     AIDA::IHistogram1D *_h_sphericity;
 
     AIDA::IDataPointSet *_h_R_Durham[6];
     AIDA::IHistogram1D *_h_y_Durham[5];
 
+    double _weightedTotalChargedPartNum;
+
   };
 
 


More information about the Rivet-svn mailing list