AliDecayer realisation for the EvtGen code and EvtGen itself.

(Fiorella Fionda)

diff --git a/TEvtGen/AliDecayerEvtGen.cxx b/TEvtGen/AliDecayerEvtGen.cxx
new file mode 100644 (file)
index 0000000..efd158f
--- /dev/null
+++ b/TEvtGen/AliDecayerEvtGen.cxx
@@ -0,0 +1,286 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+///////////////////////////////////////////////////////////////////////////
+//  Implementation of AliDecayer using EvtGen package.                   //
+//                                                                       //
+//  Giuseppe E. Bruno            &    Fiorella Fionda                    //
+//  (Giuseppe.Bruno@ba.infn.it)       (Fiorella.Fionda@ba.infn.it)       //
+///////////////////////////////////////////////////////////////////////////
+
+#include <TSystem.h>
+#include <TDatabasePDG.h>
+#include <TParticle.h>
+#include <TClonesArray.h>
+#include <TLorentzVector.h>
+
+#include "EvtGenBase/EvtStdHep.hh"
+#include "EvtGenBase/EvtRandomEngine.hh"
+#include "EvtGen/EvtGen.hh"
+#include "EvtGenBase/EvtParticle.hh"
+#include "EvtGenBase/EvtPDL.hh"
+#include "EvtGenBase/EvtParticleFactory.hh"
+#include "AliDecayerEvtGen.h"
+#include "AliLog.h"
+
+ClassImp(AliDecayerEvtGen)
+//____________________________________________________________
+AliDecayerEvtGen::AliDecayerEvtGen():
+  fRandomEngine(0x0),
+  fGenerator(0x0),
+  fEvtstdhep(0x0),
+  fDecayTablePath(0x0),
+  fParticleTablePath(0x0),
+  fDecay(kAll)
+  {
+  // Default constructor
+  fEvtstdhep = new EvtStdHep();
+  fDecayTablePath = gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DECAY.DEC"); //default decay table
+  fParticleTablePath = gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/evt.pdl"); //particle table
+  }
+//_________________________________________________________________
+AliDecayerEvtGen::AliDecayerEvtGen(const AliDecayerEvtGen &decayer):
+  AliDecayer(decayer),
+  fRandomEngine(decayer.fRandomEngine),
+  fGenerator(decayer.fGenerator),
+  fEvtstdhep(decayer.fEvtstdhep),
+  fDecayTablePath(decayer.fDecayTablePath),
+  fParticleTablePath(decayer.fParticleTablePath),
+  fDecay(decayer.fDecay)
+  {
+  // Copy Constructor
+  decayer.Copy(*this);
+  }
+//____________________________________________________________
+AliDecayerEvtGen::~AliDecayerEvtGen()
+  {
+  // Destructor
+  if(fRandomEngine) {delete fRandomEngine;}
+  fRandomEngine = 0;
+  if(fGenerator) {delete fGenerator;}
+  fGenerator = 0;
+  if(fEvtstdhep) {delete fEvtstdhep;}
+  fEvtstdhep = 0;
+  if(fDecayTablePath) {delete fDecayTablePath;}
+  fDecayTablePath = 0;
+  if(fParticleTablePath) {delete fParticleTablePath;}
+  fParticleTablePath = 0;
+  }
+
+//___________________________________________________________
+void AliDecayerEvtGen::Init()
+  {
+  //Standard AliDecayerEvtGen initializer:
+  //initialize EvtGen with default decay table (DECAY.DEC), particle table (evt.pdl)
+  //and fRandomEngine for generation of random numbers
+  //
+  if(fGenerator){
+  AliWarning(" AliDecayerEvtGen already initialized!!!!\n");
+  return;
+  }
+  fRandomEngine=new EvtNUMRandomEngine();
+  fGenerator=new EvtGen(fDecayTablePath,fParticleTablePath,fRandomEngine);
+  }
+//____________________________________________________________
+void AliDecayerEvtGen::Decay(Int_t ipart, TLorentzVector *p)
+  {
+  //
+  //Decay a particle
+  //input: pdg code and momentum of the particle to be decayed  
+  //all informations about decay products are stored in fEvtstdhep 
+  //
+  EvtId IPART=EvtPDL::evtIdFromStdHep(ipart);
+  EvtVector4R p_init=EvtVector4R::EvtVector4R(p->E(),p->Px(),p->Py(),p->Pz());
+  EvtParticle *froot_part=EvtParticleFactory::particleFactory(IPART,p_init);
+  fGenerator->generateDecay(froot_part);
+  fEvtstdhep->init();
+  froot_part->makeStdHep(*fEvtstdhep);
+  //froot_part->printTree(); //to print the decay chain 
+  froot_part->deleteTree();
+  }
+
+//____________________________________________________________
+Int_t AliDecayerEvtGen::ImportParticles(TClonesArray *particles)
+  {
+  //
+  //Input: pointer to a TClonesArray - Output(Int_t): number of decay products    
+  //Put all the informations about the decay products in the 
+  //TClonesArray particles
+  //
+  if (particles == 0) return 0;
+  TClonesArray &clonesParticles = *particles;
+  clonesParticles.Clear();
+
+  int j;
+  int istat;
+  int partnum;
+  double px,py,pz,e;
+  double x,y,z,t;
+  EvtVector4R p4,x4;
+
+  Int_t npart=fEvtstdhep->getNPart();
+  for(int i=0;i<fEvtstdhep->getNPart();i++){
+  j=i+1;
+  int jmotherfirst=fEvtstdhep->getFirstMother(i)+1;
+  int jmotherlast=fEvtstdhep->getLastMother(i)+1;
+  int jdaugfirst=fEvtstdhep->getFirstDaughter(i)+1;
+  int jdauglast=fEvtstdhep->getLastDaughter(i)+1;
+
+  partnum=fEvtstdhep->getStdHepID(i);
+  
+  //verify if all particles of decay chain are in the TDatabasePDG
+  TParticlePDG *partPDG = TDatabasePDG::Instance()->GetParticle(partnum);
+  if(!partPDG)
+    {
+    AliWarning("Particle code non known in TDatabasePDG - set pdg = 89");
+    partnum=89; //internal use for unspecified resonance data
+    }
+
+  istat=fEvtstdhep->getIStat(i);
+
+  if(istat!=1 && istat!=2) Info("ImportParticles","Attention: unknown status code!");
+  if(istat == 2) istat = 11; //status decayed
+
+  p4=fEvtstdhep->getP4(i);
+  x4=fEvtstdhep->getX4(i);
+  px=p4.get(1);
+  py=p4.get(2);
+  pz=p4.get(3);
+  e=p4.get(0);
+
+  x=x4.get(1);//[mm]
+  y=x4.get(2);//[mm]
+  z=x4.get(3);//[mm]
+  t=x4.get(0);//[mm]
+
+  AliDebug(1,Form("partnum = %d istat = %d primaMadre = %d ultimaMadre = %d primaF = %d ultimaF=%d x=%f y=%f z=%f t=%f e=%f px=%f \n",partnum,istat,jmotherfirst,jmotherlast,jdaugfirst,jdauglast,x,y,z,t,e,px));
+
+  new(clonesParticles[i]) TParticle(partnum,istat,jmotherfirst,-1,jdaugfirst,jdauglast,px,py,pz,e,x,y,z,t);
+
+  //set polarization!!!
+  }
+
+  return npart; 
+
+  }
+
+void  AliDecayerEvtGen::Copy(TObject &) const
+  {
+  //
+  // Copy *this onto AliDecayerEvtGen -- not implemented
+  //
+   Fatal("Copy","Not implemented!\n");
+  }
+
+void AliDecayerEvtGen::ForceDecay()
+  {
+  //
+  // Intupt: none - Output: none
+  // Set the decay mode to decay particles: for each case is read a 
+  // different decay table. case kAll read the default decay table only   
+  //
+  Decay_t decay = fDecay;
+  switch(decay)
+    {
+     case kAll: 
+     break;
+     case kBJpsiDiElectron:
+     SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOJPSITOELE.DEC"));
+     break;
+     case kBJpsi:
+     SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOJPSI.DEC"));
+     break;
+     case kBJpsiDiMuon:
+     SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOJPSITOMU.DEC"));
+     break;
+     case kBSemiElectronic:
+     SetDecayTablePath(gSystem->ExpandPathName("$ALICE_ROOT/TEvtGen/EvtGen/DecayTable/BTOELE.DEC"));
+     break;  
+     case kHardMuons:
+     case kChiToJpsiGammaToMuonMuon:
+     case kChiToJpsiGammaToElectronElectron:
+     case kBSemiMuonic:
+     case kSemiMuonic:
+     case kDiMuon:
+     case kSemiElectronic:
+     case kDiElectron:
+     case kBPsiPrimeDiMuon:
+     case kPiToMu:
+     case kKaToMu:
+     case kAllMuonic:
+     case kWToMuon:
+     case kWToCharm:
+     case kWToCharmToMuon:
+     case kZDiMuon:
+     case kZDiElectron:
+     case kHadronicD:
+     case kHadronicDWithout4Bodies:
+     case kPhiKK:
+     case kOmega:
+     case kLambda:      
+     case kNoDecay:
+     case kNoDecayHeavy:
+     case kNeutralPion:
+     case kBPsiPrimeDiElectron:
+     AliWarning(Form("Warning: case %s not implemented for this class!",(Char_t)decay));
+     break;
+     }
+     ReadDecayTable();
+  }
+
+Float_t AliDecayerEvtGen::GetPartialBranchingRatio(Int_t)
+  {
+   // This method is dummy
+   return  1.;
+  }
+
+Float_t AliDecayerEvtGen::GetLifetime(Int_t kf)
+  {    
+  //
+  //Input: pdg code of a particle   
+  //return lifetime in sec for a particle with particle code kf
+  //
+  EvtId IdPart=EvtPDL::evtIdFromStdHep(kf); 
+  Double_t lifetime = EvtPDL::getctau(IdPart); //c*tau (mm)
+  AliDebug(1,Form("lifetime is %f (mum) , particle id= %d",lifetime*1000,kf));
+  return lifetime*kconv; //tau (sec)
+  }
+
+void AliDecayerEvtGen::ReadDecayTable()
+    {
+     //Input none - Output none 
+     //Read the decay table that correspond to the path 
+     //fDecayTablePath
+     // 
+     TString temp = fDecayTablePath;
+     if(!temp.EndsWith("DECAY.DEC"))
+     fGenerator->readUDecay(fDecayTablePath);
+    }
+////////////////////////////////////////////////////////////
+Bool_t AliDecayerEvtGen::SetDecayTablePath(Char_t *path)
+  {
+  // 
+  //Set the path of the decay table read to force particle decays 
+  //
+  if(gSystem->AccessPathName(path)) 
+  {
+   AliWarning("Attention: This path not exist!\n");
+   return kFALSE;
+  }
+  fDecayTablePath = path;
+  return kTRUE;
+  } 
+
+
+

diff --git a/TEvtGen/AliDecayerEvtGen.h b/TEvtGen/AliDecayerEvtGen.h
new file mode 100644 (file)
index 0000000..c0c61fb
--- /dev/null
+++ b/TEvtGen/AliDecayerEvtGen.h
@@ -0,0 +1,60 @@
+#ifndef ALIDECAYEREVTGEN_H
+#define ALIDECAYEREVTGEN_H
+/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+///////////////////////////////////////////////////////////////////////////
+//  Implementation of AliDecayer using EvtGen package. It inherits       //
+//  from AliDecayer.                                                     //
+//                                                                       //
+//  Contact: Giuseppe.Bruno@ba.infn.it  &  Fiorella.Fionda@ba.infn.it    //
+///////////////////////////////////////////////////////////////////////////
+
+
+#include "AliDecayer.h"
+
+class EvtGen;
+class EvtRandomEngine;
+class EvtStdHep;
+class TLorentzVector;
+class TClonesArray;
+
+class AliDecayerEvtGen : public AliDecayer
+{
+ public:
+  AliDecayerEvtGen();
+  AliDecayerEvtGen(const AliDecayerEvtGen &decayer);
+  virtual ~AliDecayerEvtGen();
+  virtual void  Init();
+  virtual void  Decay(Int_t ipart, TLorentzVector *p);
+  virtual Int_t ImportParticles(TClonesArray *particles);
+  virtual void    SetForceDecay(Decay_t decay) {fDecay=decay;}
+  virtual void    SetForceDecay(Int_t decay){SetForceDecay((Decay_t) decay);}
+  virtual void    ForceDecay();
+  virtual Float_t GetPartialBranchingRatio(Int_t ipart);
+  virtual Float_t GetLifetime(Int_t kf);
+  Char_t*           GetDecayTablePath() {return fDecayTablePath;} 
+  virtual void    ReadDecayTable();
+  Bool_t SetDecayTablePath(Char_t *path);
+   
+  private:
+  void  Copy(TObject &decayer) const;
+  AliDecayerEvtGen &operator=(const AliDecayerEvtGen &decayer) 
+      {decayer.Copy(*this);return(*this);}
+
+  protected:
+  EvtRandomEngine *fRandomEngine;  //!pointer to EvtRandomEngine to generate random number   
+  EvtGen *fGenerator;              //!pointer to EvtGen class interface 
+  EvtStdHep *fEvtstdhep;           //!pointer to EvtGen common block
+  Char_t *fDecayTablePath;         //!pointer to decay table path
+  Char_t *fParticleTablePath;      //!pointer to particle table path
+  Decay_t     fDecay;                     // Forced decay case         
+
+  ClassDef(AliDecayerEvtGen,0)  //AliDecayer implementation using EvtGen 
+};
+
+//converts from mm/c to s
+const Double_t kconv=0.001/2.999792458e8;
+
+#endif
+
+ 

diff --git a/TEvtGen/AliGenEvtGen.cxx b/TEvtGen/AliGenEvtGen.cxx
new file mode 100644 (file)
index 0000000..08d2c58
--- /dev/null
+++ b/TEvtGen/AliGenEvtGen.cxx
@@ -0,0 +1,241 @@
+/**************************************************************************
+ * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+***************************************************************************/
+///////////////////////////////////////////////////////////////////////////
+//  Class to generate decays of particles generated by a                 //
+//  previous generator. It works as a generator, but pratically it       //
+//  performs only decays. It works with this scheme: first loops over    // 
+//  particles on the stack, selects those to be decayed, decays them     //
+//  and then pushes the decay products on the stack.                     //
+//                                                                       //
+//  Giuseppe E. Bruno            &    Fiorella Fionda                    //
+//  (Giuseppe.Bruno@ba.infn.it)       (Fiorella.Fionda@ba.infn.it)       //
+///////////////////////////////////////////////////////////////////////////
+
+#include "AliStack.h"
+#include "AliGenEvtGen.h"
+#include "AliRun.h"
+#include <TParticle.h>
+
+ClassImp(AliGenEvtGen)
+///////////////////////////////////////////////////////////////////////////
+AliGenEvtGen::AliGenEvtGen(): 
+  fStack(0x0),    
+  fDecayer(0x0),
+  fForceDecay(kAll),
+  fSwitchOff(kBeautyPart),
+  fUserDecay(kFALSE),
+  fUserDecayTablePath(0x0)
+  {
+  //
+  // Default Construction
+  //
+  }
+///////////////////////////////////////////////////////////////////////////////////////////
+AliGenEvtGen::~AliGenEvtGen() 
+  {
+  //
+  // Standard Destructor
+  //
+  if(fStack) {delete fStack;}
+  fStack = 0; 
+  if(fDecayer) {delete fDecayer;}
+  fDecayer = 0;
+  if(fUserDecayTablePath) {delete fUserDecayTablePath;}
+  fUserDecayTablePath = 0;
+  }
+///////////////////////////////////////////////////////////////////////////////////////////
+
+void AliGenEvtGen::Init() 
+  {
+  //
+  // Standard AliGenerator Initializer - no input 
+  // 1) initialize EvtGen with default decay and particle table
+  // 2) set the decay mode to force particle 
+  // 3) set a user decay table if defined
+  //
+  if(fDecayer)
+  {
+  AliWarning("AliGenEvtGen already initialized!!!");
+  return;
+  }
+  fDecayer = new AliDecayerEvtGen();
+  fDecayer->Init(); //read the default decay table DECAY.DEC and particle table
+
+  //if is set a decay mode: default decay mode is kAll  
+  fDecayer->SetForceDecay(fForceDecay);
+  fDecayer->ForceDecay();
+
+  //if is defined a user decay table
+  if(fUserDecay) 
+    {
+    fDecayer->SetDecayTablePath(fUserDecayTablePath);
+    fDecayer->ReadDecayTable();
+    }
+  }
+
+/////////////////////////////////////////////////////////////////////////////////////////////
+void AliGenEvtGen::Generate() 
+  {
+  //
+  //Generate method - Input none - Output none
+  //For each event:
+  //1)return the stack of the previous generator and select particles to be decayed by EvtGen
+  //2)decay particles selected and put the decay products on the stack
+  //
+  //
+  Float_t polar[3]= {0,0,0};  // Polarisation of daughter particles 
+  Float_t origin0[3];         // Origin of the parent particle 
+  Float_t pc[3], och[3]; // Momentum and origin of the children particles from EvtGen
+  Int_t nt;
+  Float_t tof;
+  Int_t nPrimsPythia;
+  TLorentzVector *mom=new TLorentzVector();
+  static TClonesArray *particles;
+  if(!particles) particles = new TClonesArray("TParticle",1000);
+  fStack = AliRunLoader::Instance()->Stack();
+  if(!fStack) {Info("Generate","Error: No stack found!"); return;}
+  nPrimsPythia = fStack->GetNprimary();  
+  AliDebug(1,Form("nPrimsPythia = %d \n",nPrimsPythia));
+  for (Int_t iTrack = 0; iTrack < nPrimsPythia; ++iTrack) {
+  TParticle *part = fStack->Particle(iTrack);
+  Int_t pdg=part->GetPdgCode();
+
+  AliDebug(1,Form("GetFlavour = %d e pdg = %d \n",GetFlavour(pdg),pdg));
+  
+  switch(fSwitchOff)
+  {
+  case kAllPart: 
+  break;
+  case kBeautyPart:
+  if(GetFlavour(pdg)!=5) continue;
+  break; 
+  case kCharmPart: 
+  if(GetFlavour(pdg)!=4) continue;
+  break;
+  } 
+
+  //check if particle is already decayed by Pythia  
+  if(part->GetStatusCode() != 1 || part->GetNDaughters()>0) 
+    {
+    Info("AliGenEvtGen","Attention: particle %d is already decayed by Pythia!",pdg); 
+    continue;
+    }
+  
+  part->SetStatusCode(11); //Set particle as decayed : change the status code
+  
+  mom->SetPxPyPzE(part->Px(),part->Py(),part->Pz(),part->Energy());
+  Int_t np;
+
+   do{
+     fDecayer->Decay(part->GetPdgCode(),mom);
+     np = fDecayer->ImportParticles(particles);
+     }while(np<0);
+
+  Int_t* trackIt    = new Int_t[np];
+  Int_t* pParent   = new Int_t[np];
+  AliDebug(1,Form("np = %d \n",np)); 
+  
+  for (int i = 0; i < np; i++) {
+             pParent[i]   = -1;
+             trackIt[i]   =  0;
+         }
+  //select trackable particle  
+  if (np >1) {
+  TParticle* iparticle =  (TParticle *) particles->At(0);//parent particle
+   for (int i = 1; i<np ; i++) { 
+     iparticle = (TParticle*) particles->At(i);
+     Int_t ks = iparticle->GetStatusCode();
+
+     //track last decay products
+     if(ks==1) trackIt[i]=1; 
+   
+     }//decay particles loop
+
+  }// if decay products   
+ 
+  origin0[0]=part->Vx(); //[cm]
+  origin0[1]=part->Vy(); //[cm]
+  origin0[2]=part->Vz(); //[cm]
+  //
+  // Put decay products on the stack
+  //
+  for (int i = 1; i < np; i++) {
+         TParticle* iparticle = (TParticle *) particles->At(i);
+         Int_t kf   = iparticle->GetPdgCode();
+         Int_t ksc  = iparticle->GetStatusCode();
+         Int_t jpa  = iparticle->GetFirstMother()-1; //jpa = 0 for daughters of beauty particles
+         Int_t iparent = (jpa > 0) ? pParent[jpa] : iTrack;      
+
+         och[0] = origin0[0]+iparticle->Vx()/10; //[cm]
+         och[1] = origin0[1]+iparticle->Vy()/10; //[cm]
+         och[2] = origin0[2]+iparticle->Vz()/10; //[cm]
+         pc[0]  = iparticle->Px(); //[GeV/c]
+         pc[1]  = iparticle->Py(); //[GeV/c]
+         pc[2]  = iparticle->Pz(); //[GeV/c]
+         tof = part->T()+kconv*iparticle->T(); 
+
+         AliDebug(1,Form("FirstMother = %d e indicePart = %d e pdg = %d \n",jpa,i,kf));         
+          
+         PushTrack(trackIt[i], iparent, kf, pc, och, polar,tof, kPDecay, nt, 1., ksc); 
+         if(trackIt[i]==1) AliDebug(1,Form("Trackable particles: %d e pdg %d \n",i,kf));
+         pParent[i] = nt;
+         KeepTrack(nt);
+         SetHighWaterMark(nt);
+         }// Particle loop
+         particles->Clear();
+        if (trackIt)    delete[] trackIt;
+         if (pParent)    delete[] pParent;
+       }
+  Info("Generate","AliGenEvtGen DONE");  
+  }
+
+//////////////////////////////////////////////////////////////////////////////////////////
+Int_t AliGenEvtGen::GetFlavour(Int_t pdgCode)  
+  { 
+  //
+  // return the flavour of a particle
+  // input: pdg code of the particle
+  // output: Int_t 
+  //         3 in case of strange (open and hidden)
+  //         4 in case of charm (")
+  //         5 in case of beauty (")
+  //
+  Int_t pdg = TMath::Abs(pdgCode);
+  //Resonance
+  if (pdg > 100000) pdg %= 100000;
+  if(pdg > 10000)  pdg %= 10000;
+  // meson ?
+  if(pdg > 10) pdg/=100;
+  // baryon ?
+  if(pdg > 10) pdg/=10;
+  return pdg;
+  }
+
+/////////////////////////////////////////////////////////////////////////////////////////
+Bool_t AliGenEvtGen::SetUserDecayTable(Char_t *path)
+  {
+  //
+  //Set the path of user decay table if it is defined
+  //
+  //put a comment to control if path exists 
+  if(gSystem->AccessPathName(path))
+    {
+    AliWarning("Attention: This path not exist!\n");
+    return kFALSE;
+    }
+  fUserDecayTablePath = path;
+  fUserDecay = kTRUE;
+  return kTRUE;
+  }

diff --git a/TEvtGen/AliGenEvtGen.h b/TEvtGen/AliGenEvtGen.h
new file mode 100644 (file)
index 0000000..a3d1e7c
--- /dev/null
+++ b/TEvtGen/AliGenEvtGen.h
@@ -0,0 +1,58 @@
+#ifndef ALIGENEVTGEN_H
+#define ALIGENEVTGEN_H
+
+/* Copyright(c) 2007-2009, ALICE Experiment at CERN, All rights reserved. *
+ * See cxx source for full Copyright notice                               */
+///////////////////////////////////////////////////////////////////////////
+//                                                                       //
+//  AliGenEvtGen class to performs decays of particles generated by a    //
+//  previous generator. It inherits from AliGenerator.                   //
+//                                                                       //  
+//  Origin: Giuseppe.Bruno@ba.infn.it  &  Fiorella.Fionda@ba.infn.it     //
+///////////////////////////////////////////////////////////////////////////
+
+#include "AliGenerator.h"
+#include "AliDecayerEvtGen.h"
+
+class TParticle;
+
+class AliGenEvtGen : public AliGenerator {
+
+  public:
+
+  typedef enum {kAllPart, kBeautyPart, kCharmPart} DecayOff_t;
+
+  AliGenEvtGen();
+ 
+  ~AliGenEvtGen();
+
+  virtual void Init();
+  virtual void Generate(); 
+  AliDecayer* GetDecayer() {return fDecayer;}
+  void SetForceDecay(Decay_t decay = kAll) {fForceDecay = decay;} //set a decay mode
+  Bool_t SetUserDecayTable(Char_t *path);
+  void SetParticleSwitchedOff(DecayOff_t decay) {fSwitchOff = decay;} //set particles to be decayed by EvtGen
+  
+  protected:
+  Int_t GetFlavour(Int_t pdgCode); 
+  
+ 
+  private:
+  AliStack *fStack;               //!pointer to AliStack
+  AliDecayerEvtGen  *fDecayer;     //!pointer to AliDecayerEvtGen
+  Decay_t  fForceDecay;            //!decay case
+  DecayOff_t fSwitchOff;           //!particle decay  
+  Bool_t fUserDecay;              //!TRUE if a user decay table is set 
+  Char_t    *fUserDecayTablePath;  //!pointer to path of user decay table
+
+  AliGenEvtGen(const AliGenEvtGen &EvtGen);
+  AliGenEvtGen & operator=(const AliGenEvtGen & rhs);
+  
+  
+
+  ClassDef(AliGenEvtGen,0) //AliGenerator interface to EvtGen
+};
+
+////////////////////////////////////////////////////////////////////////////////////////////////
+
+#endif

diff --git a/TEvtGen/EvtGen/DECAY.DEC b/TEvtGen/EvtGen/DECAY.DEC
new file mode 100644 (file)
index 0000000..4c85113
--- /dev/null
+++ b/TEvtGen/EvtGen/DECAY.DEC
@@ -0,0 +1,10239 @@
+# $Id: DECAY.DEC,v 1.271 2009/02/15 18:43:38 ryd Exp $
+#
+#Define the B0B0bar mass difference
+Define dm 0.489e12
+Define dgamma 0
+Define qoverp 1
+Define phaseqoverp 0
+# define the values of the CKM angles  (alpha=70, beta=40)
+Define alpha 1.365
+Define beta  0.39
+Define gamma 1.387
+Define twoBetaPlusGamma  2.167
+Define betaPlusHalfGamma 1.0835
+Define minusGamma -1.387
+Define minusTwoBeta -0.78
+
+# New definitions for psiKstar modes (Lange, July 26, 2000)
+Define PKHplus  0.159
+Define PKHzero  0.775
+Define PKHminus 0.612
+Define PKphHplus  1.563
+Define PKphHzero  0.0
+Define PKphHminus 2.712
+
+Define Aplus  0.490
+Define Azero  1.10
+Define Aminus 0.4
+#
+Define phAplus  2.5
+Define phAzero  0.0
+Define phAminus -0.17
+
+### These are for incoherent B_s and B_d mixing
+Define qoverp_incohMix_B_s0 1.0
+Define dm_incohMix_B_s0 17.1e12
+Define qoverp_incohMix_B0 1.0
+Define dm_incohMix_B0 0.489e12
+Define qoverp_incohMix_D0 1.0
+Define dm_incohMix_D0 0.0
+
+#
+# These particle aliases are used in for CP violating decays
+# in which the decay distributions depends on how the K* decayed.
+#
+#
+Alias K*L       K*0
+Alias K*S       K*0
+Alias K*BL       anti-K*0
+Alias K*BS       anti-K*0
+Alias K*0T      K*0
+Alias anti-K*0T anti-K*0
+Alias K*BR  anti-K*0
+Alias K*0R  K*0
+Alias anti-K_0*0N anti-K_0*0
+Alias K_0*0N K_0*0
+#
+ChargeConj K*L    K*BL
+ChargeConj K*S    K*BS
+ChargeConj K*0T    anti-K*0T
+ChargeConj K_0*0N  anti-K_0*0N
+ChargeConj K*0R    K*BR    
+#
+#JetSet parameter modifications
+#
+# O.Long : March 30, 2000
+#  Change from using Peterson fragmentation to Lund fragmentation
+#  with the default parameters (A=0.30, B=0.58)
+#
+#(Very important that there are no blank spaces in the parameter string!)
+# Number of flavor produced. 5 = udscb production
+JetSetPar MSTJ(104)=5
+# turn on initial state radiation (was off before)
+ #JetSetPar MSTJ(107)=1
+#use Lund fragmentation with default parameters
+JetSetPar MSTJ(11)=4
+# Lund symmetric fragmentation function parameter a 
+JetSetPar PARJ(41)=.30
+# Lund symmetric fragmentation function parameter b 
+JetSetPar PARJ(42)=.58
+# Lund symmetric fragmentation function parameter a 
+JetSetPar PARJ(43)=.50
+# Lund symmetric fragmentation function parameter b 
+JetSetPar PARJ(44)=.90
+#eps_c   from CLEO 8/97
+ #JetSetPar PARJ(54)=-.040
+# eps_b
+ #JetSetPar PARJ(55)=-.004
+#Turn of B0-B0B mixing in JetSet:
+ #JetSetPar MSTJ(26)=0
+# D* production - Lange July 19, 2002
+#050214 JetSetPar PARJ(13)=0.45
+ #JetSetPar PARJ(13)=0.4
+# eta'
+# lange - change eta' production to agree better with udsc results
+# from charmless group
+#050214 JetSetPar PARJ(26)=0.15
+ #JetSetPar PARJ(26)=0.1
+# lange - change eta production to agree better with udsc results
+# from charmless group
+#050214 JetSetPar PARJ(25)=1.0
+ #JetSetPar PARJ(25)=0.6
+# P(ssbar)/P(uubar) - has a big effect on phis. Turn down to 
+# reduce phi production in continuum.
+#050214 was no change
+ #JetSetPar PARJ(2)=0.26
+#050214JetSetPar PARJ(12)=0.2
+ #JetSetPar PARJ(12)=0.55
+ #JetSetPar PARJ(11)=0.4
+#lange/Petersen Nov 14, 2004 baryon production turned down for sp8
+ #JetSetPar PARJ(1)=0.065
+# control of L=1 mesons (in order: J1S0 J0S1 J1S1 J2S1
+ #JetSetPar PARJ(14)=0.05
+ #JetSetPar PARJ(15)=0.05
+ #JetSetPar PARJ(16)=0.05
+ #JetSetPar PARJ(17)=0.05
+#cut-off parameter used to stop fragmentation process (was .4, CLEO uses.3)
+ #JetSetPar PARJ(33)=0.3
+#
+### Parameters below are a new tuning of the continuum
+### fragmentation that softens the pion spectrum
+### but commented out untill verified to be physically
+### sound. ryd 00-02-10
+### Lund symmetric fragmentation function parameter a 
+##JetSetPar PARJ(41)=.40
+### Lund symmetric fragmentation function parameter b 
+##JetSetPar PARJ(42)=.29
+### Lund symmetric fragmentation function parameter a 
+##JetSetPar PARJ(43)=.60
+### Lund symmetric fragmentation function parameter b 
+##JetSetPar PARJ(44)=.45
+###eps_c   from CLEO 8/97
+##JetSetPar PARJ(54)=-.040
+##JetSetPar PARJ(46)=0.0
+##JetSetPar PARJ(33)=0.2
+###End of new continuum tuning
+#
+#pythia stuff for standalone generation
+# no isr
+ #JetSetPar MSTP(11)=0
+# take the qqbar energy from EvtGen
+ #JetSetPar MSTP(171)=1
+ #JetSetPar MSTP(172)=1
+
+
+# Turn on PHOTOS for all decays
+# Lange Aug19,2002
+yesPhotos
+#
+
+ModelAlias VubHybridB0 VUBHYBRID 4.62 2.27 0.22 8 8 8
+    0.28 1.40 1.60 1.80  2.00  2.50  3.00  3.50 
+    0.00 2.50 5.00 7.50 10.00 12.50 15.00 20.00 
+    0.00 0.50 1.00 1.25  1.50  1.75  2.00  2.25
+    0.00000 1.2059 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241
+    0.09098 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.47335 1.21021 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241
+    0.52695 1.21011 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241
+    0.6304 1.21096 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.33327 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.54178 1.20886 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241
+    0.59143 1.20802 1.21241 1.21241 1.21241 1.21241 1.21241 0.99995
+    0.69332 1.20698 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.66811 1.20755 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.34624 1.21241 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.23645 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.53769 1.21089 1.21241 1.21241 1.21241 1.21241 1.21241 1.21241
+    0.63768 1.20843 1.21241 1.21241 1.21241 1.21241 1.21241 0.99995
+    0.71713 1.20869 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.75566 1.20778 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.74831 1.20916 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.21312 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.62594 1.20926 1.21241 1.21241 1.21241 1.21241 1.21241 0.99995
+    0.5902 1.20516 1.21241 1.21241 1.21241 1.21241 1.21241 0.99995
+    0.68809 1.20016 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.78229 1.20459 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.82156 1.20561 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.72422 1.20485 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.00423 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.55708 1.20977 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.61225 1.20629 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.66904 1.20049 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.72567 1.20029 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.78501 1.19808 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.83037 1.19609 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.59537 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.00000 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.45962 1.20645 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.53276 1.20083 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995
+    0.57906 1.19669 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.68263 1.19458 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.74222 1.19771 1.21241 1.21241 0.99995 0.99995 0.99995 0.99995
+    0.78307 1.19479 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.58714 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.19489 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.24451 1.21241 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.14704 1.20472 1.21241 1.21241 1.21241 0.99995 0.99995 0.99995
+    0.42163 1.19849 1.21241 1.21241 0.99995 0.99995 0.99995 0.99995
+    0.37319 1.18877 1.21241 1.21241 0.99995 0.99995 0.99995 0.99995
+    0.38655 1.1995 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.34297 1.21241 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.31299 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.57067 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995
+    0.57067 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995 0.99995;
+
+# Hybrid weights for charged B -> Xu l nu decays
+ModelAlias VubHybridBp VUBHYBRID  4.62 2.27 0.22 8 8 8 
+    0.28 1.40 1.60 1.80  2.00  2.50  3.00  3.50 
+    0.00 2.50 5.00 7.50 10.00 12.50 15.00 20.00 
+    0.00 0.50 1.00 1.25  1.50  1.75  2.00  2.25
+    0.00000 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.00000 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.11287 1.26974 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.10819 1.26929 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.27133 1.26671 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.17256 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.18216 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.32091 1.26953 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.52924 1.26988 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.51586 1.26755 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.26794 1.27029 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.21166 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029 1.27029
+    0.29627 1.26757 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.52168 1.26788 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.5675 1.26647 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.62174 1.26467 1.27029 1.27029 0.9986 0.9986 0.9986 0.9986
+    0.00536 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.38235 1.2678 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.44511 1.26686 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.54476 1.26672 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.66852 1.26707 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.64859 1.26542 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.53664 1.2536 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.00000 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.22558 1.26822 1.27029 1.27029 1.27029 1.27029 1.27029 0.9986
+    0.38409 1.26531 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.52891 1.26284 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.65141 1.26648 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.70896 1.26227 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.79827 1.26101 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.4274 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.08636 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.28546 1.26121 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.40621 1.26584 1.27029 1.27029 1.27029 1.27029 0.9986 0.9986
+    0.58209 1.26637 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.61119 1.26424 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.70714 1.26297 1.27029 1.27029 0.9986 0.9986 0.9986 0.9986
+    0.73926 1.25781 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.53898 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.52061 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.00000 1.26726 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.16996 1.2649 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.23334 1.2534 1.27029 1.27029 1.27029 0.9986 0.9986 0.9986
+    0.25078 1.26036 1.27029 1.27029 0.9986 0.9986 0.9986 0.9986
+    0.27805 1.25917 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.27337 1.27029 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.18235 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.84654 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986
+    0.84654 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986 0.9986;
+
+#
+
+
+# Define model alias for the D -> K*lnu decays
+# This should correspond to MV=2.1 GeV, MA1=MA2=2.5 GeV,
+# R2=0.85 and RV=1.5
+#
+ModelAlias SLPOLE_DtoKstarlnu SLPOLE  1.00 -0.558 0.0 1.0
+                                      0.85 -0.558 0.0 1.0
+                                      1.50 -0.790 0.0 1.0
+                                      0.00 -0.558 0.0 1.0;
+#
+# Model alias for the D -> rho l nu
+# MV=2.0 GeV and MA1=MA2=2.4 GeV, R2=0.6, RV=1.5
+#
+ModelAlias SLPOLE_Dtorholnu SLPOLE 1.00 -0.605 0.0 1.0
+                                   0.60 -0.605 0.0 1.0
+                                   1.50 -0.871 0.0 1.0
+                                   0.00 -0.605 0.0 1.0;
+#
+# Model alias for the D -> omega l nu
+# MV=2.0 GeV and MA1=MA2=2.4 GeV, R2=0.6, RV=1.5
+#
+ModelAlias SLPOLE_Dtoomegalnu SLPOLE 1.00 -0.605 0.0 1.0
+                                     0.60 -0.605 0.0 1.0
+                                     1.50 -0.871 0.0 1.0
+                                     0.00 -0.605 0.0 1.0;
+#
+# Model aliases for D -> Klnu, pilenu, etalnu, eta'lnu
+#
+ModelAlias SLBKPOLE_DtoKlnu SLBKPOLE 1.0 0.25 1.0 2.112;
+ModelAlias SLBKPOLE_Dtopilnu SLBKPOLE 1.0 0.25 1.0 2.010;
+ModelAlias SLBKPOLE_Dtoetalnu SLBKPOLE 1.0 0.25 1.0 2.010;
+ModelAlias SLBKPOLE_Dtoetaplnu SLBKPOLE 1.0 0.25 1.0 2.010;
+
+
+#Decay vpho
+#1.000                     PYCONT;
+#Enddecay
+#
+# use new VSS_BMIX mixing decay model (DK,28-Oct-1999)
+Decay Upsilon(4S)
+0.508200 B+  B-                          VSS;
+0.490000 B0  anti-B0                     VSS_BMIX dm;
+0.000016    e+    e-              PHOTOS  VLL;
+0.000016    mu+   mu-             PHOTOS  VLL;
+0.000016    tau+  tau-            PHOTOS  VLL;
+0.000090         Upsilon(2S)  pi+  pi-   VVPIPI;
+0.000045         Upsilon(2S)  pi0  pi0   VVPIPI;
+0.000090         Upsilon  pi+  pi-   VVPIPI;
+0.000045         Upsilon  pi0  pi0   VVPIPI;
+0.000200        Upsilon  eta  PARTWAVE 0.0 0.0 1.0 0.0 0.0 0.0;
+
+# BF ~ (2J+1)E^3_gamma; see PRL 94, 032001
+# V-> gamma S    Partial wave (L,S)=(0,0)
+0.000094         gamma chi_b0(3P)  HELAMP 1. 0. 1. 0.; 
+# V-> gamma V    Partial wave (L,S)=(0,1)
+0.000141         gamma chi_b1(3P)  HELAMP 1. 0. 1. 0. -1. 0. -1. 0.;
+# V-> gamma T    Partial wave (L,S)=(0,1)
+0.000131         gamma chi_b2(3P)  HELAMP 2.4494897 0. 1.7320508 0. 1. 0.
+                                   1. 0. 1.7320508 0. 2.4494897 0.;
+# V-> gamma S    Partial wave (L,S)=(0,0)
+0.000003         gamma chi_b0(2P)  HELAMP 1. 0. 1. 0.; 
+# V-> gamma V    Partial wave (L,S)=(0,1)
+0.000008         gamma chi_b1(2P)  HELAMP 1. 0. 1. 0. -1. 0. -1. 0.;
+# V-> gamma T    Partial wave (L,S)=(0,1)
+0.000012         gamma chi_b2(2P)  HELAMP 2.4494897 0. 1.7320508 0. 1. 0.
+                                   1. 0. 1.7320508 0. 2.4494897 0.;
+0.000850         g       g       g       PYTHIA 4;
+0.000040         gamma   g       g       PYTHIA 4;
+Enddecay
+#
+#
+#
+Decay anti-B0
+#                   b -> c semileptonic
+#
+0.0529   D*+    e-   anti-nu_e          PHOTOS  HQET 0.77 1.33 0.92;
+0.0208   D+     e-   anti-nu_e          PHOTOS  ISGW2;
+0.0052   D_1+   e-   anti-nu_e          PHOTOS  ISGW2;
+0.0045   D_0*+   e-   anti-nu_e          PHOTOS  ISGW2;
+0.0083   D'_1+   e-   anti-nu_e          PHOTOS  ISGW2;
+0.0023   D_2*+   e-   anti-nu_e          PHOTOS  ISGW2;
+0.0003   D*+  pi0   e-   anti-nu_e      PHOTOS  GOITY_ROBERTS;
+0.0065   D*0  pi+   e-   anti-nu_e      PHOTOS  GOITY_ROBERTS;
+0.0010   D+   pi0   e-   anti-nu_e      PHOTOS  GOITY_ROBERTS;
+0.0032   D0   pi+   e-   anti-nu_e      PHOTOS  GOITY_ROBERTS;
+
+0.0529   D*+    mu-  anti-nu_mu         PHOTOS  HQET 0.77 1.33 0.92;
+0.0208   D+     mu-  anti-nu_mu         PHOTOS  ISGW2;
+0.0052   D_1+   mu-  anti-nu_mu         PHOTOS  ISGW2;
+0.0045   D_0*+   mu-  anti-nu_mu         PHOTOS  ISGW2;
+0.0083   D'_1+   mu-  anti-nu_mu         PHOTOS  ISGW2;
+0.0023   D_2*+   mu-  anti-nu_mu         PHOTOS  ISGW2;
+0.0003   D*+  pi0   mu-  anti-nu_mu    PHOTOS   GOITY_ROBERTS;
+0.0065   D*0  pi+   mu-  anti-nu_mu    PHOTOS   GOITY_ROBERTS;
+0.0010   D+   pi0   mu-  anti-nu_mu    PHOTOS   GOITY_ROBERTS;
+0.0032   D0   pi+   mu-  anti-nu_mu    PHOTOS   GOITY_ROBERTS;
+#
+#                    b -> c tau nu
+#
+0.0202   D*+    tau- anti-nu_tau                 ISGW2;
+0.0088   D+     tau- anti-nu_tau                 ISGW2;
+0.0013   D_1+   tau- anti-nu_tau                 ISGW2;
+0.0013   D_0*+   tau- anti-nu_tau                 ISGW2;
+0.0020   D'_1+   tau- anti-nu_tau                 ISGW2;
+0.0020   D_2*+   tau- anti-nu_tau                 ISGW2;
+#
+#                    b -> u l nu
+#
+# 03/28/01 S.Menke included inclusive B -> X_u l nu decays for m(X_u) > 1.26 GeV.
+#                  with some thresholds to gradually switch from excl. to incl.
+#                  modeling without big jumps keeping the better known excl.
+#                  Brs at their actual values and constrain the incl. Br and 
+#                  weights in a way that ~40% of the total BR belongs to 
+#                  the mass region > 1.5 GeV
+# October 20, 2004 - Lange - update using sl awg input
+# March 30, 2005 J. Dingfelder: Updated  B -> X_u l nu hybrid model 
+#                               (supersedes older version by S. Menke). Inclusive
+#                               decays are reweighted in three kinematic 
+#                               variables mX, q2, and El (=> 512 weights total).
+#                               The list of weights is included as ModelAlias 
+#                               once per B charge state.
+# NOTE: - The specified parameters for model VUBHYBRID are:
+#         mb, a, alpha_s, nbins_mX, nbins_q2, nbins_El, 
+#         list of threshold values for mX, q2, El in ascending order.
+#
+#         *** DO NOT CHANGE ANY OF THESE PARAMETERS! *** 
+#         Otherwise a new set of weights would be needed, since they were 
+#         generated for a certain set of values for mb, a, excl. BFs, non-res BF.
+#
+#       - If no binning (nbins or thresholds) are specified after the first 
+#         three parameters (mb,a,alpha_s) no hybrid weighting is performed.
+
+# --- Hybrid weights for neutral B -> Xu l nu decays
+# NOTE: Do NOT CHANGE any BFs without using the corresponding set of 
+#       hybrid weights (and vice versa).
+     
+0.000136   pi+    e-   anti-nu_e        PHOTOS   ISGW2;
+0.000220   rho+   e-   anti-nu_e        PHOTOS   ISGW2;
+0.000000   D(2S)+ e-   anti-nu_e        PHOTOS   ISGW2;
+0.000000   D*(2S)+ e-   anti-nu_e       PHOTOS   ISGW2;
+0.001892   Xu+ e-  anti-nu_e            PHOTOS VubHybridB0;
+#
+0.000136   pi+    mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000220   rho+   mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000000   D(2S)+ mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000000   D*(2S)+ mu-  anti-nu_mu      PHOTOS    ISGW2;
+0.001892   Xu+ mu-  anti-nu_mu          PHOTOS VubHybridB0;
+#
+0.000060   pi+    tau- anti-nu_tau                 ISGW2;
+0.000083   rho+   tau- anti-nu_tau                 ISGW2;
+0.000090   a_1+    tau- anti-nu_tau                 ISGW2;
+0.000008   b_1+    tau- anti-nu_tau                 ISGW2;
+0.000008   a_0+    tau- anti-nu_tau                 ISGW2;
+0.000008   a_2+    tau- anti-nu_tau                 ISGW2;
+0.000000   rho(2S)+ tau- anti-nu_tau                 ISGW2;
+0.000000   D(2S)+ tau- anti-nu_tau                 ISGW2;
+0.000000   D*(2S)+ tau- anti-nu_tau                 ISGW2;
+#
+
+#
+# b->u hadronic 
+# Ref. [B1]:
+# Lange Nov 14, 2004  (flip D_s K -- is opposite of D_s pi)
+0.000015   D_s-  pi+   PHSP;
+0.000028   D_s+  K-    PHSP;
+0.000029   D_s*-  pi+  SVS;
+0.000022   D_s*+  K-   SVS;
+#
+0.000016   rho+   D_s-    SVS;
+0.000016   K*-   D_s+     SVS;
+0.000016   D_s*-  rho+    SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+0.000016   D_s*+  K*-     SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+
+
+#
+# b -> s gamma
+#
+0.0000401   anti-K*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000135  anti-K_10  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000065  anti-K'_10 gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000128  anti-K'*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000166  anti-K_2*0 gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000017  anti-K''*0 gamma                   HELAMP 1.0 0.0 1.0 0.0;
+0.0003139  anti-Xsd    gamma  BTOXSGAMMA 2;
+#
+0.000000134 anti-K0   e+    e-    PHOTOS BTOSLLBALL;
+0.000001045 anti-K*0  e+    e-    PHOTOS BTOSLLBALL;
+0.000004994 anti-Xsd  e+    e-    PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000574 anti-K0   mu+   mu-   PHOTOS BTOSLLBALL;
+0.000001100 anti-K*0  mu+   mu-   PHOTOS BTOSLLBALL;
+0.000002506 anti-Xsd  mu+   mu-   PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000130 anti-K0   tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000190 anti-K*0  tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000200 anti-Xsd  tau+  tau-  PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+#
+#------------------------------------------------------------------------------
+# 2- and 3-body modes revised Feb.2005  Markus Cristinziani, SLAC
+#vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
+#    B0 -> CP eigenstates  (some exclusive b -> u)   sum=0.00036
+#
+0.0000049 pi+  pi-     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000015 pi0  pi0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000015 pi0  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  pi0  f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+# pi0 rho0 is with the 3-body modes
+0.000001  pi0  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00  a_00   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  a_00  f_0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  a_00 rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000003 f_0  f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  f_0  rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+
+# Penguin dominated modes are sin2beta, not sin2alpha
+0.0000043 phi  K_S0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0000043 phi  K_L0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000001  eta K_S0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000001  eta K_L0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000032  eta' K_S0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000032  eta' K_L0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.0000023 omega K_S0   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0000023 omega K_L0   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+#don't allow jetset to simulate the same decay but to a K0 or anti-K0
+#Lange 1/30/2003
+
+0.0  phi  anti-K0    PHSP;
+0.0  phi  K0    PHSP;
+
+0.0  eta anti-K0 PHSP;
+0.0  eta K0     PHSP;
+
+0.0  eta' anti-K0 PHSP;
+0.0  eta' K0    PHSP;
+
+0.0  omega anti-K0 PHSP;
+0.0  omega K0   PHSP;
+# phipi0 has no simple Pengun or tree so who knows - ditto phieta/eta'?
+0.00000005 phi  pi0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00000003 phi  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00000003 phi  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+#
+0.0000163 anti-K*0 eta SVS;
+0.0000038 anti-K*0 eta' SVS;
+0.000001  anti-K*0 omega SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi  omega   SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  omega omega  SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  omega rho0   SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+#
+#
+#2-body modes from Jim Olsen
+0.0000188 K-  pi+      PHSP;
+0.0000001  K+  K-       PHSP;
+# K_S0 K_L0 is 0
+# see http://babar-hn.slac.stanford.edu:5090/HyperNews/get/pubboard14/156.html
+0.0000006 K_S0 K_S0    PHSP;
+0.0000000 K_S0 K_L0    PHSP;
+0.0000006 K_L0 K_L0    PHSP;
+0.0000115 anti-K0 pi0  PHSP;
+# Put in with 0 BR to keep PYTHIA from generating it
+0.0000000 K0 anti-K0   PHSP;
+#
+# 3-body by John Back (jback@slac.stanford.edu) - Oct 15, 2002
+# JGS intersperses modes with pi0->eta,eta'
+# anti-B0 modes
+# Model generates rho(770), rho(1450) and rho(1700) resonances 
+# with interference (uses CKM angle alpha)
+0.000024     pi- pi+ pi0            BTO3PI_CP dm  alpha;
+# Put in with 0 BR to keep PYTHIA from generating these (see BTO3PI_CP above)
+0.000000     pi0 rho0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000000     rho- pi+               SVS;
+0.000000     rho+ pi-               SVS;
+0.000000     rho(2S)- pi+           SVS;
+0.000000     rho(2S)+ pi-           SVS;
+0.000000     rho(3S)- pi+           SVS;
+0.000000     rho(3S)+ pi-           SVS;
+#
+# rho0 3-body modes
+0.0000022    rho0 K_S0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.0000022    rho0 K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+#
+# rho+ 3-body modes
+0.0000085    rho+ K-                SVS;
+#
+# rho(1450) 3-body modes
+0.0000005    rho(2S)0 K_S0          SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.0000005    rho(2S)0 K_L0          SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+#
+#
+# f0,a0(980) 3-body modes
+0.0000055    f_0 anti-K0            PHSP;
+0.000001     a_00 anti-K0           PHSP;
+0.000001    a_0+ K-                PHSP;
+
+0.000001     a_0- pi+               PHSP;
+0.000003     a_0+ pi-               PHSP;
+#
+# K*0(892) 3-body modes
+0.0000017    anti-K*0 pi0           SVS;
+0.000003     anti-K*0 K0            SVS;
+#
+# K*0(1430) 3-body modes
+0.000001     anti-K_0*0 pi0         PHSP;
+0.000003     anti-K_0*0 K0          PHSP;
+#
+# K*0(1680) 3-body modes
+0.000001     anti-K''*0 pi0         SVS;
+0.000001     anti-K''*0 K0          SVS;
+
+#
+# K*-(892) 3-body modes
+0.0000085    K*- pi+                SVS;
+0.0000002     K*- K+                 SVS;
+#
+# K*-1(430) 3-body modes
+0.000033     K_0*- pi+              PHSP;
+0.0000002     K_0*- K+               PHSP;
+#
+# K*0(1680) 3-body modes
+0.000001     K''*- pi+                     SVS;
+0.0000001     K''*- K+               SVS;
+#
+#
+# Non-resonant 3-body left-overs
+# Most modes are set at 1e-6
+#
+# pi-pi+pi0: high mass + f0(400-1200) = (1+1)e-6
+####0402270.000002     pi- pi+ pi0            PHSP;
+0.000017     pi- pi+ eta            PHSP;
+0.000001     pi- pi+ eta'           PHSP;
+#
+# pi- pi+ anti-K0: 
+0.000013     pi- pi+ anti-K0        PHSP;
+#
+# K- pi+ pi0
+0.000002     K- pi+ pi0             PHSP;
+0.000020     K- pi+ eta             PHSP;
+0.000001     K- pi+ eta'            PHSP;
+#
+# pi- K+ K0. 
+# Both B0 and B0bar can decay to these => 0.5e-6 each
+0.0000005    pi- K+ anti-K0         PHSP;
+0.0000005    pi+ K- K0              PHSP;
+#
+# K- K+ anti-K0: total - phiKs = 25 -
+0.0000205    K- K+ anti-K0          PHSP;
+#
+# K- K+ pi0
+0.000001     K- K+ pi0              PHSP;
+0.000001     K- K+ eta              PHSP;
+0.000001     K- K+ eta'             PHSP;
+#
+# pi0 anti-K0 K0
+0.000002     anti-K0 K0 pi0        PHSP;
+0.000001     anti-K0 K0 eta        PHSP;
+0.000001     anti-K0 K0 eta'       PHSP;
+#
+# pi0 pi0 anti-K0
+0.000002     anti-K0 pi0 pi0      PHSP;
+0.000001     anti-K0 pi0 eta      PHSP;
+0.000001     anti-K0 pi0 eta'     PHSP;
+0.000001     anti-K0 eta eta      PHSP;
+0.000001     anti-K0 eta eta'     PHSP;
+0.000001     anti-K0 eta' eta'    PHSP;
+#
+# anti-K0 K0 anti-K0:  8*BF(KsKsKs)~=48 (8 is a guess)
+0.000048    anti-K0 K0 anti-K0     PHSP;
+#
+# pi0 pi0 pi0
+0.000001     pi0 pi0 pi0            PHSP;
+0.000001     pi0 pi0 eta            PHSP;
+0.000001     pi0 pi0 eta'           PHSP;
+0.000001     pi0 eta eta            PHSP;
+0.000001     pi0 eta eta'           PHSP;
+#
+#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# 2- and 3-body modes revised Feb.2005   Markus Cristinziani, SLAC
+#---------------------------------------------------------------------------
+#
+#4-body modes from Andrei Gritsan
+#--- 4-body rho-rho, rho-pi-pi, pi-pi-pi-pi --------------------------
+0.0000011 rho0 rho0              SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000025  rho+ rho-              SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000005  rho0 pi+ pi-           PHSP;
+0.000003  omega pi+ pi-          PHSP;
+0.000001  rho0 pi0 pi0           PHSP;
+0.000001  omega pi0 pi0          PHSP;
+0.0000005 rho0 pi0 eta           PHSP;
+0.0000002 rho0 eta eta           PHSP;
+0.0000002 rho0 pi0 eta'          PHSP;
+0.0000001 rho0 eta eta'          PHSP;
+0.0000005 omega pi0 eta          PHSP;
+0.0000002 omega eta eta          PHSP;
+0.0000002 omega pi0 eta'         PHSP;
+0.0000001 omega eta eta'         PHSP;
+
+0.000010  rho- pi+ pi0           PHSP;
+0.000005  rho- pi+ eta           PHSP;
+0.000010  rho+ pi- pi0           PHSP;
+0.000005  rho+ pi- eta           PHSP;
+0.000010  pi+ pi- pi+ pi-        PHSP;
+0.000010  pi+ pi- pi0 pi0        PHSP;
+0.000005  pi+ pi- eta pi0        PHSP;
+0.000002  pi+ pi- eta eta        PHSP;
+0.000002  pi+ pi- eta' pi0       PHSP;
+0.000001  pi+ pi- eta' eta       PHSP;
+0.0000165 a_1- pi+               SVS;
+0.0000165 a_1+ pi-               SVS;
+0.000000  b_1- pi+               SVS;
+0.000010  b_1+ pi-               SVS;
+0.000002  rho+ a_0-              SVS;
+0.000010  rho- a_0+              SVS;
+
+0.000001  f_0 pi+ pi-            PHSP;
+0.000001  f_0 pi0 pi0            PHSP;
+0.000001  a_0- pi+ pi0           PHSP;
+0.000001  a_00 pi0 pi0           PHSP;
+0.000001  a_0+ pi- pi0           PHSP;
+#--- 4-body rho-K*, rho-pi-K, K*-pi-pi, pi-pi-pi-K -------------------
+0.0000056 rho0  anti-K*0         SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  rho+  K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  pi- pi+ anti-K*0       PHSP;
+0.000010  pi+ pi0 K*-            PHSP;
+0.000005  pi+ eta K*-            PHSP;
+0.000002  pi+ eta' K*-           PHSP;
+0.000010  rho+  K- pi0           PHSP;
+0.000005  rho+  K- eta           PHSP;
+0.000002  rho+  K- eta'          PHSP;
+0.000010  rho+  anti-K0 pi-      PHSP;
+0.000010  rho-  anti-K0 pi+      PHSP;
+0.000005  rho0  anti-K0 pi0      PHSP;
+0.000002  rho0  anti-K0 eta      PHSP;
+0.000001  rho0  anti-K0 eta'     PHSP;
+0.000010  rho0  K- pi+           PHSP;
+0.000005  omega anti-K0 pi0      PHSP;
+0.000002  omega anti-K0 eta      PHSP;
+0.000001  omega anti-K0 eta'     PHSP;
+0.000005  omega K- pi+           PHSP;
+0.000010  pi+ pi0 pi- anti-K0    PHSP;
+0.000005  pi+ eta pi- anti-K0    PHSP;
+0.000010  pi0 pi0 pi0 anti-K0    PHSP;
+0.000005  pi0 pi0 eta anti-K0    PHSP;
+0.000002  pi0 eta eta anti-K0    PHSP;
+0.000002  pi0 pi0 eta' anti-K0   PHSP;
+0.000001  pi0 eta eta' anti-K0   PHSP;
+0.000010  pi+ pi- pi+ K-         PHSP;
+0.000010  pi0 pi0 pi+ K-         PHSP;
+0.000005  pi0 eta pi+ K-         PHSP;
+0.000002  eta eta pi+ K-         PHSP;
+0.000002  pi0 eta' pi+ K-        PHSP;
+0.000001  eta eta' pi+ K-        PHSP;
+0.000010  rho0 anti-K_0*0        PHSP;
+0.000020  rho+ K_0*-             PHSP;
+0.000010  pi- pi+ anti-K_0*0     PHSP;
+0.000010  pi0 pi0 anti-K_0*0     PHSP;
+0.000010  pi+ pi0 K_0*-          PHSP;
+0.000002  anti-K*0 f_0           SVS;
+0.000010  a_10 anti-K0           SVS;
+0.000010  a_1+ K-                SVS;
+0.0000066 b_10 anti-K0           SVS;
+0.000010  b_1+ K-                SVS;
+0.000005  anti-K*0 a_00          SVS;
+0.000005  K*- a_0+               SVS;
+#--- 4-body K*-K*, rho-K-K, K*-K-pi, phi-pi-pi, pi-pi-K-K ------------
+0.00000105  K*0   anti-K*0         SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  K*+   K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi   rho0             SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi   pi- pi+          PHSP;
+0.000001  K+ K- rho0             PHSP;
+0.000001  K0 anti-K0 rho0        PHSP;
+0.000001  K+ K- omega            PHSP;
+0.000001  K0 anti-K0 omega       PHSP;
+0.000001  K+ anti-K0 rho-        PHSP;
+0.000001  K- K0 rho+             PHSP;
+0.000001  K+ K- pi- pi+          PHSP;
+0.000001  K+ K- pi0 pi0          PHSP;
+0.0000005 K+ K- pi0 eta          PHSP;
+0.0000002 K+ K- eta eta          PHSP;
+0.0000002 K+ K- pi0 eta'         PHSP;
+0.0000001 K+ K- eta eta'         PHSP;
+0.000001  K0 anti-K0 pi- pi+     PHSP;
+0.000001  K0 anti-K0 pi0 pi0     PHSP;
+0.0000005 K0 anti-K0 pi0 eta     PHSP;
+0.0000002 K0 anti-K0 eta eta     PHSP;
+0.0000002 K0 anti-K0 pi0 eta'    PHSP;
+0.0000001 K0 anti-K0 eta eta'    PHSP;
+0.000001  K*- K+ pi0             PHSP;
+0.000001  K*+ K- pi0             PHSP;
+0.000001  K0 anti-K*0 pi0        PHSP;
+0.000001  K*0 anti-K0 pi0        PHSP;
+0.000001  K- K*0 pi+             PHSP;
+0.000001  K*- K0 pi+             PHSP;
+0.000001  K+ anti-K*0 pi-        PHSP;
+0.000001  K*+ anti-K0 pi-        PHSP;
+
+#--- 4-body phi-K*, phi-K-pi, K-K-K*, pi-K-K-K -----------------------
+0.0000095 phi   anti-K*0         SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  phi K- pi+             PHSP;
+0.000010  phi anti-K0 pi0        PHSP;
+0.000007  phi anti-K0 eta        PHSP;
+0.000010  K- K+ anti-K*0         PHSP; 
+0.000010  K- K*+ anti-K0         PHSP;
+0.000010  K*- K+ anti-K0         PHSP;
+0.000010  K+ K- K- pi+           PHSP;
+0.000010  K+ K- anti-K0 pi0      PHSP;
+0.000005  K+ K- anti-K0 eta      PHSP;
+0.000010  K0 anti-K0 anti-K0 pi0 PHSP;
+0.000005  K0 anti-K0 anti-K0 eta PHSP;
+0.0000046 phi   anti-K_0*0       PHSP;
+0.000010  K- K+ anti-K_0*0       PHSP;
+0.000010  K- K_0*+ anti-K0       PHSP;
+0.000010  K_0*- K+ anti-K0       PHSP;
+0.000002  phi   anti-K'_10       PHSP;
+0.000002  K- K+ anti-K'_10       PHSP;
+#--- 4-body phi-phi, phi-K-K, K-K-K-K --------------------------------
+0.0000001 phi phi                SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi K+ K-              PHSP;
+0.0000001 phi K0 anti-K0         PHSP;
+0.0000001 K+ K- K+ K-            PHSP;
+0.0000001 K0 anti-K0 K+ K-       PHSP;
+0.0000001 K0 anti-K0 K0 anti-K0  PHSP;
+#5-body modes from Jim Smith - 1/30/05
+#--- 5-body rho-rho-pi, rho-pi-pi-pi, pi-pi-pi-pi-pi -------------------
+0.000010  rho0 rho0 pi0          PHSP;
+0.000005  rho0 rho0 eta          PHSP;
+0.000002  rho0 rho0 eta'         PHSP;
+0.000010  omega rho0 pi0         PHSP;
+0.000005  omega rho0 eta         PHSP;
+0.000002  omega rho0 eta'        PHSP;
+0.000010  rho+ rho- pi0          PHSP;
+0.000005  rho+ rho- eta          PHSP;
+0.000002  rho+ rho- eta'         PHSP;
+0.000010  rho0 rho-  pi+          PHSP;
+0.0000066 rho0 rho+ pi-          PHSP; 
+0.000010  omega rho- pi+         PHSP;
+0.0000066 omega rho+ pi-         PHSP;
+0.000002  rho0 pi+ pi- pi0       PHSP;
+0.000001  rho0 pi+ pi- eta       PHSP;
+0.000001  rho0 pi+ pi- eta'      PHSP;
+0.000002  rho0 pi0 pi0 pi0       PHSP;
+0.000001  rho0 pi0 eta pi0       PHSP;
+0.000001  rho0 pi0 eta' pi0      PHSP;
+0.000002  omega pi+ pi- pi0      PHSP;
+0.000001  omega pi+ pi- eta      PHSP;
+0.000001  omega pi+ pi- eta'     PHSP;
+0.000002  omega pi0 pi0 pi0      PHSP;
+0.000001  omega pi0 eta pi0      PHSP;
+0.000001  omega pi0 eta' pi0     PHSP;
+
+0.000002  rho- pi+ pi+ pi-       PHSP;
+0.000002  rho- pi+ pi0 pi0       PHSP;
+0.000001  rho- pi+ pi0 eta       PHSP;
+0.000001  rho- pi+ pi0 eta'      PHSP;
+0.000002  rho+ pi- pi+ pi-       PHSP;
+0.000002  rho+ pi- pi0 pi0       PHSP;
+0.000001  rho+ pi- pi0 eta       PHSP;
+0.000001  rho+ pi- pi0 eta'      PHSP;
+0.000001  pi+ pi- pi+ pi- pi0    PHSP;
+0.0000005 pi+ pi- pi+ pi- eta    PHSP;
+0.0000002 pi+ pi- pi+ pi- eta'   PHSP;
+0.000001  pi+ pi- pi0 pi0 pi0    PHSP;
+0.0000005 pi+ pi- pi0 pi0 eta    PHSP;
+0.0000002 pi+ pi- pi0 pi0 eta'   PHSP;
+0.0000005 pi0 pi0 pi0 pi0 pi0    PHSP;
+0.0000002 pi0 pi0 pi0 pi0 eta    PHSP;
+0.0000001 pi0 pi0 pi0 pi0 eta'   PHSP;
+
+0.000010  a_1-  rho+             PHSP;
+0.000010  a_1+  rho-             PHSP;
+0.000010  a_10  rho0             PHSP;
+0.000010  a_10  omega            PHSP;
+0.000010  a_1-  pi+ pi0          PHSP;
+0.000005  a_1-  pi+ eta          PHSP;
+0.000002  a_1-  pi+ eta'         PHSP;
+0.000010  a_10  pi+ pi-          PHSP;
+0.000010  a_10  pi0 pi0          PHSP;
+0.000005  a_10  pi0 eta          PHSP;
+0.000002  a_10  pi0 eta'         PHSP;
+0.000010  a_1+  pi- pi0          PHSP;
+0.000005  a_1+  pi- eta          PHSP;
+0.000002  a_1+  pi- eta'         PHSP;
+
+0.000010  rho- f_0 pi+           PHSP;
+0.000010  rho+ f_0 pi-           PHSP;
+0.000010  rho0 f_0 pi0           PHSP;
+0.000005  rho0 f_0 eta           PHSP;
+0.000002  rho0 f_0 eta'          PHSP;
+0.000010  omega f_0 pi0          PHSP;
+
+0.000010  rho- a_00 pi+          PHSP;
+0.000010  rho- a_0+ pi0          PHSP;
+0.000005  rho- a_0+ eta          PHSP;
+0.000002  rho- a_0+ eta'         PHSP;
+
+0.000010  rho0 a_0- pi+          PHSP;
+0.000010  rho0 a_0+ pi-          PHSP;
+0.000010  rho0 a_00 pi0          PHSP;
+0.000005  rho0 a_00 eta          PHSP;
+0.000002  rho0 a_00 eta'         PHSP;
+0.000010  rho+ a_00 pi-          PHSP;
+0.000010  rho+ a_0- pi0          PHSP;
+0.000005  rho+ a_0- eta          PHSP;
+0.000002  rho+ a_0- eta'         PHSP;
+0.000002  f_0 pi+ pi- pi0        PHSP;
+0.000001  f_0 pi+ pi- eta        PHSP;
+0.000001  f_0 pi+ pi- eta'       PHSP;
+0.000002  f_0 pi0 pi0 pi0        PHSP;
+0.000001  f_0 pi0 pi0 eta        PHSP;
+0.000001  f_0 pi0 pi0 eta'       PHSP;
+0.000001  a_00 pi+ pi- pi0       PHSP;
+0.000001  a_00 pi0 pi0 pi0       PHSP;
+
+0.000001  a_0- pi+ pi+ pi-       PHSP;
+0.000001  a_0- pi+ pi0 pi0       PHSP;
+0.000001  a_0+ pi- pi+ pi-       PHSP;
+0.000001  a_0+ pi- pi0 pi0       PHSP;
+#--- 5-body rho-K*pi, rho-2pi-K, K*-3pi, 4pi-K -------------------
+0.000010  rho0 K*- pi+           PHSP;
+0.000010  rho+ K*- pi0           PHSP;
+0.000005  rho+ K*- eta           PHSP;
+0.000002  rho+ K*- eta'          PHSP;
+0.000010  rho0 anti-K*0 pi0      PHSP;
+0.000005  rho0 anti-K*0 eta      PHSP;
+0.000002  rho0 anti-K*0 eta'     PHSP;
+
+0.000010  rho+ anti-K*0 pi-      PHSP;
+0.000010  rho- anti-K*0 pi+      PHSP;
+0.000010  omega K*- pi+          PHSP;
+0.000010  omega anti-K*0 pi0     PHSP;
+0.000005  omega anti-K*0 eta     PHSP;
+0.000002  omega anti-K*0 eta'    PHSP;
+0.000010  pi+ pi- K*- pi+        PHSP;
+0.000010  pi+ pi0 K*- pi0        PHSP;
+0.000005  pi+ pi0 K*- eta        PHSP;
+0.000002  pi+ pi0 K*- eta'       PHSP;
+0.000010  pi+ pi- anti-K*0 pi0   PHSP;
+0.000005  pi+ pi- anti-K*0 eta   PHSP;
+0.000002  pi+ pi- anti-K*0 eta'  PHSP;
+0.000010  pi0 pi0 anti-K*0 pi0   PHSP;
+0.000005  pi0 pi0 anti-K*0 eta   PHSP;
+0.000002  pi0 pi0 anti-K*0 eta'  PHSP;
+0.000010  rho-  K- pi+ pi+       PHSP;
+0.000010  rho0  K- pi+ pi0       PHSP;
+0.000005  rho0  K- pi+ eta       PHSP;
+0.000002  rho0  K- pi+ eta'      PHSP;
+0.000010  rho+  K- pi+ pi-       PHSP;
+0.000010  rho+  K- pi0 pi0       PHSP;
+0.000005  rho+  K- pi0 eta       PHSP;
+0.000002  rho+  K- pi0 eta'      PHSP;
+
+0.000010  rho0  anti-K0 pi+ pi-  PHSP;
+0.000010  rho0  anti-K0 pi0 pi0  PHSP;
+0.000005  rho0  anti-K0 pi0 eta  PHSP;
+0.000002  rho0  anti-K0 pi0 eta' PHSP;
+0.000010  rho-  anti-K0 pi+ pi0  PHSP;
+0.000005  rho-  anti-K0 pi+ eta  PHSP;
+0.000002  rho-  anti-K0 pi+ eta' PHSP;
+0.000010  rho+  anti-K0 pi- pi0  PHSP;
+0.000005  rho+  anti-K0 pi- eta  PHSP;
+0.000002  rho+  anti-K0 pi- eta' PHSP;
+0.000010  omega K- pi+ pi0       PHSP;
+0.000005  omega K- pi+ eta       PHSP;
+0.000002  omega K- pi+ eta'      PHSP;
+0.000010  omega anti-K0 pi+ pi-  PHSP;
+0.000010  omega anti-K0 pi0 pi0  PHSP;
+0.000005  omega anti-K0 pi0 eta  PHSP;
+0.000002  omega anti-K0 pi0 eta' PHSP;
+0.000002  pi+ pi- pi+ K- pi0     PHSP;
+0.000001  pi+ pi- pi+ K- eta     PHSP;
+0.000001  pi+ pi- pi+ K- eta'    PHSP;
+0.000002  pi0 pi0 pi+ K- pi0     PHSP;
+0.000001  pi0 pi0 pi+ K- eta     PHSP;
+0.000001  pi0 pi0 pi+ K- eta'    PHSP;
+0.000002  pi+ pi- pi+ anti-K0 pi- PHSP;
+0.000002  pi+ pi- pi0 anti-K0 pi0 PHSP;
+0.000001  pi+ pi- pi0 anti-K0 eta PHSP;
+0.000001  pi+ pi- pi0 anti-K0 eta' PHSP;
+0.000002  pi0 pi0 pi0 anti-K0 pi0 PHSP;
+0.000001  pi0 pi0 pi0 anti-K0 eta PHSP;
+0.000001  pi0 pi0 pi0 anti-K0 eta' PHSP;
+
+0.000010  rho- anti-K_0*0 pi+    PHSP;
+0.000010  rho0 anti-K_0*0 pi0    PHSP;
+0.000005  rho0 anti-K_0*0 eta    PHSP;
+0.000002  rho0 anti-K_0*0 eta'   PHSP;
+0.000010  rho+ anti-K_0*0 pi-    PHSP;
+0.000010  rho+ K_0*- pi0         PHSP;
+0.000005  rho+ K_0*- eta         PHSP;
+0.000002  rho+ K_0*- eta'        PHSP;
+0.000010  rho0 K_0*- pi+         PHSP;
+0.000010  pi+ pi- K_0*- pi+      PHSP;
+0.000010  pi+ pi0 K_0*- pi0      PHSP;
+0.000005  pi+ pi0 K_0*- eta      PHSP;
+0.000002  pi+ pi0 K_0*- eta'     PHSP;
+0.000010  pi+ pi- anti-K_0*0 pi0 PHSP;
+0.000005  pi+ pi- anti-K_0*0 eta PHSP;
+0.000002  pi+ pi- anti-K_0*0 eta' PHSP;
+0.000010  pi0 pi0 anti-K_0*0 pi0 PHSP;
+0.000005  pi0 pi0 anti-K_0*0 eta PHSP;
+0.000005  pi0 pi0 anti-K_0*0 eta' PHSP;
+0.000010  K*- f_0 pi+            PHSP;
+0.000010  anti-K*0 f_0 pi0       PHSP;
+0.000005  anti-K*0 f_0 eta       PHSP;
+0.000002  anti-K*0 f_0 eta'      PHSP;
+0.000010  a_1- anti-K0 pi+       PHSP;
+0.000020  a_10 anti-K*0          PHSP;
+0.000010  a_10 K- pi+            PHSP;
+0.000010  a_10 anti-K0 pi0       PHSP;
+0.000020  a_1+ K*-               PHSP;
+0.000010  a_1+ K- pi0            PHSP;
+0.000010  a_1+ anti-K0 pi-       PHSP;
+0.000005  K*- a_00 pi+           PHSP;
+0.000005  K*- a_0+ pi0           PHSP;
+0.000002  K*- a_0+ eta           PHSP;
+0.000001  K*- a_0+ eta'          PHSP;
+0.000005  anti-K*0 a_0+ pi-      PHSP;
+0.000005  anti-K*0 a_0- pi+      PHSP;
+0.000005  anti-K*0 a_00 pi0      PHSP;
+0.000002  anti-K*0 a_00 eta      PHSP;
+0.000001  anti-K*0 a_00 eta'     PHSP;
+#--- 5-body K*-K*pi,rho-K-Kpi, K*-K-2pi, phi-3pi, 3pi-K-K ------------
+0.000001  K*+ anti-K*0 pi-       PHSP;
+0.000001  K*- K*0 pi+            PHSP;
+
+0.000001  K*- K*+ pi0            PHSP;
+0.0000005 K*- K*+ eta            PHSP;
+0.0000002 K*- K*+ eta'           PHSP;
+0.000001  K*0 anti-K*0 pi0       PHSP;
+0.0000005 K*0 anti-K*0 eta       PHSP;
+0.0000002 K*0 anti-K*0 eta'      PHSP;
+
+0.000001  phi rho- pi+           PHSP;
+0.000001  phi rho+ pi-           PHSP;
+0.000001  phi rho0 pi0           PHSP;
+0.0000005 phi rho0 eta           PHSP;
+0.0000002 phi rho0 eta'          PHSP;
+0.000001  phi omega pi0          PHSP;
+0.000001  phi f_0 pi0            PHSP;
+0.000001  phi pi+ pi- pi0        PHSP;
+0.0000005 phi pi+ pi- eta        PHSP;
+0.0000002 phi pi+ pi- eta'       PHSP;
+0.000001  phi pi0 pi0 pi0        PHSP;
+0.0000005 phi pi0 pi0 eta        PHSP;
+0.0000002 phi pi0 pi0 eta'       PHSP;
+0.000001  K0 K- rho0 pi+         PHSP;
+0.000001  K0 K- rho+ pi0         PHSP;
+0.0000005 K0 K- rho+ eta         PHSP;
+0.0000002 K0 K- rho+ eta'        PHSP;
+0.000001  K0 K- omega pi+        PHSP;
+0.000001  K0 K- f_0 pi+          PHSP;
+0.000001  K+ K- rho+ pi-         PHSP;
+0.000001  K+ K- rho- pi+         PHSP;
+0.000001  K+ K- rho0 pi0         PHSP;
+0.0000005 K+ K- rho0 eta         PHSP;
+0.0000002 K+ K- rho0 eta'        PHSP;
+0.000001  K+ K- omega pi0        PHSP;
+0.0000005 K+ K- omega eta        PHSP;
+0.0000002 K+ K- omega eta'       PHSP;
+0.000001  K+ K- f_0 pi0          PHSP;
+0.000001  K0 anti-K0 rho+ pi-    PHSP;
+0.000001  K0 anti-K0 rho- pi+    PHSP;
+0.000001  K0 anti-K0 rho0 pi0    PHSP;
+0.0000005 K0 anti-K0 rho0 eta    PHSP;
+0.0000002 K0 anti-K0 rho0 eta'   PHSP;
+0.000001  K0 anti-K0 omega pi0   PHSP;
+0.000001  K0 anti-K0 f_0 pi0     PHSP;
+0.000001  anti-K0 K+ rho0 pi-    PHSP;
+0.000001  anti-K0 K+ rho- pi0    PHSP;
+0.0000005 anti-K0 K+ rho- eta    PHSP;
+0.0000002 anti-K0 K+ rho- eta'   PHSP;
+0.000001  anti-K0 K+ omega pi-   PHSP;
+0.000001  anti-K0 K+ f_0 pi-     PHSP;
+0.0000002 anti-K0 K+ pi- pi+ pi- PHSP;
+0.0000002 anti-K0 K+ pi- pi0 pi0 PHSP;
+0.0000001 anti-K0 K+ pi- pi0 eta PHSP;
+0.0000001 anti-K0 K+ pi- pi0 eta' PHSP;
+
+0.0000002 K+ K- pi+ pi- pi0      PHSP;
+0.0000001 K+ K- pi+ pi- eta      PHSP;
+0.0000001 K+ K- pi+ pi- eta'     PHSP;
+0.0000002 K+ K- pi0 pi0 pi0      PHSP;
+0.0000001 K+ K- pi0 pi0 eta      PHSP;
+0.0000001 K+ K- pi0 pi0 eta'     PHSP;
+0.0000002 K0 anti-K0 pi+ pi- pi0 PHSP;
+0.0000001 K0 anti-K0 pi+ pi- eta PHSP;
+0.0000001 K0 anti-K0 pi+ pi- eta' PHSP;
+0.0000002 K0 anti-K0 pi0 pi0 pi0 PHSP;
+0.0000001 K0 anti-K0 pi0 pi0 eta PHSP;
+0.0000001 K0 anti-K0 pi0 pi0 eta' PHSP;
+0.0000002 K0 K- pi+ pi+ pi-      PHSP;
+0.0000002 K0 K- pi+ pi0 pi0      PHSP;
+0.0000001 K0 K- pi+ pi0 eta      PHSP;
+0.0000001 K0 K- pi+ pi0 eta'     PHSP;
+
+0.000001  K+ anti-K*0 pi- pi0    PHSP;
+0.0000005 K+ anti-K*0 pi- eta    PHSP;
+0.0000002 K+ anti-K*0 pi- eta'   PHSP;
+0.000001  K*+ anti-K0 pi- pi0    PHSP;
+0.0000005 K*+ anti-K0 pi- eta    PHSP;
+0.0000002 K*+ anti-K0 pi- eta'   PHSP;
+0.000001  K*- K+ pi+ pi-         PHSP;
+0.000001  K*- K+ pi0 pi0         PHSP;
+0.0000005 K*- K+ pi0 eta         PHSP;
+0.0000002 K*- K+ pi0 eta'        PHSP;
+0.000001  K*+ K- pi+ pi-         PHSP;
+0.000001  K*+ K- pi0 pi0         PHSP;
+0.0000005 K*+ K- pi0 eta         PHSP;
+0.0000002 K*+ K- pi0 eta'        PHSP;
+
+0.000001  K0 anti-K*0 pi+ pi-    PHSP;
+0.000001  K0 anti-K*0 pi0 pi0    PHSP;
+0.0000005 K0 anti-K*0 pi0 eta    PHSP;
+0.0000002 K0 anti-K*0 pi0 eta'   PHSP;
+0.000001  K*0 anti-K0 pi+ pi-    PHSP;
+0.000001  K*0 anti-K0 pi0 pi0    PHSP;
+0.0000005 K*0 anti-K0 pi0 eta    PHSP;
+0.0000002 K*0 anti-K0 pi0 eta'   PHSP;
+0.000001  K*0 K- pi+ pi0         PHSP;
+0.0000005 K*0 K- pi+ eta         PHSP;
+0.0000002 K*0 K- pi+ eta'        PHSP;
+0.000001  K*- K0 pi+ pi0         PHSP;
+0.0000005 K*- K0 pi+ eta         PHSP;
+0.0000002 K*- K0 pi+ eta'        PHSP;
+
+#--- 5-body phi-K*-pi, phi-K-2pi, K-K-K*-pi, 2pi-K-K-K --------------------
+0.000010  phi K*- pi+            PHSP;
+0.000010  phi anti-K*0 pi0       PHSP;
+0.000005  phi anti-K*0 eta       PHSP;
+0.000002  phi anti-K*0 eta'      PHSP;
+0.000010  phi K- rho+            PHSP;
+0.000010  phi anti-K0 rho0       PHSP;
+0.000010  phi anti-K0 omega      PHSP;
+0.000010  phi anti-K0 f_0        PHSP;
+0.000010  phi anti-K0 pi+ pi-    PHSP;
+0.000010  phi anti-K0 pi0 pi0    PHSP;
+0.000005  phi anti-K0 pi0 eta    PHSP;
+0.000002  phi anti-K0 pi0 eta'   PHSP;
+0.000010  phi K- pi+ pi0         PHSP;
+0.000005  phi K- pi+ eta         PHSP;
+0.000002  phi K- pi+ eta'        PHSP;
+0.000010  K+ K- K*- pi+          PHSP;
+0.000010  K- K*+ K- pi+          PHSP;
+
+0.000010  K- K0 anti-K*0 pi+     PHSP;
+0.000010  K- K*0 anti-K0 pi+     PHSP;
+0.000010  K0 K*- anti-K0 pi+     PHSP;
+0.000010  K- K+ K- rho+          PHSP;
+0.000010  K- K+ anti-K0 rho0     PHSP;
+0.000010  K- K+ anti-K0 omega    PHSP;
+0.000010  K- K+ anti-K0 f_0      PHSP;
+0.000010  K+ K- K- pi+ pi0       PHSP;
+0.000005  K+ K- K- pi+ eta       PHSP;
+0.000002  K+ K- K- pi+ eta'      PHSP;
+0.000010  K+ K- anti-K0 pi+ pi-  PHSP;
+0.000010  K+ K- anti-K0 pi0 pi0  PHSP;
+0.000005  K+ K- anti-K0 pi0 eta  PHSP;
+0.000005  K+ K- anti-K0 pi0 eta' PHSP;
+0.000010  K0 anti-K0 K- pi+ pi0  PHSP;
+0.000005  K0 anti-K0 K- pi+ eta  PHSP;
+0.000002  K0 anti-K0 K- pi+ eta' PHSP;
+0.000010  K0 anti-K0 anti-K0 pi+ pi-  PHSP;
+0.000010  K0 anti-K0 anti-K0 pi0 pi0  PHSP;
+0.000005  K0 anti-K0 anti-K0 pi0 eta  PHSP;
+0.000002  K0 anti-K0 anti-K0 pi0 eta' PHSP;
+0.000010  phi   K_0*- pi+        PHSP;
+0.000010  phi   anti-K_0*0 pi0   PHSP;
+0.000010  K+ K- K_0*- pi+        PHSP;
+0.000010  K+ K- anti-K_0*0 pi0   PHSP;
+0.000005  K+ K- anti-K_0*0 eta   PHSP;
+0.000002  K+ K- anti-K_0*0 eta'  PHSP;
+0.000010  K- K- K_0*+ pi+        PHSP;
+0.000010  K- anti-K_0*0 K0 pi+   PHSP;
+0.000010  K- K_0*0 anti-K0 pi+   PHSP;
+0.000010  K0 anti-K0 K_0*- pi+   PHSP;
+0.000010  K0 anti-K0 anti-K_0*0 pi0  PHSP;
+0.000005  K0 anti-K0 anti-K_0*0 eta  PHSP;
+0.000002  K0 anti-K0 anti-K_0*0 eta' PHSP;
+0.000010  anti-K0 anti-K0 K_0*0 pi0  PHSP;
+0.000005  anti-K0 anti-K0 K_0*0 eta  PHSP;
+0.000002  anti-K0 anti-K0 K_0*0 eta' PHSP;
+0.000002  phi   K'_1- pi+        PHSP;
+0.000002  phi   anti-K'_10 pi0   PHSP;
+0.000001  phi   anti-K'_10 eta   PHSP;
+0.000001  phi   anti-K'_10 eta'  PHSP;
+0.000002  K- K+ K'_1- pi+        PHSP;
+0.000002  K- K+ anti-K'_10 pi0   PHSP;
+0.000001  K- K+ anti-K'_10 eta   PHSP;
+0.000001  K- K+ anti-K'_10 eta'  PHSP;
+0.000002  K- K- K'_1+ pi+        PHSP;
+0.000002  K0 anti-K0 anti-K'_10 pi0  PHSP;
+0.000001  K0 anti-K0 anti-K'_10 eta  PHSP;
+0.000001  K0 anti-K0 anti-K'_10 eta' PHSP;
+0.000002  anti-K0 anti-K0 K'_10 pi0  PHSP;
+0.000001  anti-K0 anti-K0 K'_10 eta  PHSP;
+0.000001  anti-K0 anti-K0 K'_10 eta' PHSP;
+#--- 5-body phi-phi-K, phi-phi-pi, phi-3K, phi-K-K-pi, 5K, 4K-pi ------------
+0.0000041 phi phi anti-K0        PHSP;
+0.0000001 phi phi pi0            PHSP;
+0.000001  phi K+ K- anti-K0      PHSP;
+0.000001  phi K0 anti-K0 anti-K0 PHSP;
+0.0000001 phi K+ K- pi0          PHSP;
+0.0000001 phi K0 anti-K0 pi0     PHSP;
+0.0000001 phi anti-K0 K+ pi-     PHSP;
+0.0000001 phi K0 K- pi+          PHSP;
+0.0000001 K+ K- K+ K- anti-K0    PHSP;
+0.0000001 K0 anti-K0 K+ K- anti-K0 PHSP;
+0.0000001 K+ K- K+ K- pi0        PHSP;
+0.0000001 K+ K- K0 anti-K0 pi0   PHSP;
+0.0000001 K+ K- K+ anti-K0 pi-   PHSP;
+0.0000001 K+ K- K0 K- pi+        PHSP;
+0.0000001 K0 anti-K0 K0 anti-K0 pi0 PHSP;
+0.0000001 K0 anti-K0 anti-K0 K+ pi- PHSP;
+0.0000001 K0 anti-K0 K0 K- pi+   PHSP;
+#--- 6-body phi-phi-K*, phi-3Kpi, phi-K-K-2pi ------------
+0.000003  phi phi anti-K*0       PHSP;
+0.000001  phi K+ K- anti-K*0     PHSP;
+0.000001  phi K0 anti-K0 anti-K*0 PHSP;
+0.000001  phi anti-K0 K- K*+     PHSP;
+0.000001  phi anti-K0 anti-K0 K*0 PHSP;
+0.0000001 phi K+ K- K- pi+       PHSP;
+0.0000001 phi K+ K- anti-K0 pi0  PHSP;
+0.0000001 phi K0 anti-K0 K- pi+  PHSP;
+0.0000001 phi K0 anti-K0 anti-K0 pi0 PHSP;
+#--- 6-body a1a1 ------------
+0.000050  a_10  a_10             PHSP;
+0.000050  a_1+  a_1-             PHSP;
+#
+#       B -> cc= s          
+#
+#   Charmonium states - updated from Lange's recommendations (august 23,2000)
+# Based on new BABAR results I'm making the following changes (Lange, March 13, 2001
+#       J/psi K0   was 10, now 9 x 10^-4
+#       J/psi pi0   was 2.5, now 2.0 x 10^-4
+#       J/psi Kstar was 15, now 13
+#       Fix chic1 KS CP eigenstate
+#       adding J/psi rho and omega - for lack of better thing will use the Kstar helicity amplitudes.
+#       Psi2sKs  30 ->31
+#       chic1 Kstar: 12 ->6
+
+
+0.000436    J/psi  K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000436    J/psi  K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.0002215   J/psi  K*S               SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.0002215   J/psi  K*L               SVV_CP beta  dm -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000886     J/psi  anti-K*0T          SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.0000205    J/psi  pi0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000016     J/psi  rho0              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000030     J/psi  omega             SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+
+0.0012     J/psi  K-  pi+           PHSP;
+0.0001     J/psi  anti-K0   pi0           PHSP;
+#rl0.0007     J/psi  anti-K0   pi+  pi-      PHSP;
+#rl0.00035     J/psi  anti-K0   pi0  pi0      PHSP;
+#rl0.00035     J/psi  K-  pi+  pi0      PHSP;
+0.0013     J/psi  anti-K_10              SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0001     J/psi  anti-K'_10             SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0005     J/psi  anti-K_2*0              PHSP;
+0.000094    J/psi  phi anti-K0        PHSP;
+#
+0.00031     psi(2S)  K_S0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00031     psi(2S)  K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00012    psi(2S)  K*S              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00012    psi(2S)  K*L              SVV_CP beta  dm  -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00048     psi(2S)  anti-K*0T             SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+
+0.0004     psi(2S)  K-  pi+          PHSP;
+0.0002     psi(2S)  anti-K0   pi0          PHSP;
+0.0002     psi(2S)  anti-K0   pi+  pi-     PHSP;
+0.0001     psi(2S)  anti-K0   pi0  pi0     PHSP;
+0.0001     psi(2S)  K-  pi+  pi0     PHSP;
+0.0004     psi(2S)  anti-K_10              PHSP;
+
+#
+0.0005    eta_c K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0005    eta_c K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00027   K*S  eta_c              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00027   K*L  eta_c              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00108    anti-K*0T eta_c              SVS;
+0.0002    eta_c  K-  pi+          PHSP;
+0.0001    eta_c  anti-K0   pi0          PHSP;
+0.0002    eta_c  anti-K0   pi+  pi-     PHSP;
+0.0001    eta_c  anti-K0   pi0  pi0     PHSP;
+0.0001    eta_c  K-  pi+  pi0     PHSP;
+#
+0.00024    eta_c(2S) K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00024    eta_c(2S) K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00011   K*S  eta_c(2S)               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00011   K*L  eta_c(2S)               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00044   anti-K*0T eta_c(2S)          SVS;
+0.00008    eta_c(2S)   K-  pi+          PHSP;
+0.00005    eta_c(2S)   anti-K0   pi0          PHSP;
+0.00008    eta_c(2S)   anti-K0   pi+  pi-     PHSP;
+0.00005    eta_c(2S)   anti-K0   pi0  pi0     PHSP;
+0.00005    eta_c(2S)   K-  pi+  pi0     PHSP;
+#
+0.00005     chi_c0 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00005     chi_c0 K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00005    K*S  chi_c0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00005    K*L  chi_c0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0002     anti-K*0T chi_c0              SVS;
+0.0002     chi_c0  K-  pi+          PHSP;
+0.0001     chi_c0  anti-K0   pi0          PHSP;
+0.0002     chi_c0  anti-K0   pi+  pi-     PHSP;
+0.0001     chi_c0  anti-K0   pi0  pi0     PHSP;
+0.0001     chi_c0  K-  pi+  pi0     PHSP;
+#
+0.00020     chi_c1 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00020     chi_c1 K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000053    K*S  chi_c1              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000053    K*L  chi_c1              SVV_CP beta  dm -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000213   chi_c1  anti-K*0T             SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.0004     chi_c1  K-  pi+          PHSP;
+0.0002     chi_c1  anti-K0   pi0          PHSP;
+0.0004     chi_c1  anti-K0   pi+  pi-     PHSP;
+0.0002     chi_c1  anti-K0   pi0  pi0     PHSP;
+0.0002     chi_c1  K-  pi+  pi0     PHSP;
+#
+0.00001     chi_c2 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00001     chi_c2 K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000005    K*S  chi_c2              PHSP;
+0.000005    K*L  chi_c2              PHSP;
+0.00002     chi_c2  anti-K*0T             PHSP;
+0.0002     chi_c2  K-  pi+          PHSP;
+0.0001     chi_c2  anti-K0   pi0          PHSP;
+0.0002     chi_c2  anti-K0   pi+  pi-     PHSP;
+0.0001    chi_c2  anti-K0   pi0  pi0     PHSP;
+0.0001     chi_c2  K-  pi+  pi0     PHSP;
+#
+0.00024     psi(3770)  K_S0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00024     psi(3770)  K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00008     psi(3770)  K*S              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00008     psi(3770)  K*L              SVV_CP beta  dm  -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00032     psi(3770)  anti-K*0T             SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.00014     psi(3770)  K-  pi+                PHSP;
+0.00014     psi(3770)  anti-K0   pi0          PHSP;
+0.00014     psi(3770)  anti-K0   pi+  pi-     PHSP;
+0.00007     psi(3770)  anti-K0   pi0  pi0     PHSP;
+0.00007     psi(3770)  K-  pi+  pi0           PHSP;
+0.00029     psi(3770)  anti-K_10              PHSP;
+#
+
+
+
+# b-> c (dc=)
+#
+0.00023     D-   D+              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+#  See Ref [B1] for the next 3:
+0.000610    D*+  D-              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000610    D*-  D+              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00082     D*+ D*-              SVV_CP beta dm  1  0.56 0.0 0.96 0.0 0.47 0.0;
+
+#
+# b -> c (sc=)  -->  D Ds X   Sum = 10%
+#
+0.0065   D+  D_s-                        PHSP;
+0.0080   D*+  D_s-                       SVS;
+0.0074   D_s*- D+                        SVS;
+0.0177   D*+  D_s*-              SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+0.0006   D'_1+  D_s-             SVS;
+0.0012   D'_1+  D_s*-            SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0012   D_1+  D_s-              SVS;
+0.0024   D_1+  D_s*-             SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0042   D_2*+  D_s-                      STS;
+0.0040   D_2*+  D_s*-                    PHSP;
+#
+0.0018   D_s-  D+  pi0                   PHSP;
+0.0037   D_s-  D0  pi+                   PHSP;
+0.0018   D_s*- D+  pi0                   PHSP;
+0.0037   D_s*- D0  pi+                   PHSP;
+0.0030   D_s-  D+  pi+  pi-              PHSP;
+0.0022   D_s-  D+  pi0  pi0              PHSP;
+0.0022   D_s-  D0  pi+  pi0              PHSP;
+0.0030   D_s*- D+  pi+  pi-              PHSP;
+0.0022   D_s*- D+  pi0  pi0              PHSP;
+0.0022   D_s*- D0  pi+  pi0              PHSP;
+#
+# b -> c (sc=)  -->  D D= K  X   Sum = 8%
+# Update: Ref. [B1]
+#
+# External W-emission amplitude
+0.0017   D+  anti-D0  K-                         PHSP;
+0.0046   D+  anti-D*0  K-                        PHSP; 
+0.0031   D*+ anti-D0  K-                         PHSP;
+0.0118   D*+ anti-D*0  K-                        PHSP;
+# External+internal W-emission amplitude
+0.0015   D+  D-  anti-K0                         PHSP; 
+0.0018   D*+ D-  anti-K0                         PHSP; 
+0.0047   D+  D*-  anti-K0                        PHSP; 
+0.0088   D*+ D*-  anti-K0                        PHSP; 
+# Internal W-emission amplitude (color suppressed modes)
+0.0005   D0  anti-D0  anti-K0                    PHSP;
+
+0.0005   D0  anti-D*0  anti-K0                   PHSP; 
+0.0005   D*0  anti-D0  anti-K0                   PHSP; 
+0.0005   D*0  anti-D*0  anti-K0                  PHSP; 
+
+0.0025   D+  anti-D0  K*-                        PHSP;
+0.0025   D*+ anti-D0  K*-                        PHSP;
+0.0025   D+  anti-D*0  K*-                       PHSP;
+0.0050   D*+ anti-D*0  K*-                       PHSP;
+#
+0.0025   D+  D-  anti-K*0                        PHSP;
+0.0025   D*+ D-  anti-K*0                        PHSP;
+0.0025   D+  D*-  anti-K*0                       PHSP;
+0.0050   D*+ D*-  anti-K*0                       PHSP;
+#
+0.0005   anti-D0  D0  anti-K*0                   PHSP;
+0.0005   anti-D0  D*0  anti-K*0                  PHSP;
+0.0005   anti-D*0  D0  anti-K*0                  PHSP;
+0.0010   anti-D*0  D*0  anti-K*0                 PHSP;
+#
+#
+#  B -> D(*) X  Exclusive modes
+#
+0.00276  D*+ pi-                        SVS;
+0.00268  D+  pi-                        PHSP;
+0.0075   rho- D+                        SVS;
+#
+# D* rho HELAMP parameters taken from ICHEP 98-852.
+0.0068   rho- D*+                 SVV_HELAMP 0.152 1.47 0.936 0 0.317 0.19; 
+0.0005   D+  pi-  pi0                   PHSP;
+0.0082  D*+ pi-  pi0                   PHSP;
+0.0005   D0  pi+  pi-                   PHSP;
+0.00062  D*0 pi+  pi-                   PHSP;
+0.0005   D*0 pi0  pi0                   PHSP;
+#
+# D a1 updated Ref. [B1]
+0.00834  a_1- D+                         SVS;
+0.0011   D+  rho0  pi-                  PHSP;
+0.0011   D+  rho-  pi0                  PHSP;
+0.0039   D+  pi+  pi-  pi-              PHSP;
+0.0010   D+  pi0  pi-  pi0              PHSP;
+0.0010   D0  pi+  pi-  pi0              PHSP;
+0.0001   D0  pi0  pi0  pi0              PHSP;
+#
+# SVV_HELAMP for D* a1 taken from factorization. Recommandation 
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/168.html:
+# updated Ref. [B1]
+# October 26, 2004 Lange update
+0.0120   D*+  a_1-                     SVV_HELAMP 0.200 0.0 0.866 0.0 0.458 0.0;
+#
+0.0025   D*+  rho0  pi-                 PHSP;
+0.0010   D*+  rho-  pi0                 PHSP;
+0.0000   D*+  pi+  pi-  pi-             PHSP;
+0.0010   D*+  pi0  pi-  pi0             PHSP;
+0.0010   D*0  pi+  pi-  pi0             PHSP;
+0.0001   D*0  pi0  pi0  pi0             PHSP;
+#
+# B->D** pi and D** rho, recommendation
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html
+#further updates October 26, 2004 - Lange
+# October 26, 2004 Lange update
+#
+0.0006   D_1+ pi-                      SVS;
+0.0001   D'_1+ pi-                     SVS;
+0.0002   D_0*+ pi- PHSP;
+0.0009   D_2*+ pi-                     STS;
+0.0004   D_1+ rho-                     PHSP;  
+0.0013   D'_1+ rho-                    PHSP;  
+0.0022   D_2*+ rho-                    PHSP;  
+#
+#
+# B->DK, recommendation 
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html:
+# update: Ref. [B1]:
+#
+0.00021  D*+ K-                        SVS;
+0.00020  D+  K-                        PHSP;
+0.00033  D*+ K*-               SVV_HELAMP 0.228 0.0 0.932 0.0 0.283 0.0;
+0.00045  K*- D+                SVS;
+#more DK modes - October 26, 2004 -Lange
+0.00004 D*0 K0 SVS;
+0.00004 D0 K0 PHSP;
+0.000042 anti-K*0 D0 SVS;
+0.00001 anti-K*0 anti-D0 SVS;
+0.00004 D*0 anti-K*0 SVV_HELAMP 1. 0. 1. 0. 1. 0.;
+0.00001 anti-D*0 anti-K*0 SVV_HELAMP 1. 0. 1. 0. 1. 0.;
+#
+#
+# Color-suppressed modes. Br's are expectations from 
+# Phys. Rev. D57 (1998), 5363-5369, except for the eta' modes, which are set
+# equal to the eta modes. SVV_HELAMP parameters are the same as in D*0 rho-,
+# with phases set to 0. Recommendation http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/169.html:
+# update: See Ref [B1]:
+# further update October 26, 2004 Lange
+#
+0.000261 D0 pi0                  PHSP;
+0.00017  D*0 pi0                 SVS;
+#
+0.00029  rho0 D0                 SVS;
+0.00029  D*0 rho0                SVV_HELAMP 0.228 0.0 0.932 0.0 0.283 0.0;
+#
+0.000202 D0 eta                  PHSP;
+0.00018  D*0 eta                 SVS;
+#
+0.000125 D0 eta'                 PHSP;
+0.000123 D*0 eta'                SVS;
+#
+0.000259 omega D0                SVS;
+0.00027  D*0 omega               SVV_HELAMP 0.228 0.0  0.932 0.0 0.283 0.0;
+#
+#
+#October 26, 2004 - Lange
+0.00097 D_s0*- D+ PHSP;
+0.00097 D*+ D_s0*- SVS;
+0.0015 D_s1- D+ SVS;
+0.0093 D*+ D_s1- SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.;
+#
+ 0.0066        cd_0       anti-ud_0         PYTHIA 22;
+ 0.0134        cd_1       anti-ud_1         PYTHIA 22;
+ 0.004         cd_0       anti-ud_0         PYTHIA 23;
+ 0.008         cd_1       anti-ud_1         PYTHIA 23;
+ 0.001         cs_0       anti-ud_0         PYTHIA 23;
+ 0.002         cs_1       anti-ud_1         PYTHIA 23;
+ 0.005         cs_0       anti-cd_0         PYTHIA 22;
+ 0.010         cs_1       anti-cd_1         PYTHIA 22;
+ 0.001         cs_0       anti-cd_0         PYTHIA 23;
+ 0.002         cs_1       anti-cd_1         PYTHIA 23;
+
+0.25885        anti-u       d  c       anti-d       PYTHIA         48;
+0.04600         anti-u       d  c       anti-d       PYTHIA         13;
+0.02100         anti-u       s  c       anti-d       PYTHIA         13;
+0.01800         anti-u       c  d       anti-d       PYTHIA         48;
+0.00090         anti-u       c  s       anti-d       PYTHIA         48;
+#lange - try to crank up the psi production....
+0.07400         anti-c       s  c       anti-d       PYTHIA         13;
+0.00600         anti-c       d  c       anti-d       PYTHIA         13;
+0.00300         anti-u       d  u       anti-d       PYTHIA         48;
+0.00400         anti-c       s  u       anti-d       PYTHIA         48;
+# JGS 11/5/02 This and similar a few lines above have been divided by two
+# to solve a double-counting problem for this channel 
+0.00205         anti-u     u  d       anti-d       PYTHIA         48;
+0.00007         anti-d     d  d       anti-d       PYTHIA         48;
+0.00009         anti-s     s  d       anti-d       PYTHIA         48;
+0.00220         anti-u     u  s       anti-d       PYTHIA         48;
+0.00180         anti-d     d  s       anti-d       PYTHIA         48;
+0.00150         anti-s     s  s       anti-d       PYTHIA         48;
+0.00500                    s  g       anti-d       PYTHIA         33;
+Enddecay
+#
+#
+
+Decay B0
+#                   b -> c semileptonic
+# =10.4% total  non-D and D* = 3.4% (changing to 2.7)
+0.0529   D*-    e+   nu_e           PHOTOS HQET 0.77 1.33 0.92;
+0.0208   D-     e+   nu_e           PHOTOS ISGW2;
+0.0052   D_1-   e+   nu_e           PHOTOS ISGW2;
+0.0045   D_0*-   e+   nu_e           PHOTOS ISGW2;
+0.0083   D'_1-   e+   nu_e           PHOTOS ISGW2;
+0.0023   D_2*-   e+   nu_e           PHOTOS ISGW2;
+0.0003   D*-  pi0   e+   nu_e       PHOTOS GOITY_ROBERTS;
+0.0065   anti-D*0  pi-   e+   nu_e       PHOTOS GOITY_ROBERTS;
+0.0010   D-   pi0   e+   nu_e       PHOTOS GOITY_ROBERTS;
+0.0032   anti-D0  pi-   e+   nu_e       PHOTOS GOITY_ROBERTS;
+0.0529   D*-    mu+  nu_mu           PHOTOS HQET 0.77 1.33 0.92;
+0.0208   D-     mu+  nu_mu           PHOTOS ISGW2;
+0.0052   D_1-   mu+  nu_mu           PHOTOS ISGW2;
+0.0045   D_0*-   mu+  nu_mu          PHOTOS  ISGW2;
+0.0083   D'_1-   mu+  nu_mu          PHOTOS  ISGW2;
+0.0023   D_2*-   mu+  nu_mu          PHOTOS  ISGW2;
+0.0003   D*-  pi0   mu+  nu_mu       PHOTOS GOITY_ROBERTS;
+0.0065   anti-D*0  pi-   mu+  nu_mu   PHOTOS     GOITY_ROBERTS;
+0.0010   D-   pi0   mu+  nu_mu       PHOTOS GOITY_ROBERTS;
+0.0032   anti-D0  pi-   mu+  nu_mu    PHOTOS    GOITY_ROBERTS;
+#
+#                    b -> c tau nu
+#
+0.0202   D*-    tau+ nu_tau                  ISGW2; # taken from Belle- Phys. Rev. Lett. 99, 191807 (2007) 
+0.0088   D-     tau+ nu_tau                  ISGW2;
+0.0013   D_1-   tau+ nu_tau                  ISGW2;
+0.0013   D_0*-   tau+ nu_tau                  ISGW2;
+0.0020   D'_1-   tau+ nu_tau                  ISGW2;
+0.0020   D_2*-   tau+ nu_tau                  ISGW2;
+#
+#                    b -> u l nu
+#
+0.000136   pi-    e+   nu_e        PHOTOS  ISGW2;
+0.000220   rho-   e+   nu_e        PHOTOS  ISGW2;
+0.000000   D(2S)- e+   nu_e        PHOTOS  ISGW2;
+0.000000   D*(2S)- e+   nu_e       PHOTOS  ISGW2;
+0.001892   Xu-  e+  nu_e           PHOTOS  VubHybridB0;
+#
+0.000136   pi-    mu+  nu_mu       PHOTOS  ISGW2;
+0.000220   rho-   mu+  nu_mu       PHOTOS  ISGW2;
+0.000000   D(2S)- mu+  nu_mu       PHOTOS  ISGW2;
+0.000000   D*(2S)- mu+  nu_mu      PHOTOS  ISGW2;
+0.001892   Xu-  mu+  nu_mu         PHOTOS  VubHybridB0;
+#
+0.000060   pi-    tau+ nu_tau                ISGW2;
+0.000083   rho-   tau+ nu_tau                ISGW2;
+0.000090   a_1-    tau+ nu_tau                ISGW2;
+0.000008   b_1-    tau+ nu_tau                ISGW2;
+0.000008   a_0-    tau+ nu_tau                ISGW2;
+0.000008   a_2-    tau+ nu_tau                ISGW2;
+0.000000   rho(2S)- tau+ nu_tau                ISGW2;
+0.000000   D(2S)- tau+ nu_tau                ISGW2;
+0.000000   D*(2S)- tau+ nu_tau                ISGW2;
+#
+
+#
+# b->u hadronic 
+# Ref. [B1]:
+# Lange Nov14,2004 (flip D_s K --opposite of D_s pi)
+0.000015   D_s+  pi-   PHSP;
+0.000028   D_s-  K+    PHSP;
+0.000029   D_s*+  pi-  SVS;
+0.000022   D_s*-  K+   SVS;
+#
+0.000016   rho-   D_s+    SVS; 
+0.000016   K*+   D_s-     SVS;
+0.000016   D_s*+  rho-    SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+0.000016   D_s*-  K*+     SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+
+
+#
+# b -> s gamma
+#
+0.0000401   K*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000135  K_10  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000065  K'_10  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000128  K'*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000166  K_2*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000017  K''*0  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+0.0003139   Xsd    gamma  BTOXSGAMMA 2;
+#
+0.000000134 K0   e+    e-    PHOTOS BTOSLLBALL;
+0.000001045 K*0  e+    e-    PHOTOS BTOSLLBALL;
+0.000004994 Xsd  e+    e-    PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000574 K0   mu+   mu-   PHOTOS BTOSLLBALL;
+0.000001100 K*0  mu+   mu-   PHOTOS BTOSLLBALL;
+0.000002506 Xsd  mu+   mu-   PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000130 K0   tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000190 K*0  tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000200 Xsd  tau+  tau-  PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+#
+#------------------------------------------------------------------------------
+# 2- and 3-body modes revised Feb.2005  Markus Cristinziani, SLAC
+#vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
+#
+#    B0 -> CP eigenstates  (some exclusive b -> u)   sum=0.00036
+#
+0.0000049 pi+  pi-     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000015 pi0  pi0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000015 pi0  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  pi0  f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+# pi0 rho0 is with the 3-body modes
+0.000001  pi0  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  pi0  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta  rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  eta' rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  eta' b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 a_00    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  a_00 f_0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  a_00 rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  a_00 b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.0000003  f_0   f_0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+0.000001  f_0  rho0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  omega   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  a_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000001  f_0  b_10    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+
+# Penguin dominated modes are sin2beta, not sin2alpha
+
+0.0000043 phi  K_S0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0000043 phi  K_L0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000001  eta K_S0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000001  eta K_L0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000032  eta' K_S0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000032  eta' K_L0    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.0000023 omega K_S0   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0000023 omega K_L0   SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+#don't allow jetset to simulate the same decay but to a K0 or anti-K0
+#Lange 1/30/2003
+0.0  phi  K0    PHSP;
+0.0  phi  anti-K0    PHSP;
+0.0  eta K0     PHSP;
+0.0  eta anti-K0 PHSP;
+0.0  eta' K0    PHSP;
+0.0  eta' anti-K0 PHSP;
+0.0  omega K0   PHSP;
+0.0  omega anti-K0 PHSP;
+# phipi0 has no simple Pengiun or tree so who knows - ditto phieta/eta'?
+
+0.00000005 phi  pi0     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; 
+0.00000003 phi  eta     SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #upper bound in PDG 6*10^-7    
+0.00000003 phi  eta'    SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; #bound 10^-6
+#
+0.0000163  K*0 eta      SVS;
+0.0000038  K*0 eta'     SVS;
+
+0.000001  K*0 omega    SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi  omega   SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  omega omega  SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  omega rho0   SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+#
+#2-body modes from Jim Olsen
+0.0000188  K+  pi-      PHSP;
+0.0000001      K+  K-       PHSP;# upper bound 3.7 10^-7
+0.0000006 K_S0 K_S0    PHSP;
+0.0000000 K_S0 K_L0    PHSP;
+0.0000006 K_L0 K_L0    PHSP;
+0.0000115 K0 pi0       PHSP;
+# Put in with 0 BR to keep PYTHIA from generating it
+0.0000000 K0 anti-K0   PHSP;
+#
+# 3-body John Back (jback@slac.stanford.edu) - Oct 15, 2002
+# JGS intersperses modes with pi0->eta,eta'
+# B0 modes
+# Model generates rho(770), rho(1450) and rho(1700) resonances 
+# with interference (uses CKM angle alpha)
+
+0.000024     pi- pi+ pi0            BTO3PI_CP dm  alpha;
+# Put in with 0 BR to keep PYTHIA from generating these (see BTO3PI_CP above)
+0.000000     rho0 pi0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.000000     rho+ pi-               SVS;
+0.000000     rho- pi+               SVS;
+0.000000     rho(2S)+ pi-           SVS;
+0.000000     rho(2S)- pi+           SVS;
+0.000000     rho(3S)+ pi-           SVS;
+0.000000     rho(3S)- pi+           SVS;
+#
+# rho0 3-body modes
+0.0000022    rho0 K_S0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.0000022    rho0 K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+#
+# rho- 3-body modes
+0.0000085    rho- K+                SVS;
+#
+# rho(1450) 3-body modes
+0.0000005    rho(2S)0 K_S0          SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma -1.0 minusGamma;
+0.0000005    rho(2S)0 K_L0          SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 gamma 1.0 minusGamma;
+#
+#
+# f0,a0(980) 3-body modes
+0.0000055    f_0 K0                 PHSP;
+0.000001     a_00 K0                PHSP;
+0.000001     a_0- K+                PHSP;
+0.000001     a_0+ pi-               PHSP;
+0.000003     a_0- pi+               PHSP;
+#
+# K*0(892) 3-body modes
+0.0000017    K*0 pi0                SVS; 
+0.000003     K*0 anti-K0            SVS;
+#
+# K*0(1430) 3-body modes
+0.000001     K_0*0 pi0              PHSP;
+0.000003     K_0*0 anti-K0          PHSP;
+#
+# K*0(1680) 3-body modes
+0.000001     K''*0 pi0              SVS;
+0.000001     K''*0 anti-K0               SVS;
+
+#
+# K*+(892) 3-body modes
+0.0000085    K*+ pi-                SVS;
+0.0000002     K*+ K-                 SVS;
+#
+# K*+(1430) 3-body modes
+0.000033     K_0*+ pi-              PHSP; #rescalled by 2/3 the Belle result PhysRevD 75 012006
+0.0000002     K_0*+ K-               PHSP;
+#
+# K*0(1680) 3-body modes
+0.000001     K''*+ pi-              SVS;
+0.0000001     K''*+ K-               SVS;
+#
+#
+# Non-resonant 3-body left-overs
+# Most modes are set at 1e-6
+#
+# pi+pi-pi0: high mass + f0(400-1200) = (1+1)e-6
+#0402270.000002     pi+ pi- pi0            PHSP;
+0.000017     pi+ pi- eta            PHSP;
+0.000001     pi+ pi- eta'           PHSP;
+#
+# pi+ pi- K0: 
+0.000013     pi+ pi- K0             PHSP;
+
+#
+# K+ pi- pi0
+0.000002     K+ pi- pi0             PHSP;
+0.000020     K+ pi- eta             PHSP;
+0.000001     K+ pi- eta'            PHSP;
+#
+# pi+ K- K0. 
+# Both B0 and B0bar can decay to these => 0.5e-6 each
+0.0000005    pi+ K- K0              PHSP;
+0.0000005    pi- K+ anti-K0         PHSP;
+
+# K+ K- K0:  total-phiK  = 24.7-4.2  
+0.0000205    K+ K- K0               PHSP;
+#
+# K+ K- pi0
+0.000001     K+ K- pi0              PHSP;
+0.000001     K+ K- eta              PHSP;
+0.000001     K+ K- eta'             PHSP;
+#
+# pi0 K0 anti-K0
+0.000002     K0 anti-K0 pi0         PHSP;
+0.000001     K0 anti-K0 eta         PHSP;
+0.000001     K0 anti-K0 eta'        PHSP;
+#
+# pi0 pi0 K0
+0.000002     K0 pi0 pi0             PHSP;
+0.000001     K0 pi0 eta             PHSP;
+0.000001     K0 pi0 eta'            PHSP;
+0.000001     K0 eta eta             PHSP;
+0.000001     K0 eta eta'            PHSP;
+0.000001     K0 eta' eta'           PHSP;
+#
+# K0 anti-K0 K0: 8*BF(KsKsKs)~=48 (8 is a guess)
+0.000048     K0 anti-K0 K0          PHSP;
+#
+# pi0 pi0 pi0
+0.000001     pi0 pi0 pi0            PHSP;
+0.000001     pi0 pi0 eta            PHSP;
+0.000001     pi0 pi0 eta'           PHSP;
+0.000001     pi0 eta eta            PHSP;
+0.000001     pi0 eta eta'           PHSP;
+#
+#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# 2- and 3-body modes revised Feb.2005   Markus Cristinziani, SLAC
+#---------------------------------------------------------------------------
+#
+#4-body modes from Andrei Gritsan
+#--- 4-body rho-rho, rho-pi-pi, pi-pi-pi-pi --------------------------
+0.0000011 rho0 rho0              SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000025  rho+ rho-              SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000005  rho0 pi+ pi-           PHSP;
+0.000003  omega pi+ pi-          PHSP;
+0.000001  rho0 pi0 pi0           PHSP;
+0.000001  omega pi0 pi0          PHSP;
+0.0000005 rho0 pi0 eta           PHSP;
+0.0000002 rho0 eta eta           PHSP;
+0.0000002 rho0 pi0 eta'          PHSP;
+0.0000001 rho0 eta eta'          PHSP;
+0.0000005 omega pi0 eta          PHSP;
+0.0000002 omega eta eta          PHSP;
+0.0000002 omega pi0 eta'         PHSP;
+0.0000001 omega eta eta'         PHSP;
+0.000010  rho+ pi- pi0           PHSP;
+0.000005  rho+ pi- eta           PHSP;
+0.000010  rho- pi+ pi0           PHSP;
+0.000005  rho- pi+ eta           PHSP;
+0.000010  pi+ pi- pi+ pi-        PHSP;
+0.000010  pi+ pi- pi0 pi0        PHSP;
+0.000005  pi+ pi- eta pi0        PHSP;
+0.000002  pi+ pi- eta eta        PHSP;
+0.000002  pi+ pi- eta' pi0       PHSP;
+0.000001  pi+ pi- eta' eta       PHSP;
+0.0000165 a_1+ pi-               SVS;
+0.0000165 a_1- pi+               SVS;
+0.000000  b_1+ pi-               SVS;
+0.000010  b_1- pi+               SVS;
+0.000002  rho- a_0+              SVS;
+0.000010  rho+ a_0-              SVS;
+0.000001  f_0 pi+ pi-            PHSP;
+0.000001  f_0 pi0 pi0            PHSP;
+0.000001  a_0+ pi- pi0           PHSP;
+0.000001  a_00 pi0 pi0           PHSP;
+0.000001  a_0- pi+ pi0           PHSP;
+#--- 4-body rho-K*, rho-pi-K, K*-pi-pi, pi-pi-pi-K -------------------
+0.0000056 rho0  K*0              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  rho-  K*+              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0; 
+0.000010  pi- pi+ K*0            PHSP;
+0.000010  pi- pi0 K*+            PHSP;
+0.000005  pi- eta K*+            PHSP;
+0.000002  pi- eta' K*+           PHSP;
+0.000010  rho-  K+ pi0           PHSP;
+0.000005  rho-  K+ eta           PHSP;
+0.000002  rho-  K+ eta'          PHSP;
+0.000010  rho-  K0 pi+           PHSP;
+0.000010  rho+  K0 pi-           PHSP;
+0.000005  rho0  K0 pi0           PHSP;
+0.000002  rho0  K0 eta           PHSP;
+0.000001  rho0  K0 eta'          PHSP;
+0.000010  rho0  K+ pi-           PHSP;
+0.000005  omega K0 pi0           PHSP;
+0.000002  omega K0 eta           PHSP;
+0.000001  omega K0 eta'          PHSP;
+0.000005  omega K+ pi-           PHSP;
+0.000010  pi+ pi0 pi- K0         PHSP;
+0.000005  pi+ eta pi- K0         PHSP;
+0.000010  pi0 pi0 pi0 K0         PHSP;
+0.000005  pi0 pi0 eta K0         PHSP;
+0.000002  pi0 eta eta K0         PHSP;
+0.000002  pi0 pi0 eta' K0        PHSP;
+0.000001  pi0 eta eta' K0        PHSP;
+0.000010  pi+ pi- pi- K+         PHSP;
+0.000010  pi0 pi0 pi- K+         PHSP;
+0.000005  pi0 eta pi- K+         PHSP;
+0.000002  eta eta pi- K+         PHSP;
+0.000002  pi0 eta' pi- K+        PHSP;
+0.000001  eta eta' pi- K+        PHSP;
+0.000010  rho0 K_0*0             PHSP;
+0.000020  rho- K_0*+             PHSP;
+0.000010  pi- pi+ K_0*0          PHSP;
+0.000010  pi0 pi0 K_0*0          PHSP;
+0.000010  pi- pi0 K_0*+          PHSP;
+0.000002  K*0 f_0                SVS;
+0.000010  a_10 K0                SVS;
+0.000010  a_1- K+                SVS; 
+0.0000066 b_10 K0                SVS;
+0.000010  b_1- K+                SVS;
+0.000005  K*0 a_00               SVS;
+0.000005  K*+ a_0-               SVS;
+#--- 4-body K*-K*, rho-K-K, K*-K-pi, phi-pi-pi, pi-pi-K-K ------------
+0.00000105  K*0   anti-K*0         SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  K*+   K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi   rho0             SVV_CP alpha dm 1 1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi   pi- pi+          PHSP;
+0.000001  K+ K- rho0             PHSP;
+0.000001  K0 anti-K0 rho0        PHSP;
+0.000001  K+ K- omega            PHSP;
+0.000001  K0 anti-K0 omega       PHSP;
+
+0.000001  K- K0 rho+             PHSP;
+0.000001  K+ anti-K0 rho-        PHSP;
+0.000001  K+ K- pi- pi+          PHSP;
+0.000001  K+ K- pi0 pi0          PHSP;
+0.0000005 K+ K- pi0 eta          PHSP;
+0.0000002 K+ K- eta eta          PHSP;
+0.0000002 K+ K- pi0 eta'         PHSP;
+0.0000001 K+ K- eta eta'         PHSP;
+0.000001  K0 anti-K0 pi- pi+     PHSP;
+0.000001  K0 anti-K0 pi0 pi0     PHSP;
+0.0000005 K0 anti-K0 pi0 eta     PHSP;
+0.0000002 K0 anti-K0 eta eta     PHSP;
+0.0000002 K0 anti-K0 pi0 eta'    PHSP;
+0.0000001 K0 anti-K0 eta eta'    PHSP;
+0.000001  K*+ K- pi0             PHSP;
+0.000001  K*- K+ pi0             PHSP;
+0.000001  K*0 anti-K0 pi0        PHSP;
+0.000001  K0 anti-K*0 pi0        PHSP;
+0.000001  K+ anti-K*0 pi-        PHSP;
+0.000001  K*+ anti-K0 pi-        PHSP;
+0.000001  K- K*0 pi+             PHSP;
+0.000001  K*- K0 pi+             PHSP;
+#--- 4-body phi-K*, phi-K-pi, K-K-K*, pi-K-K-K -----------------------
+0.0000095 phi   K*0              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  phi K+ pi-             PHSP;
+0.000010  phi K0 pi0             PHSP;
+0.000007  phi K0 eta             PHSP;
+0.000010  K- K+ K*0              PHSP;
+0.000010  K*- K+ K0              PHSP;
+0.000010  K- K*+ K0              PHSP;
+
+0.000010  K+ K- K+ pi-           PHSP;
+0.000010  K+ K- K0 pi0           PHSP;
+0.000005  K+ K- K0 eta           PHSP;
+0.000010  K0 anti-K0 K0 pi0      PHSP;
+0.000005  K0 anti-K0 K0 eta      PHSP;
+0.0000046 phi   K_0*0            PHSP;
+0.000010  K- K+ K_0*0            PHSP;
+0.000010  K_0*- K+ K0            PHSP;
+0.000010  K- K_0*+ K0            PHSP;
+
+0.000002  phi   K'_10            PHSP;
+0.000002  K- K+ K'_10            PHSP;
+#--- 4-body phi-phi, phi-K-K, K-K-K-K --------------------------------
+0.0000001 phi phi                SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi K+ K-              PHSP;
+0.0000001 phi K0 anti-K0         PHSP;
+0.0000001 K+ K- K+ K-            PHSP;
+0.0000001 K0 anti-K0 K+ K-       PHSP;
+0.0000001 K0 anti-K0 K0 anti-K0  PHSP;
+#5-body modes from Jim Smith - 1/30/05
+#--- 5-body rho-rho-pi, rho-pi-pi-pi, pi-pi-pi-pi-pi -------------------
+0.000010  rho0 rho0 pi0          PHSP;
+0.000005  rho0 rho0 eta          PHSP;
+0.000002  rho0 rho0 eta'         PHSP;
+0.000010  omega rho0 pi0         PHSP;
+0.000005  omega rho0 eta         PHSP;
+0.000002  omega rho0 eta'        PHSP;
+0.000010  rho+ rho- pi0          PHSP;
+0.000005  rho+ rho- eta          PHSP;
+0.000002  rho+ rho- eta'         PHSP;
+0.000010  rho0 rho+ pi-          PHSP; 
+0.0000066 rho0 rho-  pi+          PHSP;
+0.000010  omega rho+ pi-         PHSP;
+0.0000066 omega rho- pi+         PHSP;
+0.000002  rho0 pi+ pi- pi0       PHSP;
+0.000001  rho0 pi+ pi- eta       PHSP;
+0.000001  rho0 pi+ pi- eta'      PHSP;
+0.000002  rho0 pi0 pi0 pi0       PHSP;
+0.000001  rho0 pi0 eta pi0       PHSP;
+0.000001  rho0 pi0 eta' pi0      PHSP;
+0.000002  omega pi+ pi- pi0      PHSP;
+0.000001  omega pi+ pi- eta      PHSP;
+0.000001  omega pi+ pi- eta'     PHSP;
+0.000002  omega pi0 pi0 pi0      PHSP;
+0.000001  omega pi0 eta pi0      PHSP;
+0.000001  omega pi0 eta' pi0     PHSP;
+0.000002  rho+ pi- pi+ pi-       PHSP;
+0.000002  rho+ pi- pi0 pi0       PHSP;
+0.000001  rho+ pi- pi0 eta       PHSP;
+0.000001  rho+ pi- pi0 eta'      PHSP;
+0.000002  rho- pi+ pi+ pi-       PHSP;
+0.000002  rho- pi+ pi0 pi0       PHSP;
+0.000001  rho- pi+ pi0 eta       PHSP;
+0.000001  rho- pi+ pi0 eta'      PHSP;
+0.000001  pi+ pi- pi+ pi- pi0    PHSP; 
+0.0000005 pi+ pi- pi+ pi- eta    PHSP;
+0.0000002 pi+ pi- pi+ pi- eta'   PHSP;
+0.000001  pi+ pi- pi0 pi0 pi0    PHSP;
+0.0000005 pi+ pi- pi0 pi0 eta    PHSP;
+0.0000002 pi+ pi- pi0 pi0 eta'   PHSP;
+0.0000005 pi0 pi0 pi0 pi0 pi0    PHSP;
+0.0000002 pi0 pi0 pi0 pi0 eta    PHSP;
+0.0000001 pi0 pi0 pi0 pi0 eta'   PHSP;
+0.000010  a_1+  rho-             PHSP; 
+0.000010  a_1-  rho+             PHSP;
+0.000010  a_10  rho0             PHSP; 
+0.000010  a_10  omega            PHSP;
+0.000010  a_1+  pi- pi0          PHSP;
+0.000005  a_1+  pi- eta          PHSP;
+0.000002  a_1+  pi- eta'         PHSP;
+0.000010  a_10  pi+ pi-          PHSP;
+0.000010  a_10  pi0 pi0          PHSP;
+0.000005  a_10  pi0 eta          PHSP;
+0.000002  a_10  pi0 eta'         PHSP;
+0.000010  a_1-  pi+ pi0          PHSP;
+0.000005  a_1-  pi+ eta          PHSP;
+0.000002  a_1-  pi+ eta'         PHSP;
+0.000010  rho+ f_0 pi-           PHSP;
+0.000010  rho- f_0 pi+           PHSP;
+0.000010  rho0 f_0 pi0           PHSP;
+0.000005  rho0 f_0 eta           PHSP;
+0.000002  rho0 f_0 eta'          PHSP;
+0.000010  omega f_0 pi0          PHSP;
+0.000010  rho+ a_00 pi-          PHSP;
+0.000010  rho+ a_0- pi0          PHSP;
+0.000005  rho+ a_0- eta          PHSP;
+0.000002  rho+ a_0- eta'         PHSP;
+0.000010  rho0 a_0+ pi-          PHSP;
+0.000010  rho0 a_0- pi+          PHSP;
+0.000010  rho0 a_00 pi0          PHSP;
+0.000005  rho0 a_00 eta          PHSP;
+0.000002  rho0 a_00 eta'         PHSP;
+0.000010  rho- a_00 pi+          PHSP;
+0.000010  rho- a_0+ pi0          PHSP;
+0.000005  rho- a_0+ eta          PHSP;
+0.000002  rho- a_0+ eta'         PHSP;
+0.000002  f_0 pi+ pi- pi0        PHSP;
+0.000001  f_0 pi+ pi- eta        PHSP;
+0.000001  f_0 pi+ pi- eta'       PHSP;
+0.000002  f_0 pi0 pi0 pi0        PHSP;
+0.000001  f_0 pi0 pi0 eta        PHSP;
+0.000001  f_0 pi0 pi0 eta'       PHSP;
+0.000001  a_00 pi+ pi- pi0       PHSP;
+0.000001  a_00 pi0 pi0 pi0       PHSP;
+0.000001  a_0+ pi- pi+ pi-       PHSP;
+0.000001  a_0+ pi- pi0 pi0       PHSP;
+0.000001  a_0- pi+ pi+ pi-       PHSP;
+0.000001  a_0- pi+ pi0 pi0       PHSP;
+#--- 5-body rho-K*pi, rho-2pi-K, K*-3pi, 4pi-K -------------------
+0.000010  rho0 K*+ pi-           PHSP;
+0.000010  rho- K*+ pi0           PHSP;
+0.000005  rho- K*+ eta           PHSP;
+0.000002  rho- K*+ eta'          PHSP;
+0.000010  rho0 K*0 pi0           PHSP;
+0.000005  rho0 K*0 eta           PHSP;
+0.000002  rho0 K*0 eta'          PHSP;
+0.000010  rho- K*0 pi+           PHSP;
+0.000010  rho+ K*0 pi-           PHSP;
+0.000010  omega K*+ pi-          PHSP;
+0.000010  omega K*0 pi0          PHSP;
+0.000005  omega K*0 eta          PHSP;
+0.000002  omega K*0 eta'         PHSP;
+0.000010  pi+ pi- K*+ pi-        PHSP;
+0.000010  pi- pi0 K*+ pi0        PHSP;
+0.000005  pi- pi0 K*+ eta        PHSP;
+0.000002  pi- pi0 K*+ eta'       PHSP;
+0.000010  pi+ pi- K*0 pi0        PHSP;
+0.000005  pi+ pi- K*0 eta        PHSP;
+0.000002  pi+ pi- K*0 eta'       PHSP;
+0.000010  pi0 pi0 K*0 pi0        PHSP;
+0.000005  pi0 pi0 K*0 eta        PHSP;
+0.000002  pi0 pi0 K*0 eta'       PHSP;
+0.000010  rho+  K+ pi- pi-       PHSP;
+0.000010  rho0  K+ pi- pi0       PHSP;
+0.000005  rho0  K+ pi- eta       PHSP;
+0.000002  rho0  K+ pi- eta'      PHSP;
+0.000010  rho-  K+ pi+ pi-       PHSP;
+0.000010  rho-  K+ pi0 pi0       PHSP;
+0.000005  rho-  K+ pi0 eta       PHSP;
+0.000002  rho-  K+ pi0 eta'      PHSP;
+0.000010  rho0  K0 pi+ pi-       PHSP;
+0.000010  rho0  K0 pi0 pi0       PHSP;
+0.000005  rho0  K0 pi0 eta       PHSP;
+0.000002  rho0  K0 pi0 eta'      PHSP;
+0.000010  rho+  K0 pi- pi0       PHSP;
+0.000005  rho+  K0 pi- eta       PHSP;
+0.000002  rho+  K0 pi- eta'      PHSP;
+0.000010  rho-  K0 pi+ pi0       PHSP;
+0.000005  rho-  K0 pi+ eta       PHSP;
+0.000002  rho-  K0 pi+ eta'      PHSP;
+0.000010  omega K+ pi- pi0       PHSP;
+0.000005  omega K+ pi- eta       PHSP;
+0.000002  omega K+ pi- eta'      PHSP;
+0.000010  omega K0 pi+ pi-       PHSP;
+0.000010  omega K0 pi0 pi0       PHSP;
+0.000005  omega K0 pi0 eta       PHSP;
+0.000002  omega K0 pi0 eta'      PHSP;
+0.000002  pi+ pi- pi- K+ pi0     PHSP;
+0.000001  pi+ pi- pi- K+ eta     PHSP;
+0.000001  pi+ pi- pi- K+ eta'    PHSP;
+0.000002  pi0 pi0 pi- K+ pi0     PHSP;
+0.000001  pi0 pi0 pi- K+ eta     PHSP;
+0.000001  pi0 pi0 pi- K+ eta'    PHSP;
+0.000002  pi+ pi- pi+ K0 pi-     PHSP;
+0.000002  pi+ pi- pi0 K0 pi0     PHSP;
+0.000001  pi+ pi- pi0 K0 eta     PHSP;
+0.000001  pi+ pi- pi0 K0 eta'    PHSP;
+0.000002  pi0 pi0 pi0 K0 pi0     PHSP;
+0.000001  pi0 pi0 pi0 K0 eta     PHSP;
+0.000001  pi0 pi0 pi0 K0 eta'    PHSP;
+0.000010  rho+ K_0*0 pi-         PHSP;
+0.000010  rho0 K_0*0 pi0         PHSP;
+0.000005  rho0 K_0*0 eta         PHSP;
+0.000002  rho0 K_0*0 eta'        PHSP;
+0.000010  rho- K_0*0 pi+         PHSP;
+0.000010  rho- K_0*+ pi0         PHSP;
+0.000005  rho- K_0*+ eta         PHSP;
+0.000002  rho- K_0*+ eta'        PHSP;
+0.000010  rho0 K_0*+ pi-         PHSP;
+0.000010  pi+ pi- K_0*+ pi-      PHSP;
+0.000010  pi- pi0 K_0*+ pi0      PHSP;
+0.000005  pi- pi0 K_0*+ eta      PHSP;
+0.000002  pi- pi0 K_0*+ eta'     PHSP;
+0.000010  pi+ pi- K_0*0 pi0      PHSP;
+0.000005  pi+ pi- K_0*0 eta      PHSP;
+0.000002  pi+ pi- K_0*0 eta'     PHSP;
+0.000010  pi0 pi0 K_0*0 pi0      PHSP;
+0.000005  pi0 pi0 K_0*0 eta      PHSP;
+0.000005  pi0 pi0 K_0*0 eta'     PHSP;
+0.000010  K*+ f_0 pi-            PHSP;
+0.000010  K*0 f_0 pi0            PHSP;
+0.000005  K*0 f_0 eta            PHSP;
+0.000002  K*0 f_0 eta'           PHSP;
+0.000010  a_1+ K0 pi-            PHSP;
+0.000020  a_10 K*0               PHSP;
+0.000010  a_10 K+ pi-            PHSP;
+0.000010  a_10 K0 pi0            PHSP;
+0.000020  a_1- K*+               PHSP;
+0.000010  a_1- K+ pi0            PHSP;
+0.000010  a_1- K0 pi+            PHSP;
+0.000005  K*+ a_00 pi-           PHSP;
+0.000005  K*+ a_0- pi0           PHSP;
+0.000002  K*+ a_0- eta           PHSP;
+0.000001  K*+ a_0- eta'          PHSP;
+0.000005  K*0 a_0- pi+           PHSP;
+0.000005       K*0 a_0+ pi-      PHSP;
+0.000005  K*0 a_00 pi0           PHSP;
+0.000002  K*0 a_00 eta           PHSP;
+0.000001  K*0 a_00 eta'          PHSP;
+#--- 5-body K*-K*pi,rho-K-Kpi, K*-K-2pi, phi-3pi, 3pi-K-K ------------
+0.000001  K*- K*0 pi+            PHSP;
+0.000001  K*+ anti-K*0 pi-       PHSP;
+0.000001  K*- K*+ pi0            PHSP;
+0.0000005 K*- K*+ eta            PHSP;
+0.0000002 K*- K*+ eta'           PHSP;
+0.000001  K*0 anti-K*0 pi0       PHSP;
+0.0000005 K*0 anti-K*0 eta       PHSP;
+0.0000002 K*0 anti-K*0 eta'      PHSP;
+0.000001  phi rho+ pi-           PHSP;
+0.000001  phi rho- pi+           PHSP;
+0.000001  phi rho0 pi0           PHSP;
+0.0000005 phi rho0 eta           PHSP;
+0.0000002 phi rho0 eta'          PHSP;
+0.000001  phi omega pi0          PHSP;
+0.000001  phi f_0 pi0            PHSP;
+0.000001  phi pi+ pi- pi0        PHSP;
+0.0000005 phi pi+ pi- eta        PHSP;
+0.0000002 phi pi+ pi- eta'       PHSP;
+0.000001  phi pi0 pi0 pi0        PHSP;
+0.0000005 phi pi0 pi0 eta        PHSP;
+0.0000002 phi pi0 pi0 eta'       PHSP;
+0.000001  anti-K0 K+ rho0 pi-    PHSP;
+0.000001  anti-K0 K+ rho- pi0    PHSP;
+0.0000005 anti-K0 K+ rho- eta    PHSP;
+0.0000002 anti-K0 K+ rho- eta'   PHSP;
+0.000001  anti-K0 K+ omega pi-   PHSP;
+0.000001  anti-K0 K+ f_0 pi-     PHSP;
+
+
+0.000001  K+ K- rho- pi+         PHSP;
+0.000001  K+ K- rho+ pi-         PHSP;
+0.000001  K+ K- rho0 pi0         PHSP;
+0.0000005 K+ K- rho0 eta         PHSP;
+0.0000002 K+ K- rho0 eta'        PHSP;
+0.000001  K+ K- omega pi0        PHSP;
+0.0000005 K+ K- omega eta        PHSP;
+0.0000002 K+ K- omega eta'       PHSP;
+0.000001  K+ K- f_0 pi0          PHSP;
+0.000001  K0 anti-K0 rho- pi+    PHSP;
+0.000001  K0 anti-K0 rho+ pi-    PHSP;
+0.000001  K0 anti-K0 rho0 pi0    PHSP;
+0.0000005 K0 anti-K0 rho0 eta    PHSP;
+0.0000002 K0 anti-K0 rho0 eta'   PHSP;
+0.000001  K0 anti-K0 omega pi0   PHSP;
+0.000001  K0 anti-K0 f_0 pi0     PHSP;
+0.000001  K0 K- rho0 pi+         PHSP;
+0.000001  K0 K- rho+ pi0         PHSP;
+0.0000005 K0 K- rho+ eta         PHSP;
+0.0000002 K0 K- rho+ eta'        PHSP;
+0.000001  K0 K- omega pi+        PHSP;
+0.000001  K0 K- f_0 pi+          PHSP;
+0.0000002 K0 K- pi+ pi+ pi-      PHSP;
+0.0000002 K0 K- pi+ pi0 pi0      PHSP;
+0.0000001 K0 K- pi+ pi0 eta      PHSP;
+0.0000001 K0 K- pi+ pi0 eta'     PHSP;
+0.0000002 K+ K- pi+ pi- pi0      PHSP;
+0.0000001 K+ K- pi+ pi- eta      PHSP;
+0.0000001 K+ K- pi+ pi- eta'     PHSP;
+0.0000002 K+ K- pi0 pi0 pi0      PHSP;
+0.0000001 K+ K- pi0 pi0 eta      PHSP;
+0.0000001 K+ K- pi0 pi0 eta'     PHSP;
+0.0000002 K0 anti-K0 pi+ pi- pi0 PHSP;
+0.0000001 K0 anti-K0 pi+ pi- eta PHSP;
+0.0000001 K0 anti-K0 pi+ pi- eta' PHSP;
+0.0000002 K0 anti-K0 pi0 pi0 pi0 PHSP;
+0.0000001 K0 anti-K0 pi0 pi0 eta PHSP;
+0.0000001 K0 anti-K0 pi0 pi0 eta' PHSP;
+0.0000002 anti-K0 K+ pi- pi+ pi- PHSP;
+0.0000002 anti-K0 K+ pi- pi0 pi0 PHSP;
+0.0000001 anti-K0 K+ pi- pi0 eta PHSP;
+0.0000001 anti-K0 K+ pi- pi0 eta' PHSP;
+0.000001  K*0 K- pi+ pi0         PHSP;
+0.0000005 K*0 K- pi+ eta         PHSP;
+0.0000002 K*0 K- pi+ eta'        PHSP;
+0.000001  K*- K0 pi+ pi0         PHSP;
+0.0000005 K*- K0 pi+ eta         PHSP;
+0.0000002 K*- K0 pi+ eta'        PHSP;
+0.000001  K*+ K- pi+ pi-         PHSP;
+0.000001  K*+ K- pi0 pi0         PHSP;
+0.0000005 K*+ K- pi0 eta         PHSP;
+0.0000002 K*+ K- pi0 eta'        PHSP;
+0.000001  K*- K+ pi+ pi-         PHSP;
+0.000001  K*- K+ pi0 pi0         PHSP;
+0.0000005 K*- K+ pi0 eta         PHSP;
+0.0000002 K*- K+ pi0 eta'        PHSP;
+0.000001  K*0 anti-K0 pi+ pi-    PHSP;
+0.000001  K*0 anti-K0 pi0 pi0    PHSP;
+0.0000005 K*0 anti-K0 pi0 eta    PHSP;
+0.0000002 K*0 anti-K0 pi0 eta'   PHSP;
+0.000001  K0 anti-K*0 pi+ pi-    PHSP;
+0.000001  K0 anti-K*0 pi0 pi0    PHSP;
+0.0000005 K0 anti-K*0 pi0 eta    PHSP;
+0.0000002 K0 anti-K*0 pi0 eta'   PHSP;
+0.000001  K+ anti-K*0 pi- pi0    PHSP;
+0.0000005 K+ anti-K*0 pi- eta    PHSP;
+0.0000002 K+ anti-K*0 pi- eta'   PHSP;
+0.000001  K*+ anti-K0 pi- pi0    PHSP;
+0.0000005 K*+ anti-K0 pi- eta    PHSP;
+0.0000002 K*+ anti-K0 pi- eta'   PHSP;
+#--- 5-body phi-K*-pi, phi-K-2pi, K-K-K*-pi, 2pi-K-K-K -------------------- 
+0.000010  phi K*+ pi-            PHSP;
+0.000010  phi K*0 pi0            PHSP;
+0.000005  phi K*0 eta            PHSP;
+0.000002  phi K*0 eta'           PHSP;
+0.000010  phi K+ rho-            PHSP;
+0.000010  phi K0 rho0            PHSP;
+0.000010  phi K0 omega           PHSP;
+0.000010  phi K0 f_0             PHSP;
+0.000010  phi K0 pi+ pi-         PHSP;
+0.000010  phi K0 pi0 pi0         PHSP;
+0.000005  phi K0 pi0 eta         PHSP;
+0.000002  phi K0 pi0 eta'        PHSP;
+0.000010  phi K+ pi- pi0         PHSP;
+0.000005  phi K+ pi- eta         PHSP;
+0.000002  phi K+ pi- eta'        PHSP;
+0.000010  K- K+ K*+ pi-          PHSP;
+0.000010  K+ K*- K+ pi-          PHSP;
+0.000010  K+ K*0 anti-K0 pi-     PHSP;
+0.000010  K+ K0 anti-K*0 pi-     PHSP;
+0.000010  K0 K*+ anti-K0 pi-     PHSP;
+0.000010  K- K+ K+ rho-          PHSP;
+0.000010  K- K+ K0 rho0          PHSP;
+0.000010  K- K+ K0 omega         PHSP;
+0.000010  K- K+ K0 f_0           PHSP;
+0.000010  K+ K- K+ pi- pi0       PHSP;
+0.000005  K+ K- K+ pi- eta       PHSP;
+0.000002  K+ K- K+ pi- eta'      PHSP;
+0.000010  K+ K- K0 pi+ pi-       PHSP;
+0.000010  K+ K- K0 pi0 pi0       PHSP;
+0.000005  K+ K- K0 pi0 eta       PHSP;
+0.000005  K+ K- K0 pi0 eta'      PHSP;
+0.000010  K0 anti-K0 K+ pi- pi0  PHSP;
+0.000005  K0 anti-K0 K+ pi- eta  PHSP;
+0.000002  K0 anti-K0 K+ pi- eta' PHSP;
+0.000010  K0 anti-K0 K0 pi+ pi-  PHSP;
+0.000010  K0 anti-K0 K0 pi0 pi0  PHSP;
+0.000005  K0 anti-K0 K0 pi0 eta  PHSP;
+0.000002  K0 anti-K0 K0 pi0 eta' PHSP;
+0.000010  phi   K_0*+ pi-        PHSP;
+0.000010  phi   K_0*0 pi0        PHSP;
+0.000010  K+ K- K_0*+ pi-        PHSP;
+0.000010  K+ K- K_0*0 pi0        PHSP;
+0.000005  K+ K- K_0*0 eta        PHSP;
+0.000002  K+ K- K_0*0 eta'       PHSP;
+0.000010  K+ K+ K_0*- pi-        PHSP;
+
+0.000010  K+ K_0*0 anti-K0 pi-   PHSP;
+0.000010  K+ anti-K_0*0 K0 pi-   PHSP;
+0.000010  K0 anti-K0 K_0*+ pi-   PHSP;
+0.000010  K0 anti-K0 K_0*0 pi0   PHSP;
+0.000005  K0 anti-K0 K_0*0 eta   PHSP;
+0.000002  K0 anti-K0 K_0*0 eta'  PHSP;
+0.000010  K0 K0 anti-K_0*0 pi0   PHSP;
+0.000005  K0 K0 anti-K_0*0 eta   PHSP;
+0.000002  K0 K0 anti-K_0*0 eta'  PHSP;
+0.000002  phi   K'_1+ pi-        PHSP;
+0.000002  phi   K'_10 pi0        PHSP;
+0.000001  phi   K'_10 eta        PHSP;
+0.000001  phi   K'_10 eta'       PHSP;
+0.000002  K- K+ K'_1+ pi-        PHSP;
+0.000002  K- K+ K'_10 pi0        PHSP;
+0.000001  K- K+ K'_10 eta        PHSP;
+0.000001  K- K+ K'_10 eta'       PHSP;
+0.000002  K+ K+ K'_1- pi-        PHSP;
+0.000002  K0 anti-K0 K'_10 pi0   PHSP;
+0.000001  K0 anti-K0 K'_10 eta   PHSP;
+0.000001  K0 anti-K0 K'_10 eta'  PHSP;
+0.000002  K0 K0 anti-K'_10 pi0   PHSP;
+0.000001  K0 K0 anti-K'_10 eta   PHSP;
+0.000001  K0 K0 anti-K'_10 eta'  PHSP;
+#--- 5-body phi-phi-K, phi-phi-pi, phi-3K, phi-K-K-pi, 5K, 4K-pi ------------
+0.0000041 phi phi K0             PHSP;
+0.0000001 phi phi pi0            PHSP;
+0.000001  phi K+ K- K0           PHSP;
+0.000001  phi K0 anti-K0 K0      PHSP;
+0.0000001 phi K+ K- pi0          PHSP;
+0.0000001 phi K0 anti-K0 pi0     PHSP;
+0.0000001 phi K0 K- pi+          PHSP;
+0.0000001 phi anti-K0 K+ pi-     PHSP;
+
+0.0000001 K+ K- K+ K- K0         PHSP;
+0.0000001 K0 anti-K0 K+ K- K0    PHSP;
+0.0000001 K+ K- K+ K- pi0        PHSP;
+0.0000001 K+ K- K0 anti-K0 pi0   PHSP;
+0.0000001 K+ K- K0 K- pi+        PHSP;
+0.0000001 K+ K- K+ anti-K0 pi-   PHSP;
+
+0.0000001 K0 anti-K0 K0 anti-K0 pi0 PHSP;
+0.0000001 K0 anti-K0 K0 K- pi+   PHSP;
+0.0000001 K0 anti-K0 anti-K0 K+ pi- PHSP;
+
+#--- 6-body phi-phi-K*, phi-3Kpi, phi-K-K-2pi ------------  not in PDG
+0.000003  phi phi K*0            PHSP;
+0.000001  phi K+ K- K*0          PHSP;
+0.000001  phi K0 anti-K0 K*0     PHSP;
+0.000001  phi K0 K+ K*-          PHSP;
+0.000001  phi K0 K0 anti-K*0     PHSP;
+0.0000001 phi K+ K- K+ pi-       PHSP;
+0.0000001 phi K+ K- K0 pi0       PHSP;
+0.0000001 phi K0 anti-K0 K+ pi-  PHSP;
+0.0000001 phi K0 anti-K0 K0 pi0  PHSP;
+#--- 6-body a1a1 ------------
+0.000050  a_10  a_10             PHSP;
+0.000050  a_1+  a_1-             PHSP;
+#
+#       B -> cc= s          
+0.000436    J/psi  K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000436    J/psi  K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.0002215   J/psi  K*S               SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.0002215   J/psi  K*L               SVV_CP beta  dm -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000886     J/psi  K*0T              SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus;
+0.0000205    J/psi  pi0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000016     J/psi  rho0             SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00003     J/psi  omega            SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.0012     J/psi  K+  pi-           PHSP;
+0.0001     J/psi  K0  pi0           PHSP;
+#rl0.0007     J/psi  K0  pi-  pi+      PHSP;
+#rl0.00035     J/psi  K0  pi0  pi0      PHSP;
+#rl0.00035     J/psi  K+  pi-  pi0      PHSP;
+0.0013     J/psi  K_10              SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0001     J/psi  K'_10             SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0005     J/psi  K_2*0              PHSP;
+0.000094    J/psi  phi K0        PHSP;
+#
+0.00031     psi(2S)  K_S0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00031     psi(2S)  K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00012    psi(2S)  K*S              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00012    psi(2S)  K*L              SVV_CP beta  dm  -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00048     psi(2S)  K*0T             SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus;
+
+0.0004     psi(2S)  K+  pi-          PHSP;
+0.0002     psi(2S)  K0  pi0          PHSP;
+0.0002     psi(2S)  K0  pi-  pi+     PHSP;
+0.0001     psi(2S)  K0  pi0  pi0     PHSP;
+0.0001     psi(2S)  K+  pi-  pi0     PHSP;
+0.0004     psi(2S)  K_10              PHSP;
+#
+0.0005    eta_c K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0005    eta_c K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00027   K*S  eta_c              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00027   K*L  eta_c              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00108    K*0T eta_c              SVS;
+0.0002    eta_c  K+  pi-          PHSP;
+0.0001    eta_c  K0  pi0          PHSP;
+0.0002    eta_c  K0  pi-  pi+     PHSP;
+0.0001    eta_c  K0  pi0  pi0     PHSP;
+0.0001    eta_c  K+  pi-  pi0     PHSP;
+#
+0.00024    eta_c(2S) K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00024    eta_c(2S) K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00011   K*S  eta_c(2S)               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00011   K*L  eta_c(2S)               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00044   K*0T eta_c(2S)               SVS;
+0.00008    eta_c(2S)   K+  pi-          PHSP;
+0.00005    eta_c(2S)   K0  pi0          PHSP;
+0.00008    eta_c(2S)   K0  pi-  pi+     PHSP;
+0.00005    eta_c(2S)   K0  pi0  pi0     PHSP;
+0.00005    eta_c(2S)   K+  pi-  pi0     PHSP;
+#
+0.00005     chi_c0 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00005     chi_c0 K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00005    K*S  chi_c0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00005    K*L  chi_c0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.0002     K*0T chi_c0              SVS;
+0.0002     chi_c0  K+  pi-          PHSP;
+0.0001     chi_c0  K0  pi0          PHSP;
+0.0002     chi_c0  K0  pi-  pi+     PHSP;
+0.0001     chi_c0  K0  pi0  pi0     PHSP;
+0.0001     chi_c0  K+  pi-  pi0     PHSP;
+#
+0.00020     chi_c1 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00020     chi_c1 K_L0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.000053    K*S  chi_c1              SVV_CP beta  dm  1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000053    K*L  chi_c1              SVV_CP beta  dm -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.000213    chi_c1  K*0T             SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus;
+0.0004     chi_c1  K+  pi-          PHSP;
+0.0002     chi_c1  K0  pi0          PHSP;
+0.0004     chi_c1  K0  pi-  pi+     PHSP;
+0.0002     chi_c1  K0  pi0  pi0     PHSP;
+0.0002     chi_c1  K+  pi-  pi0     PHSP;
+#
+0.00001     chi_c2 K_S0               SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0; 
+0.00001     chi_c2 K_L0              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.000005    K*S  chi_c2              PHSP;
+0.000005    K*L  chi_c2              PHSP;
+0.00002     chi_c2  K*0T             PHSP;
+0.0002     chi_c2  K+  pi-          PHSP;
+0.0001     chi_c2  K0  pi0          PHSP;
+0.0002     chi_c2  K0  pi-  pi+     PHSP;
+0.0001     chi_c2  K0  pi0  pi0     PHSP;
+0.0001     chi_c2  K+  pi-  pi0     PHSP;
+
+0.00024    psi(3770)  K_S0             SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 -1.0 0.0;
+0.00024    psi(3770)  K_L0             SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00008    psi(3770)  K*S              SVV_CP beta  dm 1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00008    psi(3770)  K*L              SVV_CP beta  dm  -1 Aplus phAplus Azero phAzero Aminus phAminus;
+0.00032    psi(3770)  K*0T             SVV_HELAMP PKHminus PKphHminus PKHzero PKphHzero PKHplus PKphHplus;
+0.00014    psi(3770)  K+  pi-          PHSP;
+0.00014    psi(3770)  K0  pi0          PHSP;
+0.00014    psi(3770)  K0  pi-  pi+     PHSP;
+0.00007    psi(3770)  K0  pi0  pi0     PHSP;
+0.00007    psi(3770)  K+  pi-  pi0     PHSP;
+0.00029    psi(3770)  K_10             PHSP;
+#
+
+
+
+#
+# b-> c (dc=)  
+#
+0.00023     D-   D+              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+#  See Ref [B1] for the next 3:
+0.000610    D*-  D+              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00061     D*+  D-              SSD_CP dm 0.0 1.0 minusTwoBeta 1.0 0.0 1.0 0.0;
+0.00082     D*- D*+              SVV_CP beta dm  1  0.56 0.0 0.96 0.0 0.47 0.0;
+#
+# b -> c (sc=)  -->  D Ds X   Sum = 10%
+#
+0.0065   D-   D_s+                       PHSP;
+0.0080   D*-  D_s+                       SVS;
+0.0074   D_s*+ D-                        SVS;
+0.0177   D_s*+ D*-                 SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+0.0006   D'_1- D_s+                SVS;
+0.0012   D'_1-  D_s*+              SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0012   D_1-  D_s+                SVS;
+0.0024   D_1-  D_s*+               SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0042   D_2*-  D_s+                      STS;
+0.0040   D_2*-  D_s*+                    PHSP;
+#
+0.0018   D_s+  D-  pi0                   PHSP;
+0.0037   D_s+  anti-D0 pi-                   PHSP;
+0.0018   D_s*+ D-  pi0                   PHSP;
+0.0037   D_s*+ anti-D0 pi-                   PHSP;
+0.0030   D_s+  D-  pi-  pi+              PHSP;
+0.0022   D_s+  D-  pi0  pi0              PHSP;
+0.0022   D_s+  anti-D0 pi-  pi0              PHSP;
+0.0030   D_s*+ D-  pi-  pi+              PHSP;
+0.0022   D_s*+ D-  pi0  pi0              PHSP;
+0.0022   D_s*+ anti-D0 pi-  pi0              PHSP;
+#
+# b -> c (sc=)  -->  D D= K  X   Sum = 8%
+# Update: Ref. [B1]
+# October 26, 2004 Lange -per breco awg
+# External W-emission amplitude
+0.0017   D-  D0  K+                         PHSP;
+0.0046   D-  D*0  K+                        PHSP;
+0.0031   D*- D0  K+                         PHSP;
+0.0118   D*- D*0  K+                        PHSP;
+# External+internal W-emission amplitude
+0.0015   D-  D+  K0                         PHSP; 
+0.0018   D*- D+  K0                         PHSP; 
+0.0047   D-  D*+  K0                        PHSP; 
+0.0088   D*- D*+  K0                        PHSP;
+# Internal W-emission amplitude (color suppressed modes) 
+0.0005   D0  anti-D0  K0                    PHSP;
+0.0005   D*0  anti-D0  K0                   PHSP; 
+0.0005   D0  anti-D*0  K0                   PHSP; 
+0.0005   D*0  anti-D*0  K0                  PHSP; 
+##########
+# K* modes
+###########  
+0.0025   D-  D0  K*+                        PHSP;
+0.0025   D*- D0  K*+                        PHSP;
+0.0025   D-  D*0  K*+                       PHSP;
+0.0050   D*- D*0  K*+                       PHSP;
+#
+0.0025   D-  D+  K*0                        PHSP;
+0.0025   D*- D+  K*0                        PHSP;
+0.0025   D-  D*+  K*0                       PHSP;
+0.0050   D*- D*+  K*0                       PHSP;
+#
+0.0005   D0  anti-D0  K*0                   PHSP;
+0.0005   D0  anti-D*0  K*0                  PHSP;
+0.0005   D*0  anti-D0  K*0                  PHSP;
+0.0010   D*0  anti-D*0  K*0                 PHSP;
+#
+#  B -> D(*) X  Exclusive modes
+#
+0.00276    D*- pi+                        SVS;
+0.00268  D-  pi+                        PHSP;
+0.0075   rho+ D-                        SVS;
+#
+# D* rho HELAMP parameters taken from ICHEP 98-852.
+0.0068   rho+ D*-                 SVV_HELAMP  0.317 0.19 0.936 0.0 0.152 1.47; 
+0.0005   D-  pi+  pi0                   PHSP;
+0.0082   D*- pi+  pi0                   PHSP;
+0.0005   anti-D0 pi-  pi+                   PHSP;
+0.00062   anti-D*0 pi-  pi+                   PHSP;
+0.0005   anti-D*0 pi0  pi0                   PHSP;
+#
+# D a1 updated Ref. [B1]
+0.00834  a_1+ D-                         SVS; 
+0.0011   D-  rho0  pi+                  PHSP;
+0.0011   D-  rho+  pi0                  PHSP;
+0.0039   D-  pi-  pi+  pi+              PHSP;
+0.0010   D-  pi0  pi+  pi0              PHSP;
+0.0010   anti-D0  pi-  pi+  pi0             PHSP;
+0.0001   anti-D0  pi0  pi0  pi0             PHSP;
+#
+# SVV_HELAMP for D* a1 taken from factorization. Recommandation 
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/168.html:
+# updated Ref. [B1]
+# October 26, 2004 Lange update
+#
+0.0120   D*-  a_1+                     SVV_HELAMP 0.458 0.0 0.866 0.0 0.458 0.0;
+#
+0.0025    D*-  rho0  pi+                 PHSP;
+0.0010   D*-  rho+  pi0                 PHSP;
+0.0000   D*-  pi-  pi+  pi+             PHSP;
+0.0010   D*-  pi0  pi+  pi0             PHSP;
+0.0010   anti-D*0  pi-  pi+  pi0             PHSP;
+0.0001   anti-D*0  pi0  pi0  pi0             PHSP;
+#
+# B->D** pi and D** rho, recommendation
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html
+#
+0.0006   D_1- pi+                      SVS;
+0.0001   D'_1- pi+                     SVS;
+0.0002   D_0*- pi+                      PHSP;
+0.0009   D_2*- pi+                     STS;
+0.0004   D_1- rho+                     PHSP;  
+0.0013   D'_1- rho+                    PHSP;  
+0.0022   D_2*- rho+                    PHSP;  
+#
+#
+# B->DK, recommendation 
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html:
+# update: Ref. [B1]:
+#
+0.00021  D*- K+                        SVS;
+0.00020  D-  K+                        PHSP;
+0.00033  D*- K*+               SVV_HELAMP  0.283 0.0 0.932 0.0 0.228 0.0;
+0.00045  K*+ D-                SVS;
+#more DK modes - October 26, 2004 -Lange
+0.00004 anti-D*0 anti-K0 SVS;
+0.00004 anti-D0 anti-K0 PHSP;
+0.000042 K*0 anti-D0 SVS;
+0.00001 K*0 D0 SVS;
+0.00004 anti-D*0 K*0 SVV_HELAMP 1. 0. 1. 0. 1. 0.;
+0.00001 D*0 K*0 SVV_HELAMP 1. 0. 1. 0. 1. 0.;
+#
+#
+# Color-suppressed modes. Br's are expectations from 
+# Phys. Rev. D57 (1998), 5363-5369, except for the eta' modes, which are set
+# equal to the eta modes. SVV_HELAMP parameters are the same as in D*0 rho-,
+# with phases set to 0. Recommendation http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/169.html:
+# update: See Ref [B1]:
+# further update October 26, 2004 Lange
+#
+0.000261  anti-D0 pi0                  PHSP;
+0.00017  anti-D*0 pi0                 SVS;
+#
+0.00029  rho0 anti-D0                 SVS;
+0.00029  anti-D*0 rho0                SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
+#
+0.000202 anti-D0 eta                  PHSP;
+0.00018  anti-D*0 eta                 SVS;
+#
+0.000125 anti-D0 eta'                 PHSP;
+0.000123 anti-D*0 eta'                SVS;
+#
+0.000259 omega anti-D0                SVS;
+0.00027  anti-D*0 omega               SVV_HELAMP 0.283 0.0 0.932 0.0 0.228 0.0;
+#
+       #October 26, 2004 - Lange 
+0.00097 D_s0*+ D- PHSP;
+0.00097  D*- D_s0*+ SVS;
+0.0015 D_s1+ D- SVS;
+0.0093 D*- D_s1+ SVV_HELAMP 0.4904 0. 0.7204 0. 0.4904 0.;
+#
+0.0066        anti-cd_0       ud_0         PYTHIA 22;
+0.0134        anti-cd_1       ud_1         PYTHIA 22;
+0.004         anti-cd_0       ud_0         PYTHIA 23;
+0.008         anti-cd_1       ud_1         PYTHIA 23;
+0.001         anti-cs_0       ud_0         PYTHIA 23;
+0.002         anti-cs_1       ud_1         PYTHIA 23;
+0.005         anti-cs_0       cd_0         PYTHIA 22;
+0.010         anti-cs_1       cd_1         PYTHIA 22;
+0.001         anti-cs_0       cd_0         PYTHIA 23;
+0.002         anti-cs_1       cd_1         PYTHIA 23;
+#
+#
+0.25885         u       anti-d  anti-c          d       PYTHIA         48;
+0.04600         u       anti-d  anti-c          d       PYTHIA         13;
+0.02100         u       anti-s  anti-c          d       PYTHIA         13;
+0.01800         u       anti-c  anti-d          d       PYTHIA         48;
+0.00090         u       anti-c  anti-s          d       PYTHIA         48;
+#lange - try to crank up the psi production....
+0.07400         c       anti-s  anti-c          d       PYTHIA         13;
+0.00600         c       anti-d  anti-c          d       PYTHIA         13;
+0.00300         u       anti-d  anti-u          d       PYTHIA         48;
+0.00400         c       anti-s  anti-u          d       PYTHIA         48;
+# JGS 11/5/02 This and similar a few lines above have been divided by two
+# to solve a double-counting problem for this channel 
+0.00205         u     anti-u  anti-d       d       PYTHIA         48;
+0.00007         d     anti-d  anti-d       d       PYTHIA         48;
+0.00009         s     anti-s  anti-d       d       PYTHIA         48;
+0.00220         u     anti-u  anti-s       d       PYTHIA         48;
+0.00180         d     anti-d  anti-s       d       PYTHIA         48;
+0.00150         s     anti-s  anti-s       d       PYTHIA         48;
+0.00500              anti-s  g            d       PYTHIA         33;
+Enddecay
+#
+#
+#
+Decay B-
+#
+#                  b -> c semileptonic
+#
+0.0650   D*0    e-   anti-nu_e         PHOTOS  HQET 0.77 1.33 0.92;
+0.0215   D0     e-   anti-nu_e         PHOTOS  ISGW2;
+0.0056   D_10   e-   anti-nu_e         PHOTOS  ISGW2;
+0.0049   D_0*0   e-   anti-nu_e         PHOTOS  ISGW2;
+0.009   D'_10   e-   anti-nu_e         PHOTOS  ISGW2;
+0.003   D_2*0   e-   anti-nu_e         PHOTOS  ISGW2;
+0.0063   D*+  pi-   e-   anti-nu_e     PHOTOS  GOITY_ROBERTS;
+0.0003   D*0  pi0   e-   anti-nu_e     PHOTOS  GOITY_ROBERTS;
+0.0052   D+   pi-   e-   anti-nu_e     PHOTOS  GOITY_ROBERTS;
+0.0010   D0   pi0   e-   anti-nu_e     PHOTOS  GOITY_ROBERTS;
+#
+0.0650   D*0   mu-  anti-nu_mu         PHOTOS   HQET 0.77 1.33 0.92;
+0.0215   D0    mu-  anti-nu_mu         PHOTOS   ISGW2;
+0.0056   D_10  mu-  anti-nu_mu         PHOTOS   ISGW2;
+0.0049   D_0*0  mu-  anti-nu_mu        PHOTOS    ISGW2;
+0.009   D'_10  mu-  anti-nu_mu        PHOTOS    ISGW2;
+0.003   D_2*0  mu-  anti-nu_mu        PHOTOS    ISGW2;
+0.0063   D*+   pi-   mu-   anti-nu_mu  PHOTOS   GOITY_ROBERTS;
+0.0003   D*0   pi0   mu-   anti-nu_mu  PHOTOS   GOITY_ROBERTS;
+0.0052   D+    pi-   mu-   anti-nu_mu  PHOTOS   GOITY_ROBERTS;
+0.0010   D0    pi0   mu-   anti-nu_mu  PHOTOS   GOITY_ROBERTS;
+#
+#                   b -> c tau nu
+#
+0.0208   D*0    tau- anti-nu_tau                 ISGW2;# equel to B0-> D*- tau+ nu
+0.0088   D0     tau- anti-nu_tau                 ISGW2;
+0.0013   D_10   tau- anti-nu_tau                 ISGW2;
+0.0013   D_0*0   tau- anti-nu_tau                 ISGW2;
+0.0020   D'_10   tau- anti-nu_tau                 ISGW2;
+0.0020   D_2*0   tau- anti-nu_tau                 ISGW2;
+#
+#                   b -> u l nu
+#
+# NOTE: Do NOT CHANGE any BFs without using the corresponding set of 
+#       hybrid weights (and vice versa). 
+
+0.000077   pi0    e-   anti-nu_e        PHOTOS   ISGW2;
+0.000080   eta    e-   anti-nu_e        PHOTOS   ISGW2;
+0.000128   rho0   e-   anti-nu_e        PHOTOS   ISGW2;
+0.000130   omega   e-   anti-nu_e       PHOTOS   ISGW2;
+0.000084   eta'   e-   anti-nu_e        PHOTOS   ISGW2;
+0.000000   D(2S)0 e-   anti-nu_e        PHOTOS   ISGW2;
+0.000000   D*(2S)0 e-   anti-nu_e       PHOTOS   ISGW2;
+0.001948   Xu0   e-  anti-nu_e          PHOTOS   VubHybridBp;
+#
+0.000077   pi0    mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000080   eta    mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000128   rho0   mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000130   omega   mu-  anti-nu_mu      PHOTOS   ISGW2;
+0.000084   eta'   mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000000   D(2S)0 mu-  anti-nu_mu       PHOTOS   ISGW2;
+0.000000   D*(2S)0 mu-  anti-nu_mu      PHOTOS   ISGW2;
+0.001948   Xu0   mu-  anti-nu_mu        PHOTOS   VubHybridBp;
+#
+0.000030   pi0    tau- anti-nu_tau                 ISGW2;
+0.000012   eta    tau- anti-nu_tau                 ISGW2;
+0.000042   rho0   tau- anti-nu_tau                 ISGW2;
+0.000042   omega   tau- anti-nu_tau                 ISGW2;
+0.000018   eta'   tau- anti-nu_tau                 ISGW2;
+0.000046   a_10    tau- anti-nu_tau                 ISGW2;
+0.000027   b_10    tau- anti-nu_tau                 ISGW2;
+0.000004   a_00    tau- anti-nu_tau                 ISGW2;
+0.000002   f_0     tau- anti-nu_tau                 ISGW2;
+0.000002   f'_0    tau- anti-nu_tau                 ISGW2;
+0.000023   f_1     tau- anti-nu_tau                 ISGW2;
+0.000023   f'_1    tau- anti-nu_tau                 ISGW2;
+0.000013   h_1     tau- anti-nu_tau                 ISGW2;
+0.000013   h'_1    tau- anti-nu_tau                 ISGW2;
+0.000020   f_2     tau- anti-nu_tau                 ISGW2;
+0.000020   f'_2    tau- anti-nu_tau                 ISGW2;
+0.000000   D(2S)0 tau- anti-nu_tau                 ISGW2;
+0.000000   D*(2S)0 tau- anti-nu_tau                 ISGW2;
+#
+#   b -> l nu
+#
+0.000180   tau- anti-nu_tau                        SLN;
+0.000000   mu-  anti-nu_mu               PHOTOS   SLN;
+0.000000   e-   anti-nu_e                PHOTOS  SLN;
+#
+#
+#
+# b->u hadronic 
+# Ref. [B1]:
+#
+0.000016   D_s-  pi0   PHSP;
+0.000020   D_s*-  pi0  SVS;
+#
+0.000028   rho0  D_s-     SVS;
+0.000028   D_s*-  rho0    SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+
+#
+# b-> s gamma
+#
+0.0000403   K*-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000135  K_1-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000065  K'_1-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000128  K'*-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000166  K_2*-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000017  K''*-  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+0.0003137   anti-Xsu    gamma  BTOXSGAMMA 2;
+#
+0.00000049  K-        e+    e-    PHOTOS BTOSLLBALL;
+0.000000700 K*-       e+    e-    PHOTOS BTOSLLBALL;
+0.000004994 anti-Xsu  e+    e-    PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.00000039  K-        mu+   mu-   PHOTOS BTOSLLBALL;
+0.000000700 K*-       mu+   mu-   PHOTOS BTOSLLBALL;
+0.000002506 anti-Xsu  mu+   mu-   PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000130 K-        tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000190 K*-       tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000200 anti-Xsu  tau+  tau-  PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+#
+#
+#------------------------------------------------------------------------------
+# 2- and 3-body modes revised Feb.2005  Markus Cristinziani, SLAC
+#vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
+# Exclusive hadronic b->u
+#
+0.000026   pi-   pi0   pi0                   CB3PI-P00 alpha;
+0.0000121  pi-   pi-   pi+                   CB3PI-MPP alpha;
+####
+0.0000049 eta pi-                PHSP;
+0.0000026 eta K-                 PHSP;
+0.0000193 K*- eta                SVS;
+0.0000084 rho- eta               SVS;
+0.0000026 eta' pi-               PHSP;
+0.0000697 eta' K-                PHSP;
+0.0000049 K*- eta'               SVS;
+0.0000087 rho- eta'              SVS;
+0.0000067 omega pi-              SVS;
+0.0000068 omega K-               SVS;
+0.000001  omega K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000106 omega rho-             SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001  phi pi-                SVS;
+0.0000083 phi K-                 SVS;
+0.000023  anti-K0 pi-            PHSP;
+0.00000128 K0 K-                  PHSP;
+0.0000055 pi- pi0                PHSP;
+0.0000121  K- pi0                 PHSP;
+#
+# 3-body John Back (jback@slac.stanford.edu) - Oct 15, 2002
+# JGS intersperses modes with pi0->eta,eta'
+# B- modes
+#
+# rho0 3-body modes
+### AJR zeroed rho0pi- in favor of the CB3PI-P00 decay above
+0.0000000    rho0 pi-               SVS;
+0.0000042    rho0 K-                SVS;
+#
+# rho- 3-body modes
+### AJR zeroed rho-pi0 in favor of the CB3PI-MPP decay above
+0.000000     rho- pi0               SVS;
+0.000001     rho- anti-K0           SVS;
+
+
+#
+# rho(1450) 3-body modes
+0.0000022    rho(2S)0 pi-           SVS;
+0.000001     rho(2S)0 K-            SVS;
+#
+# f0,a0(980) 3-body modes
+0.000001     f_0 pi-                PHSP;
+0.0000177    f_0 K-                 PHSP;
+0.000001     a_0- anti-K0           PHSP;
+0.000001     a_00 K-                PHSP;
+0.000001     a_00 pi-               PHSP;
+0.000001     a_0- pi0               PHSP;
+
+0.0000082 f_2 pi-                   PHSP;
+
+#
+# K*(1430) 3-body modes
+0.000047     anti-K_0*0 pi-         PHSP;
+0.000002     K_0*0 K-               PHSP;
+0.000047     K_0*- pi0              PHSP;
+0.000002     K_0*- K0               PHSP;
+#
+# K* (892) 3-body modes
+0.0000109    anti-K*0 pi-           SVS;
+0.0000069    K*- pi0                SVS;
+0.000003     K*- K0                 SVS;
+#
+#
+#
+# Non-resonant 3-body left-overs
+# Most modes are set at 1e-6
+#
+# pi-pi+pi-: high mass + f0(400-1200) = (1+1)e-6
+#0402270.000002     pi- pi+ pi-            PHSP;
+#
+# K- pi+ pi-: high mass + f0(400-1200) = (5+5)e-6
+0.0000100     K- pi+ pi-             PHSP;
+#
+# K- K+ pi-: just non-resonant
+0.000001     K- K+ pi-              PHSP;
+#
+# K- K+ K-: high mass structure near 1500 + non-res: total - phiK 
+0.0000296    K- K+ K-               PHSP;
+#
+# 
+# K-K-pi+: suppressed mode (1e-7)
+0.0000001    K- K- pi+              PHSP;
+# K+pi-pi-: suppressed mode (1e-7)
+0.0000001    K+ pi- pi-             PHSP;
+#
+#
+# K- K0 pi0
+0.000001     K- K0 pi0              PHSP;
+0.000001     K- K0 eta              PHSP;
+0.000001     K- K0 eta'             PHSP;
+#
+# pi- anti-K0 pi0
+0.000001     pi- anti-K0 pi0        PHSP;
+0.000001     pi- anti-K0 eta        PHSP;
+0.000001     pi- anti-K0 eta'       PHSP;
+#
+# pi- pi0 pi0
+#0402270.000001     pi- pi0 pi0            PHSP;
+0.000001     pi- pi0 eta            PHSP;
+0.000001     pi- pi0 eta'           PHSP;
+0.000001     pi- eta eta            PHSP;
+0.000001     pi- eta eta'           PHSP;
+#
+# K- pi0 pi0
+0.000001     K- pi0 pi0             PHSP;
+0.000001     K- pi0 eta             PHSP;
+0.000001     K- pi0 eta'            PHSP;
+0.000001     K- eta eta             PHSP;
+0.000001     K- eta eta'            PHSP;
+0.000001     K- eta' eta'           PHSP;
+#
+# K- anti-K0 K0:  4*BF(K+KsKs) (4 is a guess)
+0.000046     K- anti-K0 K0          PHSP;
+#
+# pi- anti-K0 K0
+0.000001     pi- anti-K0 K0         PHSP;
+#
+#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# 2- and 3-body modes revised Feb.2005   Markus Cristinziani, SLAC
+#---------------------------------------------------------------------------
+#
+#4-body modes from Andrei Gritsan
+#--- 4-body rho-rho, rho-pi-pi, pi-pi-pi-pi --------------------------
+0.000018  rho-  rho0             SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  rho0 pi- pi0           PHSP;
+0.000005  rho0 pi- eta           PHSP;
+0.000002  rho0 pi- eta'          PHSP;
+0.000010  omega pi- pi0          PHSP;
+0.000005  omega pi- eta          PHSP;
+0.000002  omega pi- eta'         PHSP;
+0.000010  rho- pi- pi+           PHSP;
+0.000010  rho- pi0 pi0           PHSP;
+0.000005  rho- pi0 eta           PHSP;
+0.000002  rho- eta eta           PHSP;
+0.000002  rho- pi0 eta'          PHSP;
+0.000010  pi+ pi- pi- pi0        PHSP;
+0.000005  pi+ pi- pi- eta        PHSP;
+0.000002  pi+ pi- pi- eta'       PHSP;
+0.000010  pi- pi0 pi0 pi0        PHSP;
+0.000005  pi- eta pi0 pi0        PHSP;
+0.000002  pi- eta eta pi0        PHSP;
+0.000002  pi- eta' pi0 pi0       PHSP;
+0.000001  pi- eta' eta pi0       PHSP;
+0.0000204 a_10  pi-              SVS;
+0.0000264 a_1-  pi0              SVS;
+0.0000067 b_10  pi-              SVS;
+0.000001  b_1-  pi0              SVS;
+0.000002  rho- f_0               SVS;
+0.000010  rho- a_00              SVS;
+0.000002  rho0 a_0-              SVS;
+0.000010  f_0 pi- pi0            PHSP;
+0.000001  a_00 pi- pi0           PHSP;
+0.000001  a_0- pi0 pi0           PHSP;
+0.000001  a_0+ pi- pi-           PHSP;
+#--- 4-body rho-K*, rho-pi-K, K*-pi-pi, pi-pi-pi-K -------------------
+0.0000092 rho-  anti-K*0         SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000005  rho0  K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  pi- pi0 anti-K*0       PHSP;
+0.000005  pi- eta anti-K*0       PHSP;
+0.000067  pi+ pi- K*-            PHSP;
+0.000010  pi0 pi0 K*-            PHSP;
+0.000005  pi0 eta K*-            PHSP;
+0.000002  eta eta K*-            PHSP;
+0.000002  pi0 eta' K*-           PHSP;
+0.000001  eta eta' K*-           PHSP;
+0.000010  rho-  anti-K0 pi0      PHSP;
+0.000005  rho-  anti-K0 eta      PHSP;
+0.000002  rho-  anti-K0 eta'     PHSP;
+0.000010  rho-  K- pi+           PHSP;
+0.000010  rho+  K- pi-           PHSP;
+0.000005  rho0  anti-K0 pi-      PHSP;
+0.000005  rho0  K- pi0           PHSP;
+0.000002  rho0  K- eta           PHSP;
+0.000001  rho0  K- eta'          PHSP;
+0.000005  omega anti-K0 pi-      PHSP;
+0.000005  omega K- pi0           PHSP;
+0.000002  omega K- eta           PHSP;
+0.000001  omega K- eta'          PHSP;
+0.000010  pi+ pi- pi- anti-K0    PHSP;
+0.000010  pi- pi0 pi0 anti-K0    PHSP;
+0.000005  pi- eta pi0 anti-K0    PHSP;
+0.000002  pi- eta eta anti-K0    PHSP;
+0.000010  pi+ pi- pi0 K-         PHSP;
+0.000005  pi+ pi- eta K-         PHSP;
+0.000010  pi0 pi0 pi0 K-         PHSP;
+0.000005  pi0 pi0 eta K-         PHSP;
+0.000002  pi0 eta eta K-         PHSP;
+0.000002  pi0 pi0 eta' K-        PHSP;
+0.000001  pi0 eta eta' K-        PHSP;
+0.000010  rho- anti-K_0*0        PHSP;
+0.000010  rho0 K_0*-             PHSP;
+0.000010  pi- pi0 anti-K_0*0     PHSP;
+0.000010  pi0 pi0 K_0*-          PHSP;
+0.0000052 K*- f_0                SVS;
+0.000010  a_10 K-                SVS;
+0.0000349 a_1- anti-K0           SVS;
+0.0000091 b_10 K-                SVS;
+0.000013  b_1- anti-K0           SVS;
+0.000005  K*- a_00               SVS;
+0.000005  anti-K*0 a_0-          SVS;
+#--- 4-body K*-K*, rho-K-K, K*-K-pi, phi-pi-pi, pi-pi-K-K ------------
+0.000001  K*0   K*-              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi   rho-             SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi   pi- pi0          PHSP;
+0.0000005 phi   pi- eta          PHSP;
+0.0000002 phi   pi- eta'         PHSP;
+0.000001  K+ K- rho-             PHSP;
+0.000001  K0 anti-K0 rho-        PHSP;
+0.000001  K0 K- rho0             PHSP;
+0.000001  K0 K- omega            PHSP;
+0.000001  K+ K- pi- pi0          PHSP;
+0.0000005 K+ K- pi- eta          PHSP;
+0.0000002 K+ K- pi- eta'         PHSP;
+0.000001  K0 anti-K0 pi- pi0     PHSP;
+0.0000005 K0 anti-K0 pi- eta     PHSP;
+0.0000002 K0 anti-K0 pi- eta'    PHSP;
+0.000001  anti-K0 K+ pi- pi-     PHSP;#2007 added      
+0.000001  K0 K- pi+ pi-          PHSP;
+0.00000025  K0 K- pi0 pi0          PHSP;
+0.0000005 K0 K- pi0 eta          PHSP;
+0.00000025 K0 K- eta eta          PHSP;#missing from B+
+0.0000002 K0 K- pi0 eta'         PHSP;
+0.0000001 K0 K- eta eta'         PHSP;
+0.000001  K*- K+ pi-             PHSP;
+0.000001  K*+ K- pi-             PHSP;
+
+0.000001  K0 anti-K*0 pi-        PHSP;
+0.000001  K*0 anti-K0 pi-        PHSP;
+0.000001  K- K*0 pi0             PHSP;
+0.0000005 K- K*0 eta             PHSP;
+0.0000002 K- K*0 eta'            PHSP;
+0.000001  K*- K0 pi0             PHSP;
+0.0000005 K*- K0 eta             PHSP;
+0.0000002 K*- K0 eta'            PHSP;
+#--- 4-body phi-K*, phi-K-pi, K-K-K*, pi-K-K-K -----------------------
+0.0000096 phi K*-                SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  phi anti-K0 pi-        PHSP;
+0.000010  phi K- pi0             PHSP;
+0.000005  phi K- eta             PHSP;
+0.000002  phi K- eta'            PHSP;
+0.000031  K- K+ K*-              PHSP;#total -phiK*
+0.000010  K- K*+ K-              PHSP;
+
+0.000010  K- K0 anti-K*0         PHSP;
+0.000010  K- K*0 anti-K0         PHSP;
+0.000010  K0 K*- anti-K0         PHSP;
+0.000010  K+ K- K- pi0           PHSP;
+0.000005  K+ K- K- eta           PHSP;
+0.000002  K+ K- K- eta'          PHSP;
+0.000010  K+ K- anti-K0 pi-      PHSP;
+0.000010  K0 anti-K0 K- pi0      PHSP;
+0.000005  K0 anti-K0 K- eta      PHSP;
+0.000002  K0 anti-K0 K- eta'     PHSP;
+0.000010  K0 anti-K0 anti-K0 pi- PHSP;
+0.000010  phi   K_0*-            PHSP;
+0.000010  K- K+ K_0*-            PHSP;
+0.000010  K- K- K_0*+            PHSP;
+0.000010  K- anti-K_0*0 K0       PHSP;
+0.000010  K- K_0*0 anti-K0       PHSP;
+0.000002  phi   K'_1-            PHSP;
+0.000002  K- K+ K'_1-            PHSP;
+#--- 4-body phi-phi, phi-K-K, K-K-K-K --------------------------------
+0.000001  phi K0 K-              PHSP;
+0.000001  K+ K- K- K0            PHSP;
+0.000001  K0 anti-K0 K0 K-       PHSP;
+#5-body modes from Jim Smith - 1/30/05
+#--- 5-body rho-rho-pi, rho-pi-pi-pi, pi-pi-pi-pi-pi -------------------
+0.000010  rho0 rho0 pi-          PHSP;
+0.000010  rho- rho0 pi0          PHSP;
+0.000005  rho- rho0 eta          PHSP;
+0.000002  rho- rho0 eta'         PHSP;
+0.000010  rho+ rho- pi-          PHSP;
+0.000010  omega rho0 pi-         PHSP;
+0.000010  omega rho- pi0         PHSP;
+0.000005  omega rho- eta         PHSP;
+0.000002  omega rho- eta'        PHSP;
+0.000002  rho0 pi+ pi- pi-       PHSP;
+0.000002  rho0 pi- pi0 pi0       PHSP;
+0.000001  rho0 pi- eta pi0       PHSP;
+0.000001  rho0 pi- eta' pi0      PHSP;
+0.000002  omega pi- pi0 pi0      PHSP;
+0.000001  omega pi- eta pi0      PHSP;
+0.000001  omega pi- eta' pi0     PHSP;
+0.000002  rho+ pi- pi- pi0       PHSP;
+0.000001  rho+ pi- pi- eta       PHSP;
+0.000001  rho+ pi- pi- eta'      PHSP;
+0.000002  rho- pi+ pi- pi0       PHSP;
+0.000001  rho- pi+ pi- eta       PHSP;
+0.000001  rho- pi+ pi- eta'      PHSP;
+0.000002  rho- pi0 pi0 pi0       PHSP;
+0.000001  rho- pi0 eta pi0       PHSP;
+0.000001  rho- eta eta pi0       PHSP;
+0.000001  rho- pi0 eta' pi0      PHSP;
+0.000001  pi+ pi- pi- pi0 pi0    PHSP;
+0.000001  pi+ pi- pi- eta pi0    PHSP;
+0.000001  pi+ pi- pi- eta' pi0   PHSP;
+0.000001  pi- pi0 pi0 pi0 pi0    PHSP;
+0.000001  pi- eta pi0 pi0 pi0    PHSP;
+0.000001  pi- eta' pi0 pi0 pi0   PHSP;
+0.000010  a_1-  rho0             PHSP;
+0.000010  a_1-  omega            PHSP;
+0.000010  a_1-  f_0              PHSP;
+0.000010  a_10  pi- pi0          PHSP;
+0.000005  a_10  pi- eta          PHSP;
+0.000002  a_10  pi- eta'         PHSP;
+0.000010  a_1+  pi- pi-          PHSP;
+0.000010  a_1-  pi+ pi-          PHSP;
+0.000010  a_1-  pi0 pi0          PHSP;
+0.000005  a_1-  pi0 eta          PHSP;
+0.000002  a_1-  pi0 eta'         PHSP;
+0.000010  rho- f_0 pi0           PHSP;
+0.000005  rho- f_0 eta           PHSP;
+0.000002  rho- f_0 eta'          PHSP;
+0.000010  rho0 f_0 pi-           PHSP;
+0.000010  omega f_0 pi-          PHSP;
+0.000010  rho- a_00 pi0          PHSP;
+0.000005  rho- a_00 eta          PHSP;
+0.000002  rho- a_00 eta'         PHSP;
+0.000010  rho+ a_0- pi-          PHSP;
+0.000010  rho0 a_0- pi0          PHSP;
+0.000005  rho0 a_0- eta          PHSP;
+0.000002  rho0 a_0- eta'         PHSP;
+0.000002  f_0 pi+ pi- pi-        PHSP;
+0.000002  f_0 pi- pi0 pi0        PHSP;
+0.000001  f_0 pi- pi0 eta        PHSP;
+0.000001  f_0 pi- pi0 eta'       PHSP;
+0.000001  a_00 pi- pi0 pi0       PHSP;
+0.000001  a_0- pi0 pi0 pi0       PHSP;
+0.000001  a_0+ pi- pi- pi0       PHSP;
+#--- 5-body rho-K*pi, rho-2pi-K, K*-3pi, 4pi-K -------------------
+0.000010  rho- K*- pi+           PHSP;
+0.000010  rho+ K*- pi-           PHSP;
+0.000010  rho- anti-K*0 pi0      PHSP;
+0.000005  rho- anti-K*0 eta      PHSP;
+0.000002  rho- anti-K*0 eta'     PHSP;
+0.000010  rho0 anti-K*0 pi-      PHSP;
+0.000010  rho0 K*- pi0           PHSP;
+0.000005  rho0 K*- eta           PHSP;
+0.000002  rho0 K*- eta'          PHSP;
+0.000010  omega anti-K*0 pi-     PHSP;
+0.000010  omega K*- pi0          PHSP;
+0.000005  omega K*- eta          PHSP;
+0.000002  omega K*- eta'         PHSP;
+0.000010  pi+ pi- anti-K*0 pi-   PHSP;
+0.000010  pi- pi0 anti-K*0 pi0   PHSP;
+0.000005  pi- eta anti-K*0 pi0   PHSP;
+0.000002  pi- eta' anti-K*0 pi0  PHSP;
+0.000010  pi- pi+ K*- pi0        PHSP;
+0.000005  pi- pi+ K*- eta        PHSP;
+0.000002  pi- pi+ K*- eta'       PHSP;
+0.000010  pi0 pi0 K*- pi0        PHSP;
+0.000005  pi0 eta K*- pi0        PHSP;
+0.000002  eta eta K*- pi0        PHSP;
+0.000002  pi0 eta' K*- pi0       PHSP;
+0.000001  eta eta' K*- pi0       PHSP;
+0.000010  rho-  anti-K0 pi+ pi-  PHSP;
+0.000010  rho-  anti-K0 pi0 pi0  PHSP;
+0.000005  rho-  anti-K0 eta pi0  PHSP;
+0.000002  rho-  anti-K0 eta' pi0 PHSP;
+0.000010  rho-  K- pi+ pi0       PHSP;
+0.000010  rho+  K- pi- pi0       PHSP;
+0.000010  rho+  anti-K0 pi- pi-  PHSP;
+0.000010  rho0  anti-K0 pi- pi0  PHSP;
+0.000010  rho0  K- pi0 pi0       PHSP;
+0.000005  rho0  K- eta pi0       PHSP;
+0.000002  rho0  K- eta' pi0      PHSP;
+0.000010  omega anti-K0 pi- pi0  PHSP;
+0.000010  omega K- pi0 pi0       PHSP;
+0.000005  omega K- eta pi0       PHSP;
+0.000002  omega K- eta' pi0      PHSP;
+0.000002  pi- pi- pi+ anti-K0 pi0 PHSP;
+0.000001  pi- pi- pi+ anti-K0 eta PHSP;
+0.000001  pi- pi- pi+ anti-K0 eta' PHSP;
+0.000002  pi- pi0 pi0 anti-K0 pi0 PHSP;
+0.000001  pi- eta pi0 anti-K0 pi0 PHSP;
+0.000001  pi- eta eta anti-K0 pi0 PHSP;
+0.000002  pi+ pi- pi+ K- pi-     PHSP;
+0.000002  pi+ pi- pi0 K- pi0     PHSP;
+0.000001  pi+ pi- eta K- pi0     PHSP;
+0.000001  pi+ pi- eta' K- pi0    PHSP;
+0.000002  pi0 pi0 pi0 K- pi0     PHSP;
+0.000001  pi0 pi0 eta K- pi0     PHSP;
+0.000001  pi0 eta eta K- pi0     PHSP;
+0.000001  pi0 pi0 eta' K- pi0    PHSP;
+0.000001  pi0 eta eta' K- pi0    PHSP;
+0.000010  rho0 anti-K_0*0 pi-    PHSP;
+0.000010  rho- anti-K_0*0 pi0    PHSP;
+0.000005  rho- anti-K_0*0 eta    PHSP;
+0.000002  rho- anti-K_0*0 eta'   PHSP;
+0.000010  rho- K_0*- pi+         PHSP;
+0.000010  rho0 K_0*- pi0         PHSP;
+0.000005  rho0 K_0*- eta         PHSP;
+0.000002  rho0 K_0*- eta'        PHSP;
+0.000010  rho+ K_0*- pi-        PHSP;#added 2007 not in PDG exist in B+
+0.000010  pi- pi- anti-K_0*0 pi+ PHSP;
+0.000010  pi- pi0 anti-K_0*0 pi0 PHSP;
+0.000005  pi- pi0 anti-K_0*0 eta PHSP;
+0.000002  pi- pi0 anti-K_0*0 eta' PHSP;
+0.000010  pi0 pi0 K_0*- pi0      PHSP;
+0.000005  pi0 pi0 K_0*- eta      PHSP;
+0.000002  pi0 pi0 K_0*- eta'     PHSP;
+0.000010  anti-K*0 f_0 pi-       PHSP;
+0.000010  K*- f_0 pi0            PHSP;
+0.000005  K*- f_0 eta            PHSP;
+0.000002  K*- f_0 eta'           PHSP;
+0.000020  a_10 K*-               PHSP;
+0.000010  a_10 K- pi0            PHSP;
+0.000010  a_10 anti-K0 pi-       PHSP;
+0.000020  a_1- anti-K*0          PHSP;
+0.000010  a_1- K- pi+            PHSP;
+0.000010  a_1- anti-K0 pi0       PHSP;
+0.000005  K*- a_00 pi0           PHSP;
+0.000005  anti-K*0 a_0- pi0      PHSP;
+#--- 5-body K*-K*pi,rho-K-Kpi, K*-K-2pi, phi-3pi, 3pi-K-K ------------
+0.000001  K*0 anti-K*0 pi-       PHSP;
+0.000001  K*0 K*- pi0       PHSP;
+0.0000005 K*0 K*- eta       PHSP;
+0.0000002 K*0 K*- eta'      PHSP;
+0.000001  phi rho0 pi-           PHSP;
+0.000001  phi omega pi-          PHSP;
+0.000001  phi f_0 pi-            PHSP;
+0.000001  phi rho- pi0           PHSP;
+0.0000005 phi rho- eta           PHSP;
+0.0000002 phi rho- eta'          PHSP;
+0.000001  phi pi+ pi- pi-        PHSP;
+0.000001  phi pi- pi0 pi0        PHSP;
+0.0000005 phi pi- eta pi0        PHSP;
+0.0000002 phi pi- eta' pi0       PHSP;
+0.000001  K0 K- rho+ pi-         PHSP;
+0.000001  K0 K- rho- pi+         PHSP;
+0.000001  K0 K- rho0 pi0         PHSP;
+0.0000005 K0 K- rho0 eta         PHSP;
+0.0000002 K0 K- rho0 eta'        PHSP;
+0.000001  K0 K- omega pi0        PHSP;
+0.000001  K0 K- f_0 pi0          PHSP;
+0.000001  K+ K- rho0 pi-         PHSP;
+0.000001  K+ K- rho- pi0         PHSP;
+0.0000005 K+ K- rho- eta         PHSP;
+0.0000002 K+ K- rho- eta'        PHSP;
+0.000001  K+ K- omega pi-        PHSP;
+0.000001  K0 anti-K0 rho0 pi-    PHSP;
+0.000001  K0 anti-K0 omega pi-   PHSP;
+0.000001  K0 anti-K0 f_0 pi-     PHSP;
+0.000001  K0 anti-K0 rho- pi0    PHSP;
+0.0000005 K0 anti-K0 rho- eta    PHSP;
+0.0000002 K0 anti-K0 rho- eta'   PHSP;
+0.000001  anti-K0 K+ rho- pi-    PHSP;
+0.0000002 anti-K0 K+ pi- pi- pi0 PHSP;
+0.0000001 anti-K0 K+ pi- pi- eta PHSP;
+0.0000001 anti-K0 K+ pi- pi- eta' PHSP;
+0.0000002 K+ K- pi+ pi- pi-      PHSP;
+0.0000002 K+ K- pi- pi0 pi0      PHSP;
+0.0000001 K+ K- pi- eta pi0      PHSP;
+0.0000001 K+ K- pi- eta' pi0     PHSP;
+0.0000002 K0 anti-K0 pi+ pi- pi- PHSP;
+0.0000002 K0 anti-K0 pi- pi0 pi0 PHSP;
+0.0000001 K0 anti-K0 pi- eta pi0 PHSP;
+0.0000001 K0 anti-K0 pi- eta' pi0 PHSP;
+0.0000002 K0 K- pi+ pi- pi0      PHSP;
+0.0000001 K0 K- pi+ pi- eta      PHSP;
+0.0000001 K0 K- pi+ pi- eta'     PHSP;
+0.0000002 K0 K- pi0 pi0 pi0      PHSP;
+0.0000001 K0 K- pi0 eta pi0      PHSP;
+0.0000001 K0 K- pi0 eta' pi0     PHSP;
+0.000001  anti-K*0 K+ pi- pi-    PHSP;
+0.000001  K*+ anti-K0 pi- pi-    PHSP;
+0.000001  K*- K+ pi- pi0         PHSP;
+0.000001  K*+ K- pi- pi0         PHSP;
+0.000001  K*0 anti-K0 pi- pi0    PHSP;
+0.000001  K0 anti-K*0 pi- pi0    PHSP;
+0.000001  K- K*0 pi0 pi0         PHSP;
+0.0000005 K- K*0 eta pi0         PHSP;
+0.0000002 K- K*0 eta' pi0        PHSP;
+0.000001  K*- K0 pi0 pi0         PHSP;
+0.0000005 K*- K0 eta pi0         PHSP;
+0.0000002 K*- K0 eta' pi0        PHSP;
+#--- 5-body phi-K*-pi, phi-K-2pi, K-K-K*-pi, 2pi-K-K-K --------------------
+0.000010  phi K*- pi0            PHSP;
+0.000010  phi K- rho0            PHSP;
+0.000010  phi K- omega           PHSP;
+0.000010  phi K- f_0             PHSP;
+0.000010  phi anti-K0 rho-       PHSP;
+0.000010  phi anti-K0 pi- pi0    PHSP;
+0.000010  phi K- pi+ pi-         PHSP;
+0.000010  phi K- pi0 pi0         PHSP;
+0.000005  phi K- eta pi0         PHSP;
+0.000002  phi K- eta' pi0        PHSP;
+0.000010  K+ K- K*- pi0          PHSP;
+0.000010  K- K*+ K- pi0          PHSP;
+0.000010  K- K0 anti-K*0 pi0     PHSP;
+0.000010  K- K*0 anti-K0 pi0     PHSP;
+
+0.000010  K0 K*- anti-K0 pi0     PHSP;
+0.000010  K+ K- K- rho0          PHSP;
+0.000010  K+ K- K- omega         PHSP;
+0.000010  K+ K- K- f_0           PHSP;
+0.000010  K- K+ anti-K0 rho-     PHSP;
+0.000010  K- K- K+ pi0 pi0       PHSP;
+0.000005  K- K- K+ eta pi0       PHSP;
+0.000002  K- K- K+ eta' pi0      PHSP;
+0.000010  K+ K- anti-K0 pi- pi0  PHSP;
+0.000010  K0 anti-K0 K- pi0 pi0  PHSP;
+0.000005  K0 anti-K0 K- eta pi0  PHSP;
+0.000002  K0 anti-K0 K- eta' pi0 PHSP;
+0.000010  K0 anti-K0 anti-K0 pi- pi0 PHSP;
+0.000010  phi   anti-K_0*0 pi-   PHSP;
+0.000010  phi   K_0*- pi0        PHSP;
+0.000005  phi   K_0*- eta        PHSP;
+0.000002  phi   K_0*- eta'       PHSP;
+0.000010  K- K+ K_0*- pi0        PHSP;
+0.000005  K- K+ K_0*- eta        PHSP;
+0.000002  K- K+ K_0*- eta'       PHSP;
+0.000010  K- K- K_0*+ pi0        PHSP;
+0.000005  K- K- K_0*+ eta        PHSP;
+0.000002  K- K- K_0*+ eta'       PHSP;
+0.000010  K- K+ anti-K_0*0 pi-   PHSP;
+0.000010  K- K- K_0*0 pi+        PHSP;
+0.000010  anti-K0 anti-K0 K_0*0 pi- PHSP;
+0.000010  K- anti-K_0*0 K0 pi0   PHSP;
+0.000005  K- anti-K_0*0 K0 eta   PHSP;
+0.000002  K- anti-K_0*0 K0 eta'  PHSP;
+0.000010  K- K_0*0 anti-K0 pi0   PHSP;
+0.000005  K- K_0*0 anti-K0 eta   PHSP;
+0.000002  K- K_0*0 anti-K0 eta'  PHSP;
+0.000002  phi   anti-K'_10 pi-   PHSP;
+0.000002  phi   K'_1- pi0        PHSP;
+0.000001  phi   K'_1- eta        PHSP;
+0.000001  phi   K'_1- eta'       PHSP;
+0.000002  K- K+ anti-K'_10 pi-   PHSP;
+0.000002  K- K+ K'_1- pi0        PHSP;
+0.000001  K- K+ K'_1- eta        PHSP;
+0.000001  K- K+ K'_1- eta'       PHSP;
+0.000002  K- K- K'_1+ pi0        PHSP;
+0.000002  K- K- K'_10 pi+   PHSP; #2007 changed anti-K'_10 to K'_10 
+0.000002  K0 anti-K0 anti-K'_10 pi- PHSP;
+0.000002  K0 anti-K0 K'_1- pi0   PHSP;
+0.000001  K0 anti-K0 K'_1- eta   PHSP;
+0.000001  K0 anti-K0 K'_1- eta'  PHSP;
+0.000002  anti-K0 anti-K0 K'_10 pi-        PHSP;
+0.000002  K0 K0 K'_1- pi0        PHSP;
+0.000001  K0 K0 K'_1- eta        PHSP;
+0.000001  K0 K0 K'_1- eta'       PHSP;
+
+0.000010  K- K'_10 anti-K0 pi0   PHSP;
+0.000005  K- K'_10 anti-K0 eta   PHSP;
+0.000002  K- K'_10 anti-K0 eta'  PHSP;
+0.000010  K- K0 anti-K'_10 pi0   PHSP;
+0.000005  K- K0 anti-K'_10 eta   PHSP;
+0.000002  K- K0 anti-K'_10 eta'  PHSP;
+#--- 5-body phi-phi-K, phi-phi-pi, phi-3K, phi-K-K-pi, 5K, 4K-pi ------------
+0.0000049 phi phi K-             PHSP;
+0.000001  phi phi pi-            PHSP;
+0.000001  phi K+ K- K-           PHSP;
+0.000001  phi K0 anti-K0 K-      PHSP;
+0.000001  phi K0 anti-K0 pi-     PHSP;
+0.000001  phi K0 K- pi0          PHSP;
+0.000001  K+ K- K+ K- pi-        PHSP;
+0.000001  K+ K- K0 anti-K0 pi-   PHSP;
+0.000001  K+ K- K- K0 pi0        PHSP;
+0.000001  K0 anti-K0 K0 K- pi0   PHSP;
+#--- 6-body phi-phi-K*, phi-3Kpi, phi-K-K-2pi ------------
+0.000003  phi phi K*-            PHSP;
+0.000001  phi K+ K- K*-          PHSP;
+0.000001  phi K0 anti-K0 K*-     PHSP;
+0.000001  phi K0 K- anti-K*0          PHSP; ## 2007 changed K*0 to anti-K*0
+0.000001  phi K+ K- K- pi0       PHSP;
+0.000001  phi K+ K- anti-K0 pi-  PHSP;
+0.000001  phi K0 anti-K0 K- pi0  PHSP;
+0.000001  phi K0 anti-K0 anti-K0 pi- PHSP;
+0.000001  phi K0 K- K- pi+       PHSP;
+#--- 6-body a1a1 ------------
+0.000050  a_1-  a_10             PHSP;
+#
+#   B -> cc= s          sum = 1.92%
+#
+0.001007  J/psi K-                         SVS;
+0.00141   J/psi K*-                        SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.000049  J/psi pi-                       SVS;
+0.00006   J/psi rho-                      SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.0002   J/psi anti-K0   pi-                    PHSP;
+0.0001   J/psi K-  pi0                    PHSP;
+#rl0.0007   J/psi K-  pi+  pi-               PHSP;
+#rl0.00035   J/psi K-  pi0  pi0               PHSP;
+#rl0.00035   J/psi anti-K0   pi-  pi0               PHSP;
+0.0018     J/psi  K_1-              SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0001   J/psi K'_1-                       SVV_HELAMP 0.5 0.0 1.0 0.0 0.5 0.0;
+0.0005   J/psi K_2*-                       PHSP;
+
+0.000052    J/psi  phi K-        PHSP;
+#
+0.000648  psi(2S) K-                        SVS;
+0.00067  psi(2S) K*-                       SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.0004   psi(2S) anti-K0   pi-                   PHSP;
+0.0002   psi(2S) K-  pi0                   PHSP;
+0.0019   psi(2S) K-  pi+  pi-              PHSP;
+0.0001   psi(2S) K-  pi0  pi0              PHSP;
+0.0001   psi(2S) anti-K0   pi-  pi0              PHSP;
+0.0004     psi(2S)  K_1-              PHSP;
+#
+0.00091   eta_c  K-                      PHSP;
+0.0012    K*-   eta_c                    SVS;
+0.0002    eta_c  anti-K0   pi-                 PHSP;
+0.0001    eta_c  K-  pi0                 PHSP;
+0.0002    eta_c  K-  pi+  pi-            PHSP;
+0.0001    eta_c  K-  pi0  pi0            PHSP;
+0.0001    eta_c  anti-K0   pi-  pi0            PHSP;
+#
+0.00034    eta_c(2S) K-                      PHSP;
+0.00048    K*- eta_c(2S)                     SVS;
+0.00008    eta_c(2S)   anti-K0   pi-                 PHSP;
+0.00005    eta_c(2S)   K-  pi0                 PHSP;
+0.00008    eta_c(2S)   K-  pi+  pi-            PHSP;
+0.00005    eta_c(2S)   K-  pi0  pi0            PHSP;
+0.00005    eta_c(2S)   anti-K0   pi-  pi0            PHSP;
+#
+0.00014   chi_c0  K-                      PHSP;
+0.0004    K*-   chi_c0                    SVS;
+0.0002    chi_c0  anti-K0   pi-                 PHSP;
+0.0001    chi_c0  K-  pi0                 PHSP;
+0.0002    chi_c0  K-  pi+  pi-            PHSP;
+0.0001    chi_c0  K-  pi0  pi0            PHSP;
+0.0001    chi_c0  anti-K0   pi-  pi0            PHSP;
+#
+0.00049    chi_c1  K-                      SVS;
+0.00036   chi_c1  K*-                     SVV_HELAMP PKHplus PKphHplus PKHzero PKphHzero PKHminus PKphHminus;
+0.0004    chi_c1  anti-K0   pi-                 PHSP;
+0.0002    chi_c1  K-  pi0                 PHSP;
+0.0004    chi_c1  K-  pi+  pi-            PHSP;
+0.0002    chi_c1  K-  pi0  pi0            PHSP;
+0.0002    chi_c1  anti-K0   pi-  pi0            PHSP;
+#
+0.00002    chi_c2  K-                      STS;
+0.00001    chi_c2  K*-                     PHSP;
+0.0002    chi_c2  anti-K0   pi-                 PHSP;
+0.0001    chi_c2  K-  pi0                 PHSP;
+0.0002    chi_c2  K-  pi+  pi-            PHSP;
+0.0001    chi_c2  K-  pi0  pi0            PHSP;
+0.0001    chi_c2  anti-K0   pi-  pi0            PHSP;
+#
+0.00049   psi(3770) K-                    SVS;   
+0.0005    psi(3770) K*-                   PHSP;
+0.0003    psi(3770) anti-K0   pi-         PHSP;
+0.0002    psi(3770) K-  pi0               PHSP;
+0.0002    psi(3770) K-  pi+  pi-          PHSP;
+0.0001    psi(3770) K-  pi0  pi0          PHSP;
+0.0001    psi(3770) anti-K0   pi-  pi0    PHSP;
+0.0003    psi(3770) K_1-                  PHSP;
+#
+#   b -> c (sc=) -> D Ds X    Sum = 10%
+#
+0.0100    D0   D_s-                      PHSP;
+0.0082    D*0  D_s-                      SVS;
+0.0076    D_s*- D0                       SVS;
+0.0171    D_s*- D*0                  SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0006   D'_10  D_s-                     SVS;
+0.0012   D'_10  D_s*-                    SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0012   D_10  D_s-                      SVS;
+0.0024   D_10  D_s*-                     SVV_HELAMP 0.48 0.0 0.734 0.0 0.48 0.0;
+0.0042   D_2*0  D_s-                      STS;
+0.0040   D_2*0  D_s*-                    PHSP;
+#
+#
+# SVV_HELAMP from factorization, recommendation
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/168.html:
+#
+0.0036    D_s-  D+  pi-                  PHSP;
+0.0018    D_s-  D0  pi0                  PHSP;
+0.0037    D_s*- D+  pi-                  PHSP;
+0.0018    D_s*- D0  pi0                  PHSP;
+0.0033    D_s-  D+  pi-  pi0             PHSP;
+0.0033    D_s-  D0  pi-  pi+             PHSP;
+0.0008    D_s-  D0  pi0  pi0             PHSP;
+0.0033    D_s*- D+  pi-  pi0             PHSP;
+0.0033    D_s*- D0  pi-  pi+             PHSP;
+0.0008    D_s*- D0  pi0  pi0             PHSP;
+#
+#  b -> c (sc=) -> D D= K X    Sum = 8%
+# update: Ref. [B1]:
+#
+0.0017   D0  D-  anti-K0                    PHSP; 
+0.0052   D0  D*-  anti-K0                   PHSP; 
+0.0031   D*0  D-  anti-K0                   PHSP; 
+0.0078   D*0  D*-  anti-K0                  PHSP;
+# External+internal W-emission amplitude
+0.00137  D0  anti-D0  K-                    PHSP; 
+0.0018   D*0  anti-D0  K-                   PHSP;
+0.0047   D0  anti-D*0  K-                   PHSP; 
+0.0053   D*0  anti-D*0  K-                  PHSP; 
+# Internal W-emission amplitude (color suppressed modes)
+0.0003    D-  D+  K-                        PHSP;
+0.0005    D*-  D+  K-                       PHSP;
+0.0015    D-  D*+  K-                       PHSP;
+0.0015    D*-  D*+  K-                      PHSP;
+#
+0.0025   D0  D-  anti-K*0                   PHSP;
+0.0025   D*0  D-  anti-K*0                  PHSP;
+0.0025   D0  D*-  anti-K*0                  PHSP;
+0.0050   D*0  D*-  anti-K*0                 PHSP;
+#
+0.0025   D0  anti-D0  K*-                   PHSP;
+0.0025   D*0  anti-D0  K*-                  PHSP;
+0.0025   D0  anti-D*0  K*-                  PHSP;
+0.0050   D*0  anti-D*0  K*-                 PHSP;
+#
+0.0005   D-  D+  K*-                        PHSP;
+0.0005   D*-  D+  K*-                       PHSP;
+0.0005   D-  D*+  K*-                       PHSP;
+0.0010   D*-  D*+  K*-                      PHSP;
+#
+# B->D(*)D(*). See Ref [B1]:
+#
+0.00042     D-   D0               PHSP;
+0.00039     D*-  D0               SVS;
+0.00063     D*0  D-               SVS;
+0.00081     D*0  D*-              SVV_HELAMP 0.47 0.0 0.96 0.0 0.56 0.0;
+
+#
+#  B -> D(*) X  Exclusive Modes
+#
+0.00519  D*0 pi-                        SVS;
+0.00484  D0   pi-                       PHSP;
+0.0134   rho-   D0                      SVS;
+# D* rho HELAMP parameters taken from ICHEP 98-852.
+0.0098   D*0  rho-               SVV_HELAMP 0.228 0.95 0.932 0.0 0.283 1.13;
+#
+0.0005   D0  pi0  pi-                   PHSP;
+0.0005   D*0 pi0  pi-                   PHSP;
+0.00102   D+  pi-  pi-                   PHSP;
+0.00125   D*+ pi-  pi-                   PHSP;
+#
+# D a1 updated Ref. [B1]:
+0.00887  a_1- D0                         SVS;
+0.0021   D0  rho0 pi-                   PHSP;
+0.0050   D0  pi+  pi-  pi-              PHSP;
+#
+# SVV_HELAMP from factorization, recommendation
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/168.html:
+# updated Ref. [B1]:
+#
+0.01900  D*0 a_1-                      SVV_HELAMP 0.200 0.0 0.866 0.0 0.458 0.0;
+#
+0.0005   D*0 rho0 pi-                   PHSP;
+0.00045  D*0 pi+  pi-  pi-              PHSP;
+#
+0.0020   D+  rho-  pi-                  PHSP;
+0.0020   D+  pi0   pi-  pi-             PHSP;
+0.0020   D*+ rho-  pi-                  PHSP;
+0.0150   D*+ pi0   pi-  pi-             PHSP;
+0.0005   D*0 rho-  pi0                  PHSP;
+0.0005   D*0 pi-   pi0  pi0             PHSP;
+#
+# B->D** pi and D** rho, recommendation
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html
+# further update October 26, 2004 Lange
+#
+0.0013   D_10 pi-                       SVS;
+0.00075   D'_10 pi-                      SVS;
+0.00086   D_2*0 pi-                      STS;
+0.00091 D_0*0 pi- PHSP;
+0.0007   D_10 rho-                      PHSP;
+0.0022   D'_10 rho-                    PHSP;
+0.0038   D_2*0 rho-                    PHSP;
+
+#
+#
+# B->DK, recommendation 
+# http://babar-hn.slac.stanford.edu:5090/HyperNews/get/event_gen/151.html:
+# update: Ref [B1]:
+#
+0.000402 D0  K-                                PHSP;
+0.000416 D*0 K-                                SVS;
+0.00053  K*- D0                        SVS;
+0.00081  D*0 K*-                       SVV_HELAMP 0.228 0.0  0.932 0.0 0.283 0.0;
+#October 26, 2004 Lange
+0.00000054 D- pi0  PHSP;
+0.00000054 D*- pi0  SVS;
+0.0000004 D- anti-K0 PHSP;#2007 changed BR  
+0.0000006 D*- anti-K0 SVS;#2007 changed BR 
+#
+#October 26, 2004 - Lange
+0.00073 D0 D_s0*- PHSP;
+0.0009 D*0 D_s0*- SVS;
+0.0031 D_s1- D0 SVS;
+0.012 D*0 D_s1- SVV_HELAMP 0.4904 0.0 0.7204 0.0 0.4904 0.0;
+#
+0.020        cd_1       anti-uu_1         PYTHIA 22;
+0.012        cd_1       anti-uu_1         PYTHIA 23;
+0.003        cs_1       anti-uu_1         PYTHIA 23;
+0.005        cs_0       anti-cu_0         PYTHIA 22;
+0.010        cs_1       anti-cu_1         PYTHIA 22;
+0.001        cs_0       anti-cu_0         PYTHIA 23;
+0.002        cs_1       anti-cu_1         PYTHIA 23;
+#
+0.214266        anti-u       d  c          anti-u       PYTHIA         48;
+0.04300         anti-u       d  c          anti-u       PYTHIA         13;
+0.02250         anti-u       s  c          anti-u       PYTHIA         13;
+#lange - try to crank up the psi production....
+0.07400         anti-c       s  c          anti-u       PYTHIA         13;
+0.00400         anti-c       d  c          anti-u       PYTHIA         13;
+0.00300         anti-u       d  u          anti-u       PYTHIA         48;
+0.00400         anti-c       s  u          anti-u       PYTHIA         48;
+# JGS 11/5/02 This and similar a few lines above have been divided by two
+# to solve a double-counting problem for this channel 
+0.00205         anti-u     u  d       anti-u       PYTHIA         48;
+0.00007         anti-d     d  d       anti-u       PYTHIA         48;
+0.00009         anti-s     s  d       anti-u       PYTHIA         48;
+0.00220         anti-u     u  s       anti-u       PYTHIA         48;
+0.00180         anti-d     d  s       anti-u       PYTHIA         48;
+0.00150         anti-s     s  s       anti-u       PYTHIA         48;
+0.00500                    s  g     anti-u       PYTHIA         33;
+#
+####
+
+Enddecay
+#
+
+Decay B+
+#
+#                  b -> c semileptonic
+0.0650   anti-D*0    e+   nu_e        PHOTOS  HQET 0.77 1.33 0.92;
+0.0215   anti-D0    e+   nu_e        PHOTOS  ISGW2;
+0.0056   anti-D_10   e+   nu_e        PHOTOS  ISGW2;
+0.0049   anti-D_0*0   e+   nu_e        PHOTOS  ISGW2;
+0.009   anti-D'_10   e+   nu_e        PHOTOS  ISGW2;
+0.003   anti-D_2*0   e+   nu_e        PHOTOS  ISGW2;
+0.0063   D*-  pi+   e+   nu_e    PHOTOS  GOITY_ROBERTS;
+0.0003   anti-D*0  pi0   e+   nu_e    PHOTOS  GOITY_ROBERTS;
+0.0052   D-   pi+   e+   nu_e    PHOTOS  GOITY_ROBERTS;
+0.0010   anti-D0  pi0   e+   nu_e    PHOTOS  GOITY_ROBERTS;
+
+0.0650   anti-D*0    mu+  nu_mu       PHOTOS   HQET 0.77 1.33 0.92;
+0.0215   anti-D0    mu+  nu_mu        PHOTOS  ISGW2;
+0.0056   anti-D_10   mu+  nu_mu       PHOTOS   ISGW2;
+0.0049   anti-D_0*0   mu+  nu_mu      PHOTOS    ISGW2;
+0.009   anti-D'_10   mu+  nu_mu      PHOTOS    ISGW2;
+0.003   anti-D_2*0   mu+  nu_mu      PHOTOS    ISGW2;
+0.0063   D*-  pi+   mu+   nu_mu       PHOTOS    GOITY_ROBERTS;
+0.0003   anti-D*0  pi0   mu+   nu_mu  PHOTOS         GOITY_ROBERTS;
+0.0052   D-   pi+   mu+   nu_mu       PHOTOS    GOITY_ROBERTS;
+0.0010   anti-D0  pi0   mu+   nu_mu   PHOTOS        GOITY_ROBERTS;
+#
+#                   b -> c tau nu
+#
+0.0208   anti-D*0    tau+ nu_tau                ISGW2; 
+0.0088   anti-D0    tau+ nu_tau                ISGW2;
+0.0013   anti-D_10   tau+ nu_tau                ISGW2;
+0.0013   anti-D_0*0   tau+ nu_tau                ISGW2;
+0.0020   anti-D'_10   tau+ nu_tau                ISGW2;
+0.0020   anti-D_2*0   tau+ nu_tau                ISGW2;
+#
+#                   b -> u l nu
+#
+0.000077   pi0    e+   nu_e        PHOTOS  ISGW2;
+0.000080   eta    e+   nu_e        PHOTOS  ISGW2;
+0.000128   rho0   e+   nu_e        PHOTOS  ISGW2;
+0.000130   omega   e+   nu_e       PHOTOS  ISGW2;
+0.000084   eta'   e+   nu_e        PHOTOS  ISGW2;
+0.000000   anti-D(2S)0 e+   nu_e   PHOTOS  ISGW2;
+0.000000   anti-D*(2S)0 e+   nu_e  PHOTOS  ISGW2;
+0.001948   Xu0   e+  nu_e          PHOTOS  VubHybridBp;
+#
+0.000077   pi0    mu+  nu_mu       PHOTOS  ISGW2;
+0.000080   eta    mu+  nu_mu       PHOTOS  ISGW2;
+0.000128   rho0   mu+  nu_mu       PHOTOS  ISGW2;
+0.000130   omega   mu+  nu_mu      PHOTOS  ISGW2;
+0.000084   eta'   mu+  nu_mu       PHOTOS  ISGW2;
+0.000000   anti-D(2S)0 mu+  nu_mu  PHOTOS  ISGW2;
+0.000000   anti-D*(2S)0 mu+  nu_mu PHOTOS  ISGW2;
+0.001948   Xu0   mu+  nu_mu        PHOTOS  VubHybridBp;
+#
+0.000030   pi0    tau+ nu_tau              ISGW2;
+0.000012   eta    tau+ nu_tau              ISGW2;
+0.000042   rho0   tau+ nu_tau              ISGW2;
+0.000042   omega   tau+ nu_tau              ISGW2;
+0.000018   eta'   tau+ nu_tau              ISGW2;
+0.000046   a_10    tau+ nu_tau              ISGW2;
+0.000027   b_10    tau+ nu_tau              ISGW2;
+0.000004   a_00    tau+ nu_tau              ISGW2;
+0.000002   f_0     tau+ nu_tau              ISGW2;
+0.000002   f'_0    tau+ nu_tau              ISGW2;
+0.000023   f_1     tau+ nu_tau              ISGW2;
+0.000023   f'_1    tau+ nu_tau              ISGW2;
+0.000013   h_1     tau+ nu_tau              ISGW2;
+0.000013   h'_1    tau+ nu_tau              ISGW2;
+0.000020   f_2     tau+ nu_tau              ISGW2;
+0.000020   f'_2    tau+ nu_tau              ISGW2;
+0.000000   anti-D(2S)0 tau+ nu_tau              ISGW2;
+0.000000   anti-D*(2S)0 tau+ nu_tau              ISGW2;
+#
+#
+#   b -> l nu
+#
+0.000180   tau+ nu_tau                     SLN;
+0.000000   mu+  nu_mu          PHOTOS     SLN;
+0.000000   e+   nu_e           PHOTOS    SLN;
+
+#
+# b->u hadronic 
+# Ref. [B1]:
+#
+0.000016   D_s+  pi0   PHSP;
+0.000020   D_s*+  pi0  SVS;
+#
+0.000028   rho0  D_s+     SVS;
+0.000028   D_s*+  rho0    SVV_HELAMP 1.0 0.0 1.0 0.0 1.0 0.0;
+
+#
+# b-> s gamma
+#
+0.0000403   K*+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000135  K_1+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000065  K'_1+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000128  K'*+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000166  K_2*+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+#0.0000017  K''*+  gamma                   HELAMP 1.0 0.0 1.0 0.0;
+0.0003137   Xsu    gamma  BTOXSGAMMA 2;
+#
+0.00000049  K+   e+    e-    PHOTOS BTOSLLBALL;
+0.000000700 K*+  e+    e-    PHOTOS BTOSLLBALL;
+0.000004994 Xsu  e+    e-    PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.00000039  K+   mu+   mu-   PHOTOS BTOSLLBALL;
+0.000000700 K*+  mu+   mu-   PHOTOS BTOSLLBALL;
+0.000002506 Xsu  mu+   mu-   PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+0.000000130 K+   tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000190 K*+  tau+  tau-  PHOTOS BTOSLLBALL;
+0.000000200 Xsu  tau+  tau-  PHOTOS BTOXSLL 4.8 0.2 0.0 0.41;
+#
+#------------------------------------------------------------------------------
+# 2- and 3-body modes revised Feb.2005  Markus Cristinziani, SLAC
+#vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv
+#
+# Exclusive hadronic b->u
+#
+0.000026   pi+   pi0   pi0                   CB3PI-P00 alpha;
+0.0000121  pi+   pi+   pi-                   CB3PI-MPP alpha;
+####
+0.0000049 eta pi+                PHSP;
+0.0000026 eta K+                 PHSP;
+0.0000193 K*+ eta                SVS;
+0.0000084 rho+ eta               SVS;
+0.0000026 eta' pi+               PHSP;
+0.0000697 eta' K+                PHSP;
+0.0000049 K*+ eta'               SVS;
+0.0000087 rho+ eta'              SVS;
+0.0000067 omega pi+              SVS;
+0.0000068 omega K+               SVS;
+0.000001  omega K*+              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000106 omega rho+             SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.0000001 phi pi+                SVS;
+0.0000083 phi K+                 SVS;
+0.000023  K0 pi+                 PHSP;
+0.00000128 anti-K0 K+             PHSP;
+0.0000055 pi+ pi0                PHSP;
+0.0000121 K+ pi0                 PHSP;
+#
+# 3-body John Back (jback@slac.stanford.edu) - Oct 15, 2002
+# JGS intersperses modes with pi0->eta,eta'
+# B+ modes
+#
+# rho0 3-body modes
+### AJR zeroed rho0pi- in favor of the CB3PI-P00 decay above
+0.0000000    rho0 pi+               SVS;
+0.0000042    rho0 K+                SVS;
+#
+# rho+ 3-body modes
+### AJR zeroed rho-pi0 in favor of the CB3PI-MPP decay above
+0.0000000    rho+ pi0               SVS;
+0.000001     rho+ K0                SVS;
+#
+# rho(1450) 3-body modes
+0.0000022    rho(2S)0 pi+           SVS;
+0.000001     rho(2S)0 K+            SVS;
+#
+# f0(980) 3-body modes
+0.000001     f_0 pi+                PHSP;
+0.0000177    f_0 K+                 PHSP;
+0.000001     a_0+ K0                PHSP;
+0.000001     a_00 K+                PHSP;
+0.000001     a_00 pi+               PHSP;
+0.000001     a_0+ pi0               PHSP;
+
+0.0000082 f_2 pi+                   PHSP;
+
+#
+# K*0(1430) 3-body modes 
+0.000047     K_0*0 pi+              PHSP;
+0.000002     anti-K_0*0 K+          PHSP;
+0.000047     K_0*+ pi0              PHSP;
+0.000002     K_0*+ anti-K0          PHSP;
+#
+# K*(892) 3-body modes
+0.0000109    K*0 pi+               SVS;        
+0.0000069    K*+ pi0                SVS;
+0.000003     K*+ anti-K0            SVS;
+
+#
+#
+# Non-resonant 3-body left-overs
+# Most modes are set at 1e-6
+#
+# pi+pi-pi+: high mass + f0(400-1200) = (1+1)e-6
+#0402270.000002     pi+ pi- pi+            PHSP;
+#
+# K+ pi- pi+:total-K*pi  
+0.0000100    K+ pi- pi+             PHSP;
+#
+# K+ K- pi+: just non resonant  
+0.000001     K+ K- pi+              PHSP;
+#
+# K+K-K+: high mass structure near 1500 + non-res: total - phiK = 
+0.0000296    K+ K- K+               PHSP;
+#
+# 
+# K+K+pi-: suppressed mode (1e-7)
+0.0000001    K+ K+ pi-              PHSP;
+# K-pi+pi+: suppressed mode (1e-7)
+0.0000001    K- pi+ pi+             PHSP;
+#
+#
+# K+ anti-K0 pi0
+0.000001     K+ anti-K0 pi0         PHSP;
+0.000001     K+ anti-K0 eta         PHSP;
+0.000001     K+ anti-K0 eta'        PHSP;
+#
+# pi+ K0 pi0
+0.000001     pi+ K0 pi0             PHSP;
+0.000001     pi+ K0 eta             PHSP;
+0.000001     pi+ K0 eta'            PHSP;
+#
+# pi+ pi0 pi0          
+#0402270.000001     pi+ pi0 pi0            PHSP;
+0.000001     pi+ pi0 eta            PHSP;
+0.000001     pi+ pi0 eta'           PHSP;
+0.000001     pi+ eta eta            PHSP;
+0.000001     pi+ eta eta'           PHSP;
+#
+# K+ pi0 pi0 
+0.000001     K+ pi0 pi0             PHSP;
+0.000001     K+ pi0 eta             PHSP;
+0.000001     K+ pi0 eta'            PHSP;
+0.000001     K+ eta eta             PHSP;
+0.000001     K+ eta eta'            PHSP;
+0.000001     K+ eta' eta'           PHSP;
+#
+#
+# K+ K0 anti-K0:  4*BF(K+KsKs) (4 is a guess)
+0.000046     K+ K0 anti-K0          PHSP;
+#
+# pi+ K0 anti-K0
+0.000001     pi+ K0 anti-K0         PHSP;
+#
+#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
+# 2- and 3-body modes revised Feb.2005   Markus Cristinziani, SLAC
+#---------------------------------------------------------------------------
+#
+#4-body modes from Andrei Gritsan
+#--- 4-body rho-rho, a1-pi, rho-pi-pi, pi-pi-pi-pi --------------------------
+0.000018  rho+  rho0             SVV_HELAMP 1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  rho0 pi+ pi0           PHSP;
+0.000005  rho0 pi+ eta           PHSP;
+0.000002  rho0 pi+ eta'          PHSP;
+0.000010  omega pi+ pi0          PHSP;
+0.000005  omega pi+ eta          PHSP;
+0.000002  omega pi+ eta'         PHSP;
+0.000010  rho+ pi+ pi-           PHSP;
+0.000010  rho+ pi0 pi0           PHSP;
+0.000005  rho+ pi0 eta           PHSP;
+0.000002  rho+ eta eta           PHSP;
+0.000002  rho+ pi0 eta'          PHSP;
+0.000010  pi- pi+ pi+ pi0        PHSP;
+0.000005  pi- pi+ pi+ eta        PHSP;
+0.000002  pi- pi+ pi+ eta'       PHSP;
+0.000010  pi+ pi0 pi0 pi0        PHSP;
+0.000005  pi+ eta pi0 pi0        PHSP;
+0.000002  pi+ eta eta pi0        PHSP;
+0.000002  pi+ eta' pi0 pi0       PHSP;
+0.000001  pi+ eta' eta pi0       PHSP;
+0.0000204 a_10  pi+              SVS;
+0.0000264 a_1+  pi0              SVS;
+0.0000067 b_10  pi+              SVS;
+0.000001  b_1+  pi0              SVS;
+0.0000020  rho+ f_0               SVS; 
+0.000010  rho+ a_00              SVS;
+0.000002  rho0 a_0+              SVS;
+0.000010  f_0 pi+ pi0            PHSP;
+0.000001  a_00 pi+ pi0           PHSP;
+0.000001  a_0+ pi0 pi0           PHSP;
+0.000001  a_0- pi+ pi+           PHSP;
+#--- 4-body rho-K*, a1-K, rho-pi-K, K*-pi-pi, pi-pi-pi-K -------------------
+0.0000092 rho+  K*0              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000005  rho0  K*+              SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  pi+ pi0 K*0            PHSP;
+0.000005  pi+ eta K*0            PHSP;
+0.000067  pi- pi+ K*+            PHSP;
+0.000010  pi0 pi0 K*+            PHSP;
+0.000005  pi0 eta K*+            PHSP;
+0.000002  eta eta K*+            PHSP;
+0.000002  pi0 eta' K*+           PHSP;
+0.000001  eta eta' K*+           PHSP;
+0.000010  rho+  K0 pi0           PHSP;
+0.000005  rho+  K0 eta           PHSP;
+0.000002  rho+  K0 eta'          PHSP;
+0.000010  rho+  K+ pi-           PHSP;
+0.000010  rho-  K+ pi+           PHSP;
+0.000005  rho0  K0 pi+           PHSP;
+0.000005  rho0  K+ pi0           PHSP;
+0.000002  rho0  K+ eta           PHSP;
+0.000001  rho0  K+ eta'          PHSP;
+0.000005  omega K0 pi+           PHSP;
+0.000005  omega K+ pi0           PHSP;
+0.000002  omega K+ eta           PHSP;
+0.000001  omega K+ eta'          PHSP;
+0.000010  pi+ pi+ pi- K0         PHSP;
+0.000010  pi+ pi0 pi0 K0         PHSP;
+0.000005  pi+ eta pi0 K0         PHSP;
+0.000002  pi+ eta eta K0         PHSP;
+0.000010  pi+ pi- pi0 K+         PHSP;
+0.000005  pi+ pi- eta K+         PHSP;
+0.000010  pi0 pi0 pi0 K+         PHSP;
+0.000005  pi0 pi0 eta K+         PHSP;
+0.000002  pi0 eta eta K+         PHSP;
+0.000002  pi0 pi0 eta' K+        PHSP;
+0.000001  pi0 eta eta' K+        PHSP;
+0.000010  rho+ K_0*0             PHSP;
+0.000010  rho0 K_0*+             PHSP;
+0.000010  pi+ pi0 K_0*0          PHSP;
+0.000010  pi0 pi0 K_0*+          PHSP;
+0.0000052 K*+ f_0                SVS; 
+0.000010  a_10 K+                SVS;
+0.0000349 a_1+ K0                SVS;
+0.0000091 b_10 K+                SVS;
+0.000013  b_1+ K0                SVS;
+0.000005  K*+ a_00               SVS;
+0.000005  K*0 a_0+               SVS;
+#--- 4-body K*-K*, rho-K-K, K*-K-pi, phi-pi-pi, pi-pi-K-K ------------
+0.000001  anti-K*0 K*+           SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi   rho+             SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000001  phi   pi+ pi0          PHSP;
+0.0000005 phi   pi+ eta          PHSP;
+0.0000002 phi   pi+ eta'         PHSP;
+0.000001  K+ K- rho+             PHSP;
+0.000001  K0 anti-K0 rho+        PHSP;
+0.000001  anti-K0 K+ rho0        PHSP;
+0.000001  anti-K0 K+ omega       PHSP;
+0.000001  K+ K- pi+ pi0          PHSP;
+0.0000005 K+ K- pi+ eta          PHSP;
+0.0000002 K+ K- pi+ eta'         PHSP;
+0.000001  K0 anti-K0 pi+ pi0     PHSP;
+0.0000005 K0 anti-K0 pi+ eta     PHSP;
+0.0000002 K0 anti-K0 pi+ eta'    PHSP;
+0.000001  K0 K- pi+ pi+          PHSP;# 2007 added the CC to B- 
+0.000001  anti-K0 K+ pi+ pi-     PHSP;
+0.00000025 anti-K0 K+ pi0 pi0     PHSP;
+0.0000005 anti-K0 K+ pi0 eta     PHSP;
+0.00000025 anti-K0 K+ eta eta          PHSP;#added 2007 
+0.0000002 anti-K0 K+ pi0 eta'    PHSP;
+0.0000001 anti-K0 K+ eta eta'    PHSP;
+0.000001  K*+ K- pi+             PHSP;
+0.000001  K*- K+ pi+             PHSP;
+0.000001  K*0 anti-K0 pi+        PHSP;
+0.000001  K0 anti-K*0 pi+        PHSP;
+0.000001  K+ anti-K*0 pi0        PHSP;
+0.0000005 K+ anti-K*0 eta        PHSP;
+0.0000002 K+ anti-K*0 eta'       PHSP;
+0.000001  K*+ anti-K0 pi0        PHSP;
+0.0000005 K*+ anti-K0 eta        PHSP;
+0.0000002 K*+ anti-K0 eta'       PHSP;
+#--- 4-body phi-K*, phi-K-pi, K-K-K*, pi-K-K-K -----------------------
+0.0000096 phi K*+                SVV_HELAMP  1.0 0.0 1.7 0.0 1.0 0.0;
+0.000010  phi K0 pi+             PHSP;
+0.000010  phi K+ pi0             PHSP;
+0.000005  phi K+ eta             PHSP;
+0.000002  phi K+ eta'            PHSP;
+0.000031  K- K+ K*+              PHSP; #total -phiK*
+0.000010  K+ K*- K+    &n