Analysis of first-order reversal curves (FORCs) is a powerful tool to probe irreversible switching events in nanomagnet assemblies. As in essence switching events are related to the intrinsic properties of the constituents and their interactions, the resulting FORC diagrams contain much information that can be cross-linked and complex to deconvolute. In order to quantify the relevant parameters that drive the FORC diagrams of arrays of perpendicularly magnetized nanomagnets, we present step-by-step simulations of assemblies of hysterons to determine the specific signatures related to different known inputs. While we explored the consequences of dipolar interactions using either mean field or magnetostatic approaches, we completed by taking the hysteron switching field distribution (SFD) as either normal or lognormal. We demonstrated that the transition between FORC diagrams composed of an isolated interaction field distribution (IFD) and a wishbone shape operates via the SFD deviation, σHsw, in the presence of a weakly dispersed interaction field. In the presence of a magnetostatic interaction field, the IFD profile is peaked and a coercive field distribution (CFD) sums to the IFD as σHsw increases. A transition between IFD + CFD and wishbone shapes is clearly demonstrated as a function of the interaction field deviation σHint. In addition, we demonstrate that whatever the considered cases, σHswcan be quantitatively extracted from the FORC diagrams within an error inferior to 10%. These findings are of interest for dipolar coupled perpendicularly magnetized nanomagnets, as in assemblies of magnetic nanowires and nanopillars, as well as bit patterned media.