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

@@ -44,16 +44,17 @@ ThickTracker::ThickTracker(const Beamline &beamline,
   Tracker(beamline, reference, revBeam, revTrack),
   hamiltonian_m(1),
   itsOpalBeamline_m(beamline.getOrigin3D(), beamline.getCoordTransformationTo()),
-  pathLength_m(0.0),
-  zStop_m(),
+  zstart_m(0.0),
   dtCurrentTrack_m(0.0),
-  dtAllTracks_m(),
-  localTrackSteps_m(),
   truncOrder_m(1), // linear
   mapCreation_m(IpplTimings::getTimer("mapCreation")),
   mapCombination_m(IpplTimings::getTimer("mapCombination")),
   mapTracking_m(IpplTimings::getTimer("mapTracking"))
-{ }
+{
+  CoordinateSystemTrafo labToRef(beamline.getOrigin3D(),
+				 beamline.getCoordTransformationTo());
+  referenceToLabCSTrafo_m = labToRef.inverted();
+}
 
 
 ThickTracker::ThickTracker(const Beamline &beamline,
@@ -69,15 +70,10 @@ ThickTracker::ThickTracker(const Beamline &beamline,
                const int& truncOrder):
 	  Tracker(beamline, bunch, reference, revBeam, revTrack),
           hamiltonian_m(truncOrder),
-	  itsBeam_m(&beam),
   	  itsDataSink_m(&ds),
   	  itsOpalBeamline_m(beamline.getOrigin3D(), beamline.getCoordTransformationTo()),
-  	  pathLength_m(0.0),
   	  zstart_m(zstart),
-  	  zStop_m(),
   	  dtCurrentTrack_m(0.0),
-  	  dtAllTracks_m(),
-  	  localTrackSteps_m(),
   	  truncOrder_m(truncOrder),
       mapCreation_m(IpplTimings::getTimer("mapCreation")),
       mapCombination_m(IpplTimings::getTimer("mapCombination")),
@@ -516,13 +512,9 @@ void ThickTracker::execute() {
     bool consSC = itsBunch_m->hasFieldSolver();
     std::string fieldSolver = itsBunch_m->getFieldSolverType();
 
-    double minStepsize=0.01; // TODO better way to implement this value
-
-
     setconsiderSC(consSC);
     if ( fieldSolver == "SC" ){
         setconsiderSC(true); //TODO which way to set true?
-        setminStepsizeIfSC(minStepsize);
         msg << "SpaceCharges considered" << endl;
     } else {
         msg << "SpaceCharges neglected" << endl;