diff -ru 9.1/EvtGenBase/EvtDecayAmp.hh 9.1/EvtGenBase/EvtDecayAmp.hh
--- 9.1/EvtGenBase/EvtDecayAmp.hh	2003-10-02 18:25:52.000000000 +0100
+++ 9.1/EvtGenBase/EvtDecayAmp.hh	2009-01-09 09:15:16.000000000 +0000
@@ -28,7 +28,7 @@
 
 public:
 
-  void makeDecay(EvtParticle* p);
+  void makeDecay(EvtParticle* p,bool recursive=true);
   inline void setWeight(double weight) {_weight=weight;}
 
   /**
@@ -60,6 +60,13 @@
   void vertex(int *i1, const EvtComplex& amp)
   { _amp2.vertex(i1,amp);}
 
+  /**
+   *  Provide access to the amplitude
+   */
+  const EvtAmp & amplitude() const 
+  { return _amp2;}
+
+
   virtual ~EvtDecayAmp() {}
 
 protected:
diff -ru 9.1/EvtGenBase/EvtDecayBase.hh 9.1/EvtGenBase/EvtDecayBase.hh
--- 9.1/EvtGenBase/EvtDecayBase.hh	2003-10-02 18:25:52.000000000 +0100
+++ 9.1/EvtGenBase/EvtDecayBase.hh	2009-01-09 09:15:40.000000000 +0000
@@ -37,7 +37,7 @@
   //by any derived class
   virtual void getName(std::string& name)=0;
   virtual void decay(EvtParticle *p)=0;
-  virtual void makeDecay(EvtParticle *p)=0;
+  virtual void makeDecay(EvtParticle *p,bool recursive=true)=0;
   virtual EvtDecayBase* clone()=0;
 
 
diff -ru 9.1/EvtGenBase/EvtDecayIncoherent.hh 9.1/EvtGenBase/EvtDecayIncoherent.hh
--- 9.1/EvtGenBase/EvtDecayIncoherent.hh	2003-10-02 18:25:52.000000000 +0100
+++ 9.1/EvtGenBase/EvtDecayIncoherent.hh	2009-01-09 09:15:48.000000000 +0000
@@ -30,7 +30,7 @@
 
 public:
 
-  void makeDecay(EvtParticle* p);
+  void makeDecay(EvtParticle* p,bool recursive=true);
 
   virtual ~EvtDecayIncoherent() {}
 
diff -ru 9.1/EvtGenBase/EvtDecayProb.hh 9.1/EvtGenBase/EvtDecayProb.hh
--- 9.1/EvtGenBase/EvtDecayProb.hh	2003-10-02 18:25:52.000000000 +0100
+++ 9.1/EvtGenBase/EvtDecayProb.hh	2009-01-09 09:15:57.000000000 +0000
@@ -29,7 +29,7 @@
 
 public:
 
-  void makeDecay(EvtParticle* p);
+  void makeDecay(EvtParticle* p,bool recursive=true);
 
   void setProb(double prob) { _prob=prob;}
   double getProb() { return _prob;}
diff -ru 9.1/EvtGenBase/EvtDiracParticle.hh 9.1/EvtGenBase/EvtDiracParticle.hh
--- 9.1/EvtGenBase/EvtDiracParticle.hh	2003-10-02 18:25:52.000000000 +0100
+++ 9.1/EvtGenBase/EvtDiracParticle.hh	2009-01-09 09:16:43.000000000 +0000
@@ -36,6 +36,9 @@
   EvtDiracParticle();
   virtual ~EvtDiracParticle();
   void init(EvtId part_n,const EvtVector4R& p4);
+  void init(EvtId part_n,const EvtVector4R& p4,
+	    const EvtDiracSpinor &,const EvtDiracSpinor &,
+	    const EvtDiracSpinor &,const EvtDiracSpinor &);
   EvtDiracSpinor spParent(int) const;
   EvtDiracSpinor sp(int) const;   
   EvtSpinDensity rotateToHelicityBasis() const;
diff -ru 9.1/EvtGenBase/EvtParticle.hh 9.1/EvtGenBase/EvtParticle.hh
--- 9.1/EvtGenBase/EvtParticle.hh	2003-10-10 10:21:43.000000000 +0100
+++ 9.1/EvtGenBase/EvtParticle.hh	2009-01-09 09:17:38.000000000 +0000
@@ -33,6 +33,7 @@
 class EvtTensor4C;
 class EvtStdHep;
 class EvtSecondary;
+class EvtRaritaSchwinger;
 
 const int MAX_DAUG =100;
 const int MAX_LEVEL=10;
@@ -107,6 +108,18 @@
   * Returns tensor in the particles own restframe for a spin 2 particle.
   */
   virtual EvtTensor4C epsTensor(int i) const; 
+   
+  /**
+   * Returns Rarita-Schwinger spinor in the parents restframe for a 
+   * Rarita-Schwinger particle.
+   */
+  virtual EvtRaritaSchwinger spRSParent(int) const;
+  
+  /**
+   * Returns Rarita-Schwinger spinor in the particles own restframe for a 
+   * Rarita-Schwinger particle.
+   */
+  virtual EvtRaritaSchwinger spRS(int) const;
   
 
   /**
diff -ru 9.1/EvtGenBase/EvtPDL.hh 9.1/EvtGenBase/EvtPDL.hh
--- 9.1/EvtGenBase/EvtPDL.hh	2007-10-14 19:50:07.000000000 +0100
+++ 9.1/EvtGenBase/EvtPDL.hh	2009-01-09 09:27:12.000000000 +0000
@@ -69,6 +69,7 @@
   static void alias(EvtId num,const std::string& newname);
   static void aliasChgConj(EvtId a,EvtId abar);
   static int entries() ;
+  static EvtId getEntry(int i);
   static void reSetMass(EvtId i, double mass) ;
   static void reSetWidth(EvtId i, double width);
   static void reSetMassMin(EvtId i, double mass);
diff -ru 9.1/EvtGenBase/EvtRaritaSchwingerParticle.hh 9.1/EvtGenBase/EvtRaritaSchwingerParticle.hh
--- 9.1/EvtGenBase/EvtRaritaSchwingerParticle.hh	2003-10-02 18:25:55.000000000 +0100
+++ 9.1/EvtGenBase/EvtRaritaSchwingerParticle.hh	2009-01-09 09:32:54.000000000 +0000
@@ -35,6 +35,11 @@
   EvtRaritaSchwingerParticle();
   virtual ~EvtRaritaSchwingerParticle();
   void init(EvtId id,const EvtVector4R& p4);
+  void init(EvtId id,const EvtVector4R& p4,
+	    const EvtRaritaSchwinger &,const EvtRaritaSchwinger &,
+	    const EvtRaritaSchwinger &,const EvtRaritaSchwinger &,
+	    const EvtRaritaSchwinger &,const EvtRaritaSchwinger &,
+	    const EvtRaritaSchwinger &,const EvtRaritaSchwinger &);
   EvtRaritaSchwinger spRSParent(int) const;
   EvtRaritaSchwinger spRS(int) const;   
   EvtSpinDensity rotateToHelicityBasis() const;
diff -ru 9.1/EvtGenBase/EvtTensorParticle.hh 9.1/EvtGenBase/EvtTensorParticle.hh
--- 9.1/EvtGenBase/EvtTensorParticle.hh	2003-10-02 18:25:57.000000000 +0100
+++ 9.1/EvtGenBase/EvtTensorParticle.hh	2009-01-09 09:19:40.000000000 +0000
@@ -33,7 +33,10 @@
 
   void init(EvtId part_n,double e,double px,double py,double pz);
   void init(EvtId part_n,const EvtVector4R& p4);
-//Returns polarization tensors.
+  void init(EvtId part_n,const EvtVector4R& p4,
+	    const EvtTensor4C&,const EvtTensor4C&,const EvtTensor4C&,
+	    const EvtTensor4C&,const EvtTensor4C&);
+  //Returns polarization tensors.
   EvtTensor4C epsTensorParent(int i) const; 
   EvtTensor4C epsTensor(int i) const; 
 
diff -ru 9.1/EvtGenBase/EvtVectorParticle.hh 9.1/EvtGenBase/EvtVectorParticle.hh
--- 9.1/EvtGenBase/EvtVectorParticle.hh	2003-10-02 18:25:57.000000000 +0100
+++ 9.1/EvtGenBase/EvtVectorParticle.hh	2009-01-09 09:19:58.000000000 +0000
@@ -35,6 +35,8 @@
 
   void init(EvtId part_n,double e,double px,double py,double pz);
   void init(EvtId part_n,const EvtVector4R& p);
+  void init(EvtId part_n,const EvtVector4R& p,
+	    const EvtVector4C&,const EvtVector4C&,const EvtVector4C&);
   EvtVector4C epsParent(int i) const; 
   EvtVector4C eps(int i) const; 
   EvtSpinDensity rotateToHelicityBasis() const;
diff -ru 9.1/src/EvtDecayAmp.cpp 9.1/src/EvtDecayAmp.cpp
--- 9.1/src/EvtDecayAmp.cpp	2004-07-12 17:13:28.000000000 +0100
+++ 9.1/src/EvtDecayAmp.cpp	2009-01-09 09:22:40.000000000 +0000
@@ -30,7 +30,7 @@
 #include "EvtGenBase/EvtReport.hh"
 
 
-void EvtDecayAmp::makeDecay(EvtParticle* p){
+void EvtDecayAmp::makeDecay(EvtParticle* p,bool recursive){
 
   int ntimes=10000;
 
@@ -153,6 +153,7 @@
   }
 
   //report(INFO,"EvtGen") << "Found " << p->getNDaug() << " daughters\n";
+  if(recursive) {
   for(i=0;i<p->getNDaug();i++){
 
     rho.SetDim(_amp2.dstates[i]);
@@ -206,6 +207,7 @@
     report(ERROR,"EvtGen") << p->getSpinDensityBackward();
 
   }
+  }
 
   if (getPHOTOS() || EvtRadCorr::alwaysRadCorr()) {
     EvtRadCorr::doRadCorr(p);
diff -ru 9.1/src/EvtDecayIncoherent.cpp 9.1/src/EvtDecayIncoherent.cpp
--- 9.1/src/EvtDecayIncoherent.cpp	2004-07-12 17:13:28.000000000 +0100
+++ 9.1/src/EvtDecayIncoherent.cpp	2009-01-09 09:35:55.000000000 +0000
@@ -26,7 +26,7 @@
 #include "EvtGenBase/EvtPDL.hh"
 
 
-void EvtDecayIncoherent::makeDecay(EvtParticle* p){
+void EvtDecayIncoherent::makeDecay(EvtParticle* p, bool recursive){
 
   int i;
   //initialize this the hard way..
@@ -45,6 +45,8 @@
     EvtRadCorr::doRadCorr(p);
   }
 
+  if(!recursive) return;
+  
   //Now decay the daughters.
   for(i=0;i<p->getNDaug();i++){
     //Need to set the spin density of the daughters to be
diff -ru 9.1/src/EvtDecayProb.cpp 9.1/src/EvtDecayProb.cpp
--- 9.1/src/EvtDecayProb.cpp	2004-07-12 17:13:29.000000000 +0100
+++ 9.1/src/EvtDecayProb.cpp	2009-01-09 09:38:44.000000000 +0000
@@ -26,7 +26,7 @@
 #include "EvtGenBase/EvtPDL.hh"
 #include "EvtGenBase/EvtReport.hh"
 
-void EvtDecayProb::makeDecay(EvtParticle* p){
+void EvtDecayProb::makeDecay(EvtParticle* p, bool recursive){
 
   int ntimes=10000;
 
@@ -70,6 +70,7 @@
     EvtRadCorr::doRadCorr(p);
   }
 
+  if(!recursive) return;
 
   //Now decay the daughters.
   int i;
diff -ru 9.1/src/EvtDiracParticle.cpp 9.1/src/EvtDiracParticle.cpp
--- 9.1/src/EvtDiracParticle.cpp	2005-04-05 22:55:03.000000000 +0100
+++ 9.1/src/EvtDiracParticle.cpp	2009-01-09 09:23:40.000000000 +0000
@@ -79,6 +79,26 @@
   setLifetime();
 }
 
+void EvtDiracParticle::init(EvtId part_n,const EvtVector4R& p4,
+			    const EvtDiracSpinor & prod1, const EvtDiracSpinor & prod2,
+			    const EvtDiracSpinor & rest1, const EvtDiracSpinor & rest2){
+  
+  _validP4=true; 
+  setp(p4);
+  setpart_num(part_n);
+  
+  if (EvtPDL::getStdHep(part_n)==0){
+    report(ERROR,"EvtGen") << "Error in EvtDiracParticle::init, part_n="
+                           << part_n.getId()<<std::endl;
+    ::abort();
+  }
+  spinor1_rest=rest1;
+  spinor2_rest=rest2;
+  spinor1_parent=prod1;
+  spinor2_parent=prod2;
+  
+  setLifetime();
+}
 
 EvtDiracSpinor EvtDiracParticle::spParent(int i) const {
 
diff -ru 9.1/src/EvtParticle.cpp 9.1/src/EvtParticle.cpp
--- 9.1/src/EvtParticle.cpp	2004-07-12 17:13:31.000000000 +0100
+++ 9.1/src/EvtParticle.cpp	2009-01-09 09:39:09.000000000 +0000
@@ -39,6 +39,7 @@
 #include "EvtGenBase/EvtScalarParticle.hh"
 #include "EvtGenBase/EvtVectorParticle.hh"
 #include "EvtGenBase/EvtTensorParticle.hh"
+#include "EvtGenBase/EvtRaritaSchwingerParticle.hh"
 #include "EvtGenBase/EvtPhotonParticle.hh"
 #include "EvtGenBase/EvtNeutrinoParticle.hh"
 #include "EvtGenBase/EvtStringParticle.hh"
@@ -688,6 +689,32 @@
   return temp;
 }
 
+EvtRaritaSchwinger EvtParticle::spRSParent(int i) const {
+  EvtRaritaSchwinger tempD;
+  int temp;
+  temp = i;
+  printParticle();
+  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
+			 <<"th Rarita-Schwinger spinor."
+			 <<" I.e. you thought it was a"
+			 <<" RaritaSchwinger particle!" << std::endl;
+  ::abort();
+  return tempD;
+}
+
+EvtRaritaSchwinger EvtParticle::spRS(int i) const {
+  EvtRaritaSchwinger tempD;
+  int temp;
+  temp = i;
+  printParticle();
+  report(ERROR,"EvtGen") << "and you have asked for the:"<<i
+			 <<"th Rarita-Schwinger spinor."
+			 <<" I.e. you thought it was a"
+			 <<" RaritaSchwinger particle!" << std::endl;
+  ::abort();
+  return tempD;
+}
+
 EvtVector4R EvtParticle::getP4Restframe() {
 
   return EvtVector4R(mass(),0.0,0.0,0.0);
diff -ru 9.1/src/EvtPDL.cpp 9.1/src/EvtPDL.cpp
--- 9.1/src/EvtPDL.cpp	2007-10-14 19:50:07.000000000 +0100
+++ 9.1/src/EvtPDL.cpp	2009-01-09 10:27:23.000000000 +0000
@@ -379,6 +379,10 @@
   return _partlist.size();
 }
 
+EvtId EvtPDL::getEntry(int i) {
+  return _partlist[i].getId();
+}
+
 void EvtPDL::reSetMass(EvtId i, double mass) {
   _partlist[i.getId()].reSetMass(mass);
 }
diff -ru 9.1/src/EvtRaritaSchwingerParticle.cpp 9.1/src/EvtRaritaSchwingerParticle.cpp
--- 9.1/src/EvtRaritaSchwingerParticle.cpp	2004-07-12 17:13:32.000000000 +0100
+++ 9.1/src/EvtRaritaSchwingerParticle.cpp	2009-01-09 09:25:08.000000000 +0000
@@ -116,7 +116,37 @@
 
 }
 
-
+void EvtRaritaSchwingerParticle::init(EvtId id,const EvtVector4R& p4,
+				      const EvtRaritaSchwinger & prod1,
+				      const EvtRaritaSchwinger & prod2,
+				      const EvtRaritaSchwinger & prod3,
+				      const EvtRaritaSchwinger & prod4,
+				      const EvtRaritaSchwinger & rest1,
+				      const EvtRaritaSchwinger & rest2,
+				      const EvtRaritaSchwinger & rest3,
+				      const EvtRaritaSchwinger & rest4){
+
+  _validP4=true;
+  setp(p4);
+  setpart_num(id);
+
+  if (EvtPDL::getStdHep(id)==0){
+    report(ERROR,"EvtGen") << "Error in EvtRaritaSchwingerParticle::init, part_n="
+                           << id.getId()<<std::endl;
+    ::abort();
+  }
+  _spinorRest[0]=rest1;
+  _spinorRest[1]=rest2;
+  _spinorRest[2]=rest3;
+  _spinorRest[3]=rest4;
+    
+  _spinor[0]=prod1;
+  _spinor[1]=prod2;
+  _spinor[2]=prod3;
+  _spinor[3]=prod4;
+  
+  setLifetime();
+}
 
 EvtSpinDensity EvtRaritaSchwingerParticle::rotateToHelicityBasis() const{
   
diff -ru 9.1/src/EvtTensorParticle.cpp 9.1/src/EvtTensorParticle.cpp
--- 9.1/src/EvtTensorParticle.cpp	2004-07-12 17:13:34.000000000 +0100
+++ 9.1/src/EvtTensorParticle.cpp	2009-01-09 09:25:44.000000000 +0000
@@ -93,6 +93,25 @@
   
 }
 
+void EvtTensorParticle::init(EvtId part_n,const EvtVector4R& p4,
+			     const EvtTensor4C& epsin1, const EvtTensor4C& epsin2,
+			     const EvtTensor4C& epsin3, const EvtTensor4C& epsin4,
+			     const EvtTensor4C& epsin5){
+ 
+  _validP4=true;
+  setp(p4);
+  setpart_num(part_n);
+
+  eps1=epsin1;
+  eps2=epsin2;
+  eps3=epsin3;
+  eps4=epsin4;
+  eps5=epsin5;
+
+  setLifetime();
+  
+}
+
 EvtTensor4C EvtTensorParticle::epsTensorParent(int i) const {
 
   EvtTensor4C temp;
diff -ru 9.1/src/EvtVectorParticle.cpp 9.1/src/EvtVectorParticle.cpp
--- 9.1/src/EvtVectorParticle.cpp	2004-07-12 17:13:35.000000000 +0100
+++ 9.1/src/EvtVectorParticle.cpp	2009-01-09 09:26:35.000000000 +0000
@@ -56,6 +56,22 @@
   setLifetime();
 }
 
+void EvtVectorParticle::init(EvtId part_n,const EvtVector4R& p4,
+			     const EvtVector4C & epsin1,
+			     const EvtVector4C & epsin2,
+			     const EvtVector4C & epsin3){
+
+  _validP4=true;
+  setp(p4);
+  setpart_num(part_n);
+  
+  eps1=epsin1;
+  eps2=epsin2;
+  eps3=epsin3;
+
+  setLifetime();
+}
+
 EvtVector4C EvtVectorParticle::epsParent(int i) const  {
 
   EvtVector4C temp;