Decay/FormFactors/PiPiAnisovichKMatrix.h

0001 // -*- C++ -*-
0002 #ifndef Herwig_PiPiAnisovichKMatrix_H
0003 #define Herwig_PiPiAnisovichKMatrix_H
0004 //
0005 // This is the declaration of the PiPiAnisovichKMatrix class.
0006 //
0007
0008 #include "KMatrix.h"
0009
0010 namespace Herwig {
0011
0012 using namespace ThePEG;
0013
0014 /**
0015  * The PiPiAnisovichKMatrix class implements the K-matrix for $\pi\pi$ scattering from Eur.Phys.J.A 16 (2003) 229-258.
0016  *
0017  * @see \ref PiPiAnisovichKMatrixInterfaces "The interfaces"
0018  * defined for PiPiAnisovichKMatrix.
0019  */
0020 class PiPiAnisovichKMatrix: public KMatrix {
0021
0022 public:
0023
0024   /**
0025    * The default constructor.
0026    */
0027   PiPiAnisovichKMatrix();
0028
0029   /**
0030    *  Compute the K-matrix for a given scale
0031    * @param s The scale
0032    * @param Whether or not to multiply by \f$\prod_i(1-s/m^2_i\f$ to regularise the poles
0033    */
0034   virtual boost::numeric::ublas::matrix<double> K(Energy2 s, bool multiplyByPoles=false) const;
0035
0036 public:
0037
0038   /** @name Functions used by the persistent I/O system. */
0039   //@{
0040   /**
0041    * Function used to write out object persistently.
0042    * @param os the persistent output stream written to.
0043    */
0044   void persistentOutput(PersistentOStream & os) const;
0045
0046   /**
0047    * Function used to read in object persistently.
0048    * @param is the persistent input stream read from.
0049    * @param version the version number of the object when written.
0050    */
0051   void persistentInput(PersistentIStream & is, int version);
0052   //@}
0053
0054   /**
0055    * The standard Init function used to initialize the interfaces.
0056    * Called exactly once for each class by the class description system
0057    * before the main function starts or
0058    * when this class is dynamically loaded.
0059    */
0060   static void Init();
0061
0062 protected:
0063
0064   /** @name Clone Methods. */
0065   //@{
0066   /**
0067    * Make a simple clone of this object.
0068    * @return a pointer to the new object.
0069    */
0070   virtual IBPtr clone() const;
0071
0072   /** Make a clone of this object, possibly modifying the cloned object
0073    * to make it sane.
0074    * @return a pointer to the new object.
0075    */
0076   virtual IBPtr fullclone() const;
0077   //@}
0078
0079 protected:
0080
0081   /** @name Standard Interfaced functions. */
0082   //@{
0083   /**
0084    * Initialize this object after the setup phase before saving an
0085    * EventGenerator to disk.
0086    * @throws InitException if object could not be initialized properly.
0087    */
0088   virtual void doinit();
0089   //@}
0090
0091 private:
0092
0093   /**
0094    * The assignment operator is private and must never be called.
0095    * In fact, it should not even be implemented.
0096    */
0097   PiPiAnisovichKMatrix & operator=(const PiPiAnisovichKMatrix &) = delete;
0098
0099 private :
0100
0101   /**
0102    *  Parameters for the \f$K\f$-matrix
0103    */
0104   //{@
0105   /**
0106    *   \f$s_0^{\rm scat}
0107    */
0108   Energy2 s0Scatt_;
0109
0110   /**
0111    *   Constants
0112    */
0113   vector<double> f1a_;
0114
0115   /**
0116    *  \f$s_A\f$
0117    */
0118   double sA_;
0119
0120   /**
0121    *  \f$s_{A0}\f$
0122    */
0123   double sA0_;
0124
0125   /**
0126    *  Pion mass
0127    */
0128   Energy mPi_;
0129   //@}
0130 };
0131
0132 }
0133
0134 #endif /* Herwig_PiPiAnisovichKMatrix_H */