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();
}
pulsesequence()
{
char
autocal[MAXSTR],
shname1[MAXSTR],
shname2[MAXSTR],
ipap_flg[MAXSTR],
ab_flg[MAXSTR],
f1180[MAXSTR],
SE[MAXSTR], /*Use Sensitivity Enhancement */
c13refoc[MAXSTR],
refpat[MAXSTR], /* pulse shape pattern refocus. pulse*/
hetsofast_flg[MAXSTR],
grad3_flg[MAXSTR]; /*gradient flag */
int
t1_counter,
phase;
double d2_init=0.0,
pwS,pwS1,pwS2,
tpwrsf = getval("tpwrsf"),
adjust = getval("adjust"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gstab = getval("gstab"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gt3 = getval("gt3"),
shlvl1 = getval("shlvl1"),
shlvl2 = getval("shlvl2"),
shdmf2 = getval("shdmf2"),
shbw = getval("shbw"),
shpw1 = getval("shpw1"),
shpw2 = getval("shpw2"),
shpw3 = 0.0,
shofs = getval("shofs"),
flipangle = getval("flipangle"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
dpwr2 = getval("dpwr2"),
d2 = getval("d2"),
tau1,
taunh = 1.0/(2.0*getval("JNH"));
void compo_pulse();
void compo1_pulse();
void makeshape_pc9();
void makeshape_refoc();
void makeshape_ndec();
getstr("autocal",autocal);
getstr("f1180",f1180);
getstr("c13refoc",c13refoc);
getstr("hetsofast_flg",hetsofast_flg);
getstr("refpat", refpat); /* pulse pattern refocussing pulse */
getstr("ipap_flg",ipap_flg);
getstr("grad3_flg",grad3_flg);
getstr("ab_flg",ab_flg);
getstr("SE",SE);
getstr("shname1",shname1);
getstr("shname2",shname2);
pwS = c_shapedpw("sech",200.0,0.0,zero, 0.0, 0.0);
pwS1 = hn_simshapedpw(refpat,shbw,shofs-4.8,"isnob3",50.0,0.0, zero, zero, 0.0, 0.0);
pwS2 = h_shapedpw(refpat,shbw,3.5,zero, 0.0, 0.0);
if (hetsofast_flg[0] == 'a')
shpw3 = h_shapedpw("isnob5",4.0,-3.0,two, 2.0e-6, 2.0e-6);
if (hetsofast_flg[0] == 'b')
shpw3 = h_shapedpw("gaus180",0.015,0.0,two, 2.0e-6, 2.0e-6);
phase = (int) (getval("phase") + 0.5);
settable(t1,2,phi1);
settable(t2,2,phi2);
if (autocal[0] == 'y')
{
(void) makeshape_pc9(flipangle, shbw, shofs); /*create pc9 pulse*/
sh90 = getRsh("hn_pc9");
shpw1 = sh90.pw;
shlvl1 = sh90.pwr;
sprintf(shname1,"hn_pc9");
(void) makeshape_refoc(refpat, shbw, shofs); /* create refocussing pulse */
shref = getDsh("hn_refoc");
shpw2 = shref.pw;
shlvl2 = shref.pwr;
shdmf2 = shref.dmf;
sprintf(shname2,"hn_refoc");
if (dmm2[2] == 'p') /* waveform capability present on channel 3 */
{
(void) makeshape_ndec(); /* create n15 wurst decoupling */
shdec = getDsh("hncompdec");
dpwr2 = shdec.pwr;
dmf2 = shdec.dmf;
dres2 = shdec.dres;
sprintf(dseq2,"hncompdec");
}
}
if (tpwrsf <4095.0) shlvl2 = shlvl2+6.0;
if (phase == 1) ;
if (phase == 2) tsadd(t1,1,4);
tau1 = d2;
if((f1180[A] == 'y') && (ni > 1.0))
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1;
if (f1180[0] == 'y') tau1 = tau1-pwN*4.0/3.0;
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t2,2,4); }
status(A);
decoffset(dof);
dec2power(pwNlvl);
dec2offset(dof2);
obsoffset(tof);
obspower(tpwr);
/**********************************************/
if (hetsofast_flg[0] != 'n')
{
if (hetsofast_flg[0] == 'a')
h_shapedpulse("isnob5",4.0,-3.0,two, 2.0e-6, 2.0e-6);
if (hetsofast_flg[0] == 'b')
h_shapedpulse("gaus180",0.015,0.0,two, 2.0e-6, 2.0e-6);
zgradpulse(gzlvl2, gt2);
delay(gstab);
delay(d1-gt2-shpw3);
lk_hold();
}
else
{
zgradpulse(gzlvl2, gt2);
delay(gstab);
delay(d1-gt2);
lk_hold();
}
obspower(shlvl1);
shaped_pulse(shname1,shpw1,zero,2.0e-4,2.0e-6);
if (SE[0] == 'y')
{
if (ipap_flg[0] == 'y')
{
if ((tau1+pwN*2.0) < pwS2) delay((pwS2*0.5-tau1*0.5-pwN)*0.5);
if (grad3_flg[0]== 'y')
delay(taunh*0.5-shpw1*0.533-pwS1*0.5+(gt3*2.0+2.0*gstab+pwN*3.0));
else
delay(taunh*0.5-shpw1*0.533-pwS1*0.5+pwN);
hn_simshapedpulse(refpat,shbw,shofs-4.8,"bip720_50_20",40.0,0.0, zero, zero, 0.0, 0.0);
obspower(shlvl2);
obspwrf(4095.0);
compo1_pulse(shname2,shpw2,pwN,shdmf2,t1,tau1,c13refoc,pwS,taunh,pwS1,pwS2,gt1,gt3,gzlvl3,grad3_flg,gstab);
obspower(shlvl2);
obspwrf(4095.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
h_sim3shapedpulse(refpat,shbw,shofs-4.8,0.0,2.0*pwN, one, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh*0.5-gt1-gstab-POWER_DELAY-pwS1*0.5);
if ((tau1+pwN*2.0) < pwS2) delay((pwS2*0.5-tau1*0.5-pwN)*0.5);
if (ab_flg[0] == 'a')
dec2rgpulse(pwN,one,0.0,0.0);
else
dec2rgpulse(pwN,three,0.0,0.0);
if (grad3_flg[0]== 'y')
{
delay(gt3+gstab);
dec2rgpulse(pwN*2.0,zero,0.0,0.0);
zgradpulse(gzlvl3,gt3);
delay(gstab);
}
}
else
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
if ((tau1+pwN*2.0) < pwS2)
delay(taunh-gt1-gstab-shpw1*0.533-adjust-(pwS2*0.5-tau1*0.5-pwN)*0.5);
else
delay(taunh-gt1-gstab-shpw1*0.533-adjust);
obspower(shlvl2);
obspwrf(tpwrsf);
compo_pulse(shname2,shpw2,pwN,shdmf2,t1,tau1,c13refoc,pwS);
obspower(shlvl1);
obspwrf(4095.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
if ((tau1+pwN*2.0) < pwS2)
delay(taunh-gt1-gstab-POWER_DELAY-(pwS2*0.5-tau1*0.5-pwN)*0.5);
else
delay(taunh-gt1-gstab-POWER_DELAY);
}
}
else
{
if ((ni == 0) || (ni == 1))
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh-gt1-gstab-WFG_START_DELAY+pwN*4.0-shpw1*0.533-adjust);
obspwrf(tpwrsf);
obspower(shlvl2);
xmtrphase(zero);
xmtron();
obsunblank();
obsprgon(shname2,1/shdmf2,9.0);
delay(shpw2);
obsprgoff();
obsblank();
xmtroff();
obspower(shlvl2);
obspwrf(4095.0);
dec2rgpulse(pwN,t1,0.0,0.0);
dec2rgpulse(pwN*2.0,zero,0.0,0.0);
dec2rgpulse(pwN,zero,0.0,0.0);
}
else
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh-gt1-gstab-shpw1*0.533-adjust);
obspower(shlvl2);
obspwrf(tpwrsf);
compo_pulse(shname2,shpw2,pwN,shdmf2,t1,tau1,c13refoc,pwS);
obspower(shlvl1);
obspwrf(4095.0);
}
if (ipap_flg[0] == 'y')
{
if (ab_flg[0] == 'b')
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh-gt1-gstab-pwN-POWER_DELAY);
dec2rgpulse(pwN,one,0.0,0.0);
}
else
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh*0.5-gt1-pwN-gstab);
dec2rgpulse(pwN*2.0,zero,0.0,0.0);
delay(taunh*0.5-pwN-POWER_DELAY);
}
}
else
{
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(taunh-gt1-gstab-POWER_DELAY);
}
}
dec2power(dpwr2);
lk_sample();
setreceiver(t2);
status(C);
}
pulsesequence()
{
double gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
satdly = getval("satdly"),
mix = getval("mix"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
pwClvl = getval("pwClvl"),
pwC = getval("pwC"),
scubad = getval("scubad");
int iphase = (int) (getval("phase") + 0.5);
char sspul[MAXSTR],mfsat[MAXSTR],scuba[MAXSTR],
C13refoc[MAXSTR],N15refoc[MAXSTR],CNrefoc[MAXSTR];
getstr("C13refoc",C13refoc);
getstr("N15refoc",N15refoc);
getstr("CNrefoc",CNrefoc);
getstr("mfsat",mfsat);
getstr("sspul", sspul);
getstr("scuba",scuba);
loadtable("tnnoesy");
sub(ct,ssctr,v12);
getelem(t1,v12,v1);
getelem(t2,v12,v2);
getelem(t3,v12,v5);
getelem(t4,v12,v8);
getelem(t5,v12,v9);
getelem(t6,v12,oph);
assign(zero,v3);
assign(one,v4);
if (iphase == 2)
{incr(v1); incr(v2); incr(v3); incr(v4);}
/* HYPERCOMPLEX MODE USES REDFIELD TRICK TO MOVE AXIAL PEAKS TO EDGE */
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v6);
if ((iphase==1)||(iphase==2))
{add(v1,v6,v1); add(v2,v6,v2); add(oph,v6,oph);
add(v3,v6,v3); add(v4,v6,v4);}
/* CHECK CONDITIONS */
if ((dm[A]=='y') || (dm[B] == 'y') || (dm[C] == 'y'))
{
abort_message("Set dm to be nnnn or nnny");
}
if ((dm2[A]=='y') || (dm2[B] == 'y') || (dm2[C] == 'y'))
{
abort_message("Set dm2 to be nnnn or nnny");
}
/* BEGIN THE ACTUAL PULSE SEQUENCE */
decpower(pwClvl); dec2power(pwNlvl);
if (CNrefoc[A] == 'y')
{
decpower(pwClvl-3.0); pwC=1.4*pwC;
dec2power(pwNlvl-3.0); pwN=1.4*pwN;
}
if (sspul[A] == 'y')
{
zgradpulse(gzlvl1,gt1);
delay(1.0e-4);
rgpulse(pw,zero,rof1,rof1);
zgradpulse(gzlvl1,gt1);
delay(1.0e-4);
}
status(A);
if (satmode[A] == 'y')
{
if (d1 > satdly) delay(d1 - satdly);
if (mfsat[A] == 'y')
{obsunblank(); mfpresat_on(); delay(satdly); mfpresat_off(); obsblank();}
else
{
obspower(satpwr);
rgpulse(satdly,v1,rof1,rof1);
}
obspower(tpwr);
if (scuba[0] == 'y')
{
rgpulse(pw,v3,2.0e-6,0.0);
rgpulse(2.0*pw,v4,2.0e-6,0.0);
rgpulse(pw,v3,2.0e-6,0.0);
delay(scubad/2.0);
rgpulse(pw,v3,2.0e-6,0.0);
rgpulse(2.0*pw,v4,2.0e-6,0.0);
rgpulse(pw,v3,2.0e-6,0.0);
delay(scubad/2.0);
}
}
else delay(d1);
obsstepsize(45.0);
initval(7.0,v7);
xmtrphase(v7);
status(B);
rgpulse(pw,v2,rof1,0.0);
xmtrphase(zero);
if (d2>0.0)
{
if ((C13refoc[A] == 'n') && (N15refoc[A] == 'n') && (CNrefoc[A] == 'n'))
{
delay(d2-4.0*pw/PI-SAPS_DELAY-rof1);
}
else if ((C13refoc[A] == 'n') && (N15refoc[A] == 'n') && (CNrefoc[A] == 'y'))
{
if (pwN > 2.0*pwC)
{
if (d2/2.0 > (pwN +0.64*pw+rof1))
{
delay(d2/2.0-pwN-0.64*pw-SAPS_DELAY);
dec2rgpulse(pwN-2.0*pwC,zero,0.0,0.0);
sim3pulse(0.0,pwC,pwC, zero,zero,zero, 0.0, 0.0);
sim3pulse(0.0,2.0*pwC,2.0*pwC, zero,one,zero, 0.0, 0.0);
sim3pulse(0.0,pwC,pwC, zero,zero,zero, 0.0, 0.0);
dec2rgpulse(pwN-2.0*pwC,zero,0.0,0.0);
delay(d2/2.0-pwN-0.64*pw-rof1);
}
else
delay(d2-4.0*pw/PI-SAPS_DELAY-rof1);
}
else
{
if (d2/2.0 > (pwN +pwC+ 0.64*pw+rof1))
{
delay(d2/2.0-pwN-pwC-0.64*pw-SAPS_DELAY);
decrgpulse(pwC,zero,0.0,0.0);
sim3pulse(0.0,2.0*pwC,2.0*pwN, zero,one,zero, 0.0, 0.0);
decrgpulse(pwC,zero,0.0,0.0);
delay(d2/2.0-pwN-pwC-0.64*pw-rof1);
}
else
delay(d2-4.0*pw/PI-SAPS_DELAY-rof1);
}
}
else if ((C13refoc[A] == 'n') && (N15refoc[A] == 'y') && (CNrefoc[A] == 'n'))
{
if (d2/2.0 > (pwN + 0.64*pw+rof1))
{
delay(d2/2.0-pwN-0.64*pw-SAPS_DELAY);
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(d2/2.0-pwN-0.64*pw-rof1);
}
else
delay(d2-1.28*pw-SAPS_DELAY-rof1);
}
else if ((C13refoc[A] == 'y') && (N15refoc[A] == 'n') && (CNrefoc[A] == 'n'))
{
if (d2/2.0 > (2.0*pwC + 0.64*pw +rof1))
{
delay(d2/2.0-2.0*pwC-0.64*pw-SAPS_DELAY);
decrgpulse(pwC,zero,0.0,0.0);
decrgpulse(2.0*pwC, one, 0.0, 0.0);
decrgpulse(pwC,zero,0.0,0.0);
delay(d2/2.0-2.0*pwC-0.64*pw-rof1);
}
else
delay(d2-1.28*pw-SAPS_DELAY -rof1);
}
else
{
abort_message("C13refoc, N15refoc, and CNrefoc must be nnn, nny, nyn, or ynn");
}
}
rgpulse(pw,v5,rof1,1.0e-6);
status(C);
decpower(dpwr); dec2power(dpwr2);
if (satmode[C] == 'y')
{
if (mfsat[C] == 'y')
{obsunblank(); mfpresat_on(); delay(mix*0.7); mfpresat_off(); obsblank();
zgradpulse(gzlvl2,gt2);
obsunblank(); mfpresat_on(); delay(mix*0.3-gt2); mfpresat_off(); obsblank();}
else
{
obspower(satpwr);
rgpulse(mix*0.7,v8,rof1,rof1);
zgradpulse(gzlvl2,gt2);
rgpulse(mix*0.3-gt2,v8,rof1,rof1);
}
obspower(tpwr);
}
else
{
delay(mix*0.7);
zgradpulse(gzlvl2,gt2);
delay(mix*0.3-gt2);
}
status(D);
rgpulse(pw,v9,rof1,rof2);
/* Phase cycle: .satdly(t1)..pw(t2)..d2..pw(t3)..mix(t4)..pw(t5)..at(t6)
(for phase=1; for phase=2 incr t1 + t2; for TPPI add two to t1,t2,t6)
t1 = 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3
t2 = 2 0 2 0 2 0 2 0 3 1 3 1 3 1 3 1 2 0 2 0 2 0 2 0 3 1 3 1 3 1 3 1
t3 = 0 0 0 0 0 0 0 0 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 3 3 3 3 3 3 3 3
t4 = 3 3 1 1 0 0 2 2 0 0 2 2 1 1 3 3 3 3 1 1 0 0 2 2 0 0 2 2 1 1 3 3
t5 = 0 0 2 2 1 1 3 3 1 1 3 3 2 2 0 0 0 0 2 2 1 1 3 3 1 1 3 3 2 2 0 0
t6 = 0 2 2 0 1 3 3 1 1 3 3 1 2 0 0 2 2 0 0 2 3 1 1 3 3 1 1 3 0 2 2 0 */
}
void pulsesequence()
{
char c1d[MAXSTR]; /* option to record only 1D C13 spectrum */
int ncyc;
double tau1 = 0.002, /* t1 delay */
post_del = 0.0,
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
compC = getval("compC"),
compH = getval("compH"),
mixpwr = getval("mixpwr"),
jCH = getval("jCH"),
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gzlvl0 = getval("gzlvl0"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
grec = getval("grec"),
phase = getval("phase");
getstr("c1d",c1d);
ncyc=1;
if(jCH > 0.0)
tau1 = 0.25/jCH;
dbl(ct, v1); /* v1 = 0 */
mod4(v1,oph);
hlv(ct,v2);
add(v2,v1,v2);
if (phase > 1.5)
incr(v1); /* hypercomplex phase increment */
initval(2.0*(double)((int)(d2*getval("sw1")+0.5)%2),v10);
add(v1,v10,v1);
add(oph,v10,oph);
mod4(v1,v1); mod4(v2,v2); mod4(oph,oph);
assign(zero,v3);
if(FIRST_FID)
{
HHmix = pbox_mix("HHmix", "DIPSI2", mixpwr, pw*compH, tpwr);
if(c1d[A] == 'n')
{
opx("CHdec"); setwave("WURST2 30k/1.2m"); pbox_par("steps","600"); cpx(pwC*compC, pwClvl);
CHdec = getDsh("CHdec");
}
}
ncyc = (int) (at/HHmix.pw) + 1;
post_del = ncyc*HHmix.pw - at;
/* BEGIN PULSE SEQUENCE */
status(A);
zgradpulse(gzlvl0, gt0);
rgpulse(pw, zero, 0.0, 0.0); /* destroy H-1 magnetization*/
zgradpulse(gzlvl0, gt0);
delay(1.0e-4);
obspower(tpwr);
txphase(v1);
decphase(zero);
dec2phase(zero);
presat();
obspower(tpwr);
delay(1.0e-5);
status(B);
if(c1d[A] == 'y')
{
rgpulse(pw,v1,0.0,0.0); /* 1H pulse excitation */
delay(d2);
rgpulse(pw,two,0.0,0.0); /* 1H pulse excitation */
assign(oph,v3);
}
else
{
decunblank(); pbox_decon(&CHdec);
rgpulse(pw,v1,0.0,0.0); /* 1H pulse excitation */
txphase(zero);
delay(d2);
pbox_decoff(); decblank();
decpower(pwClvl); decpwrf(4095.0);
delay(tau1 - POWER_DELAY);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
txphase(one); decphase(one); dec2phase(one);
delay(tau1);
simpulse(pw, pwC, one, one, 0.0, 0.0);
txphase(zero); decphase(zero); dec2phase(zero);
delay(tau1);
simpulse(2.0*pw, 2.0*pwC, zero, zero, 0.0, 0.0);
delay(tau1);
simpulse(0.0, pwC, zero, zero, 0.0, 0.0);
}
zgradpulse(gzlvl1, gt1);
delay(grec);
simpulse(0.0, pwC, zero, v3, 0.0, rof2);
txphase(v2);
obsunblank(); pbox_xmtron(&HHmix);
status(C);
setactivercvrs("ny");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
delay(post_del);
pbox_xmtroff(); obsblank();
zgradpulse(gzlvl2, gt2);
obspower(tpwr);
delay(grec);
rgpulse(pw,zero,0.0,rof2); /* 1H pulse excitation */
status(D);
setactivercvrs("yn");
startacq(alfa);
acquire(np,1.0/sw);
endacq();
}
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
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 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()
{
/* DECLARE AND LOAD VARIABLES */
char f1180[MAXSTR], /* Flag to start t1 @ halfdwell */
CT[MAXSTR],
refocN15[MAXSTR],
refocCO[MAXSTR], COshape[MAXSTR],
C180shape[MAXSTR];
int t1_counter,icosel; /* used for states tppi in t1 */
double tau1, /* t1 delay */
tauch = getval("tauch"), /* 1/4J evolution delay */
tauch1 = getval("tauch1"), /* 1/4J or 1/8JC13H evolution delay */
timeCC = getval("timeCC"), /* 13C constant-time if needed*/
corrD, corrB, /* small correction delays */
dof_dec = getval("dof_dec"), /*decoupler offset for decoupling during acq - folding */
pwClvl = getval("pwClvl"), /* power for hard C13 pulses */
pwC180lvl = getval("pwC180lvl"), /*power levels for 180 shaped pulse */
pwC180lvlF = getval("pwC180lvlF"),
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
pwC180 = getval("pwC180"), /* shaped 180 pulse on C13channl */
sw1 = getval("sw1"),
pwNlvl = getval("pwNlvl"),
pwN = getval("pwN"),
pwCOlvl = getval("pwCOlvl"),
pwCO = getval("pwCO"),
gstab = getval("gstab"),
gstab1 = getval("gstab1"), /* recovery for club sandwitch, short*/
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gt3 = getval("gt3"), /* other gradients */
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gt9 = getval("gt9"),
gzlvl0 = getval("gzlvl0"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5"),
gzlvl9 = getval("gzlvl9");
getstr("f1180",f1180);
getstr("C180shape",C180shape);
getstr("COshape",COshape);
getstr("refocCO",refocCO);
getstr("refocN15",refocN15);
getstr("CT",CT);
/* LOAD PHASE TABLE */
settable(t2,2,phi2);
settable(t3,1,phi3);
settable(t4,1,phi4);
settable(t10,1,phix);
settable(t11,4,rec);
/* INITIALIZE VARIABLES */
/* CHECK VALIDITY OF PARAMETER RANGES */
/* like in olde good times */
if((dm[A] == 'y' || dm[B] == 'y'))
{
text_error("incorrect dec1 decoupler flags! Should be 'nnn' or 'nny' ");
psg_abort(1);
}
if((dm2[A] != 'n' || dm2[B] != 'n' || dm2[C] != 'n'))
{
text_error("incorrect dec2 decoupler flags! Should be 'nnn' ");
psg_abort(1);
}
if( dpwr2 > 45 )
{
text_error("don't fry the probe, DPWR2 too large! ");
psg_abort(1);
}
if( dpwr > 50 )
{
text_error("don't fry the probe, DPWR too large! ");
psg_abort(1);
}
if( (pw > 30.0e-6) )
{
text_error("don't fry the probe, pw too high ! ");
psg_abort(1);
}
if( (pwC > 200.0e-6) )
{
text_error("don't fry the probe, pwC too high ! ");
psg_abort(1);
}
/* PHASES AND INCREMENTED TIMES */
/* Phase incrementation for hypercomplex 2D data, States-TPPI-Haberkorn element */
if (phase1 == 1) {
tsadd(t10,2,4);
icosel = +1;
}
else {
icosel = -1;
}
/* Calculate modifications to phases for States-TPPI acquisition */
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{
tsadd(t2,2,4);
tsadd(t11,2,4);
}
/* Set up f1180 */
tau1 = d2;
if((f1180[A] == 'y') && (ni > 0.0))
{
tau1 += ( 1.0 / (2.0*sw1) );
if(tau1 < 0.2e-6) tau1 = 0.0;
}
tau1 = tau1/2.0;
if(CT[A]=='y') {
refocN15[A]='y';
refocCO[A]='y';
if ( (timeCC-ni/sw1*0.5)*0.5 -pw -pwCO*0.5 -pwN-rof1 <0.2e-6 )
{
text_error("too many increments in t1 for a constant time");
psg_abort(1);
}
}
/*temporary*/
/* correction delays */
corrB=0.0;
corrD=2.0/M_PI*pwC-pw-rof1;
if (corrD < 0.0) {
corrB=-corrD;
corrD=0.0;
}
/* BEGIN PULSE SEQUENCE */
status(A);
obsoffset(tof);
obspower(tpwr);
decpower(pwClvl);
decpwrf(4095.0);
decoffset(dof);
obspwrf(4095.0);
delay(d1);
txphase(zero);
decphase(zero);
decrgpulse(pwC, zero, rof1, rof1);
zgradpulse(gzlvl0, gt0);
delay(2.0*gstab);
delay(5.0e-3);
/************ H->C13 */
txphase(zero);
decphase(zero);
rgpulse(pw,zero,rof1,rof1);
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay(gstab);
zgradpulse(gzlvl5, gt5);
delay(tauch - gt5 -gstab);
simshaped_pulse("hard",C180shape,2.0*pw,pwC180,zero,zero, 0.0, 0.0);
decpower(pwClvl);
decpwrf(4095.0);
txphase(one);
delay(tauch - gt5-gstab);
zgradpulse(gzlvl5, gt5);
delay(gstab);
rgpulse(pw, one, rof1,rof1);
/*** purge */
dec2power(pwNlvl);
dec2pwrf(4095.0);
txphase(zero);
decphase(t2);
delay(gstab);
zgradpulse(gzlvl2, gt2);
delay(3.0*gstab);
/* evolution on t1 */
/* real time here */
if(CT[A]=='n')
{
decrgpulse(pwC, t2, 0.0, 0.0);
decphase(zero);
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay(tau1);
obsunblank();
rgpulse(pw,one,rof1,0.0);
rgpulse(2.0*pw,zero,0.0,0.0);
rgpulse(pw,one,0.0,rof1);
obsblank();
if(refocN15[A]=='y') dec2rgpulse(2.0*pwN,zero,0.0,0.0); /*n15 refocusing */
if(refocCO[A]=='y')
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
delay(tau1);
/*************** CODING with CLUB sandwitches and BIPS */
zgradpulse(gzlvl3*icosel, gt3);
delay(gstab1);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
zgradpulse(-gzlvl3*icosel, gt3);
delay(gstab1 +rof1*2.0+pw*4.0 +pwC*4.0/M_PI +2.0*POWER_DELAY+2.0*PWRF_DELAY);
if(refocN15[A]=='y') delay(2.0*pwN); /*n15 refocusing */
if(refocCO[A]=='y') /* ghost CO pulse */
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
zgradpulse(-gzlvl3*icosel, gt3);
delay(gstab1);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
zgradpulse( gzlvl3*icosel, gt3);
decpower(pwClvl);
decpwrf(4095.0);
delay(gstab1);
} /* end of if bracket for real-time */
/*^^^^^^^ end of real time */
/* CONSTANT TIME VESION: */
if(CT[A]=='y')
{
decrgpulse(pwC, t2, 0.0, 0.0);
/* timeCC-t1 evolution */
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay((timeCC-tau1)*0.5 -2.0*pw -pwCO*0.5 -pwN-rof1);
obsunblank();
rgpulse(pw,one,rof1,0.0);
rgpulse(2.0*pw,zero,0.0,0.0);
rgpulse(pw,one,0.0,rof1);
obsblank();
if(refocN15[A]=='y') dec2rgpulse(2.0*pwN,zero,0.0,0.0); /*n15 refocusing */
if(refocCO[A]=='y' )
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
delay((timeCC-tau1)*0.5 -2.0*pw -pwCO*0.5 -pwN-rof1);
/* end of timeCC-t1 evolution */
/* 180 on carbons in T1 */
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
/* timeCC+t1 evolution + encoding */
delay((timeCC+tau1)*0.5 -2.0*pw -pwCO*0.5 -pwN -rof1);
obsunblank();
rgpulse(pw,one,rof1,0.0);
rgpulse(2.0*pw,zero,0.0,0.0);
rgpulse(pw,one,0.0,rof1);
obsblank();
if(refocN15[A]=='y') dec2rgpulse(2.0*pwN,zero,0.0,0.0); /*n15 refocusing */
if(refocCO[A]=='y' )
{
decpower(pwCOlvl);
decshaped_pulse(COshape,pwCO,zero, 3.0e-6, 3.0e-6);
decpower(pwC180lvl);
}
delay((timeCC+tau1)*0.5 -2.0*pw -pwCO*0.5 -pwN -rof1);
zgradpulse(-gzlvl3*icosel, gt3*2.0); /* coding */
delay(gstab +pwC*4.0/M_PI);
decshaped_pulse(C180shape,pwC180,zero, 2e-6, 2e-6); /* ghost BIP pulse */
zgradpulse(gzlvl3*icosel, gt3*2.0); /* coding */
delay(gstab);
} /*^^^^^^^ end of CONSTANT time bracket*/
/*reverse INPET */
simpulse(pw, pwC, zero, t10, 0.0, 0.0);
delay(gstab);
zgradpulse(gzlvl4,gt4);
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay(tauch - gt4 - gstab -corrD-pwC180 -POWER_DELAY-PWRF_DELAY);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
delay(corrD);
rgpulse(2.0*pw,zero,rof1,rof1);
zgradpulse(gzlvl4,gt4);
delay(tauch - gt4 - gstab -corrB-pwC180 -POWER_DELAY-PWRF_DELAY);
delay(gstab);
decphase(one);
txphase(one);
decshaped_pulse(C180shape,pwC180,zero, 0.0, 0.0);
decpower(pwClvl);
decpwrf(4095.0);
delay(corrB);
simpulse(pw, pwC, one, one, 0.0, 0.0);
decpower(pwC180lvl);
decpwrf(pwC180lvlF);
delay(gstab-POWER_DELAY-PWRF_DELAY -WFG_START_DELAY);
zgradpulse(gzlvl5,gt5);
delay(tauch1 - gt5 - gstab);
simshaped_pulse("hard",C180shape,2.0*pw,pwC180,zero,zero, 0.0, 0.0);
delay(tauch1- gt5 - gstab - WFG_STOP_DELAY);
zgradpulse(gzlvl5,gt5);
delay(gstab-rof1);
rgpulse(pw, zero,rof1,rof1);
/* echo and decoding */
delay(gstab+gt9-rof1+10.0e-6 + 2.0*POWER_DELAY+PWRF_DELAY +2.0*GRADIENT_DELAY);
decoffset(dof_dec);
rgpulse(2.0*pw, zero,rof1,rof1);
delay(10.0e-6);
zgradpulse(gzlvl9,gt9);
decpower(dpwr);
decpwrf(4095.0);
decoffset(dof_dec);
dec2power(dpwr2); /* POWER_DELAY EACH */
delay(gstab);
status(C);
setreceiver(t11);
}
void pulsesequence()
{
double p1lvl,
trim,
mix,
window,
cycles;
char sspul[MAXSTR];
/* LOAD AND INITIALIZE VARIABLES */
mix = getval("mix");
trim = getval("trim");
p1lvl = getval("p1lvl");
window=getval("window");
getstr("sspul", sspul);
/* CHECK CONDITIONS */
if ((phase1 != 3) && (arrayelements > 2))
{
fprintf(stdout, "PHASE=3 is required if MIX is arrayed!\n");
psg_abort(1);
}
if (satdly > 9.999)
{
fprintf(stdout, "Presaturation period is too long.\n");
psg_abort(1);
}
if (!newtransamp)
{
fprintf(stdout, "TOCSY requires linear amplifiers on transmitter.\n");
fprintf(stdout, "Use DECTOCSY with the appropriate re-cabling,\n");
psg_abort(1);
}
if ((p1 == 0.0) && (ix == 1))
fprintf(stdout, "Warning: P1 has a zero value.\n");
if ((rof1 < 9.9e-6) && (ix == 1))
fprintf(stdout,"Warning: ROF1 is less than 10 us\n");
if (satpwr > 40)
{
printf("satpwr too large - acquisition aborted./n");
psg_abort(1);
}
/* STEADY-STATE PHASECYCLING */
/* This section determines if the phase calculations trigger off of (SS - SSCTR)
or off of CT */
ifzero(ssctr);
hlv(ct, v13);
mod2(ct, v1);
hlv(ct, v2);
elsenz(ssctr);
sub(ssval, ssctr, v12); /* v12 = 0,...,ss-1 */
hlv(v12, v13);
mod2(v12, v1);
hlv(v12, v2);
endif(ssctr);
/* CALCULATE PHASES */
/* A 2-step cycle is performed on the first pulse (90 degrees) to suppress
axial peaks in the first instance. Second, the 2-step F2 quadrature image
suppression subcycle is added to all pulse phases and receiver phase.
Finally, a 2-step cycle is performed on the spin-lock pulses. */
mod2(v13, v13);
dbl(v1, v1);
incr(v1);
hlv(v2, v2);
mod2(v2, v2);
dbl(v2, v2);
incr(v2);
add(v13, v2, v2);
sub(v2, one, v3);
add(two, v2, v4);
add(two, v3, v5);
add(v1, v13, v1);
assign(v1, oph);
if (phase1 == 2)
incr(v1);
if (phase1 == 3)
add(v1, id2, v1);
/*HYPERCOMPLEX MODE USES REDFIELD TRICK TO MOVE AXIAL PEAKS TO EDGE*/
if ((phase1==1)||(phase1==2))
{
initval(2.0*(double)(d2_index%2),v6);
add(v1,v6,v1); add(oph,v6,oph);
}
/* CALCULATE AND INITIALIZE LOOP COUNTER */
if (pw > 0.0)
{
cycles = (mix - trim) / (64.66*pw+32*window);
cycles = 2.0*(double) (int) (cycles/2.0);
}
else
{
cycles = 0.0;
}
initval(cycles, v9); /* V9 is the MIX loop count */
/* BEGIN ACTUAL PULSE SEQUENCE CODE */
status(A);
obspower(p1lvl);
if (sspul[0] == 'y')
{
rgpulse(1000*1e-6, zero, rof1, 0.0e-6);
rgpulse(1000*1e-6, one, 0.0e-6, rof1);
}
hsdelay(d1);
if (satmode[A] == 'y')
{ obspower(satpwr);
rgpulse(satdly,zero,rof1,rof2);
obspower(p1lvl);}
status(B);
rgpulse(p1, v1, rof1, 1.0e-6);
if (satmode[B] =='y')
{
if (d2 > 0.0)
{
obspower(satpwr);
rgpulse(d2 - (2*POWER_DELAY) - 1.0e-6 - (2*p1/3.1416),zero,0.0,0.0);
}
}
else
{
if (d2 > 0.0)
delay(d2 - POWER_DELAY - 1.0e-6 - (2*p1/3.1416));
}
rcvroff();
obsunblank();
obspower(tpwr);
txphase(v3);
xmtron();
delay(trim);
if (cycles > 1.0)
{
starthardloop(v9);
mleva(window); mleva(window); mlevb(window); mlevb(window);
mlevb(window); mleva(window); mleva(window); mlevb(window);
mlevb(window); mlevb(window); mleva(window); mleva(window);
mleva(window); mlevb(window); mlevb(window); mleva(window);
rgpulse(.66*pw,v3,rof1,rof2);
endhardloop();
}
txphase(v13);
xmtroff();
/* detection */
delay(rof2);
rcvron();
obsblank();
status(C);
}
pulsesequence()
{
char Cdseq[MAXSTR], refoc[MAXSTR], sspul[MAXSTR];
int mult = (0.5 + getval("mult"));
double rg1 = 2.0e-6,
j1xh = getval("j1xh"),
gt0 = getval("gt0"),
gzlvl0 = getval("gzlvl0"),
pp = getval("pp"),
pplvl = getval("pplvl"),
compH = getval("compH"),
Cdpwr = getval("Cdpwr"),
Cdres = getval("Cdres"),
tau1 = 0.002,
tau2 = 0.002,
Cdmf = getval("Cdmf");
shape hdx;
/* Get new variables from parameter table */
getstr("Cdseq", Cdseq);
getstr("refoc", refoc);
getstr("sspul", sspul);
if (j1xh < 1.0) j1xh = 150.0;
hdx = pboxHT_F1i("gaus180", pp*compH, pplvl); /* HADAMARD stuff */
if (getval("htcal1") > 0.5) /* Optional fine power calibration */
hdx.pwr = getval("htpwr1");
if(j1xh > 0.0)
tau1 = 0.25/j1xh;
tau2 = tau1;
if (mult > 2) tau2=tau1/2.5;
else if ((mult == 0) || (mult == 2)) tau2=tau1/2.0;
dbl(ct, v1); /* v1 = 02 */
add(two,v1,oph);
mod4(oph,oph); /* oph = 02 */
/* Calculate delays */
if (dm[0] == 'y')
{
abort_message("decoupler must be set as dm=nny\n");
}
if(refoc[A]=='n' && dm[C]=='y')
{
abort_message("with refoc=n decoupler must be set as dm=nnn\n");
}
/* Relaxation delay */
status(A);
delay(0.05);
zgradpulse(gzlvl0,gt0);
if (sspul[0] == 'y')
{
obspower(tpwr); decpower(pplvl);
simpulse(pw, pp, zero, zero, rg1, rg1); /* destroy H and C magnetization */
zgradpulse(gzlvl0,gt0);
}
delay(d1);
status(B);
obspower(Cdpwr);
obsunblank(); xmtron();
obsprgon(Cdseq, 1.0/Cdmf, Cdres);
pbox_decpulse(&hdx, zero, rg1, rg1);
obsprgoff(); xmtroff(); obsblank();
obspower(tpwr); decpower(pplvl);
delay(2.0e-4);
decrgpulse(pp, v1, rg1, rg1);
delay(tau1 - POWER_DELAY);
simpulse(2.0*pw, 2.0*pp, zero, zero, rg1, rg1);
txphase(one); decphase(one);
delay(tau1);
simpulse(pw, pp, one, one, rg1, rg1);
if(refoc[A]=='y')
{
txphase(zero); decphase(zero);
delay(tau2);
simpulse(2.0*pw, 2.0*pp, zero, zero, rg1, rg1);
delay(tau2 - rof2 - POWER_DELAY);
decrgpulse(pp, zero, rg1, rof2);
}
decpower(dpwr);
status(C);
}
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() {
// 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()
{
/* DECLARE AND LOAD VARIABLES */
int t1_counter; /* used for states tppi in t1 */
double tau1, /* t1 delay */
TC = getval("TC"), /* Constant delay 1/(JCC) ~ 13.5 ms */
mix = getval("mix"), /* TOCSY mixing time */
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
rf0, /* maximum fine power when using pwC pulses */
/* 90 degree pulse at Cab (46ppm), first off-resonance null at CO (174ppm) */
pwC1, /* 90 degree pulse length on C13 at rf1 */
rf1, /* fine power for 5.1 kHz rf for 600MHz magnet */
/* 180 degree pulse at Ca (46ppm), first off-resonance null at CO(174ppm) */
pwC2, /* 180 degree pulse length at rf2 */
rf2, /* fine power for 11.4 kHz rf for 600MHz magnet */
compC = getval("compC"), /* adjustment for C13 amplifier compression */
sw1 = getval("sw1"),
gt1 = getval("gt1"),
gzlvl1 = getval("gzlvl1"),
gt2 = getval("gt2"),
gzlvl2 = getval("gzlvl2"),
gstab = getval("gstab"),
ppm,
co_ofs = getval("co_ofs"), /* offset for C' */
co_bw = getval("co_bw"), /* bandwidth for C' */
copwr = getval("copwr"), /* power for C' decoupling. Get from CO_dec.DEC*/
codmf = getval("codmf"), /* dmf for C' decoupling. Get from CO_dec.DEC */
codres = getval("codres"), /* dres for C' decoupling. Get from CO_dec.DEC */
mixbw, /* band width for mixing shape */
mixpwr = getval("mixpwr"), /* power for CC mixing. Get from ccmix.DEC*/
mixdmf= getval("mixdmf"), /* dmf for CC decoupling. Get from ccmix.DEC */
mixdres = getval("mixdres"); /* dres for CC decoupling. Get from ccmix.DEC */
/* LOAD PHASE TABLE */
settable(t1,2,phi1);
settable(t2,1,phi2);
settable(t3,1,phi3);
settable(t12,2,rec);
setautocal(); /* activate auto-calibration */
/* INITIALIZE VARIABLES */
/* maximum fine power for pwC pulses */
rf0 = 4095.0;
/* 90 degree pulse on Cab, null at CO 128ppm away */
pwC1 = sqrt(15.0)/(4.0*128.0*sfrq);
rf1 = (compC*4095.0*pwC)/pwC1;
rf1 = (int) (rf1 + 0.5);
/* 180 degree pulse on Cab, null at CO 128ppm away */
pwC2 = sqrt(3.0)/(2.0*128.0*sfrq);
rf2 = (4095.0*compC*pwC*2.0)/pwC2;
rf2 = (int) (rf2 + 0.5);
if( rf2 > 4295 )
{ printf("increase pwClvl"); psg_abort(1);}
if(( rf2 > 4095 ) && (rf2 <4296)) rf2=4095;
/* CHECK VALIDITY OF PARAMETER RANGES */
if ( 0.5*ni*1/(sw1) > TC)
{ printf(" ni is too big. Make ni equal to %d or less.\n",
((int)((TC)*2.0*sw1))); psg_abort(1);}
/* PHASES AND INCREMENTED TIMES */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase1 == 2) tsadd(t1,1,4);
tau1 = d2;
tau1 = tau1/2.0;
/* Calculate modifications to phases for States-TPPI acquisition */
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t12,2,4); }
if (autocal[0] == 'y')
{
if(FIRST_FID)
{
ppm = getval("sfrq");
ofs_check(C13ofs);
co_ofs = (C13ofs+118.0)*ppm; co_bw = 20*ppm;
CO_dec = pbox_Dsh("CO_dec", "SEDUCE1", co_bw, co_ofs, pwC*compC, pwClvl);
copwr = CO_dec.pwr; codmf = CO_dec.dmf; codres = CO_dec.dres;
mixbw = sw1;
mix_seq = pbox_Dsh("mix_seq", "FLOPSY8", mixbw, 0.0, pwC*compC, pwClvl);
mixpwr = mix_seq.pwr; mixdmf = mix_seq.dmf; mixdres = mix_seq.dres;
}
}
/* BEGIN PULSE SEQUENCE */
status(A);
delay(d1);
if ( dm3[B] == 'y' )
{ lk_hold(); lk_sampling_off();} /*freezes z0 correction, stops lock pulsing*/
rcvroff();
obspower(pwClvl);
decpower(tpwr);
dec2power(dpwr2);
dec3power(dpwr3);
obspwrf(rf1); /*fine power for Cab 90 degree pulse */
obsoffset(tof); /*13C carrier at 46 ppm */
txphase(zero);
delay(1.0e-5);
status(B);
rgpulse(pwC1, t1, 0.0,0.0);
/* xxxxxxxxxxxxxxxxxxxxxx 13Cab Constant Time Evolution xxxxxxxxxxxxxxxxxx */
obspower(copwr); obspwrf(rf0); txphase(zero);
obsunblank();
xmtron();
obsprgon("CO_dec",1.0/codmf,codres);
if ( dm3[B] == 'y' ) /* turns on 2H decoupling */
{
dec3rgpulse(1/dmf3,one,10.0e-6,2.0e-6);
dec3unblank();
dec3phase(zero);
delay(2.0e-6);
setstatus(DEC3ch, TRUE, 'w', FALSE, dmf3);
delay(TC - tau1 - pwC2/2 - 1/dmf3);
}
else
delay(TC - tau1 - pwC2/2);
obsprgoff();
xmtroff();
obsblank();
obspower(pwClvl); obspwrf(rf2);
rgpulse(pwC2, zero, 0.0, 0.0);
obspower(copwr); obspwrf(rf0); txphase(zero);
obsunblank();
xmtron();
obsprgon("CO_dec",1.0/codmf,codres);
if ( dm3[B] == 'y' ) /* turns off 2H decoupling */
{
delay(TC + tau1 - pwC2/2 - 1/dmf3);
setstatus(DEC3ch, FALSE, 'c', FALSE, dmf3);
dec3rgpulse(1/dmf3,three,2.0e-6,2.0e-6);
dec3blank();
lk_autotrig(); /* resumes lock pulsing */
}
else
delay(TC + tau1 - pwC2/2);
obsprgoff();
xmtroff();
obsblank();
obspower(pwClvl); obspwrf(rf1);
rgpulse(pwC1,t2,0.0,0.0);
status(C);
zgradpulse(gzlvl1, gt1);
delay(gstab);
/* xxxxxxxxxxxxxxxxxxxxxxxxxxxx FLOPSY 8 Spin lock for mixing xxxxxxxxxxxxxxxxxxxx */
obspower(mixpwr); obspwrf(rf0); txphase(zero);
obsunblank();
xmtron();
obsprgon("mix_seq",1.0/mixdmf,mixdres);
delay(mix-gt1-gt2);
obsprgoff();
xmtroff();
obsblank();
obspower(pwClvl); obspwrf(rf1);
zgradpulse(gzlvl2,gt2);
delay(gstab);
rgpulse(pwC1,t3,0.0,rof2);
getelem(t3,ct,v3);
add(v3,one,v3);
obspower(pwClvl); obspwrf(rf0);
delay(350e-6-rof2);
rgpulse(pwC*2.0,v3,0.0,0.0);
delay(350e-6);
status(D);
setreceiver(t12);
}
pulsesequence()
{
double p1lvl = getval("p1lvl"),
phase = getval("phase"),
mix = getval("mix"),
cycles, d2corr,
gt1 = getval("gt1"),
gzlvl1 = getval("gzlvl1");
int iphase;
char mfsat[MAXSTR],sspul[MAXSTR],T_flg[MAXSTR];
/* LOAD AND INITIALIZE PARAMETERS */
iphase = (int) (phase + 0.5);
satdly = getval("satdly");
satpwr = getval("satpwr");
satfrq = getval("satfrq");
getstr("sspul", sspul);
getstr("mfsat", mfsat);
getstr("satmode", satmode);
getstr("T_flg", T_flg);
/* CALCULATE PHASES AND INITIALIZE LOOP COUNTER FOR MIXING TIME */
settable(t1,8,phi1);
settable(t2,8,phi2);
settable(t3,8,phi3);
settable(t4,8,phi4);
settable(t5,8,phi5);
getelem(t1,ct,v1);
assign(v1,oph);
if (iphase == 2)
incr(v1); /* BC2D hypercomplex method */
/* FOR HYPERCOMPLEX, USE REDFIED TRICK TO MOVE AXIALS TO EDGE */
initval(2.0*(double)(((int)(d2*getval("sw1")+0.5)%2)),v9); /* moves axials */
if ((iphase==2)||(iphase==1)) {add(v1,v9,v1); add(oph,v9,oph);}
cycles = mix / (16.0 * (4e-6 + 2.0*pw));
initval(cycles, v10); /* mixing time cycles */
/* BEGIN ACTUAL PULSE SEQUENCE */
status(C);
obspower(p1lvl);
if (sspul[A] == 'y')
{
zgradpulse(gzlvl1,gt1);
delay(5.0e-5);
rgpulse(p1,zero,rof1,rof1);
zgradpulse(gzlvl1,gt1);
delay(5.0e-5);
}
status(A);
if (satmode[A] == 'y')
{
rcvroff();
if (d1 > satdly) delay(d1-satdly);
if (mfsat[A] == 'y')
{obsunblank(); mfpresat_on(); delay(satdly); mfpresat_off(); obsblank();}
else
{
if (tof != satfrq) obsoffset(satfrq);
obspower(satpwr);
rgpulse(satdly,zero,rof1,rof1);
if (tof != satfrq) obsoffset(tof);
}
obspower(p1lvl);
}
else
{
delay(d1);
rcvroff();
}
status(B);
obsstepsize(45.0);
initval(7.0,v4);
xmtrphase(v4);
rgpulse(p1,v1,rof1,1.0e-6);
xmtrphase(zero);
d2corr = rof1 + 1.0e-6 + (2*p1/3.1416) + SAPS_DELAY;
if (d2 > d2corr) delay(d2 - d2corr); else delay(0.0);
if ((T_flg[0] == 'y')&&(cycles > 1.5))
{
rgpulse(p1,t4,rof1,rof1);
obspower(tpwr);
{
starthardloop(v10);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
rgpulse(2.0*pw,t2,4e-6,0.0);
rgpulse(2.0*pw,t3,4e-6,0.0);
endhardloop();
}
obspower(p1lvl);
rgpulse(p1,t5,rof1,rof2);
}
/* The ROESY spin-lock unit is executed sixteen times within the
hardware loop so that it is of sufficient duration to allow
the acquisition hardware loop to be loaded in behind it on
the last pass through the spin-lock loop. */
else
{
obspower(tpwr);
rgpulse(mix,t2,rof1,rof2); /* cw spin lock */
}
rcvron();
status(C);
}
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();
}
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();
}
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();
}
pulsesequence()
{
char CT_flg[MAXSTR], /* Constant time flag */
shname1[MAXSTR],
shname2[MAXSTR],
shname3[MAXSTR],
f1180[MAXSTR],
Cdecflg[MAXSTR],
Cdecseq[MAXSTR],
grad_flg[MAXSTR]; /*gradient flag */
int t1_counter,
phase;
double d2_init=0.0,
adjust = getval("adjust"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
shlvl1 = getval("shlvl1"),
shlvl2 = getval("shlvl2"),
shlvl3 = getval("shlvl3"),
shdmf2 = getval("shdmf2"),
shpw1 = getval("shpw1"),
shpw2 = getval("shpw2"),
shpw3 = getval("shpw3"),
pwClvl = getval("pwClvl"),
pwC = getval("pwC"),
dpwr = getval("dpwr"),
CT_delay = getval("CT_delay"),
d2 = getval("d2"),
tau1 = getval("tau1"),
tauch = getval("tauch");
getstr("shname1", shname1);
getstr("shname2", shname2);
getstr("shname3", shname3);
getstr("CT_flg", CT_flg);
getstr("grad_flg",grad_flg);
getstr("f1180",f1180);
getstr("Cdecflg",Cdecflg);
getstr("Cdecseq",Cdecseq);
phase = (int) (getval("phase") + 0.5);
settable(t1,2,phi1);
settable(t2,2,phi2);
if (phase == 1) ;
if (phase == 2) tsadd(t1,1,4);
tau1 = d2;
if((f1180[A] == 'y') && (ni > 1.0))
{ tau1 += ( 1.0 / (2.0*sw1) ); if(tau1 < 0.2e-6) tau1 = 0.0; }
tau1 = tau1;
if (f1180[0] == 'y') tau1 = tau1-pwC*4.0/3.0;
if(CT_flg[0] == 'y')
{
if ( (ni/sw1) > (CT_delay-pwC*8.0/3.0))
{ text_error( " ni is too big. Make ni equal to %d or less.\n",
((int)((CT_delay-pwC*8.0/3.0)*sw1)) ); psg_abort(1); }
}
if( ix == 1) d2_init = d2;
t1_counter = (int) ( (d2-d2_init)*sw1 + 0.5 );
if(t1_counter % 2)
{ tsadd(t1,2,4); tsadd(t2,2,4); }
status(A);
decpower(pwClvl);
decoffset(dof);
obsoffset(tof);
zgradpulse(gzlvl2, gt2);
lk_sample();
delay(1.0e-4);
delay(d1-gt2);
if (ix < (int)(2.0*ni)) lk_hold(); /*force lock sampling at end of experiment*/
obspower(shlvl1);
shaped_pulse(shname1,shpw1,zero,2.0e-4,2.0e-6);
zgradpulse(gzlvl1, gt1);
delay(1.0e-4);
if (ni == 0)
{
delay(tauch-gt1-1.0e-4-WFG_START_DELAY+pwC*4.0-adjust);
obspower(shlvl2);
shaped_pulse(shname2,shpw2,zero,2.0e-6,2.0e-6);
obspower(shlvl1);
decrgpulse(pwC,t1,0.0,0.0);
decrgpulse(pwC*2.0,zero,0.0,0.0);
decrgpulse(pwC,zero,0.0,0.0);
}
else
{
delay(tauch-gt1-1.0e-4-adjust);
obspower(shlvl2);
status(B);
if(CT_flg[0] == 'y')
{
/*************************************************/
/**** CT EVOLUTION **********/
/*************************************************/
decrgpulse(pwC,t1,0.0,0.0);
delay((CT_delay-tau1)*0.25-pwC*2.0/3.0);
decpower(shlvl3);
decshaped_pulse(shname3,shpw3,zero,0.0,0.0);
delay((CT_delay+tau1)*0.25-shpw2*0.5);
shapedpulse(shname2,shpw2,zero,0.0,0.0);
delay((CT_delay+tau1)*0.25-shpw2*0.5);
decshaped_pulse(shname3,shpw3,zero,0.0,0.0);
decpower(pwClvl);
delay((CT_delay-tau1)*0.25-pwC*2.0/3.0);
decrgpulse(pwC,zero,0.0,0.0);
}
else
{
/*************************************************/
/**** REAL-TIME EVOLUTION **********/
/*************************************************/
if ((tau1) > shpw2)
{
decrgpulse(pwC,t1,0.0,0.0);
if(Cdecflg[0] == 'y')
{
decpower(getval("Cdecpwr"));
decprgon(Cdecseq, 1.0/getval("Cdecdmf"), getval("Cdecres"));
decon();
}
delay((tau1-shpw2)*0.5);
xmtrphase(zero);
xmtron();
obsunblank();
obsprgon(shname2,1/shdmf2,9.0);
delay(shpw2);
obsprgoff();
obsblank();
xmtroff();
delay((tau1-shpw2)*0.5);
if(Cdecflg[0] == 'y')
{
decoff(); decprgoff();
decpower(pwClvl);
}
decrgpulse(pwC,zero,0.0,0.0);
}
else
{
xmtrphase(zero);
xmtron();
obsunblank();
obsprgon(shname2,1/shdmf2,9.0);
delay((shpw2-tau1-pwC*2.0)*0.5);
decrgpulse(pwC,t1,0.0,0.0);
if(Cdecflg[0] == 'y')
{
decpower(getval("Cdecpwr"));
decprgon(Cdecseq, 1.0/getval("Cdecdmf"), getval("Cdecres"));
decon();
delay(tau1);
decoff(); decprgoff();
decpower(pwClvl);
}
else
delay(tau1);
decrgpulse(pwC,zero,0.0,0.0);
delay((shpw2-tau1-pwC*2.0)*0.5);
obsprgoff();
obsblank();
xmtroff();
}
}
obspower(shlvl1);
status(A);
zgradpulse(gzlvl1, gt1);
delay(1.0e-4);
delay(tauch-gt1-1.0e-4-POWER_DELAY);
decpower(dpwr);
status(C);
setreceiver(t2);
}
}
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();
}
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
//======================================================
// --------------------------------
// 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();
}
void pulsesequence()
{
/* DECLARE AND LOAD VARIABLES */
char
ch90shape[MAXSTR],
ch180shape[MAXSTR],
exp_mode[MAXSTR], /* flag to run 3D, or 2D time-shared 15N TROSY /13C HSQC-SE*/
decCACO[MAXSTR],
caco180shape[MAXSTR],
f1180[MAXSTR], /* Flag to start t1 @ halfdwell */
f2180[MAXSTR],
f3180[MAXSTR],
f4180[MAXSTR]; /* do TROSY on N15 and H1 */
int icosel, max_pcyc; /* used to get n and p type */
double
tpwrs,
ni2=getval("ni2"),
ni3=getval("ni3"),
tau1, tau1p,tau2,tau3,tau3p, /*evolution times in indirect dimensions */
tauNH=getval("tauNH"), /* 1/(4Jhn)*/
tauCH=getval("tauCH"), /* 1/(4Jch) */
tauCH1= getval("tauCH1"), /* tauCH/2.0+tauNH/2.0,*/ /* 1/(8Jch) +1/(8Jnh) */
tauCH2= getval("tauCH2"),
swC = getval("swC"), /* spectral widths in 13C methyls */
pwClvl = getval("pwClvl"), /* coarse power for C13 pulse */
pwC = getval("pwC"), /* C13 90 degree pulse length at pwClvl */
swN = getval("swN"), /* spectral widths in 15N */
pwNlvl = getval("pwNlvl"), /* power for N15 pulses */
pwN = getval("pwN"), /* N15 90 degree pulse length at pwNlvl */
ch90pwr=getval("ch90pwr"),
ch90pw=getval("ch90pw"),
ch90corr=getval("ch90corr"),
ch90dres=getval("ch90dres"),
ch90dmf=getval("ch90dmf"),
ch180pw=getval("ch180pw"),
ch180pwr=getval("ch180pwr"),
caco180pw=getval("caco180pw"),
caco180pwr=getval("caco180pwr"),
mix=getval("mix"),
tpwrsf_d = getval("tpwrsf_d"), /* fine power adustment for first soft pulse(down)*/
tpwrsf_u = getval("tpwrsf_u"), /* fine power adustment for second soft pulse(up) */
pwHs = getval("pwHs"), /* H1 90 degree pulse length at tpwrs */
compH =getval("compH"),
gstab = getval("gstab"),
gt0 = getval("gt0"),
gt1 = getval("gt1"),
gt2 = getval("gt2"),
gt3 = getval("gt3"),
gt4 = getval("gt4"),
gt5 = getval("gt5"),
gt6 = getval("gt6"),
gt7 = getval("gt7"),
gt8 = getval("gt8"),
gt9 = getval("gt9"),
gt10 = getval("gt10"),
gzlvl0 = getval("gzlvl0"),
gzlvl1 = getval("gzlvl1"),
gzlvl2 = getval("gzlvl2"),
gzlvl3 = getval("gzlvl3"),
gzlvl4 = getval("gzlvl4"),
gzlvl5 = getval("gzlvl5"),
gzlvl6 = getval("gzlvl6"),
gzlvl7 = getval("gzlvl7"),
gzlvl8 = getval("gzlvl8"),
gzlvl9 = getval("gzlvl9"),
gzlvl10 = getval("gzlvl10"),
gzlvl11 = getval("gzlvl11");
getstr("f1180",f1180);
getstr("f2180",f2180);
getstr("ch180shape",ch180shape);
getstr("ch90shape",ch90shape);
getstr("decCACO",decCACO);
getstr("caco180shape",caco180shape);
getstr("exp_mode",exp_mode);
tpwrs = tpwr - 20.0*log10(pwHs/(compH*pw*1.69)); /*needs 1.69 times more*/
tpwrs = (int) (tpwrs); /*power than a square pulse */
if (tpwrsf_d<4095.0)
tpwrs=tpwrs+6.0; /* add 6dB to let tpwrsf_d control fine power ~2048*/
if( (exp_mode[A]!='2') && (exp_mode[A]!='3') && (exp_mode[A]!='4') )
{text_error("invalid exp_mode, Should be either 2D or 3D or 4D\n "); psg_abort(1); }
/* LOAD PHASE TABLE */
settable(t1,1,phi1);
settable(t2,4,phi2);
settable(t12,4,phi2); {tsadd(t12,2,4);}
settable(t3,1,phi3);
settable(t4,2,phi4);
settable(t5,2,phi5);
settable(t6,4,phi6);
settable(t7,8,phi7);
settable(t8,8,phi8);
/* changing sign */
if( (exp_mode[A]=='4') && (exp_mode[C]=='a') )
{tsadd(t7,2,4); tsadd(t5,2,4); }
settable(t21,1,psi1); /*trosy and SE hsqc in reverse INPET */
settable(t22,1,psi2);
settable(t23,1,psi2c);
if(exp_mode[A]=='2') {settable(t31,2,rec);}
if(exp_mode[A]=='3') {settable(t31,4,rec);}
if(exp_mode[A]=='4') {settable(t31,8,rec);}
if((dm2[A] == 'y') || (dm2[B] == 'y') || (dm2[C] == 'y') || (dm2[D] == 'y'))
{ text_error("incorrect dec2 decoupler flags! Should be 'nnnn' "); psg_abort(1); }
/* special case for swapping t2 and t3 for test purposes */
if( (exp_mode[A]=='4') && (exp_mode[B]=='x') && (ni3=1) )
{
text_error("Acquiring t3 axis in ni2 dimension (instead of t2), set nt to 8! ");
tau3 = 0.5*(d3_index/swC+0.5/swC)-pw-rof1 -pwC*2.0/M_PI ; /* increment corresponds to 13C increment */
tau3p = 0.5*(d3_index/swN+0.5/swN) -pw-rof1 -pwC -pwN*2.0/M_PI -tau3;
if(d3_index % 2) { tsadd(t7,2,4); tsadd(t8,2,4); tsadd(t31,2,4); }
if (phase2 == 2) {tsadd(t7 ,1,4); tsadd(t8 ,1,4);}
tau2=0.0;
}
else {
if (phase2 == 2) {tsadd(t2 ,1,4); tsadd(t12,1,4);}
if (phase3 == 2) {tsadd(t7 ,1,4); tsadd(t8 ,1,4);}
if(d3_index % 2) { tsadd(t2,2,4); tsadd(t12,2,4); tsadd(t31,2,4); }
tau3 = 0.5*(d4_index/swC+0.5/swC)-pw -rof1 -pwC*2.0/M_PI ; /* increment corresponds to 13C increment */
tau3p = 0.5*(d4_index/swN+0.5/swN) -pw -rof1 -pwC -pwN*2.0/M_PI -tau3;
if(d4_index % 2) { tsadd(t7,2,4); tsadd(t8,2,4); tsadd(t31,2,4); }
tau2 = d3;
tau2 += 0.0*(-pw*4.0/M_PI-rof1*2.0);
if((f2180[A] == 'y') && (ni2 > 0.0)) {tau2 += ( 1.0 / (2.0*sw2) ); }
if(tau2 < 0.2e-6) {tau2 = 0.0;}
tau2 = tau2/2.0;
}
/* simultaneous Ntrosy-ChsqcSE, last part */
/* Phase incrementation for hypercomplex 2D data, States-Haberkorn element */
if (phase1 == 1) {icosel = 1; }
else { tsadd(t21,2,4); tsadd(t22,2,4); tsadd(t23,2,4); icosel = -1; }
if(d2_index % 2) { tsadd(t4,2,4); tsadd(t5,2,4); tsadd(t31,2,4); }
/* ECHO-ANTIECHO + STATES-TPPI t1, t1' in TROSY/HSQC last step */
tau1 = 1.0*d2_index/swC; /* increment corresponds to 13C increment */
tau1p = 1.0*d2_index*(1.0/swN-1.0/swC);
/* BEGIN PULSE SEQUENCE */
status(A);
obspower(tpwr);
decpower(pwClvl);
dec2power(pwNlvl);
txphase(zero);
decphase(zero);
dec2phase(zero);
delay(d1);
/* Destroy 13C magnetization*/
decrgpulse(pwC*1.0, zero, 0.0, 0.0);
zgradpulse(-gzlvl0, gt0);
delay(gstab);
/* NOESY */
if(exp_mode[A]!='2')
{ /* 3-4D */
if(exp_mode[A]=='4')
{ /* full 4D */
/* t3 evolution, the very first HSQC */
/* Hz -> HzXz INEPT */
rgpulse(pw,two,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl7, gt0); delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0); delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab); zgradpulse(gzlvl7, gt0); delay(gstab);
rgpulse(pw, one, rof1, rof1);
/* water defoc-refoc */
delay(gstab); zgradpulse(gzlvl8, gt8); delay(gstab);
/* t3 time */
if((ni3==0))
{
dec2rgpulse(pwN,t7,0.0,0.0);
dec2rgpulse(pwN,two,0.0,0.0);
decrgpulse(pwC, t8, 0.0, 0.0);
decrgpulse(pwC, two, 0.0, 0.0);
rgpulse(pw*2.0, zero, rof1, rof1);
delay(pwN*2.0+pwC*2.0);
}
else
{
dec2rgpulse(pwN,t7,0.0,0.0);
delay(tau3p);
decrgpulse(pwC, t8, 0.0, 0.0);
delay(tau3);
rgpulse(pw*2.0, zero, rof1, rof1);
delay(tau3);
decrgpulse(pwC, two, 0.0, 0.0);
delay(tau3p);
dec2rgpulse(pwN,two,0.0,0.0);
}
/* water defoc-refoc */
delay(gstab); zgradpulse(gzlvl8, gt8); delay(gstab);
/* back inept, water to +Z */
rgpulse(pw,one,rof1,rof1);
zgradpulse(gzlvl9, gt0); delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0); delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab); zgradpulse(gzlvl9, gt0); delay(gstab);
rgpulse(pw, two, rof1, rof1);
/* purge */
zgradpulse(gzlvl10, gt10); delay(2.0*gstab);
} /*end of full 4D */
/************************* t2 evolution, 1H and NOE****************************** */
/* for the case of no flipbacks in NOE part of the experiment, shift first pulse
in t2 time by 45 deg and let water bring itself back at the end
of mixing time by radiation dumping */
if( (exp_mode[D]=='t') )
{ initval(1.0, v10);
obsstepsize(45.0);
xmtrphase(v10);
}
else
{
xmtrphase(zero);
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc",pwHs,t12,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
}
rgpulse(pw,t2,rof1,rof1);
xmtrphase(zero); /* SAPS_DELAY */
delay(tau2);
decrgpulse(2.0*pwC,zero,0.0,0.0); dec2rgpulse(2.0*pwN,zero,0.0,0.0);
delay(tau2);
rgpulse(pw*2.0,zero,rof1,rof1);
delay(pwN*2.0+pwC*2.0 + SAPS_DELAY);
rgpulse(pw,zero,rof1,rof1);
if( (exp_mode[D]!='t') )
{
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc",pwHs,zero,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
}
/* NOESY period */
delay(mix-gt2-4.0*gstab );
zgradpulse(gzlvl2, gt2);
delay(4.0*gstab);
} /* end 3-4 D acquisition */
/* N-TROSY/C-HSQCse */
/* Hz -> HzXz INEPT */
rgpulse(pw,two,rof1,rof1); /* 1H pulse excitation */
zgradpulse(gzlvl0, gt0);
delay(tauCH-gt0);
decrgpulse(pwC*2.0, zero, 0.0, 0.0);
delay(tauNH -tauCH -pwC*2.0 );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
delay(tauNH - gt0 -gstab);
zgradpulse(gzlvl0, gt0);
delay(gstab);
rgpulse(pw, one, rof1, rof1);
/* on HzXz now */
/* water flipback*/
obspower(tpwrs); obspwrf(tpwrsf_u);
shaped_pulse("H2Osinc",pwHs,two,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
/* purge */
zgradpulse(gzlvl3, gt3);
dec2phase(t4);
delay(gstab*2.0);
/* t1 (C) and t1+t1'(N) evolution */
dec2rgpulse(pwN, t4, 0.0, 0.0);
delay(gt4+gstab + gt4+gstab + pwC*3.0
+2.0*(pwHs +2.0*rof1));
if(decCACO[A]=='y'){ delay(2.0*caco180pw);}
dec2rgpulse(2.0*pwN, zero, 0.0, 0.0);
delay(tau1p); /* t1 */
decrgpulse(pwC,t5,0.0,0.0);
delay(tau1*0.5);
if(decCACO[A]=='y')
{
decpower(caco180pwr);
decshaped_pulse(caco180shape,caco180pw,zero, 0.0, 0.0);
decpower(pwClvl);
}
obspower(ch180pwr); /*180 on methyls*/
shaped_pulse(ch180shape,ch180pw,zero,rof1,rof1);
obspower(tpwr);
delay(tau1*0.5);
zgradpulse(gzlvl4, gt4); /*coding */
delay(gstab + pwHs -ch180pw -2.0*GRADIENT_DELAY -2.0*POWER_DELAY -WFG_START_DELAY- WFG_STOP_DELAY);
decrgpulse(2.0*pwC,zero,0.0,0.0);
if(decCACO[A]=='y')
{
decpower(caco180pwr);
decshaped_pulse(caco180shape,caco180pw,zero, 0.0, 0.0);
decpower(pwClvl);
}
/* delay(ch180pw+2.0*rof1);*/
zgradpulse(gzlvl5, gt4);
delay(gstab - rof1 -2.0*GRADIENT_DELAY -2.0*POWER_DELAY -WFG_START_DELAY- WFG_STOP_DELAY);
/*Water flipback (flipdown actually ) */
obspower(tpwrs); obspwrf(tpwrsf_d);
shaped_pulse("H2Osinc",pwHs,three,rof1,rof1);
obspower(tpwr); obspwrf(4095.0);
/* reverse double INEPT */
sim3pulse(pw, pwC, 0.0, t21, t23, zero, rof1, rof1);
/* rgpulse(pw, t21, rof1, rof1); */
zgradpulse(gzlvl11, gt1);
delay(gstab);
delay(tauCH1 -gt1 -gstab -2.0*pwC
+ (-2.0/M_PI*pwC-0.5*(pwN-pwC) +pwN));
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -tauCH1 - 0.65*(pw + pwN)-rof1 -(pwC-pw)
-(-2.0/M_PI*pwC-0.5*(pwN-pwC) +pwN) );
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl11, gt1);
delay(gstab);
delay(tauCH1-gt1 -gstab-2.0*pwC);
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -1.3*pwN -tauCH1);
sim3pulse(pw, pwC, pwN, one, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(tauCH2-2.0*pwC-gt1-gstab);
decrgpulse(2.0*pwC,zero,0.0,0.0);
delay(tauNH -1.3*pwN-tauCH2);
sim3pulse(2.0*pw, 0.0, 2.0*pwN, zero, zero, zero, 0.0, 0.0);
zgradpulse(gzlvl1, gt1);
delay(gstab);
delay(tauNH-1.6*pwN -POWER_DELAY -ch90pw*ch90corr -gt1-gstab +pwN*0.5 -PRG_START_DELAY);
/* delay(ch90pw*0.5);
sim3pulse(0.0,0.0, pwN, one, zero, t22, 0.0, 0.0);
delay(ch90pw*0.5);*/
/* sim3shaped_pulse(ch90shape,"hard","hard",ch90pw,0.0, pwN, zero,zero, t22,0.0,0.0);*/
/* ch90corr is a fraction of ch90pw to correct for a 1H phase rollcaused by shaped 90 on
CH3 protons for a sinc pulse ch90corr=0.41 seems to be good. */
obspower(ch90pwr);
txphase(one);
obsunblank(); xmtron();
obsprgon(ch90shape,1.0/ch90dmf,ch90dres); /*PRG_START_DELAY */
delay(ch90pw*ch90corr-pwN*0.5);
dec2rgpulse(pwN, t22, 0.0, 0.0);
delay(ch90pw*(1.0-ch90corr)-pwN*0.5);
obsprgoff(); xmtroff(); obsblank(); /*PRG_STOP_DELAY */
obspower(tpwr); /*POWER_DEALY */
delay( gstab +gt6 +2.0*GRADIENT_DELAY
+2.0*POWER_DELAY -0.65*pw -POWER_DELAY
+pwN*0.5 -ch90pw*(1.0-ch90corr)
-PRG_STOP_DELAY);
rgpulse(2.0*pw, zero, rof1, rof1);
dec2power(dpwr2); decpower(dpwr); /* 2.0*POWER_DELAY */
zgradpulse(gzlvl6*icosel, gt6); /* 2.0*GRADIENT_DELAY */
delay(gstab);
status(C);
setreceiver(t31);
}
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");
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 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
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
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();
}
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();
}
