seqdiffweight.h

/***************************************************************************
                          seqdiffweight.h  -  description
                             -------------------
    begin                : Wed Feb 12 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 SEQDIFFWEIGHT_H
#define SEQDIFFWEIGHT_H
 
#include <odinseq/seqparallel.h>
#include <odinseq/seqlist.h>
#include <odinseq/seqgrad.h>
#include <odinseq/seqgradpulse.h>
#include <odinseq/seqgradconst.h>
#include <odinseq/seqsimvec.h>
 
class SeqDiffWeight : public SeqObjList, public SeqSimultanVector, public virtual SeqGradInterface {
 
 public:
 
  SeqDiffWeight(const STD_string& object_label, const fvector& bvals, float maxgradstrength,
                    const SeqObjBase& midpart, direction chan, bool stejskalTanner=false, const STD_string& nucleus = "" );
 
 
  SeqDiffWeight(const STD_string& object_label, unsigned int ndir, const fvector& bvals, float maxgradstrength,
                    const SeqObjBase& midpart, unsigned int baseline_rep=0, bool stejskalTanner=false, const STD_string& nucleus = "" );
 
 
  SeqDiffWeight(const STD_string& object_label = "unnamedSeqDiffWeight" );
 
  SeqDiffWeight(const SeqDiffWeight& sgdw) {SeqDiffWeight::operator = (sgdw);}
 
  SeqDiffWeight& operator = (const SeqDiffWeight& sgdw);
 
  const SeqGradInterface& get_grad1() const {return par1;}
 
  const SeqGradInterface& get_grad2() const {return par2;}
 
  double get_grad1_duration() const {return par1.get_duration();}
 
  double get_grad2_duration() const {return par2.get_duration();}
 
  double get_midpart_duration() const {return middle_part.get_duration();}
 
  const darray& get_b_vectors() const {return b_vectors_cache;}
 
 
 
  // 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);
 
 
  // final overrider of virtual functions from SeqClass
  void clear_container() {SeqObjList::clear(); SeqSimultanVector::clear();}
 
 
  // implemented virtual functions from SeqVector
  bool is_qualvector() const {return false;}
 
 private:
 
  void build_seq();
 
 
  SeqGradVectorPulse pfg1[n_directions];
  SeqGradVectorPulse pfg2[n_directions];
 
 
  SeqParallel par1;
  SeqParallel par2;
 
  SeqObjList middle_part;
 
  darray b_vectors_cache;
 
};
 
 
 
class SeqDiffWeightFlowComp : public SeqGradChanList, public SeqSimultanVector {
 
 public:
 
  SeqDiffWeightFlowComp(const STD_string& object_label, const fvector& bvals, float maxgradstrength,
                    direction chan, double stimdelay=1.0, const STD_string& nucleus = "" );
 
  SeqDiffWeightFlowComp(const STD_string& object_label = "unnamedSeqDiffWeightFlowComp" );
 
  SeqDiffWeightFlowComp(const SeqDiffWeightFlowComp& sgdwfc) {SeqDiffWeightFlowComp::operator = (sgdwfc);}
 
  SeqDiffWeightFlowComp& operator = (const SeqDiffWeightFlowComp& sgdwfc);
 
 
  // final overrider of virtual functions from SeqClass
  void clear_container() {SeqGradChanList::clear(); SeqSimultanVector::clear();}
 
 
  // implemented virtual functions from SeqVector
  bool is_qualvector() const {return false;}
 
 private:
 
  void build_seq();
 
  SeqGradVectorPulse pfg1;
  SeqGradVectorPulse pfg2;
  SeqGradVectorPulse pfg3;
  SeqGradDelay delay;
};
 
 
 
 
#endif