From bed file to enrichment score

Question

[GK1610]

I have a consensus peak bed file for which I have to quantify the enrichment score of a bam file against their respective control file. What is the best way in macs2 software or any other software to get enrichment scores of given bed file?

Also, is the col. 7 in macs2 narrowpeak file is the enrichment score for a given peak regions averaged over the width of the peak?

[ATpoint]

The enrichment comes from the position with the highest fragment pileup in the ChIP against the local control lambda, so no, the width is not directly of relevance for the enrichment. Check in the original paper, it is quiet intuitive to read.