def GetInteractionRate(f,tallyNum=54,src=2.3E3):
    """
    Returns the interaction rate of the mctal file
    """
    m = mctal.MCTAL(f)
    t = m.tallies[tallyNum]
    return (t.data[-1]*src,t.errors[-1]*t.data[-1]*src)