Question: peak to gene custom annotation

[JABioinf]

Dear Pando team,
I would like to use a complementary set of distal peaks to the "GREAT" peak selection methods (https://quadbio.github.io/Pando/articles/association.html). To do this, does the infer_grn() function accept both peak_to_gene_method and peak_to_gene_domains or they are exclusive?

Is there a way to extract the final peak-to-gene association used from the Seurat plus object as a grange object after running infer_grn() to be able to add in a second run the custom distal peak-gene association to it?

If not, could you explain how the GREAT method is implemented and can be reproduced manually? Does it consist just of extending the proximity test by "extend=" and intersecting the resulting peaks with proximal regions from neighboring genes (excluding exons and TSS) ? or is there additional step?

thank you for your help!