CNOK: A C++ Glauber model code for single-nucleon knockout reactions

A C++ program package was developed to calculate parallel momentum distributions (including the stripping mechanism and the diffractive dissociation mechanism) of the heavy residue (core) in single-nucleon knockout reactions induced by intermediate-energy (around and more than tens of MeV per nucleon) beams of stable and radioactive atomic nuclei. The program implements the Glauber reaction model that is based on the eikonal approximation and the sudden approximation when dealing with the scattering process, and uses a t-ρ-ρ method to build the projectile-target scattering optical potential, where the nucleonic densities of the core and the target are needed as input. The radial wavefunction of the valence nucleon is solved by a built-in bound-wave solver driven by a globally convergent search engine for optimized potential parameters of the valence nucleon to reproduce user-specified separation energies and root-mean-square (rms) radii of the valence nucleon. The YAML file format is adopted to facilitate user input. The program is expected to especially serve the interpretation of data analysis results of single-nucleon knockout from radioactive ion beams (RIBs), where the experimental practitioners may be more conversant with C++.