NCOPair.h

Go to the documentation of this file.

// Copyright 2025 Julien Lamy, ICube, Université de Strasbourg-CNRS.
// Part of msl, distributed under the terms of the MIT license.
 
#ifndef _8c24ac6b_445b_442f_a45b_1f44e54f5720
#define _8c24ac6b_445b_442f_a45b_1f44e54f5720
 
#include <cmath>
 
#include <MrGlobalDefinitions/MrResult.h>
 
#include <MrMeasSrv/SeqIF/libRT/sFREQ_PHASE.h>
#include <MrMeasSrv/SeqIF/libRT/sGRAD_PULSE.h>
#include <MrMeasSrv/SeqIF/libRT/sREADOUT.h>
#include <MrMeasSrv/SeqIF/libRT/sSLICE_POS.h>
 
#include "msl/Vector.h"
 
namespace msl
{
 
class NCOPair
{
public:
    template <typename T>
    void setTime(T && event)
    {
        this->_set.setStartTime(event.getStartTime());
        // NOTE: sREADOUT has a non-integer duration, make sure the NCO always
        // ends after the readout
        this->_reset.setStartTime(
            event.getStartTime()+std::ceil(event.getDuration()));
    }
    
    template <typename T>
    void update(sSLICE_POS const & slice, T && pulse, double additionalPhase=0)
    {
        this->_set.prepSet(slice, pulse);
        this->_set.increasePhase(additionalPhase);
        this->_reset.prepNeg(slice, pulse);
        this->_reset.decreasePhase(additionalPhase);
    }
    
    void update(
        sSLICE_POS const & slice, sREADOUT const & adc,
        sGRAD_PULSE const & gradient, Vector3l const & point,
        double additionalPhase=0);
    
    NLS_STATUS run();
private:
    sFREQ_PHASE _set, _reset;
};
 
}
 
#endif // _8c24ac6b_445b_442f_a45b_1f44e54f5720