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.h

@@ -532,18 +532,26 @@ inline void ThickTracker::visitSBend(const SBend &bend) {
 
     // Fringe fields currently not working
     //FIXME e1 not initialised
-    //double e1 = bend.getEntranceAngle();
-    //elements_m.push_back(std::make_tuple(hamiltonian_m.fringeField(e1, h),
-    //                                     1, 0.0));
-
+    //insert Entrance Fringefield
+    double e1 = bend.getEntranceAngle();
+    if (e1 > 1e-6){
+        elements_m.push_back(std::make_tuple(hamiltonian_m.fringeField(e1, h),
+                                         1, 0.0));
+    }
+    
+    //instert Dipole "body"
     elements_m.push_back(std::make_tuple(hamiltonian_m.sbend(gamma, h, k0),
                                          nSlices,
                                          arclength));
 
     //FIXME e2 not initialised
-    //double e2 = bend.getExitAngle();
-    //elements_m.push_back(std::make_tuple(hamiltonian_m.fringeField(e2, h),
-    //                                     1, 0.0));
+    //insert Exit Fringe field
+    double e2 = bend.getExitAngle();
+    if (e1 > 1e-6){
+    elements_m.push_back(std::make_tuple(hamiltonian_m.fringeField(e2, h),
+                                         1, 0.0));
+    }
+
 }