How to calculate the value of accessibility

We used the following script to calculate the value of accessibility.
First, we find the region coordinates on which we focused using Blat search and then executed the script.

# -*- coding: utf-8 -*-

"""
Calculate average of signal whose the region the user focuse on from Bedgraph format.
"""
__version__ = "1.00"
__date__ = "25 Jun 2022"

import sys

def averageBedgraph(filename, chromosome, start, end):
    """
    @function   averageBedgraph();
    Calculate average of signal whose the region the user focuse on from Bedgraph format.
    @param  {string} filename :  Bedgraph file path
    @param  {string} chromosome :   chromosome number which the user focuse
    @param  {int}    start : start point
    @param  {int}    end :   end point
    """

    total = 0
    with open(filename) as lines:
        for line in lines:
            c, s, e, v = line.split()
            if c != chromosome:
                continue
            s, e, v = int(s), int(e), int(v)
            if s <= start and end <= e:
                total += (end - start) * v
            elif s <= start < e:
                total += (e - start) * v
            elif start < s and e <= end:
                total += (e - s) * v
            elif s < end <= e:
                total += (end - s) * v
            elif end < s:
                break
    print('average : %s' % (total / (end - start)))
    print('done')

if __name__ == '__main__':
    argvs = sys.argv
    argc = len(argvs)

    if (argc != 5):   # Checking input
        print("USAGE : python3 averageBedgraph.py <INPUT_FILE> <CHROMOSOME> <START> <END>")
        quit()

    averageBedgraph(str(argvs[1]),str(argvs[2]),int(argvs[3]),int(argvs[4]))
quit()

Example of usage:

python3 averageBedgraph.py example.bedgraph chr1 10150 10250

Output of usage:

average : 12.74
done