We propose and evaluate a new scatter simulation algorithm to correct for scattered events in 3D PET. The algorithm simulates both single and double scatter. We show that single and double scatter distribution can have different shapes. This may cause difficulties if the total scatter is estimated by scaling the single scatter contribution. The whole implementation is based on the STIR library (Software for Tomographic Image Reconstruction [1]).