Rivet analyses referenceBABAR_2012_I892684Cross sections for $e^+e^-\to$ $K^+K^-\pi^+\pi^-$, $K^+K^-\pi^0\pi^0$, and $K^+K^-K^+K^-$ between 1.28 and 5.0 GeVExperiment: BABAR (PEP-II) Inspire ID: 892684 Status: VALIDATED Authors:
Beam energies: (0.6, 0.6); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.7, 0.7); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.8, 0.8); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (0.9, 0.9); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.0, 1.0); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.1, 1.1); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.2, 1.2); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.3, 1.3); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.4, 1.4); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.5, 1.5); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.6, 1.6); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.7, 1.7); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.8, 1.8); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (1.9, 1.9); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.0, 2.0); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.1, 2.1); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.2, 2.2); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.3, 2.3); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.4, 2.4); (2.5, 2.5); (2.5, 2.5); (2.5, 2.5); (2.5, 2.5) GeV Run details:
Measurement of the cross section for $e^+e^- \to$ $K^+K^-\pi^+\pi^-$, $K^+K^-\pi^0\pi^0$ and $K^+K^-K^+K^-$ via radiative return, including the identification of $K^{*0}$, $\phi$ and $f_0(980)$ mesons for energies between 1.28 and 5.0 GeV. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples. Source code: BABAR_2012_I892684.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/FinalState.hh"
4#include "Rivet/Projections/UnstableParticles.hh"
5
6
7namespace Rivet {
8
9
10 /// @brief e+e- > K+K- + pi+pi-, pi0pi0 or K+K-
11 class BABAR_2012_I892684 : public Analysis {
12 public:
13
14 /// Constructor
15 RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2012_I892684);
16
17
18 /// @name Analysis methods
19 /// @{
20
21 /// Book histograms and initialise projections before the run
22 void init() {
23
24 // Initialise and register projections
25 declare(FinalState(), "FS");
26 declare(UnstableParticles(), "UFS");
27 bool matched = false;
28 for(unsigned int ix=0;ix<8;++ix) {
29 book(_sigma[ix],1+ix,1,1);
30 for (const string& en : _sigma[ix].binning().edges<0>()) {
31 const double end = std::stod(en)*GeV;
32 if (isCompatibleWithSqrtS(end)) {
33 _ecms[ix] = en;
34 matched=true;
35 break;
36 }
37 }
38 }
39 if(!matched)
40 MSG_ERROR("Beam energy incompatible with analysis.");
41 }
42
43 void findChildren(const Particle & p,map<long,int> & nRes, int &ncount) {
44 for (const Particle &child : p.children()) {
45 if(child.children().empty()) {
46 nRes[child.pid()]-=1;
47 --ncount;
48 }
49 else
50 findChildren(child,nRes,ncount);
51 }
52 }
53
54
55 /// Perform the per-event analysis
56 void analyze(const Event& event) {
57 const FinalState& fs = apply<FinalState>(event, "FS");
58
59 map<long,int> nCount;
60 int ntotal(0);
61 for (const Particle& p : fs.particles()) {
62 nCount[p.pid()] += 1;
63 ++ntotal;
64 }
65 const FinalState& ufs = apply<FinalState>(event, "UFS");
66 for (const Particle& p : ufs.particles()) {
67 if(p.children().empty()) continue;
68 // K*0
69 if(abs(p.pid())==313) {
70 map<long,int> nRes=nCount;
71 int ncount = ntotal;
72 findChildren(p,nRes,ncount);
73 // K* K+/- pi-/+
74 if(ncount !=2 ) continue;
75 bool matched = true;
76 for(auto const & val : nRes) {
77 if(abs(val.first)==321 || abs(val.first)==211) {
78 continue;
79 }
80 else if(val.second!=0) {
81 matched = false;
82 break;
83 }
84 }
85 if(matched==false) continue;
86 if((nCount[321] == 1 && nCount[-321] ==0 &&
87 nCount[211] == 0 && nCount[-211] == 1) ||
88 (nCount[321] == 0 && nCount[-321] ==1 &&
89 nCount[211] == 1 && nCount[-211] == 0)) {
90 if(!_ecms[1].empty()) _sigma[1]->fill(_ecms[1]);
91 }
92 }
93 else if(p.pid()==333) {
94 map<long,int> nRes=nCount;
95 int ncount = ntotal;
96 findChildren(p,nRes,ncount);
97 // phi pi+pi-
98 if(ncount==2) {
99 bool matched = true;
100 for(auto const & val : nRes) {
101 if(abs(val.first)==211) {
102 if(val.second!=1) {
103 matched = false;
104 break;
105 }
106 }
107 else if(val.second!=0) {
108 matched = false;
109 break;
110 }
111 }
112 if(matched) {
113 if(!_ecms[2].empty()) _sigma[2]->fill(_ecms[2]);
114 }
115 }
116 for (const Particle& p2 : ufs.particles()) {
117 if(p2.pid()!=9010221&&p2.pid()!=9000221) continue;
118 if(p2.parents()[0].isSame(p)) continue;
119 map<long,int> nResB = nRes;
120 int ncountB = ncount;
121 findChildren(p2,nResB,ncountB);
122 if(ncountB!=0) continue;
123 bool matched2 = true;
124 for(auto const & val : nResB) {
125 if(val.second!=0) {
126 matched2 = false;
127 break;
128 }
129 }
130 if(matched2) {
131 if(p2.pid()==9010221) {
132 if(!_ecms[6].empty()) _sigma[6]->fill(_ecms[6]);
133 if(!_ecms[3].empty()) _sigma[3]->fill(_ecms[3]);
134 }
135 else {
136 if(!_ecms[4].empty()) _sigma[4]->fill(_ecms[4]);
137 }
138 }
139 }
140 }
141 }
142 if(ntotal==4) {
143 if(nCount[321]==1 && nCount[-321]==1 && nCount[211]==1 && nCount[-211]==1) {
144 if(!_ecms[0].empty()) _sigma[0]->fill(_ecms[0]);
145 }
146 else if( nCount[321]==1 && nCount[-321]==1 && nCount[111]==2) {
147 if(!_ecms[5].empty()) _sigma[5]->fill(_ecms[5]);
148 }
149 else if( nCount[321]==2 && nCount[-321]==2) {
150 if(!_ecms[7].empty()) _sigma[7]->fill(_ecms[7]);
151 }
152 }
153 }
154
155
156 /// Normalise histograms etc., after the run
157 void finalize() {
158 double fact = crossSection()/ sumOfWeights() /nanobarn;
159 for(unsigned int ix=0;ix<8;++ix)
160 scale(_sigma[ix],fact);
161 }
162
163 /// @}
164
165
166 /// @name Histograms
167 /// @{
168 BinnedHistoPtr<string> _sigma[8];
169 string _ecms[8];
170 /// @}
171
172
173 };
174
175
176 RIVET_DECLARE_PLUGIN(BABAR_2012_I892684);
177
178
179}
|