SemiLeptonicBaryonDecayer.h

// -*- C++ -*-
#ifndef HERWIG_SemiLeptonicBaryonDecayer_H
#define HERWIG_SemiLeptonicBaryonDecayer_H
//
// This is the declaration of the SemiLeptonicBaryonDecayer class.
//
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Decay/WeakCurrents/LeptonNeutrinoCurrent.h"
#include "Herwig/Decay/FormFactors/BaryonFormFactor.h"
#include "Herwig/Decay/PhaseSpaceMode.h"
#include "ThePEG/Helicity/LorentzRSSpinorBar.h"
 
namespace Herwig {
using namespace ThePEG;
 
class SemiLeptonicBaryonDecayer: public DecayIntegrator {
 
public:
 
  SemiLeptonicBaryonDecayer();
 
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const;
 
  virtual int modeNumber(bool & cc, tcPDPtr parent, 
                         const tPDVector & children) const;
 
  double me2(const int ichan,const Particle & part,
             const tPDVector & outgoing,
             const vector<Lorentz5Momentum> & momenta,
             MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
  
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
  
  static void Init();
  
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
protected:
 
  double halfHalf(const Particle & part,
                  const tPDVector & outgoing,
                  const vector<Lorentz5Momentum> & momenta,
                  MEOption meopt) const;
 
  double halfThreeHalf(const Particle & part,
                       const tPDVector & outgoing,
                       const vector<Lorentz5Momentum> & momenta,
                       MEOption meopt) const;
 
private:
 
  SemiLeptonicBaryonDecayer & operator=(const SemiLeptonicBaryonDecayer &) = delete;
 
private:
 
  LeptonNeutrinoCurrentPtr _current;
 
  BaryonFormFactorPtr _form;
 
  vector<double> _maxwgt;
 
  vector<int> _modemap; 
 
  mutable RhoDMatrix _rho;
 
  mutable vector<LorentzSpinor<SqrtEnergy> > _inHalf;
 
  mutable vector<LorentzSpinorBar<SqrtEnergy> > _inHalfBar;
 
  mutable vector<LorentzRSSpinor<SqrtEnergy> > _inThreeHalf;
 
  mutable vector<LorentzRSSpinorBar<SqrtEnergy> > _inThreeHalfBar;
 
  mutable vector<unsigned int> _constants;
 
  mutable vector<PDT::Spin> _ispin;
 
  mutable unsigned int _ibar;
};
 
}
 
 
#endif /* HERWIG_SemiLeptonicBaryonDecayer_H */