Decay/FormFactors/KPiIHalfFOCUSKMatrix.h

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