|
[Rivet-svn] r2583 - in trunk: data/plotinfo src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgTue 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 |