Decay/FormFactors/BallZwickyVectorFormFactor.h

0001 // -*- C++ -*-
0002 //
0003 // BallZwickyVectorFormFactor.h is a part of Herwig - A multi-purpose Monte Carlo event generator
0004 // Copyright (C) 2002-2019 The Herwig Collaboration
0005 //
0006 // Herwig is licenced under version 3 of the GPL, see COPYING for details.
0007 // Please respect the MCnet academic guidelines, see GUIDELINES for details.
0008 //
0009 #ifndef HERWIG_BallZwickyVectorFormFactor_H
0010 #define HERWIG_BallZwickyVectorFormFactor_H
0011 //
0012 // This is the declaration of the BallZwickyVectorFormFactor class.
0013 //
0014 #include "ScalarFormFactor.h"
0015 namespace Herwig {
0016 using namespace ThePEG;
0017
0018   /** \ingroup Decay
0019    *
0020    *  The BallZwickyVectorFormFactor class implements the form-factors
0021    *  of hep-ph/0412079 for the B meson to light vector mesons.
0022    *
0023    *  This class is one of the few which includes the penguin form factors in addition
0024    *  to the standard weak decay form factors.
0025    *
0026    * @see ScalarFormFactor
0027    * @see BallZwickyScalarFormFactor
0028    */
0029
0030 class BallZwickyVectorFormFactor: public ScalarFormFactor {
0031
0032 public:
0033
0034   /**
0035    * Default constructor
0036    */
0037   BallZwickyVectorFormFactor();
0038
0039   /** @name Form-Factors */
0040   //@{
0041   /**
0042    * The form factor for the weak decay of a scalar to a vector.
0043    * @param q2 The scale \f$q^2\f$.
0044    * @param iloc The location in the form factor list.
0045    * @param id0 The PDG code of the incoming meson.
0046    * @param id1 The PDG code of the outgoing meson.
0047    * @param m0 The mass of the incoming meson.
0048    * @param m1 The mass of the outgoing meson.
0049    * @param A0 The form factor \f$A_0\f$
0050    * @param A1 The form factor \f$A_1\f$
0051    * @param A2 The form factor \f$A_2\f$
0052    * @param V  The form factor \f$V\f$
0053    */
0054   void ScalarVectorFormFactor(Energy2 q2, unsigned int iloc, int id0, int id1,
0055                   Energy m0, Energy m1,Complex & A0,
0056                   Complex & A1,Complex & A2, Complex & V) const;
0057
0058   /**
0059    * The form factor for the weak penguin decay of a scalar meson to a vector meson.
0060    * @param q2 The scale \f$q^2\f$.
0061    * @param iloc The location in the form factor list.
0062    * @param id0 The PDG code of the incoming meson.
0063    * @param id1 The PDG code of the outgoing meson.
0064    * @param m0 The mass of the incoming meson.
0065    * @param m1 The mass of the outgoing meson.
0066    * @param T1 The form factor \f$T_1\f$.
0067    * @param T2 The form factor \f$T_2\f$.
0068    * @param T3 The form factor \f$T_3\f$.
0069    */
0070   void ScalarVectorSigmaFormFactor(Energy2 q2,unsigned int iloc,int id0,int id1,
0071                    Energy m0, Energy m1, Complex & T1,
0072                    Complex & T2, Complex & T3) const;
0073   //@}
0074
0075   /**
0076    * Output the setup information for the particle database
0077    * @param os The stream to output the information to
0078    * @param header Whether or not to output the information for MySQL
0079    * @param create Whether or not to add a statement creating the object
0080    */
0081   virtual void dataBaseOutput(ofstream & os,bool header,bool create) const;
0082
0083 public:
0084
0085   /** @name Functions used by the persistent I/O system. */
0086   //@{
0087   /**
0088    * Function used to write out object persistently.
0089    * @param os the persistent output stream written to.
0090    */
0091   void persistentOutput(PersistentOStream & os) const;
0092
0093   /**
0094    * Function used to read in object persistently.
0095    * @param is the persistent input stream read from.
0096    * @param version the version number of the object when written.
0097    */
0098   void persistentInput(PersistentIStream & is, int version);
0099   //@}
0100
0101   /**
0102    * Standard Init function used to initialize the interfaces.
0103    */
0104   static void Init();
0105
0106 protected:
0107
0108   /** @name Clone Methods. */
0109   //@{
0110   /**
0111    * Make a simple clone of this object.
0112    * @return a pointer to the new object.
0113    */
0114   virtual IBPtr clone() const {return new_ptr(*this);}
0115
0116   /** Make a clone of this object, possibly modifying the cloned object
0117    * to make it sane.
0118    * @return a pointer to the new object.
0119    */
0120   virtual IBPtr fullclone() const {return new_ptr(*this);}
0121   //@}
0122
0123 protected:
0124
0125   /** @name Standard Interfaced functions. */
0126   //@{
0127   /**
0128    * Initialize this object after the setup phase before saving and
0129    * EventGenerator to disk.
0130    * @throws InitException if object could not be initialized properly.
0131    */
0132   virtual void doinit();
0133   //@}
0134
0135 private:
0136
0137   /**
0138    * Private and non-existent assignment operator.
0139    */
0140   BallZwickyVectorFormFactor & operator=(const BallZwickyVectorFormFactor &) = delete;
0141
0142 private:
0143
0144   /** @name Coefficients for the form factors.*/
0145   //@{
0146
0147   /**
0148    * The coefficient \f$r_1\f$ for the \f$V(q^2)\f$ form factor.
0149    */
0150   vector<double> _Vr1;
0151
0152   /**
0153    * The coefficient \f$r_2\f$ for the \f$V(q^2)\f$ form factor.
0154    */
0155   vector<double> _Vr2;
0156
0157   /**
0158    * The coefficient \f$r_1\f$ for the \f$A_0(q^2)\f$ form factor.
0159    */
0160   vector<double> _A0r1;
0161
0162   /**
0163    * The coefficient \f$r_2\f$ for the \f$A_0(q^2)\f$ form factor.
0164    */
0165   vector<double> _A0r2;
0166
0167   /**
0168    * The coefficient \f$r_1\f$ for the \f$A_1(q^2)\f$ form factor.
0169    */
0170   vector<double> _A1r1;
0171
0172   /**
0173    * The coefficient \f$r_2\f$ for the \f$A_1(q^2)\f$ form factor.
0174    */
0175   vector<double> _A1r2;
0176
0177   /**
0178    * The coefficient \f$r_1\f$ for the \f$A_2(q^2)\f$ form factor.
0179    */
0180   vector<double> _A2r1;
0181
0182   /**
0183    * The coefficient \f$r_2\f$ for the \f$A_2(q^2)\f$ form factor.
0184    */
0185   vector<double> _A2r2;
0186
0187   /**
0188    * The coefficient \f$r_1\f$ for the \f$T_1(q^2)\f$ form factor.
0189    */
0190   vector<double> _T1r1;
0191
0192   /**
0193    * The coefficient \f$r_2\f$ for the \f$T_1(q^2)\f$ form factor.
0194    */
0195   vector<double> _T1r2;
0196
0197   /**
0198    * The coefficient \f$r_1\f$ for the \f$T_2(q^2)\f$ form factor.
0199    */
0200   vector<double> _T2r1;
0201
0202   /**
0203    * The coefficient \f$r_2\f$ for the \f$T_2(q^2)\f$ form factor.
0204    */
0205   vector<double> _T2r2;
0206
0207   /**
0208    * The coefficient \f$r_1\f$ for the \f$\tilde{T}_3(q^2)\f$ form factor.
0209    */
0210   vector<double> _T3r1;
0211
0212   /**
0213    * The coefficient \f$r_2\f$ for the \f$\tilde{T}_3(q^2)\f$ form factor.
0214    */
0215   vector<double> _T3r2;
0216   // the constants for the form-factors
0217   //@}
0218
0219   /** @name Masses for the form factors.*/
0220   //@{
0221
0222   /**
0223    * The mass \f$m_R^2\f$ for the \f$V(q^2)\f$ form factor.
0224    */
0225   vector<Energy2> _VmR2;
0226
0227   /**
0228    * The mass \f$m_{\rm fit}^2\f$ for the \f$V(q^2)\f$ form factor.
0229    */
0230   vector<Energy2> _Vmfit2;
0231
0232   /**
0233    * The mass \f$m_R^2\f$ for the \f$A_0(q^2)\f$ form factor.
0234    */
0235   vector<Energy2> _A0mR2;
0236
0237   /**
0238    * The mass \f$m_{\rm fit}^2\f$ for the \f$A_0(q^2)\f$ form factor.
0239    */
0240   vector<Energy2> _A0mfit2;
0241
0242   /**
0243    * The mass \f$m_R^2\f$ for the \f$A_1(q^2)\f$ form factor.
0244    */
0245   vector<Energy2> _A1mR2;
0246
0247   /**
0248    * The mass \f$m_{\rm fit}^2\f$ for the \f$A_1(q^2)\f$ form factor.
0249    */
0250   vector<Energy2> _A1mfit2;
0251
0252   /**
0253    * The mass \f$m_R^2\f$ for the \f$A_2(q^2)\f$ form factor.
0254    */
0255   vector<Energy2> _A2mR2;
0256
0257   /**
0258    * The mass \f$m_{\rm fit}^2\f$ for the \f$A_2(q^2)\f$ form factor.
0259    */
0260   vector<Energy2> _A2mfit2;
0261
0262   /**
0263    * The mass \f$m_R^2\f$ for the \f$T_1(q^2)\f$ form factor.
0264    */
0265   vector<Energy2> _T1mR2;
0266
0267   /**
0268    * The mass \f$m_{\rm fit}^2\f$ for the \f$T_1(q^2)\f$ form factor.
0269    */
0270   vector<Energy2> _T1mfit2;
0271
0272   /**
0273    * The mass \f$m_R^2\f$ for the \f$T_2(q^2)\f$ form factor.
0274    */
0275   vector<Energy2> _T2mR2;
0276
0277   /**
0278    * The mass \f$m_{\rm fit}^2\f$ for the \f$T_2(q^2)\f$ form factor.
0279    */
0280   vector<Energy2> _T2mfit2;
0281
0282   /**
0283    * The mass \f$m_R^2\f$ for the \f$\tilde{T}_3(q^2)\f$ form factor.
0284    */
0285   vector<Energy2> _T3mR2;
0286
0287   /**
0288    * The mass \f$m_{\rm fit}^2\f$ for the \f$\tilde{T}_3(q^2)\f$ form factor.
0289    */
0290   vector<Energy2> _T3mfit2;
0291   // the masses for the form-factors
0292   //@}
0293
0294   /**
0295    * Cut-off parameter for the switch to a small \f$q^2\f$ expansion for the \f$T_3\f$
0296    * form factor.
0297    */
0298   Energy2 _cutoff;
0299 };
0300
0301 }
0302
0303 #endif /* HERWIG_BallZwickyVectorFormFactor_H */