Decay/FormFactors/LightBaryonQuarkModelFormFactor.h

0001 // -*- C++ -*-
0002 #ifndef HERWIG_LightBaryonQuarkModelFormFactor_H
0003 #define HERWIG_LightBaryonQuarkModelFormFactor_H
0004 // This is the declaration of the LightBaryonQuarkModelFormFactor class.
0005
0006 #include "BaryonFormFactor.h"
0007
0008 namespace Herwig {
0009 using namespace ThePEG;
0010
0011 /** \ingroup Decay
0012  *
0013  *  The <code>LightBaryonQuarkModelFormFactor</code> class implements the
0014  *  quark model calculation of hep-ph/9409272 for the form-factors for
0015  *  the decay of baryons containing the light quarks. This can be adapted for
0016  *  other models.
0017  *
0018  * @see BaryonFormFactor
0019  *
0020  */
0021 class LightBaryonQuarkModelFormFactor: public BaryonFormFactor {
0022
0023 public:
0024
0025   /**
0026    * Default constructor
0027    */
0028   LightBaryonQuarkModelFormFactor();
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
0052 protected:
0053
0054   /** @name Clone Methods. */
0055   //@{
0056   /**
0057    * Make a simple clone of this object.
0058    * @return a pointer to the new object.
0059    */
0060   virtual IBPtr clone() const {return new_ptr(*this);}
0061
0062   /** Make a clone of this object, possibly modifying the cloned object
0063    * to make it sane.
0064    * @return a pointer to the new object.
0065    */
0066   virtual IBPtr fullclone() const {return new_ptr(*this);}
0067   //@}
0068
0069 public:
0070
0071   /** @name Form Factors */
0072   //@{
0073   /**
0074    * The form factor for the weak decay of a spin \f$\frac12\f$ baryon to a
0075    * spin \f$\frac12\f$ baryon.
0076    * @param q2 The scale \f$q^2\f$.
0077    * @param iloc The location in the form factor list.
0078    * @param id0 The PDG code of the incoming baryon.
0079    * @param id1 The PDG code of the outgoing baryon.
0080    * @param m0 The mass of the incoming baryon.
0081    * @param m1 The mass of the outgoing baryon.
0082    * @param f1v The form factor \f$F^V_1\f$.
0083    * @param f2v The form factor \f$F^V_2\f$.
0084    * @param f3v The form factor \f$F^V_3\f$.
0085    * @param f1a The form factor \f$F^A_1\f$.
0086    * @param f2a The form factor \f$F^A_2\f$.
0087    * @param f3a The form factor \f$F^A_3\f$.
0088    * @param virt Whether the q2 is space or timelike
0089    */
0090   virtual void SpinHalfSpinHalfFormFactor(Energy2 q2,int iloc, int id0, int id1,
0091                       Energy m0, Energy m1,
0092                       Complex & f1v,Complex & f2v,Complex & f3v,
0093                       Complex & f1a,Complex & f2a,Complex & f3a,
0094                       FlavourInfo flavour,
0095                       Virtuality virt=SpaceLike);
0096   //@}
0097
0098   /**
0099    * Output the setup information for the particle database
0100    * @param os The stream to output the information to
0101    * @param header Whether or not to output the information for MySQL
0102    * @param create Whether or not to add a statement creating the object
0103    */
0104   virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
0105
0106 protected:
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   LightBaryonQuarkModelFormFactor & operator=(const LightBaryonQuarkModelFormFactor &) = delete;
0124
0125 private:
0126
0127   /**
0128    * The form-factor \f$f_1\f$ at \f$q^2=0\f$.
0129    */
0130   vector<double> _f1;
0131
0132   /**
0133    * The form-factor \f$f_2\f$ at \f$q^2=0\f$.
0134    */
0135   vector<InvEnergy> _f2;
0136
0137   /**
0138    * The form-factor \f$g_1\f$ at \f$q^2=0\f$.
0139    */
0140   vector<double> _g1;
0141
0142   /**
0143    * The form-factor \f$g_2\f$ at \f$q^2=0\f$.
0144    */
0145   vector<InvEnergy> _g2;
0146
0147   /**
0148    * The mass for the energy dependence of the \f$f_1\f$ form factor.
0149    */
0150   vector<Energy> _Lambdaf1;
0151
0152   /**
0153    * The mass for the energy dependence of the \f$f_2\f$ form factor.
0154    */
0155   vector<Energy> _Lambdaf2;
0156
0157   /**
0158    * The mass for the energy dependence of the \f$g_1\f$ form factor.
0159    */
0160   vector<Energy> _Lambdag1;
0161
0162   /**
0163    * The mass for the energy dependence of the \f$g_2\f$ form factor.
0164    */
0165   vector<Energy> _Lambdag2;
0166
0167 };
0168
0169 }
0170
0171
0172 #endif /* HERWIG_LightBaryonQuarkModelFormFactor_H */