|
|
[Rivet-svn] r3622 - in branches/2011-07-aida2yoda: . bin data/anainfo data/plotinfo data/refdata include/Rivet src/Analyses src/Core src/Projections
blackhole at projects.hepforge.org
blackhole at projects.hepforge.org
Fri Mar 16 18:22:33 GMT 2012
Author: hoeth
Date: Fri Mar 16 18:22:32 2012
New Revision: 3622
Log:
merge c3574-3583 from trunk
Added:
branches/2011-07-aida2yoda/data/anainfo/LHCB_2011_I917009.info
- copied unchanged from r3583, trunk/data/anainfo/LHCB_2011_I917009.info
branches/2011-07-aida2yoda/data/plotinfo/LHCB_2011_I917009.plot
- copied, changed from r3583, trunk/data/plotinfo/LHCB_2011_I917009.plot
branches/2011-07-aida2yoda/data/refdata/LHCB_2011_I917009.aida
- copied unchanged from r3583, trunk/data/refdata/LHCB_2011_I917009.aida
branches/2011-07-aida2yoda/src/Analyses/LHCB_2011_I917009.cc
- copied, changed from r3583, trunk/src/Analyses/LHCB_2011_I917009.cc
Modified:
branches/2011-07-aida2yoda/ChangeLog
branches/2011-07-aida2yoda/bin/rivet-rescale
branches/2011-07-aida2yoda/data/anainfo/Makefile.am
branches/2011-07-aida2yoda/data/plotinfo/JADE_1998_S3612880.plot
branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
branches/2011-07-aida2yoda/data/refdata/Makefile.am
branches/2011-07-aida2yoda/include/Rivet/Analysis.hh
branches/2011-07-aida2yoda/src/Analyses/ATLAS_2011_I925932.cc
branches/2011-07-aida2yoda/src/Analyses/Makefile.am
branches/2011-07-aida2yoda/src/Core/Analysis.cc
branches/2011-07-aida2yoda/src/Projections/InitialQuarks.cc
Modified: branches/2011-07-aida2yoda/ChangeLog
==============================================================================
--- branches/2011-07-aida2yoda/ChangeLog Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/ChangeLog Fri Mar 16 18:22:32 2012 (r3622)
@@ -1,3 +1,18 @@
+2012-02-26 Andy Buckley <andy.buckley at cern.ch>
+
+ * Adding LHCB_2011_I917009 analysis, from Alex Grecu.
+
+ * src/Core/Analysis.cc, include/Rivet/Analysis.hh: Add a numeric-arg version of histoPath().
+
+2012-02-24 Holger Schulz <holger.schulz at physik.hu-berlin.de>
+
+ * Adding ATLAS Ks/Lambda analysis.
+
+2012-02-20 Andy Buckley <andy.buckley at cern.ch>
+
+ * src/Analyses/ATLAS_2011_I925932.cc: Using new overflow-aware
+ normalize() in place of counters and scale(..., 1/count)
+
2012-02-14 Andy Buckley <andy.buckley at cern.ch>
* Splitting MC_GENERIC to put the PDF and PID plotting into
Modified: branches/2011-07-aida2yoda/bin/rivet-rescale
==============================================================================
--- branches/2011-07-aida2yoda/bin/rivet-rescale Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/bin/rivet-rescale Fri Mar 16 18:22:32 2012 (r3622)
@@ -316,8 +316,12 @@
else:
newarea = tempref.getArea()
- scalefactor = newarea/oldarea
- if scalefactor != 1.0:
+ if not oldarea == 0.0:
+ scalefactor = newarea/oldarea
+ else:
+ scalefactor= 0.0
+ if scalefactor != 1.0 and scalefactor > 0.0:
+ print scalefactor
oldarea = histos[name].getArea()
newarea = histos[name].getArea() * scalefactor
Copied: branches/2011-07-aida2yoda/data/anainfo/LHCB_2011_I917009.info (from r3583, trunk/data/anainfo/LHCB_2011_I917009.info)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/anainfo/LHCB_2011_I917009.info Fri Mar 16 18:22:32 2012 (r3622, copy of r3583, trunk/data/anainfo/LHCB_2011_I917009.info)
@@ -0,0 +1,43 @@
+Name: LHCB_2011_I917009
+Year: 2011
+Summary: $V^0$ production ratios in $pp$ collisions at $\sqrt{s}= 0.9$ and 7 TeV at LHCb
+Experiment: LHCb
+Collider: LHC
+InspireID: 917009
+Status: VALIDATED
+Authors:
+ - Alex Grecu <Alex.Grecu at cern.ch>
+References:
+ - JHEP08:034,2011
+ - arXiv:1107.0882[hep-ex]
+ - doi:10.1007/JHEP08(2011)034
+RunInfo:
+ QCD events. LHCb minimum bias, Perugia 0 and Perugia NOCR tune events used for reproducing published histograms.
+NumEvents: 5000000
+Beams: [p+, p+]
+Energies: [900, 7000]
+# PtCuts: [[0.25, 2.5], [0.15,2.5]]
+Description:
+ This paper presents the production ratios for $\bar{\Lambda}/\Lambda$ and
+ $\bar{\Lambda}/K_{s}^{0}$ measured by LHCb detector in 2010 at $\sqrt{s} = 0.9$
+ TeV and 7 TeV as functions of the transverse momentum $p_\perp$ and the rapidity $y$
+ in the ranges $0.15 < p_\perp < \unit{2.50}{GeV/c}$ and $2.0 < y < 4.5$, respectively.
+ The results for the two energies are merged and represented as a function of
+ rapidity loss $\Delta y = y_{\mathrm{beam}} - y$.
+NeedsCrossSection: true
+BibKey: Aaij:2011va
+BibTeX: '@article{Aaij:2011va,
+ author = "Aaij, R. and others",
+ title = "{Measurement of $V^0$ production ratios in $pp$
+ collisions at $\sqrt{s} = 0.9$ and 7\,TeV}",
+ collaboration = "LHCb Collaboration",
+ journal = "JHEP",
+ volume = "1108",
+ pages = "034",
+ doi = "10.1007/JHEP08(2011)034",
+ year = "2011",
+ note = "* Temporary entry *",
+ eprint = "1107.0882",
+ archivePrefix = "arXiv",
+ primaryClass = "hep-ex",
+}'
Modified: branches/2011-07-aida2yoda/data/anainfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/data/anainfo/Makefile.am Fri Mar 16 18:22:32 2012 (r3622)
@@ -116,6 +116,7 @@
JADE_OPAL_2000_S4300807.info \
JADE_1998_S3612880.info \
LHCB_2010_S8758301.info \
+ LHCB_2011_I917009.info \
LHCB_2011_I919315.info \
MC_DIJET.info \
MC_DIPHOTON.info \
Modified: branches/2011-07-aida2yoda/data/plotinfo/JADE_1998_S3612880.plot
==============================================================================
--- branches/2011-07-aida2yoda/data/plotinfo/JADE_1998_S3612880.plot Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/data/plotinfo/JADE_1998_S3612880.plot Fri Mar 16 18:22:32 2012 (r3622)
@@ -58,7 +58,7 @@
YLabel=$1/\sigma \text{d}\sigma/\text{d}(d_2)$
# END PLOT
-# BEGIN PLOT /JADE_1998_S3612880/d11-x01-y01
+# BEGIN PLOT /JADE_1998_S3612880/d12-x01-y01
Title=Differential 2-Jet rate (Durham), $\sqrt{s}=22$ GeV
XLabel=$D_2$
YLabel=$1/\sigma \text{d}\sigma/\text{d}(d_2)$
Copied and modified: branches/2011-07-aida2yoda/data/plotinfo/LHCB_2011_I917009.plot (from r3583, trunk/data/plotinfo/LHCB_2011_I917009.plot)
==============================================================================
--- trunk/data/plotinfo/LHCB_2011_I917009.plot Sun Feb 26 22:39:15 2012 (r3583, copy source)
+++ branches/2011-07-aida2yoda/data/plotinfo/LHCB_2011_I917009.plot Fri Mar 16 18:22:32 2012 (r3622)
@@ -1,236 +1,236 @@
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-ErrorBars=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y02
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ Gev/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.4
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y03
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y02
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y03
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMax=0.7
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d03-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.55
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d04-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendXPos=0.05
-LegendYPos=0.35
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d05-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendYPos=0.48
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d06-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.3
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d07-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-LegendYPos=0.5
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d08-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.3
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-YMax=1.2
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y02
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-LegendYPos=0.3
-LegendXPos=0.1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y03
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-LegendYPos=0.3
-LegendXPos=0.1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-LegendXPos=0.4
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y02
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
-XLabel=Rapidity (y)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.48
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y03
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.7
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d11-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity ($y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d12-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.5$ GeV/$c$)
-XLabel=Rapidity (y)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-YMin=0.1
-YMax=0.43
-FullRange=1
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d13-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d14-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
-XLabel=$p_\perp$ [GeV/$c$]
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d15-x01-y01
-Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/\Lambda$
-LogY=0
-YMin=0.25
-LegendYPos=0.65
-RatioPlot=0
-# END PLOT
-
-# BEGIN PLOT /LHCB_2011_I917009/d16-x01-y01
-Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
-XLabel=Rapidity loss ($\Delta y$)
-YLabel=$\bar{\Lambda}/K_{s}^{0}$
-LogY=0
-LegendYPos=0.33
-RatioPlot=0
-# END PLOT
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+ErrorBars=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y02
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ Gev/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.4
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d01-x01-y03
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y02
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d02-x01-y03
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMax=0.7
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d03-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.55
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d04-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.35
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d05-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendYPos=0.48
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d06-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.3
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d07-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+LegendYPos=0.5
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d08-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.3
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+YMax=1.2
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y02
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+LegendYPos=0.3
+LegendXPos=0.1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d09-x01-y03
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+LegendYPos=0.3
+LegendXPos=0.1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+LegendXPos=0.4
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y02
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/$c$)
+XLabel=Rapidity (y)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.48
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d10-x01-y03
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.7
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d11-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity ($y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d12-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.5$ GeV/$c$)
+XLabel=Rapidity (y)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+YMin=0.1
+YMax=0.43
+FullRange=1
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d13-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d14-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)
+XLabel=$p_\perp$ [GeV/$c$]
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d15-x01-y01
+Title=$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/\Lambda$
+LogY=0
+YMin=0.25
+LegendYPos=0.65
+RatioPlot=0
+# END PLOT
+
+# BEGIN PLOT /LHCB_2011_I917009/d16-x01-y01
+Title=$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/$c$)
+XLabel=Rapidity loss ($\Delta y$)
+YLabel=$\bar{\Lambda}/K_{s}^{0}$
+LogY=0
+LegendYPos=0.33
+RatioPlot=0
+# END PLOT
Modified: branches/2011-07-aida2yoda/data/plotinfo/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/data/plotinfo/Makefile.am Fri Mar 16 18:22:32 2012 (r3622)
@@ -108,6 +108,7 @@
JADE_OPAL_2000_S4300807.plot \
JADE_1998_S3612880.plot \
LHCB_2010_S8758301.plot \
+ LHCB_2011_I917009.plot \
LHCB_2011_I919315.plot \
MC_GENERIC.plot \
MC_IDENTIFIED.plot \
Copied: branches/2011-07-aida2yoda/data/refdata/LHCB_2011_I917009.aida (from r3583, trunk/data/refdata/LHCB_2011_I917009.aida)
==============================================================================
--- /dev/null 00:00:00 1970 (empty, because file is newly added)
+++ branches/2011-07-aida2yoda/data/refdata/LHCB_2011_I917009.aida Fri Mar 16 18:22:32 2012 (r3622, copy of r3583, trunk/data/refdata/LHCB_2011_I917009.aida)
@@ -0,0 +1,510 @@
+<?xml version="1.0" encoding="utf-8"?>
+<!DOCTYPE aida
+ SYSTEM 'http://aida.freehep.org/schemas/3.2.1/aida.dtd'>
+<aida version="3.2">
+ <implementation package="LHCb raw2aida" version="1.0"/>
+ <dataPointSet dimension="2" name="d01-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="4.864977e-01" errorPlus="4.864977e-01" value="1.622000e+00"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="7.249828e-02" errorPlus="7.249828e-02" value="9.040000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="5.467175e-02" errorPlus="5.467175e-02" value="6.100000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="1.348518e-01" errorPlus="1.348518e-01" value="4.200000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d01-x01-y02" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ Gev/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="1.001699e-01" errorPlus="1.001699e-01" value="7.230000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="4.579301e-02" errorPlus="4.579301e-02" value="7.720000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="5.108816e-02" errorPlus="5.108816e-02" value="7.460000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="6.657327e-02" errorPlus="6.657327e-02" value="6.170000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d01-x01-y03" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="1.164174e-01" errorPlus="1.164174e-01" value="9.040000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="5.188449e-02" errorPlus="5.188449e-02" value="7.450000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="4.601087e-02" errorPlus="4.601087e-02" value="7.570000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="4.390900e-02" errorPlus="4.390900e-02" value="4.850000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d02-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 0.65$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="4.440721e-02" errorPlus="4.440721e-02" value="1.970000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="2.607681e-02" errorPlus="2.607681e-02" value="2.180000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.059126e-02" errorPlus="2.059126e-02" value="1.800000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.744329e-02" errorPlus="3.744329e-02" value="1.580000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d02-x01-y02" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.65 < p_\perp < 1.00$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="3.828838e-02" errorPlus="3.828838e-02" value="3.160000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="2.641969e-02" errorPlus="2.641969e-02" value="3.060000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.505993e-02" errorPlus="2.505993e-02" value="3.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.041381e-02" errorPlus="3.041381e-02" value="2.990000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d02-x01-y03" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($1.00 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="5.353504e-02" errorPlus="5.353504e-02" value="4.630000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.264966e-02" errorPlus="3.264966e-02" value="4.290000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.577709e-02" errorPlus="3.577709e-02" value="4.130000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.280244e-02" errorPlus="3.280244e-02" value="3.230000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d03-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="9.436631e-02" errorPlus="9.436631e-02" value="9.340000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.535534e-02" errorPlus="3.535534e-02" value="8.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.858756e-02" errorPlus="3.858756e-02" value="7.270000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="5.060632e-02" errorPlus="5.060632e-02" value="5.390000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d04-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="3.162278e-02" errorPlus="3.162278e-02" value="2.850000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="2.213594e-02" errorPlus="2.213594e-02" value="2.630000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.184033e-02" errorPlus="2.184033e-02" value="2.580000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="2.282542e-02" errorPlus="2.282542e-02" value="2.520000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d05-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)">
+ <dataPoint>
+ <measurement errorMinus="1.250000e-01" errorPlus="1.250000e-01" value="3.750000e-01"/>
+ <measurement errorMinus="6.095900e-02" errorPlus="6.095900e-02" value="8.060000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="5.750000e-01"/>
+ <measurement errorMinus="4.816638e-02" errorPlus="4.816638e-02" value="7.310000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="7.250000e-01"/>
+ <measurement errorMinus="4.891830e-02" errorPlus="4.891830e-02" value="7.370000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="9.000000e-01"/>
+ <measurement errorMinus="4.891830e-02" errorPlus="4.891830e-02" value="7.750000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="1.100000e+00"/>
+ <measurement errorMinus="4.104875e-02" errorPlus="4.104875e-02" value="7.010000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="6.500000e-01" errorPlus="6.500000e-01" value="1.850000e+00"/>
+ <measurement errorMinus="3.905125e-02" errorPlus="3.905125e-02" value="7.450000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d06-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($2.0 < y < 4.0$)">
+ <dataPoint>
+ <measurement errorMinus="1.250000e-01" errorPlus="1.250000e-01" value="3.750000e-01"/>
+ <measurement errorMinus="1.878829e-02" errorPlus="1.878829e-02" value="1.770000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="5.750000e-01"/>
+ <measurement errorMinus="2.012461e-02" errorPlus="2.012461e-02" value="2.180000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="7.250000e-01"/>
+ <measurement errorMinus="2.507987e-02" errorPlus="2.507987e-02" value="2.840000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="9.000000e-01"/>
+ <measurement errorMinus="2.683282e-02" errorPlus="2.683282e-02" value="3.230000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="1.100000e+00"/>
+ <measurement errorMinus="2.830194e-02" errorPlus="2.830194e-02" value="3.680000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="6.500000e-01" errorPlus="6.500000e-01" value="1.850000e+00"/>
+ <measurement errorMinus="3.176476e-02" errorPlus="3.176476e-02" value="4.420000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d07-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.120000e+00"/>
+ <measurement errorMinus="5.060632e-02" errorPlus="5.060632e-02" value="5.390000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.620000e+00"/>
+ <measurement errorMinus="3.858756e-02" errorPlus="3.858756e-02" value="7.270000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.120000e+00"/>
+ <measurement errorMinus="3.535534e-02" errorPlus="3.535534e-02" value="8.000000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.620000e+00"/>
+ <measurement errorMinus="9.436631e-02" errorPlus="9.436631e-02" value="9.340000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d08-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=0.9$ TeV ($0.25 < p_\perp < 2.50$ Gev/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.120000e+00"/>
+ <measurement errorMinus="2.282542e-02" errorPlus="2.282542e-02" value="2.520000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.620000e+00"/>
+ <measurement errorMinus="2.184033e-02" errorPlus="2.184033e-02" value="2.580000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.120000e+00"/>
+ <measurement errorMinus="2.213594e-02" errorPlus="2.213594e-02" value="2.630000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.620000e+00"/>
+ <measurement errorMinus="3.162278e-02" errorPlus="3.162278e-02" value="2.850000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d09-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="1.999725e-01" errorPlus="1.999725e-01" value="8.720000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.935734e-02" errorPlus="3.935734e-02" value="9.570000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.584690e-02" errorPlus="3.584690e-02" value="9.420000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.940812e-02" errorPlus="3.940812e-02" value="8.760000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="1.328157e-01" errorPlus="1.328157e-01" value="9.000000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d09-x01-y02" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="6.580274e-02" errorPlus="6.580274e-02" value="9.740000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="4.134005e-02" errorPlus="4.134005e-02" value="9.680000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.546830e-02" errorPlus="3.546830e-02" value="9.240000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.671512e-02" errorPlus="3.671512e-02" value="8.960000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="5.280152e-02" errorPlus="5.280152e-02" value="8.620000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d09-x01-y03" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="4.468781e-02" errorPlus="4.468781e-02" value="9.870000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.758989e-02" errorPlus="3.758989e-02" value="9.660000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.534119e-02" errorPlus="3.534119e-02" value="9.280000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.623534e-02" errorPlus="3.623534e-02" value="9.030000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="4.031129e-02" errorPlus="4.031129e-02" value="7.920000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d10-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 0.65$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="4.036087e-02" errorPlus="4.036087e-02" value="1.820000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="2.617250e-02" errorPlus="2.617250e-02" value="1.910000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.507987e-02" errorPlus="2.507987e-02" value="1.850000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="2.531798e-02" errorPlus="2.531798e-02" value="1.750000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="3.354102e-02" errorPlus="3.354102e-02" value="2.070000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d10-x01-y02" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.65 < p_\perp < 1.00$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="3.269557e-02" errorPlus="3.269557e-02" value="3.200000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.059412e-02" errorPlus="3.059412e-02" value="3.280000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.828427e-02" errorPlus="2.828427e-02" value="3.150000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="2.844293e-02" errorPlus="2.844293e-02" value="2.990000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="3.138471e-02" errorPlus="3.138471e-02" value="3.210000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d10-x01-y03" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($1.00 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="3.668787e-02" errorPlus="3.668787e-02" value="4.830000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.373426e-02" errorPlus="3.373426e-02" value="4.780000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.354102e-02" errorPlus="3.354102e-02" value="4.580000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.275668e-02" errorPlus="3.275668e-02" value="4.560000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="3.162278e-02" errorPlus="3.162278e-02" value="3.990000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d11-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="4.720169e-02" errorPlus="4.720169e-02" value="9.780000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="3.417601e-02" errorPlus="3.417601e-02" value="9.520000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="3.201562e-02" errorPlus="3.201562e-02" value="9.310000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="3.289377e-02" errorPlus="3.289377e-02" value="8.890000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="4.134005e-02" errorPlus="4.134005e-02" value="8.100000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d12-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.5$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.250000e+00"/>
+ <measurement errorMinus="2.961419e-02" errorPlus="2.961419e-02" value="2.940000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="2.750000e+00"/>
+ <measurement errorMinus="2.816026e-02" errorPlus="2.816026e-02" value="2.790000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.250000e+00"/>
+ <measurement errorMinus="2.707397e-02" errorPlus="2.707397e-02" value="2.740000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="3.750000e+00"/>
+ <measurement errorMinus="2.617250e-02" errorPlus="2.617250e-02" value="2.760000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.250000e+00"/>
+ <measurement errorMinus="2.961419e-02" errorPlus="2.961419e-02" value="2.860000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d13-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)">
+ <dataPoint>
+ <measurement errorMinus="1.750000e-01" errorPlus="1.750000e-01" value="3.250000e-01"/>
+ <measurement errorMinus="3.676955e-02" errorPlus="3.676955e-02" value="9.540000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="5.750000e-01"/>
+ <measurement errorMinus="3.584690e-02" errorPlus="3.584690e-02" value="9.300000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="7.250000e-01"/>
+ <measurement errorMinus="3.584690e-02" errorPlus="3.584690e-02" value="9.430000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="9.000000e-01"/>
+ <measurement errorMinus="3.453983e-02" errorPlus="3.453983e-02" value="9.230000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="1.100000e+00"/>
+ <measurement errorMinus="3.534119e-02" errorPlus="3.534119e-02" value="9.360000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="6.500000e-01" errorPlus="6.500000e-01" value="1.850000e+00"/>
+ <measurement errorMinus="3.289377e-02" errorPlus="3.289377e-02" value="9.190000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d14-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($2.0 < y < 4.5$)">
+ <dataPoint>
+ <measurement errorMinus="1.750000e-01" errorPlus="1.750000e-01" value="3.250000e-01"/>
+ <measurement errorMinus="2.408319e-02" errorPlus="2.408319e-02" value="1.620000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="5.750000e-01"/>
+ <measurement errorMinus="2.517936e-02" errorPlus="2.517936e-02" value="2.310000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="7.500000e-02" errorPlus="7.500000e-02" value="7.250000e-01"/>
+ <measurement errorMinus="2.716616e-02" errorPlus="2.716616e-02" value="2.880000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="9.000000e-01"/>
+ <measurement errorMinus="2.828427e-02" errorPlus="2.828427e-02" value="3.510000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="1.000000e-01" errorPlus="1.000000e-01" value="1.100000e+00"/>
+ <measurement errorMinus="3.059412e-02" errorPlus="3.059412e-02" value="4.120000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="6.500000e-01" errorPlus="6.500000e-01" value="1.850000e+00"/>
+ <measurement errorMinus="3.436568e-02" errorPlus="3.436568e-02" value="4.920000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d15-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/\Lambda$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.000000e-01" errorPlus="2.000000e-01" value="4.720000e+00"/>
+ <measurement errorMinus="4.134005e-02" errorPlus="4.134005e-02" value="8.100000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="5.170000e+00"/>
+ <measurement errorMinus="3.289377e-02" errorPlus="3.289377e-02" value="8.890000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="5.670000e+00"/>
+ <measurement errorMinus="3.201562e-02" errorPlus="3.201562e-02" value="9.310000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="6.170000e+00"/>
+ <measurement errorMinus="3.417601e-02" errorPlus="3.417601e-02" value="9.520000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="6.670000e+00"/>
+ <measurement errorMinus="4.720169e-02" errorPlus="4.720169e-02" value="9.780000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+ <dataPointSet dimension="2" name="d16-x01-y01" path="/REF/LHCB_2011_I917009" title="$\bar{\Lambda}/K_{s}^{0}$ ratio at $\sqrt{s}=7$ TeV ($0.15 < p_\perp < 2.50$ GeV/c)">
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="6.670000e+00"/>
+ <measurement errorMinus="2.961419e-02" errorPlus="2.961419e-02" value="2.940000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="6.170000e+00"/>
+ <measurement errorMinus="2.816026e-02" errorPlus="2.816026e-02" value="2.790000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="5.670000e+00"/>
+ <measurement errorMinus="2.707397e-02" errorPlus="2.707397e-02" value="2.740000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="5.170000e+00"/>
+ <measurement errorMinus="2.617250e-02" errorPlus="2.617250e-02" value="2.760000e-01"/>
+ </dataPoint>
+ <dataPoint>
+ <measurement errorMinus="2.500000e-01" errorPlus="2.500000e-01" value="4.670000e+00"/>
+ <measurement errorMinus="2.961419e-02" errorPlus="2.961419e-02" value="2.860000e-01"/>
+ </dataPoint>
+ </dataPointSet>
+</aida>
Modified: branches/2011-07-aida2yoda/data/refdata/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/data/refdata/Makefile.am Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/data/refdata/Makefile.am Fri Mar 16 18:22:32 2012 (r3622)
@@ -1,5 +1,6 @@
dist_pkgdata_DATA = \
LHCB_2010_S8758301.aida \
+ LHCB_2011_I917009.aida \
LHCB_2011_I919315.aida \
JADE_OPAL_2000_S4300807.aida \
JADE_1998_S3612880.aida \
Modified: branches/2011-07-aida2yoda/include/Rivet/Analysis.hh
==============================================================================
--- branches/2011-07-aida2yoda/include/Rivet/Analysis.hh Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/include/Rivet/Analysis.hh Fri Mar 16 18:22:32 2012 (r3622)
@@ -343,6 +343,9 @@
/// Get the canonical histogram path for the named histogram in this analysis.
const std::string histoPath(const std::string& hname) const;
+ /// Get the canonical histogram path for the numbered histogram in this analysis.
+ const std::string histoPath(size_t datasetId, size_t xAxisId, size_t yAxisId) const;
+
/// Get the internal histogram name for given d, x and y (cf. HepData)
const std::string makeAxisCode(size_t datasetId, size_t xAxisId, size_t yAxisId) const;
Modified: branches/2011-07-aida2yoda/src/Analyses/ATLAS_2011_I925932.cc
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/ATLAS_2011_I925932.cc Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/src/Analyses/ATLAS_2011_I925932.cc Fri Mar 16 18:22:32 2012 (r3622)
@@ -9,6 +9,13 @@
namespace Rivet {
+ namespace {
+ inline double calcMT(const FourMomentum& a, const FourMomentum& b) {
+ return sqrt(2.0 * a.pT() * b.pT() * (1.0 - cos(a.phi() - b.phi())) );
+ }
+ }
+
+
class ATLAS_2011_I925932 : public Analysis {
public:
@@ -31,12 +38,6 @@
WFinder wfinder_bare_mu(fs, -2.4, 2.4, 20.0*GeV, MUON, 0.0*GeV, 1000.0*GeV, 25.0*GeV, 0.0,true,false);
addProjection(wfinder_bare_mu, "WFinder_bare_mu");
- // Counters
- _sumW_e_dressed = 0;
- _sumW_e_bare = 0;
- _sumW_mu_dressed = 0;
- _sumW_mu_bare = 0;
-
// Book histograms
_hist_wpt_dressed_el = bookHisto1D(1, 1, 1);
_hist_wpt_bare_el = bookHisto1D(1, 1, 2);
@@ -45,11 +46,6 @@
}
- inline double calcMT(const FourMomentum& a, const FourMomentum& b) {
- return sqrt(2.0 * a.pT() * b.pT() * (1.0 - cos(a.phi() - b.phi())) );
- }
-
-
/// Do the analysis
void analyze(const Event& event) {
@@ -69,7 +65,6 @@
const FourMomentum el = wfinder_dressed_el.constituentLeptons()[0].momentum();
const FourMomentum nu = wfinder_dressed_el.constituentNeutrinos()[0].momentum();
if (calcMT(el, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) {
- _sumW_e_dressed += event.weight();
_hist_wpt_dressed_el->fill(wfinder_dressed_el.bosons()[0].momentum().pT()/GeV, event.weight());
}
}
@@ -79,7 +74,6 @@
const FourMomentum el = wfinder_bare_el.constituentLeptons()[0].momentum();
const FourMomentum nu = wfinder_bare_el.constituentNeutrinos()[0].momentum();
if (calcMT(el, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) {
- _sumW_e_bare += event.weight();
_hist_wpt_bare_el->fill(wfinder_bare_el.bosons()[0].momentum().pT()/GeV, event.weight());
}
}
@@ -89,7 +83,6 @@
const FourMomentum mu = wfinder_dressed_mu.constituentLeptons()[0].momentum();
const FourMomentum nu = wfinder_dressed_mu.constituentNeutrinos()[0].momentum();
if (calcMT(mu, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) {
- _sumW_mu_dressed += event.weight();
_hist_wpt_dressed_mu->fill(wfinder_dressed_mu.bosons()[0].momentum().pT()/GeV, event.weight());
}
}
@@ -99,7 +92,6 @@
const FourMomentum mu = wfinder_bare_mu.constituentLeptons()[0].momentum();
const FourMomentum nu = wfinder_bare_mu.constituentNeutrinos()[0].momentum();
if (calcMT(mu, nu) > 40.0*GeV && nu.pT() > 25.0*GeV) {
- _sumW_mu_bare += event.weight();
_hist_wpt_bare_mu->fill(wfinder_bare_mu.bosons()[0].momentum().pT()/GeV, event.weight());
}
}
@@ -109,10 +101,10 @@
// Normalize histos
void finalize() {
- scale(_hist_wpt_dressed_el, 1/_sumW_e_dressed);
- scale(_hist_wpt_bare_el, 1/_sumW_e_bare);
- scale(_hist_wpt_dressed_mu, 1/_sumW_mu_dressed);
- scale(_hist_wpt_bare_mu, 1/_sumW_mu_bare);
+ normalize(_hist_wpt_dressed_el);
+ normalize(_hist_wpt_bare_el);
+ normalize(_hist_wpt_dressed_mu);
+ normalize(_hist_wpt_bare_mu);
}
//@}
@@ -120,8 +112,6 @@
private:
- double _sumW_e_dressed, _sumW_e_bare, _sumW_mu_dressed, _sumW_mu_bare;
-
Histo1DPtr _hist_wpt_dressed_el;
Histo1DPtr _hist_wpt_bare_el;
Histo1DPtr _hist_wpt_dressed_mu;
Copied and modified: branches/2011-07-aida2yoda/src/Analyses/LHCB_2011_I917009.cc (from r3583, trunk/src/Analyses/LHCB_2011_I917009.cc)
==============================================================================
--- trunk/src/Analyses/LHCB_2011_I917009.cc Sun Feb 26 22:39:15 2012 (r3583, copy source)
+++ branches/2011-07-aida2yoda/src/Analyses/LHCB_2011_I917009.cc Fri Mar 16 18:22:32 2012 (r3622)
@@ -1,9 +1,9 @@
// -*- C++ -*-
#include "Rivet/Analysis.hh"
-#include "Rivet/RivetAIDA.hh"
+#include "Rivet/RivetYODA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/UnstableFinalState.hh"
-#include "LWH/Histogram1D.h"
+//#include "LWH/Histogram1D.h"
#include "Rivet/Math/MathUtils.hh"
#include "Rivet/Math/Constants.hh"
@@ -59,12 +59,13 @@
} else {
MSG_ERROR("Incompatible beam energy!");
}
- //pT edges are the same for all 3 histos in this suite
- const BinEdges& pt_edges = binEdges(dsShift+5,1,1);
- // booking temporary histos
- for (int i = 0; i<12; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(y_nbins, rap_min, rap_max));
- for (int i = 12; i<15; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(pt_edges));
- for (int i = 15; i<18; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(y_nbins, rap_beam - rap_max, rap_beam - rap_min));
+ /// @todo YODA
+ ////pT edges are the same for all 3 histos in this suite
+ //const BinEdges& pt_edges = binEdges(dsShift+5,1,1);
+ //// booking temporary histos
+ //for (int i = 0; i<12; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(y_nbins, rap_min, rap_max));
+ //for (int i = 12; i<15; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(pt_edges));
+ //for (int i = 15; i<18; i++ ) _tmphistos[i].reset(new LWH::Histogram1D(y_nbins, rap_beam - rap_max, rap_beam - rap_min));
addProjection(UnstableFinalState(), "UFS");
}
@@ -101,36 +102,38 @@
if ( (y < rap_min) || (y > rap_max) ) continue;
pT = sqrt((qmom.px() * qmom.px()) + (qmom.py() * qmom.py()));
if ( (pT < pt_min) || (pT > pt3_edge) ) continue;
- // Filling corresponding temporary histograms for pT intervals
- if ( (pT >= pt_min ) && (pT < pt1_edge) ) _tmphistos[partIdx*3]->fill(y, weight);
- if ( (pT >= pt1_edge) && (pT < pt2_edge) ) _tmphistos[partIdx*3+1]->fill(y, weight);
- if ( (pT >= pt2_edge) && (pT < pt3_edge) ) _tmphistos[partIdx*3+2]->fill(y, weight);
- // Fill histo in rapidity for whole pT interval
- _tmphistos[partIdx+9]->fill(y, weight);
- // Fill histo in pT for whole rapidity interval
- _tmphistos[partIdx+12]->fill(pT, weight);
- // Fill histo in rapidity loss for whole pT interval
- _tmphistos[partIdx+15]->fill(rap_beam - y, weight);
+ /// @todo YODA
+ //// Filling corresponding temporary histograms for pT intervals
+ //if ( (pT >= pt_min ) && (pT < pt1_edge) ) _tmphistos[partIdx*3]->fill(y, weight);
+ //if ( (pT >= pt1_edge) && (pT < pt2_edge) ) _tmphistos[partIdx*3+1]->fill(y, weight);
+ //if ( (pT >= pt2_edge) && (pT < pt3_edge) ) _tmphistos[partIdx*3+2]->fill(y, weight);
+ //// Fill histo in rapidity for whole pT interval
+ //_tmphistos[partIdx+9]->fill(y, weight);
+ //// Fill histo in pT for whole rapidity interval
+ //_tmphistos[partIdx+12]->fill(pT, weight);
+ //// Fill histo in rapidity loss for whole pT interval
+ //_tmphistos[partIdx+15]->fill(rap_beam - y, weight);
}
}
// Generate the ratio histograms
void finalize() {
- // needed to determine AIDA to save the file!
- tree().mkdirs(histoDir());
- int dsId = dsShift+1;
- for (int j=0; j<3; j++ ) {
- histogramFactory().divide(histoPath(dsId, 1, j+1), *_tmphistos[j], *_tmphistos[3+j]);
- histogramFactory().divide(histoPath(dsId+1, 1, j+1), *_tmphistos[j], *_tmphistos[6+j]);
- }
- dsId += 2;
- for (int j = 3; j<6; j++) {
- histogramFactory().divide(histoPath(dsId, 1, 1), *_tmphistos[3*j], *_tmphistos[3*j+1]);
- dsId ++;
- histogramFactory().divide(histoPath(dsId, 1, 1), *_tmphistos[3*j], *_tmphistos[3*j+2]);
- dsId ++;
- }
+ /// @todo YODA
+ //// needed to determine AIDA to save the file!
+ //tree().mkdirs(histoDir());
+ //int dsId = dsShift+1;
+ //for (int j=0; j<3; j++ ) {
+ // histogramFactory().divide(histoPath(dsId, 1, j+1), *_tmphistos[j], *_tmphistos[3+j]);
+ // histogramFactory().divide(histoPath(dsId+1, 1, j+1), *_tmphistos[j], *_tmphistos[6+j]);
+ //}
+ //dsId += 2;
+ //for (int j = 3; j<6; j++) {
+ // histogramFactory().divide(histoPath(dsId, 1, 1), *_tmphistos[3*j], *_tmphistos[3*j+1]);
+ // dsId ++;
+ // histogramFactory().divide(histoPath(dsId, 1, 1), *_tmphistos[3*j], *_tmphistos[3*j+2]);
+ // dsId ++;
+ //}
}
//@}
@@ -210,7 +213,7 @@
// The limits of the rapidity window
double rap_min;
double rap_max;
- // Indicates which set of histograms will be output to aida file (according to beam energy)
+ // Indicates which set of histograms will be output to yoda file (according to beam energy)
int dsShift;
// Map between PDG id and particle lifetimes in seconds
std::map<int, double> partLftMap;
@@ -218,7 +221,8 @@
static const int stablePDGIds[205];
/// @name Helper Histograms:
//@{
- shared_ptr<LWH::Histogram1D> _tmphistos[18];
+ /// @todo YODA
+ //shared_ptr<LWH::Histogram1D> _tmphistos[18];
// Histograms are defined in the following order: anti-Lambda, Lambda and K0s.
// First 3 suites of 3 histograms correspond to each particle in bins of y for the 3 pT intervals. (9 histos)
// Next 3 histograms contain the particles in y bins for the whole pT interval (3 histos)
@@ -336,7 +340,7 @@
9900024, 9900041, 9900042 };
- // Plugin hook
+ // Hook for the plugin system
DECLARE_RIVET_PLUGIN(LHCB_2011_I917009);
}
Modified: branches/2011-07-aida2yoda/src/Analyses/Makefile.am
==============================================================================
--- branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/src/Analyses/Makefile.am Fri Mar 16 18:22:32 2012 (r3622)
@@ -98,6 +98,7 @@
if ENABLE_UNVALIDATED
RivetLHCbAnalyses_la_SOURCES += \
LHCB_2010_S8758301.cc \
+ LHCB_2011_I917009.cc \
LHCB_2011_I919315.cc
endif
Modified: branches/2011-07-aida2yoda/src/Core/Analysis.cc
==============================================================================
--- branches/2011-07-aida2yoda/src/Core/Analysis.cc Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/src/Core/Analysis.cc Fri Mar 16 18:22:32 2012 (r3622)
@@ -53,6 +53,11 @@
}
+ const string Analysis::histoPath(size_t datasetId, size_t xAxisId, size_t yAxisId) const {
+ return histoDir() + "/" + makeAxisCode(datasetId, xAxisId, yAxisId);
+ }
+
+
const string Analysis::makeAxisCode(size_t datasetId, size_t xAxisId, size_t yAxisId) const {
stringstream axisCode;
axisCode << "d";
Modified: branches/2011-07-aida2yoda/src/Projections/InitialQuarks.cc
==============================================================================
--- branches/2011-07-aida2yoda/src/Projections/InitialQuarks.cc Fri Mar 16 17:51:40 2012 (r3621)
+++ branches/2011-07-aida2yoda/src/Projections/InitialQuarks.cc Fri Mar 16 18:22:32 2012 (r3622)
@@ -22,7 +22,7 @@
const GenVertex* pv = p->production_vertex();
const GenVertex* dv = p->end_vertex();
const PdgId pid = abs(p->pdg_id());
- bool passed = (inRange(pid, 1, 5));
+ bool passed = inRange((long)pid, 1, 5);
if (passed) {
if (pv != 0) {
foreach (const GenParticle* pp, particles_in(pv)) {
More information about the Rivet-svn
mailing list
|