Rivet analyses referenceBABAR_2013_I1206605Dalitz plot analysis of $D^+\to K^+K^-\pi^+$Experiment: BABAR (PEP-II) Inspire ID: 1206605 Status: VALIDATED NOHEPDATA Authors:
Beam energies: ANY Run details:
Kinematic distributions in the decay $D^+\to \to K^+K^-\pi^+$. Resolution/acceptance effects have been not unfolded. Given the agreement with the model in the paper this analysis should only be used for qualitative studies. Source code: BABAR_2013_I1206605.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 D+ -> K+K-pi+
10 class BABAR_2013_I1206605 : public Analysis {
11 public:
12
13 /// Constructor
14 RIVET_DEFAULT_ANALYSIS_CTOR(BABAR_2013_I1206605);
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==411);
24 declare(ufs, "UFS");
25 DecayedParticles DP(ufs);
26 DP.addStable(PID::PI0);
27 DP.addStable(PID::K0S);
28 DP.addStable(PID::ETA);
29 DP.addStable(PID::ETAPRIME);
30 declare(DP, "DP");
31 // histos
32 book(_h_Kppi,1,1,1);
33 book(_h_KK ,1,1,2);
34 book(_h_Kmpi,1,1,3);
35 book(_dalitz, "dalitz",50,0.,2.,50,0.9,3.1);
36 }
37
38
39 /// Perform the per-event analysis
40 void analyze(const Event& event) {
41 static const map<PdgId,unsigned int> & mode = { { 321,1},{-321,1}, { 211,1}};
42 static const map<PdgId,unsigned int> & modeCC = { { 321,1},{-321,1}, {-211,1}};
43 DecayedParticles DP = apply<DecayedParticles>(event, "DP");
44 // loop over particles
45 for(unsigned int ix=0;ix<DP.decaying().size();++ix) {
46 int sign = 1;
47 if (DP.decaying()[ix].pid()>0 && DP.modeMatches(ix,3,mode)) {
48 sign=1;
49 }
50 else if (DP.decaying()[ix].pid()<0 && DP.modeMatches(ix,3,modeCC)) {
51 sign=-1;
52 }
53 else
54 continue;
55 const Particle & Kp = DP.decayProducts()[ix].at( sign*321)[0];
56 const Particle & Km = DP.decayProducts()[ix].at(-sign*321)[0];
57 const Particle & pip= DP.decayProducts()[ix].at( sign*211)[0];
58 double mminus = (Km.momentum()+pip.momentum() ).mass2();
59 double mplus = (Kp.momentum()+pip.momentum() ).mass2();
60 double mKK = (Kp.momentum()+Km.momentum()).mass2();
61 _h_Kppi->fill(mplus);
62 _h_Kmpi->fill(mminus);
63 _h_KK ->fill(mKK);
64 _dalitz->fill(mminus,mKK);
65 }
66 }
67
68
69 /// Normalise histograms etc., after the run
70 void finalize() {
71 normalize(_h_Kppi);
72 normalize(_h_Kmpi);
73 normalize(_h_KK );
74 normalize(_dalitz);
75 }
76
77 /// @}
78
79
80 /// @name Histograms
81 /// @{
82 Histo1DPtr _h_Kppi,_h_Kmpi,_h_KK;
83 Histo2DPtr _dalitz;
84 /// @}
85
86
87 };
88
89
90 RIVET_DECLARE_PLUGIN(BABAR_2013_I1206605);
91
92}
|