def make(Name, Sc, Stps, Del, Be):
import math
ampl = [] # normalized to 0...100
ph = [] # normalized to 0...100
Scale = float(Sc)
steps = int(Stps)
Delta = float(Del)
Beta = float(Be)
pi = 3.14159265
k = (steps - 1) / 2.0
if Beta == 0.:
alpha = 0. #No division by Zero
else:
alpha = (2.0 / (steps - 1)) * math.atan(Delta / Beta)
absBeta = math.fabs(Beta)
for i in range(steps):
di = absBeta * math.tan(alpha * (float(i - k)))
di = di
ampl.append(Scale * (100. + di) / 100.)
ph.append(0.0)
TopCmds.SAVE_SHAPE(Name, "NoRotation", ampl, ph)
def make(Scale, input, name):
import math
ampl = [] # normalized to 0...100
ph = [] # normalized to 0...100
pi = 3.14159265
durat = 1000 * float(input[2])
RIRO = float(input[3])
steps = int(input[4])
Delta = float(input[5])
RIOsteps = int(steps * RIRO / durat)
Start = 100.0 - (Delta)
End = 100.0 + (Delta)
for i in range(steps):
if i < RIOsteps:
RF = 1.0 * i * Start / RIOsteps
if i >= RIOsteps:
RF = 1.0 * Start + (i - RIOsteps) * (End -
Start) / (steps -
(2 * RIOsteps))
if i > (steps - RIOsteps):
RF = 1.0 * End - (i - steps + RIOsteps) * End / RIOsteps
ampl.append(Scale * RF)
ph.append(0.0)
TopCmds.SAVE_SHAPE(name[0], "NoRotation", ampl, ph)
def make(input, name):
import math
amplitudes = [] # normalized to 0...100
phases = [] # in degrees
pi = 3.14159265
NomRF = float(input[0])
MAS = float(input[1])
steps = int(input[2])
DUO_B = MAS * float(input[3])
DUO_C = MAS * float(input[4])
duration = 2.0e06 / MAS
deltaT = duration / steps
RFmax = 0
for i in range(steps):
ti = i * deltaT
co = math.cos(2 * pi * DUO_C * ti * 1.0e-6)
si = math.sin(2 * pi * DUO_C * ti * 1.0e-6)
# Calculate RF portion
rf1 = DUO_C
rf2 = -1.0 * DUO_B * si
rf = math.sqrt(rf1 * rf1 + rf2 * rf2)
if rf > RFmax:
RFmax = rf
# Calculate Phase Portion
phcor = -1 * DUO_B * si / DUO_C
ph = (math.atan2(rf2, rf1) + phcor) * 180 / pi
if ph >= 0.0:
nph = int(ph / 360.0)
phint = ph - (nph) * 360
else:
nph = int(ph / 360.0)
phint = ph - (nph - 1) * 360
if RFmax > (NomRF * 1000):
TopCmds.ERRMSG(
message="ERROR the maximum RF is above the nominal RF",
title="WAVE ERROR")
amplitudes.append(rf / NomRF / 10)
phases.append(phint)
TopCmds.SAVE_SHAPE(name, "Excitation", amplitudes, phases)
return
def make(Ttot, RIRO, Sc, Del, Bet, Name):
import math
ampl = [] # normalized to 0...100
ph = [] # normalized to 0...100
#TopCmds.MSG(Sc+" "+Del+" "+Bet)
durat = 1000. * float(Ttot)
steps = int(1000 * float(Ttot))
Read = float(RIRO)
Delta = float(Del)
Beta = float(Bet)
Scale = float(Sc)
RIOsteps = int(steps * Read / durat)
TopCmds
tansteps = steps - 2 * RIOsteps
Start = 100. - (Delta * (Beta / math.fabs(Beta)))
End = 100. + (Delta * (Beta / math.fabs(Beta)))
k = (tansteps - 1) / 2
alpha = (2.0 / (tansteps - 1)) * math.atan(Delta / Beta)
absBeta = math.fabs(Beta)
for i in range(steps):
if i < RIOsteps:
RF = 1.0 * i * Start / RIOsteps
if i >= RIOsteps:
j = i - RIOsteps
di = absBeta * math.tan(alpha * float(j - k))
RF = 100. + di
if i > (steps - RIOsteps):
j = i - steps + RIOsteps
RF = 1.0 * End - j * End / RIOsteps
ampl.append(Scale * RF / 100)
ph.append(0.0)
TopCmds.SAVE_SHAPE(Name[0], "NoRotation", ampl, ph)
def make(input, Names):
import math
Iampl = [] # normalized to 0...100
Sampl = [] # normalized to 0...100
Iphases = [] # in degrees
Sphases = [] # in degrees
Names = []
iph = [] # in degrees
pi = 3.14159265
NomRFI = float(input[0])
NomRFS = float(input[1])
MAS = float(input[2])
steps = int(input[3])
CF = MAS * float(input[4])
BI = MAS * float(input[5])
BS = MAS * float(input[6])
Iphase = float(input[7])
duration = 2.0e06 / MAS
deltaT = duration / steps
RFImax = 0
RFSmax = 0
h = 0 - 1
for i in range(steps):
j = i + 1
k = 1 + math.fmod(j - 1, steps)
if (k == 1): h = -1 * h
ti = i * deltaT
co = math.cos(2 * pi * CF * ti * 1.0e-6)
si = math.sin(2 * pi * CF * ti * 1.0e-6)
#RF part
RFIx = CF * (h)
RFIy = BI * co
RFSx = CF * (h)
RFSy = BS * co
RFI = math.sqrt(RFIx * RFIx + RFIy * RFIy)
RFS = math.sqrt(RFSx * RFSx + RFSy * RFSy)
if RFI > RFImax:
RFImax = RFI
if RFS > RFSmax:
RFSmax = RFS
if RFImax > (NomRFI * 1000):
TopCmds.ERRMSG(message="ERROR the maximum RF is above the nominal RF",\
title="WAVE ERROR")
if RFSmax > (NomRFS * 1000):
TopCmds.ERRMSG(message="ERROR the maximum RF is above the nominal RF",\
title="WAVE ERROR")
#Phase Part
Iphcor = Iphase + BI * h * (co - 1) / (CF)
Sphcor = BS * h * (co - 1) / (CF)
Iph = (Iphcor + math.atan2(RFIy, RFIx)) * 180 / pi
Sph = (Sphcor + math.atan2(RFSy, RFSx)) * 180 / pi
if (Iph >= 0.0):
Inph = int(Iph / 360.0)
Iph_int = Iph - (Inph) * 360
else:
Inph = int(Iph / 360.0)
Iph_int = Iph - (Inph - 1) * 360
if (Sph >= 0.0):
Snph = int(Sph / 360.0)
Sph_int = Sph - (Snph) * 360
else:
Snph = int(Sph / 360.0)
Sph_int = Sph - (Snph - 1) * 360
Iampl.append(RFI / NomRFI / 10)
Sampl.append(RFS / NomRFS / 10)
Iphases.append(Iph_int)
Sphases.append(Sph_int)
TopCmds.SAVE_SHAPE(Names[0], "Excitation", Iampl, Iphases)
TopCmds.SAVE_SHAPE(Names[1], "Excitation", Sampl, Sphases)