MapAnalyser.h

Go to the documentation of this file.

//
// Class: MapAnalyser
//   Organizes the function for a linear map analysis from
//   ThickTracker.
//   Transfer map -> tunes, symplecticity and stability
//   Sigma Matrix -> (not projected) beam emittance
//
// This class is in an unfinished state.
// For some dicussion see https://gitlab.psi.ch/OPAL/src/issues/464
// Some cleanup was done in https://gitlab.psi.ch/OPAL/src/merge_requests/294
//
// Copyright (c) 2018, Philippe Ganz, ETH Zürich
// All rights reserved
//
// Implemented as part of the Master thesis
// "s-based maps from TPS & Lie-Series applied to Proton-Therapy Gantries"
//
// This file is part of OPAL.
//
// OPAL is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
 
#ifndef MAP_ANALYSER_H
#define MAP_ANALYSER_H
 
#include "Classic/FixedAlgebra/FMatrix.h"
 
#include "Utility/IpplTimings.h"
 
#include <array>
#include <complex>
 
class MapAnalyser
{
public:
    typedef FMatrix<double, 6, 6> fMatrix_t;
    typedef FMatrix<std::complex<double>, 6, 6> cfMatrix_t;
 
    MapAnalyser();
 
    void linTAnalyze(const fMatrix_t& /*tMatrix*/) {};
 
 
    void linSigAnalyze(fMatrix_t& /*sigMatrix*/) {};
private:
 
    void eigenDecomp_m(const fMatrix_t& M, cfMatrix_t& eigenVal, cfMatrix_t& eigenVec, cfMatrix_t& invEigenVec);
 
    cfMatrix_t getBlockDiagonal_m(const fMatrix_t& M, cfMatrix_t& eigenVecM, cfMatrix_t& invEigenVecM);
 
    void printPhaseShift_m(fMatrix_t& Sigma, fMatrix_t tM, cfMatrix_t& oldN);
 
    void setNMatrix_m(fMatrix_t& M, cfMatrix_t& N, cfMatrix_t& invN);
 
 
    fMatrix_t createRotMatrix_m(std::array<double, 3> phi);
 
    fMatrix_t createSkewMatrix_m();
 
    fMatrix_t realPartOfMatrix_m(cfMatrix_t cM);
 
    fMatrix_t imagPartOfMatrix_m(cfMatrix_t cM);
 
    cfMatrix_t complexTypeMatrix_m(fMatrix_t M);
 
    cfMatrix_t invertMatrix_m(const cfMatrix_t &M);
 
    void rearrangeEigen_m(cfMatrix_t& /*EigenVal*/, cfMatrix_t& /*EigenVec*/) {};
 
    void normalizeEigen_m(cfMatrix_t& eigenVec, cfMatrix_t& invEigenVec);
 
private:
    IpplTimings::TimerRef mapAnalysis_m;
    IpplTimings::TimerRef bunchAnalysis_m;
};
 
#endif