Rivet analyses referenceCLEOC_2011_I913909Dalitz plot analysis of $D^0\to K_S^0\pi^0\pi^0$Experiment: CLEOC (CESR) Inspire ID: 913909 Status: VALIDATED NOHEPDATA Authors:
Beam energies: ANY Run details:
Measurement of the mass distributions in the decay $D^0\to K_S^0\pi^0\pi^0$. The data were read from the plots in the paper, with the background given subtracted. Resolution/acceptance effects have been not unfolded and given the agreement with the model in the paper this analysis should only be used for qualitative studies. Source code: CLEOC_2011_I913909.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/UnstableParticles.hh"
4#include "Rivet/Projections/DecayedParticles.hh"
5
6namespace Rivet {
7
8
9 /// @brief D0 -> K0 pi0pi0
10 class CLEOC_2011_I913909 : public Analysis {
11 public:
12
13 /// Constructor
14 RIVET_DEFAULT_ANALYSIS_CTOR(CLEOC_2011_I913909);
15
16
17 /// @name Analysis methods
18 /// @{
19
20 /// Book histograms and initialise projections before the run
21 void init() {
22 // Initialise and register projections
23 UnstableParticles ufs = UnstableParticles(Cuts::abspid==421);
24 declare(ufs, "UFS");
25 DecayedParticles D0(ufs);
26 D0.addStable(PID::PI0);
27 D0.addStable(PID::K0S);
28 D0.addStable(PID::ETA);
29 D0.addStable(PID::ETAPRIME);
30 declare(D0, "D0");
31 // histograms
32 book(_h_Kpi ,1,1,1);
33 book(_h_pipi,1,1,2);
34 book(_dalitz, "dalitz",50,0.,1.9,50,0.3,3.1);
35 }
36
37 /// Perform the per-event analysis
38 void analyze(const Event& event) {
39 static const map<PdgId,unsigned int> & mode = { {111,2}, { 310,1}};
40 DecayedParticles D0 = apply<DecayedParticles>(event, "D0");
41 // loop over particles
42 for(unsigned int ix=0;ix<D0.decaying().size();++ix) {
43 if( ! D0.modeMatches(ix,3,mode ) ) continue;
44 const Particles & K0 = D0.decayProducts()[ix].at(310);
45 const Particles & pi0= D0.decayProducts()[ix].at(111);
46 double mpipi = (pi0[0].momentum()+pi0[1].momentum()).mass2();
47 _h_pipi->fill(mpipi);
48 for(unsigned int ix=0;ix<2;++ix) {
49 double mKpi = (K0[0].momentum()+pi0[ix].momentum()).mass2();
50 _h_Kpi->fill(mKpi);
51 _dalitz->fill(mpipi,mKpi);
52 }
53 }
54 }
55
56
57 /// Normalise histograms etc., after the run
58 void finalize() {
59 normalize(_h_pipi);
60 normalize(_h_Kpi);
61 normalize(_dalitz);
62 }
63
64 /// @}
65
66
67 /// @name Histograms
68 /// @{
69 Histo1DPtr _h_Kpi, _h_pipi;
70 Histo2DPtr _dalitz;
71 /// @}
72
73
74 };
75
76
77 RIVET_DECLARE_PLUGIN(CLEOC_2011_I913909);
78
79}
|
