Opal maps

This MR will implement Map tracking for:

• Quadrupole
• Dipole
• Drifts
• Bends

TODO list:

• Code implementation (@ganz_p): ThickTracker, Hamiltonian, MapAnalyser
• Doxygen documentation (@ganz_p)
• Documentation (@ganz_p)
• Code cleaning (@frey_m)
• Regression tests (@snuverink_j)
• Examples (@ganz_p)
• Jupyter notebook (@ganz_p)
• Use beam line order

Remaining issues:

• Output dumps in local coordinates only
• Dipole until 1st order
• Tune calculation not stable
• Space charge calculations needs to be verified
• Fringe fields for Bends

src/Algorithms/ThickTracker.cpp

@@ -241,9 +241,11 @@ void ThickTracker::checkElementOrder_m() {
     for (FieldList::iterator it = elements.begin(); it != elements.end(); ++it) {
         double pos = it->getElement()->getElementPosition();
         // we have to take this length due to dipole --> arclength
-        if ( std::abs(pos - currentEnd) > threshold_m ) {
+        if ( currentEnd - pos > threshold_m ) {
+            std::cout << pos << " - " << currentEnd << std::endl;
             throw OpalException("ThickTracker::checkOrder_m()",
                                 std::string("Elements are not in ascending order or overlap.") +
+                                " Element Name: " + it->getElement()->getName() +
                                 " ELEMEDGE position: " + std::to_string(pos) +
                                 " Beamline end " + std::to_string(currentEnd) +
                                 " element length " + std::to_string(std::get<2>(*el))) ;