seqacqspiral.h

/***************************************************************************
                          seqacqspiral.h  -  description
                             -------------------
    begin                : Tue Nov 11 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 SEQACQSPIRAL_H
#define SEQACQSPIRAL_H
 
#include <odinseq/seqlist.h>
#include <odinseq/seqparallel.h>
#include <odinseq/seqacq.h>
#include <odinseq/seqdelay.h>
#include <odinseq/seqgradspiral.h>
#include <odinseq/seqrotmatrixvector.h>
#include <odinseq/seqgradtrapez.h>
 
class SeqAcqSpiral : public virtual SeqAcqInterface, public SeqObjList  {
 
  // virtual functions of SeqFreqChanInterface are marhshalled to acq
 
 public:
 
 
  SeqAcqSpiral(const STD_string& object_label,double sweepwidth, float fov,
             unsigned int sizeRadial, unsigned int numofSegments, LDRtrajectory&  traj,
             bool inout=false, bool optimize=false, const STD_string& nucleus="",const dvector& phaselist=0);
 
 
  SeqAcqSpiral(const SeqAcqSpiral& sas);
 
 
 
  SeqAcqSpiral(const STD_string& object_label="unnamedSeqAcqSpiral");
 
  SeqAcqSpiral& operator = (const SeqAcqSpiral& sas);
 
 
  const SeqVector& get_segment_vector() const {return rotvec;}
 
 
  fvector get_ktraj(unsigned int iseg, direction channel) const;
 
  fvector get_denscomp() const;
 
 
  // marshalling virtual functions of SeqAcqInterface to acq, except:
  double get_acquisition_center() const {return preacq.get_duration()+acq.get_acquisition_center();}
  double get_acquisition_start() const {return preacq.get_duration()+acq.get_acquisition_start();}
  SeqAcqInterface& set_sweepwidth(double sw, float os_factor); // produce warning
 
 
 private:
 
  // overwriting virtual functions from SeqClass
  bool prep();
 
  void common_init();
 
  void build_seq();
 
  SeqParallel par;
 
  SeqGradSpiral spirgrad_in;
  SeqGradSpiral spirgrad_out;
 
  SeqDelay preacq;
  SeqAcq acq;
 
  SeqGradTrapezParallel gbalance;
 
  SeqRotMatrixVector rotvec;
 
  bool inout_traj;
};
 
#endif