;2.3) Program: 	HETCOR after 1H CSA filter
;Description: 	Suppresses NH/OH signals in HETCOR spectra (overlap with aromatic H)
;Requires: 	!!!!!;  
;Reference(s):	KSR & Mao, J. Magn. Reson. 157, 210-217 (2002).

;srhetcor.hcsafA (XwinNMR 2.6)
;from srhetcor1D.hcsafA
;from srhetcor.mh
; lg-480 consists of 480 t(360), with 22 steps in t(360)
; total duration is p10 = 7920 u
; 7920 u /480 = 16.5 u = t(360) 
; good because 16.5 / 22 = 0.75 u sharp, so digital increment will be exact
; t(360)/sin(54.7)/4 = t(90 true) = 5.05u, set sp1 accordingly

; lg-480 consists of 480 t(360), with 22 steps in t(360)
; total duration is p10 = 6336 u
; 6336 u /480 = 13.2 u = t(360) 
; good because 13.2 / 22 = 0.6 u sharp, so digital increment will be exact
; t(360)/sin(54.7)/4 = t(90 true) = 4.04u, set sp1 accordingly

;States-TPPI detected 1H-13C HETCOR, 
;with LG-CP to suppress spin diffusion
;MH 8-8-00

;p10: t1, set as t360*960 for lgs-960
;p22: 35.3 deg pulse
;p21: 54.7 deg pulse
;p9: t1 dwell time
;p3: 1H excitation
;p8: TPPM pulse
;p15: CP, short
;spnam1: homonuclear decoupling sequence
;sp1: FSLG decoupling during t1
;pl2: 1H excitation
;pl3: 1H LG-CP, ca. 1.8dB > HHCP
;pl1: 13C CP, fixed amplitude
;pl12: TPPM decoupling during t2

define delay swi
"swi=d23-p21"
define delay lwi
"lwi=d23+d23-p21"
define delay swix   ;at 13C 180 pulse
"swix=lwi-swi-p2/2"

;CANNOT do RAMP CP on 13C. Signals go away!!

  "d8=p21"	;HW-8 pulse
  "d25=(1s/l31)*0.1226-(p2/2)-2u"
  "d26=(1s/l31)*0.0773-p2"
  "d27=(1s/l31)*0.2236-p2"
  "d28=(1s/l31)*1.0433-p2"
  "d29=(1s/l31)*0.7744-(p2/2)-d13"

define loopcounter td1half
define delay tr
  "td1half=td1/2"			;for STATES-TPPI
  "p22=(p3*353)/900"
  "p23=(p3*547)/900"		;magic-angle pulse
  "l1=aq/(p8*2)"
  "l13=l3/2"
  "l23=l13-1"
  "tr=48*d23" ;-swi-p21+p1+p2/2"

  100m
  "p10=p9;"       

1 ze                             ;zero data
  tr   				 ;just to get tr calculated

;2D acquisition
12 d1 do:f2                      ;recycle time, decoupler off
  10u pl11:f1 			;set power level of f1
  p1:f1 ph9			;destroy 13C magnetization
  5m
  p1:f1 ph10
  5m
  p1:f1 ph9			;destroy 13C magnetization
  5m
  p1:f1 ph10
  5m pl1:f1

  10u fq1:f2		;slightly off resonance for CSA filter

  10u pl22:f2		;high power (e.g.6 dB) for HW8 selection

6  swi
  p21:f2 ph22 ;-x
  swi
  p21:f2 ph21  ;y
  lwi
  p21:f2 ph23  ;-y
  swi
  p21:f2 ph22  ;-x
  lwi
  p21:f2 ph20  ;x
  swi
  p21:f2 ph23  ;-y
  lwi
  p21:f2 ph21  ;y
  swi
  p21:f2 ph20  ;x
  swi
  d8 
 
  swi
  p21:f2 ph22 ;-x
  swi
  p21:f2 ph21  ;y
  lwi
  p21:f2 ph23  ;-y
  swi
  p21:f2 ph22  ;-x
  lwi
  p21:f2 ph20  ;x
  swi
  p21:f2 ph23  ;-y
  lwi
  p21:f2 ph21  ;y
  swi
  p21:f2 ph20  ;x
  swix
  
;tr/2

  (p2 ph18):f1  ;180 13C pulse for decoupling ------------------------

  swix
  p21:f2 ph20 ;x invert all pulses  
  swi
  p21:f2 ph23  ;-y
  lwi
  p21:f2 ph21  ;y
  swi
  p21:f2 ph20  ;x	
  lwi
  p21:f2 ph22  ;-x
  swi
  p21:f2 ph21  ;y
  lwi
  p21:f2 ph23  ;-y
  swi
  p21:f2 ph22  ;-x
  lwi 

  p21:f2 ph20 ;x invert all pulses  
  swi
  p21:f2 ph23  ;-y
  lwi
  p21:f2 ph21  ;y
  swi
  p21:f2 ph20  ;x	
  lwi
  p21:f2 ph22  ;-x
  swi
  p21:f2 ph21  ;y
  lwi
  p21:f2 ph23  ;-y
  swi
  p21:f2 ph22  ;-x
  swi		
  
  lo to 6 times l0 		;end of CSA filter

  1u fq2:f2                      ;small offset during evolution
  2u pl2:f2
;  1u pl1:f1                     ; set power level of f1

16 p3:f2 ph1			;excitation pulse 
  p23:f2 ph3       		;55 deg -y pulse,normal to eff field 

;-------------- t1 evolution  ---------------------------------
  (p10:sp1 ph4):f2		;FSLG, 1H CS evolution
  1u pl2:f2 fq4:f2            ;reset O2= on resonance

;  (p22 ph7):f2 		;35 deg -y pulse. to -z
 
  (p3 ph8):f2			;90 deg +-y pulse, to xy plane
;  (p22 ph5):f2			;35 deg pulse, +y
  1u fq3:f2 pl3:f2		;fq3=ca.-30 kHz
17 (p15 ph2):f1 (p15 ph11):f2	;LG-CP

  1u fq4:f2
  1u cw:f2 pl12:f2

   d25 pl11:f1           
  p2:f1 ph27             ;180 13C #1
   d26             
  p2:f1 ph24             ;180 13C #2
   d27           
  p2:f1 ph25             ;180 13C #3
   d28                  
  p2:f1 ph26             ;180 13C #4
   d29 pl12:f2         


  1u:f1 ph0			;dead time delay, reset RF phase to detect
  2u adc ph31			;start ADC with ph31 signal routing
18 p8:f2 ph13			;TPPM decoupling
  p8:f2 ph14
lo to 18 times l1
  aq do:f2			
  rcyc=12 			;do NS scans

;1D:
;wr #0

;2D:
19  100m wr #0 if #0 ip1 zd	;write  and increment file
lo to 12 times 2		
 "p10=p10+p13;"			;calculate next pulse length for LG evolution
  1m
lo to 12 times td1half		;number of rows

exit                       

ph0=0
ph1=0			;excitation, States-TPPI
;ph1=1			;cosine tests
ph3=3			;1H 55 deg pulse
ph4=0 			;1H FSLG phase
ph5=1			;1H 35 deg pulse, +y
ph6=2
ph7=3 			;1H 35 deg. pulse
ph8=1 3			;1H 90 deg. pulse

ph9=+x			;destroy 13C magnetization
ph10=+y			;destroy 13C magnetization

ph11= 0			;1H CP
ph2= +x +x +y +y -x -x -y -y;  0 0 2 2 1 1 3 3  
ph13=(360) 0		;TPPM1
ph14=(360) 10		;TPPM2

ph18=+x +x +x +x +x +x +x +x	;180 on 13C in center of 1H CSA filter
     -x -x -x -x -x -x -x -x
     +y +y +y +y +y +y +y +y
     -y -y -y -y -y -y -y -y

ph19=+y +y +y +y +y +y +y +y    ;90 in compos. 180 on 13C in center of 1H CSA filter
     -y -y -y -y -y -y -y -y
     +x +x +x +x +x +x +x +x	
     -x -x -x -x -x -x -x -x 
ph20=+x			;for HW-8 CSA filter on 1H (M ends up along z, so x vs 0 doesn't matter)
ph21=+y
ph22=-x
ph23=-y

ph27= +y +y -x -x -y -y +x +x    ;180 13C #1
ph24= -y -y +x +x +y +y -x -x    ;180 13C #2
ph25= +y +y -x -x -y -y +x +x    ;180 13C #3
     -y -y +x +x +y +y -x -x
ph26= -x -x +y +y +x +x -y -y    ;180 13C #4
     +x +x -y -y -x -x +y +y

ph31=0 2 1 3 2 0 3 1 ; Receiver