1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/UnstableParticles.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief e+e- -> 3(pi+pi-)
10  class CMD3_2013_I1217420 : public Analysis {
11  public:
12
13    /// Constructor
14    RIVET_DEFAULT_ANALYSIS_CTOR(CMD3_2013_I1217420);
15
16
17    /// @name Analysis methods
18    /// @{
19
20    /// Book histograms and initialise projections before the run
21    void init() {
22
23      // Initialise and register projections
24      declare(FinalState(), "FS");
25      declare(UnstableParticles(), "UFS");
26
27      // Book histograms
28
29      book(_c_all, "/TMP/all");
30
31    }
32
33    void findChildren(const Particle & p,map<long,int> & nRes, int &ncount) {
34      for(const Particle &child : p.children()) {
35	if(child.children().empty()) {
36	  --nRes[child.pid()];
37	  --ncount;
38	}
39	else
40	  findChildren(child,nRes,ncount);
41      }
42    }
43
44    /// Perform the per-event analysis
45    void analyze(const Event& event) {
46      // find the final-state particles
47      const FinalState& fs = apply<FinalState>(event, "FS");
48      map<long,int> nCount;
49      int ntotal(0);
50      for (const Particle& p : fs.particles()) {
51	nCount[p.pid()] += 1;
52	++ntotal;
53      }
54      if(ntotal==6 && nCount[211]==3 && nCount[-211]==3) {
55	_c_all->fill();
56      }
57    }
58
59
60    /// Normalise histograms etc., after the run
61    void finalize() {
62
63      double fact = crossSection()/nanobarn/sumOfWeights();
64      double sigma = _c_all->val()*fact;
65      double error = _c_all->err()*fact;
66      for (unsigned int ihist=1;ihist<4;++ihist) {
67        Estimate1DPtr  mult;
68        book(mult, 1, 1, ihist);
69        for (auto& b : mult->bins()) {
70          if (inRange(sqrtS()/MeV, b.xMin(), b.xMax())) {
71            b.set(sigma, error);
72          }
73        }
74      }
75    }
76
77    /// @}
78
79
80    /// @name Histograms
81    /// @{
82    CounterPtr _c_all;
83    /// @}
84
85
86  };
87
88
89  RIVET_DECLARE_PLUGIN(CMD3_2013_I1217420);
90
91
92}