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

@@ -53,36 +53,36 @@ ThickTracker::ThickTracker(const Beamline &beamline,
   mapTracking_m(IpplTimings::getTimer("mapTracking"))
 {
   CoordinateSystemTrafo labToRef(beamline.getOrigin3D(),
-				 beamline.getCoordTransformationTo());
+                                 beamline.getCoordTransformationTo());
   referenceToLabCSTrafo_m = labToRef.inverted();
 }
 
 
 ThickTracker::ThickTracker(const Beamline &beamline,
-               PartBunchBase<double, 3> *bunch,
-               Beam &beam,
-			   DataSink &ds,
-			   const PartData &reference,
-               bool revBeam, bool revTrack,
-			   const std::vector<unsigned long long> &maxSteps,
-			   double zstart,
-			   const std::vector<double> &zstop,
-			   const std::vector<double> &dt,
-               const int& truncOrder):
-	  Tracker(beamline, bunch, reference, revBeam, revTrack),
-          hamiltonian_m(truncOrder),
-  	  itsDataSink_m(&ds),
-  	  itsOpalBeamline_m(beamline.getOrigin3D(), beamline.getCoordTransformationTo()),
-  	  zstart_m(zstart),
-  	  threshold_m(1.0e-6),
-  	  dtCurrentTrack_m(0.0),
-  	  truncOrder_m(truncOrder),
-      mapCreation_m(IpplTimings::getTimer("mapCreation")),
-      mapCombination_m(IpplTimings::getTimer("mapCombination")),
-      mapTracking_m(IpplTimings::getTimer("mapTracking"))
+                           PartBunchBase<double, 3> *bunch,
+                           Beam &beam,
+                           DataSink &ds,
+                           const PartData &reference,
+                           bool revBeam, bool revTrack,
+                           const std::vector<unsigned long long> &maxSteps,
+                           double zstart,
+                           const std::vector<double> &zstop,
+                           const std::vector<double> &dt,
+                           const int& truncOrder)
+    : Tracker(beamline, bunch, reference, revBeam, revTrack)
+    , hamiltonian_m(truncOrder)
+    , itsDataSink_m(&ds)
+    , itsOpalBeamline_m(beamline.getOrigin3D(), beamline.getCoordTransformationTo())
+    , zstart_m(zstart)
+    , threshold_m(1.0e-6)
+    , dtCurrentTrack_m(0.0)
+    , truncOrder_m(truncOrder)
+    , mapCreation_m(IpplTimings::getTimer("mapCreation"))
+    , mapCombination_m(IpplTimings::getTimer("mapCombination"))
+    , mapTracking_m(IpplTimings::getTimer("mapTracking"))
 {
   CoordinateSystemTrafo labToRef(beamline.getOrigin3D(),
-				 beamline.getCoordTransformationTo());
+                                 beamline.getCoordTransformationTo());
   referenceToLabCSTrafo_m = labToRef.inverted();
 
   for (std::vector<unsigned long long>::const_iterator it = maxSteps.begin(); it != maxSteps.end(); ++ it) {