Decay/FormFactors/BaryonSimpleFormFactor.h

0001 // -*- C++ -*-
0002 #ifndef HERWIG_BaryonSimpleFormFactor_H
0003 #define HERWIG_BaryonSimpleFormFactor_H
0004 //
0005 // This is the declaration of the BaryonSimpleFormFactor class.
0006 //
0007 #include "BaryonFormFactor.h"
0008
0009 namespace Herwig {
0010 using namespace ThePEG;
0011
0012 /** \ingroup Decay
0013  *
0014  *  The BaryonSimpleFormFactor class is a simple model for the form-factors
0015  *  for the semi-leptonic decay of the light (i.e. uds) baryons. The form-factors
0016  *  are assumed to be constant and are taken from the quark model results
0017  *  of PRD25, 206 (1982).
0018  *
0019  * @ see BaryonFormFactor
0020  */
0021 class BaryonSimpleFormFactor: public BaryonFormFactor {
0022
0023 public:
0024
0025   /**
0026    * Default constructor
0027    */
0028   BaryonSimpleFormFactor();
0029
0030   /** @name Functions used by the persistent I/O system. */
0031   //@{
0032   /**
0033    * Function used to write out object persistently.
0034    * @param os the persistent output stream written to.
0035    */
0036   void persistentOutput(PersistentOStream & os) const;
0037
0038   /**
0039    * Function used to read in object persistently.
0040    * @param is the persistent input stream read from.
0041    * @param version the version number of the object when written.
0042    */
0043   void persistentInput(PersistentIStream & is, int version);
0044   //@}
0045
0046   /**
0047    * Standard Init function used to initialize the interfaces.
0048    */
0049   static void Init();
0050
0051 public:
0052
0053   /** @name Form Factors */
0054   //@{
0055   /**
0056    * The form factor for the weak decay of a spin \f$\frac12\f$ baryon to a
0057    * spin \f$\frac12\f$ baryon.
0058    * @param q2 The scale \f$q^2\f$.
0059    * @param iloc The location in the form factor list.
0060    * @param id0 The PDG code of the incoming baryon.
0061    * @param id1 The PDG code of the outgoing baryon.
0062    * @param m0 The mass of the incoming baryon.
0063    * @param m1 The mass of the outgoing baryon.
0064    * @param f1v The form factor \f$F^V_1\f$.
0065    * @param f2v The form factor \f$F^V_2\f$.
0066    * @param f3v The form factor \f$F^V_3\f$.
0067    * @param f1a The form factor \f$F^A_1\f$.
0068    * @param f2a The form factor \f$F^A_2\f$.
0069    * @param f3a The form factor \f$F^A_3\f$.
0070    * @param virt Whether the q2 is space or timelike
0071    */
0072   virtual void SpinHalfSpinHalfFormFactor(Energy2 q2,int iloc, int id0, int id1,
0073                       Energy m0, Energy m1,
0074                       Complex & f1v,Complex & f2v,Complex & f3v,
0075                       Complex & f1a,Complex & f2a,Complex & f3a,
0076                       FlavourInfo flavour,
0077                       Virtuality virt=SpaceLike);
0078   //@}
0079
0080   /**
0081    * Output the setup information for the particle database
0082    * @param os The stream to output the information to
0083    * @param header Whether or not to output the information for MySQL
0084    * @param create Whether or not to add a statement creating the object
0085    */
0086   virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
0087
0088 protected:
0089
0090   /** @name Clone Methods. */
0091   //@{
0092   /**
0093    * Make a simple clone of this object.
0094    * @return a pointer to the new object.
0095    */
0096   virtual IBPtr clone() const {return new_ptr(*this);}
0097
0098   /** Make a clone of this object, possibly modifying the cloned object
0099    * to make it sane.
0100    * @return a pointer to the new object.
0101    */
0102   virtual IBPtr fullclone() const {return new_ptr(*this);}
0103   //@}
0104
0105 protected:
0106
0107
0108   /** @name Standard Interfaced functions. */
0109   //@{
0110   /**
0111    * Initialize this object after the setup phase before saving and
0112    * EventGenerator to disk.
0113    * @throws InitException if object could not be initialized properly.
0114    */
0115   virtual void doinit();
0116   //@}
0117
0118 private:
0119
0120   /**
0121    * Private and non-existent assignment operator.
0122    */
0123   BaryonSimpleFormFactor & operator=(const BaryonSimpleFormFactor &) = delete;
0124
0125 private:
0126
0127   /**
0128    * The axial vector coupling in \f$\beta\f$ decay.
0129    */
0130   double _gA;
0131
0132   /**
0133    *  The \f$\frac{D}{D+F}\f$ ratio.
0134    */
0135   double _alphaD;
0136
0137   /**
0138    *  The \f$\eta_V\f$ \f$SU(3)\f$ breaking parameter.
0139    */
0140   double _etaV;
0141
0142   /**
0143    *  The \f$\eta_V\f$ \f$SU(3)\f$ breaking parameter.
0144    */
0145   double _etaA;
0146
0147   /**
0148    *  The \f$\rho_E\f$  \f$SU(3)\f$ breaking parameter for the electic dipole moment.
0149    */
0150   double _rhoE;
0151
0152   /**
0153    *  The \f$\rho_M\f$  \f$SU(3)\f$ breaking parameter for the electic dipole moment.
0154    */
0155   double _rhoM;
0156
0157   /**
0158    *  The calculated constants for the \f$f_1\f$ form factor.
0159    */
0160   vector<double> _f1;
0161
0162   /**
0163    *  The calculated constants for the \f$f_2\f$ form factor.
0164    */
0165   vector<double> _f2;
0166
0167   /**
0168    *  The calculated constants for the \f$g_1\f$ form factor.
0169    */
0170   vector<double> _g1;
0171
0172   /**
0173    *  The calculated constants for the \f$g_2\f$ form factor.
0174    */
0175   vector<double> _g2;
0176
0177 };
0178
0179 }
0180
0181 #endif /* HERWIG_BaryonSimpleFormFactor_H */