Added options + changed polynomial addition + C++Python bindings -> increased speed

I made several changes:
1. Polynomial addition with the same variables
2. C++/Python binding when getting the result from C++
3. Separated loop from MapBeamLine in two loops 
All these changes have improved the running times