seqsat.h

/***************************************************************************
                          seqsat.h  -  description
                             -------------------
    begin                : Tue Oct 19 2004
    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 SEQSAT_H
#define SEQSAT_H
 
#include <odinseq/seqpulsar.h>
#include <odinseq/seqgradpulse.h>
 
 
enum satNucleus {fat=0,water};
 
 
 
class SeqPulsarSat : public SeqPulsar {
 
 public:
 
 
  SeqPulsarSat(const STD_string& object_label="unnamedSeqPulsarSat", satNucleus nuc=fat,float bandwidth=0.3);
 
  SeqPulsarSat(const STD_string& object_label,float bandwidth,
                           double freqoffset=0.0, float flipangle=90.0, float rel_filterwidth=0.3);
 
 
  SeqPulsarSat(const SeqPulsarSat& spg);
 
  SeqPulsarSat& operator = (const SeqPulsarSat& spg);
 
};
 
 
 
 
class SeqSat : public SeqObjList, public virtual SeqPulsInterface, public virtual SeqGradInterface {
 
  // virtual functions of SeqPulsInterface are marshalled to puls
  // virtual functions of SeqFreqChanInterface are marshalled to puls
 
 public:
  SeqSat(const STD_string& object_label="unnamedSeqSat",satNucleus nuc=fat,float bandwidth=0.3, unsigned int npulses=1);
 
 
  SeqSat(const SeqSat& spg);
 
  SeqSat& operator = (const SeqSat& spg);
 
 
  // overloading virtual function from SeqGradInterface
  SeqGradInterface& set_strength(float gradstrength);
  SeqGradInterface& invert_strength();
  float get_strength() const;
  fvector get_gradintegral() const;
  double get_gradduration() const;
  SeqGradInterface& set_gradrotmatrix(const RotMatrix& matrix);
 
 private:
  void build_seq();
 
  SeqPulsarSat puls;
  SeqGradConstPulse spoiler_read_pos;
  SeqGradConstPulse spoiler_slice_neg;
  SeqGradConstPulse spoiler_read_neg;
  SeqGradConstPulse spoiler_slice_pos;
  SeqGradConstPulse spoiler_phase_pos;
 
  unsigned int npulses_cache;
 
};
 
#endif