In-beam PET is one of the imaging-based methods for monitoring carbon therapy by reconstruction of positron-emitters produced by incident particles. The accurate Monte Carlo simulation is necessary for range verification and therapeutic dose monitoring by means of the in-beam PET. GATE is capable of performing in-beam PET imaging but very few studies assess the accuracy of activity range measured from PET imaging with GATE simulated results. In this work, we present both experimental and GATE simulation studies of in-beam PET imaging of a phantom, which is irradiated by carbon ion beam, for range verification. The experiment data is acquired at Wuwei Heavy Ion Cancer Treatment Center in Gansu, China, with a dual-head plate PET prototype. The PET prototype consists of 2 panels, each arranged in 2 × 2 matrix. Each detector block is composed of 20 × 20 LYSO crystal array coupled with the position-sensitive photomultiplier tube (HAMAMATSU H8500). A homogeneous plastic phantom is irradiated with monoenergetic 131.05 MeV and 190.19 MeV 12C ion pencil beams. The irradiation and full response of PET prototype are simulated with GATE macros. A novel mathematical model is established, which is capable of calculating and revealing the variation of the positron activity. Our focus is on the reconstructed activity ranges obtained by the experimental measurement and GATE MC prediction combined with the positron activity distribution calculation mathematical model. Results show that the reconstructed activity distribution of experimental data and GATE MC prediction are in good agreement. The measured and simulated 1D activity peak and falling edge positions are within 1.0 mm in all cases. The 20%, 50% and 80% PET peak positions predicted by GATE are close to measurements. These results indicated that it is feasible to assess the accuracy of activity range measured from the dual-head plate PET system with the modeled GATE hadron-PET simulation.