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

blackhole at projects.hepforge.org blackhole at projects.hepforge.org
Sat Feb 9 11:22:52 GMT 2013


Author: richardn
Date: Sat Feb  9 11:22:51 2013
New Revision: 4134

Log:
added BaBar charmonium production analysis

Added:
   branches/2012-06-aidarivet/data/anainfo/BABAR_2003_I593379.info
      - copied, changed from r4111, branches/2012-06-aidarivet/data/anainfo/BABAR_2007_S7266081.info
   branches/2012-06-aidarivet/data/plotinfo/BABAR_2003_I593379.plot
      - copied, changed from r4111, branches/2012-06-aidarivet/data/plotinfo/BABAR_2005_S6181155.plot
   branches/2012-06-aidarivet/data/refdata/BABAR_2003_I593379.aida
   branches/2012-06-aidarivet/src/Analyses/BABAR_2003_I593379.cc
      - copied, changed from r4111, branches/2012-06-aidarivet/src/Analyses/BELLE_2001_S4598261.cc
Modified:
   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

Copied and modified: branches/2012-06-aidarivet/data/anainfo/BABAR_2003_I593379.info (from r4111, branches/2012-06-aidarivet/data/anainfo/BABAR_2007_S7266081.info)
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/BABAR_2007_S7266081.info	Fri Jan 25 18:56:20 2013	(r4111, copy source)
+++ branches/2012-06-aidarivet/data/anainfo/BABAR_2003_I593379.info	Sat Feb  9 11:22:51 2013	(r4134)
@@ -1,37 +1,37 @@
-Name: BABAR_2007_S7266081
-Summary: Measurements of Semi-Leptonic Tau Decays into Three Charged Hadrons
+Name: BABAR_2003_I593379
+Summary: Measurement of inclusive charmonium production
 Status: VALIDATED
 Authors:
  - Peter Richardson <Peter.Richardson at durham.ac.uk>
 References:
- - Phys.Rev.Lett.100:011801,2008
- - arXiv:0707.2981
- - SLAC-R-936
+ - Phys.Rev. D67 032002, 2003
+ - hep-ex/0207097
 RunInfo:
-  Tau production, can be any process but original data was in  $e^+ e^-$ at the $\Upsilon(4S)$ resonance, with CoM boost -- 8.0~GeV~($e^−$) and 3.5~GeV~($e^+$)
+  Production of charmionum at the $\Upslion(4S)$ resonace.
 NumEvents: 1000000
-SpiresID: 7266081
+InspireID: 593379
 Beams: [e+, e-]
 Energies: [[3.5, 8.0]]
 PtCuts: [0]
-Description: Measurement of tau decays to three charged hadrons using
- a data sample corresponding to an integrated luminosity of 342 $fb^{-1}$ collected
+Description: Measurement of  $J/\Psi$, $\Psi'$, $\chi_{c1}$ and 
+    $\chi_{c2}$ production using
+ a data sample corresponding to an integrated luminosity of 20.3$\mathrm{fb}^{-1}$ collected
  with the BABAR detector at the SLAC PEP-II electron-positron storage ring operating
- at a center-of-mass energy near 10.58 GeV.
-BibKey: Aubert:2007mh
-BibTeX: '@Article{Aubert:2007mh,
-     author    = "Aubert, Bernard and others",
-     collaboration = "BABAR",
-     title     = "{Exclusive branching fraction measurements of semileptonic
-                  tau decays into three charged hadrons, $\tau^- \to \phi
-                  \pi^- \nu_\tau$ and $\tau^- \to \phi K^- \nu_\tau$}",
-     journal   = "Phys. Rev. Lett.",
-     volume    = "100",
-     year      = "2008",
-     pages     = "011801",
-     eprint    = "0707.2981",
-     archivePrefix = "arXiv",
-     primaryClass  =  "hep-ex",
-     doi       = "10.1103/PhysRevLett.100.011801",
-     SLACcitation  = "%%CITATION = 0707.2981;%%"
+ at a centre-of-mass energy near 10.58 GeV.
+BibKey: Aubert:2002hc
+BibTeX: '@article{Aubert:2002hc,
+      author         = "Aubert, Bernard and others",
+      title          = "{Study of inclusive production of charmonium mesons in
+                        $B$ decay}",
+      collaboration  = "BABAR Collaboration",
+      journal        = "Phys.Rev.",
+      volume         = "D67",
+      pages          = "032002",
+      doi            = "10.1103/PhysRevD.67.032002",
+      year           = "2003",
+      eprint         = "hep-ex/0207097",
+      archivePrefix  = "arXiv",
+      primaryClass   = "hep-ex",
+      reportNumber   = "SLAC-PUB-9327, BABAR-PUB-02-04",
+      SLACcitation   = "%%CITATION = HEP-EX/0207097;%%",
 }'

Modified: branches/2012-06-aidarivet/data/anainfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/anainfo/Makefile.am	Fri Feb  8 09:49:59 2013	(r4133)
+++ branches/2012-06-aidarivet/data/anainfo/Makefile.am	Sat Feb  9 11:22:51 2013	(r4134)
@@ -74,6 +74,7 @@
   ATLAS_2012_CONF_2012_109.info \
   ATLAS_2012_CONF_2012_153.info \
   ATLAS_2013_I1217867.info \
+  BABAR_2003_I593379.info \
   BABAR_2005_S6181155.info \
   BABAR_2007_S6895344.info \
   BABAR_2007_S7266081.info \

Copied and modified: branches/2012-06-aidarivet/data/plotinfo/BABAR_2003_I593379.plot (from r4111, branches/2012-06-aidarivet/data/plotinfo/BABAR_2005_S6181155.plot)
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/BABAR_2005_S6181155.plot	Fri Jan 25 18:56:20 2013	(r4111, copy source)
+++ branches/2012-06-aidarivet/data/plotinfo/BABAR_2003_I593379.plot	Sat Feb  9 11:22:51 2013	(r4134)
@@ -1,42 +1,72 @@
-# BEGIN PLOT /BABAR_2005_S6181155/d01-x01-y01
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ with $\Xi_c^0\to\Xi^-\pi^+$ at the $\Upsilon(4S)$
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [fb/GeV]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y01
+Title=${\rm Br}(B\to J/\Psi)$ at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to J/\Psi)$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d02-x01-y01
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ with $\Xi_c^0\to\Xi^-\pi^+$ at the $\Upsilon(4S)$
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [fb/GeV]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y02
+Title=${\rm Br}(B\to J/\Psi)$~direct at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to J/\Psi)$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d02-x01-y02
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ with $\Xi_c^0\to\Xi^-\pi^+$ in the continuum region
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [fb/GeV]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y03
+Title=${\rm Br}(B\to \chi_{c1})$ at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to \chi_{c1})$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d03-x01-y01
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ with $\Xi_c^0\to\Xi^-\pi^+$ 
-XLabel=$\sqrt{s}$ [GeV]
-YLabel=$\sigma$ [fb]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y04
+Title=${\rm Br}(B\to \chi_{c1})$~direct at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to \chi_{c1})$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d04-x01-y01
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ at the $\Upsilon(4S)$
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$1/\sigma\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [$\text{GeV}^{-1}$]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y05
+Title=${\rm Br}(B\to \chi_{c2})$ at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to \chi_{c2})$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d05-x01-y01
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ at the $\Upsilon(4S)$
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$1/\sigma\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [$\text{GeV}^{-1}$]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y06
+Title=${\rm Br}(B\to \chi_{c2})$~direct at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to \chi_{c2})$ [\%]
 LogY=0
 # END PLOT
-# BEGIN PLOT /BABAR_2005_S6181155/d05-x01-y02
-Title=$\sigma(e^+e^-\to \Xi_c^0+\bar{\Xi}_c^0+X)$ in the continuum region
-XLabel=$p_{\Xi_c^0}$ [GeV]
-YLabel=$1/\sigma\text{d}\sigma/\text{d}p_{\Xi_c^0}$ [$\text{GeV}^{-1}$]
+# BEGIN PLOT /BABAR_2003_I593379/d01-x01-y07
+Title=${\rm Br}(B\to \Psi')$ at the $\Upsilon(4S)$
+XLabel=
+YLabel=${\rm Br}(B\to \Psi')$ [\%]
 LogY=0
 # END PLOT
+# BEGIN PLOT /BABAR_2003_I593379/d06-x01-y01
+Title=${\rm Br}(B\to J/\Psi)$ at the $\Upsilon(4S)$
+XLabel=$p^*_{J/\Psi}$ [GeV]
+YLabel=$\text{d}B/\text{d}p^*_{J/\Psi}$ [0.1/GeV]
+LogY=0
+# END PLOT
+# BEGIN PLOT /BABAR_2003_I593379/d07-x01-y01
+Title=${\rm Br}(B\to \chi_{c1})$ at the $\Upsilon(4S)$
+XLabel=$p^*_{\chi_{c1}}$ [GeV]
+YLabel=$\text{d}B/\text{d}p^*_{\chi_{c1}}$ [0.1/GeV]
+LogY=0
+# END PLOT
+# BEGIN PLOT /BABAR_2003_I593379/d07-x01-y02
+Title=${\rm Br}(B\to \chi_{c2})$ at the $\Upsilon(4S)$
+XLabel=$p^*_{\chi_{c2}}$ [GeV]
+YLabel=$\text{d}B/\text{d}p^*_{\chi_{c2}}$ [0.1/GeV]
+LogY=0
+# END PLOT
+# BEGIN PLOT /BABAR_2003_I593379/d08-x01-y01
+Title=${\rm Br}(B\to \Psi')$ at the $\Upsilon(4S)$
+XLabel=$p^*_{\Psi'}$ [GeV]
+YLabel=$\text{d}B/\text{d}p^*_{\Psi'}$ [0.1/GeV]
+LogY=0
+# END PLOT
+# BEGIN PLOT /BABAR_2003_I593379/d10-x01-y01
+Title=${\rm Br}(B\to J/\Psi)$~(direct) at the $\Upsilon(4S)$
+XLabel=$p^*_{J/\Psi}$ [GeV]
+YLabel=$\text{d}B/\text{d}p^*_{J/\Psi}$ [0.1/GeV]
+LogY=0
+# END PLOT
\ No newline at end of file

Modified: branches/2012-06-aidarivet/data/plotinfo/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Fri Feb  8 09:49:59 2013	(r4133)
+++ branches/2012-06-aidarivet/data/plotinfo/Makefile.am	Sat Feb  9 11:22:51 2013	(r4134)
@@ -71,6 +71,7 @@
   ATLAS_2012_CONF_2012_109.plot \
   ATLAS_2012_CONF_2012_153.plot \
   ATLAS_2013_I1217867.plot \
+  BABAR_2003_I593379.plot \
   BABAR_2005_S6181155.plot \
   BABAR_2007_S6895344.plot \
   BABAR_2007_S7266081.plot \

Added: branches/2012-06-aidarivet/data/refdata/BABAR_2003_I593379.aida
==============================================================================
--- /dev/null	00:00:00 1970	(empty, because file is newly added)
+++ branches/2012-06-aidarivet/data/refdata/BABAR_2003_I593379.aida	Sat Feb  9 11:22:51 2013	(r4134)
@@ -0,0 +1,417 @@
+<?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="d06-x01-y01" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="4.69838884761e-06" errorPlus="1.15981979842e-05" errorMinus="1.16037222509e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="6.97770129737e-05" errorPlus="1.16009601176e-05" errorMinus="1.16009601175e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000158063081602" errorPlus="1.1598197984e-05" errorMinus="1.3921152141e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000313629194645" errorPlus="1.6241344164e-05" errorMinus="1.8561536188e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000432074997445" errorPlus="1.8561536188e-05" errorMinus="2.3201920235e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000515717919893" errorPlus="2.0884490345e-05" errorMinus="2.5522112259e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00055760014805" errorPlus="2.5516587992e-05" errorMinus="2.7842304282e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000703885492999" errorPlus="2.3201920235e-05" errorMinus="2.5519350125e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000743447529134" errorPlus="2.5516587992e-05" errorMinus="3.0159734173e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000708760658382" errorPlus="2.7839542149e-05" errorMinus="2.7845066415e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000769198898461" errorPlus="2.7842304282e-05" errorMinus="2.7842304282e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000776278246266" errorPlus="2.7839542149e-05" errorMinus="2.5519350125e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000822795334204" errorPlus="2.7842304282e-05" errorMinus="3.0165258439e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000818270959758" errorPlus="2.5522112259e-05" errorMinus="2.7845066415e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000964561828974" errorPlus="3.0159734173e-05" errorMinus="3.0162496305e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000839384707172" errorPlus="2.7842304283e-05" errorMinus="3.0162496305e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000718850731551" errorPlus="2.7839542149e-05" errorMinus="2.5519350125e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00051015222117" errorPlus="2.3199158101e-05" errorMinus="2.7845066416e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00024112319391" errorPlus="2.0881728212e-05" errorMinus="2.5519350125e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-7.02134300448e-06" errorPlus="2.32019202351e-05" errorMinus="2.31991581017e-05"/>
+</dataPoint>
+  </dataPointSet>
+<dataPointSet name="d07-x01-y01" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="7.89771449333e-06" errorPlus="1.04804581394e-05" errorMinus="9.17040087199e-06"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="5.25268603263e-05" errorPlus="2.7511202616e-05" errorMinus="7.61076463525e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="1.46244744719e-05" errorPlus="3.79895845949e-05" errorMinus="-1.16991643454e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000128684579318" errorPlus="4.8472118895e-05" errorMinus="4.97821761622e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00019427464143" errorPlus="6.4192806104e-05" errorMinus="6.5504939532e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000320127337843" errorPlus="7.7293378778e-05" errorMinus="7.7293378778e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000187898752574" errorPlus="6.8122977906e-05" errorMinus="6.9433035174e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00012903337428" errorPlus="7.2053149708e-05" errorMinus="7.33632069761e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000413403000052" errorPlus="8.3841588955e-05" errorMinus="8.3841588955e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000353229640651" errorPlus="8.5151646222e-05" errorMinus="8.7773836917e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000392616481254" errorPlus="7.8601359885e-05" errorMinus="7.8603436046e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.0003547140954" errorPlus="7.8601359885e-05" errorMinus="7.8605512207e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000336458416236" errorPlus="7.0743092441e-05" errorMinus="7.2053149709e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000375847332999" errorPlus="6.9433035173e-05" errorMinus="6.9435111334e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000216109621274" errorPlus="5.7642519767e-05" errorMinus="5.7642519767e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000169032682226" errorPlus="5.7642519767e-05" errorMinus="5.8952577034e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="1.22493468745e-07" errorPlus="4.97801000018e-05" errorMinus="4.84721188947e-05"/>
+</dataPoint>
+  </dataPointSet>
+<dataPointSet name="d07-x01-y02" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="5.24017467249e-06" errorPlus="1.17903930131e-05" errorMinus="1.17903930131e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="3.27510917031e-05" errorPlus="2.48908296943e-05" errorMinus="2.62008733625e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-3.93013100437e-06" errorPlus="4.97816593887e-05" errorMinus="4.58515283842e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00011135371179" errorPlus="7.2052401747e-05" errorMinus="7.20524017463e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000168995633188" errorPlus="7.8602620087e-05" errorMinus="7.86026200876e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000205676855895" errorPlus="8.9082969433e-05" errorMinus="8.9082969432e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000144104803493" errorPlus="9.1703056769e-05" errorMinus="9.30131004362e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000117903930131" errorPlus="0.00010480349345" errorMinus="0.000100873362445"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00016768558952" errorPlus="0.000103493449781" errorMinus="0.000106113537118"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000203056768559" errorPlus="0.000102183406113" errorMinus="0.000100873362445"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000267248908297" errorPlus="0.000102183406113" errorMinus="0.000102183406114"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="3.01310043668e-05" errorPlus="0.00010480349345" errorMinus="0.000100873362445"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000184716157205" errorPlus="0.00010480349345" errorMinus="0.000106113537118"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000221397379913" errorPlus="9.1703056768e-05" errorMinus="9.0393013101e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="7.4672489083e-05" errorPlus="9.4323144105e-05" errorMinus="9.5633187773e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-7.4672489083e-05" errorPlus="9.69432314411e-05" errorMinus="0.000107423580786"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-1.70305676856e-05" errorPlus="0.000136244541485" errorMinus="0.000133624454148"/>
+</dataPoint>
+  </dataPointSet>
+<dataPointSet name="d08-x01-y01" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="6.32545081525e-06" errorPlus="4.53514739225e-06" errorMinus="4.5340675953e-06"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="2.3514739229e-05" errorPlus="1.08854335385e-05" errorMinus="1.08843537415e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="8.15192743764e-05" errorPlus="1.99546485256e-05" errorMinus="1.99546485261e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000138616780045" errorPlus="2.3583846237e-05" errorMinus="2.4489795918e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000173038548753" errorPlus="2.9931972789e-05" errorMinus="3.1744951949e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000235578231293" errorPlus="3.2653061224e-05" errorMinus="3.35611705e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000229182593672" errorPlus="3.5375229457e-05" errorMinus="3.5374149659e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000214624770543" errorPlus="3.5375229457e-05" errorMinus="3.8093078501e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000254488716121" errorPlus="3.7189288414e-05" errorMinus="3.7187128819e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000222698412698" errorPlus="3.8094158299e-05" errorMinus="3.9002267573e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000297029478458" errorPlus="3.8095238095e-05" errorMinus="3.8095238095e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000296076017709" errorPlus="3.9002267573e-05" errorMinus="4.0816326531e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000325962639024" errorPlus="3.9911456646e-05" errorMinus="4.1723356009e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000216167800454" errorPlus="3.3560090702e-05" errorMinus="3.5373069863e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000119977324263" errorPlus="3.4466040384e-05" errorMinus="3.44671201814e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-9.77324263039e-06" errorPlus="2.53968253968e-05" errorMinus="2.44897959184e-05"/>
+</dataPoint>
+  </dataPointSet>
+<dataPointSet name="d10-x01-y01" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="0.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="2.25173210162e-06" errorPlus="1.15473441109e-05" errorMinus="1.15473441109e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="1.59930715935e-05" errorPlus="1.84757505774e-05" errorMinus="1.84757505773e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="8.51589134499e-05" errorPlus="2.54069064111e-05" errorMinus="2.54014076762e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="8.04272517321e-05" errorPlus="3.69515011549e-05" errorMinus="3.46420323326e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00020733256351" errorPlus="4.1570438799e-05" errorMinus="3.9260969977e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000227999560101" errorPlus="4.6192125811e-05" errorMinus="4.6186627076e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00033411965248" errorPlus="4.3882656989e-05" errorMinus="4.3879907621e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000350173210162" errorPlus="4.8501594633e-05" errorMinus="5.0811063456e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000384697019685" errorPlus="4.8501594634e-05" errorMinus="5.0808314087e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="0.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000391509952711" errorPlus="4.8501594633e-05" errorMinus="4.849334653e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.05" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000467609699769" errorPlus="4.8501594633e-05" errorMinus="5.0808314088e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.15" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00048827669636" errorPlus="5.3120532277e-05" errorMinus="4.6186627076e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.25" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000605946882217" errorPlus="4.8498845266e-05" errorMinus="5.0808314088e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.35" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000712064225228" errorPlus="4.3885406357e-05" errorMinus="4.6186627076e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.45" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000868992631695" errorPlus="5.0811063455e-05" errorMinus="5.0802815353e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.55" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000792664687122" errorPlus="6.2358407566e-05" errorMinus="6.2352908831e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.65" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000684006928406" errorPlus="6.0048938745e-05" errorMinus="6.0043440008e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.75" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.00049451226218" errorPlus="3.9263719344e-05" errorMinus="4.1567689432e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.85" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="0.000226501154734" errorPlus="4.1570438799e-05" errorMinus="4.3879907621e-05"/>
+</dataPoint>
+<dataPoint>
+<measurement value="1.95" errorPlus="0.05" errorMinus="0.05"/>
+<measurement value="-1.61112944023e-05" errorPlus="3.00285934235e-05" errorMinus="3.23270647751e-05"/>
+</dataPoint>
+  </dataPointSet>
+<dataPointSet name="d01-x01-y01" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="1.057" errorPlus="0.0417612260356" errorMinus="0.0417612260356"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y02" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.74" errorPlus="0.0487647413609" errorMinus="0.0487647413609"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y03" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.21" errorPlus="0.0546443043693" errorMinus="0.0546443043693"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y04" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.19" errorPlus="0.0535350352573" errorMinus="0.0535350352573"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y05" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.21" errorPlus="0.0546443043693" errorMinus="0.0546443043693"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y06" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.19" errorPlus="0.0535350352573" errorMinus="0.0535350352573"/>
+</dataPoint>
+</dataPointSet>
+<dataPointSet name="d01-x01-y07" dimension="2" path="/REF/BABAR_2003_I593379" title="B" >
+<dataPoint>
+<measurement value="10.58" errorPlus="0.5" errorMinus="0.5"/>
+<measurement value="0.297" errorPlus="0.0282842712475" errorMinus="0.0282842712475"/>
+</dataPoint>
+</dataPointSet>
+</aida>

Modified: branches/2012-06-aidarivet/data/refdata/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/data/refdata/Makefile.am	Fri Feb  8 09:49:59 2013	(r4133)
+++ branches/2012-06-aidarivet/data/refdata/Makefile.am	Sat Feb  9 11:22:51 2013	(r4134)
@@ -62,6 +62,7 @@
   ATLAS_2012_I1188891.aida \
   ATLAS_2012_CONF_2012_001.aida \
   ATLAS_2013_I1217867.aida \
+  BABAR_2003_I593379.aida \
   BABAR_2005_S6181155.aida \
   BABAR_2006_S6511112.aida \
   BABAR_2007_S6895344.aida \

Copied and modified: branches/2012-06-aidarivet/src/Analyses/BABAR_2003_I593379.cc (from r4111, branches/2012-06-aidarivet/src/Analyses/BELLE_2001_S4598261.cc)
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/BELLE_2001_S4598261.cc	Fri Jan 25 18:56:20 2013	(r4111, copy source)
+++ branches/2012-06-aidarivet/src/Analyses/BABAR_2003_I593379.cc	Sat Feb  9 11:22:51 2013	(r4134)
@@ -9,21 +9,19 @@
 namespace Rivet {
 
 
-  /// @brief BELLE pi0 spectrum at Upsilon(4S)
+  /// @brief Babar chamronium spectra
   /// @author Peter Richardson
-  class BELLE_2001_S4598261 : public Analysis {
+  class BABAR_2003_I593379 : public Analysis {
   public:
 
-    BELLE_2001_S4598261() 
-      : Analysis("BELLE_2001_S4598261"), _weightSum(0.)
+    BABAR_2003_I593379() 
+      : Analysis("BABAR_2003_I593379"), _weightSum(0.)
     { }
 
 
     void analyze(const Event& e) {
       const double weight = e.weight();
-
       const UnstableFinalState& ufs = applyProjection<UnstableFinalState>(e, "UFS");
-
       // find the upsilons
       ParticleVector upsilons;
       // first in unstable final state
@@ -51,33 +49,79 @@
       // find an upsilons
       foreach (const Particle& p, upsilons) {
 	_weightSum += weight;
-	// find the neutral pions from the decay
-	vector<GenParticle *> pions;
-	findDecayProducts(p.genParticle(),pions);
+	// find the charmonium resonances
+	vector<GenParticle *> allJpsi,primaryJpsi,Psiprime,
+	  all_chi_c1,all_chi_c2,primary_chi_c1,primary_chi_c2;
+	findDecayProducts(p.genParticle(),allJpsi,primaryJpsi,Psiprime,
+			  all_chi_c1,all_chi_c2,primary_chi_c1,primary_chi_c2);
 	LorentzTransform cms_boost(-p.momentum().boostVector());
-	for(unsigned int ix=0;ix<pions.size();++ix) {
+	for(unsigned int ix=0;ix<allJpsi.size();++ix) {
+	  double pcm = 
+	    cms_boost.transform(FourMomentum(allJpsi[ix]->momentum())).vector3().mod();
+	  _hist_all_Jpsi->fill(pcm,weight);
+	}
+	_mult_JPsi->fill(10.58,weight*double(allJpsi.size()));
+	for(unsigned int ix=0;ix<primaryJpsi.size();++ix) {
+	  double pcm = 
+	    cms_boost.transform(FourMomentum(primaryJpsi[ix]->momentum())).vector3().mod();
+	  _hist_primary_Jpsi->fill(pcm,weight);
+	}
+	_mult_JPsi_direct->fill(10.58,weight*double(primaryJpsi.size()));
+	for(unsigned int ix=0;ix<Psiprime.size();++ix) {
+	  double pcm = 
+	    cms_boost.transform(FourMomentum(Psiprime[ix]->momentum())).vector3().mod();
+	  _hist_Psi_prime->fill(pcm,weight);
+	}
+	_mult_Psi2S->fill(10.58,weight*double(Psiprime.size()));
+	for(unsigned int ix=0;ix<all_chi_c1.size();++ix) {
 	  double pcm = 
-	    cms_boost.transform(FourMomentum(pions[ix]->momentum())).vector3().mod();
-	  _histdSigDp->fill(pcm,weight);
+	    cms_boost.transform(FourMomentum(all_chi_c1[ix]->momentum())).vector3().mod();
+	  _hist_chi_c1->fill(pcm,weight);
 	}
-	_histMult->fill(0.,double(pions.size())*weight);
+	_mult_chi_c1->fill(10.58,weight*double(all_chi_c1.size()));
+	_mult_chi_c1_direct->fill(10.58,weight*double(primary_chi_c1.size()));
+	for(unsigned int ix=0;ix<all_chi_c2.size();++ix) {
+	  double pcm = 
+	    cms_boost.transform(FourMomentum(all_chi_c2[ix]->momentum())).vector3().mod();
+	  _hist_chi_c2->fill(pcm,weight);
+	}
+	_mult_chi_c2->fill(10.58,weight*double(all_chi_c2.size()));
+	_mult_chi_c2_direct->fill(10.58,weight*double(primary_chi_c2.size()));
       }
     } // analyze
 
     void finalize() {
 
-      scale(_histdSigDp, 1./_weightSum);
-      scale(_histMult  , 1./_weightSum);
+      scale(_hist_all_Jpsi    ,0.5*0.1/_weightSum);
+      scale(_hist_chi_c1      ,0.5*0.1/_weightSum);
+      scale(_hist_chi_c2      ,0.5*0.1/_weightSum);
+      scale(_hist_Psi_prime   ,0.5*0.1/_weightSum);
+      scale(_hist_primary_Jpsi,0.5*0.1/_weightSum);
+      scale(_mult_JPsi         ,0.5*100./_weightSum);
+      scale(_mult_JPsi_direct  ,0.5*100./_weightSum);
+      scale(_mult_chi_c1       ,0.5*100./_weightSum);
+      scale(_mult_chi_c1_direct,0.5*100./_weightSum);
+      scale(_mult_chi_c2       ,0.5*100./_weightSum);
+      scale(_mult_chi_c2_direct,0.5*100./_weightSum);
+      scale(_mult_Psi2S        ,0.5*100./_weightSum);
     } // finalize
 
 
     void init() {
       addProjection(UnstableFinalState(), "UFS");
 
-      // spectrum
-      _histdSigDp = bookHistogram1D(1, 1, 1);
-      // multiplicity
-      _histMult   = bookHistogram1D(2, 1, 1);
+      _mult_JPsi          = bookHistogram1D(1, 1, 1);
+      _mult_JPsi_direct   = bookHistogram1D(1, 1, 2);
+      _mult_chi_c1        = bookHistogram1D(1, 1, 3);
+      _mult_chi_c1_direct = bookHistogram1D(1, 1, 4);
+      _mult_chi_c2        = bookHistogram1D(1, 1, 5);
+      _mult_chi_c2_direct = bookHistogram1D(1, 1, 6);
+      _mult_Psi2S         = bookHistogram1D(1, 1, 7);
+      _hist_all_Jpsi      = bookHistogram1D(6, 1, 1);
+      _hist_chi_c1        = bookHistogram1D(7, 1, 1);
+      _hist_chi_c2        = bookHistogram1D(7, 1, 2);
+      _hist_Psi_prime     = bookHistogram1D(8, 1, 1);
+      _hist_primary_Jpsi  = bookHistogram1D(10, 1, 1);
     } // init
 
   private:
@@ -86,28 +130,68 @@
     // count of weights
     double _weightSum;
     /// Histograms
-    AIDA::IHistogram1D* _histdSigDp;
-    AIDA::IHistogram1D* _histMult;
+    AIDA::IHistogram1D* _hist_all_Jpsi;
+    AIDA::IHistogram1D* _hist_chi_c1;
+    AIDA::IHistogram1D* _hist_chi_c2;
+    AIDA::IHistogram1D* _hist_Psi_prime;
+    AIDA::IHistogram1D* _hist_primary_Jpsi;
+
+    AIDA::IHistogram1D* _mult_JPsi;
+    AIDA::IHistogram1D* _mult_JPsi_direct;
+    AIDA::IHistogram1D* _mult_chi_c1;
+    AIDA::IHistogram1D* _mult_chi_c1_direct;
+    AIDA::IHistogram1D* _mult_chi_c2;
+    AIDA::IHistogram1D* _mult_chi_c2_direct;
+    AIDA::IHistogram1D* _mult_Psi2S;
     //@}
 
     void findDecayProducts(const GenParticle & p,
-			   vector<GenParticle *> & pions) {
+			   vector<GenParticle *> & allJpsi,
+			   vector<GenParticle *> & primaryJpsi,
+			   vector<GenParticle *> & Psiprime,
+			   vector<GenParticle *> & all_chi_c1,
+			   vector<GenParticle *> & all_chi_c2,
+			   vector<GenParticle *> & primary_chi_c1,
+			   vector<GenParticle *> & primary_chi_c2) {
       const GenVertex* dv = p.end_vertex();
+      bool isOnium(false);
+      for (GenVertex::particles_in_const_iterator pp = dv->particles_in_const_begin() ;
+	   pp != dv->particles_in_const_end() ; ++pp) {
+	int id = (*pp)->pdg_id();
+	id = id%1000;
+	id -= id%10;
+	id /= 10;
+	if(id==44) isOnium = true;
+      }
       for (GenVertex::particles_out_const_iterator
 	     pp = dv->particles_out_const_begin();
 	   pp != dv->particles_out_const_end(); ++pp) {
 	int id = (*pp)->pdg_id();
-	if(id==111) {
-	  pions.push_back(*pp);
+	if(id==100443) {
+	  Psiprime.push_back(*pp);
+	}
+	else if(id==20443) {
+	  all_chi_c1.push_back(*pp);
+	  if(!isOnium) primary_chi_c1.push_back(*pp);
+	}
+	else if(id==445) {
+	  all_chi_c2.push_back(*pp);
+	  if(!isOnium) primary_chi_c2.push_back(*pp);
+	}
+	else if(id==443) {
+	  allJpsi.push_back(*pp);
+	  if(!isOnium) primaryJpsi.push_back(*pp);
+	}
+	if((*pp)->end_vertex()) {
+	  findDecayProducts(**pp,allJpsi,primaryJpsi,Psiprime,
+			    all_chi_c1,all_chi_c2,primary_chi_c1,primary_chi_c2);
 	}
-	else if((*pp)->end_vertex())
-	  findDecayProducts(**pp,pions);
       }
     }
   };
 
 
   // The hook for the plugin system
-  DECLARE_RIVET_PLUGIN(BELLE_2001_S4598261);
+  DECLARE_RIVET_PLUGIN(BABAR_2003_I593379);
 
 }

Modified: branches/2012-06-aidarivet/src/Analyses/Makefile.am
==============================================================================
--- branches/2012-06-aidarivet/src/Analyses/Makefile.am	Fri Feb  8 09:49:59 2013	(r4133)
+++ branches/2012-06-aidarivet/src/Analyses/Makefile.am	Sat Feb  9 11:22:51 2013	(r4134)
@@ -332,6 +332,7 @@
     ARGUS_1993_S2653028.cc \
     ARGUS_1993_S2669951.cc \
     ARGUS_1993_S2789213.cc \
+    BABAR_2003_I593379.cc \
     BABAR_2005_S6181155.cc \
     BABAR_2007_S6895344.cc \
     BABAR_2007_S7266081.cc \


More information about the Rivet-svn mailing list