def HDec(units): Stuff = [] p90 = pul.GetPar('pH90', "") amp = pul.GetPar('aH', units) CPD = pul.GetPar('prgHDec', "") MaxB1 = 1000000. / 4. / p90 if CPD == "mlev" or CPD == "None": pul.SetPar('prgHDec', "tppm15", "") TopCmds.XCMD(pul.xcmd_name(pul.pulDict['prgHDec'])) CPD = pul.GetPar('prgHDec', "") TopCmds.MSG(str(CPD) + " CPD") Stuff = CPDtools.CPDparse(CPD, "1H") TopCmds.MSG(str(Stuff) + " Stuff") amp0 = CPDtools.Find_old_pl(Stuff[0], units) TopCmds.MSG(str(amp0) + " amp0") if units == "W": amp = WtodB(amp) amp0 = WtodB(amp0) decpw0 = CPDtools.Find_old_pw(Stuff[1], "1H") B1_0 = MaxB1 * (math.pow(10, (amp - amp0) / 20.)) / 1000. if B1_0 > 1.: B1out = '% .1f' % B1_0 if B1_0 > MaxB1 / 1000.: B1out = '75.0' if B1_0 <= 1.: B1out = '75.0' index = TopCmds.INPUT_DIALOG("Mr Setup Input", "Decoupling Window", \ ["Desired 1H Decoupling Amplitude","File"],\ [B1out,CPD],["kHz",""],["1","1",],\ ["Accept","Close"], ['a','c'], 10) pul.SetPar('prgHDec', index[1], "") adjust = 20 * (math.log10(1000. * float(index[0]) / MaxB1)) Hamp = amp - adjust decpwH = (MaxB1 / 1000. / float(index[0])) * (170. / 180.) * 2 * p90 if units == "W": Hamp = dBtoW(Hamp) value = TopCmds.SELECT("Adjusting the H decoupling parameters:",\ "This will set\n 1H power ("+ Stuff[0] +") to: "+ str('%.2f' %Hamp)+" "+ units+"\n \ Pulse width (" + Stuff[1] +") to: " +str('%3.2f' %decpwH)+" us",["Update", "Keep Previous"]) if value != 1: pl = "" if Stuff[0] != "": pl = pul.pp_2_xcmd(Stuff[0], "") if pl == pul.pulDict['aHdec']: pul.SetPar('aHdec', Hamp, units) elif pl == pul.pulDict['aHdec2']: pul.SetPar('aHdec2', Hamp, units) elif pl == pul.pulDict['aHdec3']: pul.SetPar('aHdec3', Hamp, units) if Stuff[1] == 'pcpd': pul.SetPar("PCPD 2", decpwH, "") elif Stuff[1] == 'p31': pul.SetPar("P31", decpwH, "") pul.SetPar("P30", decpwH, "") elif Stuff[1] == 'p62': pul.SetPar("P61", decpwH, "") pul.SetPar("P62", decpwH, "")
def DDec(units): Stuff = [] p90 = pul.GetPar('pD90', "") amp = pul.GetPar('aD', units) CPD = pul.GetPar('prgDDec', "") MaxB1 = 1000000. / 4. / p90 if units == "W": amp = WtodB(amp) if CPD == "mlev" or CPD == "None": TopCmds.XCMD(pul.xcmd_name(pul.pulDict['prgDDec'])) CPD = pul.GetPar('prgDDec', "") Stuff = CPDtools.CPDparse(CPD, "2H") amp0 = CPDtools.Find_old_pl(Stuff[0]) decpw0 = CPDtools.Find_old_pw(Stuff[1], "2H") B1_0 = MaxB1 * (math.pow(10, (amp - amp0) / 20.)) / 1000. if B1_0 > 1.: B1out = '% .1f' % B1_0 if B1_0 <= 1.: B1out = '15.0' index = TopCmds.INPUT_DIALOG("Mr Setup Input", "Decoupling Window", \ ["Desired 2H Decoupling Amplitude","File"],\ [B1out,CPD],["kHz",""],["1","1",],\ ["Accept","Close"], ['a','c'], 10) pul.SetPar('aDdec', index[1], "") adjust = 20 * (math.log10(1000. * float(index[0]) / MaxB1)) ampli = amp - adjust decpw = (MaxB1 / 1000. / float(index[0])) * 2 * p90 if units == "W": ampli = dBtoW(ampli) value = TopCmds.SELECT("Adjusting the 2H decoupling parameters:",\ "This will set\n 2H power ("+ Stuff[0] +") to: "+ str('%3.2f' %ampli)+" "+ units+"\n \ Pulse width (" + Stuff[1] +"= 180deg) to: " +str('%3.2f' %decpw) + " us",["Update", "Keep Previous"]) if value != 1: pl = "" if Stuff[0] != "": pl = pul.pp_2_xcmd(Stuff[0], "") if pl == "": pul.SetPar('aDdec', ampli, units) elif pl == pul.pulDict['aDdec']: pul.SetPar('aDdec', ampli, units) elif pl == pul.pulDict['aHdec']: DecError('aHdec', "13C") elif pl == pul.pulDict['aHdec2']: DecError('aHdec2', "13C") elif pl == pul.pulDict['aHdec3']: DecError('aHdec3', "13C") if Stuff[1] == 'pcpd': pul.SetPar("PCPD 3", decpw, "") elif Stuff[1] == 'p31': pul.SetPar("P31", decpw, "") pul.SetPar("P30", decpw, "") elif Stuff[1] == 'p62': pul.SetPar("P61", decpw, "") pul.SetPar("P62", decpw, "")
def CalDec(p90, amp, cpd, nuc, units, dfltWave, dfltB0, dfltPH, aOption, aFaults): """ p90 : Dict key for Hard Pulse of Decoupled Nucleus amp : Dict key for Hard Pulse Amplitude cpd : Dict key for CPD file nuc : Decoupled Nucleus units : Watts (W) or Decibels (dB) dfltWave: Default CPD dfltB0 : Default field dfltPH : Default tip angle used in CPD aOption : List of accepted amplitude dict keys aFaults : List of amplitude dict keys that will cause a PLEASE CONFIRM Message """ Stuff = [] P90 = pul.GetPar(p90, "") Amp = pul.GetPar(amp, units) CPD = pul.GetPar(cpd, "") MaxB1 = 1000000. / 4. / P90 if CPD == "mlev" or CPD == "None" or CPD == None or CPD == "": pul.SetPar(cpd, dfltWave, "") TopCmds.XCMD(pul.xcmd_name(pul.pulDict[cpd])) CPD = pul.GetPar(cpd, "") Stuff = CPDtools.CPDparse(CPD, nuc) Amp0 = CPDtools.Find_old_pl(Stuff[0], units) if units == "W": Amp = WtodB(Amp) Amp0 = WtodB(Amp0) decpw0 = CPDtools.Find_old_pw(Stuff[1], nuc) B1_0 = MaxB1 * (math.pow(10, (Amp - Amp0) / 20.)) / 1000. if B1_0 > 1.: B1out = '% .3f' % B1_0 if B1_0 > MaxB1 / 1000.: B1out = '% .3f' % dfltB0 if B1_0 <= 1.: B1out = '% .3f' % dfltB0 index = TopCmds.INPUT_DIALOG("Mr Setup Input", "Decoupling Window", \ ["Desired "+nuc+" Decoupling Amplitude","File"],\ [B1out,CPD],["kHz",""],["1","1"],\ ["Accept","Close"], [spc,ret], 10) if index == None: TopCmds.EXIT() pul.SetPar(cpd, index[1], "") pl = pul.pp_2_xcmd(Stuff[0], "") matched = 0 for a in aOption: if pl == pul.pulDict[a]: matched = 1 Hamp = DecSafely( float(index[0]) * 1000, a, MaxB1, 2 * 1000. * dfltB0, Amp, units) if matched == 0: Hamp = DecSafely( float(index[0]) * 1000, aOption[0], MaxB1, 2 * 1000. * dfltB0, Amp, units) decpw = (MaxB1 / 1000. / float(index[0])) * (dfltPH / 180.) * 2 * P90 if units == "W": Hamp = dBtoW(Hamp) value = TopCmds.SELECT("Adjusting the "+nuc+" decoupling parameters:",\ "This will set\n "+nuc+" power ("+ Stuff[0] +") to: "+ str('%.2f' %Hamp)+" "+ units+"\n \ Pulse width (" + Stuff[1] +") to: " +str('%3.2f' %decpw)+" us",["Update", "Keep Previous"]) if value != 1: if Stuff[0] == "": pul.SetPar(aOption[0], Hamp, units) for i in aFaults: if pl == pul.pulDict[i]: DecError(f, nuc) for a in aOption: if pl == pul.pulDict[a]: pul.SetPar(a, Hamp, units) return Stuff[1], decpw