Gamma ray tracking is an efficient detection technique in studying exotic nuclei which lies far from beta stability line. To achieve very powerful and extraordinary resolution ability, new detectors based on gamma ray tracking are currently being developed. To reach this achievement, the neutron-gamma discrimination in these detectors is also an important task. In this paper, by suitable layered feedforward neural networks (LFNNs), we have constructed novel and consistent empirical physical formulas (EPFs) for some highly nonlinear detector counts measured in neutron-gamma discrimination. The detector counts data used in the discrimination was actually borrowed from our previous paper. The counts used here had been originally measured versus the following parameters: energy deposited in the first interaction points, difference in the incoming direction of initial gamma rays, and finally figure of merit values of the clusters determined by tracking. The LFNN-EPFs are of explicit mathematical functional form. Therefore, by various suitable operations of mathematical analysis, these LENN-EPFs can be used to derivate further physical functions which might be potentially relevant to neutron-gamma discrimination performance of gamma ray tracking. (C) 2012 Elsevier Ltd. All rights reserved.