Pith. sign in

REVIEW 3 major objections 6 minor 32 references

C-Code Reader of Form Factors

T0 review · 3 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A C reimplementation of the lattice QCD form-factor reader for $\Lambda_b \to \Lambda_c \ell^- \bar{\nu}_\ell$ matches the Mathematica output while running about 64 times faster per data point.

desk verdict A plausible software port with a useful speedup, but the 'identical results' claim rests on visual agreement only and needs numerical validation before the tool can be trusted. read the letter →

arxiv 2412.11437 v1 pith:V5BJV5CN submitted 2024-12-16 hep-ph

Ahmed Rashed This is my paper · ORCID
classification hep-ph
keywords CcodelatticeQCDformfactorsLambda_btoLambda_cdecaydifferentialrateMathematicacomparisonnewphysicscouplingswebapplicationerrorpropagation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that a C-language program can read the lattice QCD form factors for the decay $\Lambda_b \to \Lambda_c \ell^- \bar{\nu}_\ell$ and generate the same differential decay-rate data as a previous Mathematica implementation, but about 64 times faster per data point. The program implements the full set of scalar, vector/axial-vector, and tensor form factors and the effective Hamiltonian with new-physics couplings, including error propagation from the lattice covariance matrices. The reported agreement with the earlier code is based on plotted curves; if the agreement is exact, the practical gain is that new-physics parameter scans and fits become much cheaper and the analysis is easier to embed in other tools.

What carries the argument

The carrying object is the C program itself, organized as main.c with supporting modules for constants, calculations, and file manipulation. It maps the data flow of the Mathematica code into custom read-file functions for the lattice covariance matrices, then computes nominal and higher-order form factors and the hadronic helicity amplitudes entering the decay rate. Error propagation is handled by numerical differentiation, and the output is a tabulated file of $q^2$, differential decay rate, and error. A companion web interface runs the same C program on a server through a graphical input page.

What would settle it

Run both codes on the same covariance files and compare every printed value of $q^2$, the differential decay rate, and its error using a numerical tolerance such as a relative difference below $10^{-12}$; any point above that tolerance would show that the C code does not reproduce the Mathematica results, making the speed comparison less meaningful.

Watch

Extended reading notes

Core claim

The central claim is that the C implementation is a faithful, much faster replacement for the Mathematica reader used in Ref. [1] for the lattice QCD form factors of $\Lambda_b \to \Lambda_c$. The code reads the nominal and higher-order covariance files, reconstructs the ten form factors $F_0,F_+,F_\perp,G_0,G_+,G_\perp,h_+,h_\perp,\tilde{h}_+,\tilde{h}_\perp$, computes the differential decay rate for the Standard Model and for new-physics scenarios such as $g_L=-2.2$, and propagates uncertainties. The paper reports "identical results" at the level of the plotted differential distributions and a wall-clock ratio of 1:64.2 per data point in favor of the C code.

Load-bearing premise

The verification of "identical results" rests on visual agreement between plotted curves, so a discrepancy too small to see in the figures would not be caught, and the central speedup claim is meaningful only if the numerical outputs actually match.

Editorial extensions

If this is right

  • New-physics scans over the couplings $g_L,g_R,g_S,g_P,g_T$ become roughly two orders of magnitude cheaper, allowing thousands of $q^2$ points per parameter choice.
  • The analysis becomes reproducible without a Mathematica license, since the C code and its data files are distributed openly and the web interface hides the compilation step.
  • The same structured reader can be adapted to other lattice form-factor data sets that publish covariance matrices, extending the speed gain to related decays.
  • If the planned LaTeX-expression interface is completed, the tool generalizes from this decay to arbitrary user-provided expressions and data files.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • An implication the author leaves implicit is that the speedup makes Monte Carlo or profile-likelihood fits to LHCb's $\Lambda_b$ decay ratios practical, because a single likelihood evaluation drops from seconds to tens of milliseconds.
  • A testable extension is to add a numerical regression test to the released code that compares its output against a stored reference file from the Mathematica code with a fixed tolerance, converting the visual check into a reproducible assertion.
  • The architecture is portable to other flavor transitions: replacing the form-factor parametrization and hadron masses would give the same speed advantage, though the Hamiltonian and helicity amplitudes would need to be rewritten.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The manuscript reports a C-language implementation that reads the lattice QCD form-factor data for Lambda_b -> Lambda_c l nu from Ref. [1] and computes differential decay rates for the tau mode. The code is a port of the Mathematica code used in Ref. [1]. The paper claims that the C code produces identical results to the Mathematica code and that it is about 64.2 times faster per data point (Section 3, Section 4, Table 1). It also describes a prototype Node.js web interface and plans for a generalized expression parser. The physics content is limited to reproducing the existing calculation; the central claims are correctness of the port and computational speed.

Significance. If the claims are substantiated, the code would be a useful practical tool for the heavy-flavor phenomenology community: lattice form-factor analyses of Lambda_b -> Lambda_c tau nu often require repeated evaluations, and a fast, open-source reader would lower the barrier to using the Ref. [1] data files. The paper's strengths are that the code is released on GitHub under GPL, the modular structure is described, and the speed comparison is performed at four different data-set sizes. However, the scientific value of the paper depends entirely on establishing that the C code reproduces the Mathematica results to numerical accuracy and that the speed measurement is meaningful. The current evidence for both is thin, so the significance is conditional on additional validation.

major comments (3)
  1. [Section 4 and Section 6] The claim that the C code gives 'identical results' to the Mathematica code is not supported by quantitative evidence. Section 4 says only that the graphs 'align closely' and 'confirm the correctness of the C implementation,' and Section 6 repeats 'identical results.' There is no numerical comparison: no table of q2 values, no maximum absolute difference, no root-mean-square difference, and no tolerance for the decay rates or their error bars. Because the appendix states that errors are computed using numerical differentiation and error propagation, small differences in interpolation or covariance handling could produce offsets that are invisible at the plotted line width. Please provide a quantitative equivalence test, for example a machine-readable diff of the C and Mathematica outputs at all q2 points, together with the maximum and mean absolute differences and a stated tolerance.
  2. [Section 4, Table 1, and Figure 4] The 64.2x speed ratio is derived from slopes of wall-clock-time plots, but the measurement is weakly documented. The text says five trials were averaged, but no variance, standard deviation, or individual trial times are reported, and Figure 4 shows no error bars. No hardware specifications, operating system, compiler version, compiler optimization flags, Mathematica version, or timing methodology are given, so the measurement cannot be reproduced. In addition, it is not stated whether the 10-point measurement includes Mathematica kernel startup, which would affect the slope interpretation. Please report per-trial timings, mean and standard deviation, the exact versions of all software, and a description of how wall-clock time was measured.
  3. [Abstract, Section 3, and Conclusion] The code availability statement is internally inconsistent: the abstract gives the URL as https://github.com/darkfiresmith96/Lattice_QCD while Section 3 gives https://github.com/darkfiresmith96/LatticeQCD (with visible spacing artifacts in the rendered text), and no version or commit hash is provided. Since the paper's central deliverable is a code release, a single working URL with a specific commit identifier, and a note on how to reproduce the benchmark from that version, is necessary. The current discrepancy makes the code effectively unreproducible as published.
minor comments (6)
  1. [Table 1 and Figure 4 captions] The phrase 'Results of of the speed test' appears in both the table caption and the Figure 4 caption; delete the duplicated 'of.'
  2. [Figure 4 caption] The caption says 'Left and right panels are the linear and log plots, consequently,' but 'consequently' should be 'respectively.'
  3. [Figures 2 and 3] The panels are described as 'Standard Model results' and 'new physics results,' but the new physics point (g_L = -2.2) is only stated in the figure captions; the ranges of the axes and the number of q2 points plotted would be clearer if stated in the text.
  4. [Acknowledgments] The acknowledgments refer to 'B.S.' without identifying this person in the author list, and the grant attribution is not linked to a named author; this should be clarified.
  5. [Section 5] The web application is described as a prototype hosted on a free-tier Heroku instance, but no URL or repository location is given for the web code. Please provide access information or state that the web interface is not part of the released repository.
  6. [References] Reference [6] is a conference talk and does not include a DOI or arXiv identifier; please verify that this citation is complete.

Circularity Check

0 steps flagged · score 1.0 of 10

No circular derivation is present; the C code is a port of the authors' own Mathematica reader, so the reference comparison is internal, but the speed claim is an independent measurement and no result is fitted or defined in terms of its input.

full rationale

This paper makes no physics prediction and contains no fitted parameter that is later relabeled as a prediction. The C code is a port of the Mathematica reader from Ref. [1] (which shares author Rashed), and the claimed verification is that the C output 'align[s] closely' with the Mathematica output in Figs. 2 and 3 and is 'consistent with Ref. [1]' (Section 4). For a porting task, agreement with the source implementation is the appropriate reference check, and the 64.2x speedup is an independent wall-clock measurement, not a quantity derived from the Mathematica output by construction. The underlying lattice form factors come from external lattice QCD work (Refs. [27,28]), and no equation in the paper has the target result equal to an input by definition. The absence of numerical tolerances in the visual comparison is a validation rigor concern, not circularity. Under the rule that self-citation is not circularity unless load-bearing, no specific circular step can be exhibited.

Assumptions & free parameters 0 free parameters · 3 assumptions · 0 invented entities

The paper's physics content is entirely imported from prior work: the effective Hamiltonian, the lattice form factors, and the Mathematica reference code. The only new artifact is the C implementation, which adds no free parameters or invented entities. The key assumptions are the correctness of the imported physics and of the Mathematica baseline.

assumptions (3)
  • domain assumption The effective Hamiltonian in Section 2 correctly describes b -> c l nu transitions with NP scalar, vector, and tensor couplings.
    The paper uses this Hamiltonian to define the decay rate calculation; it is taken from Refs. [29,30,31].
  • domain assumption The lattice QCD form factors and covariance matrices from Refs. [1,27] are accurate inputs.
    The code reads these supplemental data files and the paper does not re-derive or validate the lattice results.
  • domain assumption The Mathematica code of Ref. [1] is a correct reference implementation.
    Correctness of the C port is asserted by matching this code, but the Mathematica code itself is not independently verified in this paper.

how reviews work

0 comments
Cite this review

Pith. "Pith review of C-Code Reader of Form Factors." pith.science (2026). https://pith.science/paper/V5BJV5CN

@misc{pith2026241211437,
  author       = {Pith},
  title        = {Pith review of: C-Code Reader of Form Factors},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/V5BJV5CN}},
  note         = {Machine review of arXiv:2412.11437}
}
abstract

The process $\Lambda_b \rightarrow \Lambda_c \ell^- \bar{\nu}_\ell$ serves as a tool for exploring new physics, with contributions from scalar, vector, and tensor hadronic currents in various models. These form factors are derived from the quark model or lattice QCD. This work introduces a C-code for efficiently reading lattice QCD form factors for these currents, significantly outperforming a previous Mathematica-based approach, with a speed improvement of over 64 times per data point. The code, available on GitHub $(https://github.com/darkfiresmith96/Lattice_QCD)$, also features a web interface for user inputs.

Figures

Figures reproduced from arXiv: 2412.11437 by the authors.

Figure 1
Figure 1. Example of data from code 4 6 8 10 0.0 0.5 1.0 1.5 q2 [GeV2 ] d Γ/dq 2 [x10 -15 GeV -1 ] 4 6 8 10 0.0 0.5 1.0 1.5 2.0 q2 [GeV2 ] d Γ/dq 2 [x10 -15 GeV -1 ] [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. The differential cross section distribution from the Mathematica code. The left panel is the [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. The differential cross section distribution from the c-code. The left panel is the Standard [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (3 more)
Figure 4
Figure 4. Figure 4: Results of of the speed test. This shows that the c-code (red) performs faster than the [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]
Figure 5
Figure 5. Figure 5: Web interface for user inputs [PITH_FULL_IMAGE:figures/full_fig_p007_5.png]
Figure 6
Figure 6. Figure 6: Web application which can parse user’s expression [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

32 extracted references · 13 canonical work pages

  1. [1]

    Datta, S

    A. Datta, S. Kamali, S. Meinel and A. Rashed, JHEP 1708, 131 (2017) doi:10.1007/JHEP08(2017)131 [arXiv:1702.02243 [hep-ph]]

    arXiv 2017

  2. [2]

    Non-relativistic Model for the Semileptonic \Lambda_b \to \Lambda_c Decay

    D. Chakraverty, T. De, B. Dutta-Roy and K. S. Gupta, Mod. Phys. Lett. A 12, 195 (1997) doi:10.1142/S0217732397000194 [hep-ph/9612369]

    work page Pith review arXiv 1997

  3. [3]

    Combined nonrelativistic constituent quark model and heavy quark effective theory study of semileptonic decays of \Lambda_b and \Xi_b baryons

    C. Albertus, E. Hernandez and J. Nieves, Phys. Rev. D 71, 014012 (2005) doi:10.1103/PhysRevD.71.014012 [nucl-th/0412006]

    work page Pith review arXiv 2005

  4. [4]

    Shivashankara, W

    S. Shivashankara, W. Wu and A. Datta, Phys. Rev. D 91, no. 11, 115003 (2015) doi:10.1103/PhysRevD.91.115003 [arXiv:1502.07230 [hep-ph]]

    arXiv 2015

  5. [5]

    J. P. Lees et al. [BaBar Collaboration], Phys. Rev. D 88, no. 7, 072012 (2013) doi:10.1103/PhysRevD.88.072012 [arXiv:1303.0571 [hep-ex]]

    arXiv 2013

  6. [6]

    Ciezarek [LHCb Collaboration], talk presented atFlavor Physics & CP violation 2015(Nagoya, Japan,25-29 May 2015), [fpcp2015.hepl.phys.nagoya-u.ac.jp]

    G. Ciezarek [LHCb Collaboration], talk presented atFlavor Physics & CP violation 2015(Nagoya, Japan,25-29 May 2015), [fpcp2015.hepl.phys.nagoya-u.ac.jp]

    work page 2015

  7. [7]

    Bhattacharya, A

    B. Bhattacharya, A. Datta, D. London and S. Shivashankara, Phys. Lett. B 742, 370 (2015) doi:10.1016/j.physletb.2015.02.011 [arXiv:1412.7164 [hep-ph]]

    arXiv 2015

  8. [8]

    Bhattacharya, A

    B. Bhattacharya, A. Datta, J. P. Guevin, D. London and R. Watanabe, JHEP 1701, 015 (2017) doi:10.1007/JHEP01(2017)015 [arXiv:1609.09078 [hep-ph]]

    arXiv 2017

Show all 32 references
  1. [9]

    Aaij et al

    R. Aaij et al. [LHCb Collaboration], Phys. Rev. Lett. 113, 151601 (2014) doi:10.1103/PhysRevLett.113.151601 [arXiv:1406.6482 [hep-ex]]

    2014 arXiv

  2. [10]

    A. K. Alok, B. Bhattacharya, A. Datta, D. Kumar, J. Kumar and D. Lon- don, Phys. Rev. D 96, no. 9, 095009 (2017) doi:10.1103/PhysRevD.96.095009 [arXiv:1704.07397 [hep-ph]]

    2017 arXiv

  3. [11]

    Datta, J

    A. Datta, J. Kumar and D. London, Phys. Lett. B 797, 134858 (2019) doi:10.1016/j.physletb.2019.134858 [arXiv:1903.10086 [hep-ph]]

    2019

  4. [12]

    Datta, J

    A. Datta, J. Liao and D. Marfatia, Phys. Lett. B 768, 265 (2017) doi:10.1016/j.physletb.2017.02.058 [arXiv:1702.01099 [hep-ph]]

    2017 arXiv

  5. [13]

    Datta, J

    A. Datta, J. Kumar, J. Liao and D. Marfatia, Phys. Rev. D 97, no. 11, 115038 (2018) doi:10.1103/PhysRevD.97.115038 [arXiv:1705.08423 [hep-ph]]

    2018 arXiv

  6. [14]

    Datta, hep-ph/9504429

    A. Datta, hep-ph/9504429

    arXiv

  7. [15]

    Aaij et al

    R. Aaij et al. [LHCb Collaboration], Nature Phys. 11, 743 (2015) doi:10.1038/nphys3415 [arXiv:1504.01568 [hep-ex]]

    2015 arXiv

  8. [16]

    Cardarelli and S

    F. Cardarelli and S. Simula, Phys. Lett. B 421, 295 (1998) doi:10.1016/S0370- 2693(97)01581-5 [hep-ph/9711207]

    1998 arXiv

  9. [17]

    H. G. Dosch, E. Ferreira, M. Nielsen and R. Rosenfeld, Phys. Lett. B 431, 173 (1998) doi:10.1016/S0370-2693(98)00566-8 [hep-ph/9712350]. 8

    1998 arXiv

  10. [18]

    C. S. Huang, C. F. Qiao and H. G. Yan, Phys. Lett. B 437, 403 (1998) doi:10.1016/S0370-2693(98)00909-5 [hep-ph/9805452]

    1998 arXiv

  11. [19]

    R. S. Marques de Carvalho, F. S. Navarra, M. Nielsen, E. Ferreira and H. G. Dosch, Phys. Rev. D 60, 034009 (1999) doi:10.1103/PhysRevD.60.034009 [hep-ph/9903326]

    1999 arXiv

  12. [20]

    M. q. Huang and D. W. Wang, Phys. Rev. D 69, 094003 (2004) doi:10.1103/PhysRevD.69.094003 [hep-ph/0401094]

    2004 arXiv

  13. [21]

    Pervin, W

    M. Pervin, W. Roberts and S. Capstick, Phys. Rev. C 72, 035201 (2005) doi:10.1103/PhysRevC.72.035201 [nucl-th/0503030]

    2005 arXiv

  14. [22]

    H. W. Ke, X. Q. Li and Z. T. Wei, Phys. Rev. D 77, 014020 (2008) doi:10.1103/PhysRevD.77.014020 [arXiv:0710.1927 [hep-ph]]

    2008 arXiv

  15. [23]

    Y. M. Wang, Y. L. Shen and C. D. Lu, Phys. Rev. D 80, 074012 (2009) doi:10.1103/PhysRevD.80.074012 [arXiv:0907.4008 [hep-ph]]

    2009 arXiv

  16. [24]

    Azizi, M

    K. Azizi, M. Bayar, Y. Sarac and H. Sundu, Phys. Rev. D 80, 096007 (2009) doi:10.1103/PhysRevD.80.096007 [arXiv:0908.1758 [hep-ph]]

    2009 arXiv

  17. [25]

    Khodjamirian, C

    A. Khodjamirian, C. Klein, T. Mannel and Y.-M. Wang, JHEP 1109, 106 (2011) doi:10.1007/JHEP09(2011)106 [arXiv:1108.2971 [hep-ph]]

    2011 arXiv

  18. [26]

    Gutsche, M

    T. Gutsche, M. A. Ivanov, J. G. K¨ orner, V. E. Lyubovitskij and P. Santorelli, Phys. Rev. D 90, no. 11, 114033 (2014) Erratum: [Phys. Rev. D 94, no. 5, 059902 (2016)] doi:10.1103/PhysRevD.90.114033, 10.1103/PhysRevD.94.059902 [arXiv:1410.6043 [hep-ph]]

    2014 arXiv

  19. [27]

    Detmold, C

    W. Detmold, C. Lehner and S. Meinel, Phys. Rev. D 92, no. 3, 034503 (2015) doi:10.1103/PhysRevD.92.034503 [arXiv:1503.01421 [hep-lat]]

    2015 arXiv

  20. [28]

    Detmold and S

    W. Detmold and S. Meinel, Phys. Rev. D 93, no. 7, 074501 (2016) doi:10.1103/PhysRevD.93.074501 [arXiv:1602.01399 [hep-lat]]

    2016 arXiv

  21. [29]

    C. H. Chen and C. Q. Geng, Phys. Rev. D 71, 077501 (2005) doi:10.1103/PhysRevD.71.077501 [hep-ph/0503123]

    2005 arXiv

  22. [30]

    Bhattacharya, V

    T. Bhattacharya, V. Cirigliano, S. D. Cohen, A. Filipuzzi, M. Gonzalez-Alonso, M. L. Graesser, R. Gupta and H. W. Lin, Phys. Rev. D 85, 054512 (2012) doi:10.1103/PhysRevD.85.054512 [arXiv:1110.6448 [hep-ph]]

    2012 arXiv

  23. [31]

    Datta, M

    A. Datta, M. Duraisamy and D. Ghosh, Phys. Rev. D 86, 034027 (2012) doi:10.1103/PhysRevD.86.034027 [arXiv:1206.3760 [hep-ph]]

    2012 arXiv

  24. [32]

    Feldmann and M

    T. Feldmann and M. W. Y. Yip, Phys. Rev. D 85, 014035 (2012) Erratum: [Phys. Rev. D 86, 079901 (2012)] doi:10.1103/PhysRevD.85.014035, 10.1103/phys- revd.86.079901 [arXiv:1111.1844 [hep-ph]]. 9

    2012 arXiv

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.