[Rivet-svn] r4084 - in branches/2012-06-aidarivet: . data/anainfo data/plotinfo data/refdata src/Analyses

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Fri Dec 14 13:57:51 GMT 2012


Author: hoeth
Date: Fri Dec 14 13:57:50 2012
New Revision: 4084

Log:
Added LHCB_2012_I1119400 analysis

Added:
   branches/2012-06-aidarivet/data/anainfo/LHCB_2012_I1119400.info
   branches/2012-06-aidarivet/data/plotinfo/LHCB_2012_I1119400.plot
   branches/2012-06-aidarivet/data/refdata/LHCB_2012_I1119400.aida
   branches/2012-06-aidarivet/src/Analyses/LHCB_2012_I1119400.cc
Modified:
   branches/2012-06-aidarivet/ChangeLog
   branches/2012-06-aidarivet/data/anainfo/Makefile.am
   branches/2012-06-aidarivet/data/plotinfo/Makefile.am
   branches/2012-06-aidarivet/data/refdata/Makefile.am
   branches/2012-06-aidarivet/src/Analyses/Makefile.am

Modified: branches/2012-06-aidarivet/ChangeLog
==============================================================================
--- branches/2012-06-aidarivet/ChangeLog	Fri Dec 14 13:38:42 2012	(r4083)
+++ branches/2012-06-aidarivet/ChangeLog	Fri Dec 14 13:57:50 2012	(r4084)
@@ -2,6 +2,8 @@
 
 	* Added CMS_2012_PAS_FWD_11_003 analysis
 
+	* Added LHCB_2012_I1119400 analysis
+
 2012-12-12  Andy Buckley  <andy.buckley at cern.ch>
 
 	* Correction to jet acceptance in CMS_2011_S9120041, from Sercan Sen: thanks!

Added: branches/2012-06-aidarivet/data/anainfo/LHCB_2012_I1119400.info
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ branches/2012-06-aidarivet/data/anainfo/LHCB_2012_I1119400.info	Fri Dec 14 13:57:50 2012	(r4084)
@@ -0,0 +1,41 @@
+Name: LHCB_2012_I1119400
+Year: 2012
+Summary: Measurement of prompt hadron production ratios in $pp$ collisions at $\sqrt{s} = $ 0.9 and 7 TeV
+Experiment: LHCb
+Collider: LHC
+SpireID: 9635220
+InspireID: 1119400
+Status: UNVALIDATED
+Authors:
+ - Andrea Contu <Andrea.Contu at cern.ch>
+ - Alex Grecu <Alex.Grecu at cern.ch>
+References:
+ - arXiv:1206.5160 [hep-ex]
+RunInfo:
+   Minimum bias events at sqrt(s) = 0.9 and 7 TeV.
+NumEvents: 500000
+NeedCrossSection: no
+Beams: [p+, p+]
+Energies: [900, 7000]
+PtCuts: [0]
+Description:
+  'Measurement of the production ratios of prompt charged particles (protons,
+   pions and kaons). Promptness is defined as originating from the primary
+   interaction, either directly, or through the subsequent decay of a
+   resonance.'
+BibKey: Aaij:2012ut
+BibTeX: '@article{Aaij:2012ut,
+      author         = "Aaij, R and others",
+      title          = "{Measurement of prompt hadron production ratios in $pp$
+                        collisions at $\sqrt{s} = $ 0.9 and 7 TeV}",
+      collaboration  = "LHCb Collaboration",
+      journal        = "Eur.Phys.J.",
+      volume         = "C72",
+      pages          = "2168",
+      year           = "2012",
+      eprint         = "1206.5160",
+      archivePrefix  = "arXiv",
+      primaryClass   = "hep-ex",
+      reportNumber   = "CERN-PH-EP-2012-171, LHCB-PAPER-2011-037",
+      SLACcitation   = "%%CITATION = ARXIV:1206.5160;%%"
+}'

Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/Makefile.am	Fri Dec 14 13:38:42 2012	(r4083)
+++ branches/2012-06-aidarivet/data/anainfo/Makefile.am	Fri Dec 14 13:57:50 2012	(r4084)
@@ -172,6 +172,7 @@
   LHCB_2010_I867355.info \
   LHCB_2011_I917009.info \
   LHCB_2011_I919315.info \
+  LHCB_2012_I1119400.info \
   LHCF_2012_I1115479.info \
   MC_DIJET.info \
   MC_DIPHOTON.info \

Added: branches/2012-06-aidarivet/data/plotinfo/LHCB_2012_I1119400.plot
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ branches/2012-06-aidarivet/data/plotinfo/LHCB_2012_I1119400.plot	Fri Dec 14 13:57:50 2012	(r4084)
@@ -0,0 +1,180 @@
+# BEGIN PLOT /LHCB_2012_I1119400/d
+XLabel=$\eta$
+LogY=0
+LegendXPos=0.05
+LegendYPos=0.4
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d0[1-2]-
+YLabel=$\bar{p}/p$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d01-x01-y01
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d01-x01-y02
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d01-x01-y03
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d02-x01-y01
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d02-x01-y02
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d02-x01-y03
+Title=$\bar{p}/p$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d0[3-4]-
+YLabel=$K^-/K^+$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d03-x01-y01
+Title=$K^-/K^+$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d03-x01-y02
+Title=$K^-/K^+$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d03-x01-y03
+Title=$K^-/K^+$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d04-x01-y01
+Title=$K^-/K^+$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d04-x01-y02
+Title=$K^-/K^+$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d04-x01-y03
+Title=$K^-/K^+$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d0[5-6]-
+YLabel=$\pi^-/\pi^+$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d05-x01-y01
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d05-x01-y02
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d05-x01-y03
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d06-x01-y01
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d06-x01-y02
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d06-x01-y03
+Title=$\pi^-/\pi^+$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d0[7-8]-
+YLabel=$(p+\bar{p})/(\pi^+ + \pi^-)$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d07-x01-y01
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d07-x01-y02
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d07-x01-y03
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d08-x01-y01
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d08-x01-y02
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d08-x01-y03
+Title=$(p+\bar{p})/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d(09|10)-
+YLabel=$(K^+ + K^-)/(\pi^+ + \pi^-)$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d09-x01-y01
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d09-x01-y02
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d09-x01-y03
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d10-x01-y01
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d10-x01-y02
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d10-x01-y03
+Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+
+# BEGIN PLOT /LHCB_2012_I1119400/d1[1-2]-
+YLabel=$(p+\bar{p})/(K^+ + K^-)$
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d11-x01-y01
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d11-x01-y02
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=0.9$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d11-x01-y03
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=0.9$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d12-x01-y01
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=7$ TeV ($0.0 < p_\perp < 0.8$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d12-x01-y02
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=7$ TeV ($0.8 < p_\perp < 1.2$ GeV/$c$)
+# END PLOT
+
+# BEGIN PLOT /LHCB_2012_I1119400/d12-x01-y03
+Title=$(p+\bar{p})/(K^+ + K^-)$ ratio at $\sqrt{s}=7$ TeV ($p_\perp > 1.2$ GeV/$c$)
+# END PLOT

Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Fri Dec 14 13:38:42 2012	(r4083)
+++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Fri Dec 14 13:57:50 2012	(r4084)
@@ -169,6 +169,7 @@
   LHCB_2010_I867355.plot \
   LHCB_2011_I917009.plot \
   LHCB_2011_I919315.plot \
+  LHCB_2012_I1119400.plot \
   LHCF_2012_I1115479.plot \
   MC_GENERIC.plot \
   MC_IDENTIFIED.plot \

Added: branches/2012-06-aidarivet/data/refdata/LHCB_2012_I1119400.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ branches/2012-06-aidarivet/data/refdata/LHCB_2012_I1119400.aida	Fri Dec 14 13:57:50 2012	(r4084)
@@ -0,0 +1,605 @@
+<?xml version="1.0" encoding="utf-8" ?>
+<!DOCTYPE aida SYSTEM "http://aida.freehep.org/schemas/3.3/aida.dtd">
+<aida version="3.3">
+  <implementation version="1.0" package="HepData"/>
+  <dataPointSet name="d01-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.751" errorPlus="0.16337074401495513" errorMinus="0.16337074401495513"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.729" errorPlus="0.040607881008493905" errorMinus="0.040607881008493905"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.66" errorPlus="0.04687216658103186" errorMinus="0.04687216658103186"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.107" errorPlus="0.3495725961799637" errorMinus="0.3495725961799637"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.684" errorPlus="0.05000999900019995" errorMinus="0.05000999900019995"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.576" errorPlus="0.0327566787083184" errorMinus="0.0327566787083184"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.451" errorPlus="0.03905124837953327" errorMinus="0.03905124837953327"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d01-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.794" errorPlus="0.0902551937563706" errorMinus="0.0902551937563706"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.614" errorPlus="0.04805205510693585" errorMinus="0.04805205510693585"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.456" errorPlus="0.03420526275297414" errorMinus="0.03420526275297414"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.328" errorPlus="0.05000999900019995" errorMinus="0.05000999900019995"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.734" errorPlus="0.12448694710691559" errorMinus="0.12448694710691559"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.015" errorPlus="0.038078865529319536" errorMinus="0.038078865529319536"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.957" errorPlus="0.05197114584074512" errorMinus="0.05197114584074512"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.181" errorPlus="0.1960229578391266" errorMinus="0.1960229578391266"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.942" errorPlus="0.037643060449437424" errorMinus="0.037643060449437424"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.916" errorPlus="0.02308679276123039" errorMinus="0.02308679276123039"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.906" errorPlus="0.04026164427839479" errorMinus="0.04026164427839479"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d02-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PBAR)/SIG(Q=P)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.88" errorPlus="0.040024992192379" errorMinus="0.040024992192379"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.905" errorPlus="0.027202941017470887" errorMinus="0.027202941017470887"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.903" errorPlus="0.02435159132377184" errorMinus="0.02435159132377184"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.831" errorPlus="0.050990195135927854" errorMinus="0.050990195135927854"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.834" errorPlus="0.06935416353759881" errorMinus="0.06935416353759881"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.001" errorPlus="0.0643816744112795" errorMinus="0.0643816744112795"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.002" errorPlus="0.09326306878931231" errorMinus="0.09326306878931231"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.87" errorPlus="0.2671872002922296" errorMinus="0.2671872002922296"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.847" errorPlus="0.041" errorMinus="0.041"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.792" errorPlus="0.029410882339705485" errorMinus="0.029410882339705485"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.68" errorPlus="0.04220189569201839" errorMinus="0.04220189569201839"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d03-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.936" errorPlus="0.07021395872616784" errorMinus="0.07021395872616784"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.783" errorPlus="0.03860051813123756" errorMinus="0.03860051813123756"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.723" errorPlus="0.03324154027718932" errorMinus="0.03324154027718932"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.506" errorPlus="0.05192301994298868" errorMinus="0.05192301994298868"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d04-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.992" errorPlus="0.06328506932918695" errorMinus="0.06328506932918695"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.062" errorPlus="0.04031128874149275" errorMinus="0.04031128874149275"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.161" errorPlus="0.055226805085936304" errorMinus="0.055226805085936304"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d04-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.995" errorPlus="0.10131633629380803" errorMinus="0.10131633629380803"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.966" errorPlus="0.019924858845171277" errorMinus="0.019924858845171277"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.948" errorPlus="0.015231546211727818" errorMinus="0.015231546211727818"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.898" errorPlus="0.025709920264364885" errorMinus="0.025709920264364885"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d04-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=K-)/SIG(Q=K+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.991" errorPlus="0.02213594362117866" errorMinus="0.02213594362117866"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.999" errorPlus="0.017088007490635063" errorMinus="0.017088007490635063"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.93" errorPlus="0.01802775637731995" errorMinus="0.01802775637731995"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.958" errorPlus="0.049819674828324605" errorMinus="0.049819674828324605"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d05-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.979" errorPlus="0.011180339887498949" errorMinus="0.011180339887498949"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.968" errorPlus="0.011704699910719625" errorMinus="0.011704699910719625"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.928" errorPlus="0.01746424919657298" errorMinus="0.01746424919657298"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d05-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.987" errorPlus="0.016401219466856725" errorMinus="0.016401219466856725"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.971" errorPlus="0.014866068747318505" errorMinus="0.014866068747318505"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.951" errorPlus="0.015620499351813309" errorMinus="0.015620499351813309"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.971" errorPlus="0.02483948469674844" errorMinus="0.02483948469674844"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d05-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.97" errorPlus="0.0212602916254693" errorMinus="0.0212602916254693"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.926" errorPlus="0.022022715545545243" errorMinus="0.022022715545545243"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.871" errorPlus="0.022472205054244233" errorMinus="0.022472205054244233"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.816" errorPlus="0.038288379438153296" errorMinus="0.038288379438153296"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d06-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.011" errorPlus="0.007211102550927979" errorMinus="0.007211102550927979"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.002" errorPlus="0.00670820393249937" errorMinus="0.00670820393249937"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.976" errorPlus="0.00670820393249937" errorMinus="0.00670820393249937"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d06-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.002" errorPlus="0.009219544457292887" errorMinus="0.009219544457292887"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.998" errorPlus="0.00806225774829855" errorMinus="0.00806225774829855"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.003" errorPlus="0.008544003745317531" errorMinus="0.008544003745317531"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.998" errorPlus="0.012041594578792296" errorMinus="0.012041594578792296"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d06-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PI-)/SIG(Q=PI+)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.015" errorPlus="0.011180339887498949" errorMinus="0.011180339887498949"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.998" errorPlus="0.010770329614269008" errorMinus="0.010770329614269008"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.0" errorPlus="0.012083045973594572" errorMinus="0.012083045973594572"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.974" errorPlus="0.020808652046684813" errorMinus="0.020808652046684813"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d07-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.086" errorPlus="0.02102379604162864" errorMinus="0.02102379604162864"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.062" errorPlus="0.00806225774829855" errorMinus="0.00806225774829855"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.076" errorPlus="0.01004987562112089" errorMinus="0.01004987562112089"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d07-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.328" errorPlus="0.10423531071570709" errorMinus="0.10423531071570709"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.208" errorPlus="0.01649242250247064" errorMinus="0.01649242250247064"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.175" errorPlus="0.011401754250991379" errorMinus="0.011401754250991379"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.233" errorPlus="0.022803508501982758" errorMinus="0.022803508501982758"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d07-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.3" errorPlus="0.03492849839314596" errorMinus="0.03492849839314596"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.272" errorPlus="0.024041630560342614" errorMinus="0.024041630560342614"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.252" errorPlus="0.02118962010041709" errorMinus="0.02118962010041709"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.301" errorPlus="0.04876474136094643" errorMinus="0.04876474136094643"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d08-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.085" errorPlus="0.017029386365926404" errorMinus="0.017029386365926404"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.069" errorPlus="0.00806225774829855" errorMinus="0.00806225774829855"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.051" errorPlus="0.007071067811865475" errorMinus="0.007071067811865475"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d08-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.235" errorPlus="0.03920459156782532" errorMinus="0.03920459156782532"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.174" errorPlus="0.009219544457292887" errorMinus="0.009219544457292887"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.156" errorPlus="0.006324555320336759" errorMinus="0.006324555320336759"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.184" errorPlus="0.01044030650891055" errorMinus="0.01044030650891055"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d08-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.262" errorPlus="0.014560219778561038" errorMinus="0.014560219778561038"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.245" errorPlus="0.011401754250991379" errorMinus="0.011401754250991379"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.242" errorPlus="0.01044030650891055" errorMinus="0.01044030650891055"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.244" errorPlus="0.01746424919657298" errorMinus="0.01746424919657298"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d09-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.18" errorPlus="0.026076809620810593" errorMinus="0.026076809620810593"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.171" errorPlus="0.023021728866442676" errorMinus="0.023021728866442676"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.173" errorPlus="0.02501999200639361" errorMinus="0.02501999200639361"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d09-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.184" errorPlus="0.05608029957123981" errorMinus="0.05608029957123981"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.267" errorPlus="0.015524174696260024" errorMinus="0.015524174696260024"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.247" errorPlus="0.011704699910719625" errorMinus="0.011704699910719625"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.268" errorPlus="0.018973665961010275" errorMinus="0.018973665961010275"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d09-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.351" errorPlus="0.029120439557122076" errorMinus="0.029120439557122076"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.319" errorPlus="0.01969771560359221" errorMinus="0.01969771560359221"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.314" errorPlus="0.019235384061671346" errorMinus="0.019235384061671346"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.281" errorPlus="0.03324154027718932" errorMinus="0.03324154027718932"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d10-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.181" errorPlus="0.024020824298928628" errorMinus="0.024020824298928628"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.173" errorPlus="0.02102379604162864" errorMinus="0.02102379604162864"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.131" errorPlus="0.016031219541881397" errorMinus="0.016031219541881397"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d10-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.224" errorPlus="0.024083189157584593" errorMinus="0.024083189157584593"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.263" errorPlus="0.01044030650891055" errorMinus="0.01044030650891055"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.262" errorPlus="0.009486832980505138" errorMinus="0.009486832980505138"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.275" errorPlus="0.011401754250991379" errorMinus="0.011401754250991379"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d10-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=KAON)/SIG(Q=PION)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.371" errorPlus="0.014560219778561038" errorMinus="0.014560219778561038"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.357" errorPlus="0.012649110640673518" errorMinus="0.012649110640673518"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.367" errorPlus="0.01392838827718412" errorMinus="0.01392838827718412"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.328" errorPlus="0.020615528128088305" errorMinus="0.020615528128088305"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d11-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.481" errorPlus="0.13923002549737612" errorMinus="0.13923002549737612"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.363" errorPlus="0.06612110101926616" errorMinus="0.06612110101926616"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.433" errorPlus="0.08628441342444185" errorMinus="0.08628441342444185"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d11-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.831" errorPlus="0.8229246624084127" errorMinus="0.8229246624084127"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.779" errorPlus="0.07433034373659253" errorMinus="0.07433034373659253"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.709" errorPlus="0.05629387178015028" errorMinus="0.05629387178015028"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.865" errorPlus="0.10022474744293448" errorMinus="0.10022474744293448"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d11-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.855" errorPlus="0.12066896867049125" errorMinus="0.12066896867049125"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.851" errorPlus="0.08612200647918046" errorMinus="0.08612200647918046"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.799" errorPlus="0.07884795495128584" errorMinus="0.07884795495128584"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.067" errorPlus="0.20500000000000002" errorMinus="0.20500000000000002"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d12-x01-y01" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.465" errorPlus="0.11128791488746656" errorMinus="0.11128791488746656"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.398" errorPlus="0.06711929677819935" errorMinus="0.06711929677819935"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.379" errorPlus="0.06811754546370562" errorMinus="0.06811754546370562"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d12-x01-y02" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="1.051" errorPlus="0.2049780476051033" errorMinus="0.2049780476051033"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.66" errorPlus="0.040024992192379" errorMinus="0.040024992192379"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.593" errorPlus="0.032557641192199414" errorMinus="0.032557641192199414"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.671" errorPlus="0.0478539444560216" errorMinus="0.0478539444560216"/>
+    </dataPoint>
+  </dataPointSet>
+  <dataPointSet name="d12-x01-y03" dimension="2" path="/REF/LHCB_2012_I1119400" title="SIG(Q=PROTON)/SIG(Q=KAON)" >
+    <dataPoint>
+      <measurement value="2.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.705" errorPlus="0.0478539444560216" errorMinus="0.0478539444560216"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.682" errorPlus="0.03863935817272331" errorMinus="0.03863935817272331"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="3.75" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.659" errorPlus="0.037656340767525456" errorMinus="0.037656340767525456"/>
+    </dataPoint>
+    <dataPoint>
+      <measurement value="4.25" errorPlus="0.25" errorMinus="0.25"/>
+      <measurement value="0.744" errorPlus="0.06987131027825369" errorMinus="0.06987131027825369"/>
+    </dataPoint>
+  </dataPointSet>
+</aida>

Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/refdata/Makefile.am	Fri Dec 14 13:38:42 2012	(r4083)
+++ branches/2012-06-aidarivet/data/refdata/Makefile.am	Fri Dec 14 13:57:50 2012	(r4084)
@@ -3,6 +3,7 @@
   LHCB_2010_I867355.aida \
   LHCB_2011_I917009.aida \
   LHCB_2011_I919315.aida \
+  LHCB_2012_I1119400.aida \
   LHCF_2012_I1115479.aida \
   JADE_OPAL_2000_S4300807.aida \
   JADE_1998_S3612880.aida \

Added: branches/2012-06-aidarivet/src/Analyses/LHCB_2012_I1119400.cc
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ branches/2012-06-aidarivet/src/Analyses/LHCB_2012_I1119400.cc	Fri Dec 14 13:57:50 2012	(r4084)
@@ -0,0 +1,376 @@
+// -*- C++ -*-
+#include "Rivet/Analysis.hh"
+#include "Rivet/RivetAIDA.hh"
+#include "Rivet/Tools/Logging.hh"
+#include "Rivet/Tools/ParticleIdUtils.hh"
+#include "Rivet/Projections/ChargedFinalState.hh"
+#include "LWH/Histogram1D.h"
+#include "Rivet/Math/MathUtils.hh"
+#include "Rivet/Math/Constants.hh"
+
+#include "HepMC/GenEvent.h"
+#include "HepMC/GenParticle.h"
+#include "HepMC/GenVertex.h"
+#include "HepMC/SimpleVector.h"
+
+#include <iostream>
+#include <sstream>
+#include <string>
+#include <vector>
+#include <map>
+#include <utility>
+
+namespace Rivet {
+
+
+  class LHCB_2012_I1119400 : public Analysis {
+    public:
+
+    /// @name Constructors etc.
+    //@{
+
+    /// Constructor
+    LHCB_2012_I1119400() : Analysis("LHCB_2012_I1119400"),
+      _p_min(5.0),
+      _pt_min(0.0),_pt1_edge(0.8), _pt2_edge(1.2),
+      _eta_nbins(4),
+      _eta_min(2.5),
+      _eta_max(4.5)
+    {   }
+
+    //@}
+
+
+    public:
+
+    /// @name Analysis methods
+    //@{
+
+    /// Book histograms and initialise projections before the run
+    void init() {
+      fillMap(_partLftMap);
+
+      int id_shift = 0;
+      if (fuzzyEquals(sqrtS(), 7*TeV)) id_shift = 1;
+      // define ratios if second pdgid in pair is -1, it means that is a antiparticle/particle ratio
+      _ratiotypes["pbarp"] = make_pair(2212, -1);
+      _ratiotypes_hepdataid["pbarp"] = 1 + id_shift;
+      _ratiotypes["kminuskplus"] = make_pair(321, -1);
+      _ratiotypes_hepdataid["kminuskplus"] = 3 + id_shift;
+      _ratiotypes["piminuspiplus"] = make_pair(211, -1);
+      _ratiotypes_hepdataid["piminuspiplus"] = 5 + id_shift;
+      _ratiotypes["pk"] = make_pair(2212, 321);
+      _ratiotypes_hepdataid["pk"] = 11 + id_shift;
+      _ratiotypes["ppi"] = make_pair(2212, 211);
+      _ratiotypes_hepdataid["ppi"] = 7 + id_shift;
+      _ratiotypes["kpi"] = make_pair(321, 211);
+      _ratiotypes_hepdataid["kpi"] = 9 + id_shift;
+
+      std::map<std::string, std::pair<int, int> >::iterator it;
+
+      // booking histograms
+      for (it=_ratiotypes.begin(); it!=_ratiotypes.end(); it++) {
+        _h_ratio_lowpt[it->first]  = new LWH::Histogram1D(_eta_nbins-1, _eta_min+0.5, _eta_max);
+        _h_ratio_midpt[it->first]  = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+        _h_ratio_highpt[it->first] = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+        _h_num_lowpt[it->first]    = new LWH::Histogram1D(_eta_nbins-1, _eta_min+0.5, _eta_max);
+        _h_num_midpt[it->first]    = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+        _h_num_highpt[it->first]   = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+        _h_den_lowpt[it->first]    = new LWH::Histogram1D(_eta_nbins-1, _eta_min+0.5, _eta_max);
+        _h_den_midpt[it->first]    = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+        _h_den_highpt[it->first]   = new LWH::Histogram1D(_eta_nbins,   _eta_min,     _eta_max);
+      }
+
+      addProjection(ChargedFinalState(_eta_min, _eta_max, _pt_min*GeV), "CFS");
+    }
+
+
+    // Perform the per-event analysis
+    void analyze(const Event& event) {
+      const double weight = event.weight();
+      const ChargedFinalState& cfs = applyProjection<ChargedFinalState>(event, "CFS");
+
+      foreach (const Particle& p, cfs.particles()) {
+        int id = p.pdgId();
+        // continue if particle not a proton, a kaon or a pion
+        if ( !( (abs(id) == 211) || (abs(id) == 321) || (abs(id) == 2212)))  {
+          continue;
+        }
+
+        // cut in momentum
+        const FourMomentum& qmom = p.momentum();
+        if (qmom.p().mod() < _p_min) continue;
+
+        // Lifetime cut: ctau sum of all particle ancestors < 10^-9 m according to the paper (see eq. 5)
+        const double MAX_CTAU = 1.0e-9; // [m]
+        double ancestor_lftsum = getMotherLifeTimeSum(p);
+        if ( (ancestor_lftsum < 0.0) || (ancestor_lftsum > MAX_CTAU) ) continue;
+
+        double eta = qmom.eta();
+        double pT  = qmom.pT();
+
+        std::map<std::string, std::pair<int, int> >::iterator it;
+
+        for (it=_ratiotypes.begin(); it!=_ratiotypes.end(); it++) {
+          // check what type of ratio is
+          if ((it->second.second)==-1) {
+            // check ptbin
+            if (pT < _pt1_edge) {
+              // filling histos for numerator and denominator
+              if (id == -abs(it->second.first)) _h_num_lowpt[it->first]->fill(eta, weight);
+              if (id ==  abs(it->second.first)) _h_den_lowpt[it->first]->fill(eta, weight);
+            }
+            else if (pT < _pt2_edge) {
+              // filling histos for numerator and denominator
+              if (id == -abs(it->second.first)) _h_num_midpt[it->first]->fill(eta, weight);
+              if (id ==  abs(it->second.first)) _h_den_midpt[it->first]->fill(eta, weight);
+            }
+            else {
+              // filling histos for numerator and denominator
+              if (id == -abs(it->second.first)) _h_num_highpt[it->first]->fill(eta, weight);
+              if (id ==  abs(it->second.first)) _h_den_highpt[it->first]->fill(eta, weight);
+            }
+          }
+          else {
+            // check what type of ratio is
+            if (pT < _pt1_edge) {
+              // filling histos for numerator and denominator
+              if (abs(id) == abs(it->second.first))  _h_num_lowpt[it->first]->fill(eta, weight);
+              if (abs(id) == abs(it->second.second)) _h_den_lowpt[it->first]->fill(eta, weight);
+            }
+            else if (pT < _pt2_edge) {
+              // filling histos for numerator and denominator
+              if (abs(id) == abs(it->second.first))  _h_num_midpt[it->first]->fill(eta, weight);
+              if (abs(id) == abs(it->second.second)) _h_den_midpt[it->first]->fill(eta, weight);
+            }
+            else {
+              // filling histos for numerator and denominator
+              if (abs(id) == abs(it->second.first))  _h_num_highpt[it->first]->fill(eta, weight);
+              if (abs(id) == abs(it->second.second)) _h_den_highpt[it->first]->fill(eta, weight);
+            }
+          }
+        }
+      }
+    }
+
+
+    // Generate the ratio histograms
+    void finalize() {
+    // needed to determine AIDA to save the file!
+      tree().mkdirs(histoDir());
+
+      std::map<std::string, std::pair<int, int> >::iterator it;
+
+      // booking histograms
+      for (it=_ratiotypes.begin(); it!=_ratiotypes.end(); it++) {
+        int dsId=_ratiotypes_hepdataid[it->first];
+        histogramFactory().divide(histoPath(dsId, 1, 1), *(_h_num_lowpt[it->first]),  *(_h_den_lowpt[it->first]));
+        histogramFactory().divide(histoPath(dsId, 1, 2), *(_h_num_midpt[it->first]),  *(_h_den_midpt[it->first]));
+        histogramFactory().divide(histoPath(dsId, 1, 3), *(_h_num_highpt[it->first]), *(_h_den_highpt[it->first]));
+      }
+    }
+
+    //@}
+
+
+    private:
+
+
+    // Get particle lifetime from hardcoded data
+    double getLifeTime(int pid) {
+      pid = abs(pid);
+      double lft = -1.0;
+      map<int, double>::iterator pPartLft = _partLftMap.find(pid);
+      // search stable particle list
+      if (pPartLft == _partLftMap.end()) {
+        if (pid <= 100) return 0.0;
+        for (size_t i=0; i < sizeof(_stablePDGIds)/sizeof(unsigned int); i++) {
+          if (pid == _stablePDGIds[i]) {
+            lft = 0.0;
+            break;
+          }
+        }
+      } else {
+        lft = (*pPartLft).second;
+      }
+      if (lft < 0.0 && PID::isHadron(pid)) {
+        MSG_ERROR("Could not determine lifetime for particle with PID " << pid
+            << "... This will be considered unprompt!");
+      }
+      return lft;
+    }
+
+    // Data members like post-cuts event weight counters go here
+    const double getMotherLifeTimeSum(const Particle& p) {
+      if ( !p.hasGenParticle() ) return -1.;
+      double lftSum = 0.;
+      double plft = 0.;
+      const GenParticle* part = &(p.genParticle());
+      GenVertex* ivtx = part->production_vertex();
+      while(ivtx)
+      {
+        if (ivtx->particles_in_size() < 1) { lftSum = -1.; break; };
+        const HepMC::GenVertex::particles_in_const_iterator iPart_invtx = ivtx->particles_in_const_begin();
+        part = (*iPart_invtx);
+        if ( !(part) ) { lftSum = -1.; break; };
+        ivtx = part->production_vertex();
+        if ( (part->pdg_id() == 2212) || !(ivtx) ) break; // reached beam
+        plft = getLifeTime(part->pdg_id());
+        if (plft < 0.) { lftSum = -1.; break; };
+        lftSum += plft;
+      };
+      return (lftSum * c_light);
+    }
+
+    /// @name Private variables
+    // Momentum threshold
+    double _p_min;
+
+    // The edges of the intervals of transversal momentum
+    double _pt_min;
+    double _pt1_edge;
+    double _pt2_edge;
+
+    // The limits of the pseudorapidity window
+    int    _eta_nbins;
+    double _eta_min;
+    double _eta_max;
+
+    // Map between PDG id and particle lifetimes in seconds
+    std::map<int, double> _partLftMap;
+
+    // Set of PDG Ids for stable particles (PDG Id <= 100 are considered stable)
+    static const int _stablePDGIds[205];
+
+    // Define histograms
+    // ratio
+    std::map<std::string, LWH::Histogram1D * > _h_ratio_lowpt;
+    std::map<std::string, LWH::Histogram1D * > _h_ratio_midpt;
+    std::map<std::string, LWH::Histogram1D * > _h_ratio_highpt;
+    // numerator
+    std::map<std::string, LWH::Histogram1D * > _h_num_lowpt;
+    std::map<std::string, LWH::Histogram1D * > _h_num_midpt;
+    std::map<std::string, LWH::Histogram1D * > _h_num_highpt;
+    // denominator
+    std::map<std::string, LWH::Histogram1D * > _h_den_lowpt;
+    std::map<std::string, LWH::Histogram1D * > _h_den_midpt;
+    std::map<std::string, LWH::Histogram1D * > _h_den_highpt;
+
+    // Map of ratios and IDs of numerator and denominator
+    std::map<string, pair<int,int> > _ratiotypes;
+    std::map<string, int > _ratiotypes_hepdataid;
+
+    // Fill the PDG Id to Lifetime[seconds] map
+    // Data was extracted from LHCb Particle Table through LHCb::ParticlePropertySvc
+    bool fillMap(map<int, double> &m) {
+      m[6] =  4.707703E-25;  m[11] =  1.E+16;  m[12] =  1.E+16;
+      m[13] =  2.197019E-06;  m[14] =  1.E+16;  m[15] =  2.906E-13;  m[16] =  1.E+16;  m[22] =  1.E+16;
+      m[23] =  2.637914E-25;  m[24] =  3.075758E-25;  m[25] =  9.4E-26;  m[35] =  9.4E-26;
+      m[36] =  9.4E-26;  m[37] =  9.4E-26;  m[84] =  3.335641E-13;  m[85] =  1.290893E-12;
+      m[111] =  8.4E-17;  m[113] =  4.405704E-24;  m[115] =  6.151516E-24;  m[117] =  4.088275E-24;
+      m[119] =  2.102914E-24;  m[130] =  5.116E-08;  m[150] =  1.525E-12;  m[211] =  2.6033E-08;
+      m[213] =  4.405704E-24;  m[215] =  6.151516E-24;  m[217] =  4.088275E-24;  m[219] =  2.102914E-24;
+      m[221] =  5.063171E-19;  m[223] =  7.752794E-23;  m[225] =  3.555982E-24;  m[227] =  3.91793E-24;
+      m[229] =  2.777267E-24;  m[310] =  8.953E-11;  m[313] =  1.308573E-23;  m[315] =  6.038644E-24;
+      m[317] =  4.139699E-24;  m[319] =  3.324304E-24;  m[321] =  1.238E-08;  m[323] =  1.295693E-23;
+      m[325] =  6.682357E-24;  m[327] =  4.139699E-24;  m[329] =  3.324304E-24;  m[331] =  3.210791E-21;
+      m[333] =  1.545099E-22;  m[335] =  9.016605E-24;  m[337] =  7.565657E-24;  m[350] =  1.407125E-12;
+      m[411] =  1.04E-12;  m[413] =  6.856377E-21;  m[415] =  1.778952E-23;  m[421] =  4.101E-13;
+      m[423] =  1.000003E-19;  m[425] =  1.530726E-23;  m[431] =  5.E-13;  m[433] =  1.000003E-19;
+      m[435] =  3.291061E-23;  m[441] =  2.465214E-23;  m[443] =  7.062363E-21;  m[445] =  3.242425E-22;
+      m[510] =  1.525E-12;  m[511] =  1.525E-12;  m[513] =  1.000019E-19;  m[515] =  1.31E-23;
+      m[521] =  1.638E-12;  m[523] =  1.000019E-19;  m[525] =  1.31E-23;  m[530] =  1.536875E-12;
+      m[531] =  1.472E-12;  m[533] =  1.E-19;  m[535] =  1.31E-23;  m[541] =  4.5E-13;
+      m[553] =  1.218911E-20;  m[1112] =  4.539394E-24;  m[1114] =  5.578069E-24;  m[1116] =  1.994582E-24;
+      m[1118] =  2.269697E-24;  m[1212] =  4.539394E-24;  m[1214] =  5.723584E-24;  m[1216] =  1.994582E-24;
+      m[1218] =  1.316424E-24;  m[2112] =  8.857E+02;  m[2114] =  5.578069E-24;  m[2116] =  4.388081E-24;
+      m[2118] =  2.269697E-24;  m[2122] =  4.539394E-24;  m[2124] =  5.723584E-24;  m[2126] =  1.994582E-24;
+      m[2128] =  1.316424E-24;  m[2212] =  1.E+16;  m[2214] =  5.578069E-24;  m[2216] =  4.388081E-24;
+      m[2218] =  2.269697E-24;  m[2222] =  4.539394E-24;  m[2224] =  5.578069E-24;  m[2226] =  1.994582E-24;
+      m[2228] =  2.269697E-24;  m[3112] =  1.479E-10;  m[3114] =  1.670589E-23;  m[3116] =  5.485102E-24;
+      m[3118] =  3.656734E-24;  m[3122] =  2.631E-10;  m[3124] =  4.219309E-23;  m[3126] =  8.227653E-24;
+      m[3128] =  3.291061E-24;  m[3212] =  7.4E-20;  m[3214] =  1.828367E-23;  m[3216] =  5.485102E-24;
+      m[3218] =  3.656734E-24;  m[3222] =  8.018E-11;  m[3224] =  1.838582E-23;  m[3226] =  5.485102E-24;
+      m[3228] =  3.656734E-24;  m[3312] =  1.639E-10;  m[3314] =  6.648608E-23;  m[3322] =  2.9E-10;
+      m[3324] =  7.233101E-23;  m[3334] =  8.21E-11;  m[4112] =  2.991874E-22;  m[4114] =  4.088274E-23;
+      m[4122] =  2.E-13;  m[4132] =  1.12E-13;  m[4212] =  3.999999E-22;  m[4214] =  3.291061E-22;
+      m[4222] =  2.951624E-22;  m[4224] =  4.417531E-23;  m[4232] =  4.42E-13;  m[4332] =  6.9E-14;
+      m[4412] =  3.335641E-13;  m[4422] =  3.335641E-13;  m[4432] =  3.335641E-13;  m[5112] =  1.E-19;
+      m[5122] =  1.38E-12;  m[5132] =  1.42E-12;  m[5142] =  1.290893E-12;  m[5212] =  1.E-19;
+      m[5222] =  1.E-19;  m[5232] =  1.42E-12;  m[5242] =  1.290893E-12;  m[5312] =  1.E-19;
+      m[5322] =  1.E-19;  m[5332] =  1.55E-12;  m[5342] =  1.290893E-12;  m[5442] =  1.290893E-12;
+      m[5512] =  1.290893E-12;  m[5522] =  1.290893E-12;  m[5532] =  1.290893E-12;  m[5542] =  1.290893E-12;
+      m[10111] =  2.48382E-24;  m[10113] =  4.635297E-24;  m[10115] =  2.54136E-24;  m[10211] =  2.48382E-24;
+      m[10213] =  4.635297E-24;  m[10215] =  2.54136E-24;  m[10223] =  1.828367E-24;  m[10225] =  3.636531E-24;
+      m[10311] =  2.437823E-24;  m[10313] =  7.313469E-24;  m[10315] =  3.538775E-24;
+      m[10321] =  2.437823E-24;  m[10323] =  7.313469E-24;  m[10325] =  3.538775E-24;
+      m[10331] =  4.804469E-24;  m[10411] =  4.38E-24;  m[10413] =  3.29E-23;  m[10421] =  4.38E-24;
+      m[10423] =  3.22653E-23;  m[10431] =  6.5821E-22;  m[10433] =  6.5821E-22;  m[10441] =  6.453061E-23;
+      m[10511] =  4.39E-24;  m[10513] =  1.65E-23;  m[10521] =  4.39E-24;  m[10523] =  1.65E-23;
+      m[10531] =  4.39E-24;  m[10533] =  1.65E-23;  m[11114] =  2.194041E-24;  m[11116] =  1.828367E-24;
+      m[11212] =  1.880606E-24;  m[11216] =  1.828367E-24;  m[12112] =  2.194041E-24;
+      m[12114] =  2.194041E-24;  m[12116] =  5.063171E-24;  m[12126] =  1.828367E-24;
+      m[12212] =  2.194041E-24;  m[12214] =  2.194041E-24;  m[12216] =  5.063171E-24;
+      m[12224] =  2.194041E-24;  m[12226] =  1.828367E-24;  m[13112] =  6.582122E-24;  m[13114] =  1.09702E-23;
+      m[13116] =  5.485102E-24;  m[13122] =  1.316424E-23;  m[13124] =  1.09702E-23;  m[13126] =  6.928549E-24;
+      m[13212] =  6.582122E-24;  m[13214] =  1.09702E-23;  m[13216] =  5.485102E-24;  m[13222] =  6.582122E-24;
+      m[13224] =  1.09702E-23;  m[13226] =  5.485102E-24;  m[13314] =  2.742551E-23;
+      m[13324] =  2.742551E-23;  m[14122] =  1.828367E-22;  m[20022] =  1.E+16;  m[20113] =  1.567172E-24;
+      m[20213] =  1.567172E-24;  m[20223] =  2.708692E-23;  m[20313] =  3.782829E-24;
+      m[20315] =  2.384827E-24;  m[20323] =  3.782829E-24;  m[20325] =  2.384827E-24;
+      m[20333] =  1.198929E-23;  m[20413] =  2.63E-24;  m[20423] =  2.63E-24;  m[20433] =  6.5821E-22;
+      m[20443] =  7.395643E-22;  m[20513] =  2.63E-24;  m[20523] =  2.63E-24;  m[20533] =  2.63E-24;
+      m[21112] =  2.632849E-24;  m[21114] =  3.291061E-24;  m[21212] =  2.632849E-24;
+      m[21214] =  6.582122E-24;  m[22112] =  4.388081E-24;  m[22114] =  3.291061E-24;
+      m[22122] =  2.632849E-24;  m[22124] =  6.582122E-24;  m[22212] =  4.388081E-24;
+      m[22214] =  3.291061E-24;  m[22222] =  2.632849E-24;  m[22224] =  3.291061E-24;
+      m[23112] =  7.313469E-24;  m[23114] =  2.991874E-24;  m[23122] =  4.388081E-24;
+      m[23124] =  6.582122E-24;  m[23126] =  3.291061E-24;  m[23212] =  7.313469E-24;
+      m[23214] =  2.991874E-24;  m[23222] =  7.313469E-24;  m[23224] =  2.991874E-24;
+      m[30113] =  2.632849E-24;  m[30213] =  2.632849E-24;  m[30221] =  1.880606E-24;
+      m[30223] =  2.089563E-24;  m[30313] =  2.056913E-24;  m[30323] =  2.056913E-24;
+      m[30443] =  2.419898E-23;  m[31114] =  1.880606E-24;  m[31214] =  3.291061E-24;
+      m[32112] =  3.989164E-24;  m[32114] =  1.880606E-24;  m[32124] =  3.291061E-24;
+      m[32212] =  3.989164E-24;  m[32214] =  1.880606E-24;  m[32224] =  1.880606E-24;
+      m[33122] =  1.880606E-23;  m[42112] =  6.582122E-24;  m[42212] =  6.582122E-24;
+      m[43122] =  2.194041E-24;  m[53122] =  4.388081E-24;  m[100111] =  1.645531E-24;
+      m[100113] =  1.64553E-24;  m[100211] =  1.645531E-24;  m[100213] =  1.64553E-24;
+      m[100221] =  1.196749E-23;  m[100223] =  3.061452E-24;  m[100313] =  2.837122E-24;
+      m[100323] =  2.837122E-24;  m[100331] =  4.459432E-25;  m[100333] =  4.388081E-24;
+      m[100441] =  4.701516E-23;  m[100443] =  2.076379E-21;  m[100553] =  2.056913E-20;
+      m[200553] =  3.242425E-20;  m[300553] =  3.210791E-23;  m[9000111] =  8.776163E-24;
+      m[9000211] =  8.776163E-24;  m[9000443] =  8.227652E-24;  m[9000553] =  5.983747E-24;
+      m[9010111] =  3.164482E-24;  m[9010211] =  3.164482E-24;  m[9010221] =  9.403031E-24;
+      m[9010443] =  8.438618E-24;  m[9010553] =  8.3318E-24;  m[9020443] =  1.061633E-23;
+      m[9030221] =  6.038644E-24;  m[9042413] =  2.07634E-21;  m[9050225] =  1.394517E-24;
+      m[9060225] =  3.291061E-24;  m[9080225] =  4.388081E-24;  m[9090225] =  2.056913E-24;
+      m[9910445] =  2.07634E-21;  m[9920443] =  2.07634E-21;
+      return true;
+    }
+
+  };
+
+
+  const int LHCB_2012_I1119400::_stablePDGIds[205] = {
+    311, 543, 545, 551, 555, 557, 1103, 2101, 2103, 2203, 3101, 3103, 3201, 3203, 3303,
+    4101, 4103, 4124, 4201, 4203, 4301, 4303, 4312, 4314, 4322, 4324, 4334, 4403, 4414,
+    4424, 4434, 4444, 5101, 5103, 5114, 5201, 5203, 5214, 5224, 5301, 5303, 5314, 5324,
+    5334, 5401, 5403, 5412, 5414, 5422, 5424, 5432, 5434, 5444, 5503, 5514, 5524, 5534,
+    5544, 5554, 10022, 10333, 10335, 10443, 10541, 10543, 10551, 10553, 10555, 11112,
+    12118, 12122, 12218, 12222, 13316, 13326, 20543, 20553, 20555, 23314, 23324, 30343,
+    30353, 30363, 30553, 33314, 33324, 41214, 42124, 52114, 52214, 100311, 100315, 100321,
+    100325, 100411, 100413, 100421, 100423, 100551, 100555, 100557, 110551, 110553, 110555,
+    120553, 120555, 130553, 200551, 200555, 210551, 210553, 220553, 1000001, 1000002,
+    1000003, 1000004, 1000005, 1000006, 1000011, 1000012, 1000013, 1000014, 1000015,
+    1000016, 1000021, 1000022, 1000023, 1000024, 1000025, 1000035, 1000037, 1000039,
+    2000001, 2000002, 2000003, 2000004, 2000005, 2000006, 2000011, 2000012, 2000013,
+    2000014, 2000015, 2000016, 3000111, 3000113, 3000211, 3000213, 3000221, 3000223,
+    3000331, 3100021, 3100111, 3100113, 3200111, 3200113, 3300113, 3400113, 4000001,
+    4000002, 4000011, 4000012, 5000039, 9000221, 9900012, 9900014, 9900016, 9900023,
+    9900024, 9900041, 9900042 };
+
+
+  // Plugin hook
+  DECLARE_RIVET_PLUGIN(LHCB_2012_I1119400);
+
+}

Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/Makefile.am	Fri Dec 14 13:38:42 2012	(r4083)
+++ branches/2012-06-aidarivet/src/Analyses/Makefile.am	Fri Dec 14 13:57:50 2012	(r4084)
@@ -145,7 +145,8 @@
 RivetLHCbAnalyses_la_SOURCES += \
     LHCB_2010_S8758301.cc \
     LHCB_2011_I917009.cc \
-    LHCB_2011_I919315.cc
+    LHCB_2011_I919315.cc \
+    LHCB_2012_I1119400.cc
 endif
 
 


More information about the Rivet-svn mailing list