seqgradspiral.h

/***************************************************************************
                          seqgradspiral.h  -  description
                             -------------------
    begin                : Tue Aug 13 2002
    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 SEQGRADSPIRAL_H
#define SEQGRADSPIRAL_H
 
#include <tjutils/tjnumeric.h>
 
#include <odinseq/odinpulse.h>
#include <odinseq/seqgradchan.h>
#include <odinseq/seqgradwave.h>
#include <odinseq/seqgradconst.h>
#include <odinseq/seqgradchanparallel.h>
 
 
class SeqGradSpiral: public SeqGradChanParallel, public MinimizationFunction {
 
 public:
  SeqGradSpiral(const STD_string& object_label, LDRtrajectory& traj, double dt, float resolution, unsigned int sizeRadial, unsigned int numofSegments, bool inwards=false, bool optimize=false, const STD_string& nucleus="");
 
  SeqGradSpiral(const STD_string& object_label = "unnamedSeqGradSpiral" );
 
  SeqGradSpiral(const SeqGradSpiral& sgs);
 
  SeqGradSpiral& operator = (const SeqGradSpiral& sgs);
 
 
  unsigned int spiral_size() const {return denscomp.size();}
 
  double get_ramp_duration() const {return (get_gradduration()-spiral_dur);}
 
 
  SeqGradSpiral& set_predelay_duration(double dur);
 
  const fvector& get_denscomp() const {return denscomp;}
 
  fvector get_ktraj(direction channel) const;
 
 
 private:
  void common_init() {traj_cache=0;}
  void build_seq();
 
  float readout_npts() const;
 
  // implement virtual functions of MinimizationFunction to minimize readout duration
  float evaluate(const fvector& spirpar) const;
  unsigned int numof_fitpars() const {return 1;}
  
  SeqGradWave gx;
  SeqGradWave gy;
 
  SeqGradDelay gxdelay;
  SeqGradDelay gydelay;
 
  fvector kx;
  fvector ky;
 
  fvector denscomp;
 
  double spiral_dur;
  double predelay;
 
  // cache values for readout_npts/evaluate
  LDRtrajectory* traj_cache;
  double dt_cache;
  float resolution_cache;
  unsigned int sizeRadial_cache;
  float gamma_cache;
 
};
 
 
 
 
#endif