Rivet analyses
MC analysis of Dalitz plots in three-body ηc decays
Experiment: ()
Status: VALIDATED NOHEPDATA
Authors: - Peter Richardson
References: none listed
Beams: * *
Beam energies: ANY
Run details: - Any process producing eta_c mesons
Monte Carlo analysis of common three-body decays of etac mesons. The π+π−η, π0π0η, π+π−η′, π0π0η′, K+K−η, KS, L0KS, L0η, K+K−η′, KS, L0KS, L0η′, K+K−π0, KS, L0KS, L0π0, KS, L0K+π−, KS, L0K−π+ modes are included.
Source
code:MC_DALITZ_ETAC.cc
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/UnstableParticles.hh"
#include "Rivet/Projections/DecayedParticles.hh"
namespace Rivet {
/// @brief eta_c Dalitz plots
class MC_DALITZ_ETAC : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(MC_DALITZ_ETAC);
/// @name Analysis methods
/// @{
/// Book histograms and initialise projections before the run
void init() {
// Initialise and register projections
UnstableParticles ufs = UnstableParticles(Cuts::pid== 441);
declare(ufs, "UFS");
DecayedParticles ETAC(ufs);
ETAC.addStable(PID::PI0);
ETAC.addStable(PID::K0S);
ETAC.addStable(PID::ETA);
ETAC.addStable(PID::ETAPRIME);
declare(ETAC,"ETAC");
// histos
// pi+pi- eta
book(_h1_pippim,"h01_pippim",100,0.2,2.5);
book(_h1_pipeta,"h01_pipeta",100,0.5,3.0);
book(_h1_pimeta,"h01_pimeta",100,0.5,3.0);
book(_dalitz[0], "dalitz01",50,0.,9.,50,0.,9.);
// pi0 pi0 eta
book(_h2_pi0pi0,"h02_pi0pi0",100,0.2,2.5);
book(_h2_pi0eta,"h02_pi0eta",100,0.5,3.0);
book(_dalitz[1], "dalitz02",50,0.,9.,50,0.,9.);
// pi+pi- eta'
book(_h3_pippim,"h03_pippim",100,0.2,2.2);
book(_h3_pipeta,"h03_pipeta",100,1.0,3.0);
book(_h3_pimeta,"h03_pimeta",100,1.0,3.0);
book(_dalitz[2], "dalitz03",50,0.,9.,50,0.,9.);
// pi0 pi0 eta'
book(_h4_pi0pi0,"h04_pi0pi0",100,0.2,2.2);
book(_h4_pi0eta,"h04_pi0eta",100,1.0,3.0);
book(_dalitz[3], "dalitz04",50,0.,9.,50,0.,9.);
// K+K- eta
book(_h5_KpKm ,"h05_KpKm" ,100,0.5,2.5);
book(_h5_Kpeta,"h05_Kpeta",100,1.0,3.0);
book(_h5_Kmeta,"h05_Kmeta",100,1.0,3.0);
book(_dalitz[4], "dalitz05",50,1.,7.,50,1.,7.);
// KS0 KS0 eta
book(_h6_KS0KS0,"h06_KS0KS0",100,0.5,2.5);
book(_h6_KS0eta,"h06_KS0eta",100,1.0,3.0);
book(_dalitz[5], "dalitz06",50,1.,7.,50,1.,7.);
// KL0 KL0 eta
book(_h7_KL0KL0,"h07_KL0KL0",100,0.5,2.5);
book(_h7_KL0eta,"h07_KL0eta",100,1.0,3.0);
book(_dalitz[6], "dalitz07",50,1.,7.,50,1.,7.);
// K+K- eta'
book(_h8_KpKm ,"h08_KpKm" ,100,0.9,2.2);
book(_h8_Kpeta,"h08_Kpeta",100,1.3,2.7);
book(_h8_Kmeta,"h08_Kmeta",100,1.3,2.7);
book(_dalitz[7], "dalitz08",50,1.5,6.5,50,1.5,6.5);
// KS0 KS0 eta'
book(_h9_KS0KS0,"h09_KS0KS0",100,0.9,2.2);
book(_h9_KS0eta,"h09_KS0eta",100,1.3,2.7);
book(_dalitz[8], "dalitz09",50,1.5,6.5,50,1.5,6.5);
// KL0 KL0 eta'
book(_h10_KL0KL0,"h10_KL0KL0",100,0.9,2.2);
book(_h10_KL0eta,"h10_KL0eta",100,1.3,2.7);
book(_dalitz[9] , "dalitz10",50,1.5,6.5,50,1.5,6.5);
// K+K- pi0
book(_h11_KpKm ,"h11_KpKm" ,100,0.9,2.8);
book(_h11_Kppi0,"h11_Kppi0",100,0.6,2.8);
book(_h11_Kmpi0,"h11_Kmpi0",100,0.6,2.8);
book(_dalitz[10],"dalitz11",50,0.3,6.5,50,0.3,6.5);
// KS0 KS0 pi0
book(_h12_KS0KS0,"h12_KS0KS0",100,0.9,2.8);
book(_h12_KS0pi0,"h12_KS0pi0",100,0.6,2.8);
book(_dalitz[11] ,"dalitz12",50,0.3,6.5,50,0.3,6.5);
// KL0 KL0 pi0
book(_h13_KL0KL0,"h13_KL0KL0",100,0.9,2.8);
book(_h13_KL0pi0,"h13_KL0pi0",100,0.6,2.8);
book(_dalitz[12] , "dalitz13",50,0.3,6.5,50,0.3,6.5);
// KS0 K+ pi-
book(_h14_KpKS0 ,"h14_KpKS0" ,100,0.9,2.8);
book(_h14_Kppim ,"h14_Kppim" ,100,0.6,2.8);
book(_h14_KS0pim,"h14_KS0pim",100,0.6,2.8);
book(_dalitz[13] ,"dalitz14",50,0.3,6.5,50,0.3,6.5);
// KS0 K- pi+
book(_h15_KmKS0 ,"h15_KmKS0" ,100,0.9,2.8);
book(_h15_Kmpip ,"h15_Kmpip" ,100,0.6,2.8);
book(_h15_KS0pip,"h15_KS0pip",100,0.6,2.8);
book(_dalitz[14] , "dalitz15",50,0.3,6.5,50,0.3,6.5);
// KL0 K+ pi-
book(_h16_KpKL0 ,"h16_KpKL0" ,100,0.9,2.8);
book(_h16_Kppim ,"h16_Kppim" ,100,0.6,2.8);
book(_h16_KL0pim,"h16_KL0pim",100,0.6,2.8);
book(_dalitz[15] ,"dalitz16",50,0.3,6.5,50,0.3,6.5);
// KL0 K- pi+
book(_h17_KmKL0 ,"h17_KmKL0" ,100,0.9,2.8);
book(_h17_Kmpip ,"h17_Kmpip" ,100,0.6,2.8);
book(_h17_KL0pip,"h17_KL0pip",100,0.6,2.8);
book(_dalitz[16] ,"dalitz17",50,0.3,6.5,50,0.3,6.5);
}
/// Perform the per-event analysis
void analyze(const Event& event) {
static const map<PdgId,unsigned int> & mode1 = { { 211,1}, {-211,1}, { 221,1}};
static const map<PdgId,unsigned int> & mode2 = { { 111,2}, { 221,1}};
static const map<PdgId,unsigned int> & mode3 = { { 211,1}, {-211,1}, { 331,1}};
static const map<PdgId,unsigned int> & mode4 = { { 111,2}, { 331,1}};
static const map<PdgId,unsigned int> & mode5 = { { 321,1}, {-321,1}, { 221,1}};
static const map<PdgId,unsigned int> & mode6 = { { 310,2}, { 221,1}};
static const map<PdgId,unsigned int> & mode7 = { { 130,2}, { 221,1}};
static const map<PdgId,unsigned int> & mode8 = { { 321,1}, {-321,1}, { 331,1}};
static const map<PdgId,unsigned int> & mode9 = { { 310,2}, { 331,1}};
static const map<PdgId,unsigned int> & mode10 = { { 130,2}, { 331,1}};
static const map<PdgId,unsigned int> & mode11 = { { 321,1}, {-321,1}, { 111,1}};
static const map<PdgId,unsigned int> & mode12 = { { 310,2}, { 111,1}};
static const map<PdgId,unsigned int> & mode13 = { { 130,2}, { 111,1}};
static const map<PdgId,unsigned int> & mode14 = { { 321,1}, {-211,1}, { 310,1}};
static const map<PdgId,unsigned int> & mode15 = { {-321,1}, { 211,1}, { 310,1}};
static const map<PdgId,unsigned int> & mode16 = { { 321,1}, {-211,1}, { 130,1}};
static const map<PdgId,unsigned int> & mode17 = { {-321,1}, { 211,1}, { 130,1}};
DecayedParticles ETAC = apply<DecayedParticles>(event, "ETAC");
// loop over particles
for(unsigned int ix=0;ix<ETAC.decaying().size();++ix) {
if(ETAC.decaying()[ix].mass()<2.93 || ETAC.decaying()[ix].mass()>3.03) continue;
// pi+ pi- eta
if (ETAC.modeMatches(ix,3,mode1)) {
const Particle & pip = ETAC.decayProducts()[ix].at( 211)[0];
const Particle & pim = ETAC.decayProducts()[ix].at(-211)[0];
const Particle & eta = ETAC.decayProducts()[ix].at( 221)[0];
double mplus = (pip.momentum()+eta.momentum()).mass2();
double mminus = (pim.momentum()+eta.momentum()).mass2();
double mpipi = (pip.momentum()+pim.momentum()).mass2();
_h1_pippim->fill(sqrt(mpipi));
_h1_pipeta->fill(sqrt(mplus));
_h1_pimeta->fill(sqrt(mminus));
_dalitz[0]->fill(mplus,mminus);
}
// pi0 pi0 eta
else if (ETAC.modeMatches(ix,3,mode2)) {
const Particles & pi0 = ETAC.decayProducts()[ix].at( 111);
const Particle & eta = ETAC.decayProducts()[ix].at( 221)[0];
double m1 = (pi0[0].momentum()+eta.momentum()).mass2();
double m2 = (pi0[1].momentum()+eta.momentum()).mass2();
double mpipi = (pi0[0].momentum()+pi0[1].momentum()).mass2();
_h2_pi0pi0->fill(sqrt(mpipi));
_h2_pi0eta->fill(sqrt(m1));
_h2_pi0eta->fill(sqrt(m2));
_dalitz[1]->fill(m1,m2);
_dalitz[1]->fill(m2,m1);
}
// pi+ pi- eta'
else if (ETAC.modeMatches(ix,3,mode3)) {
const Particle & pip = ETAC.decayProducts()[ix].at( 211)[0];
const Particle & pim = ETAC.decayProducts()[ix].at(-211)[0];
const Particle & etaP = ETAC.decayProducts()[ix].at( 331)[0];
double mplus = (pip.momentum()+etaP.momentum()).mass2();
double mminus = (pim.momentum()+etaP.momentum()).mass2();
double mpipi = (pip.momentum()+pim.momentum()).mass2();
_h3_pippim->fill(sqrt(mpipi));
_h3_pipeta->fill(sqrt(mplus));
_h3_pimeta->fill(sqrt(mminus));
_dalitz[2]->fill(mplus,mminus);
}
// pi0 pi0 eta'
else if (ETAC.modeMatches(ix,3,mode4)) {
const Particles & pi0 = ETAC.decayProducts()[ix].at( 111);
const Particle & etaP = ETAC.decayProducts()[ix].at( 331)[0];
double m1 = (pi0[0].momentum()+etaP.momentum()).mass2();
double m2 = (pi0[1].momentum()+etaP.momentum()).mass2();
double mpipi = (pi0[0].momentum()+pi0[1].momentum()).mass2();
_h4_pi0pi0->fill(sqrt(mpipi));
_h4_pi0eta->fill(sqrt(m1));
_h4_pi0eta->fill(sqrt(m2));
_dalitz[3]->fill(m1,m2);
_dalitz[3]->fill(m2,m1);
}
// K+K- eta
else if (ETAC.modeMatches(ix,3,mode5)) {
const Particle & Kp = ETAC.decayProducts()[ix].at( 321)[0];
const Particle & Km = ETAC.decayProducts()[ix].at(-321)[0];
const Particle & eta = ETAC.decayProducts()[ix].at( 221)[0];
double mplus = (Kp.momentum()+eta.momentum()).mass2();
double mminus = (Km.momentum()+eta.momentum()).mass2();
double mKK = (Kp.momentum()+Km .momentum()).mass2();
_h5_KpKm->fill(sqrt(mKK));
_h5_Kpeta->fill(sqrt(mplus));
_h5_Kmeta->fill(sqrt(mminus));
_dalitz[4]->fill(mplus,mminus);
}
// KS0 KS0 eta
else if (ETAC.modeMatches(ix,3,mode6)) {
const Particles & KS0 = ETAC.decayProducts()[ix].at( 310);
const Particle & eta = ETAC.decayProducts()[ix].at( 221)[0];
double m1 = (KS0[0].momentum()+eta.momentum()).mass2();
double m2 = (KS0[1].momentum()+eta.momentum()).mass2();
double mKSKS = (KS0[0].momentum()+KS0[1].momentum()).mass2();
_h6_KS0KS0->fill(sqrt(mKSKS));
_h6_KS0eta->fill(sqrt(m1));
_h6_KS0eta->fill(sqrt(m2));
_dalitz[5]->fill(m1,m2);
_dalitz[5]->fill(m2,m1);
}
// KL0 KL0 eta
else if (ETAC.modeMatches(ix,3,mode7)) {
const Particles & KL0 = ETAC.decayProducts()[ix].at( 130);
const Particle & eta = ETAC.decayProducts()[ix].at( 221)[0];
double m1 = (KL0[0].momentum()+eta.momentum()).mass2();
double m2 = (KL0[1].momentum()+eta.momentum()).mass2();
double mKLKL = (KL0[0].momentum()+KL0[1].momentum()).mass2();
_h7_KL0KL0->fill(sqrt(mKLKL));
_h7_KL0eta->fill(sqrt(m1));
_h7_KL0eta->fill(sqrt(m2));
_dalitz[6]->fill(m1,m2);
_dalitz[6]->fill(m2,m1);
}
// K+K- eta'
else if (ETAC.modeMatches(ix,3,mode8)) {
const Particle & Kp = ETAC.decayProducts()[ix].at( 321)[0];
const Particle & Km = ETAC.decayProducts()[ix].at(-321)[0];
const Particle & etaP = ETAC.decayProducts()[ix].at( 331)[0];
double mplus = (Kp.momentum()+etaP.momentum()).mass2();
double mminus = (Km.momentum()+etaP.momentum()).mass2();
double mKK = (Kp.momentum()+Km .momentum()).mass2();
_h8_KpKm->fill(sqrt(mKK));
_h8_Kpeta->fill(sqrt(mplus));
_h8_Kmeta->fill(sqrt(mminus));
_dalitz[7]->fill(mplus,mminus);
}
// KS0 KS0 eta'
else if (ETAC.modeMatches(ix,3,mode9)) {
const Particles & KS0 = ETAC.decayProducts()[ix].at( 310);
const Particle & etaP = ETAC.decayProducts()[ix].at( 331)[0];
double m1 = (KS0[0].momentum()+etaP.momentum()).mass2();
double m2 = (KS0[1].momentum()+etaP.momentum()).mass2();
double mKSKS = (KS0[0].momentum()+KS0[1].momentum()).mass2();
_h9_KS0KS0->fill(sqrt(mKSKS));
_h9_KS0eta->fill(sqrt(m1));
_h9_KS0eta->fill(sqrt(m2));
_dalitz[8]->fill(m1,m2);
_dalitz[8]->fill(m2,m1);
}
// KL0 KL0 eta'
else if (ETAC.modeMatches(ix,3,mode10)) {
const Particles & KL0 = ETAC.decayProducts()[ix].at( 130);
const Particle & etaP = ETAC.decayProducts()[ix].at( 331)[0];
double m1 = (KL0[0].momentum()+etaP.momentum()).mass2();
double m2 = (KL0[1].momentum()+etaP.momentum()).mass2();
double mKLKL = (KL0[0].momentum()+KL0[1].momentum()).mass2();
_h10_KL0KL0->fill(sqrt(mKLKL));
_h10_KL0eta->fill(sqrt(m1));
_h10_KL0eta->fill(sqrt(m2));
_dalitz[9]->fill(m1,m2);
_dalitz[9]->fill(m2,m1);
}
// K+ K- pi0
else if (ETAC.modeMatches(ix,3,mode11)) {
const Particle & Kp = ETAC.decayProducts()[ix].at( 321)[0];
const Particle & Km = ETAC.decayProducts()[ix].at(-321)[0];
const Particle & pi0 = ETAC.decayProducts()[ix].at( 111)[0];
double mplus = (Kp.momentum()+pi0.momentum()).mass2();
double mminus = (Km.momentum()+pi0.momentum()).mass2();
double mKK = (Kp.momentum()+Km .momentum()).mass2();
_h11_KpKm->fill(sqrt(mKK));
_h11_Kppi0->fill(sqrt(mplus));
_h11_Kmpi0->fill(sqrt(mminus));
_dalitz[10]->fill(mplus,mminus);
}
// KS0 KS0 pi0
else if (ETAC.modeMatches(ix,3,mode12)) {
const Particles & KS0 = ETAC.decayProducts()[ix].at( 310);
const Particle & pi0 = ETAC.decayProducts()[ix].at( 111)[0];
double mplus = (KS0[0].momentum()+pi0.momentum()).mass2();
double mminus = (KS0[1].momentum()+pi0.momentum()).mass2();
double mKK = (KS0[0].momentum()+KS0[1].momentum()).mass2();
_h12_KS0KS0->fill(sqrt(mKK));
_h12_KS0pi0->fill(sqrt(mplus));
_h12_KS0pi0->fill(sqrt(mminus));
_dalitz[11]->fill(mplus,mminus);
}
// KL0 KL0 pi0
else if (ETAC.modeMatches(ix,3,mode13)) {
const Particles & KL0 = ETAC.decayProducts()[ix].at( 130);
const Particle & pi0 = ETAC.decayProducts()[ix].at( 111)[0];
double mplus = (KL0[0].momentum()+pi0.momentum()).mass2();
double mminus = (KL0[1].momentum()+pi0.momentum()).mass2();
double mKK = (KL0[0].momentum()+KL0[1].momentum()).mass2();
_h13_KL0KL0->fill(sqrt(mKK));
_h13_KL0pi0->fill(sqrt(mplus));
_h13_KL0pi0->fill(sqrt(mminus));
_dalitz[12]->fill(mplus,mminus);
}
// KS0 K+ pi-
else if (ETAC.modeMatches(ix,3,mode14)) {
const Particle & KS0 = ETAC.decayProducts()[ix].at( 310)[0];
const Particle & Kp = ETAC.decayProducts()[ix].at( 321)[0];
const Particle & pim = ETAC.decayProducts()[ix].at(-211)[0];
double mplus = (Kp.momentum()+pim.momentum()).mass2();
double mminus = (KS0.momentum()+pim.momentum()).mass2();
double mKK = (KS0.momentum()+Kp.momentum()).mass2();
_h14_KpKS0->fill(sqrt(mKK));
_h14_Kppim->fill(sqrt(mplus));
_h14_KS0pim->fill(sqrt(mminus));
_dalitz[13]->fill(mplus,mminus);
}
// KS0 K- pi+
else if (ETAC.modeMatches(ix,3,mode15)) {
const Particle & KS0 = ETAC.decayProducts()[ix].at( 310)[0];
const Particle & Km = ETAC.decayProducts()[ix].at(-321)[0];
const Particle & pip = ETAC.decayProducts()[ix].at( 211)[0];
double mplus = (Km.momentum()+pip.momentum()).mass2();
double mminus = (KS0.momentum()+pip.momentum()).mass2();
double mKK = (KS0.momentum()+Km.momentum()).mass2();
_h15_KmKS0->fill(sqrt(mKK));
_h15_Kmpip->fill(sqrt(mplus));
_h15_KS0pip->fill(sqrt(mminus));
_dalitz[14]->fill(mplus,mminus);
}
// KL0 K+ pi-
else if (ETAC.modeMatches(ix,3,mode16)) {
const Particle & KL0 = ETAC.decayProducts()[ix].at( 130)[0];
const Particle & Kp = ETAC.decayProducts()[ix].at( 321)[0];
const Particle & pim = ETAC.decayProducts()[ix].at(-211)[0];
double mplus = (Kp.momentum()+pim.momentum()).mass2();
double mminus = (KL0.momentum()+pim.momentum()).mass2();
double mKK = (KL0.momentum()+Kp.momentum()).mass2();
_h16_KpKL0->fill(sqrt(mKK));
_h16_Kppim->fill(sqrt(mplus));
_h16_KL0pim->fill(sqrt(mminus));
_dalitz[15]->fill(mplus,mminus);
}
// KL0 K- pip
else if (ETAC.modeMatches(ix,3,mode17)) {
const Particle & KL0 = ETAC.decayProducts()[ix].at( 130)[0];
const Particle & Km = ETAC.decayProducts()[ix].at(-321)[0];
const Particle & pip = ETAC.decayProducts()[ix].at( 211)[0];
double mplus = (Km.momentum()+pip.momentum()).mass2();
double mminus = (KL0.momentum()+pip.momentum()).mass2();
double mKK = (KL0.momentum()+Km.momentum()).mass2();
_h17_KmKL0->fill(sqrt(mKK));
_h17_Kmpip->fill(sqrt(mplus));
_h17_KL0pip->fill(sqrt(mminus));
_dalitz[16]->fill(mplus,mminus);
}
}
}
/// Normalise histograms etc., after the run
void finalize() {
normalize(_h1_pippim);
normalize(_h1_pipeta);
normalize(_h1_pimeta);
normalize(_h2_pi0pi0);
normalize(_h2_pi0eta);
normalize(_h3_pippim);
normalize(_h3_pipeta);
normalize(_h3_pimeta);
normalize(_h4_pi0pi0);
normalize(_h4_pi0eta);
normalize(_h5_KpKm);
normalize(_h5_Kpeta);
normalize(_h5_Kmeta);
normalize(_h6_KS0KS0);
normalize(_h6_KS0eta);
normalize(_h7_KL0KL0);
normalize(_h7_KL0eta);
normalize(_h8_KpKm);
normalize(_h8_Kpeta);
normalize(_h8_Kmeta);
normalize(_h9_KS0KS0);
normalize(_h9_KS0eta);
normalize(_h10_KL0KL0);
normalize(_h10_KL0eta);
normalize(_h11_KpKm );
normalize(_h11_Kppi0);
normalize(_h11_Kmpi0);
normalize(_h12_KS0KS0);
normalize(_h12_KS0pi0);
normalize(_h13_KL0KL0);
normalize(_h13_KL0pi0);
normalize(_h14_KpKS0);
normalize(_h14_Kppim);
normalize(_h14_KS0pim);
normalize(_h15_KmKS0);
normalize(_h15_Kmpip);
normalize(_h15_KS0pip);
normalize(_h16_KpKL0);
normalize(_h16_Kppim);
normalize(_h16_KL0pim);
normalize(_h17_KmKL0);
normalize(_h17_Kmpip);
normalize(_h17_KL0pip);
for(unsigned int ix=0;ix<17;++ix)
normalize(_dalitz[ix]);
}
/// @}
/// @name Histograms
/// @{
Histo1DPtr _h1_pippim,_h1_pipeta,_h1_pimeta;
Histo1DPtr _h2_pi0pi0,_h2_pi0eta;
Histo1DPtr _h3_pippim,_h3_pipeta,_h3_pimeta;
Histo1DPtr _h4_pi0pi0,_h4_pi0eta;
Histo1DPtr _h5_KpKm,_h5_Kpeta,_h5_Kmeta;
Histo1DPtr _h6_KS0KS0,_h6_KS0eta;
Histo1DPtr _h7_KL0KL0,_h7_KL0eta;
Histo1DPtr _h8_KpKm,_h8_Kpeta,_h8_Kmeta;
Histo1DPtr _h9_KS0KS0,_h9_KS0eta;
Histo1DPtr _h10_KL0KL0,_h10_KL0eta;
Histo1DPtr _h11_KpKm,_h11_Kppi0,_h11_Kmpi0;
Histo1DPtr _h12_KS0KS0,_h12_KS0pi0;
Histo1DPtr _h13_KL0KL0,_h13_KL0pi0;
Histo1DPtr _h14_KpKS0,_h14_Kppim,_h14_KS0pim;
Histo1DPtr _h15_KmKS0,_h15_Kmpip,_h15_KS0pip;
Histo1DPtr _h16_KpKL0,_h16_Kppim,_h16_KL0pim;
Histo1DPtr _h17_KmKL0,_h17_Kmpip,_h17_KL0pip;
Histo2DPtr _dalitz[17];
/// @}
};
RIVET_DECLARE_PLUGIN(MC_DALITZ_ETAC);
}