Handwritten derivations

Numerical model for Steve's elegant delta_epsilon/epsilon integration

This script is used for a MC simulation of the relative efficiency
changes due to mp displacement by Rd, given an impact parameter resolution function Fb(b).
A high statistics calculation with 1E9 events is saved in a root file.

The generated 2D functions are shown in mc1.png

The slices in Rd and b are displayed in mc2.png:

The results agree with Steve's, apart from the first two bcut bins (0 < b < 10 and 10 < b < 20).