Neutrinos are elementary particles in the Standard Model. Neutrino oscillation is a quantum mechanical phenomenon beyond the Standard Model. Neutrino oscillation can be described by two independent mass-squared differences △m21^2, △m31^2 （or △m32^2） and a 3× 3 unitary matrix, containing three mixing angles θ12, θ23, θ13, and one charge-parity （CP） phase. θ12 is about 34° and determined by solar neutrino experiments and the reactor neutrino experiment KamLAND. θ23 is about 45° and determined by atmospheric neutrino experiments and accelerator neutrino experiments. 013 can be measured by either accelerator or reactor neutrino experiments. On Mar. 8, 2012, the Daya Bay Reactor Neutrino Experiment reported the first observation of non-zero 013 with 5.2 standard deviations. In June, with 2.5× previous data, Daya Bay improved the measurement of sin2 2013 = 0.089 ± 0.010（stat）± 0.005（syst）.
In reactor neutrino experiments, the analysis of time correlations between different physical events is an important task. Such analysis can help to understand the physical mechanisms of the signal and background events as well as the details of event selection and background estimation. This study investigates a ＂sampling and mixing＂ method used for producing large MC data samples for the Daya Bay reactor neutrino experiment. We designed a simple, generic mixing algorithm and generated large MC data samples for physics analysis from several samples according to their respective event rates. Basic plots based on the mixed data are shown.
