TEvtGen/EvtGenModels/EvtVubdGamma.hh

   1 //-----------------------------------------------------------------------
   2 // File and Version Information:
   3 //
   4 // Copyright Information: See EvtGen/COPYRIGHT
   5 //
   6 //
   7 // Description:
   8 //    3                                         2                2
   9 //   d Gamma                 /   _      _ _2  mb           _2  mb
  10 //  ---------- = 12 Gamma   | (1+x-z)(z-x-p ) -- W  + (1-z+p ) -- W
  11 //         _ 2           0   \                 2  1             2  2
  12 //  dx dz dp                                   2
  13 //                                _   _  _2  mb                 2   \.
  14 //                             + [x(z-x)-p ] -- (W + 2mb W  + mb W ) |
  15 //                                            4   3       4       5 /
  16 //
  17 //   with
  18 //        2 E           2
  19 //           l    _2   p        2 v.p    _
  20 //   x = ------ , p = --- , z = ------ , x = 1-x
  21 //         mb           2         mb
  22 //                    mb
  23 //
  24 //   the triple differential decay rate according to
  25 //   hep-ph/9905351 v2
  26 //
  27 // Environment:
  28 //      Software developed for the BaBar Detector at the SLAC B-Factory.
  29 //
  30 // Author List:
  31 //      Sven Menke
  32 //
  33 //-----------------------------------------------------------------------
  34
  35 #ifndef EVTVUBDGAMMA_HH
  36 #define EVTVUBDGAMMA_HH
  37
  38 class EvtVubdGamma {
  39
  40 public:
  41
  42   // Constructors
  43
  44   EvtVubdGamma(const double &alphas);
  45
  46   // Destructor
  47
  48   virtual ~EvtVubdGamma( );
  49
  50   // Operators
  51
  52   // Selectors
  53
  54   // Modifiers
  55
  56   // Methods
  57
  58   double getdGdxdzdp(const double &x, const double &z, const double &p2);
  59
  60 protected:
  61
  62   // Helper functions
  63
  64   double delta(const double &x, const double &xmin, const double &xmax);
  65
  66   double getW1nodelta(const double &x, const double &z, const double &p2);
  67
  68   double getW2nodelta(const double &x, const double &z, const double &p2);
  69
  70   double getW3nodelta(const double &x, const double &z, const double &p2);
  71
  72   double getW4nodelta(const double &x, const double &z, const double &p2);
  73
  74   double getW5nodelta(const double &x, const double &z, const double &p2);
  75
  76   double getW1delta(const double &x, const double &z);
  77
  78   double getW4plus5delta(const double &x, const double &z);
  79
  80 private:
  81
  82   // Friends
  83
  84   // Data members
  85
  86   double _alphas;
  87   double _epsilon1;
  88   double _epsilon2;
  89   double _epsilon3;
  90 };
  91
  92
  93 #endif // EVTVUBDGAMMA_HH
  94
  95