diff --git a/Decay/HeavyMeson/HQETStrongDecayer.cc b/Decay/HeavyMeson/HQETStrongDecayer.cc
--- a/Decay/HeavyMeson/HQETStrongDecayer.cc
+++ b/Decay/HeavyMeson/HQETStrongDecayer.cc
@@ -1,283 +1,364 @@
 // -*- C++ -*-
 //
 // This is the implementation of the non-inlined, non-templated member
 // functions of the HQETStrongDecayer class.
 //
 
 #include "HQETStrongDecayer.h"
 #include "ThePEG/Interface/ClassDocumentation.h"
 #include "ThePEG/Interface/ParVector.h"
 #include "ThePEG/EventRecord/Particle.h"
 #include "ThePEG/Utilities/DescribeClass.h"
 #include "ThePEG/Persistency/PersistentOStream.h"
 #include "ThePEG/Persistency/PersistentIStream.h"
+#include "ThePEG/Helicity/WaveFunction/TensorWaveFunction.h"
+#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
 #include "ThePEG/Helicity/WaveFunction/ScalarWaveFunction.h"
-#include "ThePEG/Helicity/WaveFunction/VectorWaveFunction.h"
-#include "ThePEG/Helicity/WaveFunction/TensorWaveFunction.h"
 #include "Herwig/Decay/TwoBodyDecayMatrixElement.h"
 #include "Herwig/Utilities/Kinematics.h"
+#include "ThePEG/Helicity/epsilon.h"
+#include "ThePEG/Helicity/HelicityFunctions.h"
 
 using namespace Herwig;
 
 HQETStrongDecayer::HQETStrongDecayer()
-  : fPi_(130.2*MeV), g_(0.565), deltaEta_(1./43.7),
-    incoming_ ({413,413,423,433}),
-    outgoingH_({421,411,421,431}),
-    outgoingL_({211,111,111,111}),
-    type_     ({  1,  1,  1, -1}),
-    maxWeight_({ 1., 1., 1., 1.})
+  : fPi_(130.2*MeV), g_(0.565), h_(0.565), deltaEta_(1./43.7), Lambda_(1.*GeV),
+    incoming_ ({413,413,423,433,                   //D* decay modes   [0-3],   VtoSS
+                415,415,425,425,435,435,           //D*_2 decay modes [4-9],   TtoSS
+                435,435,415,415,425,425,435,435}), //D*_2 decay modes [10-17], TtoVS
+    outgoingH_({421,411,421,431,
+                411,421,411,421,411,421,
+                413,423,413,423,413,423,413,423}),
+    outgoingL_({211,111,111,111,
+                111,211,-211,111,311,321,
+                311,321,111,211,-211,111,311,321}),
+    type_     ({1,  1,  1, -1,
+                2,  2,  2,  2,  2,  2,
+                3,  3,  3,  3,  3,  3,  3,  3}),
+    maxWeight_({1., 1., 1., 1.,
+                1., 1., 1., 1., 1., 1., 1., 1.,
+                1., 1., 1., 1., 1., 1.})
 {}
 
 void HQETStrongDecayer::doinit() {
   DecayIntegrator::doinit();
   // check consistence of the parameters
   unsigned int isize=incoming_.size();
   if(isize!=outgoingH_.size()||isize!=outgoingL_.size()||
      isize!=maxWeight_.size()||isize!=type_     .size())
-    throw InitException() << "Inconsistent parameters HQETStrongDecayer" 
+    throw InitException() << "Inconsistent parameters in HQETStrongDecayer"
   			  << Exception::abortnow;
   // set up the integration channels
-  vector<double> wgt(0);
   PhaseSpaceModePtr mode;
   for(unsigned int ix=0;ix<incoming_.size();++ix) {
      tPDPtr     in = getParticleData(incoming_[ix]);
      tPDVector out = {getParticleData(outgoingH_[ix]),
   		      getParticleData(outgoingL_[ix])};
-    if(in&&out[0]&&out[1]) 
+    if(in&&out[0]&&out[1])
       mode=new_ptr(PhaseSpaceMode(in,out,maxWeight_[ix]));
     else
       mode=PhaseSpaceModePtr();
     addMode(mode);
   }
 }
 
 void HQETStrongDecayer::doinitrun() {
   DecayIntegrator::doinitrun();
   if(initialize()) {
     for(unsigned int ix=0;ix<incoming_.size();++ix)
       if(mode(ix)) maxWeight_[ix] = mode(ix)->maxWeight();
   }
 }
 
 IBPtr HQETStrongDecayer::clone() const {
   return new_ptr(*this);
 }
 
 IBPtr HQETStrongDecayer::fullclone() const {
   return new_ptr(*this);
 }
 
 void HQETStrongDecayer::persistentOutput(PersistentOStream & os) const {
-  os << ounit(fPi_,MeV) << g_ << deltaEta_ << maxWeight_;
+  os << ounit(fPi_,MeV) << g_ << h_ << deltaEta_ << ounit(Lambda_,GeV) << maxWeight_;
 }
 
 void HQETStrongDecayer::persistentInput(PersistentIStream & is, int) {
-  is >> iunit(fPi_,MeV) >> g_ >> deltaEta_ >> maxWeight_;
+  is >> iunit(fPi_,MeV) >> g_ >> h_ >> deltaEta_ >> iunit(Lambda_,GeV) >> maxWeight_;
 }
 
 // The following static variable is needed for the type
-// description system in ThePEG. 
+// description system in ThePEG.
 DescribeClass<HQETStrongDecayer,DecayIntegrator>
 describeHerwigHQETStrongDecayer("Herwig::HQETStrongDecayer", "HwHMDecay.so");
 
 void HQETStrongDecayer::Init() {
 
   static ClassDocumentation<HQETStrongDecayer> documentation
     ("The HQETStrongDecayer class performs the strong decays of excited heavy mesons using HQET results.");
-  
+
   static Parameter<HQETStrongDecayer,Energy> interfacefPi
     ("fPi",
      "The pion decay constant",
      &HQETStrongDecayer::fPi_, MeV, 130.2*MeV, 100.0*MeV, 200.0*MeV,
      false, false, Interface::limited);
 
   static Parameter<HQETStrongDecayer,double> interfaceg
     ("g",
      "The coupling for 1- to 0- decays",
      &HQETStrongDecayer::g_, 0.565, 0.0, 1.0,
      false, false, Interface::limited);
 
+  static Parameter<HQETStrongDecayer,double> interfaceh
+    ("h",
+    "The coupling for D*_2 and D*_2s decays",
+    &HQETStrongDecayer::h_, 0.565, 0.0, 1.0,
+    false, false, Interface::limited);
+
   static ParVector<HQETStrongDecayer,double> interfaceMaxWeight
     ("MaxWeight",
      "The maximum weight for the decay mode",
      &HQETStrongDecayer::maxWeight_,
      0, 0, 0, 0., 100000., false, false, true);
-  
+
   static Parameter<HQETStrongDecayer,double> interfaceDeltaEta
     ("DeltaEta",
      "The mixing parameter for eta-pi0 of isospin violating decays",
      &HQETStrongDecayer::deltaEta_, 1./43.7, 0.0, 1.,
      false, false, Interface::limited);
+
+  static Parameter<HQETStrongDecayer,Energy> interfacefLambda
+    ("Lambda",
+    "D*_2 and D*_2s strong decays momentum scale",
+    &HQETStrongDecayer::Lambda_, GeV, 1.*GeV, .1*GeV, 2.*GeV,
+    false, false, Interface::limited);
 }
 
 int HQETStrongDecayer::modeNumber(bool & cc,tcPDPtr parent,
 				  const tPDVector & children) const {
   if(children.size()!=2) return -1;
   int id(parent->id());
   int idbar = parent->CC() ? parent->CC()->id() : id;
   int id1(children[0]->id());
   int id1bar = children[0]->CC() ? children[0]->CC()->id() : id1;
   int id2(children[1]->id());
   int id2bar = children[1]->CC() ? children[1]->CC()->id() : id2;
   int imode(-1);
   unsigned int ix(0);
   cc=false;
   do {
     if(id   ==incoming_[ix]) {
       if((id1   ==outgoingH_[ix]&&id2   ==outgoingL_[ix])||
   	 (id2   ==outgoingH_[ix]&&id1   ==outgoingL_[ix])) imode=ix;
     }
     if(idbar==incoming_[ix]) {
       if((id1bar==outgoingH_[ix]&&id2bar==outgoingL_[ix])||
   	 (id2bar==outgoingH_[ix]&&id1bar==outgoingL_[ix])) {
   	imode=ix;
   	cc=true;
       }
     }
     ++ix;
   }
   while(ix<incoming_.size()&&imode<0);
   return imode;
 }
 
 void HQETStrongDecayer::
 constructSpinInfo(const Particle & part, ParticleVector decay) const {
   switch (part.dataPtr()->iSpin()) {
   case PDT::Spin1:
     Helicity::VectorWaveFunction::constructSpinInfo(vecIn_,const_ptr_cast<tPPtr>(&part),
 						    Helicity::incoming,true,false);
     break;
+  case PDT::Spin2:
+    Helicity::TensorWaveFunction::constructSpinInfo(tensorIn_,const_ptr_cast<tPPtr>(&part),
+                Helicity::incoming,true,false);
+    break;
   default:
     assert(false);
   }
   // set up the spin information for the decay products
   for(unsigned int ix=0;ix<decay.size();++ix) {
     switch (decay[ix]->dataPtr()->iSpin()) {
+    case PDT::Spin1:
+      Helicity::VectorWaveFunction::constructSpinInfo(vecOut_,decay[ix],
+              Helicity::outgoing,true,false);
+      break;
     case PDT::Spin0:
       Helicity::ScalarWaveFunction::constructSpinInfo(decay[ix],Helicity::outgoing,true);
       break;
     default:
       assert(false);
     }
   }
 }
 
 // matrix elememt for the process
 double HQETStrongDecayer::me2(const int, const Particle & part,
 			      const tPDVector & outgoing,
 			      const vector<Lorentz5Momentum> & momenta,
 			      MEOption meopt) const {
   if(!ME()) {
     if(abs(type_[imode()])==1) {
       ME(new_ptr(TwoBodyDecayMatrixElement(PDT::Spin1,PDT::Spin0,PDT::Spin0)));
     }
+    if(abs(type_[imode()])==2) {
+      ME(new_ptr(TwoBodyDecayMatrixElement(PDT::Spin2,PDT::Spin0,PDT::Spin0)));
+    }
+    if(abs(type_[imode()])==3) {
+      ME(new_ptr(TwoBodyDecayMatrixElement(PDT::Spin2,PDT::Spin1,PDT::Spin0)));
+    }
   }
   // stuff for incoming particle
   if(meopt==Initialize) {
     if(abs(type_[imode()])==1) {
       rho_ = RhoDMatrix(PDT::Spin1);
       Helicity::VectorWaveFunction::calculateWaveFunctions(vecIn_,rho_,const_ptr_cast<tPPtr>(&part),
 							   Helicity::incoming,false);
     }
+    else if(abs(type_[imode()])==2 || abs(type_[imode()])==3) {
+      rho_ = RhoDMatrix(PDT::Spin2);
+      Helicity::TensorWaveFunction::calculateWaveFunctions(tensorIn_,rho_,const_ptr_cast<tPPtr>(&part),
+                 Helicity::incoming,false);
+    }
     else {
       cerr << part << "\n";
       cerr << "testing " << imode() << " " << type_[imode()] << "\n";
       assert(false);
     }
   }
   double output(0.);
   // calculate the matrix element
-  // 1- to 0-
+  InvEnergy fact;
+  Energy pcm = Kinematics::pstarTwoBodyDecay(part.mass(),momenta[0].mass(),
+            momenta[1].mass()); //test subject
+
+  double test(0.), ratio(0.);
+  // HeavyVectorMeson to PScalarMeson + PScalarMeson
   if(abs(type_[imode()])==1) {
-    InvEnergy fact = -2.*g_/fPi_*sqrt(momenta[0].mass()/part.mass());
+    fact = -2.*g_/fPi_*sqrt(momenta[0].mass()/part.mass());
     if(abs(outgoing[1]->id())==111) {
-      fact *= type_[imode()]>0 ? sqrt(0.5) : 0.5*deltaEta_*sqrt(0.5) ;
+      fact *= type_[imode()]>0 ? 0.5 : 0.125*deltaEta_*sqrt(0.5) ;
     }
-    for(unsigned int ix=0;ix<3;++ix) { 
+    for(unsigned int ix=0;ix<3;++ix) {
       (*ME())(ix,0,0) = fact*(vecIn_[ix]*momenta[1]);
     }
-    output = ME()->contract(rho_).real();
     // analytic test of the answer
-    Energy pcm = Kinematics::pstarTwoBodyDecay(part.mass(),momenta[0].mass(),
-					       momenta[1].mass());
-    double test = 4.*sqr(g_)*momenta[0].mass()*sqr(pcm)/3./sqr(fPi_)/part.mass();
+    test = 4.*sqr(g_)*momenta[0].mass()*sqr(pcm)/3./sqr(fPi_)/part.mass();
     if(abs(outgoing[1]->id())==111) {
-      test *= type_[imode()]>0 ? 0.5 : 0.125*sqr(deltaEta_) ;
+      test *= type_[imode()]>0 ? 0.5 : 0.5*deltaEta_*sqrt(0.5) ;
     }
-    double ratio = (output-test)/(output+test);
-    generator()->log() << "testing matrix element for " << part.PDGName() << " -> " 
-		       << outgoing[0]->PDGName() << " " << outgoing[1]->PDGName() << " " 
-		       << output << " " << test << " " << ratio << endl;
   }
-  // for(unsigned int ix=0;ix<5;++ix) {
-  //   (*ME())(ix,0,0) = Complex(_coupling[imode()]/part.mass()*
-  // 			      ((_tensors[ix]*momenta[1])*momenta[0]));
-  // }
-  // double output = ME()->contract(_rho).real();
-  // // test of the answer
-  // // Energy pcm = Kinematics::pstarTwoBodyDecay(part.mass(),momenta[0].mass(),
-  // // 					     momenta[1].mass());
-  // // double test = Energy4(pow<4,1>(2*pcm))*sqr( _coupling[imode()]/part.mass())/120.;
-  // // cout << "testing matrix element for " << part.PDGName() << " -> " 
-  // //      << outgoing[0]->PDGName() << " " << outgoing[1]->PDGName() << " " 
-  // //      << output << " " << test << " " << (output-test)/(output+test) << endl;
+  // HeavyTensorMeson to PScalarMeson + PScalarMeson
+  else if(abs(type_[imode()])==2) {
+    fact = -2.*h_/fPi_*sqrt(momenta[0].mass()/part.mass());
+    if(abs(outgoing[1]->id())==111) {
+      fact *= sqrt(0.5); // ??
+    }
+    for(unsigned int ix=0;ix<5;++ix) {
+      (*ME())(ix,0,0) = (fact/Lambda_)*((tensorIn_[ix]*momenta[1])*momenta[0]);
+    }
+    // analytic test of the answer
+    test = 8.*sqr(h_)*momenta[0].mass()*sqr(sqr(pcm))/10./sqr(fPi_)/sqr(Lambda_)/part.mass();
+    if(abs(outgoing[1]->id())==111) {
+      test *= sqrt(0.5);
+    }
+  }
+  // HeavyTensorMeson to VectorMeson + PScalarMeson
+  else if(abs(type_[imode()])==3) {
+    //get the polarization vectors
+    vecOut_.resize(3);
+    vecOut_={
+        HelicityFunctions::polarizationVector(momenta[0],0,Helicity::outgoing),
+        HelicityFunctions::polarizationVector(momenta[0],1,Helicity::outgoing),
+        HelicityFunctions::polarizationVector(momenta[0],2,Helicity::outgoing)};
+    fact = -2.*h_/fPi_*sqrt(momenta[0].mass()/part.mass());
+    if(abs(outgoing[1]->id())==111) {
+      fact *= sqrt(0.5); // ??
+    }
+    for(unsigned int ix=0;ix<5;++ix) {
+      for(unsigned int iy=0;iy<3;++iy) {
+        if(iy==1) (*ME())(ix,iy,0)=0.;
+        else{
+          LorentzVector<complex<InvEnergy> > vtemp =
+                  (fact/sqr(Lambda_))*epsilon(momenta[0],vecOut_[iy],momenta[1]); ///!!!
+          (*ME())(ix,iy,0) = (momenta[1]*tensorIn_[ix]).dot(vtemp);
+        }
+      }
+    }
+    // analytic test of the answer
+    test = 4.*sqr(h_)*momenta[0].mass()*sqr(sqr(pcm))/5./sqr(fPi_)/sqr(Lambda_)/part.mass(); ///!!!
+    if(abs(outgoing[1]->id())==111) {
+      test *= sqrt(0.5);
+    }
+  }
+  else {
+    assert(false);
+  }
+  output = ME()->contract(rho_).real();
+  //testing
+  ratio = (output-test)/(output+test);
+  generator()->log() << "testing matrix element for " << part.PDGName() << " -> "
+        << outgoing[0]->PDGName() << " " << outgoing[1]->PDGName() << " "
+        << output << " " << test << " " << ratio << endl;
+
   // return the answer
   return output;
 }
 
 bool HQETStrongDecayer::twoBodyMEcode(const DecayMode & dm,int & mecode,
 					       double & coupling) const {
   // int imode(-1);
   // int id(dm.parent()->id());
   // int idbar = dm.parent()->CC() ? dm.parent()->CC()->id() : id;
   // ParticleMSet::const_iterator pit(dm.products().begin());
   // int id1((**pit).id());
   // int id1bar = (**pit).CC() ? (**pit).CC()->id() : id1;
   // ++pit;
   // int id2((**pit).id());
   // int id2bar = (**pit).CC() ? (**pit).CC()->id() : id2;
   // unsigned int ix(0); bool order(false);
   // do {
   //   if(id   ==incoming_[ix]) {
   //     if(id1==outgoingH_[ix]&&id2==outgoingL_[ix]) {
   // 	imode=ix;
   // 	order=true;
   //     }
   //     if(id2==outgoingH_[ix]&&id1==outgoingL_[ix]) {
   // 	imode=ix;
   // 	order=false;
   //     }
   //   }
   //   if(idbar==incoming_[ix]&&imode<0) {
   //     if(id1bar==outgoingH_[ix]&&id2bar==outgoingL_[ix]) {
   // 	imode=ix;
   // 	order=true;
   //     }
   //     if(id2bar==outgoingH_[ix]&&id1bar==outgoingL_[ix]) {
   // 	imode=ix;
   // 	order=false;
   //     }
   //   }
   //   ++ix;
   // }
   // while(ix<incoming_.size()&&imode<0);
   // coupling=_coupling[imode]*dm.parent()->mass();
   // mecode=7;
   // return order;
 }
 
 void HQETStrongDecayer::dataBaseOutput(ofstream & output,
 				       bool header) const {
   if(header) output << "update decayers set parameters=\"";
   // parameters for the DecayIntegrator base class
   DecayIntegrator::dataBaseOutput(output,false);
   // the rest of the parameters
   // couplings
   output << "newdef " << name() << ":fPi " << fPi_/MeV << "\n";
   output << "newdef " << name() << ":g   " << g_      << "\n";
-  if(header) output << "\n\" where BINARY ThePEGName=\"" 
+  if(header) output << "\n\" where BINARY ThePEGName=\""
 		    << fullName() << "\";" << endl;
   for(unsigned int ix=0;ix<incoming_.size();++ix) {
     output << "newdef " << name() << ":MaxWeight " << ix << " " << maxWeight_[ix] << "\n";
   }
 }
diff --git a/Decay/HeavyMeson/HQETStrongDecayer.h b/Decay/HeavyMeson/HQETStrongDecayer.h
--- a/Decay/HeavyMeson/HQETStrongDecayer.h
+++ b/Decay/HeavyMeson/HQETStrongDecayer.h
@@ -1,217 +1,239 @@
 // -*- C++ -*-
 #ifndef Herwig_HQETStrongDecayer_H
 #define Herwig_HQETStrongDecayer_H
 //
 // This is the declaration of the HQETStrongDecayer class.
 //
 
 #include "Herwig/Decay/DecayIntegrator.h"
 #include "Herwig/Decay/PhaseSpaceMode.h"
 #include "ThePEG/Helicity/LorentzPolarizationVector.h"
 
 namespace Herwig {
 
 using namespace ThePEG;
 
 /** \ingroup Decay
  * The HQETStrongDecayer class implements the strong decays of excited heavy, i.e. bottom and charm mesons, using
  * heavy quark effective theory results
  *
  * @see \ref HQETStrongDecayerInterfaces "The interfaces"
  * defined for HQETStrongDecayer.
  */
 class HQETStrongDecayer: public DecayIntegrator {
 
 public:
 
   /**
    * The default constructor.
    */
   HQETStrongDecayer();
 
   /**
    * Which of the possible decays is required
    * @param cc Is this mode the charge conjugate
    * @param parent The decaying particle
    * @param children The decay products
    */
-  virtual int modeNumber(bool & cc, tcPDPtr parent, 
+  virtual int modeNumber(bool & cc, tcPDPtr parent,
 			 const tPDVector & children) const;
 
   /**
    * Return the matrix element squared for a given mode and phase-space channel.
-   * @param ichan The channel we are calculating the matrix element for. 
+   * @param ichan The channel we are calculating the matrix element for.
    * @param part The decaying Particle.
    * @param outgoing The particles produced in the decay
    * @param momenta  The momenta of the particles produced in the decay
    * @param meopt Option for the calculation of the matrix element
    * @return The matrix element squared for the phase-space configuration.
    */
   double me2(const int ichan,const Particle & part,
 	     const tPDVector & outgoing,
 	     const vector<Lorentz5Momentum> & momenta,
 	     MEOption meopt) const;
 
   /**
    *   Construct the SpinInfos for the particles produced in the decay
    */
   virtual void constructSpinInfo(const Particle & part,
 				 ParticleVector outgoing) const;
 
   /**
    * Specify the \f$1\to2\f$ matrix element to be used in the running width calculation.
    * @param dm The DecayMode
    * @param mecode The code for the matrix element as described
    *               in the GenericWidthGenerator class, in this case 7.
    * @param coupling The coupling for the matrix element.
    * @return True or False if this mode can be handled.
    */
   bool twoBodyMEcode(const DecayMode & dm, int & mecode, double & coupling) const;
 
   /**
    * Output the setup information for the particle database
    * @param os The stream to output the information to
    * @param header Whether or not to output the information for MySQL
    */
   virtual void dataBaseOutput(ofstream & os,bool header) const;
-  
+
 public:
 
   /** @name Functions used by the persistent I/O system. */
   //@{
   /**
    * Function used to write out object persistently.
    * @param os the persistent output stream written to.
    */
   void persistentOutput(PersistentOStream & os) const;
 
   /**
    * Function used to read in object persistently.
    * @param is the persistent input stream read from.
    * @param version the version number of the object when written.
    */
   void persistentInput(PersistentIStream & is, int version);
   //@}
 
   /**
    * The standard Init function used to initialize the interfaces.
    * Called exactly once for each class by the class description system
    * before the main function starts or
    * when this class is dynamically loaded.
    */
   static void Init();
 
 protected:
 
   /** @name Clone Methods. */
   //@{
   /**
    * Make a simple clone of this object.
    * @return a pointer to the new object.
    */
   virtual IBPtr clone() const;
 
   /** Make a clone of this object, possibly modifying the cloned object
    * to make it sane.
    * @return a pointer to the new object.
    */
   virtual IBPtr fullclone() const;
   //@}
-  
+
 protected:
-  
+
   /** @name Standard Interfaced functions. */
   //@{
   /**
    * Initialize this object after the setup phase before saving and
    * EventGenerator to disk.
    * @throws InitException if object could not be initialized properly.
    */
   virtual void doinit();
 
   /**
    * Initialize this object to the begining of the run phase.
    */
   virtual void doinitrun();
   //@}
 
 private:
 
   /**
    * The assignment operator is private and must never be called.
    * In fact, it should not even be implemented.
    */
   HQETStrongDecayer & operator=(const HQETStrongDecayer &) = delete;
 
 private:
 
   /**
    *  Coupings for the decays
    */
   //@{
   /**
    *  Pion decay constant
    */
   Energy fPi_;
 
   /**
    *  Coupling for decays within the \f$(0^-,1^-)\f$ multiplet
    */
   double g_;
 
   /**
-   *   \f$\eta-\pi^0$ mixing for isospin violating decays
+   *  Coupling for decays within the \f$(1^+ ,2^+)\f$ multiplet
+   */
+  double h_;
+
+  /**
+   *   \f$\eta-\pi^0 \f$ mixing for isospin violating decays
    */
   double deltaEta_;
   //@}
 
   /**
+   *   A momentum scale characterising the convergence of the
+   *   derivative expansion. We expect Lambda_ ~ 1 GeV.
+   */
+  Energy Lambda_;
+  //@}
+
+  /**
    *  Particles etc for thye different modes
    */
   //@{
   /**
    * the PDG codes for the incoming particles
    */
   vector<int> incoming_;
 
   /**
    * the PDG codes for outgoing heavy meson
    */
   vector<int> outgoingH_;
 
   /**
    * the PDG codes for outgoing light meson
    */
   vector<int> outgoingL_;
 
   /**
    *  Type of decay
    */
   vector<int> type_;
 
   /**
    * the maximum weight for the decay
    */
   vector<double> maxWeight_;
   //@}
 
   /**
    *  Storage of wavefunctions etx
    */
   //@{
   /**
    *   Storage of the \f$\rho\f$ matrix
    */
   mutable RhoDMatrix rho_;
-  
+
   /**
    *  Storage of polarization vectors of the decaying particle
    */
   mutable vector<Helicity::LorentzPolarizationVector> vecIn_;
 
+  /**
+   *  Storage of polarization tensors of the decaying particle
+   */
+  mutable vector<Helicity::LorentzTensor<double> > tensorIn_;
+
+  /**
+   *  Storage of polarization vectors of the decay product
+   */
+  mutable vector<Helicity::LorentzPolarizationVector> vecOut_;
   //@}
+
 };
 
 }
 
 #endif /* Herwig_HQETStrongDecayer_H */
diff --git a/src/defaults/meson_decays.in b/src/defaults/meson_decays.in
--- a/src/defaults/meson_decays.in
+++ b/src/defaults/meson_decays.in
@@ -1,1882 +1,1882 @@
 # -*- ThePEG-repository -*-
 
 #
 # file containing the decays for the mesons
 #
 #
 # the 1^1S_0 multiplet
 #
 #
 decaymode pi+->mu+,nu_mu; 0.999877 1 /Herwig/Decays/PScalarlnu
 decaymode pi+->e+,nu_e; 0.000123 1 /Herwig/Decays/PScalarlnu
 #
 decaymode pi0->gamma,gamma; 0.987989 1 /Herwig/Decays/DecayME0
 decaymode pi0->gamma,e-,e+; 0.01198 1 /Herwig/Decays/PVff
 decaymode pi0->e-,e+,e-,e+; 3.1E-5 1 /Herwig/Decays/PScalar4f
 #
 decaymode pi-->mu-,nu_mubar; 0.999877 1 /Herwig/Decays/PScalarlnu
 decaymode pi-->e-,nu_ebar; 0.000123 1 /Herwig/Decays/PScalarlnu
 #
 decaymode eta->gamma,gamma; 0.393818 1 /Herwig/Decays/PVV
 decaymode eta->pi0,pi0,pi0; 0.325114 1 /Herwig/Decays/Etapipipi
 decaymode eta->pi+,pi-,pi0; 0.22701 1 /Herwig/Decays/Etapipipi
 decaymode eta->pi+,pi-,gamma; 0.046902 1 /Herwig/Decays/Etapipigamma
 decaymode eta->e-,e+,gamma; 0.006 1 /Herwig/Decays/PVff
 decaymode eta->pi0,gamma,gamma; 0.00044 1 /Herwig/Decays/Etapigammagamma
 decaymode eta->pi+,pi-,e-,e+; 0.0004 1 /Herwig/Decays/Mambo
 decaymode eta->mu-,mu+,gamma; 0.00031 1 /Herwig/Decays/PVff
 decaymode eta->mu-,mu+; 6.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode eta'->pi+,pi-,eta; 0.444682 1 /Herwig/Decays/Etapipipi
 decaymode eta'->pi+,pi-,gamma; 0.29379 1 /Herwig/Decays/Etapipigamma
 decaymode eta'->pi0,pi0,eta; 0.207852 1 /Herwig/Decays/Etapipipi
 decaymode eta'->omega,gamma; 0.0303 1 /Herwig/Decays/PVV
 decaymode eta'->gamma,gamma; 0.0212 1 /Herwig/Decays/PVV
 decaymode eta'->pi0,pi0,pi0; 0.00155 1 /Herwig/Decays/Etapipipi
 decaymode eta'->e-,e+,gamma; 0.000522 1 /Herwig/Decays/PVff
 decaymode eta'->mu-,mu+,gamma; 0.000104 1 /Herwig/Decays/PVff
 #
 decaymode eta_c->g,g; 0.720947 1 /Herwig/Decays/QDecayME0
 decaymode eta_c->eta,pi+,pi-; 0.032667 1 /Herwig/Decays/DecayME0
 decaymode eta_c->eta',pi+,pi-; 0.027333 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K-,K0,pi+; 0.023333 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K+,Kbar0,pi-; 0.023333 1 /Herwig/Decays/DecayME0
 decaymode eta_c->eta,pi0,pi0; 0.016333 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*0,K*bar0,pi+,pi-; 0.015 1 /Herwig/Decays/Mambo
 decaymode eta_c->eta',pi0,pi0; 0.013667 1 /Herwig/Decays/DecayME0
 decaymode eta_c->rho+,rho-; 0.013333 1 /Herwig/Decays/PVV
 decaymode eta_c->K0,Kbar0,pi0; 0.011667 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K+,K-,pi0; 0.011667 1 /Herwig/Decays/DecayME0
 decaymode eta_c->f_2,f_2; 0.01 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*-,K0,pi+; 0.00847 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*bar0,K+,pi-; 0.00847 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*0,K-,pi+; 0.00847 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*+,Kbar0,pi-; 0.00847 1 /Herwig/Decays/DecayME0
 decaymode eta_c->rho0,rho0; 0.006667 1 /Herwig/Decays/PVV
 decaymode eta_c->pi+,pi-,pi+,pi-; 0.005333 1 /Herwig/Decays/Mambo
 decaymode eta_c->K*0,K*bar0; 0.0046 1 /Herwig/Decays/PVV
 decaymode eta_c->K*+,K*-; 0.0046 1 /Herwig/Decays/PVV
 decaymode eta_c->K*-,K+,pi0; 0.00423 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*bar0,K0,pi0; 0.00423 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*0,Kbar0,pi0; 0.00423 1 /Herwig/Decays/DecayME0
 decaymode eta_c->K*+,K-,pi0; 0.00423 1 /Herwig/Decays/DecayME0
 decaymode eta_c->phi,phi; 0.0027 1 /Herwig/Decays/PVV
 decaymode eta_c->phi,K0,Kbar0; 0.00157 1 /Herwig/Decays/DecayME0
 decaymode eta_c->phi,K+,K-; 0.00157 1 /Herwig/Decays/DecayME0
 decaymode eta_c->n0,nbar0; 0.0013 1 /Herwig/Decays/DecayME0
 decaymode eta_c->p+,pbar-; 0.0013 1 /Herwig/Decays/DecayME0
 decaymode eta_c->gamma,gamma; 0.00028 1 /Herwig/Decays/PVV
 #
 decaymode eta_b->g,g; 1 1 /Herwig/Decays/QDecayME0
 #
 decaymode K+->mu+,nu_mu; 0.634277 1 /Herwig/Decays/PScalarlnu
 decaymode K+->pi+,pi0; 0.209159 1 /Herwig/Decays/DecayME0
 decaymode K+->pi+,pi+,pi-; 0.0559 1 /Herwig/Decays/DecayME0
 decaymode K+->pi0,e+,nu_e; 0.0498 1 /Herwig/Decays/DecayME101
 decaymode K+->pi0,mu+,nu_mu; 0.0332 1 /Herwig/Decays/DecayME101
 decaymode K+->pi+,pi0,pi0; 0.01757 1 /Herwig/Decays/DecayME0
 decaymode K+->pi+,pi-,e+,nu_e; 4.1E-5 1 /Herwig/Decays/Mambo
 decaymode K+->pi0,pi0,e+,nu_e; 2.2E-5 1 /Herwig/Decays/Mambo
 decaymode K+->e+,nu_e; 1.6E-5 1 /Herwig/Decays/PScalarlnu
 decaymode K+->pi+,pi-,mu+,nu_mu; 1.4E-5 1 /Herwig/Decays/Mambo
 decaymode K+->pi+,gamma,gamma; 1.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode K0->K_L0; 0.5 1 /Herwig/Decays/DecayME0
 decaymode K0->K_S0; 0.5 1 /Herwig/Decays/DecayME0
 #
 decaymode K-->mu-,nu_mubar; 0.634277 1 /Herwig/Decays/PScalarlnu
 decaymode K-->pi-,pi0; 0.209159 1 /Herwig/Decays/DecayME0
 decaymode K-->pi-,pi-,pi+; 0.0559 1 /Herwig/Decays/DecayME0
 decaymode K-->pi0,e-,nu_ebar; 0.0498 1 /Herwig/Decays/DecayME101
 decaymode K-->pi0,mu-,nu_mubar; 0.0332 1 /Herwig/Decays/DecayME101
 decaymode K-->pi-,pi0,pi0; 0.01757 1 /Herwig/Decays/DecayME0
 decaymode K-->pi-,pi+,e-,nu_ebar; 4.1E-5 1 /Herwig/Decays/Mambo
 decaymode K-->pi0,pi0,e-,nu_ebar; 2.2E-5 1 /Herwig/Decays/Mambo
 decaymode K-->e-,nu_ebar; 1.6E-5 1 /Herwig/Decays/PScalarlnu
 decaymode K-->pi-,pi+,mu-,nu_mubar; 1.4E-5 1 /Herwig/Decays/Mambo
 decaymode K-->pi-,gamma,gamma; 1.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode Kbar0->K_L0; 0.5 1 /Herwig/Decays/DecayME0
 decaymode Kbar0->K_S0; 0.5 1 /Herwig/Decays/DecayME0
 #
 decaymode B_c+->B_s*0,rho+; 0.202 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s0,pi+; 0.164 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s0,rho+; 0.072 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s*0,pi+; 0.065 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s*0,mu+,nu_mu; 0.0506 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B_s*0,e+,nu_e; 0.0506 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->c,sbar; 0.049 1 /Herwig/Decays/QDecayME0
 decaymode B_c+->B_s0,mu+,nu_mu; 0.0403 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B_s0,e+,nu_e; 0.0403 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B*0,rho+; 0.0257 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B+,Kbar0; 0.0198 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,mu+,nu_mu; 0.019 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->Jpsi,e+,nu_e; 0.019 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B*+,K*bar0; 0.0167 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->tau+,nu_tau; 0.016 1 /Herwig/Decays/PScalarlnu
 decaymode B_c+->B*+,Kbar0; 0.016 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->dbar=u,cbar=c; 0.011134 1 /Herwig/Decays/WeakPartonic100
 decaymode B_c+->B0,pi+; 0.0106 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s0,K+; 0.0106 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B0,rho+; 0.0096 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B*0,pi+; 0.0095 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,mu+,nu_mu; 0.0075 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c,e+,nu_e; 0.0075 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->Jpsi,D_s*+; 0.0067 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B*0,mu+,nu_mu; 0.0058 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B*0,e+,nu_e; 0.0058 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->Jpsi,tau+,nu_tau; 0.0048 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B+,K*bar0; 0.0043 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,rho+; 0.0042 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,rho+; 0.004 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B_s*0,K+; 0.0037 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B0,mu+,nu_mu; 0.0034 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B0,e+,nu_e; 0.0034 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c,D_s+; 0.0028 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,D_s*+; 0.0027 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,tau+,nu_tau; 0.0023 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c,pi+; 0.002 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,D_s+; 0.0017 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,pi+; 0.0013 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->psi(2S),mu+,nu_mu; 0.00094 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->psi(2S),e+,nu_e; 0.00094 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->B*+,rho0; 0.0009 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B0,K+; 0.0007 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B*0,K*+; 0.00058 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B*0,K+; 0.00055 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B+,pi0; 0.00037 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B+,rho0; 0.00034 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D*+,D*bar0; 0.00033 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B*+,pi0; 0.00033 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,D*+; 0.00028 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,K*+; 0.00022 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,K*+; 0.0002 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c(2S),mu+,nu_mu; 0.0002 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c(2S),e+,nu_e; 0.0002 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D*0,mu+,nu_mu; 0.00018 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D*0,e+,nu_e; 0.00018 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c,D+; 0.00015 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->B0,K*+; 0.00015 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,K+; 0.00013 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,K+; 0.00011 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->eta_c,D*+; 0.0001 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->Jpsi,D+; 9.0E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D*0,tau+,nu_tau; 8.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->psi(2S),tau+,nu_tau; 8.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D+,D*bar0; 7.5E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D+,Dbar0; 5.3E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D*+,Dbar0; 4.9E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D0,mu+,nu_mu; 4.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D0,e+,nu_e; 4.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D_s*+,D*bar0; 2.6E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D0,tau+,nu_tau; 2.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->eta_c(2S),tau+,nu_tau; 1.6E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c+->D_s+,D*bar0; 7.0E-6 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D_s+,Dbar0; 5.0E-6 1 /Herwig/Decays/BcHadronic
 decaymode B_c+->D_s*+,Dbar0; 5.0E-6 1 /Herwig/Decays/BcHadronic
 #
 decaymode B_c-->B_s*bar0,rho-; 0.202 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B_sbar0,pi-; 0.164 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B_sbar0,rho-; 0.072 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B_s*bar0,pi-; 0.065 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B_s*bar0,e-,nu_ebar; 0.0506 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B_s*bar0,mu-,nu_mubar; 0.0506 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->cbar,s; 0.049 1 /Herwig/Decays/QDecayME0
 decaymode B_c-->B_sbar0,e-,nu_ebar; 0.0403 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B_sbar0,mu-,nu_mubar; 0.0403 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B*bar0,rho-; 0.0257 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B-,K0; 0.0198 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,e-,nu_ebar; 0.019 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Jpsi,mu-,nu_mubar; 0.019 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B*-,K*0; 0.0167 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*-,K0; 0.016 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->tau-,nu_taubar; 0.016 1 /Herwig/Decays/PScalarlnu
 decaymode B_c-->d=ubar,c=cbar; 0.011134 1 /Herwig/Decays/WeakPartonic100
 decaymode B_c-->B_sbar0,K-; 0.0106 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Bbar0,pi-; 0.0106 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Bbar0,rho-; 0.0096 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*bar0,pi-; 0.0095 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,e-,nu_ebar; 0.0075 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->eta_c,mu-,nu_mubar; 0.0075 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Jpsi,D_s*-; 0.0067 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*bar0,e-,nu_ebar; 0.0058 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B*bar0,mu-,nu_mubar; 0.0058 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Jpsi,tau-,nu_taubar; 0.0048 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B-,K*0; 0.0043 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,rho-; 0.0042 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,rho-; 0.004 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B_s*bar0,K-; 0.0037 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Bbar0,e-,nu_ebar; 0.0034 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Bbar0,mu-,nu_mubar; 0.0034 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->eta_c,D_s-; 0.0028 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,D_s*-; 0.0027 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,tau-,nu_taubar; 0.0023 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->eta_c,pi-; 0.002 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,D_s-; 0.0017 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,pi-; 0.0013 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->psi(2S),e-,nu_ebar; 0.00094 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->psi(2S),mu-,nu_mubar; 0.00094 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->B*-,rho0; 0.0009 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Bbar0,K-; 0.0007 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*bar0,K*-; 0.00058 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*bar0,K-; 0.00055 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B-,pi0; 0.00037 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B-,rho0; 0.00034 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->B*-,pi0; 0.00033 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D*-,D*0; 0.00033 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,D*-; 0.00028 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,K*-; 0.00022 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,K*-; 0.0002 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c(2S),e-,nu_ebar; 0.0002 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->eta_c(2S),mu-,nu_mubar; 0.0002 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->D*bar0,e-,nu_ebar; 0.00018 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->D*bar0,mu-,nu_mubar; 0.00018 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Bbar0,K*-; 0.00015 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,D-; 0.00015 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,K-; 0.00013 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,K-; 0.00011 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->eta_c,D*-; 0.0001 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Jpsi,D-; 9.0E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D*bar0,tau-,nu_taubar; 8.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->psi(2S),tau-,nu_taubar; 8.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->D-,D*0; 7.5E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D-,D0; 5.3E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D*-,D0; 4.9E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Dbar0,e-,nu_ebar; 4.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->Dbar0,mu-,nu_mubar; 4.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->D_s*-,D*0; 2.6E-5 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->Dbar0,tau-,nu_taubar; 2.0E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->eta_c(2S),tau-,nu_taubar; 1.6E-5 1 /Herwig/Decays/KiselevBc
 decaymode B_c-->D_s-,D*0; 7.0E-6 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D_s*-,D0; 5.0E-6 1 /Herwig/Decays/BcHadronic
 decaymode B_c-->D_s-,D0; 5.0E-6 1 /Herwig/Decays/BcHadronic
 #
 # the 1^3S_1 multiplet
 #
 #
 decaymode rho+->pi+,pi0; 0.99955 1 /Herwig/Decays/Vector2Meson
 decaymode rho+->pi+,gamma; 0.00045 1 /Herwig/Decays/VectorVP
 #
 decaymode rho0->pi+,pi-; 0.998847 1 /Herwig/Decays/Vector2Meson
 decaymode rho0->pi0,gamma; 0.0006 1 /Herwig/Decays/VectorVP
 decaymode rho0->eta,gamma; 0.000295 1 /Herwig/Decays/VectorVP
 decaymode rho0->pi+,pi-,pi0; 0.000101 1 /Herwig/Decays/DecayME0
 decaymode rho0->e-,e+; 4.7E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode rho0->mu-,mu+; 4.7E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode rho0->pi0,pi0,gamma; 4.5E-5 1 /Herwig/Decays/DecayME0
 decaymode rho0->pi+,pi-,pi+,pi-; 1.8E-5 1 /Herwig/Decays/Mambo
 #
 decaymode rho-->pi-,pi0; 0.99955 1 /Herwig/Decays/Vector2Meson
 decaymode rho-->pi-,gamma; 0.00045 1 /Herwig/Decays/VectorVP
 #
 decaymode omega->pi+,pi-,pi0; 0.891174 1 /Herwig/Decays/Vector3Pion
 decaymode omega->pi0,gamma; 0.09025 1 /Herwig/Decays/VectorVP
 decaymode omega->pi+,pi-; 0.017 1 /Herwig/Decays/Vector2Meson
 decaymode omega->pi0,e-,e+; 0.000797 1 /Herwig/Decays/VectorVPff
 decaymode omega->eta,gamma; 0.00049 1 /Herwig/Decays/VectorVP
 decaymode omega->pi0,mu-,mu+; 7.8E-5 1 /Herwig/Decays/VectorVPff
 decaymode omega->e-,e+; 7.2E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode omega->mu-,mu+; 7.2E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode omega->pi0,pi0,gamma; 6.7E-5 1 /Herwig/Decays/DecayME0
 #
 decaymode phi->K+,K-; 0.491536 1 /Herwig/Decays/Vector2Meson
 decaymode phi->K_L0,K_S0; 0.339679 1 /Herwig/Decays/Vector2Meson
 decaymode phi->pi+,pi-,pi0; 0.152856 1 /Herwig/Decays/Vector3Pion
 decaymode phi->eta,gamma; 0.01301 1 /Herwig/Decays/VectorVP
 decaymode phi->pi0,gamma; 0.00125 1 /Herwig/Decays/VectorVP
 decaymode phi->f_0,gamma; 0.00044 1 /Herwig/Decays/VVSDecayer
 decaymode phi->e-,e+; 0.000297 1 /Herwig/Decays/Vector2Leptons
 decaymode phi->mu-,mu+; 0.000296 1 /Herwig/Decays/Vector2Leptons
 decaymode phi->pi+,pi-; 0.000114 1 /Herwig/Decays/Vector2Meson
 decaymode phi->eta,e-,e+; 0.000114 1 /Herwig/Decays/VectorVPff
 decaymode phi->pi0,pi0,gamma; 0.000109 1 /Herwig/Decays/DecayME0
 decaymode phi->eta,pi0,gamma; 8.3E-5 1 /Herwig/Decays/DecayME0
 decaymode phi->a_00,gamma; 7.6E-5 1 /Herwig/Decays/VVSDecayer
 decaymode phi->eta',gamma; 6.2E-5 1 /Herwig/Decays/VectorVP
 decaymode phi->omega,pi0; 5.2E-5 1 /Herwig/Decays/VectorVP
 decaymode phi->pi0,e-,e+; 1.3E-5 1 /Herwig/Decays/VectorVPff
 decaymode phi->eta,mu-,mu+; 6.0E-6 1 /Herwig/Decays/VectorVPff
 decaymode phi->pi+,pi-,pi+,pi-; 4.0E-6 1 /Herwig/Decays/Mambo
 decaymode phi->pi0,mu-,mu+; 3.0E-6 1 /Herwig/Decays/VectorVPff
 #
 decaymode Jpsi->g,g,g; 0.486225 1 /Herwig/Decays/QDecayME130
 decaymode Jpsi->u,ubar; 0.0791 1 /Herwig/Decays/QDecayME0
 decaymode Jpsi->e-,e+; 0.05935 1 /Herwig/Decays/Vector2Leptons
 decaymode Jpsi->mu-,mu+; 0.05935 1 /Herwig/Decays/Vector2Leptons
 decaymode Jpsi->g,g,gamma; 0.05 1 /Herwig/Decays/QDecayME130
 decaymode Jpsi->d,dbar; 0.0198 1 /Herwig/Decays/QDecayME0
 decaymode Jpsi->s,sbar; 0.0198 1 /Herwig/Decays/QDecayME0
 decaymode Jpsi->gamma,eta_c; 0.013 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->omega,pi+,pi-,pi+,pi-; 0.0085 1 /Herwig/Decays/Mambo
 decaymode Jpsi->rho-,pi+; 0.00563 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->rho0,pi0; 0.00563 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->rho+,pi-; 0.00563 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->gamma,pi+,pi-,pi0,pi0; 0.005 1 /Herwig/Decays/Mambo
 decaymode Jpsi->omega,pi+,pi-; 0.00475 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,eta'; 0.00471 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->n0,nbar0,pi+,pi-; 0.0045 1 /Herwig/Decays/Mambo
 decaymode Jpsi->p+,pbar-,pi+,pi-; 0.0045 1 /Herwig/Decays/Mambo
 decaymode Jpsi->f_2,omega; 0.0043 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->a_20,rho0; 0.003634 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->a_2-,rho+; 0.003633 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->a_2+,rho-; 0.003633 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,pi+,pi-,pi0; 0.0036 1 /Herwig/Decays/Mambo
 decaymode Jpsi->K*-,K*_2+; 0.00335 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K*bar0,K*_20; 0.00335 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K*0,K*_2bar0; 0.00335 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K*+,K*_2-; 0.00335 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pi+,pi-,pi0; 0.0033 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,eta,pi+,pi-; 0.0031 1 /Herwig/Decays/Mambo
 decaymode Jpsi->gamma,rho+,rho-; 0.003 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,eta(1475); 0.0028 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K-,K*+; 0.0023 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->Kbar0,K*0; 0.0023 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->K0,K*bar0; 0.0023 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->K+,K*-; 0.0023 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->p+,pbar-,pi+,pi-,pi0; 0.0023 1 /Herwig/Decays/Mambo
 decaymode Jpsi->omega,pi0,pi0; 0.00217 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,nbar0; 0.00217 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->p+,pbar-; 0.00217 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->gamma,K+,K-,pi+,pi-; 0.0021 1 /Herwig/Decays/Mambo
 decaymode Jpsi->n0,nbar0,eta; 0.00209 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,pbar-,eta; 0.00209 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pbar-,n0,pi+; 0.00204 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,nbar0,pi-; 0.00204 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,K*0,K*bar0; 0.002 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,K*+,K*-; 0.002 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K'_1-,K+; 0.0019 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->K'_1bar0,K0; 0.0019 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->K'_10,Kbar0; 0.0019 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->K'_1+,K-; 0.0019 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->omega,eta; 0.00174 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->phi,pi+,pi-,pi+,pi-; 0.00166 1 /Herwig/Decays/Mambo
 decaymode Jpsi->b_1-,pi+; 0.00162 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->b_10,pi0; 0.00162 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->b_1+,pi-; 0.00162 1 /Herwig/Decays/VAPDecayer
 decaymode Jpsi->gamma,eta,pi0,pi0; 0.0016 1 /Herwig/Decays/Mambo
 decaymode Jpsi->gamma,omega,omega; 0.00159 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Lambda0,Lambdabar0; 0.00154 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->gamma,rho0,rho0; 0.0015 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pi+,pi-; 0.00147 1 /Herwig/Decays/Vector2Meson
 decaymode Jpsi->gamma,f_2; 0.00138 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,f_0(1710); 0.00134 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->omega,K*-,K+; 0.001325 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,K*bar0,K0; 0.001325 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,K*0,Kbar0; 0.001325 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,K*+,K-; 0.001325 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma-,Sigmabar+; 0.00131 1 /Herwig/Decays/Vector2Leptons
 decaymode Jpsi->Sigma0,Sigmabar0; 0.00131 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->Sigma+,Sigmabar-; 0.00131 1 /Herwig/Decays/Vector2Leptons
 decaymode Jpsi->n0,nbar0,omega; 0.0013 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,pbar-,omega; 0.0013 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Delta-,Deltabar+; 0.0011 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Delta0,Deltabar0; 0.0011 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Delta+,Deltabar-; 0.0011 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Delta++,Deltabar--; 0.0011 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pbar-,K0,Sigma*+; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pbar-,p+,pi0; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->nbar0,K+,Sigma*-; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,nbar0,pi0; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,K-,Sigma*bar+; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,Kbar0,Sigma*bar-; 0.00102 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,pi+,pi-,pi+,pi-; 0.001 1 /Herwig/Decays/Mambo
 decaymode Jpsi->gamma,eta; 0.0009657 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->eta,e+,e-; 1.43e-5 1   /Herwig/Decays/VectorVPff
 decaymode Jpsi->gamma,f_2,f_2; 0.00095 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->eta',n0,nbar0; 0.0009 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->eta',p+,pbar-; 0.0009 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Xi-,Xibar+; 0.0009 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->Xi0,Xibar0; 0.0009 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->pbar-,K+,Lambda0; 0.00089 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->nbar0,K0,Lambda0; 0.00089 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,Kbar0,Lambdabar0; 0.00089 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,K-,Lambdabar0; 0.00089 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,f'_2; 0.0008 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,f'_1; 0.00079 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,pi+,pi-; 0.00078 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->f_0(1710),omega; 0.00076 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->phi,eta; 0.00074 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->phi,K0,Kbar0; 0.00074 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,K+,K-; 0.00074 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,K0,Kbar0; 0.00071 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,K+,K-; 0.00071 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,f'_1; 0.00068 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,f_1; 0.00061 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->pbar-,K0,Sigma+; 0.00058 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->nbar0,K+,Sigma-; 0.00058 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,K-,Sigmabar+; 0.00058 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,Kbar0,Sigmabar-; 0.00058 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,f_0(1500); 0.00057 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->f_0(1710),phi; 0.00057 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->Sigmabar-,pi+,Lambda0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigmabar0,pi0,Lambda0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigmabar+,pi-,Lambda0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma-,pi+,Lambdabar0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma0,pi0,Lambdabar0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma+,pi-,Lambdabar0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*-,Sigma*bar+; 0.000515 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->Sigma*0,Sigma*bar0; 0.000515 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->Sigma*+,Sigma*bar-; 0.000515 1 /Herwig/Decays/Vector2Baryons
 decaymode Jpsi->pbar-,K+,Sigma*0; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->nbar0,K0,Sigma*0; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K-,K*+,phi; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Kbar0,K*0,phi; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K0,K*bar0,phi; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K+,K*-,phi; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,Kbar0,Sigma*bar0; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,K-,Sigma*bar0; 0.00051 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,f'_2; 0.00045 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,eta(1405); 0.00045 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->omega,pi0; 0.00045 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->Xi*bar0,Xi0; 0.000446 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Xi*bar+,Xi-; 0.000446 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Xi*-,Xibar+; 0.000446 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Xi*0,Xibar0; 0.000446 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,phi,phi; 0.0004 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,eta'; 0.0004 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->gamma,n0,nbar0; 0.00038 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,p+,pbar-; 0.00038 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,f_0; 0.00032 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->pbar-,K+,Sigma0; 0.00029 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->nbar0,K0,Sigma0; 0.00029 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->n0,Kbar0,Sigmabar0; 0.00029 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->p+,K-,Sigmabar0; 0.00029 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,f_1; 0.00026 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K-,K+; 0.000237 1 /Herwig/Decays/Vector2Meson
 decaymode Jpsi->Lambda0,Lambdabar0,pi0; 0.00022 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->rho0,eta; 0.000193 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->omega,eta'; 0.000182 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->Sigma*bar-,Sigma+; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*bar0,Sigma0; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*bar+,Sigma-; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*-,Sigmabar+; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*0,Sigmabar0; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->Sigma*+,Sigmabar-; 0.000155 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->K_L0,K_S0; 0.000146 1 /Herwig/Decays/Vector2Meson
 decaymode Jpsi->omega,f_0; 0.00014 1 /Herwig/Decays/VVSDecayer
 decaymode Jpsi->rho0,eta'; 0.000105 1 /Herwig/Decays/VectorVP
 decaymode Jpsi->phi,n0,nbar0; 4.5E-5 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->phi,p+,pbar-; 4.5E-5 1 /Herwig/Decays/DecayME0
 decaymode Jpsi->gamma,pi0; 3.3E-5 1 /Herwig/Decays/VectorVP
 #
 decaymode Upsilon->g,g,g; 0.814282 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon->u,ubar; 0.036 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon->c,cbar; 0.036 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon->e-,e+; 0.0249 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon->mu-,mu+; 0.0249 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon->tau-,tau+; 0.0247 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon->g,g,gamma; 0.02099 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon->d,dbar; 0.009 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon->s,sbar; 0.009 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon->gamma,f_2; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode Upsilon->gamma,pi+,pi-; 6.3E-5 1 /Herwig/Decays/DecayME0
 decaymode Upsilon->gamma,f'_2; 3.7E-5 1 /Herwig/Decays/DecayME0
 decaymode Upsilon->gamma,pi0,pi0; 1.7E-5 1 /Herwig/Decays/DecayME0
 decaymode Upsilon->gamma,K+,K-; 1.1E-5 1 /Herwig/Decays/DecayME0
 #
 decaymode K*+->K0,pi+; 0.666007 1 /Herwig/Decays/Vector2Meson
 decaymode K*+->K+,pi0; 0.333003 1 /Herwig/Decays/Vector2Meson
 decaymode K*+->K+,gamma; 0.00099 1 /Herwig/Decays/VectorVP
 #
 decaymode K*0->K+,pi-; 0.665127 1 /Herwig/Decays/Vector2Meson
 decaymode K*0->K0,pi0; 0.332563 1 /Herwig/Decays/Vector2Meson
 decaymode K*0->K0,gamma; 0.00231 1 /Herwig/Decays/VectorVP
 #
 decaymode K*-->Kbar0,pi-; 0.666007 1 /Herwig/Decays/Vector2Meson
 decaymode K*-->K-,pi0; 0.333003 1 /Herwig/Decays/Vector2Meson
 decaymode K*-->K-,gamma; 0.00099 1 /Herwig/Decays/VectorVP
 #
 decaymode K*bar0->K-,pi+; 0.665127 1 /Herwig/Decays/Vector2Meson
 decaymode K*bar0->Kbar0,pi0; 0.332563 1 /Herwig/Decays/Vector2Meson
 decaymode K*bar0->Kbar0,gamma; 0.00231 1 /Herwig/Decays/VectorVP
 #
 decaymode D*0->D0,pi0; 0.619 1 /Herwig/Decays/HQETMeson
 decaymode D*0->D0,gamma; 0.381 1 /Herwig/Decays/VectorVP
 #
 decaymode D*+->D0,pi+; 0.677 1 /Herwig/Decays/HQETMeson
 decaymode D*+->D+,pi0; 0.307 1 /Herwig/Decays/HQETMeson
 decaymode D*+->D+,gamma; 0.016 1 /Herwig/Decays/VectorVP
 #
 decaymode D*bar0->Dbar0,pi0; 0.619 1 /Herwig/Decays/HQETMeson
 decaymode D*bar0->Dbar0,gamma; 0.381 1 /Herwig/Decays/VectorVP
 #
 decaymode D*-->Dbar0,pi-; 0.677 1 /Herwig/Decays/HQETMeson
 decaymode D*-->D-,pi0; 0.307 1 /Herwig/Decays/HQETMeson
 decaymode D*-->D-,gamma; 0.016 1 /Herwig/Decays/VectorVP
 #
 decaymode D_s*+->D_s+,pi0; 0.942 1 /Herwig/Decays/Vector2Meson
 decaymode D_s*+->D_s+,gamma; 0.058 1 /Herwig/Decays/VectorVP
 #
 decaymode D_s*-->D_s-,gamma; 0.942 1 /Herwig/Decays/VectorVP
 decaymode D_s*-->D_s-,pi0; 0.058 1 /Herwig/Decays/Vector2Meson
 #
 decaymode B*0->B0,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B*+->B+,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B*bar0->Bbar0,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B*-->B-,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B_s*0->B_s0,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B_s*bar0->B_sbar0,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B_c*+->B_c+,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 decaymode B_c*-->B_c-,gamma; 1 1 /Herwig/Decays/VectorVP
 #
 # the 1^1P_1 multiplet
 #
 #
 decaymode b_1+->omega,pi+; 0.9984 1 /Herwig/Decays/PVectorVP
 decaymode b_1+->pi+,gamma; 0.0016 1 /Herwig/Decays/PVectorVP
 #
 decaymode b_10->omega,pi0; 1 1 /Herwig/Decays/PVectorVP
 #
 decaymode b_1-->omega,pi-; 0.9984 1 /Herwig/Decays/PVectorVP
 decaymode b_1-->pi-,gamma; 0.0016 1 /Herwig/Decays/PVectorVP
 #
 decaymode h_1->rho0,pi0; 0.333334 1 /Herwig/Decays/PVectorVP
 decaymode h_1->rho-,pi+; 0.333333 1 /Herwig/Decays/PVectorVP
 decaymode h_1->rho+,pi-; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode h'_1->K*-,K+; 0.25 1 /Herwig/Decays/PVectorVP
 decaymode h'_1->K*bar0,K0; 0.25 1 /Herwig/Decays/PVectorVP
 decaymode h'_1->K*0,Kbar0; 0.25 1 /Herwig/Decays/PVectorVP
 decaymode h'_1->K*+,K-; 0.25 1 /Herwig/Decays/PVectorVP
 #
 decaymode h_c->g,g,g; 0.6 1 /Herwig/Decays/QDecayME130
 decaymode h_c->eta_c,gamma; 0.4 1 /Herwig/Decays/PVectorVP
 #
 decaymode h_b->g,g,g; 0.568 1 /Herwig/Decays/QDecayME130
 decaymode h_b->eta_b,gamma; 0.414 1 /Herwig/Decays/PVectorVP
 decaymode h_b->g,g,gamma; 0.018 1 /Herwig/Decays/QDecayME130
 #
 decaymode K_1+->rho+,K0; 0.298667 1 /Herwig/Decays/PVectorVP
 decaymode K_1+->K*_00,pi+; 0.199111 1 /Herwig/Decays/Vector2Meson
 decaymode K_1+->rho0,K+; 0.149333 1 /Herwig/Decays/PVectorVP
 decaymode K_1+->K*0,pi+; 0.122222 1 /Herwig/Decays/PVectorVP
 decaymode K_1+->K*_0+,pi0; 0.099556 1 /Herwig/Decays/Vector2Meson
 decaymode K_1+->K+,omega; 0.07 1 /Herwig/Decays/PVectorVP
 decaymode K_1+->K*+,pi0; 0.061111 1 /Herwig/Decays/PVectorVP
 #
 decaymode K_10->rho-,K+; 0.298412 1 /Herwig/Decays/PVectorVP
 decaymode K_10->K*_0+,pi-; 0.198941 1 /Herwig/Decays/Vector2Meson
 decaymode K_10->rho0,K0; 0.149206 1 /Herwig/Decays/PVectorVP
 decaymode K_10->K*+,pi-; 0.122009 1 /Herwig/Decays/PVectorVP
 decaymode K_10->K*_00,pi0; 0.09947 1 /Herwig/Decays/Vector2Meson
 decaymode K_10->K0,omega; 0.07 1 /Herwig/Decays/PVectorVP
 decaymode K_10->K*0,pi0; 0.061005 1 /Herwig/Decays/PVectorVP
 decaymode K_10->K0,gamma; 0.000957 1 /Herwig/Decays/PVectorVP
 #
 decaymode K_1-->rho-,Kbar0; 0.298667 1 /Herwig/Decays/PVectorVP
 decaymode K_1-->K*_0bar0,pi-; 0.199111 1 /Herwig/Decays/Vector2Meson
 decaymode K_1-->rho0,K-; 0.149333 1 /Herwig/Decays/PVectorVP
 decaymode K_1-->K*bar0,pi-; 0.122222 1 /Herwig/Decays/PVectorVP
 decaymode K_1-->K*_0-,pi0; 0.099556 1 /Herwig/Decays/Vector2Meson
 decaymode K_1-->K-,omega; 0.07 1 /Herwig/Decays/PVectorVP
 decaymode K_1-->K*-,pi0; 0.061111 1 /Herwig/Decays/PVectorVP
 #
 decaymode K_1bar0->rho+,K-; 0.298412 1 /Herwig/Decays/PVectorVP
 decaymode K_1bar0->K*_0-,pi+; 0.198941 1 /Herwig/Decays/Vector2Meson
 decaymode K_1bar0->rho0,Kbar0; 0.149206 1 /Herwig/Decays/PVectorVP
 decaymode K_1bar0->K*-,pi+; 0.122009 1 /Herwig/Decays/PVectorVP
 decaymode K_1bar0->K*_0bar0,pi0; 0.09947 1 /Herwig/Decays/Vector2Meson
 decaymode K_1bar0->Kbar0,omega; 0.07 1 /Herwig/Decays/PVectorVP
 decaymode K_1bar0->K*bar0,pi0; 0.061005 1 /Herwig/Decays/PVectorVP
 decaymode K_1bar0->Kbar0,gamma; 0.000957 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_10->D*+,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D_10->D*0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_1+->D*0,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D_1+->D*+,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_1bar0->D*-,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D_1bar0->D*bar0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_1-->D*bar0,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D_1-->D*-,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D'_s1+->D_s*+,pi0; 0.656068 1 /Herwig/Decays/PVectorVP
 decaymode D'_s1+->D_s+,gamma; 0.206139 1 /Herwig/Decays/PVectorVP
 decaymode D'_s1+->D_s+,pi+,pi-; 0.091862 1 /Herwig/Decays/DecayME0
 decaymode D'_s1+->D_s+,pi0,pi0; 0.045931 1 /Herwig/Decays/DecayME0
 #
 decaymode D'_s1-->D_s*-,pi0; 0.656068 1 /Herwig/Decays/PVectorVP
 decaymode D'_s1-->D_s-,gamma; 0.206139 1 /Herwig/Decays/PVectorVP
 decaymode D'_s1-->D_s-,pi-,pi+; 0.091862 1 /Herwig/Decays/DecayME0
 decaymode D'_s1-->D_s-,pi0,pi0; 0.045931 1 /Herwig/Decays/DecayME0
 #
 decaymode B_10->B*+,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B_10->B*0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_1+->B*0,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B_1+->B*+,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_1bar0->B*-,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B_1bar0->B*bar0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_1-->B*bar0,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B_1-->B*-,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_s10->B*0,Kbar0; 0.5 1 /Herwig/Decays/PVectorVP
 decaymode B_s10->B*+,K-; 0.5 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_s1bar0->B*-,K+; 0.5 1 /Herwig/Decays/PVectorVP
 decaymode B_s1bar0->B*bar0,K0; 0.5 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_c1+->B_c*+,gamma; 0.82 1 /Herwig/Decays/DecayME0
 decaymode B_c1+->B_c+,gamma; 0.18 1 /Herwig/Decays/PVectorVP
 #
 decaymode B_c1-->B_c*-,gamma; 0.82 1 /Herwig/Decays/DecayME0
 decaymode B_c1-->B_c-,gamma; 0.18 1 /Herwig/Decays/PVectorVP
 #
 # the 1^3P_0 multiplet
 #
 #
 decaymode a'_0+->omega,rho+; 0.827532 1 /Herwig/Decays/SVV
 decaymode a'_0+->pi+,eta; 0.07734 1 /Herwig/Decays/ScalarSS
 decaymode a'_0+->K+,Kbar0; 0.068059 1 /Herwig/Decays/ScalarSS
 decaymode a'_0+->pi+,eta'; 0.027069 1 /Herwig/Decays/ScalarSS
 #
 decaymode a'_00->omega,rho0; 0.827533 1 /Herwig/Decays/SVV
 decaymode a'_00->pi0,eta; 0.07734 1 /Herwig/Decays/ScalarSS
 decaymode a'_00->K-,K+; 0.034029 1 /Herwig/Decays/ScalarSS
 decaymode a'_00->Kbar0,K0; 0.034029 1 /Herwig/Decays/ScalarSS
 decaymode a'_00->pi0,eta'; 0.027069 1 /Herwig/Decays/ScalarSS
 #
 decaymode a'_0-->omega,rho-; 0.827532 1 /Herwig/Decays/SVV
 decaymode a'_0-->pi-,eta; 0.07734 1 /Herwig/Decays/ScalarSS
 decaymode a'_0-->K-,K0; 0.068059 1 /Herwig/Decays/ScalarSS
 decaymode a'_0-->pi-,eta'; 0.027069 1 /Herwig/Decays/ScalarSS
 #
 decaymode f'_0->sigma,sigma; 0.438181 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->rho+,rho-; 0.150788 1 /Herwig/Decays/SVV
 decaymode f'_0->rho0,rho0; 0.075394 1 /Herwig/Decays/SVV
 decaymode f'_0->pi+,pi-; 0.052606 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->pi'-,pi+; 0.050424 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->pi'0,pi0; 0.050424 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->pi'+,pi-; 0.050424 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->K0,Kbar0; 0.026364 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->K+,K-; 0.026364 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->pi0,pi0; 0.026303 1 /Herwig/Decays/ScalarSS
 decaymode f'_0->a_1-,pi+; 0.017091 1 /Herwig/Decays/PPVDecayer
 decaymode f'_0->a_10,pi0; 0.017091 1 /Herwig/Decays/PPVDecayer
 decaymode f'_0->a_1+,pi-; 0.017091 1 /Herwig/Decays/PPVDecayer
 decaymode f'_0->eta,eta; 0.001455 1 /Herwig/Decays/ScalarSS
 #
 decaymode f_0(1710)->K0,Kbar0; 0.316456 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1710)->K+,K-; 0.316456 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1710)->eta,eta; 0.303797 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1710)->pi+,pi-; 0.042194 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1710)->pi0,pi0; 0.021097 1 /Herwig/Decays/ScalarSS
 #
 decaymode chi_c0->g,g; 0.788568 1 /Herwig/Decays/QDecayME0
 decaymode chi_c0->f'_0,f_0(1710); 0.026 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->rho-,pi+,pi0; 0.016 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->rho0,pi+,pi-; 0.016 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->rho+,pi-,pi0; 0.016 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->Jpsi,gamma; 0.013 1 /Herwig/Decays/SVV
 decaymode chi_c0->K_1-,K+; 0.01 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K_1bar0,K0; 0.01 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K_10,Kbar0; 0.01 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K_1+,K-; 0.01 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->pi-,pi+,pi-,pi+; 0.00739 1 /Herwig/Decays/Mambo
 decaymode chi_c0->K0,Kbar0; 0.0054 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K+,K-; 0.0054 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->pi+,pi-; 0.0048 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K-,K*0,pi+; 0.0044 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->Kbar0,K*+,pi-; 0.0044 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K0,K*-,pi+; 0.0044 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K+,K*bar0,pi-; 0.0044 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->pi+,pi-,K+,K-; 0.00277 1 /Herwig/Decays/Mambo
 decaymode chi_c0->pi0,pi0; 0.0024 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K*_00,K*_0bar0; 0.002363 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K*_0+,K*_0-; 0.002363 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->omega,omega; 0.0023 1 /Herwig/Decays/SVV
 decaymode chi_c0->K*-,K+,pi0; 0.0022 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K*bar0,K0,pi0; 0.0022 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K*0,Kbar0,pi0; 0.0022 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->K*+,K-,pi0; 0.0022 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->f_0,f_0; 0.00215 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K+,K-,K+,K-; 0.0021 1 /Herwig/Decays/Mambo
 decaymode chi_c0->K*_0-,K*_2+; 0.001912 1 /Herwig/Decays/STSDecayer
 decaymode chi_c0->K*_0bar0,K*_20; 0.001912 1 /Herwig/Decays/STSDecayer
 decaymode chi_c0->K*_00,K*_2bar0; 0.001912 1 /Herwig/Decays/STSDecayer
 decaymode chi_c0->K*_0+,K*_2-; 0.001912 1 /Herwig/Decays/STSDecayer
 decaymode chi_c0->pi+,pi-,p+,pbar-; 0.0019 1 /Herwig/Decays/Mambo
 decaymode chi_c0->eta,eta; 0.0019 1 /Herwig/Decays/ScalarSS
 decaymode chi_c0->K*0,K*bar0; 0.0018 1 /Herwig/Decays/SVV
 decaymode chi_c0->K*+,K*-; 0.0018 1 /Herwig/Decays/SVV
 decaymode chi_c0->K+,K-,K_S0,K_S0; 0.0015 1 /Herwig/Decays/Mambo
 decaymode chi_c0->phi,phi; 0.0009 1 /Herwig/Decays/SVV
 decaymode chi_c0->Lambda0,Lambdabar0; 0.00044 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->gamma,gamma; 0.000276 1 /Herwig/Decays/SVV
 decaymode chi_c0->n0,nbar0; 0.000216 1 /Herwig/Decays/DecayME0
 decaymode chi_c0->p+,pbar-; 0.000216 1 /Herwig/Decays/DecayME0
 #
 decaymode chi_b0->g,g; 0.968 1 /Herwig/Decays/QDecayME0
 decaymode chi_b0->Upsilon,gamma; 0.032 1 /Herwig/Decays/DecayME0
 #
 decaymode K*_0+->K0,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode K*_0+->K+,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode K*_00->K+,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode K*_00->K0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode K*_0-->Kbar0,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode K*_0-->K-,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode K*_0bar0->K-,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode K*_0bar0->Kbar0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_0*+->D0,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode D_0*+->D+,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_0*0->D+,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode D_0*0->D0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_0*-->Dbar0,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode D_0*-->D-,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_0*bar0->D-,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode D_0*bar0->Dbar0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_s0+->D_s+,pi0; 1 1 /Herwig/Decays/ScalarSS
 #
 decaymode D_s0-->D_s-,pi0; 1 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_00->B+,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode B*_00->B0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_0+->B0,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode B*_0+->B+,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_0bar0->B-,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode B*_0bar0->Bbar0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_0-->Bbar0,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode B*_0-->B-,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_s00->B0,Kbar0; 0.5 1 /Herwig/Decays/ScalarSS
 decaymode B*_s00->B+,K-; 0.5 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_s0bar0->B-,K+; 0.5 1 /Herwig/Decays/ScalarSS
 decaymode B*_s0bar0->Bbar0,K0; 0.5 1 /Herwig/Decays/ScalarSS
 #
 decaymode B*_c0+->B_c+,gamma; 1 1 /Herwig/Decays/DecayME0
 #
 decaymode B*_c0-->B_c-,gamma; 1 1 /Herwig/Decays/DecayME0
 #
 # the 1^3P_1 multiplet
 #
 #
 decaymode a_1+->pi+,pi0,pi0; 0.499099 1 /Herwig/Decays/a1Simple
 decaymode a_1+->pi+,pi+,pi-; 0.4881 1 /Herwig/Decays/a1Simple
 decaymode a_1+->K+,K*bar0; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_1+->K*+,Kbar0; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_1+->pi+,gamma; 1.0E-6 1 /Herwig/Decays/PVectorVP
 #
 decaymode a_10->pi+,pi-,pi0; 0.987199 1 /Herwig/Decays/a1Simple
 decaymode a_10->K*0,Kbar0; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_10->K*+,K-; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_10->pi0,gamma; 1.0E-6 1 /Herwig/Decays/PVectorVP
 #
 decaymode a_1-->pi-,pi0,pi0; 0.499099 1 /Herwig/Decays/a1Simple
 decaymode a_1-->pi-,pi-,pi+; 0.4881 1 /Herwig/Decays/a1Simple
 decaymode a_1-->K*-,K0; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_1-->K-,K*0; 0.0064 1 /Herwig/Decays/PVectorVP
 decaymode a_1-->pi-,gamma; 1.0E-6 1 /Herwig/Decays/PVectorVP
 #
 decaymode f_1->a_0-,pi+; 0.14257 1 /Herwig/Decays/Vector2Meson
 decaymode f_1->a_00,pi0; 0.14257 1 /Herwig/Decays/Vector2Meson
 decaymode f_1->a_0+,pi-; 0.14257 1 /Herwig/Decays/Vector2Meson
 decaymode f_1->rho-,pi+,pi0; 0.110471 1 /Herwig/Decays/DecayME0
 decaymode f_1->rho0,pi+,pi-; 0.110471 1 /Herwig/Decays/DecayME0
 decaymode f_1->rho+,pi-,pi0; 0.110471 1 /Herwig/Decays/DecayME0
 decaymode f_1->eta,pi+,pi-; 0.107123 1 /Herwig/Decays/DecayME0
 decaymode f_1->rho0,gamma; 0.055235 1 /Herwig/Decays/DecayME0
 decaymode f_1->eta,pi0,pi0; 0.053562 1 /Herwig/Decays/DecayME0
 decaymode f_1->K-,K0,pi+; 0.008071 1 /Herwig/Decays/DecayME0
 decaymode f_1->K+,Kbar0,pi-; 0.008071 1 /Herwig/Decays/DecayME0
 decaymode f_1->K0,Kbar0,pi0; 0.004036 1 /Herwig/Decays/DecayME0
 decaymode f_1->K+,K-,pi0; 0.004036 1 /Herwig/Decays/DecayME0
 decaymode f_1->phi,gamma; 0.000743 1 /Herwig/Decays/DecayME0
 #
 decaymode f'_1->K*-,K+; 0.1886 1 /Herwig/Decays/PVectorVP
 decaymode f'_1->K*bar0,K0; 0.1886 1 /Herwig/Decays/PVectorVP
 decaymode f'_1->K*0,Kbar0; 0.1886 1 /Herwig/Decays/PVectorVP
 decaymode f'_1->K*+,K-; 0.1886 1 /Herwig/Decays/PVectorVP
 decaymode f'_1->K-,K0,pi+; 0.0709 1 /Herwig/Decays/DecayME0
 decaymode f'_1->K+,Kbar0,pi-; 0.0709 1 /Herwig/Decays/DecayME0
 decaymode f'_1->K0,Kbar0,pi0; 0.0355 1 /Herwig/Decays/DecayME0
 decaymode f'_1->K+,K-,pi0; 0.0355 1 /Herwig/Decays/DecayME0
 decaymode f'_1->a_0-,pi+; 0.01 1 /Herwig/Decays/Vector2Meson
 decaymode f'_1->a_00,pi0; 0.01 1 /Herwig/Decays/Vector2Meson
 decaymode f'_1->a_0+,pi-; 0.01 1 /Herwig/Decays/Vector2Meson
 decaymode f'_1->phi,gamma; 0.0028 1 /Herwig/Decays/DecayME0
 #
 decaymode chi_c1->g,g; 0.608569 1 /Herwig/Decays/QDecayME0
 decaymode chi_c1->Jpsi,gamma; 0.356 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->rho-,pi+,pi0; 0.0039 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->rho0,pi+,pi-; 0.0039 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->rho+,pi-,pi0; 0.0039 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->pi+,pi-,pi+,pi-; 0.0037 1 /Herwig/Decays/Mambo
 decaymode chi_c1->pi+,pi-,K+,K-; 0.00237 1 /Herwig/Decays/Mambo
 decaymode chi_c1->K-,K0,pi+; 0.0023 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K+,Kbar0,pi-; 0.0023 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*0,K*bar0; 0.0016 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*+,K*-; 0.0016 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K0,Kbar0,pi0; 0.00115 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K+,K-,pi0; 0.00115 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*-,K0,pi+; 0.00107 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*bar0,K+,pi-; 0.00107 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*0,K-,pi+; 0.00107 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*+,Kbar0,pi-; 0.00107 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*0,Kbar0,pi0; 0.000537 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*-,K+,pi0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*bar0,K0,pi0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->K*+,K-,pi0; 0.00053 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->pi+,pi-,p+,pbar-; 0.00039 1 /Herwig/Decays/Mambo
 decaymode chi_c1->K+,K-,K+,K-; 0.00039 1 /Herwig/Decays/Mambo
 decaymode chi_c1->Lambda0,Lambdabar0; 0.00024 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->n0,nbar0; 6.7E-5 1 /Herwig/Decays/DecayME0
 decaymode chi_c1->p+,pbar-; 6.7E-5 1 /Herwig/Decays/DecayME0
 #
 decaymode chi_b1->g,g; 0.65 1 /Herwig/Decays/QDecayME0
 decaymode chi_b1->Upsilon,gamma; 0.35 1 /Herwig/Decays/DecayME0
 #
 decaymode K'_1+->K*0,pi+; 0.634667 1 /Herwig/Decays/PVectorVP
 decaymode K'_1+->K*+,pi0; 0.317333 1 /Herwig/Decays/PVectorVP
 decaymode K'_1+->K0,rho+; 0.024 1 /Herwig/Decays/PVectorVP
 decaymode K'_1+->K+,rho0; 0.012 1 /Herwig/Decays/PVectorVP
 decaymode K'_1+->K+,omega; 0.012 1 /Herwig/Decays/PVectorVP
 #
 decaymode K'_10->K*+,pi-; 0.634005 1 /Herwig/Decays/PVectorVP
 decaymode K'_10->K*0,pi0; 0.317002 1 /Herwig/Decays/PVectorVP
 decaymode K'_10->K+,rho-; 0.023669 1 /Herwig/Decays/PVectorVP
 decaymode K'_10->K0,rho0; 0.011834 1 /Herwig/Decays/PVectorVP
 decaymode K'_10->K0,omega; 0.011834 1 /Herwig/Decays/PVectorVP
 decaymode K'_10->K0,gamma; 0.001656 1 /Herwig/Decays/PVectorVP
 #
 decaymode K'_1-->K*bar0,pi-; 0.634667 1 /Herwig/Decays/PVectorVP
 decaymode K'_1-->K*-,pi0; 0.317333 1 /Herwig/Decays/PVectorVP
 decaymode K'_1-->Kbar0,rho-; 0.024 1 /Herwig/Decays/PVectorVP
 decaymode K'_1-->K-,rho0; 0.012 1 /Herwig/Decays/PVectorVP
 decaymode K'_1-->K-,omega; 0.012 1 /Herwig/Decays/PVectorVP
 #
 decaymode K'_1bar0->K*-,pi+; 0.634005 1 /Herwig/Decays/PVectorVP
 decaymode K'_1bar0->K*bar0,pi0; 0.317002 1 /Herwig/Decays/PVectorVP
 decaymode K'_1bar0->K-,rho+; 0.023669 1 /Herwig/Decays/PVectorVP
 decaymode K'_1bar0->Kbar0,rho0; 0.011834 1 /Herwig/Decays/PVectorVP
 decaymode K'_1bar0->Kbar0,omega; 0.011834 1 /Herwig/Decays/PVectorVP
 decaymode K'_1bar0->Kbar0,gamma; 0.001656 1 /Herwig/Decays/PVectorVP
 #
 decaymode D'_10->D*+,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D'_10->D*0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D'_1+->D*0,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D'_1+->D*+,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D'_1bar0->D*-,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D'_1bar0->D*bar0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D'_1-->D*bar0,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode D'_1-->D*-,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_s1+->D*0,K+; 0.559471 1 /Herwig/Decays/PVectorVP
 decaymode D_s1+->D*+,K0; 0.440529 1 /Herwig/Decays/PVectorVP
 #
 decaymode D_s1-->D*bar0,K-; 0.559471 1 /Herwig/Decays/PVectorVP
 decaymode D_s1-->D*-,Kbar0; 0.440529 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_10->B*+,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B'_10->B*0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_1+->B*0,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B'_1+->B*+,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_1bar0->B*-,pi+; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B'_1bar0->B*bar0,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_1-->B*bar0,pi-; 0.666667 1 /Herwig/Decays/PVectorVP
 decaymode B'_1-->B*-,pi0; 0.333333 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_s10->B*0,Kbar0; 0.5 1 /Herwig/Decays/PVectorVP
 decaymode B'_s10->B*+,K-; 0.5 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_s1bar0->B*-,K+; 0.5 1 /Herwig/Decays/PVectorVP
 decaymode B'_s1bar0->B*bar0,K0; 0.5 1 /Herwig/Decays/PVectorVP
 #
 decaymode B'_c1+->B_c+,gamma; 0.88 1 /Herwig/Decays/PVectorVP
 decaymode B'_c1+->B_c*+,gamma; 0.12 1 /Herwig/Decays/DecayME0
 #
 decaymode B'_c1-->B_c-,gamma; 0.88 1 /Herwig/Decays/PVectorVP
 decaymode B'_c1-->B_c*-,gamma; 0.12 1 /Herwig/Decays/DecayME0
 #
 # the 1^3P_2 multiplet
 #
 #
 decaymode a_2+->pi+,rho0; 0.348989 1 /Herwig/Decays/TensorVP
 decaymode a_2+->rho+,pi0; 0.348989 1 /Herwig/Decays/TensorVP
 decaymode a_2+->pi+,eta; 0.143657 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2+->omega,rho+; 0.102417 1 /Herwig/Decays/TVV
 decaymode a_2+->K+,Kbar0; 0.048104 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2+->eta',pi+; 0.005199 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2+->pi+,gamma; 0.002645 1 /Herwig/Decays/TensorVP
 #
 decaymode a_20->rho-,pi+; 0.349434 1 /Herwig/Decays/TensorVP
 decaymode a_20->rho+,pi-; 0.349434 1 /Herwig/Decays/TensorVP
 decaymode a_20->pi0,eta; 0.144053 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_20->omega,rho0; 0.103473 1 /Herwig/Decays/TVV
 decaymode a_20->K0,Kbar0; 0.024184 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_20->K+,K-; 0.024184 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_20->eta',pi0; 0.005229 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_20->gamma,gamma; 9.0E-6 1 /Herwig/Decays/TVV
 #
 decaymode a_2-->rho-,pi0; 0.348989 1 /Herwig/Decays/TensorVP
 decaymode a_2-->pi-,rho0; 0.348989 1 /Herwig/Decays/TensorVP
 decaymode a_2-->pi-,eta; 0.143657 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2-->omega,rho-; 0.102417 1 /Herwig/Decays/TVV
 decaymode a_2-->K-,K0; 0.048104 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2-->eta',pi-; 0.005199 1 /Herwig/Decays/Tensor2PScalar
 decaymode a_2-->pi-,gamma; 0.002645 1 /Herwig/Decays/TensorVP
 #
 decaymode f_2->pi+,pi-; 0.570829 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->pi0,pi0; 0.285414 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->rho+,rho-; 0.046681 1 /Herwig/Decays/TVV
 decaymode f_2->K0,Kbar0; 0.023516 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->K+,K-; 0.023516 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->rho0,rho0; 0.02334 1 /Herwig/Decays/TVV
 decaymode f_2->sigma,sigma; 0.022649 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->eta,eta; 0.004041 1 /Herwig/Decays/Tensor2PScalar
 decaymode f_2->gamma,gamma; 1.4E-5 1 /Herwig/Decays/TVV
 #
 decaymode f'_2->K0,Kbar0; 0.444355 1 /Herwig/Decays/Tensor2PScalar
 decaymode f'_2->K+,K-; 0.444355 1 /Herwig/Decays/Tensor2PScalar
 decaymode f'_2->eta,eta; 0.103083 1 /Herwig/Decays/Tensor2PScalar
 decaymode f'_2->pi+,pi-; 0.00547 1 /Herwig/Decays/Tensor2PScalar
 decaymode f'_2->pi0,pi0; 0.002736 1 /Herwig/Decays/Tensor2PScalar
 decaymode f'_2->gamma,gamma; 1.0E-6 1 /Herwig/Decays/TVV
 #
 decaymode chi_c2->g,g; 0.700639 1 /Herwig/Decays/QDecayME0
 decaymode chi_c2->Jpsi,gamma; 0.202 1 /Herwig/Decays/TVV
 decaymode chi_c2->K0,Kbar0; 0.0093 1 /Herwig/Decays/Tensor2PScalar
 decaymode chi_c2->K+,K-; 0.0093 1 /Herwig/Decays/Tensor2PScalar
 decaymode chi_c2->K0,K-,pi0,pi+; 0.008 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K0,Kbar0,pi+,pi-; 0.008 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K+,K-,pi+,pi-; 0.008 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K+,Kbar0,pi0,pi-; 0.008 1 /Herwig/Decays/Mambo
 decaymode chi_c2->rho-,pi+,pi0; 0.007 1 /Herwig/Decays/DecayME0
 decaymode chi_c2->rho0,pi+,pi-; 0.007 1 /Herwig/Decays/DecayME0
 decaymode chi_c2->rho+,pi-,pi0; 0.007 1 /Herwig/Decays/DecayME0
 decaymode chi_c2->pi+,pi-,pi+,pi-; 0.0053 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K*0,K*bar0; 0.0038 1 /Herwig/Decays/TVV
 decaymode chi_c2->K*+,K*-; 0.0038 1 /Herwig/Decays/TVV
 decaymode chi_c2->omega,omega; 0.002 1 /Herwig/Decays/TVV
 decaymode chi_c2->K0,Kbar0,pi0,pi0; 0.002 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K+,K-,pi0,pi0; 0.002 1 /Herwig/Decays/Mambo
 decaymode chi_c2->phi,phi; 0.0019 1 /Herwig/Decays/TVV
 decaymode chi_c2->pi+,pi-; 0.001427 1 /Herwig/Decays/Tensor2PScalar
 decaymode chi_c2->pi+,pi-,p+,pbar-; 0.00121 1 /Herwig/Decays/Mambo
 decaymode chi_c2->K+,K-,K+,K-; 0.00095 1 /Herwig/Decays/Mambo
 decaymode chi_c2->pi0,pi0; 0.000713 1 /Herwig/Decays/Tensor2PScalar
 decaymode chi_c2->Lambda0,Lambdabar0; 0.00027 1 /Herwig/Decays/DecayME0
 decaymode chi_c2->gamma,gamma; 0.000259 1 /Herwig/Decays/TVV
 decaymode chi_c2->n0,nbar0; 6.6E-5 1 /Herwig/Decays/DecayME0
 decaymode chi_c2->p+,pbar-; 6.6E-5 1 /Herwig/Decays/DecayME0
 #
 decaymode chi_b2->g,g; 0.78 1 /Herwig/Decays/QDecayME0
 decaymode chi_b2->Upsilon,gamma; 0.22 1 /Herwig/Decays/TVV
 #
 decaymode K*_2+->K0,pi+; 0.3327 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2+->K+,pi0; 0.16635 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2+->K*0,pi+; 0.164683 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K*0,rho+; 0.089342 1 /Herwig/Decays/TVV
 decaymode K*_2+->K*+,pi0; 0.082342 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K0,rho+; 0.058006 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K*+,rho0; 0.044671 1 /Herwig/Decays/TVV
 decaymode K*_2+->K+,rho0; 0.029003 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K+,omega; 0.029003 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K+,gamma; 0.0024 1 /Herwig/Decays/TensorVP
 decaymode K*_2+->K+,eta; 0.0015 1 /Herwig/Decays/Tensor2PScalar
 #
 decaymode K*_20->K+,pi-; 0.3335 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_20->K0,pi0; 0.16675 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_20->K*+,pi-; 0.165079 1 /Herwig/Decays/TensorVP
 decaymode K*_20->K*+,rho-; 0.089557 1 /Herwig/Decays/TVV
 decaymode K*_20->K*0,pi0; 0.08254 1 /Herwig/Decays/TensorVP
 decaymode K*_20->K+,rho-; 0.058145 1 /Herwig/Decays/TensorVP
 decaymode K*_20->K*0,rho0; 0.044779 1 /Herwig/Decays/TVV
 decaymode K*_20->K0,rho0; 0.029073 1 /Herwig/Decays/TensorVP
 decaymode K*_20->K0,omega; 0.029073 1 /Herwig/Decays/TensorVP
 decaymode K*_20->K0,eta; 0.001504 1 /Herwig/Decays/Tensor2PScalar
 #
 decaymode K*_2-->Kbar0,pi-; 0.3327 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2-->K-,pi0; 0.16635 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2-->K*bar0,pi-; 0.164683 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->K*bar0,rho-; 0.089342 1 /Herwig/Decays/TVV
 decaymode K*_2-->K*-,pi0; 0.082342 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->Kbar0,rho-; 0.058006 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->K*-,rho0; 0.044671 1 /Herwig/Decays/TVV
 decaymode K*_2-->K-,rho0; 0.029003 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->K-,omega; 0.029003 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->K-,gamma; 0.0024 1 /Herwig/Decays/TensorVP
 decaymode K*_2-->K-,eta; 0.0015 1 /Herwig/Decays/Tensor2PScalar
 #
 decaymode K*_2bar0->K-,pi+; 0.3335 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2bar0->Kbar0,pi0; 0.16675 1 /Herwig/Decays/Tensor2PScalar
 decaymode K*_2bar0->K*-,pi+; 0.165079 1 /Herwig/Decays/TensorVP
 decaymode K*_2bar0->K*-,rho+; 0.089557 1 /Herwig/Decays/TVV
 decaymode K*_2bar0->K*bar0,pi0; 0.08254 1 /Herwig/Decays/TensorVP
 decaymode K*_2bar0->K-,rho+; 0.058145 1 /Herwig/Decays/TensorVP
 decaymode K*_2bar0->K*bar0,rho0; 0.044779 1 /Herwig/Decays/TVV
 decaymode K*_2bar0->Kbar0,rho0; 0.029073 1 /Herwig/Decays/TensorVP
 decaymode K*_2bar0->Kbar0,omega; 0.029073 1 /Herwig/Decays/TensorVP
 decaymode K*_2bar0->Kbar0,eta; 0.001504 1 /Herwig/Decays/Tensor2PScalar
 #
-decaymode D*_20->D+,pi-; 0.448087 1 /Herwig/Decays/Tensor2PScalar
-decaymode D*_20->D0,pi0; 0.224044 1 /Herwig/Decays/Tensor2PScalar
-decaymode D*_20->D*+,pi-; 0.218579 1 /Herwig/Decays/TensorVP
-decaymode D*_20->D*0,pi0; 0.10929 1 /Herwig/Decays/TensorVP
+decaymode D*_20->D+,pi-; 0.448087 1 /Herwig/Decays/HQETMeson
+decaymode D*_20->D0,pi0; 0.224044 1 /Herwig/Decays/HQETMeson
+decaymode D*_20->D*+,pi-; 0.218579 1 /Herwig/Decays/HQETMeson
+decaymode D*_20->D*0,pi0; 0.10929 1 /Herwig/Decays/HQETMeson
 #
-decaymode D*_2+->D0,pi+; 0.448087 1 /Herwig/Decays/Tensor2PScalar
-decaymode D*_2+->D+,pi0; 0.224044 1 /Herwig/Decays/Tensor2PScalar
-decaymode D*_2+->D*0,pi+; 0.218579 1 /Herwig/Decays/TensorVP
-decaymode D*_2+->D*+,pi0; 0.10929 1 /Herwig/Decays/TensorVP
+decaymode D*_2+->D0,pi+; 0.448087 1 /Herwig/Decays/HQETMeson
+decaymode D*_2+->D+,pi0; 0.224044 1 /Herwig/Decays/HQETMeson
+decaymode D*_2+->D*0,pi+; 0.218579 1 /Herwig/Decays/HQETMeson
+decaymode D*_2+->D*+,pi0; 0.10929 1 /Herwig/Decays/HQETMeson
 #
 decaymode D*_2bar0->D-,pi+; 0.448087 1 /Herwig/Decays/Tensor2PScalar
 decaymode D*_2bar0->Dbar0,pi0; 0.224044 1 /Herwig/Decays/Tensor2PScalar
 decaymode D*_2bar0->D*-,pi+; 0.218579 1 /Herwig/Decays/TensorVP
 decaymode D*_2bar0->D*bar0,pi0; 0.10929 1 /Herwig/Decays/TensorVP
 #
 decaymode D*_2-->Dbar0,pi-; 0.448087 1 /Herwig/Decays/Tensor2PScalar
 decaymode D*_2-->D-,pi0; 0.224044 1 /Herwig/Decays/Tensor2PScalar
 decaymode D*_2-->D*bar0,pi-; 0.218579 1 /Herwig/Decays/TensorVP
 decaymode D*_2-->D*-,pi0; 0.10929 1 /Herwig/Decays/TensorVP
 #
-decaymode D_s2+->D0,K+; 0.453903 1 /Herwig/Decays/Tensor2PScalar
-decaymode D_s2+->D+,K0; 0.410175 1 /Herwig/Decays/Tensor2PScalar
-decaymode D_s2+->D*0,K+; 0.077088 1 /Herwig/Decays/TensorVP
-decaymode D_s2+->D*+,K0; 0.058834 1 /Herwig/Decays/TensorVP
+decaymode D_s2+->D0,K+; 0.453903 1 /Herwig/Decays/HQETMeson
+decaymode D_s2+->D+,K0; 0.410175 1 /Herwig/Decays/HQETMeson
+decaymode D_s2+->D*0,K+; 0.077088 1 /Herwig/Decays/HQETMeson
+decaymode D_s2+->D*+,K0; 0.058834 1 /Herwig/Decays/HQETMeson
 #
-decaymode D_s2-->Dbar0,K-; 0.453903 1 /Herwig/Decays/Tensor2PScalar
-decaymode D_s2-->D-,Kbar0; 0.410175 1 /Herwig/Decays/Tensor2PScalar
-decaymode D_s2-->D*bar0,K-; 0.077088 1 /Herwig/Decays/TensorVP
-decaymode D_s2-->D*-,Kbar0; 0.058834 1 /Herwig/Decays/TensorVP
+decaymode D_s2-->Dbar0,K-; 0.453903 1 /Herwig/Decays/HQETMeson
+decaymode D_s2-->D-,Kbar0; 0.410175 1 /Herwig/Decays/HQETMeson
+decaymode D_s2-->D*bar0,K-; 0.077088 1 /Herwig/Decays/HQETMeson
+decaymode D_s2-->D*-,Kbar0; 0.058834 1 /Herwig/Decays/HQETMeson
 #
 decaymode B_20->B+,pi-; 0.345404 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_20->B*+,pi-; 0.321262 1 /Herwig/Decays/TensorVP
 decaymode B_20->B0,pi0; 0.172703 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_20->B*0,pi0; 0.160631 1 /Herwig/Decays/TensorVP
 #
 decaymode B_2+->B0,pi+; 0.345404 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2+->B*0,pi+; 0.321262 1 /Herwig/Decays/TensorVP
 decaymode B_2+->B+,pi0; 0.172703 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2+->B*+,pi0; 0.160631 1 /Herwig/Decays/TensorVP
 #
 decaymode B_2bar0->B-,pi+; 0.345404 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2bar0->B*-,pi+; 0.321262 1 /Herwig/Decays/TensorVP
 decaymode B_2bar0->Bbar0,pi0; 0.172703 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2bar0->B*bar0,pi0; 0.160631 1 /Herwig/Decays/TensorVP
 #
 decaymode B_2-->Bbar0,pi-; 0.345404 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2-->B*bar0,pi-; 0.321262 1 /Herwig/Decays/TensorVP
 decaymode B_2-->B-,pi0; 0.172703 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_2-->B*-,pi0; 0.160631 1 /Herwig/Decays/TensorVP
 #
 decaymode B_s20->B0,Kbar0; 0.47 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_s20->B+,K-; 0.47 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_s20->B*0,Kbar0; 0.03 1 /Herwig/Decays/TensorVP
 decaymode B_s20->B*+,K-; 0.03 1 /Herwig/Decays/TensorVP
 #
 decaymode B_s2bar0->B-,K+; 0.47 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_s2bar0->Bbar0,K0; 0.47 1 /Herwig/Decays/Tensor2PScalar
 decaymode B_s2bar0->B*-,K+; 0.03 1 /Herwig/Decays/TensorVP
 decaymode B_s2bar0->B*bar0,K0; 0.03 1 /Herwig/Decays/TensorVP
 #
 decaymode B_c2+->B_c+,gamma; 1 1 /Herwig/Decays/TensorVP
 #
 decaymode B_c2-->B_c-,gamma; 1 1 /Herwig/Decays/TensorVP
 #
 # the 1^1D_2 multiplet
 #
 #
 decaymode pi_2+->pi+,f_2; 0.546654 1 /Herwig/Decays/DecayME0
 decaymode pi_2+->rho0,pi+; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_2+->rho+,pi0; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_2+->pi+,sigma; 0.084624 1 /Herwig/Decays/Tensor2PScalar
 decaymode pi_2+->rho+,omega; 0.02629 1 /Herwig/Decays/DecayME0
 decaymode pi_2+->K*bar0,K+; 0.020448 1 /Herwig/Decays/DecayME0
 decaymode pi_2+->K*+,Kbar0; 0.020448 1 /Herwig/Decays/DecayME0
 #
 decaymode pi_20->pi0,f_2; 0.546654 1 /Herwig/Decays/DecayME0
 decaymode pi_20->rho-,pi+; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_20->rho+,pi-; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_20->pi0,sigma; 0.084624 1 /Herwig/Decays/Tensor2PScalar
 decaymode pi_20->rho0,omega; 0.02629 1 /Herwig/Decays/DecayME0
 decaymode pi_20->K*-,K+; 0.010224 1 /Herwig/Decays/DecayME0
 decaymode pi_20->K*bar0,K0; 0.010224 1 /Herwig/Decays/DecayME0
 decaymode pi_20->K*0,Kbar0; 0.010224 1 /Herwig/Decays/DecayME0
 decaymode pi_20->K*+,K-; 0.010224 1 /Herwig/Decays/DecayME0
 #
 decaymode pi_2-->pi-,f_2; 0.546654 1 /Herwig/Decays/DecayME0
 decaymode pi_2-->rho-,pi0; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_2-->rho0,pi-; 0.150768 1 /Herwig/Decays/DecayME0
 decaymode pi_2-->pi-,sigma; 0.084624 1 /Herwig/Decays/Tensor2PScalar
 decaymode pi_2-->rho-,omega; 0.02629 1 /Herwig/Decays/DecayME0
 decaymode pi_2-->K*-,K0; 0.020448 1 /Herwig/Decays/DecayME0
 decaymode pi_2-->K*0,K-; 0.020448 1 /Herwig/Decays/DecayME0
 #
 decaymode eta_2->a_20,pi0; 0.290634 1 /Herwig/Decays/DecayME0
 decaymode eta_2->a_2-,pi+; 0.290633 1 /Herwig/Decays/DecayME0
 decaymode eta_2->a_2+,pi-; 0.290633 1 /Herwig/Decays/DecayME0
 decaymode eta_2->a_0-,pi+; 0.022356 1 /Herwig/Decays/Tensor2PScalar
 decaymode eta_2->a_00,pi0; 0.022356 1 /Herwig/Decays/Tensor2PScalar
 decaymode eta_2->a_0+,pi-; 0.022356 1 /Herwig/Decays/Tensor2PScalar
 decaymode eta_2->K*-,K+; 0.015258 1 /Herwig/Decays/DecayME0
 decaymode eta_2->K*bar0,K0; 0.015258 1 /Herwig/Decays/DecayME0
 decaymode eta_2->K*0,Kbar0; 0.015258 1 /Herwig/Decays/DecayME0
 decaymode eta_2->K*+,K-; 0.015258 1 /Herwig/Decays/DecayME0
 #
 decaymode eta'_2->a_2-,pi+; 0.278395 1 /Herwig/Decays/DecayME0
 decaymode eta'_2->a_20,pi0; 0.278395 1 /Herwig/Decays/DecayME0
 decaymode eta'_2->a_2+,pi-; 0.278395 1 /Herwig/Decays/DecayME0
 decaymode eta'_2->f_2,eta; 0.139195 1 /Herwig/Decays/DecayME0
 decaymode eta'_2->a_0-,pi+; 0.00854 1 /Herwig/Decays/Tensor2PScalar
 decaymode eta'_2->a_00,pi0; 0.00854 1 /Herwig/Decays/Tensor2PScalar
 decaymode eta'_2->a_0+,pi-; 0.00854 1 /Herwig/Decays/Tensor2PScalar
 #
 decaymode eta_b2->h_b,gamma; 0.9 1 /Herwig/Decays/DecayME0
 decaymode eta_b2->g,g; 0.1 1 /Herwig/Decays/QDecayME0
 #
 decaymode K_2(1770)+->K+,f_2; 0.74 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)+->K*0,pi+; 0.153333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)+->K*+,pi0; 0.076667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)+->K*_20,pi+; 0.02 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)+->K*_2+,pi0; 0.01 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1770)0->K0,f_2; 0.74 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)0->K*+,pi-; 0.153333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)0->K*0,pi0; 0.076667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)0->K*_2+,pi-; 0.02 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)0->K*_20,pi0; 0.01 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1770)-->K-,f_2; 0.74 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)-->K*bar0,pi-; 0.153333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)-->K*-,pi0; 0.076667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)-->K*_2bar0,pi-; 0.02 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)-->K*_2-,pi0; 0.01 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1770)bar0->Kbar0,f_2; 0.74 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)bar0->K*-,pi+; 0.153333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)bar0->K*bar0,pi0; 0.076667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)bar0->K*_2-,pi+; 0.02 1 /Herwig/Decays/DecayME0
 decaymode K_2(1770)bar0->K*_2bar0,pi0; 0.01 1 /Herwig/Decays/DecayME0
 #
 # the 1^3D_1 multiplet
 #
 #
 decaymode rho''+->rho+,sigma; 0.209974 1 /Herwig/Decays/DecayME0
 decaymode rho''+->h_1,pi+; 0.127484 1 /Herwig/Decays/VAPDecayer
 decaymode rho''+->pi+,pi0; 0.119985 1 /Herwig/Decays/Vector2Meson
 decaymode rho''+->pi'0,pi+; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''+->pi'+,pi0; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''+->rho+,eta; 0.092239 1 /Herwig/Decays/VectorVP
 decaymode rho''+->rho+,rho0; 0.067492 1 /Herwig/Decays/VectorVV
 decaymode rho''+->a_10,pi+; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''+->a_1+,pi0; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''+->K*bar0,K+; 0.017998 1 /Herwig/Decays/VectorVP
 decaymode rho''+->K*+,Kbar0; 0.017998 1 /Herwig/Decays/VectorVP
 decaymode rho''+->K+,Kbar0; 0.001872 1 /Herwig/Decays/Vector2Meson
 #
 decaymode rho''0->rho0,sigma; 0.209974 1 /Herwig/Decays/DecayME0
 decaymode rho''0->h_1,pi0; 0.127484 1 /Herwig/Decays/VAPDecayer
 decaymode rho''0->pi+,pi-; 0.119985 1 /Herwig/Decays/Vector2Meson
 decaymode rho''0->pi'-,pi+; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''0->pi'+,pi-; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''0->rho0,eta; 0.092239 1 /Herwig/Decays/VectorVP
 decaymode rho''0->rho+,rho-; 0.067492 1 /Herwig/Decays/VectorVV
 decaymode rho''0->a_1-,pi+; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''0->a_1+,pi-; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''0->K*-,K+; 0.008999 1 /Herwig/Decays/VectorVP
 decaymode rho''0->K*bar0,K0; 0.008999 1 /Herwig/Decays/VectorVP
 decaymode rho''0->K*0,Kbar0; 0.008999 1 /Herwig/Decays/VectorVP
 decaymode rho''0->K*+,K-; 0.008999 1 /Herwig/Decays/VectorVP
 decaymode rho''0->K-,K+; 0.000936 1 /Herwig/Decays/Vector2Meson
 decaymode rho''0->Kbar0,K0; 0.000936 1 /Herwig/Decays/Vector2Meson
 #
 decaymode rho''-->rho-,sigma; 0.209974 1 /Herwig/Decays/DecayME0
 decaymode rho''-->h_1,pi-; 0.127484 1 /Herwig/Decays/VAPDecayer
 decaymode rho''-->pi-,pi0; 0.119985 1 /Herwig/Decays/Vector2Meson
 decaymode rho''-->pi'-,pi0; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''-->pi'0,pi-; 0.112486 1 /Herwig/Decays/Vector2Meson
 decaymode rho''-->rho-,eta; 0.092239 1 /Herwig/Decays/VectorVP
 decaymode rho''-->rho-,rho0; 0.067492 1 /Herwig/Decays/VectorVV
 decaymode rho''-->a_1-,pi0; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''-->a_10,pi-; 0.059993 1 /Herwig/Decays/VAPDecayer
 decaymode rho''-->K*-,K0; 0.017998 1 /Herwig/Decays/VectorVP
 decaymode rho''-->K*0,K-; 0.017998 1 /Herwig/Decays/VectorVP
 decaymode rho''-->K-,K0; 0.001872 1 /Herwig/Decays/Vector2Meson
 #
 decaymode omega''->omega,pi+,pi-; 0.333333 1 /Herwig/Decays/DecayME0
 decaymode omega''->rho-,pi+; 0.166667 1 /Herwig/Decays/VectorVP
 decaymode omega''->rho+,pi-; 0.166667 1 /Herwig/Decays/VectorVP
 decaymode omega''->omega,pi0,pi0; 0.166667 1 /Herwig/Decays/DecayME0
 decaymode omega''->rho0,pi0; 0.166666 1 /Herwig/Decays/VectorVP
 #
 decaymode psi(3770)->D0,Dbar0; 0.595969 1 /Herwig/Decays/Vector2Meson
 decaymode psi(3770)->D+,D-; 0.400064 1 /Herwig/Decays/Vector2Meson
 decaymode psi(3770)->Jpsi,pi+,pi-; 0.00193 1 /Herwig/Decays/OniumPiPi
 decaymode psi(3770)->Jpsi,eta; 0.0009 1 /Herwig/Decays/VectorVP
 decaymode psi(3770)->Jpsi,pi0,pi0; 0.0008 1 /Herwig/Decays/OniumPiPi
 decaymode psi(3770)->phi,eta; 0.00031 1 /Herwig/Decays/VectorVP
 decaymode psi(3770)->e-,e+; 1.1E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(3770)->mu-,mu+; 1.1E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(3770)->tau-,tau+; 5.0E-6 1 /Herwig/Decays/Vector2Leptons
 #
 decaymode Upsilon_1(1D)->chi_b0,gamma; 0.602 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(1D)->chi_b1,gamma; 0.318 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(1D)->g,g,g; 0.062 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon_1(1D)->chi_b2,gamma; 0.016 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(1D)->Upsilon,pi+,pi-; 0.001333 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(1D)->Upsilon,pi0,pi0; 0.000667 1 /Herwig/Decays/DecayME0
 #
 decaymode K''*+->K0,pi+; 0.258 1 /Herwig/Decays/Vector2Meson
 decaymode K''*+->K0,rho+; 0.209333 1 /Herwig/Decays/VectorVP
 decaymode K''*+->K*0,pi+; 0.199333 1 /Herwig/Decays/VectorVP
 decaymode K''*+->K+,pi0; 0.129 1 /Herwig/Decays/Vector2Meson
 decaymode K''*+->K+,rho0; 0.104667 1 /Herwig/Decays/VectorVP
 decaymode K''*+->K*+,pi0; 0.099667 1 /Herwig/Decays/VectorVP
 #
 decaymode K''*0->K+,pi-; 0.258 1 /Herwig/Decays/Vector2Meson
 decaymode K''*0->K+,rho-; 0.209333 1 /Herwig/Decays/VectorVP
 decaymode K''*0->K*+,pi-; 0.199333 1 /Herwig/Decays/VectorVP
 decaymode K''*0->K0,pi0; 0.129 1 /Herwig/Decays/Vector2Meson
 decaymode K''*0->K0,rho0; 0.104667 1 /Herwig/Decays/VectorVP
 decaymode K''*0->K*0,pi0; 0.099667 1 /Herwig/Decays/VectorVP
 #
 decaymode K''*-->Kbar0,pi-; 0.258 1 /Herwig/Decays/Vector2Meson
 decaymode K''*-->Kbar0,rho-; 0.209333 1 /Herwig/Decays/VectorVP
 decaymode K''*-->K*bar0,pi-; 0.199333 1 /Herwig/Decays/VectorVP
 decaymode K''*-->K-,pi0; 0.129 1 /Herwig/Decays/Vector2Meson
 decaymode K''*-->K-,rho0; 0.104667 1 /Herwig/Decays/VectorVP
 decaymode K''*-->K*-,pi0; 0.099667 1 /Herwig/Decays/VectorVP
 #
 decaymode K''*bar0->K-,pi+; 0.258 1 /Herwig/Decays/Vector2Meson
 decaymode K''*bar0->K-,rho+; 0.209333 1 /Herwig/Decays/VectorVP
 decaymode K''*bar0->K*-,pi+; 0.199333 1 /Herwig/Decays/VectorVP
 decaymode K''*bar0->Kbar0,pi0; 0.129 1 /Herwig/Decays/Vector2Meson
 decaymode K''*bar0->Kbar0,rho0; 0.104667 1 /Herwig/Decays/VectorVP
 decaymode K''*bar0->K*bar0,pi0; 0.099667 1 /Herwig/Decays/VectorVP
 #
 # the 1^3D_2 multiplet
 #
 #
 decaymode Upsilon_2(1D)->chi_b0,gamma; 0.786 1 /Herwig/Decays/TensorVP
 decaymode Upsilon_2(1D)->chi_b2,gamma; 0.203 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_2(1D)->g,g,g; 0.009 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon_2(1D)->Upsilon,pi+,pi-; 0.001333 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_2(1D)->Upsilon,pi0,pi0; 0.000667 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1820)+->K*_20,pi+; 0.513333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)+->K*_2+,pi0; 0.256667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)+->K+,f_2; 0.18 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)+->K*0,pi+; 0.033333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)+->K*+,pi0; 0.016667 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1820)0->K*_2+,pi-; 0.513333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)0->K*_20,pi0; 0.256667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)0->K0,f_2; 0.18 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)0->K*+,pi-; 0.033333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)0->K*0,pi0; 0.016667 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1820)-->K*_2bar0,pi-; 0.513333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)-->K*_2-,pi0; 0.256667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)-->K-,f_2; 0.18 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)-->K*bar0,pi-; 0.033333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)-->K*-,pi0; 0.016667 1 /Herwig/Decays/DecayME0
 #
 decaymode K_2(1820)bar0->K*_2-,pi+; 0.513333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)bar0->K*_2bar0,pi0; 0.256667 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)bar0->Kbar0,f_2; 0.18 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)bar0->K*-,pi+; 0.033333 1 /Herwig/Decays/DecayME0
 decaymode K_2(1820)bar0->K*bar0,pi0; 0.016667 1 /Herwig/Decays/DecayME0
 #
 # the 1^3D_3 multiplet
 #
 #
 decaymode rho_3+->pi+,pi0; 0.23581 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->a_20,pi+; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->a_2+,pi0; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->pi+,omega; 0.16 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->rho+,rho0; 0.081547 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->K*bar0,K+; 0.018985 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->K*+,Kbar0; 0.018985 1 /Herwig/Decays/DecayME0
 decaymode rho_3+->K+,Kbar0; 0.015787 1 /Herwig/Decays/DecayME0
 #
 decaymode rho_30->pi+,pi-; 0.235811 1 /Herwig/Decays/DecayME0
 decaymode rho_30->a_2-,pi+; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_30->a_2+,pi-; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_30->pi0,omega; 0.16 1 /Herwig/Decays/DecayME0
 decaymode rho_30->rho+,rho-; 0.081547 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K*-,K+; 0.009492 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K*bar0,K0; 0.009492 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K*0,Kbar0; 0.009492 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K*+,K-; 0.009492 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K0,Kbar0; 0.007894 1 /Herwig/Decays/DecayME0
 decaymode rho_30->K+,K-; 0.007894 1 /Herwig/Decays/DecayME0
 #
 decaymode rho_3-->pi-,pi0; 0.23581 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->a_2-,pi0; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->a_20,pi-; 0.234443 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->pi-,omega; 0.16 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->rho-,rho0; 0.081547 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->K*-,K0; 0.018985 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->K*0,K-; 0.018985 1 /Herwig/Decays/DecayME0
 decaymode rho_3-->K-,K0; 0.015787 1 /Herwig/Decays/DecayME0
 #
 decaymode omega_3->rho0,pi0; 0.194933 1 /Herwig/Decays/DecayME0
 decaymode omega_3->rho-,pi+; 0.194932 1 /Herwig/Decays/DecayME0
 decaymode omega_3->rho+,pi-; 0.194932 1 /Herwig/Decays/DecayME0
 decaymode omega_3->b_1-,pi+; 0.138401 1 /Herwig/Decays/DecayME0
 decaymode omega_3->b_10,pi0; 0.138401 1 /Herwig/Decays/DecayME0
 decaymode omega_3->b_1+,pi-; 0.138401 1 /Herwig/Decays/DecayME0
 #
 decaymode phi_3->K0,Kbar0; 0.322582 1 /Herwig/Decays/DecayME0
 decaymode phi_3->K+,K-; 0.322582 1 /Herwig/Decays/DecayME0
 decaymode phi_3->K*-,K+; 0.088709 1 /Herwig/Decays/DecayME0
 decaymode phi_3->K*bar0,K0; 0.088709 1 /Herwig/Decays/DecayME0
 decaymode phi_3->K*0,Kbar0; 0.088709 1 /Herwig/Decays/DecayME0
 decaymode phi_3->K*+,K-; 0.088709 1 /Herwig/Decays/DecayME0
 #
 decaymode Upsilon_3(1D)->chi_b2,gamma; 0.954 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_3(1D)->g,g,g; 0.043 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon_3(1D)->Upsilon,pi+,pi-; 0.002 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_3(1D)->Upsilon,pi0,pi0; 0.001 1 /Herwig/Decays/DecayME0
 #
 decaymode K_3*+->K+,eta; 0.300601 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K0,rho+; 0.207081 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K*0,pi+; 0.133601 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K0,pi+; 0.125585 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K+,rho0; 0.10354 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K*+,pi0; 0.0668 1 /Herwig/Decays/DecayME0
 decaymode K_3*+->K+,pi0; 0.062792 1 /Herwig/Decays/DecayME0
 #
 decaymode K_3*0->K0,eta; 0.300601 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K+,rho-; 0.207081 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K*+,pi-; 0.133601 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K+,pi-; 0.125585 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K0,rho0; 0.10354 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K*0,pi0; 0.0668 1 /Herwig/Decays/DecayME0
 decaymode K_3*0->K0,pi0; 0.062792 1 /Herwig/Decays/DecayME0
 #
 decaymode K_3*-->K-,eta; 0.300601 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->Kbar0,rho-; 0.207081 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->K*bar0,pi-; 0.133601 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->Kbar0,pi-; 0.125585 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->K-,rho0; 0.10354 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->K*-,pi0; 0.0668 1 /Herwig/Decays/DecayME0
 decaymode K_3*-->K-,pi0; 0.062792 1 /Herwig/Decays/DecayME0
 #
 decaymode K_3*bar0->Kbar0,eta; 0.300601 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->K-,rho+; 0.207081 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->K*-,pi+; 0.133601 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->K-,pi+; 0.125585 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->Kbar0,rho0; 0.10354 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->K*bar0,pi0; 0.0668 1 /Herwig/Decays/DecayME0
 decaymode K_3*bar0->Kbar0,pi0; 0.062792 1 /Herwig/Decays/DecayME0
 #
 # the 1^3F_4 multiplet
 #
 #
 # the 2^1S_0 multiplet
 #
 #
 decaymode pi'+->rho0,pi+; 0.5 1 /Herwig/Decays/PPVDecayer
 decaymode pi'+->rho+,pi0; 0.5 1 /Herwig/Decays/PPVDecayer
 #
 decaymode pi'0->rho-,pi+; 0.5 1 /Herwig/Decays/PPVDecayer
 decaymode pi'0->rho+,pi-; 0.5 1 /Herwig/Decays/PPVDecayer
 #
 decaymode pi'-->rho-,pi0; 0.5 1 /Herwig/Decays/PPVDecayer
 decaymode pi'-->rho0,pi-; 0.5 1 /Herwig/Decays/PPVDecayer
 #
 decaymode eta(1295)->sigma,eta; 0.349999 1 /Herwig/Decays/ScalarSS
 decaymode eta(1295)->a_0-,pi+; 0.216667 1 /Herwig/Decays/ScalarSS
 decaymode eta(1295)->a_00,pi0; 0.216667 1 /Herwig/Decays/ScalarSS
 decaymode eta(1295)->a_0+,pi-; 0.216667 1 /Herwig/Decays/ScalarSS
 #
 decaymode eta(1475)->K*-,K+; 0.25 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1475)->K*bar0,K0; 0.25 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1475)->K*0,Kbar0; 0.25 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1475)->K*+,K-; 0.25 1 /Herwig/Decays/PPVDecayer
 #
 decaymode eta_c(2S)->g,g; 0.9752 1 /Herwig/Decays/QDecayME0
 decaymode eta_c(2S)->eta_c,pi+,pi-; 0.0125 1 /Herwig/Decays/DecayME0
 decaymode eta_c(2S)->eta_c,pi0,pi0; 0.0063 1 /Herwig/Decays/DecayME0
 decaymode eta_c(2S)->h_c,gamma; 0.006 1 /Herwig/Decays/DecayME0
 #
 decaymode eta_b(2S)->g,g; 0.9943 1 /Herwig/Decays/QDecayME0
 decaymode eta_b(2S)->h_b,gamma; 0.0057 1 /Herwig/Decays/DecayME0
 #
 decaymode K'+->K*_00,pi+; 0.3 1 /Herwig/Decays/ScalarSS
 decaymode K'+->K*0,pi+; 0.279487 1 /Herwig/Decays/PPVDecayer
 decaymode K'+->K*_0+,pi0; 0.15 1 /Herwig/Decays/ScalarSS
 decaymode K'+->K*+,pi0; 0.139744 1 /Herwig/Decays/PPVDecayer
 decaymode K'+->K0,rho+; 0.087179 1 /Herwig/Decays/PPVDecayer
 decaymode K'+->K+,rho0; 0.04359 1 /Herwig/Decays/PPVDecayer
 #
 decaymode K'0->K*_0+,pi-; 0.3 1 /Herwig/Decays/ScalarSS
 decaymode K'0->K*+,pi-; 0.279487 1 /Herwig/Decays/PPVDecayer
 decaymode K'0->K*_00,pi0; 0.15 1 /Herwig/Decays/ScalarSS
 decaymode K'0->K*0,pi0; 0.139744 1 /Herwig/Decays/PPVDecayer
 decaymode K'0->K+,rho-; 0.087179 1 /Herwig/Decays/PPVDecayer
 decaymode K'0->K0,rho0; 0.04359 1 /Herwig/Decays/PPVDecayer
 #
 decaymode K'-->K*_0bar0,pi-; 0.3 1 /Herwig/Decays/ScalarSS
 decaymode K'-->K*bar0,pi-; 0.279487 1 /Herwig/Decays/PPVDecayer
 decaymode K'-->K*_0-,pi0; 0.15 1 /Herwig/Decays/ScalarSS
 decaymode K'-->K*-,pi0; 0.139744 1 /Herwig/Decays/PPVDecayer
 decaymode K'-->Kbar0,rho-; 0.087179 1 /Herwig/Decays/PPVDecayer
 decaymode K'-->K-,rho0; 0.04359 1 /Herwig/Decays/PPVDecayer
 #
 decaymode K'bar0->K*_0-,pi+; 0.3 1 /Herwig/Decays/ScalarSS
 decaymode K'bar0->K*-,pi+; 0.279487 1 /Herwig/Decays/PPVDecayer
 decaymode K'bar0->K*_0bar0,pi0; 0.15 1 /Herwig/Decays/ScalarSS
 decaymode K'bar0->K*bar0,pi0; 0.139744 1 /Herwig/Decays/PPVDecayer
 decaymode K'bar0->K-,rho+; 0.087179 1 /Herwig/Decays/PPVDecayer
 decaymode K'bar0->Kbar0,rho0; 0.04359 1 /Herwig/Decays/PPVDecayer
 #
 # the 2^3S_1 multiplet
 #
 #
 decaymode rho'+->pi+,pi0; 0.270072 1 /Herwig/Decays/Vector2Meson
 decaymode rho'+->pi'0,pi+; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'+->pi'+,pi0; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'+->rho+,sigma; 0.124088 1 /Herwig/Decays/VVSDecayer
 decaymode rho'+->a_10,pi+; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'+->a_1+,pi0; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'+->rho+,rho0; 0.080292 1 /Herwig/Decays/VectorVV
 decaymode rho'+->h_1,pi+; 0.058394 1 /Herwig/Decays/VAPDecayer
 #
 decaymode rho'0->pi+,pi-; 0.270072 1 /Herwig/Decays/Vector2Meson
 decaymode rho'0->pi'-,pi+; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'0->pi'+,pi-; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'0->rho0,sigma; 0.124088 1 /Herwig/Decays/VVSDecayer
 decaymode rho'0->a_1-,pi+; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'0->a_1+,pi-; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'0->rho+,rho-; 0.080292 1 /Herwig/Decays/VectorVV
 decaymode rho'0->h_1,pi0; 0.058394 1 /Herwig/Decays/VAPDecayer
 #
 decaymode rho'-->pi-,pi0; 0.270072 1 /Herwig/Decays/Vector2Meson
 decaymode rho'-->pi'-,pi0; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'-->pi'0,pi-; 0.135037 1 /Herwig/Decays/Vector2Meson
 decaymode rho'-->rho-,sigma; 0.124088 1 /Herwig/Decays/VVSDecayer
 decaymode rho'-->a_1-,pi0; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'-->a_10,pi-; 0.09854 1 /Herwig/Decays/VAPDecayer
 decaymode rho'-->rho-,rho0; 0.080292 1 /Herwig/Decays/VectorVV
 decaymode rho'-->h_1,pi-; 0.058394 1 /Herwig/Decays/VAPDecayer
 #
 decaymode omega'->rho-,pi+; 0.233 1 /Herwig/Decays/VectorVP
 decaymode omega'->rho0,pi0; 0.233 1 /Herwig/Decays/VectorVP
 decaymode omega'->rho+,pi-; 0.233 1 /Herwig/Decays/VectorVP
 decaymode omega'->omega,sigma; 0.112876 1 /Herwig/Decays/DecayME0
 decaymode omega'->b_1-,pi+; 0.062708 1 /Herwig/Decays/VAPDecayer
 decaymode omega'->b_10,pi0; 0.062708 1 /Herwig/Decays/VAPDecayer
 decaymode omega'->b_1+,pi-; 0.062708 1 /Herwig/Decays/VAPDecayer
 #
 decaymode phi'->K*-,K+; 0.233645 1 /Herwig/Decays/VectorVP
 decaymode phi'->K*bar0,K0; 0.233645 1 /Herwig/Decays/VectorVP
 decaymode phi'->K*0,Kbar0; 0.233645 1 /Herwig/Decays/VectorVP
 decaymode phi'->K*+,K-; 0.233645 1 /Herwig/Decays/VectorVP
 decaymode phi'->K_L0,K_S0; 0.03271 1 /Herwig/Decays/Vector2Meson
 decaymode phi'->K+,K-; 0.03271 1 /Herwig/Decays/Vector2Meson
 #
 decaymode psi(2S)->Jpsi,pi+,pi-; 0.318 1 /Herwig/Decays/OniumPiPi
 decaymode psi(2S)->g,g,g; 0.1830929 1 /Herwig/Decays/QDecayME130
 decaymode psi(2S)->Jpsi,pi0,pi0; 0.1646 1 /Herwig/Decays/OniumPiPi
 decaymode psi(2S)->chi_c0,gamma; 0.092 1 /Herwig/Decays/VVSDecayer
 decaymode psi(2S)->chi_c1,gamma; 0.087 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->chi_c2,gamma; 0.081 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Jpsi,eta; 0.0309 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->e-,e+; 0.00735 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(2S)->mu-,mu+; 0.00734 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(2S)->tau-,tau+; 0.00285 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(2S)->eta_c,gamma; 0.0026 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->Jpsi,pi0; 0.00126 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->K_1-,K+; 0.001 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->K_1bar0,K0; 0.001 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->K_10,Kbar0; 0.001 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->K_1+,K-; 0.001 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->h_c,pi0; 0.001 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->eta,pi+,pi-,pi0; 0.00095 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->p+,pbar-,pi+,pi-,pi0; 0.00073 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->eta',pi+,pi-,pi0; 0.00045 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->pi+,pi-,p+,pbar-; 0.0004 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->K*-,K0,rho+; 0.00035 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*bar0,K+,rho-; 0.00035 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*0,K-,rho+; 0.00035 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*+,Kbar0,rho-; 0.00035 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*-,K0,pi+; 0.000335 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*bar0,K+,pi-; 0.000335 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*0,K-,pi+; 0.000335 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*+,Kbar0,pi-; 0.000335 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->b_1-,pi+; 0.0003 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->b_10,pi0; 0.0003 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->b_1+,pi-; 0.0003 1 /Herwig/Decays/VAPDecayer
 decaymode psi(2S)->Lambda0,Lambdabar0,pi+,pi-; 0.00028 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->n0,pbar-,pi+; 0.000266 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,nbar0,pi-; 0.000266 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->n0,nbar0; 0.000265 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->p+,pbar-; 0.000265 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Lambda0,Lambdabar0; 0.00025 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Sigma-,Sigmabar+; 0.00023 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(2S)->Sigma0,Sigmabar0; 0.00023 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Sigma+,Sigmabar-; 0.00023 1 /Herwig/Decays/Vector2Leptons
 decaymode psi(2S)->rho-,pi+,pi0; 0.00022 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->rho0,pi+,pi-; 0.00022 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->rho+,pi-,pi0; 0.00022 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->f_2,gamma; 0.00021 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,f_2; 0.0002 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Xi-,Xibar+; 0.0002 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Xi0,Xibar0; 0.0002 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->K*-,K*_2+; 0.00019 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*bar0,K*_20; 0.00019 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,f_0(1710); 0.00019 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*0,K*_2bar0; 0.00019 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*+,K*_2-; 0.00019 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,K0,Kbar0; 0.000185 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,K+,K-; 0.000185 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Lambdabar0,p+,K-,pi+,pi-; 0.00018 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->Lambda0,pbar-,K+,pi+,pi-; 0.00018 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->K*-,K+,rho0; 0.000175 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*bar0,K0,rho0; 0.000175 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*0,Kbar0,rho0; 0.000175 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*+,K-,rho0; 0.000175 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*-,K+,pi0; 0.000168 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*bar0,K0,pi0; 0.000168 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*0,Kbar0,pi0; 0.000168 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K*+,K-,pi0; 0.000168 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->eta',gamma; 0.00015 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->pi+,pi-,pi0; 0.000136 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->n0,nbar0,pi0; 0.000133 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,pbar-,pi0; 0.000133 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Delta-,Deltabar+; 0.000128 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Delta0,Deltabar0; 0.000128 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Delta+,Deltabar-; 0.000128 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Delta++,Deltabar--; 0.000128 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K+,K-,K+,K-,pi0; 0.00011 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->Sigma*-,Sigma*bar+; 0.00011 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Sigma*0,Sigma*bar0; 0.00011 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Sigma*+,Sigma*bar-; 0.00011 1 /Herwig/Decays/Vector2Baryons
 decaymode psi(2S)->Kbar0,K*0; 0.000109 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->K0,K*bar0; 0.000109 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->Lambdabar0,n0,Kbar0; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Lambdabar0,p+,K-; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->phi,f_0; 0.0001 1 /Herwig/Decays/VVSDecayer
 decaymode psi(2S)->n0,pbar-,rho+; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,nbar0,rho-; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Lambda0,pbar-,K+; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->Lambda0,nbar0,K0; 0.0001 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->f_0(1710),gamma; 0.0001 1 /Herwig/Decays/VVSDecayer
 decaymode psi(2S)->rho-,a_2+; 8.7E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->rho0,a_20; 8.7E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->rho+,a_2-; 8.7E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->pi+,pi-; 8.0E-5 1 /Herwig/Decays/Vector2Meson
 decaymode psi(2S)->phi,K0,Kbar0; 7.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->phi,K+,K-; 7.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,n0,nbar0; 6.9E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,p+,pbar-; 6.9E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K+,K-,K+,K-; 6.0E-5 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->n0,nbar0,eta; 6.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,pbar-,eta; 6.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->K_L0,K_S0; 5.3E-5 1 /Herwig/Decays/Vector2Meson
 decaymode psi(2S)->K+,K-; 5.3E-5 1 /Herwig/Decays/Vector2Meson
 decaymode psi(2S)->phi,pi+,pi-; 5.3E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->n0,nbar0,rho0; 5.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,pbar-,rho0; 5.0E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->phi,f'_2; 4.4E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->omega,eta'; 3.2E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->phi,eta'; 3.1E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->phi,eta; 2.8E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->phi,pi0,pi0; 2.7E-5 1 /Herwig/Decays/DecayME0
 decaymode psi(2S)->p+,pbar-,K+,K-; 2.7E-5 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->rho0,eta; 2.2E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->omega,pi0; 2.1E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->pi+,pi-,pi+,pi-; 2.0E-5 1 /Herwig/Decays/Mambo
 decaymode psi(2S)->rho0,eta'; 1.9E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->K-,K*+; 1.7E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->K+,K*-; 1.7E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->rho-,pi+; 1.1E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->rho0,pi0; 1.1E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->rho+,pi-; 1.1E-5 1 /Herwig/Decays/VectorVP
 decaymode psi(2S)->Xi*-,Xi*bar+; 11.51e-5 1 /Herwig/Decays/Vector2Baryons
 #
 decaymode Upsilon(2S)->g,g,g; 0.4077 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon(2S)->Upsilon,pi+,pi-; 0.188 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(2S)->Upsilon,pi0,pi0; 0.09 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(2S)->chi_b2,gamma; 0.0715 1 /Herwig/Decays/DecayME0
 decaymode Upsilon(2S)->chi_b1,gamma; 0.069 1 /Herwig/Decays/DecayME0
 decaymode Upsilon(2S)->chi_b0,gamma; 0.038 1 /Herwig/Decays/VVSDecayer
 decaymode Upsilon(2S)->u,ubar; 0.0271 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(2S)->c,cbar; 0.0271 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(2S)->e-,e+; 0.0192 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(2S)->mu-,mu+; 0.0192 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(2S)->tau-,tau+; 0.019 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(2S)->g,g,gamma; 0.0106 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon(2S)->d,dbar; 0.0068 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(2S)->s,sbar; 0.0068 1 /Herwig/Decays/QDecayME0
 #
 decaymode K'*+->K*0,pi+; 0.622667 1 /Herwig/Decays/VectorVP
 decaymode K'*+->K*+,pi0; 0.311333 1 /Herwig/Decays/VectorVP
 decaymode K'*+->K0,pi+; 0.044 1 /Herwig/Decays/Vector2Meson
 decaymode K'*+->K+,pi0; 0.022 1 /Herwig/Decays/Vector2Meson
 #
 decaymode K'*0->K*+,pi-; 0.622667 1 /Herwig/Decays/VectorVP
 decaymode K'*0->K*0,pi0; 0.311333 1 /Herwig/Decays/VectorVP
 decaymode K'*0->K+,pi-; 0.044 1 /Herwig/Decays/Vector2Meson
 decaymode K'*0->K0,pi0; 0.022 1 /Herwig/Decays/Vector2Meson
 #
 decaymode K'*-->K*bar0,pi-; 0.622667 1 /Herwig/Decays/VectorVP
 decaymode K'*-->K*-,pi0; 0.311333 1 /Herwig/Decays/VectorVP
 decaymode K'*-->Kbar0,pi-; 0.044 1 /Herwig/Decays/Vector2Meson
 decaymode K'*-->K-,pi0; 0.022 1 /Herwig/Decays/Vector2Meson
 #
 decaymode K'*bar0->K*-,pi+; 0.622667 1 /Herwig/Decays/VectorVP
 decaymode K'*bar0->K*bar0,pi0; 0.311333 1 /Herwig/Decays/VectorVP
 decaymode K'*bar0->K-,pi+; 0.044 1 /Herwig/Decays/Vector2Meson
 decaymode K'*bar0->Kbar0,pi0; 0.022 1 /Herwig/Decays/Vector2Meson
 #
 # the 2^1P_1 multiplet
 #
 #
 decaymode h_b(2P)->g,g,g; 0.628 1 /Herwig/Decays/QDecayME130
 decaymode h_b(2P)->eta_b(2S),gamma; 0.193 1 /Herwig/Decays/DecayME0
 decaymode h_b(2P)->eta_b,gamma; 0.125 1 /Herwig/Decays/DecayME0
 decaymode h_b(2P)->eta_b2,gamma; 0.034 1 /Herwig/Decays/DecayME0
 decaymode h_b(2P)->g,g,gamma; 0.02 1 /Herwig/Decays/QDecayME130
 #
 # the 2^3P_0 multiplet
 #
 #
 decaymode chi_b0(2P)->g,g; 0.9435 1 /Herwig/Decays/QDecayME0
 decaymode chi_b0(2P)->Upsilon(2S),gamma; 0.046 1 /Herwig/Decays/DecayME0
 decaymode chi_b0(2P)->Upsilon,gamma; 0.009 1 /Herwig/Decays/DecayME0
 decaymode chi_b0(2P)->Upsilon_1(1D),gamma; 0.0015 1 /Herwig/Decays/DecayME0
 #
 # the 2^3P_1 multiplet
 #
 #
 decaymode chi_b1(2P)->g,g; 0.6468 1 /Herwig/Decays/QDecayME0
 decaymode chi_b1(2P)->Upsilon(2S),gamma; 0.21 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->Upsilon,gamma; 0.085 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->Upsilon_2(1D),gamma; 0.0236 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->omega,Upsilon; 0.0163 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->Upsilon_1(1D),gamma; 0.0097 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->chi_b1,pi+,pi-; 0.0057 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(2P)->chi_b1,pi0,pi0; 0.0029 1 /Herwig/Decays/DecayME0
 #
 # the 2^3P_2 multiplet
 #
 #
 decaymode chi_c2(2P)->D0,Dbar0; 0.5291 1 /Herwig/Decays/Tensor2PScalar
 decaymode chi_c2(2P)->D+,D-; 0.4709 1 /Herwig/Decays/Tensor2PScalar
 #
 decaymode chi_b2(2P)->g,g; 0.73267 1 /Herwig/Decays/QDecayME0
 decaymode chi_b2(2P)->Upsilon(2S),gamma; 0.162 1 /Herwig/Decays/TVV
 decaymode chi_b2(2P)->Upsilon,gamma; 0.071 1 /Herwig/Decays/TVV
 decaymode chi_b2(2P)->Upsilon_3(1D),gamma; 0.0142 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(2P)->omega,Upsilon; 0.011 1 /Herwig/Decays/TVV
 decaymode chi_b2(2P)->chi_b2,pi+,pi-; 0.004 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(2P)->Upsilon_2(1D),gamma; 0.0029 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(2P)->chi_b2,pi0,pi0; 0.002 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(2P)->Upsilon_1(1D),gamma; 0.00023 1 /Herwig/Decays/DecayME0
 #
 # the 2^1D_2 multiplet
 #
 #
 # the 2^3D_1 multiplet
 #
 #
 decaymode Upsilon_1(2D)->chi_b0(2P),gamma; 0.5381 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(2D)->chi_b1(2P),gamma; 0.2728 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(2D)->chi_b0,gamma; 0.1046 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(2D)->chi_b1,gamma; 0.0663 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(2D)->chi_b2(2P),gamma; 0.0144 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_1(2D)->chi_b2,gamma; 0.0038 1 /Herwig/Decays/DecayME0
 #
 # the 2^3D_2 multiplet
 #
 #
 decaymode Upsilon_2(2D)->chi_b1(2P),gamma; 0.6382 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_2(2D)->chi_b2(2P),gamma; 0.1721 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_2(2D)->chi_b1,gamma; 0.1461 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_2(2D)->chi_b2,gamma; 0.0436 1 /Herwig/Decays/DecayME0
 #
 # the 2^3D_3 multiplet
 #
 #
 decaymode Upsilon_3(2D)->chi_b2(2P),gamma; 0.8053 1 /Herwig/Decays/DecayME0
 decaymode Upsilon_3(2D)->chi_b2,gamma; 0.1947 1 /Herwig/Decays/DecayME0
 #
 # the 3^1S_0 multiplet
 #
 #
 decaymode eta_b(3S)->g,g; 0.9983 1 /Herwig/Decays/QDecayME0
 decaymode eta_b(3S)->h_b(2P),gamma; 0.0012 1 /Herwig/Decays/DecayME0
 decaymode eta_b(3S)->h_b,gamma; 0.0005 1 /Herwig/Decays/DecayME0
 #
 # the 3^3S_1 multiplet
 #
 #
 decaymode Upsilon(3S)->g,g,g; 0.3649 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon(3S)->chi_b2(2P),gamma; 0.131 1 /Herwig/Decays/DecayME0
 decaymode Upsilon(3S)->chi_b1(2P),gamma; 0.126 1 /Herwig/Decays/DecayME0
 decaymode Upsilon(3S)->chi_b0(2P),gamma; 0.059 1 /Herwig/Decays/VVSDecayer
 decaymode Upsilon(3S)->Upsilon(2S),gamma,gamma; 0.05 1 /Herwig/Decays/DecayME0
 decaymode Upsilon(3S)->Upsilon,pi+,pi-; 0.0448 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(3S)->u,ubar; 0.031 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(3S)->c,cbar; 0.031 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(3S)->Upsilon(2S),pi+,pi-; 0.028 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(3S)->e-,e+; 0.0218 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(3S)->mu-,mu+; 0.0218 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(3S)->tau-,tau+; 0.0216 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(3S)->Upsilon,pi0,pi0; 0.0206 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(3S)->Upsilon(2S),pi0,pi0; 0.02 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(3S)->g,g,gamma; 0.0095 1 /Herwig/Decays/QDecayME130
 decaymode Upsilon(3S)->d,dbar; 0.008 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(3S)->s,sbar; 0.008 1 /Herwig/Decays/QDecayME0
 decaymode Upsilon(3S)->chi_b0,gamma; 0.003 1 /Herwig/Decays/DecayME0
 #
 # the 3^1P_1 multiplet
 #
 #
 # the 3^3P_0 multiplet
 #
 #
 decaymode chi_b0(3P)->g,g; 0.9773 1 /Herwig/Decays/QDecayME0
 decaymode chi_b0(3P)->Upsilon(3S),gamma; 0.0105 1 /Herwig/Decays/DecayME0
 decaymode chi_b0(3P)->Upsilon,gamma; 0.005 1 /Herwig/Decays/DecayME0
 decaymode chi_b0(3P)->Upsilon(2S),gamma; 0.0046 1 /Herwig/Decays/DecayME0
 decaymode chi_b0(3P)->Upsilon_1(2D),gamma; 0.0026 1 /Herwig/Decays/DecayME0
 #
 # the 3^3P_1 multiplet
 #
 #
 decaymode chi_b1(3P)->g,g; 0.6925 1 /Herwig/Decays/QDecayME0
 decaymode chi_b1(3P)->Upsilon(3S),gamma; 0.1531 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(3P)->Upsilon,gamma; 0.0559 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(3P)->Upsilon(2S),gamma; 0.054 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(3P)->Upsilon_2(2D),gamma; 0.0315 1 /Herwig/Decays/DecayME0
 decaymode chi_b1(3P)->Upsilon_1(2D),gamma; 0.013 1 /Herwig/Decays/DecayME0
 #
 # the 3^3P_2 multiplet
 #
 #
 decaymode chi_b2(3P)->g,g; 0.8404 1 /Herwig/Decays/QDecayME0
 decaymode chi_b2(3P)->Upsilon(3S),gamma; 0.0827 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(3P)->Upsilon,gamma; 0.026 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(3P)->Upsilon(2S),gamma; 0.0258 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(3P)->Upsilon_3(2D),gamma; 0.0206 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(3P)->Upsilon_2(2D),gamma; 0.0042 1 /Herwig/Decays/DecayME0
 decaymode chi_b2(3P)->Upsilon_1(2D),gamma; 0.0003 1 /Herwig/Decays/DecayME0
 #
 # the 4^3S_1 multiplet
 #
 #
 decaymode Upsilon(4S)->B+,B-; 0.505812 1 /Herwig/Decays/Vector2Meson
 decaymode Upsilon(4S)->B0,Bbar0; 0.493811 1 /Herwig/Decays/Vector2Meson
 decaymode Upsilon(4S)->Upsilon(2S),pi+,pi-; 0.000129 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(4S)->Upsilon,pi+,pi-; 9.0E-5 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(4S)->Upsilon(2S),pi0,pi0; 6.5E-5 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(4S)->Upsilon,pi0,pi0; 4.5E-5 1 /Herwig/Decays/OniumPiPi
 decaymode Upsilon(4S)->e-,e+; 1.6E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(4S)->mu-,mu+; 1.6E-5 1 /Herwig/Decays/Vector2Leptons
 decaymode Upsilon(4S)->tau-,tau+; 1.6E-5 1 /Herwig/Decays/Vector2Leptons
 #
 # the 5^3S_1 multiplet
 #
 #
 # The meson which are not part of a multiplet
 #
 #
 decaymode f_0(1500)->pi+,pi-; 0.232667 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->sigma,sigma; 0.127426 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->pi0,pi0; 0.116333 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->pi'-,pi+; 0.081683 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->pi'0,pi0; 0.081683 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->pi'+,pi-; 0.081683 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->eta,eta; 0.051 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->K0,Kbar0; 0.043 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->K+,K-; 0.043 1 /Herwig/Decays/ScalarSS
 decaymode f_0(1500)->rho+,rho-; 0.042475 1 /Herwig/Decays/SVV
 decaymode f_0(1500)->rho0,rho0; 0.021238 1 /Herwig/Decays/SVV
 decaymode f_0(1500)->a_1-,pi+; 0.019604 1 /Herwig/Decays/PPVDecayer
 decaymode f_0(1500)->a_10,pi0; 0.019604 1 /Herwig/Decays/PPVDecayer
 decaymode f_0(1500)->a_1+,pi-; 0.019604 1 /Herwig/Decays/PPVDecayer
 decaymode f_0(1500)->eta,eta'; 0.019 1 /Herwig/Decays/ScalarSS
 #
 decaymode sigma->pi+,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode sigma->pi0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode f_0->pi+,pi-; 0.493994 1 /Herwig/Decays/ScalarSS
 decaymode f_0->pi0,pi0; 0.247 1 /Herwig/Decays/ScalarSS
 decaymode f_0->K+,K-; 0.139 1 /Herwig/Decays/ScalarSS
 decaymode f_0->K0,Kbar0; 0.12 1 /Herwig/Decays/ScalarSS
 decaymode f_0->gamma,gamma; 6.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode a_00->eta,pi0; 0.906997 1 /Herwig/Decays/ScalarSS
 decaymode a_00->K+,K-; 0.049 1 /Herwig/Decays/ScalarSS
 decaymode a_00->K0,Kbar0; 0.044 1 /Herwig/Decays/ScalarSS
 decaymode a_00->gamma,gamma; 3.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode a_0+->eta,pi+; 0.907 1 /Herwig/Decays/ScalarSS
 decaymode a_0+->K+,Kbar0; 0.093 1 /Herwig/Decays/ScalarSS
 #
 decaymode a_0-->eta,pi-; 0.907 1 /Herwig/Decays/ScalarSS
 decaymode a_0-->K-,K0; 0.093 1 /Herwig/Decays/ScalarSS
 #
 decaymode eta(1405)->eta,sigma; 0.479845 1 /Herwig/Decays/ScalarSS
 decaymode eta(1405)->a_0-,pi+; 0.159949 1 /Herwig/Decays/ScalarSS
 decaymode eta(1405)->a_00,pi0; 0.159949 1 /Herwig/Decays/ScalarSS
 decaymode eta(1405)->a_0+,pi-; 0.159949 1 /Herwig/Decays/ScalarSS
 decaymode eta(1405)->K*-,K+; 0.010077 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1405)->K*bar0,K0; 0.010077 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1405)->K*0,Kbar0; 0.010077 1 /Herwig/Decays/PPVDecayer
 decaymode eta(1405)->K*+,K-; 0.010077 1 /Herwig/Decays/PPVDecayer
 #
 decaymode K_S0->pi+,pi-; 0.69224 1 /Herwig/Decays/DecayME0
 decaymode K_S0->pi0,pi0; 0.307006 1 /Herwig/Decays/DecayME0
 decaymode K_S0->pi-,e+,nu_e; 0.000352 1 /Herwig/Decays/DecayME101
 decaymode K_S0->pi+,e-,nu_ebar; 0.000352 1 /Herwig/Decays/DecayME101
 decaymode K_S0->pi+,pi-,e-,e+; 4.7E-5 1 /Herwig/Decays/Mambo
 decaymode K_S0->gamma,gamma; 3.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode K_L0->pi-,e+,nu_e; 0.202618 1 /Herwig/Decays/DecayME101
 decaymode K_L0->pi+,e-,nu_ebar; 0.202618 1 /Herwig/Decays/DecayME101
 decaymode K_L0->pi0,pi0,pi0; 0.195569 1 /Herwig/Decays/DecayME0
 decaymode K_L0->pi-,mu+,nu_mu; 0.135079 1 /Herwig/Decays/DecayME101
 decaymode K_L0->pi+,mu-,nu_mubar; 0.135079 1 /Herwig/Decays/DecayME101
 decaymode K_L0->pi+,pi-,pi0; 0.125581 1 /Herwig/Decays/DecayME0
 decaymode K_L0->pi+,pi-; 0.001976 1 /Herwig/Decays/DecayME0
 decaymode K_L0->pi0,pi0; 0.000869 1 /Herwig/Decays/DecayME0
 decaymode K_L0->gamma,gamma; 0.000548 1 /Herwig/Decays/DecayME0
 decaymode K_L0->pi0,pi-,e+,nu_e; 2.6E-5 1 /Herwig/Decays/Mambo
 decaymode K_L0->pi0,pi+,e-,nu_ebar; 2.6E-5 1 /Herwig/Decays/Mambo
 decaymode K_L0->e-,e+,gamma; 1.0E-5 1 /Herwig/Decays/DecayME0
 decaymode K_L0->pi0,gamma,gamma; 1.0E-6 1 /Herwig/Decays/DecayME0
 #
 decaymode kappa0->K+,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode kappa0->K0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode kappabar0->K-,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode kappabar0->Kbar0,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode kappa-->Kbar0,pi-; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode kappa-->K-,pi0; 0.333333 1 /Herwig/Decays/ScalarSS
 #
 decaymode kappa+->K0,pi+; 0.666667 1 /Herwig/Decays/ScalarSS
 decaymode kappa+->K+,pi0; 0.333333 1 /Herwig/Decays/ScalarSS