In the present work we are interested in understanding various properties of quark matter subject to strong magnetic fields described by the Nambu-Jona-Lasinio model with Polyakov loop. We start by analysing the differences arising from two different vector interactions in the Lagrangian densities, at zero temperature, and apply the results to stellar matter. We then investigate the position of the critical end point for different chemical potential and density scenarios.