We have upgraded the new Tibet ASgamma experiment in China since 2014 to measure the chemical composition of cosmic rays around the knee. This hybrid experiment consist of an air-shower-core detector array (YAC-II) to detect high energy electromagnetic component, the Tibet air-shower array (Tibet-III) and a large underground water-Cherenkov muon-detector array (MD). We have carried out a detailed air-shower Monte Carlo (MC) simulation to study the performance of the hybrid detectors by using CORSIKA (version 7.5000), which includes EPOS-LHC, QGSJETII-04, SIBYLL2.1 and SIBYLL2.3 hadronic interaction models. The preliminary results of the interaction model checking above 50 TeV energy region are reported in this paper, and the primary proton and helium spectra in the energy range 50 TeV to 1015 eV was derived from YAC-I data and is smoothly connected with direct observation data at lower energies and also with our previously reported works at higher energies within statistical errors. The knee of the (P+He) spectra is located around 400 TeV. The interaction model dependence in deriving the primary (P+He) spectra is found to be small (less than 25% in absolute intensity, 10% in position of the knee), and the composition model dependence is less than 10% in absolute intensity.