#params: <peak start> <peak stop> <shielded> <unshielded>

shielded_total=0
unshielded_total=0

for i in `seq $1 $2`; do
  shielded=`grep -E "^$i	" "$3" | cut -f 2`
  unshielded=`grep -E "^$i	" "$4" | cut -f 2`
  func=`echo -n $shielded_total + $shielded | tr -d "\r\l"`
  shielded_total=$(( $func ))
  func=`echo -n $unshielded_total + $unshielded | tr -d "\r\l"`
  unshielded_total=$(( $func ))
  echo -ne "\r$i"
done

echo -ne "\r"
echo "Shielded total: $shielded_total"
echo "Unshielded total: $unshielded_total"
echo -n "Linear coefficient μ: "
qalc "-1/0.009*ln($shielded_total/$unshielded_total)"