diff --git a/analyses/pluginMC/MC_D_Dalitz.cc b/analyses/pluginMC/MC_D_Dalitz.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_D_Dalitz.cc @@ -0,0 +1,271 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/FinalState.hh" +#include "Rivet/Projections/FastJets.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class MC_D_Dalitz : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(MC_D_Dalitz); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(UnstableParticles(), "UFS"); + + // Book histograms + _h_plus1 = bookHisto1D("h_plus1" ,200,0.,3. ); + _h_minus1 = bookHisto1D("h_minus1" ,200,0.,3.2 ); + _h_pipi1 = bookHisto1D("h_pipi1" ,200,0.,2. ); + _h_minus2 = bookHisto1D("h_minus2" ,200,0.,3.2 ); + _h_neutral2 = bookHisto1D("h_neutral2",200,0.,3.2 ); + _h_pipi2 = bookHisto1D("h_pipi2" ,200,0.,2. ); + _h_Kpilow3 = bookHisto1D("h_Kpilow3" ,200,0.,2. ); + _h_Kpihigh3 = bookHisto1D("h_Kpihigh3" ,200,0.,3.2 ); + _h_Kpiall3 = bookHisto1D("h_Kpiall3" ,200,0.,3. ); + _h_pipi3 = bookHisto1D("h_pipi3" ,200,0.,2. ); + _h_Kpip4 = bookHisto1D("h_Kpip4" ,200,0.,3.2 ); + _h_pipi4 = bookHisto1D("h_pipi4" ,200,0.,2. ); + _h_Kpi04 = bookHisto1D("h_Kpi04" ,200,0.,3.2); + _h_kppim5 = bookHisto1D("h_kppim5" ,200,0.,3. ); + _h_kppip5 = bookHisto1D("h_kppip5" ,200,0.,3.1 ); + _h_pippim5 = bookHisto1D("h_pippim5" ,200,0.,2. ); + _h_kppim6 = bookHisto1D("h_kppim6" ,200,0.,3.5); + _h_kppip6 = bookHisto1D("h_kppip6" ,200,0.,3.5); + _h_pippim6 = bookHisto1D("h_pippim6" ,200,0.,2.5); + _dalitz1 = bookHisto2D("dalitz1" ,50,0.3,3.2,50,0.3,3.2); + _dalitz2 = bookHisto2D("dalitz2" ,50,0.3,3. ,50,0.3,3. ); + _dalitz3 = bookHisto2D("dalitz3" ,50,0.3,2. ,50,0.07,2. ); + _dalitz4 = bookHisto2D("dalitz4" ,50,0.3,3.1 ,50,0.07,2. ); + _dalitz5 = bookHisto2D("dalitz5" ,50,0.,3. ,50,0.,2. ); + _dalitz6 = bookHisto2D("dalitz6" ,50,0.3,3.5,50,0.07,2.5); + } + + void findDecayProducts(const Particle & mother, unsigned int & nstable, + Particles & pip , Particles & pim , Particles & pi0 , + Particles & Kp , Particles & Km , Particles & K0) { + for(const Particle & p : mother.children()) { + int id = p.pdgId(); + if ( id == PID::KPLUS ) { + Kp.push_back(p); + ++nstable; + } + else if (id == PID::KMINUS ) { + Km.push_back(p); + ++nstable; + } + else if (id == PID::PIPLUS) { + pip.push_back(p); + ++nstable; + } + else if (id == PID::PIMINUS) { + pim.push_back(p); + ++nstable; + } + else if (id == PID::PI0) { + pi0.push_back(p); + ++nstable; + } + else if (id == PID::K0S||id == PID::K0L) { + K0.push_back(p); + ++nstable; + } + else if ( !p.children().empty() ) { + findDecayProducts(p, nstable, pip, pim, pi0, Kp , Km, K0); + } + else + ++nstable; + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + double weight = event.weight(); + for(const Particle& meson : apply(event, "UFS"). + particles(Cuts::abspid== 411 ||Cuts::abspid== 421 ||Cuts::abspid== 431 )) { + unsigned int nstable(0); + Particles pip, pim, pi0, Kp , Km, K0; + findDecayProducts(meson, nstable, pip, pim, pi0, Kp , Km, K0); + if(meson.pdgId()<0) { + swap(pim,pip); + swap(Kp,Km); + } + if(abs(meson.pdgId())==421) { + if(pim.size()==1&&pip.size()==1&&K0.size()==1) { + double mminus = (pim[0].momentum()+K0[0].momentum() ).mass2(); + double mplus = (pip[0].momentum()+K0[0].momentum() ).mass2(); + double mpipi = (pip[0].momentum()+pim[0].momentum()).mass2(); + _h_plus1 ->fill(mplus,weight); + _h_minus1 ->fill(mminus,weight); + _h_pipi1 ->fill(mpipi,weight); + _dalitz1 ->fill(mplus,mminus,weight); + } + else if (pip.size()==1&&Km.size()==1&&pi0.size()==1) { + double mneut = (Km[0].momentum()+pip[0].momentum()).mass2(); + double mminus = (Km[0].momentum()+pi0[0].momentum()).mass2(); + double mpipi = (pip[0].momentum()+pi0[0].momentum()).mass2(); + _h_neutral2 ->fill(mneut,weight); + _h_minus2 ->fill(mminus,weight); + _h_pipi2 ->fill(mpipi,weight); + _dalitz2->fill(mminus,mneut,weight); + } + } + else if(abs(meson.pdgId())==411) { + if(pip.size()==2&&Km.size()==1) { + double mplus = (Km[0].momentum() +pip[0].momentum()).mass2(); + double mminus = (Km[0].momentum() +pip[1].momentum()).mass2(); + double mpipi = (pip[0].momentum()+pip[1].momentum()).mass2(); + if(mplusfill( mminus,weight); + _h_Kpihigh3->fill( mplus,weight); + _h_Kpiall3 ->fill( mminus,weight); + _h_Kpiall3 ->fill( mplus,weight); + _h_pipi3 ->fill( mpipi,weight); + _dalitz3->fill(mminus,mpipi, weight); + } + else if (pip.size()==1&&pi0.size()==1&&K0.size()==1) { + double mminus = (K0[0].momentum()+pip[0].momentum()).mass2(); + double mplus = (K0[0].momentum()+pi0[0].momentum()).mass2(); + double mpipi = (pip[0].momentum()+pi0[0].momentum()).mass2(); + _h_Kpip4 ->fill( mminus, weight); + _h_pipi4 ->fill( mpipi , weight); + _h_Kpi04 ->fill( mplus , weight); + _dalitz4->fill(mplus,mpipi, weight); + } + else if (pim.size()==1&&Kp.size()==1&&pip.size()==1) { + double mplus = (Kp [0].momentum()+pip[0].momentum()).mass2(); + double mminus = (Kp [0].momentum()+pim[0].momentum()).mass2(); + double mpipi = (pip[0].momentum()+pim[0].momentum()).mass2(); + _h_kppim5 ->fill(mminus,weight); + _h_kppip5 ->fill(mplus ,weight); + _h_pippim5->fill(mpipi ,weight); + _dalitz5->fill(mminus,mpipi, weight); + } + } + else if(abs(meson.pdgId())==431) { + if (pim.size()==1&&Kp.size()==1&&pip.size()==1) { + double mplus = (Kp [0].momentum()+pip[0].momentum()).mass2(); + double mminus = (Kp [0].momentum()+pim[0].momentum()).mass2(); + double mpipi = (pip[0].momentum()+pim[0].momentum()).mass2(); + _h_kppim6 ->fill(mminus,weight); + _h_kppip6 ->fill(mplus ,weight); + _h_pippim6->fill(mpipi ,weight); + _dalitz6->fill(mminus,mpipi, weight); + } + } + } + } + + /// Normalise histograms etc., after the run + void finalize() { + normalize(_h_plus1); + normalize(_h_minus1); + normalize(_h_pipi1); + normalize(_dalitz1); + normalize(_h_minus2); + normalize(_h_pipi2); + normalize(_h_neutral2); + normalize(_dalitz2); + normalize(_h_Kpilow3); + normalize(_h_Kpihigh3); + normalize(_h_Kpiall3); + normalize(_h_pipi3); + normalize(_dalitz3); + normalize(_h_Kpip4); + normalize(_h_pipi4); + normalize(_h_Kpi04); + normalize(_dalitz4); + normalize(_h_kppim5); + normalize(_h_kppip5); + normalize(_h_pippim5); + normalize(_dalitz5); + normalize(_h_kppim6); + normalize(_h_kppip6); + normalize(_h_pippim6); + normalize(_dalitz6); + } + //@} + + /// @name Histograms + //@{ + // Histograms for D^0\to \bar{K}^0\pi^+\pi^- + //m^2_+ + Histo1DPtr _h_plus1; + //m^2_+ + Histo1DPtr _h_minus1; + //m^2_{\pi\pi} + Histo1DPtr _h_pipi1; + // Dalitz plot + Histo2DPtr _dalitz1; + + // Histograms for D^0\to K^-\pi^+\pi^0 + // Histogram for the K^-\pi^+ mass + Histo1DPtr _h_minus2; + // Histogram for the \pi^+\pi^0 mass + Histo1DPtr _h_pipi2; + // Histogram for the K^-\pi^0 mass + Histo1DPtr _h_neutral2; + // Dalitz plot + Histo2DPtr _dalitz2; + + // Histograms for D^+\to K^-\pi^+\pi^+ + // Histogram for K^-\pi^+ low + Histo1DPtr _h_Kpilow3; + // Histogram for K^-\pi^+ high + Histo1DPtr _h_Kpihigh3; + // Histogram for K^-\pi^+ all + Histo1DPtr _h_Kpiall3; + // Histogram for \pi^+\pi^- + Histo1DPtr _h_pipi3; + // Dalitz plot + Histo2DPtr _dalitz3; + + // Histograms for D^+\to\bar{K}^0\pi^+\pi^0 + // Histogram for the \bar{K}^0\pi^+ mass + Histo1DPtr _h_Kpip4; + // Histogram for the \pi^+\pi^0 mass + Histo1DPtr _h_pipi4; + // Histogram for the \bar{K}^0\pi^0 mass + Histo1DPtr _h_Kpi04; + // Dalitz plot + Histo2DPtr _dalitz4; + + // Histograms for D^+\to K^+\pi^-\pi^+ + // Histogram for K^+\pi^- + Histo1DPtr _h_kppim5; + // Histogram for K^+\pi^+ + Histo1DPtr _h_kppip5; + // Histogram for \pi^+\pi^- + Histo1DPtr _h_pippim5; + // Dalitz plot + Histo2DPtr _dalitz5; + + // Histograms for D_s^+\to K^+\pi^-\pi^+ + // Histogram for K^+\pi^- + Histo1DPtr _h_kppim6; + // Histogram for K^+\pi^+ + Histo1DPtr _h_kppip6; + // Histogram for \pi^+\pi^- + Histo1DPtr _h_pippim6; + // Dalitz plot + Histo2DPtr _dalitz6; + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(MC_D_Dalitz); + + +} diff --git a/analyses/pluginMC/MC_D_Dalitz.info b/analyses/pluginMC/MC_D_Dalitz.info new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_D_Dalitz.info @@ -0,0 +1,8 @@ +Name: MC_D_Dalitz +Summary: MC analysis of Dalitz plots in three-body $D$-meson decays +Status: VALIDATED +Authors: + - Peter Richardson +RunInfo: Any process producing D mesons +Description: + 'Monte Carlo analysis of $D^0\to \bar{K}^0\pi^+\pi^-$, $D^0\to K^-\pi^+\pi^0$, $D^+\to K^-\pi^+\pi^+$, $D^+\to\bar{K}^0\pi^+\pi^0$, $D^+\to K^+\pi^-\pi^+$ and $D_s^+\to K^+\pi^-\pi^+$ including kinematic distributions and Dalitz plots.' diff --git a/analyses/pluginMC/MC_D_Dalitz.plot b/analyses/pluginMC/MC_D_Dalitz.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_D_Dalitz.plot @@ -0,0 +1,136 @@ +BEGIN PLOT /MC_D_Dalitz/h_plus1 +Title=$\bar{K}^0\pi^+$ mass distribution in $D^0\to \bar{K}^0\pi^+\pi^-$ +XLabel=$m^2_{\bar{K}^0\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{\bar{K}^0\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_minus1 +Title=$\bar{K}^0\pi^-$ mass distribution in $D^0\to \bar{K}^0\pi^+\pi^-$ +XLabel=$m^2_{\bar{K}^0\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{\bar{K}^0\pi^-}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pipi1 +Title=$\pi^+\pi^-$ mass distribution in $D^0\to \bar{K}^0\pi^+\pi^-$ +XLabel=$m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz1 +Title=Dalitz plot for $D^0\to \bar{K}^0\pi^+\pi^-$ +XLabel=$m^2_{\bar{K}^0\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{\bar{K}^0\pi^-}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}m^2_{\bar{K}^0\pi^+}/{\rm d}m^2_{\bar{K}^0\pi^-}$ [$\rm{GeV}^{-4}$] +END PLOT + +BEGIN PLOT /MC_D_Dalitz/h_minus2 +Title=$K^-\pi^0$ mass distribution in $D^0\to K^-\pi^+\pi^0$ +XLabel=$m^2_{K^-\pi^0}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{K^-\pi^0}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_neutral2 +Title=$K^-\pi^+$ mass distribution in $D^0\to K^-\pi^+\pi^0$ +XLabel=$m^2_{K^-\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{K^-\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pipi2 +Title=$\pi^+\pi^0$ mass distribution in $D^0\to K^-\pi^+\pi^0$ +XLabel=$m^2_{\pi^+\pi^0}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{\pi^+\pi^0}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz2 +Title=Dalitz plot for $D^0\to K^-\pi^+\pi^0$ +XLabel=$m^2_{K^-\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{K^-\pi^0}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}m^2_{K^-\pi^+}/{\rm d}m^2_{K^-\pi^0}$ [$\rm{GeV}^{-4}$] +END PLOT + +BEGIN PLOT /MC_D_Dalitz/h_Kpilow3 +Title=Lower $K^-\pi^+$ mass distribution in $D^+\to K^-\pi^+\pi^+$ +XLabel=${m^{\rm low}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}{m^{\rm low}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_Kpihigh3 +Title=Higher $K^-\pi^+$ mass distribution in $D^+\to K^-\pi^+\pi^+$ +XLabel=${m^{\rm high}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}{m^{\rm high}_{K^-\pi^+}}^$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_Kpiall3 +Title=All $K^-\pi^+$ mass distribution in $D^+\to K^-\pi^+\pi^+$ +XLabel=${m^{\rm all}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}{m^{\rm all}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pipi3 +Title=$\pi^+\pi^+$ mass distribution in $D^+\to K^-\pi^+\pi^+$ +XLabel=$m^2_{\pi^+\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m^2_{\pi^+\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz3 +Title=Dalitz plot for $D^+\to K^-\pi^+\pi^+$ +XLabel=${m^{\rm low}_{K^-\pi^+}}^2$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{\pi^+\pi^+}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}{m^{\rm low}_{K^-\pi^+}}^2/{\rm d}m^2_{\pi^+\pi^+}$ [$\rm{GeV}^{-4}$] +END PLOT + +BEGIN PLOT /MC_D_Dalitz/h_Kpip4 +Title=$\bar{K}^0\pi^+$ mass distribution in $D^+\to \bar{K}^0\pi^+\pi^0$ +XLabel=$m^2_{\bar{K}^0\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{\bar{K}^0\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pipi4 +Title=$\pi^0\pi^+$ mass distribution in $D^+\to \bar{K}^0\pi^+\pi^0$ +XLabel=$m^2_{\pi^0\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{\pi^0\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_Kpi04 +Title=$\bar{K}^0\pi^0$ mass distribution in $D^+\to \bar{K}^0\pi^+\pi^0$ +XLabel=$m^2_{\bar{K}^0\pi^0}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{\bar{K}^0\pi^0}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz4 +Title=Dalitz plot for $D^+\to \bar{K}^0\pi^+\pi^0$ +XLabel=$m^2_{\bar{K}^0\pi^0}$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{\pi^+\pi^0}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}m^2_{\bar{K}^0\pi^0}/{\rm d}m^2_{\pi^+\pi^0}$ [$\rm{GeV}^{-4}$] +END PLOT + +BEGIN PLOT /MC_D_Dalitz/h_kppim5 +Title=$K^+\pi^-$ mass distribution in $D^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_kppip5 +Title=$K^+\pi^+$ mass distribution in $D^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{K^+\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pippim5 +Title=$\pi^-\pi^+$ mass distribution in $D^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{\pi^-\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{\pi^-\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz5 +Title=Dalitz plot for $D^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}m^2_{K^+\pi^-}/{\rm d}m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-4}$] +END PLOT + +BEGIN PLOT /MC_D_Dalitz/h_kppim6 +Title=$K^+\pi^-$ mass distribution in $D_s^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_kppip6 +Title=$K^+\pi^+$ mass distribution in $D_s^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{K^+\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/h_pippim6 +Title=$\pi^-\pi^+$ mass distribution in $D_s^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{\pi^-\pi^+}$ [$\rm{GeV}^{-2}$] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}^{\rm low}_{\pi^-\pi^+}$ [$\rm{GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_D_Dalitz/dalitz6 +Title=Dalitz plot for $D_s^+\to K^+\pi^-\pi^+$ +XLabel=$m^2_{K^+\pi^-}$ [$\rm{GeV}^{-2}$] +YLabel=$m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-2}$] +ZLabel=$1/\Gamma{\rm d}^2 \Gamma/{\rm d}m^2_{K^+\pi^-}/{\rm d}m^2_{\pi^+\pi^-}$ [$\rm{GeV}^{-4}$] +END PLOT diff --git a/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.cc b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.cc @@ -0,0 +1,248 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/UnstableParticles.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class MC_OmegaPhia1_3Pion_Decay : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(MC_OmegaPhia1_3Pion_Decay); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(UnstableParticles(), "UFS"); + // book histograms a_1 + // Histograms for a_10 -> pi0pi0pi0 + _hist0 = bookHisto1D("hist0",200,0.2,1.5); + // // dalitz plot + _dalitz0 = bookHisto2D("dalitz0",50,0.2,1.5,50,0.2,1.5); + // Histograms for a_1+ -> pi0pi0pi+ + // Mass of the pi0pi0 pair + _hist1A = bookHisto1D("hist1A",200,0.2,1.5); + // Mass of the pi0pi+ pair + _hist1B = bookHisto1D("hist1B",200,0.2,1.5); + // dalitz plot + _dalitz1 = bookHisto2D("dalitz1",50,0.2,1.5,50,0.2,1.5); + // Histograms for a_10 -> pi+pi-pi0 + // Mass of the pi+pi- pair + _hist2A = bookHisto1D("hist2A",200,0.2,1.5); + // Mass of the pi+pi0 pair + _hist2B = bookHisto1D("hist2B",200,0.2,1.5); + // Mass of the pi-pi0 pair + _hist2C = bookHisto1D("hist2C",200,0.2,1.5); + // dalitz plot + _dalitz2 = bookHisto2D("dalitz2",50,0.2,1.5,50,0.2,1.5); + // Histograms for a_1+ -> pi+pi+pi- + // Mass of the pi+pi+ pair + _hist3A = bookHisto1D("hist3A",200,0.2,1.5); + // Mass of the pi+pi- pair + _hist3B = bookHisto1D("hist3B",200,0.2,1.5); + // dalitz plot + _dalitz3 = bookHisto2D("dalitz3",50,0.2,1.5,50,0.2,1.5); + + // Book histograms omega/phi + for(unsigned int ix=0;ix<2;++ix) { + double mmax = ix==0 ? 0.8 : 1.0; + ostringstream title1; title1 << "xhist_" << ix+1; + _h_xhist .push_back(bookHisto1D(title1.str(),200,-300.,300. )); + ostringstream title2; title2 << "yhist_" << ix+1; + _h_yhist .push_back(bookHisto1D(title2.str(),200,0. ,400. )); + ostringstream title3; title3 << "mplus_" << ix+1; + _h_mplus .push_back(bookHisto1D(title3.str(),200,200.,mmax*1000.)); + ostringstream title4; title4 << "mminus_" << ix+1; + _h_mminus .push_back(bookHisto1D(title4.str(),200,200.,mmax*1000.)); + ostringstream title5; title5 << "m0_" << ix+1; + _h_m0 .push_back(bookHisto1D(title5.str(),200,200.,mmax*1000.)); + ostringstream title6; title6 << "dalitz_" << ix+1; + _h_dalitz.push_back(bookHisto2D(title6.str(),50,0.2,mmax,50,0.2,mmax)); + } + } + + + void findDecayProducts(const Particle & mother, unsigned int & nstable, + Particles & pip , Particles & pim , Particles & pi0) { + for(const Particle & p : mother.children()) { + int id = p.pdgId(); + if (id == PID::PIPLUS) { + pip.push_back(p); + ++nstable; + } + else if (id == PID::PIMINUS) { + pim.push_back(p); + ++nstable; + } + else if (id == PID::PI0) { + pi0.push_back(p); + ++nstable; + } + else if ( !p.children().empty() ) { + findDecayProducts(p, nstable, pip, pim, pi0); + } + else + ++nstable; + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + double weight = event.weight(); + for(const Particle& meson : + apply(event, "UFS").particles(Cuts::pid==PID::PHI || Cuts::pid==PID::OMEGA || + Cuts::abspid==20213 || Cuts::pid==20113 )) { + unsigned int nstable(0); + Particles pip, pim, pi0; + findDecayProducts(meson, nstable, pip, pim, pi0); + if(nstable !=3) continue; + if(meson.pdgId()<0) { + swap(pim,pip); + } + if(meson.pdgId()== PID::PHI || meson.pdgId()==PID::OMEGA) { + if(pip.size()!=1 || pim.size()!=1 || pi0.size()!=1) continue; + unsigned int iloc = meson.pdgId() == PID::OMEGA ? 0 : 1; + LorentzTransform boost = LorentzTransform::mkFrameTransformFromBeta(meson.momentum().betaVec()); + FourMomentum pp = boost.transform(pip[0].momentum()); + FourMomentum pm = boost.transform(pim[0].momentum()); + FourMomentum p0 = boost.transform(pi0[0].momentum()); + double mp = (pp+p0).mass(), mm = (pm+pp).mass(); + _h_mplus [iloc]->fill(mp/MeV,weight); + _h_mminus[iloc]->fill((pm+p0).mass()/MeV,weight); + _h_m0 [iloc]->fill(mm/MeV,weight); + double x = pp.t()-pm.t(); + double y = p0.t()-p0.mass(); + _h_xhist[iloc]->fill(x/MeV,weight); + _h_yhist[iloc]->fill(y/MeV,weight); + _h_dalitz[iloc]->fill(mp,mm,weight); + } + else { + // a_1+ -> pi+pi+pi- + if(pip.size()==2&&pim.size()==1) { + _hist3A->fill((pip[0].momentum()+pip[1].momentum()).mass(),weight); + _hist3B->fill((pip[0].momentum()+pim[0].momentum()).mass(),weight); + _hist3B->fill((pip[1].momentum()+pim[0].momentum()).mass(),weight); + _dalitz3->fill((pip[0].momentum()+pim[0].momentum()).mass(),(pip[1].momentum()+pim[0].momentum()).mass(),weight); + _dalitz3->fill((pip[1].momentum()+pim[0].momentum()).mass(),(pip[0].momentum()+pim[0].momentum()).mass(),weight); + } + // a_1+ -> pi0pi0pi+ + else if(pip.size()==1&&pi0.size()==2) { + _hist1A->fill((pi0[0].momentum()+pi0[1].momentum()).mass(),weight); + _hist1B->fill((pip[0].momentum()+pi0[0].momentum()).mass(),weight); + _hist1B->fill((pip[0].momentum()+pi0[1].momentum()).mass(),weight); + _dalitz1->fill((pip[0].momentum()+pi0[0].momentum()).mass(),(pip[0].momentum()+pi0[1].momentum()).mass(),weight); + _dalitz1->fill((pip[0].momentum()+pi0[1].momentum()).mass(),(pip[0].momentum()+pi0[0].momentum()).mass(),weight); + } + // a_10 -> pi0pi0pi0 + else if(pi0.size()==3) { + _hist0->fill((pi0[0].momentum()+pi0[1].momentum()).mass(),weight); + _hist0->fill((pi0[0].momentum()+pi0[2].momentum()).mass(),weight); + _hist0->fill((pi0[1].momentum()+pi0[2].momentum()).mass(),weight); + _dalitz0->fill((pi0[0].momentum()+pi0[1].momentum()).mass(),(pi0[0].momentum()+pi0[2].momentum()).mass(),weight); + _dalitz0->fill((pi0[0].momentum()+pi0[1].momentum()).mass(),(pi0[1].momentum()+pi0[2].momentum()).mass(),weight); + _dalitz0->fill((pi0[0].momentum()+pi0[2].momentum()).mass(),(pi0[1].momentum()+pi0[2].momentum()).mass(),weight); + _dalitz0->fill((pi0[0].momentum()+pi0[2].momentum()).mass(),(pi0[0].momentum()+pi0[1].momentum()).mass(),weight); + _dalitz0->fill((pi0[1].momentum()+pi0[2].momentum()).mass(),(pi0[0].momentum()+pi0[1].momentum()).mass(),weight); + _dalitz0->fill((pi0[1].momentum()+pi0[2].momentum()).mass(),(pi0[0].momentum()+pi0[2].momentum()).mass(),weight); + } + // a_10 -> pi+pi-pi0 + else if(pi0.size()==1&&pip.size()==1&&pim.size()==1) { + _hist2A->fill((pim[0].momentum()+pip[0].momentum()).mass(),weight); + _hist2B->fill((pip[0].momentum()+pi0[0].momentum()).mass(),weight); + _hist2C->fill((pim[0].momentum()+pi0[0].momentum()).mass(),weight); + _dalitz2->fill((pim[0].momentum()+pi0[0].momentum()).mass(),(pip[0].momentum()+pi0[0].momentum()).mass(),weight); + } + } + } + } + + + /// Normalise histograms etc., after the run + void finalize() { + // a_1 + normalize(_hist0); + normalize(_dalitz0); + normalize(_hist1A); + normalize(_hist1B); + normalize(_dalitz1); + normalize(_hist2A); + normalize(_hist2B); + normalize(_hist2C); + normalize(_dalitz2); + normalize(_hist3A); + normalize(_hist3B); + normalize(_dalitz3); + // omega/phi + for(unsigned int ix=0;ix<2;++ix) { + normalize(_h_xhist [ix]); + normalize(_h_yhist [ix]); + normalize(_h_mplus [ix]); + normalize(_h_mminus[ix]); + normalize(_h_m0 [ix]); + normalize(_h_dalitz[ix]); + } + } + + //@} + + /// @name Histograms a_1 + //@{ + // Histograms for a_10 -> pi0pi0pi0 + Histo1DPtr _hist0; + // dalitz plot + Histo2DPtr _dalitz0; + // Histograms for a_1+ -> pi0pi0pi+ + // Mass of the pi0pi0 pair + Histo1DPtr _hist1A; + // Mass of the pi0pi+ pair + Histo1DPtr _hist1B; + // dalitz plot + Histo2DPtr _dalitz1; + // Histograms for a_10 -> pi+pi-pi0 + // Mass of the pi+pi- pair + Histo1DPtr _hist2A; + // Mass of the pi+pi0 pair + Histo1DPtr _hist2B; + // Mass of the pi-pi0 pair + Histo1DPtr _hist2C; + // dalitz plot + Histo2DPtr _dalitz2; + // Histograms for a_1+ -> pi+pi+pi- + // Mass of the pi+pi+ pair + Histo1DPtr _hist3A; + // Mass of the pi+pi- pair + Histo1DPtr _hist3B; + // dalitz plot + Histo2DPtr _dalitz3; + //@} + + /// @name Histograms omega/phi + //@{ + // Histogram for the x-values + vector _h_xhist; + // Histogram for the y-values + vector _h_yhist; + // The mass of the \rho^+ + vector _h_mplus; + // The mass of the \rho^- + vector _h_mminus; + // The mass of the \rho^0 + vector _h_m0; + // Dalitz plot + vector _h_dalitz; + //@} + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(MC_OmegaPhia1_3Pion_Decay); + + +} diff --git a/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.info b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.info new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.info @@ -0,0 +1,10 @@ +Name: MC_OmegaPhia1_3Pion_Decay +Summary: MC analysis of $\omega$, $\phi$ and $a_1\to3\pi$ decays +Status: VALIDATED +Authors: + - Peter Richardson +RunInfo: Any process producing a1, omega or phi mesons +NeedCrossSection: no +Description: + 'Analysis of the mass distributions and Dalitz plots in $a_1\to3\pi$ and $\omega,\phi\to\pi^+\pi^-\pi^0$ decays. + Based on old Herwig++ internal analysis.' diff --git a/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.plot b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_OmegaPhia1_3Pion_Decay.plot @@ -0,0 +1,130 @@ +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist0 +Title=$\pi^0\pi^0$ mass distribution in $a_1^0\to3\pi^0$ decays +XLabel=$m_{\pi^0\pi^0}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^0\pi^0}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz0 +Title=Dalitz plot for $a_1^0\to3\pi^0$ decays +XLabel=$m_{\pi^0\pi^0}$ [GeV] +YLabel=$m_{\pi^0\pi^0}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^0\pi^0}/{\rm d}m_{\pi^0\pi^0}$ [${\rm GeV}^{-2}$] +END PLOT + +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist1A +Title=$\pi^0\pi^0$ mass distribution in $a_1^+\to\pi^+2\pi^0$ decays +XLabel=$m_{\pi^0\pi^0}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^0\pi^0}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist1B +Title=$\pi^0\pi^+$ mass distribution in $a_1^+\to\pi^+2\pi^0$ decays +XLabel=$m_{\pi^0\pi^+}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^0\pi^+}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz1 +Title=Dalitz plot for $a_1^+\to\pi^+2\pi^0$ decays +XLabel=$m_{\pi^+\pi^0}$ [GeV] +YLabel=$m_{\pi^+\pi^0}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^0}/{\rm d}m_{\pi^+\pi^0}$ [${\rm GeV}^{-2}$] +END PLOT + +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist2A +Title=$\pi^+\pi^-$ mass distribution in $a_1^0\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^-}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^-}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist2B +Title=$\pi^+\pi^0$ mass distribution in $a_1^0\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^0}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^0}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist2C +Title=$\pi^-\pi^0$ mass distribution in $a_1^0\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^-\pi^0}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^-\pi^0}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz2 +Title=Dalitz plot for $a_1^0\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^-\pi^0}$ [GeV] +YLabel=$m_{\pi^+\pi^0}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^0}/{\rm d}m_{\pi^+\pi^0}$ [${\rm GeV}^{-2}$] +END PLOT + +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist3A +Title=$\pi^+\pi^+$ mass distribution in $a_1^+\to2\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^+}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^+}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/hist3B +Title=$\pi^-\pi^+$ mass distribution in $a_1^+\to2\pi^+\pi^-$ decays +XLabel=$m_{\pi^-\pi^+}$ [GeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^-\pi^+}$ [${\rm GeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz3 +Title=Dalitz plot for $a_1^+\to2\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^-}$ [GeV] +YLabel=$m_{\pi^+\pi^-}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d}m_{\pi^+\pi^0}/{\rm d}m_{\pi^+\pi^0}$ [${\rm GeV}^{-2}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/xhist_1 +Title=$x$ distribution in $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$x$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} x$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/yhist_1 +Title=$y$ distribution in $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$y$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} y$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/mminus_1 +Title=$\pi^-\pi^0$ mass distribution in $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^-\pi^0}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^-\pi^0}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/mplus_1 +Title=$\pi^+\pi^0$ mass distribution in $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^0}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^0}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/m0_1 +Title=$\pi^+\pi^-$ mass distribution in $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^-}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^-}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz_1 +Title=Dalitz plot for $\omega\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^0}$ [GeV] +YLabel=$m_{\pi^+-\pi^0}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^0}/{\rm d} m_{\pi^-\pi^0}$ [$\rm{GeV}^{-2}$] +END PLOT + +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/xhist_2 +Title=$x$ distribution in $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$x$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} x$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/yhist_2 +Title=$y$ distribution in $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$y$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} y$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/mminus_2 +Title=$\pi^-\pi^0$ mass distribution in $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^-\pi^0}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^-\pi^0}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/mplus_2 +Title=$\pi^+\pi^0$ mass distribution in $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^0}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^0}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/m0_2 +Title=$\pi^+\pi^-$ mass distribution in $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^-}$ [MeV] +YLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^-}$ [$\rm{MeV}^{-1}$] +END PLOT +BEGIN PLOT /MC_OmegaPhia1_3Pion_Decay/dalitz_2 +Title=Dalitz plot for $\phi\to\pi^0\pi^+\pi^-$ decays +XLabel=$m_{\pi^+\pi^0}$ [GeV] +YLabel=$m_{\pi^+-\pi^0}$ [GeV] +ZLabel=$1/\Gamma{\rm d}\Gamma/{\rm d} m_{\pi^+\pi^0}/{\rm d} m_{\pi^-\pi^0}$ [$\rm{GeV}^{-2}$] +END PLOT diff --git a/analyses/pluginMC/MC_Semi_Leptonic_Decay.cc b/analyses/pluginMC/MC_Semi_Leptonic_Decay.cc new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_Semi_Leptonic_Decay.cc @@ -0,0 +1,180 @@ +// -*- C++ -*- +#include "Rivet/Analysis.hh" +#include "Rivet/Projections/UnstableParticles.hh" + +namespace Rivet { + + + /// @brief Add a short analysis description here + class MC_Semi_Leptonic_Decay : public Analysis { + public: + + /// Constructor + DEFAULT_RIVET_ANALYSIS_CTOR(MC_Semi_Leptonic_Decay); + + + /// @name Analysis methods + //@{ + + /// Book histograms and initialise projections before the run + void init() { + + // Initialise and register projections + declare(UnstableParticles(), "UFS"); + + } + + void findDecayProducts(const Particle & mother, + unsigned int & nstable, + Particles& lp, Particles& lm, + Particles& nu, Particles& nub, + Particles& out) { + for(const Particle & p : mother.children()) { + int id = p.pdgId(); + if ( id == PID::EMINUS || id == PID::MUON ) { + lm .push_back(p); + ++nstable; + } + else if (id == PID::EPLUS || id == PID::ANTIMUON) { + lp .push_back(p); + ++nstable; + } + else if ( id == PID::NU_E || id == PID::NU_EBAR ) { + nu .push_back(p); + ++nstable; + } + else if (id == PID::NU_MU || id == PID::NU_MUBAR ) { + nub.push_back(p); + ++nstable; + } + else if (PID::isMeson(id)) { + out.push_back(p); + ++nstable; + } + else if ( !p.children().empty() ) { + findDecayProducts(p,nstable,lp,lm,nu,nub,out); + } + else + ++nstable; + } + } + + /// Perform the per-event analysis + void analyze(const Event& event) { + double weight = event.weight(); + // loop over unstable particles + for(const Particle& meson : apply(event, "UFS").particles()) { + int id = meson.pdgId(); + // spin 0 mesons + if(!PID::isMeson(id)) continue; + if(abs(id)%10!=1) continue; + unsigned int nstable(0); + Particles lp, lm, nu, nub, out; + findDecayProducts(meson,nstable,lp,lm,nu,nub,out); + if(nstable!=3 || out.size()!=1) continue; + int ilep=0; + FourMomentum plep,pmnu=out[0].momentum(); + double me2(0.); + if( lp.size()==1 && nu.size()==1 && out.size()==1 ) { + if(nu[0].pdgId() != -lp[0].pdgId()+1) continue; + ilep = lp[0].pdgId(); + plep = nu[0].momentum()+lp[0].momentum(); + pmnu += nu[0].momentum(); + me2 = lp[0].mass2(); + } + else if( lm.size()==1 && nub.size()==1 && out.size()==1 ) { + if(nub[0].pdgId() != -lm[0].pdgId()-1) continue; + ilep = lm[0].pdgId(); + plep = nub[0].momentum()+lm[0].momentum(); + pmnu += nub[0].momentum(); + me2 = lm[0].mass2(); + } + else + continue; + // check if histos already exist + unsigned int iloc=0; bool found(false); + while(!found&&iloc<_incoming.size()) { + if(_incoming[iloc] == id && + _outgoing[iloc] == out[0].pdgId() && + ilep==_outgoingL[iloc]) found=true; + else ++iloc; + } + if(!found) { + iloc=_incoming.size(); + _incoming.push_back(id); + _outgoing.push_back(out[0].pdgId()); + _outgoingL.push_back(ilep); + ostringstream title; + title << "h_" << abs(id); + if(id>0) title << "p"; + else title << "m"; + title << "_" << abs(out[0].pdgId()); + if(out[0].pdgId()>0) title << "p"; + else title << "m"; + title << "_" << abs(ilep); + if(ilep>0) title << "p"; + else title << "m"; + title << "_"; + _energy.push_back(bookHisto1D(title.str()+"energy", + 200,0.0,0.5*meson.mass()/MeV)); + _scale .push_back(bookHisto1D(title.str()+"scale", + 200,0.0,meson.mass()/MeV)); + } + // add the results to the histogram + _scale[iloc]->fill(plep.mass()/MeV,weight); + double ee = 0.5/meson.mass()*(meson.mass2()-pmnu.mass2()+me2); + _energy[iloc]->fill(ee/MeV,weight); + } + } + + /// Normalise histograms etc., after the run + void finalize() { + for(unsigned int ix=0;ix<_energy.size();++ix) { + normalize(_energy); + normalize(_scale ); + } + } + + //@} + + + /// @name Histograms + //@{ + /** + * PDG codes of the decaying mesons + */ + vector _incoming; + + /** + * PDG codes of the decay products + */ + vector _outgoing; + + /** + * Identidies of the leptons + */ + vector _outgoingL; + + /** + * Histograms + */ + //@{ + /** + * The lepton energy + */ + vector _energy; + + /** + * The \f$q\f$ value + */ + vector _scale; + //@} + + }; + + + // The hook for the plugin system + DECLARE_RIVET_PLUGIN(MC_Semi_Leptonic_Decay); + + +} diff --git a/analyses/pluginMC/MC_Semi_Leptonic_Decay.info b/analyses/pluginMC/MC_Semi_Leptonic_Decay.info new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_Semi_Leptonic_Decay.info @@ -0,0 +1,9 @@ +Name: MC_Semi_Leptonic_Decay +Summary: MC analysis of the Kinematics in Semi-Leptonic Meson decays +Status: VALIDATED +Authors: + - Peter Richardson +RunInfo: Any type of process producing mesons which decay semi-leptonically +NeedCrossSection: no +Description: + 'Analysis of the kinematics in semileptonic meson decays.' diff --git a/analyses/pluginMC/MC_Semi_Leptonic_Decay.plot b/analyses/pluginMC/MC_Semi_Leptonic_Decay.plot new file mode 100644 --- /dev/null +++ b/analyses/pluginMC/MC_Semi_Leptonic_Decay.plot @@ -0,0 +1,972 @@ +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_13m_energy +XLabel=$E_{\mu^+}$ [Me]V +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}E_{\mu^+}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_13m_scale +XLabel=$m_{\mu^+\nu_\mu}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m_{\mu^+\nu_\mu}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_13p_energy +XLabel=$E_{\mu^-}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}E_{\mu^-}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_13p_scale +XLabel=$m_{\mu^-\bar\nu_\mu}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m_{\mu^-\bar\nu_\mu}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_11m_energy +XLabel=$E_{e^+}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}E_{e^+}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_11m_scale +XLabel=$m_{e^+\nu_e}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m_{e^+\nu_e}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_11p_energy +XLabel=$E_{e^-}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}E_{e^-}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_.*_11p_scale +XLabel=$m_{e^-\bar\nu_e}$ [MeV] +YLabel=$1/\Gamma{\rm d} \Gamma/{\rm d}m_{e^-\bar\nu_e}$ [${\rm MeV}^{-1}$] +LogY=0 +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_10313p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to K_1(1270)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_10313p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to K_1(1270)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_10313p_11p_energy +Title=$e^-$ energy in the decay $D^-\to K_1(1270)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_10313p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to K_1(1270)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_10313m_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \bar{K}_1(1270)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_10313m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \bar{K}_1(1270)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_10313m_11m_energy +Title=$e^+$ energy in the decay $D^+\to \bar{K}_1(1270)^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_10313m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \bar{K}_1(1270)^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_313p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to K^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_313p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to K^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_313p_11p_energy +Title=$e^-$ energy in the decay $D^-\to K^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_313p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to K^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_313m_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \bar{K}^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_313m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \bar{K}^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_313m_11m_energy +Title=$e^+$ energy in the decay $D^+\to \bar{K}^{*0} e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_313m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \bar{K}^{*0} e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_311p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to K^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_311p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to K^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_311p_11p_energy +Title=$e^-$ energy in the decay $D^-\to K^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_311p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to K^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_311m_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \bar{K}^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_311m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \bar{K}^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_311m_11m_energy +Title=$e^+$ energy in the decay $D^+\to \bar{K}^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_311m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \bar{K}^0 e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_111p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to \pi^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_111p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to \pi^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_111p_11p_energy +Title=$e^-$ energy in the decay $D^-\to \pi^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_111p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to \pi^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_111p_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \pi^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_111p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \pi^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_111p_11m_energy +Title=$e^+$ energy in the decay $D^+\to \pi^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_111p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \pi^0 e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_221p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_221p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_221p_11p_energy +Title=$e^-$ energy in the decay $D^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_221p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_221p_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_221p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_221p_11m_energy +Title=$e^+$ energy in the decay $D^+\to \eta e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_221p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \eta e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_223p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_223p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_223p_11p_energy +Title=$e^-$ energy in the decay $D^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_223p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_223p_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_223p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_223p_11m_energy +Title=$e^+$ energy in the decay $D^+\to \omega e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_223p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \omega e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_113p_13p_energy +Title=$\mu^-$ energy in the decay $D^-\to \rho^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_113p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^-\to \rho^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_113p_11p_energy +Title=$e^-$ energy in the decay $D^-\to \rho^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411m_113p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^-\to \rho^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_113p_13m_energy +Title=$\mu^+$ energy in the decay $D^+\to \rho^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_113p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D^+\to \rho^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_113p_11m_energy +Title=$e^+$ energy in the decay $D^+\to \rho^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_411p_113p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D^+\to \rho^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_10323p_13p_energy +Title=$\mu^-$ energy in the decay $D^0\to K_1(1270)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_10323p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^0\to K_1(1270)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_10323p_11p_energy +Title=$e^-$ energy in the decay $D^0\to K_1(1270)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_10323p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^0\to K_1(1270)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_10323m_13m_energy +Title=$\mu^+$ energy in the decay $\bar{D}^0\to K_1(1270)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_10323m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $\bar{D}^0\to K_1(1270)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_10323m_11m_energy +Title=$e^+$ energy in the decay $\bar{D}^0\to K_1(1270)^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_10323m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $\bar{D}^0\to K_1(1270)^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_323p_13p_energy +Title=$\mu^-$ energy in the decay $D^0\to K^{*+} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_323p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^0\to K^{*+} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_323p_11p_energy +Title=$e^-$ energy in the decay $D^0\to K^{*+} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_323p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^0\to K^{*+} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_323m_13m_energy +Title=$\mu^+$ energy in the decay $\bar{D}^0\to K^{*-} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_323m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $\bar{D}^0\to K^{*-} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_323m_11m_energy +Title=$e^+$ energy in the decay $\bar{D}^0\to K^{*-} e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_323m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $\bar{D}^0\to K^{*-} e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_321p_13p_energy +Title=$\mu^-$ energy in the decay $D^0\to K^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_321p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^0\to K^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_321p_11p_energy +Title=$e^-$ energy in the decay $D^0\to K^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_321p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^0\to K^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_321m_13m_energy +Title=$\mu^+$ energy in the decay $\bar{D}^0\to K^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_321m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $\bar{D}^0\to K^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_321m_11m_energy +Title=$e^+$ energy in the decay $\bar{D}^0\to K^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_321m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $\bar{D}^0\to K^- e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_211p_13p_energy +Title=$\mu^-$ energy in the decay $D^0\to \pi^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_211p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^0\to \pi^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_211p_11p_energy +Title=$e^-$ energy in the decay $D^0\to \pi^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_211p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^0\to \pi^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_211m_13m_energy +Title=$\mu^+$ energy in the decay $\bar{D}^0\to \pi^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_211m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $\bar{D}^0\to \pi^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_211m_11m_energy +Title=$e^+$ energy in the decay $\bar{D}^0\to \pi^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_211m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $\bar{D}^0\to \pi^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_213p_13p_energy +Title=$\mu^-$ energy in the decay $D^0\to \rho^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_213p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D^0\to \rho^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_213p_11p_energy +Title=$e^-$ energy in the decay $D^0\to \rho^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421m_213p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D^0\to \rho^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_213m_13m_energy +Title=$\mu^+$ energy in the decay $\bar{D}^0\to \rho^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_213m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $\bar{D}^0\to \rho^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_213m_11m_energy +Title=$e^+$ energy in the decay $\bar{D}^0\to \rho^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_421p_213m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $\bar{D}^0\to \rho^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_311m_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \bar{K}^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_311m_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \bar{K}^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_311m_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \bar{K}^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_311m_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \bar{K}^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_311p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to K^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_311p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to K^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_311p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to K^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_311p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to K^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_313m_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \bar{K}^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_313m_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \bar{K}^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_313m_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \bar{K}^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_313m_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \bar{K}^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_313p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to K^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_313p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to K^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_313p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to K^{*0} e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_313p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to K^{*0} e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_221p_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_221p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_221p_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_221p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_221p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_221p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_221p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to \eta e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_221p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to \eta e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_331p_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \eta^\prime \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_331p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \eta^\prime \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_331p_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \eta^\prime e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_331p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \eta^\prime e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_331p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to \eta^\prime \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_331p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to \eta^\prime \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_331p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to \eta^\prime e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_331p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to \eta^\prime e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_223p_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_223p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_223p_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_223p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_223p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_223p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_223p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to \omega e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_223p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to \omega e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_333p_13p_energy +Title=$\mu^-$ energy in the decay $D_s^-\to \phi \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_333p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $D_s^-\to \phi \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_333p_11p_energy +Title=$e^-$ energy in the decay $D_s^-\to \phi e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431m_333p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $D_s^-\to \phi e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_333p_13m_energy +Title=$\mu^+$ energy in the decay $D_s^+\to \phi \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_333p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $D_s^+\to \phi \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_333p_11m_energy +Title=$e^+$ energy in the decay $D_s^+\to \phi e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_431p_333p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $D_s^+\to \phi e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10411p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D^*_0(2400)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10411p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D^*_0(2400)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10411p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D^*_0(2400)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10411p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D^*_0(2400)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10411m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D^*_0(2400)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10411m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D^*_0(2400)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10411m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D^*_0(2400)^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10411m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D^*_0(2400)^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_20413p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D_1(2430)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_20413p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D_1(2430)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_20413p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D_1(2430)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_20413p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D_1(2430)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_20413m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D_1(2430)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_20413m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D_1(2430)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_20413m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D_1(2430)^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_20413m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D_1(2430)^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10413p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D_1(2420)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10413p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D_1(2420)^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10413p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D_1(2420)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_10413p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D_1(2420)^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10413m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D_1(2420)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10413m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D_1(2420)^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10413m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D_1(2420)^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_10413m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D_1(2420)^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_411p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_411p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_411p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_411p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_411m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_411m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_411m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_411m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_415p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D_2^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_415p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D_2^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_415p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D_2^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_415p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D_2^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_415m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D_2^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_415m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D_2^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_415m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D_2^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_415m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D_2^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_413p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to D^{*+} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_413p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to D^{*+} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_413p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to D^{*+} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_413p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to D^{*+} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_413m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to D^{*-} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_413m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to D^{*-} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_413m_11m_energy +Title=$e^+$ energy in the decay $B^0\to D^{*-} e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_413m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to D^{*-} e^+\nu_e$ +END PLOT + + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_211p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to \pi^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_211p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to \pi^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_211p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to \pi^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_211p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to \pi^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_211m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to \pi^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_211m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to \pi^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_211m_11m_energy +Title=$e^+$ energy in the decay $B^0\to \pi^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_211m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to \pi^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_213p_13p_energy +Title=$\mu^-$ energy in the decay $\bar{B}^0\to \rho^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_213p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $\bar{B}^0\to \rho^+ \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_213p_11p_energy +Title=$e^-$ energy in the decay $\bar{B}^0\to \rho^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511m_213p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $\bar{B}^0\to \rho^+ e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_213m_13m_energy +Title=$\mu^+$ energy in the decay $B^0\to \rho^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_213m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^0\to \rho^- \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_213m_11m_energy +Title=$e^+$ energy in the decay $B^0\to \rho^- e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_511p_213m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^0\to \rho^- e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10421p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D^*_0(2400)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10421p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D^*_0(2400)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10421p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D^*_0(2400)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10421p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D^*_0(2400)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10421m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}^*_0(2400)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10421m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}^*_0(2400)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10421m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}^*_0(2400)^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10421m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}^*_0(2400)^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10423p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D_1(2420)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10423p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D_1(2420)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10423p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D_1(2420)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_10423p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D_1(2420)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10423m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}_1(2420)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10423m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}_1(2420)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10423m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}_1(2420)^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_10423m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}_1(2420)^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_20423p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D_1(2430)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_20423p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D_1(2430)^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_20423p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D_1(2430)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_20423p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D_1(2430)^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_20423m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}_1(2430)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_20423m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}_1(2430)^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_20423m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}_1(2430)^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_20423m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}_1(2430)^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_421p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_421p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_421p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_421p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_421m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_421m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_421m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_421m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_425p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D_2^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_425p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D_2^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_425p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D_2^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_425p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D_2^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_425m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}_2^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_425m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}_2^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_425m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}_2^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_425m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}_2^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_423p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to D^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_423p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to D^{*0} \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_423p_11p_energy +Title=$e^-$ energy in the decay $B^-\to D^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_423p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to D^{*0} e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_423m_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \bar{D}^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_423m_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \bar{D}^{*0} \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_423m_11m_energy +Title=$e^+$ energy in the decay $B^+\to \bar{D}^{*0} e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_423m_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \bar{D}^{*0} e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_113p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to \rho^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_113p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to \rho^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_113p_11p_energy +Title=$e^-$ energy in the decay $B^-\to \rho^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_113p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to \rho^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_113p_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \rho^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_113p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \rho^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_113p_11m_energy +Title=$e^+$ energy in the decay $B^+\to \rho^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_113p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \rho^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_223p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_223p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to \omega \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_223p_11p_energy +Title=$e^-$ energy in the decay $B^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_223p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to \omega e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_223p_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_223p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \omega \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_223p_11m_energy +Title=$e^+$ energy in the decay $B^+\to \omega e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_223p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \omega e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_111p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to \pi^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_111p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to \pi^0 \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_111p_11p_energy +Title=$e^-$ energy in the decay $B^-\to \pi^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_111p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to \pi^0 e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_111p_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \pi^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_111p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \pi^0 \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_111p_11m_energy +Title=$e^+$ energy in the decay $B^+\to \pi^0 e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_111p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \pi^0 e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_221p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_221p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to \eta \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_221p_11p_energy +Title=$e^-$ energy in the decay $B^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_221p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to \eta e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_221p_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_221p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \eta \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_221p_11m_energy +Title=$e^+$ energy in the decay $B^+\to \eta e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_221p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \eta e^+\nu_e$ +END PLOT + +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_331p_13p_energy +Title=$\mu^-$ energy in the decay $B^-\to \eta^\prime \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_331p_13p_scale +Title=Mass of $\mu^-\bar\nu_\mu$ in the decay $B^-\to \eta^\prime \mu^-\bar\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_331p_11p_energy +Title=$e^-$ energy in the decay $B^-\to \eta^\prime e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521m_331p_11p_scale +Title=Mass of $e^-\bar\nu_e$ in the decay $B^-\to \eta^\prime e^-\bar\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_331p_13m_energy +Title=$\mu^+$ energy in the decay $B^+\to \eta^\prime \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_331p_13m_scale +Title=Mass of $\mu^+\nu_\mu$ in the decay $B^+\to \eta^\prime \mu^+\nu_\mu$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_331p_11m_energy +Title=$e^+$ energy in the decay $B^+\to \eta^\prime e^+\nu_e$ +END PLOT +BEGIN PLOT /MC_Semi_Leptonic_Decay/h_521p_331p_11m_scale +Title=Mass of $e^+\nu_e$ in the decay $B^+\to \eta^\prime e^+\nu_e$ +END PLOT