1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/UnstableParticles.hh"
 4
 5namespace Rivet {
 6
 7
 8  /// @brief Spectrum for Omega_c
 9  class BABAR_2007_I746745 : public Analysis {
10  public:
11
12    /// Constructor
13    RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2007_I746745);
14
15
16    /// @name Analysis methods
17    /// @{
18
19    /// Book histograms and initialise projections before the run
20    void init() {
21      declare(UnstableParticles(),"UFS");
22      book(_h_p,1,1,1);
23      book(_b  ,2,1,1);
24      book(_r  ,3,1,1);
25      book(_ups,"/TMP/ups");
26    }
27
28
29    /// Perform the per-event analysis
30    void analyze(const Event& event) {
31      static const int idOmega = 4332;
32      const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
33      bool ups = !ufs.particles(Cuts::pid==300553).empty();
34      if(ups) _ups->fill();
35      for (const Particle& p : ufs.particles(Cuts::abspid==idOmega)) {
36	_h_p->fill(p.momentum().p3().mod());
37	if(p.children().size()==2) {
38	  int sign = p.pid()/p.abspid();
39	  if((p.children()[0].pid()==sign*3334 &&
40	      p.children()[1].pid()==sign*211) ||
41	     (p.children()[1].pid()==sign*3334 &&
42	      p.children()[0].pid()==sign*211) ) {
43	    if(ups) _b->fill();
44	    else    _r->fill("10.58"s);
45	  }
46	}
47      }
48    }
49
50
51    /// Normalise histograms etc., after the run
52    void finalize() {
53      normalize(_h_p);
54      if(_ups->effNumEntries()!=0) {
55	scale(_b, 0.5/ *_ups);
56      }
57      scale(_r,crossSection()/sumOfWeights()/femtobarn);
58    }
59    /// @}
60
61
62    /// @name Histograms
63    /// @{
64    Histo1DPtr _h_p;
65    BinnedHistoPtr<string> _r;
66    CounterPtr _b, _ups;
67    /// @}
68
69
70  };
71
72
73  RIVET_DECLARE_PLUGIN(BABAR_2007_I746745);
74
75}