Caption: Flowchart of the main procedures in PISCOLA. Gaussian Process is used to fit the light curves of SNe Ia. Afterwards, an SED template is warped to match the observed SN colors, which is then corrected for redshift (K-correction) and Milky Way dust extinction.