|
[Rivet-svn] r4134 - in branches/2012-06-aidarivet: data/anainfo data/plotinfo data/refdata src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgSat 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 |