A large scale Monte Carlo production has been pursued since spring 2018 for the ILD detector optimization studies based on physics benchmark processes. A production system based on ILCDirac has been developed to produce samples in timely manner. The system and its performance are presented.