seqgradecho.h

/***************************************************************************
                          seqgradecho.h  -  description
                             -------------------
    begin                : Mon Feb 17 2003
    copyright            : (C) 2000-2021 by Thies H. Jochimsen
    email                : thies@jochimsen.de
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program 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 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#ifndef SEQGRADECHO_H
#define SEQGRADECHO_H
 
#include <odinseq/seqgradpulse.h>
#include <odinseq/seqgradphase.h>
#include <odinseq/seqacqread.h>
#include <odinseq/seqacqdeph.h>
#include <odinseq/seqpulsar.h>
#include <odinseq/seqsimvec.h>
 
 
class SeqGradEcho : public SeqObjList, public virtual SeqGradInterface, public virtual SeqAcqInterface {
 
  // virtual functions of SeqFreqChanInterface are marhshalled to acqread
 
 public:
 
  SeqGradEcho(const STD_string& object_label, SeqPulsar& exc,
              double sweepwidth, unsigned int readnpts, float FOVread,
              unsigned int phasenpts, float FOVphase, encodingScheme scheme=linearEncoding,
              reorderScheme reorder=noReorder, unsigned int nsegments=1, unsigned int reduction=1, unsigned int acl_bands=DEFAULT_ACL_BANDS,
              bool balanced=false, float partial_fourier_phase=0.0, float partial_fourier_read=0.0, bool partial_fourier_read_at_end=false,
              float os_factor=1.0, const STD_string& nucleus="");
 
 
  SeqGradEcho(const STD_string& object_label,
              unsigned int readnpts, float FOVread, unsigned int phasenpts, float FOVphase, unsigned int slicenpts, float FOVslice,
              SeqPulsar& exc, double sweepwidth,
              unsigned int reduction=1, unsigned int acl_bands=DEFAULT_ACL_BANDS,
              bool balanced=false, float partial_fourier_phase=0.0, float partial_fourier_phase3d=0.0,
              float partial_fourier_read=0.0, bool partial_fourier_read_at_end=false,
              float os_factor=1.0, const STD_string& nucleus="");
 
 
  SeqGradEcho(const STD_string& object_label="unnamedSeqGradEcho");
 
  SeqGradEcho(const SeqGradEcho& sge);
 
  SeqGradEcho& operator = (const SeqGradEcho& sge);
 
  SeqVector& get_pe_vector();
 
  SeqGradEcho& set_pe_reorder_scheme(reorderScheme scheme,unsigned int nsegments);
 
  const SeqVector& get_pe_reorder_vector() const;
 
  SeqVector& get_pe3d_vector();
 
  const SeqVector& get_exc_vector() const {return *pulsptr.get_handled();}
 
  SeqGradEcho& set_freq_reorder_scheme(reorderScheme scheme,unsigned int nsegments=1) {(*pulsptr.get_handled()).set_reorder_scheme(scheme,nsegments); return *this;}
 
 
 
  const SeqVector& get_freq_reorder_vector() const {return (*pulsptr.get_handled()).get_reorder_vector();}
 
 
  double get_echo_time() const;
 
 
  SeqGradEcho& set_midpart(const SeqObjBase& soa);
 
 
  // implementing virtual functions of SeqGradInterface
  SeqGradInterface& set_strength(float gradstrength) {return *this;}  // nothing useful
  float get_strength() const {return 0.0;}  // nothing useful
  SeqGradInterface& invert_strength();
  double get_gradduration() const {return get_duration();}
  SeqGradInterface& set_gradrotmatrix(const RotMatrix& matrix);
  fvector get_gradintegral() const;
 
 
  // marshalling virtual functions of SeqAcqInterface to acqread, except:
  double get_acquisition_center() const {return get_preacq()+acqread.get_acquisition_center();}
  double get_acquisition_start() const {return get_preacq()+acqread.get_acquisition_start();}
  SeqAcqInterface& set_template_type(templateType type);
 
 
 
 
 private:
  double get_preacq() const;
 
  void common_init(const STD_string& objlabel);
 
  void build_seq();
 
  // hide these functions from the user
  SeqAcqInterface& set_pe_vector(const SeqVector& pevec) {return *this;}
 
 
  Handler<SeqPulsNdim*> pulsptr;
  SeqPulsarReph pls_reph;
  SeqGradVector phase;
  SeqGradVector phase3d;
  SeqGradVector phase_rew;
  SeqGradVector phase3d_rew;
  SeqSimultanVector phasesim;
  SeqSimultanVector phasesim3d;
  SeqSimultanVector phasereordsim;
  SeqAcqRead acqread;
  SeqGradConst readdeph;
  SeqParallel postexcpart;
  SeqParallel postacqpart;
  SeqObjList midpart;
 
  geometryMode mode;
  bool balanced_grads;
};
 
#endif