Implement B -> D^() K^ and B_s -> D_s^(*) rho QCD factorization amplitudes #799

smeiser · 2024-02-12T15:53:21Z

Take care of requested changes
Clarify how to deal with the form-factors

following [BBNS:2000A]

mreboud

Thanks @smeiser, I left some comments, nothing big

mreboud · 2024-03-11T09:46:15Z

eos/parameters/hadronic-matrix-elements.yaml

@@ -821,35 +821,35 @@ groups:
          max:     0.0
          latex:    '$a_{T,3}^{B\to D,\mathrm{BGL97}}$'
          unit:     '1'
-      # form factor parameters for B_s->D_s inspired by [BSZ:2015A]
+       # form factor parameters for B_s->D_s inspired by [BSZ:2015A]


The new indentation looks weird

mreboud · 2024-03-11T09:51:05Z

eos/parameters/hadronic-matrix-elements.yaml

-          central:  0.0
-          min:      0.0
-          max:      0.0
+          central:  0.6657


The new policy is to set all the form factors value to zero and to prompt a one-time warning whenever default parameters are used (#791). We can discuss it.

That is true for specialised parametrisations, such as the HQE. BSZ2015 / simplified series expansion should be non-zero, shouldn't it?

eos/b-decays/Makefile.am

mreboud · 2024-03-11T09:54:37Z

eos/b-decays/bq-to-dq-vec.cc

+
+        std::shared_ptr<VectorLCDAs> lcdas;
+
+        SpecifiedOption opt_cp_conjugate;


This should be a BooleanOption

mreboud · 2024-03-11T09:56:31Z

eos/b-decays/bq-to-dq-vec.cc

+        complex<double> AVLL(const double & u, const double & z) const
+        {
+            // evaluate AVLL chosing the correct branch for physical z = m_c / m_b < 1
+            if((0.0 < z) && (z < 1.0))


coding style, missing space

mreboud · 2024-03-11T09:57:25Z

eos/b-decays/bq-to-dq-vec.cc

+            if((0.0 < z) && (z < 1.0))
+            {
+                return
+                    log(u * (1.0 - z * z)) / (1.0 - u * (1.0 - z * z)) - pow(log(u * (1.0 - z * z)), 2) - dilog(1.0 - u * (1.0 - z * z));


Can you use power_of instead of pow for all the integer powers?

Also, no space immediately after return.

mreboud · 2024-03-11T09:58:50Z

eos/b-decays/bq-to-dq-vec.cc

+            // evaluate AVLL chosing the correct branch for unphysical z = m_b / m_c > 1
+            else if (z > 1.0)
+            {
+                return


Coding style, no new line after return

mreboud · 2024-03-11T10:02:55Z

eos/b-decays/bq-to-dstarq-vec.cc

+
+        std::shared_ptr<VectorLCDAs> lcdas;
+
+        SpecifiedOption opt_cp_conjugate;


Should be a BooleanOption

mreboud · 2024-03-11T10:03:56Z

eos/b-decays/bq-to-dstarq-vec.cc

+        complex<double> AVLL(const double & u, const double & z) const
+        {
+            // evaluate AVLL chosing the correct branch for z = m_c / m_b < 1 and z = - m_c / m_b > -1
+            if((-1.0 < z) && (z < 1.0))


Coding style: space

mreboud · 2024-03-11T10:04:21Z

eos/b-decays/bq-to-dstarq-vec.cc

+            if((-1.0 < z) && (z < 1.0))
+            {
+                return
+                    log(u * (1.0 - z * z)) / (1.0 - u * (1.0 - z * z)) - pow(log(u * (1.0 - z * z)), 2) - dilog(1.0 - u * (1.0 - z * z));


integer pow -> power_of

dvandyk

Looks pretty good, with a few requests for changes.

dvandyk · 2024-03-12T11:50:21Z

eos/parameters/hadronic-matrix-elements.yaml

-          central:  0.0
-          min:      0.0
-          max:      0.0
+          central:  0.6657


That is true for specialised parametrisations, such as the HQE. BSZ2015 / simplified series expansion should be non-zero, shouldn't it?

dvandyk · 2024-03-12T11:54:29Z

eos/b-decays/bq-to-dq-vec.cc

+            if((0.0 < z) && (z < 1.0))
+            {
+                return
+                    log(u * (1.0 - z * z)) / (1.0 - u * (1.0 - z * z)) - pow(log(u * (1.0 - z * z)), 2) - dilog(1.0 - u * (1.0 - z * z));


Also, no space immediately after return.

dvandyk · 2024-03-12T11:54:49Z

eos/b-decays/bq-to-dq-vec.cc

+            // evaluate fVLL chosing the correct branch for physical z = m_c / m_b < 1
+            if ((0.0 < z) && (z < 1.0))
+            {
+                return


dvandyk · 2024-03-12T11:54:57Z

eos/b-decays/bq-to-dq-vec.cc

+            // evaluate fVLL chosing the correct branch for unphysical z = m_b / m_c > 1
+            else if (z > 1.0)
+            {
+                return


dvandyk · 2024-03-12T11:55:10Z

eos/b-decays/bq-to-dq-vec.cc

+            // evaluate ASLR chosing the correct branch for physical z = m_c / m_b < 1
+            if ((0.0 < z) && (z < 1.0))
+            {
+                return


dvandyk · 2024-03-12T11:57:55Z

eos/b-decays/bq-to-dq-vec.cc

+                -((8.0 * (4.0 * u + 3.0)) / (1.0 + z)) + (8.0 * (1.0 - z)) / (1.0 + z) * (ATLL(u, z) + ATLL(1.0 - u, z));
+        };
+
+        complex<double> phiV(const double & u, const double & mu) const


Can you document what this correspond to in the BBNS paper, eq. an equation? Isn't this one of the two-particle LCDAs?

dvandyk · 2024-03-12T11:59:00Z

eos/b-decays/bq-to-dq-vec.cc

+            const double TSLR_im = integrate<GSL::QAGS>(im_integrand_SLR, 0.0, 1.0);
+
+            const double TTLL_re = integrate<GSL::QAGS>(re_integrand_TLL, 0.0, 1.0);
+            const double TTLL_im = 0.0;


Can you comment why the identity TTLL_im = 0 holds?

dvandyk · 2024-03-12T12:01:06Z

eos/b-decays/bq-to-dstarq-vec.cc

+            // evaluate AVLL chosing the correct branch for z = m_b / m_c > 1 and z = - m_b / m_c < -1
+            else if ((z > 1.0) || (z < -1.0))
+            {
+                return


Ditto. Not indicating of all those again 😄

dvandyk · 2024-03-12T12:02:02Z

eos/parameters/hadronic-matrix-elements.yaml

+          central: 0.685068
+          min:     0.685068
+          max:     0.685068
+          latex:   '$f_+^{B->D}(q^2 = M_K^{*2})$'


Is this LaTeX label correct? q^2 = M_{K^*}^2, isn't it?

dvandyk · 2024-03-12T12:02:08Z

eos/parameters/hadronic-matrix-elements.yaml

+          central: 0.70023
+          min:     0.70023
+          max:     0.70023
+          latex:   '$a_0^{B->D^*}(q^2 = M_K^{*2})$'


smeiser force-pushed the BtoDV branch from 3018318 to dec8bd6 Compare February 19, 2024 09:24

dvandyk force-pushed the BtoDV branch from dec8bd6 to cd27824 Compare February 19, 2024 14:30

smeiser force-pushed the BtoDV branch from cd27824 to 8bf876a Compare February 21, 2024 09:52

mreboud requested a review from a team March 11, 2024 09:44

mreboud assigned smeiser Mar 11, 2024

mreboud added the enhancement Feature request label Mar 11, 2024

smeiser added 2 commits March 11, 2024 10:45

[parameters] Add default values for B_s->D_s form factors

9b835e2

[b-decays] Add B_s->D_s^(*) rho and B->D^(*) Kstar amplitude

ab4fe9e

following [BBNS:2000A]

mreboud force-pushed the BtoDV branch from 8bf876a to ab4fe9e Compare March 11, 2024 09:45

mreboud requested changes Mar 11, 2024

View reviewed changes

dvandyk requested changes Mar 12, 2024

View reviewed changes

smeiser marked this pull request as draft May 24, 2024 09:11

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Implement B -> D^() K^ and B_s -> D_s^(*) rho QCD factorization amplitudes #799

Implement B -> D^() K^ and B_s -> D_s^(*) rho QCD factorization amplitudes #799

smeiser commented Feb 12, 2024 •

edited

mreboud left a comment

mreboud Mar 11, 2024

mreboud Mar 11, 2024

dvandyk Mar 12, 2024

mreboud Mar 11, 2024

mreboud Mar 11, 2024

mreboud Mar 11, 2024

dvandyk Mar 12, 2024

mreboud Mar 11, 2024

mreboud Mar 11, 2024

mreboud Mar 11, 2024

mreboud Mar 11, 2024

dvandyk left a comment

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024

dvandyk Mar 12, 2024


		std::shared_ptr<VectorLCDAs> lcdas;

		SpecifiedOption opt_cp_conjugate;

Implement B -> D^(*) K^* and B_s -> D_s^(*) rho QCD factorization amplitudes #799

Are you sure you want to change the base?

Implement B -> D^(*) K^* and B_s -> D_s^(*) rho QCD factorization amplitudes #799

Conversation

smeiser commented Feb 12, 2024 • edited

mreboud left a comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

dvandyk left a comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Choose a reason for hiding this comment

Implement B -> D^() K^ and B_s -> D_s^(*) rho QCD factorization amplitudes #799

Implement B -> D^() K^ and B_s -> D_s^(*) rho QCD factorization amplitudes #799

smeiser commented Feb 12, 2024 •

edited