A crystal routine for collimation studies in circular proton accelerators

A routine has been developed to simulate interactions of protons with bent crystals in a version of SixTrack for collimation studies. This routine is optimized to produce high-statistics tracking simulations for a highly efficient collimation system, like the one of the CERN Large Hadron Collider (LHC). The routine has recently been reviewed and improved through a comparison with experimental data, benchmarked against other codes and updated by adding better models of low-probability interactions. In this paper, data taken with 400 GeV/c proton beams at the CERN-SPS North Area are used to verify the prediction of the routine, including the results of a more recent analysis.