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

@@ -180,18 +180,6 @@ void ThickTracker::changeDT() {
 }
 
 
-//TODO verify Space Charge Maps
-/**returns a SC map
- */
-ThickTracker::map_t ThickTracker::getSpaceChargeMap(double& length){
-
-    series_t H;
-
-//     hamiltonian::spaceCharge(H, itsBunch_m);
-    return ExpMap(- H * length, 1 );
-}
-
-
 //TODO complete and test fringe fields
 void ThickTracker::insertFringeField(SBend* pSBend, lstruct_t& mBL,
         double& beta0, double& gamma0, double& P0, double& q, std::array<double,2>& paramFringe, std::string e){
@@ -272,19 +260,6 @@ void ThickTracker::execute() {
 
     prepareSections();
     
-    
-    bool consSC = itsBunch_m->hasFieldSolver();
-    std::string fieldSolver = itsBunch_m->getFieldSolverType();
-
-    setconsiderSC(consSC);
-    if ( fieldSolver == "SC" ){
-        setconsiderSC(true); //TODO which way to set true?
-        msg << "SpaceCharges considered" << endl;
-    } else {
-        msg << "SpaceCharges neglected" << endl;
-    }
-
-
     msg << "Truncation order: " << this->truncOrder_m << endl;