Simulation of LHCb RICH Detectors using