pulsesequence() {
// Define Variables and Objects and Get Parameter Values
MPSEQ dec = getblew("blewH",0,0.0,0.0,0,1);
strncpy(dec.ch,"dec",3);
putCmd("chHblew='dec'\n");
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + getval("rd") + getval("ad") + at;
d.dutyoff = d1 + 4.0e-6;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(phRec,4,table4);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
_mpseqon(dec, phHhx);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_mpseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(phRec,4,table2);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
initval(getval("periods"),v2);
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Set Phase Tables
settable(phX90,4,table1);
settable(phRec,4,table2);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
xgate(1.0);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Apply a Rotorsync Delay
rgpulse(getval("pwX90"),phX90,0.0,0.0);
rotorsync(v2);
rgpulse(getval("pwX90"),phX90,0.0,0.0);
xgate(getval("xperiods"));
rgpulse(getval("pwX90"),phX90,0.0,0.0);
delay(10.0e-6);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence(){
//Define Variables and Get Parameter Values
double pwTune = getval("pwTune");
pwTune = pwTune*6.0;
at = pwTune*2.0;
char atval[MAXSTR];
sprintf(atval,"at = %f\n", at);
putCmd(atval);
int chTune = (int) getval("chTune");
if ((chTune < 1) || (chTune > 4)) {
abort_message("chTune(%d) must be between 1 and 4\n", chTune);
}
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Set Phase Tables
settable(phTune,4,table1);
settable(phRec,4,table2);
setreceiver(t2);
//Begin Sequence
obspwrf(getval("aTune"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2e-6); sp1off(); delay(2.0e-6);
//Begin Phase Detected Pulse
set4Tune(chTune,getval("gain"));
delay(1.0e-4);
ShapedXmtNAcquire("phtran",pwTune,phTune,6.0e-6,chTune);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aXecho = getval("aXecho");
double t1Xechoinit = getval("t1Xecho");
double pwXecho = getval("pwXecho");
double t2Xechoinit = getval("t2Xecho");
double t1Xecho = t1Xechoinit - pwXecho/2.0;
if (t1Xecho < 0.0) t1Xecho = 0.0;
double t2Xecho = t2Xechoinit - pwXecho/2.0 - getval("rd");
if (t2Xecho < 0.0) t2Xecho = 0.0;
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + pwXecho;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(phXecho,16,table4);
settable(phRec,8,table5);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
// sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
sp1on();
_cp_(hx,phHhx,phXhx);
sp1off();
// Begin Decoupling
_dseqon(dec);
// X Hahn Echo
txphase(phXecho);
obspwrf(aXecho);
delay(t1Xecho);
rgpulse(pwXecho,phXecho,0.0,0.0);
delay(t2Xecho);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ xdec = getdseq("X");
strncpy(xdec.t.ch,"obs",3); // These four statements assure
strncpy(xdec.s.ch,"obs",3); // that the X decoupling will
putCmd("chXtppm='obs'\n"); // be on X for either TPPM or
putCmd("chXspinal='obs'\n"); // SPINAL.
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = 3.0*getval("pwH90") + getval("pwHtilt") + getval("tHX");
d.dutyoff = d1 + 4.0e-6 + getval("tHmix");
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d.c3 = d.c3 + (!strcmp(xdec.seq,"tppm"));
d.c3 = d.c3 + ((!strcmp(xdec.seq,"tppm")) && (xdec.t.a > 0.0));
d.t3 = d2_;
d.c4 = d.c4 + (!strcmp(xdec.seq,"spinal"));
d.c4 = d.c4 + ((!strcmp(xdec.seq,"spinal")) && (xdec.s.a > 0.0));
d.t4 = d2_;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,16,table1);
settable(phHmix1,16,table2);
settable(phHmix2,4,table3);
settable(phHtilt,4,table4);
settable(phXhx,4,table5);
settable(phHhx,4,table6);
settable(phRec,4,table7);
if (phase1 == 2) tsadd(phH90,1,4);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H Preparation and tilt
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decunblank();
// Delay for 1H Wideline T2
_dseqon(xdec);
delay(d2);
_dseqoff(xdec);
// Mix period for Spin Diffison
decrgpulse(getval("pwH90"),phHmix1,0.0,0.0);
delay(getval("tHmix"));
decrgpulse(getval("pwH90"),phHmix2,0.0,0.0);
// Tilt Pulse and Ramped H to X Cross Polarization with LG Offset
decrgpulse(getval("pwHtilt"),phHtilt,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
obsblank(); _blank34();
_dseqon(dec);
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
//Define Variables and Objects and Get Parameter Values
double tHX3 = (getval("tHX"))/3.0; //Define MOIST CP in the Sequence
MPSEQ sd = getsammyd("smydH",0,0.0,0.0,0,1);
strncpy(sd.ch,"dec",3);
putCmd("chHsmyd='dec'\n");
MPSEQ so = getsammyo("smyoX",0,0.0,0.0,0,1);
strncpy(so.ch,"obs",3);
putCmd("chXsmyo='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
d2 = sd.nelem*sd.telem;
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("pwHlock")+ getval("tHX") + d2_;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phHlock,4,table2);
settable(phHcomp,4,table3);
settable(ph1Xhx,4,table4);
settable(ph2Xhx,4,table5);
settable(ph1Hhx,4,table6);
settable(ph2Hhx,4,table7);
settable(phHsmyd,4,table8);
settable(phXsmyo,4,table9);
settable(phHdec,4,table10);
settable(phRec,4,table11);
setreceiver(phRec);
// Begin Sequence
txphase(ph1Xhx);
decphase(phH90);
obspwrf(getval("aXhx"));
decpwrf(getval("aH90"));
obsunblank();
decunblank();
_unblank34();
delay(d1);
sp1on();
delay(2.0e-6);
sp1off();
delay(2.0e-6);
// H 90-degree Pulse
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
// Prelock with Compensation Pulse
decunblank();
decon();
if (getval("onHlock") > 0) {
decphase(phHlock);
delay(getval("pwHlock"));
decphase(phHcomp);
decpwrf(getval("aHcomp"));
delay(getval("pwHcomp"));
}
// H to X MOIST Cross Polarization
xmtron();
decphase(ph1Hhx);
decpwrf(getval("aHhx"));
delay(tHX3);
txphase(ph1Xhx);
decphase(ph1Hhx);
delay(tHX3);
txphase(ph2Xhx);
decphase(ph2Hhx);
delay(tHX3);
xmtroff();
decoff();
// SAMMY Spinlocks on X and H
_mpseqon(sd,phHsmyd);
_mpseqon(so,phXsmyo);
delay(d2);
_mpseqoff(sd);
_mpseqoff(so);
decphase(phHdec);
// Begin Acquisition
obsblank();
_blank34();
_dseqon(dec);
_decdoffset(getval("rd"),getval("ofHdec"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
decoffset(dof);
obsunblank();
decunblank();
_unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aXecho = getval("aXecho"); // define the echoX group in the sequence
double t1Xechoinit = getval("t1Xecho");
double pwXecho = getval("pwXecho");
double t2Xechoinit = getval("t2Xecho");
double t1Xecho = t1Xechoinit - pwXecho/2.0;
if (t1Xecho < 0.0) t1Xecho = 0.0;
double t2Xecho = t2Xechoinit - pwXecho/2.0 - getval("rd");
if (t2Xecho < 0.0) t2Xecho = 0.0;
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90") + pwXecho;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = t1Xecho + t2Xecho + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = t1Xecho + t2Xecho + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(phXecho,16,table2);
settable(phRec,8,table3);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H Decoupler on Before Echo
_dseqon(dec);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// X Solid-State Echo
txphase(phXecho);
obspwrf(aXecho);
delay(t1Xecho);
rgpulse(pwXecho,phXecho,0.0,0.0);
delay(t2Xecho);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
//======================================================
// Define Variables and Objects and Get Parameter Values
//======================================================
// --------------------------------
// Acquisition Decoupling
// -------------------------------
char Xseq[MAXSTR];
getstr("Xseq",Xseq);
DSEQ dec = getdseq("X");
strncpy(dec.t.ch,"dec",3);
putCmd("chXtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chXspinal='dec'\n");
//-------------------------------------
// Homonuclear Decoupling During Echo
//-------------------------------------
MPDEC homo1 = getmpdec("hdec1H",0,0.0,0.0,0,1);
strncpy(homo1.mps.ch,"obs",3);
putCmd("chHhdec1='obs'\n");
// --------------------
// H echo calculation
// --------------------
double t1Hecho = getval("t1Hecho") - getval("pwHecho")/2.0 -
((!strcmp(homo1.dm,"y"))?getval("pwHshort1")*2.:0.0);
if (t1Hecho < 0.0) t1Hecho = 0.0;
double t2Hecho = getval("t2Hecho") - getval("pwHecho")/2.0 -
((!strcmp(homo1.dm,"y"))?getval("pwHshort1")*2.:0.0) -
getval("rd")- getval("ad");
if (t2Hecho < 0.0) t2Hecho = 0.0;
double t1H_echo = 0.0;
double t2H_echo = 0.0;
double t1H_left = 0.0;
double t2H_left = 0.0;
if (!strcmp(homo1.dm,"y")) {
t2H_echo = homo1.mps.t*((int)(t2Hecho/homo1.mps.t));
t2H_left = t2Hecho - t2H_echo;
t1H_echo = t2H_echo;
t1H_left = t1Hecho - t1H_echo;
}
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
//----------------------
// Dutycycle Protection
//----------------------
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90");
d.dutyoff = d1 + 4.0e-6;
if (!strcmp(homo1.dm,"y"))
d.dutyon += t1H_echo + t2H_echo;
else
d.dutyoff += t1H_echo + t2H_echo;
d.c1 = d.c1 + (!strcmp(Xseq,"tppm"));
d.c1 = d.c1 + ((!strcmp(Xseq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(Xseq,"spinal"));
d.c2 = d.c2 + ((!strcmp(Xseq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
//------------------------
// Set Phase Tables
//-----------------------
settable(phH90,4,table1);
settable(phHecho,8,table2);
settable(phRec,4,table3);
setreceiver(phRec);
//=======================
// Begin Sequence
//=======================
txphase(phH90); decphase(zero);
obspwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
//------------------------
// H Direct Polarization
//------------------------
rgpulse(getval("pwH90"),phH90,0.0,0.0);
obsunblank(); decunblank(); _unblank34();
// -----------------------------
// H Hahn Echo
// -----------------------------
if (!strcmp(homo1.dm,"y")) {
delay (t1H_left);
if (getval("pwHshort1") > 0.0 ) {
obspwrf(getval("aHhdec1"));
rgpulse(getval("pwHshort1"),three,0.0,0.0);
obsunblank();
}
if (!strcmp(homo1.dm,"y")) _mpseqon(homo1.mps,zero);
delay(t1H_echo);
if (!strcmp(homo1.dm,"y")) _mpseqoff(homo1.mps);
if (getval("pwHshort1") > 0.0 ) {
obspwrf(getval("aHhdec1")); txphase(one);
rgpulse(getval("pwHshort1"),one,0.0,0.0);
obsunblank();
}
}
else delay(t1Hecho);
txphase(phHecho);
obspwrf(getval("aHecho"));
rgpulse(getval("pwHecho"),phHecho,0.0,0.0);
obsunblank();
if (!strcmp(homo1.dm,"y")) {
if (getval("pwHshort1") > 0.0 ) {
obspwrf(getval("aHhdec1"));
rgpulse(getval("pwHshort1"),three,0.0,0.0);
obsunblank();
}
if (!strcmp(homo1.dm,"y")) _mpseqon(homo1.mps,zero);
delay(t2H_echo);
if (!strcmp(homo1.dm,"y")) _mpseqoff(homo1.mps);
if(getval("pwHshort1")>0 ) {
obspwrf(getval("aHhdec1"));
rgpulse(getval("pwHshort1"),one,0.0,0.0);
obsunblank();
}
delay(t2H_left);
}
else delay(t2Hecho);
//====================
// Begin Acquisition
//====================
_dseqon(dec);
obsblank(); decblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pw1Hhytrap") + getval("pw2Hhytrap") + getval("tHX");
d.dutyoff = d1 + 4.0e-6 + getval("t1HYtrap") + getval("t2HYtrap");
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Hhytrap,4,table1);
settable(phYhytrap,4,table2);
settable(ph2Hhytrap,4,table3);
settable(phXhx,4,table4);
settable(phHhx,4,table5);
settable(phRec,4,table6);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(ph1Hhytrap); dec2phase(phYhytrap);
obspwrf(getval("aXhx")); decpwrf(getval("aHhytrap")); dec2pwrf(getval("aYhytrap"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// TRAPDOR on H with Y Modulation
decrgpulse(getval("pw1Hhytrap"),ph1Hhytrap,0.0,0.0);
decphase(ph2Hhytrap);
decunblank();
dec2on();
delay(getval("t1HYtrap"));
dec2off();
decrgpulse(getval("pw2Hhytrap"),ph2Hhytrap,0.0,0.0);
decphase(phHhx);
decunblank();
decphase(phHhx);
decpwrf(getval("aHhx"));
delay(getval("t2HYtrap"));
// H to X Cross Polarization
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Set the Maximum Dynamic Table and v-var Numbers
settablenumber(20);
setvvarnumber(30);
// Define Variables and Objects and Get Parameter Values
MPSEQ dumbo = getdumbogen("dumboX","dcf1X",0,0.0,0.0,0,1);
strncpy(dumbo.ch,"obs",3);
putCmd("chXdumbo='obs'\n");
MPSEQ c7 = getpostc7("c7X",0,0.0,0.0,0,1);
MPSEQ c7ref = getpostc7("c7X",c7.iSuper,c7.phAccum,c7.phInt,1,1);
strncpy(c7.ch,"obs",3);
putCmd("chXc7='obs'\n");
WMPA wdumbo = getwdumbogen("wdumboX","dcfX");
strncpy(wdumbo.ch,"obs",3);
putCmd("chXwdumbo='obs'\n");
double tXzfinit = getval("tXzf"); //Define the Z-filter delay in the sequence
double tXzf = tXzfinit - 5.0e-6 - wdumbo.r1;
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = c7.t + getval("pwXtilt") + d2_ + getval("pwXtilt") + c7ref.t + getval("pwX90") +
+ wdumbo.q*wdumbo.cycles*wdumbo.pw;
d.dutyoff = 4.0e-6 + d1 + tXzfinit + wdumbo.r2 + at - wdumbo.q*wdumbo.cycles*wdumbo.pw;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Xc7,4,table1);
settable(phXdumbo,4,table2);
settable(ph2Xc7,4,table3);
settable(phX90,16,table4);
settable(phXwdumbo,4,table5);
settable(phRec,16,table6);
settable(ph1Xtilt,4,table7);
settable(ph2Xtilt,4,table8);
// Set the Small-Angle Prep Phase
double obsstep = 360.0/(PSD*8192);
obsstepsize(obsstep);
int phfX90 = initphase(0.0, obsstep);
//Add STATES Quadrature Phase
if (phase1 == 2)
initval((45.0/obsstep),v1);
else
initval(0.0,v1);
initval((d2*c7.of[0]*360.0/obsstep),v2);
initval(0.0,v3);
obsstepsize(obsstep);
setreceiver(phRec);
// Begin Sequence
xmtrphase(v1); txphase(ph1Xc7);
obspwrf(getval("aXc7"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// C7 Recoupling of 2Q coherence
_mpseq(c7, ph1Xc7);
// F1 Evolution With DUMBO
xmtrphase(v3);
if (!getval("scXdcf1")){
obspwrf(getval("aX90"));
rgpulse(getval("pwXtilt"),ph1Xtilt,0.0,0.0);
}
obspwrf(getval("aXdumbo"));
obsunblank();
_mpseqon(dumbo,phXdumbo);
delay(d2);
_mpseqoff(dumbo);
if (!getval("scXdcf1")){
obspwrf(getval("aX90"));
rgpulse(getval("pwXtilt"),ph2Xtilt,0.0,0.0);
}
obspwrf(getval("aX90"));
obsunblank();
// C7 Transfer to 1Q Coherence
xmtrphase(v2);
_mpseq(c7ref, ph2Xc7);
// Z-filter Delay
delay(tXzf);
// Detection Pulse
txphase(phX90);
obspwrf(getval("aX90"));
startacq(5.0e-6);
rcvroff();
delay(wdumbo.r1);
rgpulse(getval("pwX90"), phX90, 0.0, 0.0);
obsunblank();
xmtrphase(v3);
delay(wdumbo.r2);
// Apply WPMLG Cycles During Acqusition
decblank(); _blank34();
_wdumbo(wdumbo,phXwdumbo);
endacq();
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Set the Maximum Dynamic Table and v-var Numbers
settablenumber(10);
setvvarnumber(30);
// Define Variables and Objects and Get Parameter Values
double aXprep = getval("aXprep");
double pwXprep = getval("pwXprep");
double phvXprep = getval("phXprep");
WMPA wsam = getwsamn("wsamX");
strncpy(wsam.ch,"obs",3);
putCmd("chXwsam='obs'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwXprep") + wsam.cycles1*wsam.cycles*wsam.pw;
d.dutyoff = d1 + 4.0e-6 + 5.0e-6 + wsam.r1 + wsam.r2 +
at - wsam.cycles1*wsam.cycles*wsam.pw;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phXprep,4,table1);
settable(phXwsam,4,table2);
settable(phRec,4,table3);
setreceiver(phRec);
// Set the Small-Angle Step
double obsstep = 360.0/(PSD*8192);
obsstepsize(obsstep);
int phfXprep = initphase(phvXprep,obsstep);
int phXzero = initphase(0.0,obsstep);
// Begin Sequence
xmtrphase(phfXprep); txphase(phXprep);
obspwrf(aXprep);
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Standard 90-degree prepX pulse for SAM
startacq(5.0e-6);
rcvroff();
delay(wsam.r1);
rgpulse(pwXprep, phXprep, 0.0, 0.0);
xmtrphase(phXzero);
delay(wsam.r2);
// Apply WSAM Cycles
decblank(); _blank34();
_wsamn(wsam, phXwsam);
endacq();
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hy = getcp("HY",0.0,0.0,0,1);
strncpy(hy.fr,"dec",3);
strncpy(hy.to,"dec2",4);
putCmd("frHY='dec'\n");
putCmd("toHY='dec2'\n");
CP yx = getcp("YX",0.0,0.0,0,1);
strncpy(yx.fr,"dec2",4);
strncpy(yx.to,"obs",3);
putCmd("frYX='dec2'\n");
putCmd("toYX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwY90") + getval("pwH90") + getval("tHY") + getval("tYX")
+ 2.0*getval("pwX90");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Create Phasetables
settable(phH90,4,table1);
settable(phHhy,4,table2);
settable(phY90,4,table3);
settable(phYhy,4,table4);
settable(phHyx,4,table5);
settable(phYyx,4,table6);
settable(phXyx,8,table7);
settable(phRec,8,table8);
// Begin Sequence
setreceiver(phRec);
txphase(phXyx); decphase(phH90); dec2phase(phY90);
obspwrf(getval("aXyx")); decpwrf(getval("aH90")); dec2pwrf(getval("aY90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to Y Cross Polarization with a Y Prepulse
dec2rgpulse(getval("pwY90"),phY90,0.0,0.0);
dec2phase(phYhy);
dec2pwrf(getval("aYyx"));
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhy);
_cp_(hy,phHhy,phYhy);
// Y to X Cross Polarization
decphase(phHyx); dec2phase(phYyx);
decpwrf(getval("aHyx"));
decunblank(); decon();
_cp_(yx,phYyx,phXyx);
decphase(phHhy);
dec2blank();
decoff();
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hy = getcp("HY",0.0,0.0,0,1);
strncpy(hy.fr,"dec",3);
strncpy(hy.to,"dec2",4);
putCmd("frHY='dec'\n");
putCmd("toHY='dec2'\n");
CP yx = getcp("YX",0.0,0.0,0,1);
strncpy(yx.fr,"dec2",4);
strncpy(yx.to,"obs",3);
putCmd("frYX='dec2'\n");
putCmd("toYX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
// Set Constant-time Period for d3.
if (d3_index == 0) d3_init = getval("d3");
double d3_ = (ni - 1)/sw1 + d3_init;
putCmd("d3acqret = %f\n",roundoff(d3_,12.5e-9));
putCmd("d3dwret = %f\n",roundoff(1.0/sw2,12.5e-9));
// Set Mixing Period to N Rotor Cycles
double taur,mix,srate;
mix = getval("tXmix");
srate = getval("srate");
taur = 0.0;
if (srate >= 500.0)
taur = roundoff((1.0/srate), 0.125e-6);
else {
printf("ABORT: Spin Rate (srate) must be greater than 500\n");
psg_abort(1);
}
mix = roundoff(mix,taur);
mix = mix - getval("pwX90");
if (mix < 0.0) mix = 0.0;
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwY90") + getval("pwH90") + getval("tHY") + getval("tYX")
+ 2.0*getval("pwX90") + mix;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = d2 + d3 + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = d2 + d3 + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Create Phasetables
settable(phH90,4,table1);
settable(phHhy,4,table2);
settable(phY90,4,table3);
settable(phYhy,4,table4);
settable(phHyx,4,table5);
settable(phYyx,4,table6);
settable(phXyx,8,table7);
settable(phXmix1,8,table8);
settable(phXmix2,8,table9);
settable(phRec,8,table10);
if (phase2 == 2)
tsadd(phXyx,1,4);
if (phase1 == 2)
tsadd(phYhy,1,4);
// Begin Sequence
setreceiver(phRec);
txphase(phXyx); decphase(phH90); dec2phase(phY90);
obspwrf(getval("aXyx")); decpwrf(getval("aH90")); dec2pwrf(getval("aY90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to Y Cross Polarization with a Y Prepulse
dec2rgpulse(getval("pwY90"),phY90,0.0,0.0);
dec2phase(phYhy);
dec2pwrf(getval("aYyx"));
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhy);
_cp_(hy,phHhy,phYhy);
// F1 Indirect Period For Y
_dseqon(dec);
delay(d2);
_dseqoff(dec);
// Y to X Cross Polarization
decphase(phHyx); dec2phase(phYyx);
decpwrf(getval("aHyx"));
decunblank(); decon();
_cp_(yx,phYyx,phXyx);
decphase(phHhy);
decoff();
// F2 Indirect Period for X
txphase(phXmix1);
obspwrf(getval("aX90"));
_dseqon(dec);
delay(d3);
_dseqoff(dec);
// RAD(DARR) Mixing For X
decpwrf(getval("aHmix"));
decunblank(); decon();
rgpulse(getval("pwX90"),phXmix1,0.0,0.0);
txphase(phXmix2);
obsunblank();
delay(mix);
rgpulse(getval("pwX90"),phXmix2,0.0,0.0);
decoff();
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double pwX90 = getval("pwX90");
double d22 = d2/2.0 - pwX90;
if (d22 < 0.0) d22 = 0.0;
MPSEQ fh = getfslg("fslgH",0,0.0,0.0,0,1);
strncpy(fh.ch,"dec",3);
putCmd("chHfslg='dec'\n");
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = 3.0*getval("pwHtilt") + d2_ + getval("pwH90") + getval("tHX");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phHtilt,4,table1);
settable(phHfslg,4,table2);
settable(phHtilt2,4,table3);
settable(phH90,4,table4);
settable(phHtilt3,4,table5);
settable(phXhx,4,table6);
settable(phHhx,4,table7);
settable(phRec,4,table8);
//Add STATES TPPI ("States with "FAD")
tsadd(phRec,2*d2_index,4);
if (phase1 == 2) {
tsadd(phHtilt3,2*d2_index+3,4);
tsadd(phHhx,2*d2_index+3,4);
}
else {
tsadd(phHtilt3,2*d2_index,4);
tsadd(phHhx,2*d2_index,4);
}
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phHtilt);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank();decunblank();_unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H Preparation with a Tilt Pulse
decrgpulse(getval("pwHtilt"),phHtilt,0.0,0.0);
// FSLG spinlock on H and Reverse Tilt to Zed
_mpseqon(fh,phHfslg);
delay(d22);
rgpulse(2.0*pwX90,zero,0.0,0.0);
txphase(phXhx);
delay(d22);
_mpseqoff(fh);
decpwrf(getval("aH90"));
decrgpulse(getval("pwHtilt"),phHtilt2,0.0,0.0);
// H 90 and Ramped H to X Cross Polarization with LG Offset
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decrgpulse(getval("pwHtilt"),phHtilt3,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
obsblank(); _blank34();
_dseqon(dec);
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
PBOXPULSE shca = getpboxpulse("shcaX",0,1);
strncpy(shca.ch,"obs",3);
putCmd("chXshca ='obs'\n");
PBOXPULSE shco = getpboxpulse("shcoX",0,1);
strncpy(shco.ch,"obs",3);
putCmd("chXshco ='obs'\n");
PBOXPULSE shcaco = combine_PBOXPULSE(shca,shco,0,1);
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ mix = getdseq("Hmix");
strncpy(mix.t.ch,"dec",3);
putCmd("chHmixtppm='dec'\n");
strncpy(mix.s.ch,"dec",3);
putCmd("chHmixspinal='dec'\n");
double shcacolen = (shcaco.pw + 2.0*shcaco.t2)/2.0;
double d22 = d2/2.0;
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + 2.0*shcaco.pw +
4.0*getval("pwX90");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = d2_ + getval("tZF") - getval("pwX90") + getval("rd") +
getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = d2_ + getval("tZF") - getval("pwX90") + getval("rd") +
getval("ad") + at;
d.c3 = d.c3 + (!strcmp(mix.seq,"tppm"));
d.c3 = d.c3 + ((!strcmp(mix.seq,"tppm")) && (mix.t.a > 0.0));
d.t3 = 2.0*getval("taua") + 2.0*getval("taub") + 2.0*shcaco.t1 +
2.0*shcaco.t2 - 4.0*shcacolen - 2.5*getval("pwX90") - 6.0e-6;
d.c4 = d.c4 + (!strcmp(mix.seq,"spinal"));
d.c4 = d.c4 + ((!strcmp(mix.seq,"spinal")) && (mix.s.a > 0.0));
d.t4 = 2.0*getval("taua") + 2.0*getval("taub") + 2.0*shcaco.t1 +
2.0*shcaco.t2 - 4.0*shcacolen - 2.5*getval("pwX90") - 6.0e-6;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Up 2D
int errval = (int) ((getval("taua") - shcacolen)*2.0*sw1);
if ((getval("taua") - ni/(2.0*sw1) - shcacolen) < 0.0) {
text_error("Error:ni is too large. Make ni equal to %d or less.\n",errval);
psg_abort(1);
}
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(ph1Xshcaco,4,table4);
settable(ph1X90,4,table5);
settable(ph2X90,4,table6);
settable(ph2Xshcaco,4,table7);
settable(ph3X90,4,table8);
settable(ph4X90,4,table9);
settable(phRec,4,table10);
setreceiver(phRec);
// States Acquisition
if (phase1 == 2)
tsadd(phXhx,3,4);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
decphase(zero);
// Begin F1 Refocused INEPT
_dseqon(mix);
txphase(ph1Xshcaco);
obspwrf(shcaco.a);
delay(getval("taua") - d22 - shcacolen);
_pboxpulse(shcaco,ph1Xshcaco);
obsblank();
obspower(getval("tpwr"));
delay(3.0e-6);
obsunblank();
obsunblank();
txphase(ph1X90);
obspwrf(getval("aX90"));
delay(getval("taua") + d22 - shcacolen - getval("pwX90") - 3.0e-6);
rgpulse(getval("pwX90"),ph1X90,0.0,0.0);
rgpulse(getval("pwX90"),ph2X90,0.0,0.0);
txphase(ph2Xshcaco);
obspwrf(shcaco.a);
obsunblank();
delay(getval("taub") - shcacolen - getval("pwX90"));
_pboxpulse(shcaco,ph2Xshcaco);
obsblank();
obspower(getval("tpwr"));
delay(3.0e-6);
obsunblank();
txphase(ph3X90);
obspwrf(getval("aX90"));
obsunblank();
delay(getval("taub") - shcacolen - getval("pwX90")/2.0 - 3.0e-6);
rgpulse(getval("pwX90"),ph3X90,0.0,0.0);
txphase(ph4X90);
obsunblank();
_dseqoff(mix);
_dseqon(dec);
delay(getval("tZF") - getval("pwX90"));
rgpulse(getval("pwX90"),ph4X90,0.0,0.0);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
MPSEQ fh = getfslg("fslgH",0,0.0,0.0,0,1);
strncpy(fh.ch,"dec",3);
putCmd("chHfslg='dec'\n");
MPSEQ fx = getfslg("fslgX",0,0.0,0.0,0,1);
strncpy(fx.ch,"obs",3);
putCmd("chXfslg='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + d2_;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(phHtilt,4,table4);
settable(phHfslg,4,table5);
settable(phXlock,4,table6);
settable(phRec,4,table7);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Tilt Pulse on H with Continued X Spinlock
xmtron(); decon();
decpwrf(getval("aH90")); obspwrf(getval("aXhx"));
decrgpulse(getval("pwHtilt"),phHtilt,0.0,0.0);
// FSLG Spinlocks on X and H
xmtron(); decon();
_mpseqon(fh,phHfslg); _mpseqon(fx,phXlock);
delay(d2);
_mpseqoff(fh); _mpseqoff(fx);
// Begin Acquisition
obsblank(); _blank34();
_dseqon(dec);
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
SHAPE p1 = getpulse("90H",0.0,0.0,1,0);
strncpy(p1.pars.ch,"dec",3);
putCmd("chH90='dec'\n");
p1.pars.array = disarry("xx", p1.pars.array);
p1 = update_shape(p1,0.0,0.0,1);
MPSEQ ph = getpmlgxmx("pmlgH",0,0.0,0.0,1,0);
strncpy(ph.ch,"dec",3);
putCmd("chHpmlg='dec'\n");
double pwHpmlg = getval("pwHpmlg");
ph.nelem = (int) (d2/(2.0*pwHpmlg) + 0.1);
ph.array = disarry("xx", ph.array);
ph = update_mpseq(ph,0,p1.pars.phAccum,p1.pars.phInt,1);
SHAPE p2 = getpulse("90H",0.0,0.0,2,0);
strncpy(p2.pars.ch,"dec",3);
putCmd("chH90='dec'\n");
p2.pars.array = disarry("xx", p2.pars.array);
p2 = update_shape(p2,ph.phAccum,ph.phInt,2);
double pwX180 = getval("pwX180");
double d22 = ph.t/2.0 - pwX180/2.0;
if (d22 < 0.0) d22 = 0.0;
// CP hx and DSEQ dec Return to the Reference Phase
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = p1.pars.t + d2_ + p2.pars.t + getval("pwH90") + getval("pwHtilt") +
getval("tHX");
d.dutyoff = d1 + 4.0e-6 + getval("tHmix");
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1H90,4,table1);
settable(phHpmlg,4,table2);
settable(ph2H90,4,table3);
settable(ph3H90,4,table4);
settable(phHtilt,4,table5);
settable(phXhx,4,table6);
settable(phHhx,4,table7);
settable(phRec,4,table8);
//Add STATES TPPI ("States with "FAD")
tsadd(phRec,2*d2_index,4);
if (phase1 == 2) {
tsadd(ph2H90,2*d2_index+3,4);
}
else {
tsadd(ph2H90,2*d2_index,4);
}
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(ph1H90);
obspwrf(getval("aX180")); decpwrf(getval("aH90"));
obsunblank();decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Offset H Preparation with a Tilt Pulse
_shape(p1,ph1H90);
// Offset SAMn Spinlock on H During F1 with Optional pwX180
_mpseqon(ph,phHpmlg);
delay(d22);
rgpulse(pwX180,zero,0.0,0.0);
obspwrf(getval("aX90"));
txphase(phXhx);
delay(d22);
_mpseqoff(ph);
// Offset 90-degree Pulse to Zed and Spin-Diffusion Mix
_shape(p2,ph2H90);
decpwrf(getval("aH90"));
delay(getval("tHmix"));
// H90, 35-degree Tilt and H-to-X Cross Polarization with LG Offset
decrgpulse(getval("pwH90"),ph3H90,0.0,0.0);
decunblank();
decrgpulse(getval("pwHtilt"),phHtilt,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin Acquisition
obsblank(); _blank34();
_dseqon(dec);
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
MPSEQ spc5 = getspc5("spc5X",0,0.0,0.0,0,1);
MPSEQ spc5ref = getspc5("spc5X",spc5.iSuper,spc5.phAccum,spc5.phInt,1,1);
strncpy(spc5.ch,"obs",3);
putCmd("chXspc5='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + getval("pwX90") +
spc5.t + spc5ref.t;
d.dutyoff = d1 + 4.0e-6 + 2.0*getval("tZF");
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = d2_ + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = d2_ + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Create Phasetables
settable(phH90,4,table1);
settable(phHhx,4,table2);
settable(phXhx,4,table3);
settable(phXmix1,4,table4);
settable(phXmix2,4,table5);
settable(phRec,4,table6);
setreceiver(phRec);
if (phase1 == 2)
tsadd(phXhx,1,4);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// F2 Indirect Period for X
obspwrf(getval("aX90"));
_dseqon(dec);
delay(d2);
_dseqoff(dec);
// Mixing with SPC5 Recoupling
rgpulse(getval("pwX90"),phXmix1,0.0,0.0);
obspwrf(getval("aXspc5"));
xmtrphase(v1); txphase(phXmix1);
delay(getval("tZF"));
decpwrf(getval("aHmix"));
decon();
_mpseq(spc5, phXmix1);
xmtrphase(v2); txphase(phXmix2);
_mpseq(spc5ref, phXmix2);
decoff();
obspwrf(getval("aX90"));
xmtrphase(zero); txphase(phXmix2);
delay(getval("tZF"));
rgpulse(getval("pwX90"),phXmix2,0.0,0.0);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
MPSEQ fh = getfslg("fslgH",0,0.0,0.0,0,1);
strncpy(fh.ch,"dec",3);
putCmd("chHfslg='dec'\n");
double tHXhmqc = getval("tHXhmqc"); //parameters for hmqcHX implemented
double pwHhxhmqc = getval("pwHhxhmqc"); //directly in the pulse sequence
double pmHhxhmqc = getval("pmHhxhmqc");
double pwXhxhmqc = getval("pwXhxhmqc");
double aXhxhmqc = getval("aXhxhmqc");
double aHhxhmqc = getval("aHhxhmqc");
double d2init = getval("d2");
d2init = d2init - pwXhxhmqc;
if (d2init < 0.0) d2init = 0.0;
double d22 = d2init/2.0;
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX")+ 2.0*tHXhmqc + d2_ +
4.0*pmHhxhmqc + 2.0*pwHhxhmqc;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,4,table1);
settable(phXhx,4,table2);
settable(phHhx,4,table3);
settable(ph1Hhxhmqc,4,table4);
settable(ph2Hhxhmqc,4,table5);
settable(ph3Hhxhmqc,4,table6);
settable(ph4Hhxhmqc,4,table7);
settable(ph5Hhxhmqc,4,table8);
settable(phXhxhmqc,32,table9);
settable(ph6Hhxhmqc,4,table10);
settable(ph7Hhxhmqc,8,table11);
settable(ph8Hhxhmqc,4,table12);
settable(ph9Hhxhmqc,4,table13);
settable(phRec,16,table14);
// Add STATES TPPI (States with FAD)
tsadd(ph3Hhxhmqc,2*d2_index,4);
tsadd(ph4Hhxhmqc,2*d2_index,4);
tsadd(ph5Hhxhmqc,2*d2_index,4);
tsadd(ph6Hhxhmqc,2*d2_index,4);
tsadd(phRec,2*d2_index,4);
if (phase1 == 2) {
tsadd(ph3Hhxhmqc,1,4);
tsadd(ph4Hhxhmqc,1,4);
tsadd(ph5Hhxhmqc,1,4);
tsadd(ph6Hhxhmqc,1,4);
}
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// Begin hmqcHX with fh (FSLG) Between Pulses
obspwrf(aXhxhmqc);
_mpseqon(fh,ph1Hhxhmqc); // First "tau" period for J evolution
delay(tHXhmqc);
_mpseqoff(fh);
decpwrf(aHhxhmqc);
decrgpulse(pmHhxhmqc, ph2Hhxhmqc, 0.0, 0.0);// Create HX double-quantum coherence
decrgpulse(pwHhxhmqc, ph3Hhxhmqc, 0.0, 0.0);
decrgpulse(pmHhxhmqc, ph4Hhxhmqc, 0.0, 0.0);
_mpseqon(fh,ph5Hhxhmqc); // Begin F1 evolution with FSLG
delay(d22);
rgpulse(pwXhxhmqc, phXhxhmqc, 0.0,0.0);
delay(d22);
_mpseqoff(fh); // End F1 evolution with FSLG
decpwrf(aHhxhmqc);
decrgpulse(pmHhxhmqc, ph6Hhxhmqc, 0.0, 0.0);// Refocus HX double quantum coherence
decrgpulse(pwHhxhmqc, ph7Hhxhmqc, 0.0, 0.0);
decrgpulse(pmHhxhmqc, ph8Hhxhmqc, 0.0, 0.0);
_mpseqon(fh,ph9Hhxhmqc); // Second "tau" period for J evolution
delay(tHXhmqc);
_mpseqoff(fh);
// Begin Acquisition
decphase(phHhx);
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aYxy8 = getval("aYxy8");
double pwYxy8 = getval("pwYxy8");
double nYxy8 = getval("nYxy8");
int cycles = (int) nYxy8/2.0;
nYxy8 = 2.0*cycles;
int counter = (int) (nYxy8 - 1.0);
initval((nYxy8 - 1.0),v8);
double onYxy8 = getval("onYxy8");
double srate = getval("srate");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ mix = getdseq("Hmix");
strncpy(mix.t.ch,"dec",3);
putCmd("chHmixtppm='mix'\n");
strncpy(mix.s.ch,"dec",3);
putCmd("chHmixspinal='mix'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90") + 4.0*nYxy8*pwYxy8 + getval("pwX180");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d.c3 = d.c3 + (!strcmp(mix.seq,"tppm"));
d.c3 = d.c3 + ((!strcmp(mix.seq,"tppm")) && (mix.t.a > 0.0));
d.t3 = 2.0*nYxy8*(1.0/srate - 2.0*pwYxy8) + 1.0/srate - getval("pwX180");
d.c4 = d.c4 + (!strcmp(mix.seq,"spinal"));
d.c4 = d.c4 + ((!strcmp(mix.seq,"spinal")) && (mix.s.a > 0.0));
d.t4 = 2.0*nYxy8*(1.0/srate - 2.0*pwYxy8) + 1.0/srate - getval("pwX180");
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(ph1Yxy8,8,table2);
settable(ph2Yxy8,4,table3);
settable(phX180,4,table4);
settable(phRec,4,table5);
if (counter < 0) tsadd(phRec,2,4);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// X Single Pulse
rgpulse(getval("pwX90"),phX90,0.0,0.0);
// xy8Y Period One
obspwrf(getval("aX180"));
txphase(phX180);
if (counter >= 0) {
_dseqon(mix);
delay(pwYxy8/2.0);
dec2pwrf(aYxy8);
sub(v1,v1,v1);
if (counter >= 1) {
if (counter > 1) loop(v8,v9);
getelem(ph1Yxy8,v1,v4);
incr(v1);
getelem(ph2Yxy8,ct,v2);
add(v4,v2,v2);
dec2phase(v2);
delay(0.5/srate - pwYxy8);
if (onYxy8 == 2)
dec2rgpulse(pwYxy8,v2,0.0,0.0);
else
delay(pwYxy8);
if (counter > 1) endloop(v9);
}
// X Refocussing Pulse
delay(0.5/srate - pwYxy8/2.0 - getval("pwX180")/2.0);
rgpulse(getval("pwX180"),phX180,0.0,0.0);
dec2pwrf(aYxy8);
delay(0.5/srate - pwYxy8/2.0 - getval("pwX180")/2.0);
// xy8Y Period Two
if (counter >= 1) {
if (counter > 1) loop(v8,v9);
if (onYxy8 == 2)
dec2rgpulse(pwYxy8,v2,0.0,0.0);
else
delay(pwYxy8);
getelem(ph1Yxy8,v1,v4);
incr(v1);
getelem(ph2Yxy8,ct,v2);
add(v4,v2,v2);
dec2phase(v2);
delay(0.5/srate - pwYxy8);
if (counter > 1) endloop(v9);
}
delay(pwYxy8/2.0);
_dseqoff(mix);
}
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Set the Maximum Dynamic Table Number
settablenumber(10);
setvvarnumber(30);
// Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
WMPA toss = gettoss5("tossX");
strncpy(toss.ch,"obs",3);
putCmd("chXtoss='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + 5.0*toss.pw;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = toss.rtau - 5.0*toss.pw + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = toss.rtau - 5.0*toss.pw + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,22,table1);
settable(phfXhx,11,table2);
settable(phXhx,44,table3);
settable(phHhx,4,table4);
settable(phHdec,4,table5);
settable(phXtoss,44,table6);
settable(phRec,44,table7);
setreceiver(phRec);
obsstepsize(360.0/(PSD*8192));
// Begin Sequence
xmtrphase(phfXhx); txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization - Shifted by -6COG11
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
decphase(phHdec);
// TOSS5 Sideband Suppression with included
// (-6,-5,-6,-5,-6)COG11 Cycle
_dseqon(dec);
_toss5(toss, phXtoss);
// Begin Acquisition with Quadrature Phase
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
//
// Set the Maximum Dynamic Table Number
//
settablenumber(10);
setvvarnumber(30);
//Define Variables and Objects and Get Parameter Values
CP hx = getcp("HX",0.0,0.0,0,1);
strncpy(hx.fr,"dec",3);
strncpy(hx.to,"obs",3);
putCmd("frHX='dec'\n");
putCmd("toHX='obs'\n");
WMPA cpmg = getcpmg("cpmgX");
strncpy(cpmg.ch,"obs",3);
putCmd("chXcpmg='obs'\n");
double aXecho = getval("aXecho"); // define the echoX group in the sequence
double t1Xechoinit = getval("t1Xecho");
double pwXecho = getval("pwXecho");
double t2Xechoinit = getval("t2Xecho");
double t1Xecho = t1Xechoinit - pwXecho/2.0 - getval("pwX90")/2.0;
if (t1Xecho < 0.0) t1Xecho = 0.0;
double t2Xecho = t2Xechoinit - pwXecho/2.0 - cpmg.r1 - cpmg.t2 - getval("ad");
if (t2Xecho < 0.0) t2Xecho = 0.0;
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHX") + pwXecho + (cpmg.cycles - 1)*cpmg.pw;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = t1Xecho + t2Xecho + getval("rd") + getval("ad") +
at - (cpmg.cycles - 1)*cpmg.pw;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = t1Xecho + t2Xecho + getval("rd") + getval("ad") +
at - (cpmg.cycles - 1)*cpmg.pw;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,64,table1);
settable(phXhx,64,table2);
settable(phHhx,64,table3);
settable(phXecho,64,table4);
settable(phXcpmg,64,table5);
settable(phRec,64,table6);
setreceiver(phRec);
// Begin Sequence
txphase(phXhx); decphase(phH90);
obspwrf(getval("aXhx")); decpwrf(getval("aH90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhx);
_cp_(hx,phHhx,phXhx);
// H Decoupling On
decphase(zero);
_dseqon(dec);
// X Hahn Echo
txphase(phXecho);
obspwrf(aXecho);
delay(t1Xecho);
rgpulse(pwXecho,phXecho,0.0,0.0);
delay(t2Xecho);
// Apply CPMG Cycles
obsblank(); _blank34();
delay(cpmg.r1);
startacq(getval("ad"));
_cpmg(cpmg,phXcpmg);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Set the Maximum Dynamic Table and v-var Numbers
settablenumber(10);
setvvarnumber(30);
// Define Variables and Objects and Get Parameter Values
double aXprep1 = getval("aXprep1"); // Define Tilted Pulses using "prep1X".
double pw1Xprep1 = getval("pw1Xprep1");
double pw2Xprep1 = getval("pw2Xprep1");
double phXprep1 = getval("phXprep1");
WMPA wpmlg = getwpmlg("wpmlgX");
strncpy(wpmlg.ch,"obs",3);
putCmd("chXwpmlg='obs'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pw1Xprep1") + getval("pw2Xprep1") + 2.0*wpmlg.q*wpmlg.cycles*wpmlg.pw;
d.dutyoff = d1 + 4.0e-6 + 5.0e-6 + wpmlg.r1 + wpmlg.r2 +
at - 2.0*wpmlg.q*wpmlg.cycles*wpmlg.pw;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Xprep1,4,table1);
settable(ph2Xprep1,4,table2);
settable(phXwpmlg,4,table3);
settable(phRec,4,table4);
setreceiver(phRec);
// Set the Small-Angle Step
double obsstep = 360.0/(PSD*8192);
obsstepsize(obsstep);
int phfXprep1 = initphase(phXprep1, obsstep);
int phXzero = initphase(0.0, obsstep);
// Begin Sequence
xmtrphase(phfXprep1); txphase(ph1Xprep1);
obspwrf(aXprep1);
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Tilted Preparation Pulse for FSLG or PMLG "prep1X"
startacq(5.0e-6);
rcvroff();
delay(wpmlg.r1);
rgpulse(pw1Xprep1, ph1Xprep1, 0.0, 0.0);
rgpulse(pw2Xprep1, ph2Xprep1, 0.0, 0.0);
xmtrphase(phXzero);
delay(wpmlg.r2);
// Apply WPMLG Cycles
decblank(); _blank34();
_wpmlg(wpmlg, phXwpmlg);
endacq();
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
PBOXPULSE shp1 = getpboxpulse("sft1A",0,1);
PBOXPULSE shp2 = getpboxpulse("sft2A",0,1);
PBOXPULSE shp3 = getpboxpulse("sft3A",0,1);
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwX90") + shp1.pw + shp2.pw + shp3.pw;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phX90,4,table1);
settable(phAsft1,4,table2);
settable(phAsft2,4,table3);
settable(phAsft3,4,table4);
settable(phRec,4,table5);
setreceiver(phRec);
// Begin Sequence
txphase(phX90); decphase(zero);
obspwrf(getval("aX90"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// X Direct Polarization
rgpulse(getval("pwX90"),phX90,0.0,0.0);
delay(20.0e-6);
// X Shaped Pulse
_pboxpulse(shp1, phAsft1);
delay(20.0e-6);
// X Simultaneous Shaped Pulse
_pboxsimpulse(shp1,shp2,phAsft1,phAsft2);
delay(20.0e-6);
// X 3-channel Simultaneous Shaped Pulse
delay(20.0e-6);
_pboxsim3pulse(shp1,shp2,shp3,phAsft1,phAsft2,phAsft3);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double pw1Xstmas = getval("pw1Xstmas");
double pw2Xstmas = getval("pw2Xstmas");
double tXzfselinit = getval("tXzfsel");
double tXzfsel = tXzfselinit - 3.0e-6;
if (tXzfsel < 0.0) tXzfsel = 0.0;
double d2init = getval("d2");
double d2 = d2init - pw1Xstmas/2.0 - pw2Xstmas/2.0;
if (d2 < 0.0) d2 = 0.0;
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pw1Xstmas") + getval("pw2Xstmas") + getval("pwXzfsel");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = d2_ + tXzfsel + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = d2_ + tXzfsel + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Xstmas,4,table1);
settable(ph2Xstmas,4,table2);
settable(phXzfsel,8,table3);
settable(phRec,8,table4);
if (phase1 == 2) {
tsadd(ph1Xstmas,1,4);
}
setreceiver(phRec);
// Begin Sequence
txphase(ph1Xstmas); decphase(zero);
obspower(getval("tpwr"));
obspwrf(getval("aXstmas"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H Decoupler on Before STMAS
_dseqon(dec);
// Two-Pulse STMAS
rgpulse(getval("pw1Xstmas"),ph1Xstmas,0.0,0.0);
txphase(ph2Xstmas);
delay(d2);
rgpulse(getval("pw2Xstmas"),ph2Xstmas,0.0,0.0);
// Z-filter Pulse
txphase(phXzfsel);
obsblank();
obspower(getval("dbXzfsel"));
obspwrf(getval("aXzfsel"));
delay(3.0e-6);
obsunblank();
delay(tXzfsel);
rgpulse(getval("pwXzfsel"),phXzfsel,0.0,0.0);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aXfam2 = getval("aXfam2");
double pw1Xfam2 = getval("pw1Xfam2");
double pw2Xfam2 = getval("pw2Xfam2");
double pw3Xfam2 = getval("pw3Xfam2");
double pw4Xfam2 = getval("pw4Xfam2");
double nXfam2 = getval("nXfam2");
initval(nXfam2,v4);
putCmd("pw2Xmqmas=pwXfam1"); // Sequence uses pwXfam1 and sets pw2Xmqmas
double d2init = getval("d2"); // Define the Split d2 in the Pulse Sequence
double ival = getval("ival");
double d20 = 1.0;
double d21 = 0.0;
double d22 = 0.0;
if (ival == 1.5) {
d20 = 9.0*d2init/16.0;
d21 = 7.0*d2init/16.0;
d22 = 0.0;
}
else if (ival == 2.5) {
d20 = 12.0*d2init/31.0;
d21 = 0.0*d2init/31.0;
d22 = 19.0*d2init/31.0;
}
else {
d20 = 1.0*d2init;
d21 = 0.0*d2init;
d22 = 0.0*d2init;
}
double tXechselinit = getval("tXechsel"); // Adjust the selective echo delay for the
double tXechsel = tXechselinit - 3.0e-6; // attenuator switch time.
if (tXechsel < 0.0) tXechsel = 0.0;
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
// Set Constant-time Period for d2.
if (d2_index == 0) d2_init = getval("d2");
double d2_ = (ni - 1)/sw1 + d2_init;
putCmd("d2acqret = %f\n",roundoff(d2_,12.5e-9));
putCmd("d2dwret = %f\n",roundoff(1.0/sw1,12.5e-9));
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pw1Xmqmas") + nXfam2*(pw1Xfam2 + pw2Xfam2 + pw3Xfam2 +pw4Xfam2) +
getval("pwXechsel");
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = d2_ + tXechselinit + getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = d2_ + tXechselinit + getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
if (phase1 == 0) {
settable(phf1Xmqmas,12,table1);
settable(ph1Xfam2,6,table2);
settable(ph2Xfam2,6,table3);
settable(phfXechsel,96,table4);
settable(phRec,48,table5);
}
else {
settable(phf1Xmqmas,6,table6);
settable(ph1Xfam2,6,table7);
settable(ph2Xfam2,6,table8);
settable(phfXechsel,48,table9);
settable(phRec,24,table10);
if (phase1 == 2) {
tsadd(phf1Xmqmas,30,360);
}
}
setreceiver(phRec);
obsstepsize(1.0);
// Begin Sequence
xmtrphase(phf1Xmqmas); decphase(zero);
obspower(getval("tpwr"));
obspwrf(getval("aXmqmas"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H Decoupler on Before MQMAS
_dseqon(dec);
// Two-Pulse MQMAS with DFS Conversion
rgpulse(getval("pw1Xmqmas"),zero,0.0,0.0);
xmtrphase(zero); txphase(ph1Xfam2);
obspwrf(aXfam2);
delay(d20);
// X FAM2 Pulse
loop(v4,v5);
xmtron();
delay(pw1Xfam2);
xmtroff();
txphase(ph2Xfam2);
delay(pw2Xfam2);
xmtron();
delay(pw3Xfam2);
xmtroff();
txphase(ph2Xfam2);
delay(pw4Xfam2);
endloop(v5);
// Tau Delay and Second Selective Echo Pulse
xmtrphase(phfXechsel);
obsblank();
obspower(getval("dbXechsel"));
obspwrf(getval("aXechsel"));
delay(3.0e-6);
obsunblank();
delay(d21 + tXechsel);
rgpulse(getval("pwXechsel"),zero,0.0,0.0);
delay(d22);
// Begin Acquisition
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Define Variables and Objects and Get Parameter Values
double aXhxpto2 = getval("aXhxpto2");
double pw1Xhxpto2 = getval("pw1Xhxpto2");
double pw2Xhxpto2 = getval("pw2Xhxpto2");
double t1HXpto2init = getval("t1HXpto2");
double tau1 = t1HXpto2init - pw1Xhxpto2/2.0;
double t2HXpto2init = getval("t2HXpto2");
double tau2 = t2HXpto2init - pw1Xhxpto2/2.0 - pw2Xhxpto2/2.0;
double t3HXpto2init = getval("t3HXpto2");
double tau3 = t3HXpto2init - pw2Xhxpto2/2.0;
MPSEQ r18 = getr1825("r18H",0,0.0,0.0,0,1);
MPSEQ r18ref = getr1825("r18H",r18.iSuper,r18.phAccum,r18.phInt,1,1);
strncpy(r18.ch,"dec",3);
strncpy(r18ref.ch,"dec",3);
putCmd("chHr18='dec'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = tau1 + tau2 + pw1Xhxpto2 + tau3 + pw2Xhxpto2;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(ph1Hhxpto2,4,table1);
settable(ph2Hhxpto2,4,table2);
settable(ph1Xhxpto2,4,table3);
settable(ph2Xhxpto2,4,table4);
settable(phHdec,4,table5);
settable(phRec,4,table6);
setreceiver(phRec);
// Begin Sequence
txphase(ph1Xhxpto2); decphase(ph1Hhxpto2);
obspwrf(aXhxpto2); decpwrf(r18.a);
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to X Cross Polarization with PRESTO1
_mpseqon(r18,ph1Hhxpto2);
delay(tau1);
rgpulse(pw1Xhxpto2/2.0,ph1Xhxpto2,0.0,0.0);
_mpseqoff(r18);
_mpseqon(r18ref,ph2Hhxpto2);
rgpulse(pw1Xhxpto2/2.0,ph1Xhxpto2,0.0,0.0);
delay(tau2);
_mpseqoff(r18ref);
decphase(zero);
_dseqon(dec);
rgpulse(pw2Xhxpto2,ph2Xhxpto2,0.0,0.0);
delay(tau3);
decphase(zero);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
pulsesequence() {
// Define Variables and Objects and Get Parameter Values
CP hy = getcp("HY",0.0,0.0,0,1);
strncpy(hy.fr,"dec",3);
strncpy(hy.to,"dec2",4);
putCmd("frHY='dec'\n");
putCmd("toHY='dec2'\n");
GP inept = getinept("ineptYX");
strncpy(inept.ch1,"dec2",4);
strncpy(inept.ch2,"obs",3);
putCmd("ch1YXinept='dec2'\n");
putCmd("ch2YXinept='obs'\n");
DSEQ dec = getdseq("H");
strncpy(dec.t.ch,"dec",3);
putCmd("chHtppm='dec'\n");
strncpy(dec.s.ch,"dec",3);
putCmd("chHspinal='dec'\n");
DSEQ mix = getdseq("Hmix");
strncpy(mix.t.ch,"dec",3);
putCmd("chHmixtppm='dec'\n");
strncpy(mix.s.ch,"dec",3);
putCmd("chHmixspinal='dec'\n");
// Dutycycle Protection
double simpw1 = inept.pw1;
if (inept.pw2 > inept.pw1) simpw1 = inept.pw2;
double simpw2 = inept.pw3;
if (inept.pw4 > inept.pw3) simpw2 = inept.pw4;
DUTY d = init_dutycycle();
d.dutyon = getval("pwH90") + getval("tHY") + 2.0*simpw1 + 2.0*simpw2;
d.dutyoff = d1 + 4.0e-6;
d.c1 = d.c1 + (!strcmp(dec.seq,"tppm"));
d.c1 = d.c1 + ((!strcmp(dec.seq,"tppm")) && (dec.t.a > 0.0));
d.t1 = inept.t1 + inept.t2 + inept.t3 + inept.t4 +
getval("rd") + getval("ad") + at;
d.c2 = d.c2 + (!strcmp(dec.seq,"spinal"));
d.c2 = d.c2 + ((!strcmp(dec.seq,"spinal")) && (dec.s.a > 0.0));
d.t2 = inept.t1 + inept.t2 + inept.t3 + inept.t4 +
getval("rd") + getval("ad") + at;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phH90,16,table1);
settable(phHhy,4,table2);
settable(phYhy,4,table3);
settable(ph1Yyxinept,4,table4);
settable(ph1Xyxinept,4,table5);
settable(ph2Yyxinept,4,table6);
settable(ph2Xyxinept,16,table7);
settable(ph3Yyxinept,8,table8);
settable(ph3Xyxinept,4,table9);
settable(phRec,8,table10);
setreceiver(phRec);
// Begin Sequence
txphase(ph1Xyxinept); decphase(phH90); dec2phase(phYhy);
obspwrf(getval("aXyxinept")); decpwrf(getval("aH90")); dec2pwrf(getval("aYhy"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// H to Y Cross Polarization
decrgpulse(getval("pwH90"),phH90,0.0,0.0);
decphase(phHhy);
_cp_(hy,phHhy,phYhy);
decphase(zero);
// INEPT Transfer from Y to X
_dseqon(mix);
_ineptref(inept,ph1Yyxinept,ph1Xyxinept,ph2Yyxinept,ph2Xyxinept,ph3Yyxinept,ph3Xyxinept);
_dseqoff(mix);
// Begin Acquisition
_dseqon(dec);
obsblank(); _blank34();
delay(getval("rd"));
startacq(getval("ad"));
acquire(np, 1/sw);
endacq();
_dseqoff(dec);
obsunblank(); decunblank(); _unblank34();
}
void pulsesequence() {
// Set the Maximum Dynamic Table Number
settablenumber(10);
setvvarnumber(30);
// Define Variables and Objects and Get Parameter Values
WMPA xx = getxx("xxX");
strncpy(xx.ch,"obs",3);
putCmd("chXxx='obs'\n");
//--------------------------------------
// Copy Current Parameters to Processed
//-------------------------------------
putCmd("groupcopy('current','processed','acquisition')");
// Dutycycle Protection
DUTY d = init_dutycycle();
d.dutyon = getval("pwXprep") + 2.0*xx.cycles*xx.pw;
d.dutyoff = d1 + 4.0e-6 + 5.0e-6 + xx.r1 + xx.r2 +
at - 2.0*xx.cycles*xx.pw;
d = update_dutycycle(d);
abort_dutycycle(d,10.0);
// Set Phase Tables
settable(phXprep,4,table1);
settable(phXxx,4,table2);
settable(phRec,4,table3);
setreceiver(phRec);
// Set the Small-Angle Prep Phase
double obsstep = 360.0/(PSD*8192);
obsstepsize(obsstep);
int phfXprep = initphase(getval("phXprep"), obsstep);
int phXzero = initphase(0.0, obsstep);
// Begin Sequence
xmtrphase(phfXprep); txphase(phXprep);
obspwrf(getval("aXprep"));
obsunblank(); decunblank(); _unblank34();
delay(d1);
sp1on(); delay(2.0e-6); sp1off(); delay(2.0e-6);
// Preparation Pulse with Initial Point
startacq(5.0e-6);
rcvroff();
delay(xx.r1);
rgpulse(getval("pwXprep"), phXprep, 0.0, 0.0);
xmtrphase(phXzero);
// Apply Semi-windowless WHH4 Cycles
decblank(); _blank34();
_xx(xx, phXxx);
endacq();
obsunblank(); decunblank(); _unblank34();
}