# We set T, T', S, and S'. The T'=dT/d(kappa) and S'=dS/d(kappa). # where kappa = 2*PI*frequency/(c*beta) # The T' and S' are set up as separate polynomials, because # the accuracy of calculating a derivative from the polynomial # fitting is very low. #--------------------------------------- rf_gap_ttf_arr = [] for i_gap in range(n_gaps): rf_gap_ttf = RfGapTTF() polyT = split_polynom_coeffs(lns[8+i_gap]) polyTp = split_polynom_coeffs(lns[9+i_gap]) polyS = split_polynom_coeffs(lns[10+i_gap]) polySp = split_polynom_coeffs(lns[11+i_gap]) rf_gap_ttf.setT_TTF(polyT) rf_gap_ttf.setTp_TTF(polyTp) rf_gap_ttf.setS_TTF(polyS) rf_gap_ttf.setSp_TTF(polySp) gap_length = gap_lengths[i_gap] relative_amplitude = gap_E0_amplitudes[i_gap] rf_gap_ttf.setParameters(polyT,polyTp,polyS,polySp,beta_min,beta_max,rf_freq,gap_length,relative_amplitude) rf_gap_ttf_arr.append(rf_gap_ttf) #--------directions of the cavity can be +1 or -1 directionZ = +1 if(directionZ < 0): for i_gap in range(len(gap_border_points)): gap_border_points[i_gap] = - gap_border_points[i_gap] for i_gap in range(len(gap_positions)): gap_positions[i_gap] = - gap_positions[i_gap] gap_border_points.reverse()
while(1 < 2): count += 1 rf_gap_ttf = RfGapTTF() polyT = Polynomial(4) polyT.coefficient(2,2.0) T_ttf = rf_gap_ttf.setT_TTF(polyT) polyS = Polynomial(5) polyS.coefficient(3,3.0) S_ttf = rf_gap_ttf.setS_TTF(polyS) polyT = Polynomial(4) polyT.coefficient(2,2.0) T_ttf = rf_gap_ttf.setTp_TTF(polyT) polyS = Polynomial(5) polyS.coefficient(3,3.0) S_ttf = rf_gap_ttf.setSp_TTF(polyS) T_ttf = rf_gap_ttf.getT_TTF() S_ttf = rf_gap_ttf.getS_TTF() Tp_ttf = rf_gap_ttf.getTp_TTF() Sp_ttf = rf_gap_ttf.getSp_TTF() if(count % 100000 == 0): print "count=",count
# We set T, T', S, and S'. The T'=dT/d(kappa) and S'=dS/d(kappa).
# where kappa = 2*PI*frequency/(c*beta)
# The T' and S' are set up as separate polynomials, because
# the accuracy of calculating a derivative from the polynomial
# fitting is very low.
#---------------------------------------
rf_gap_ttf_arr = []
for i_gap in range(n_gaps):
rf_gap_ttf = RfGapTTF()
polyT = split_polynom_coeffs(lns[8 + i_gap])
polyTp = split_polynom_coeffs(lns[9 + i_gap])
polyS = split_polynom_coeffs(lns[10 + i_gap])
polySp = split_polynom_coeffs(lns[11 + i_gap])
rf_gap_ttf.setT_TTF(polyT)
rf_gap_ttf.setTp_TTF(polyTp)
rf_gap_ttf.setS_TTF(polyS)
rf_gap_ttf.setSp_TTF(polySp)
gap_length = gap_lengths[i_gap]
relative_amplitude = gap_E0_amplitudes[i_gap]
rf_gap_ttf.setParameters(polyT, polyTp, polyS, polySp, beta_min, beta_max,
rf_freq, gap_length, relative_amplitude)
rf_gap_ttf_arr.append(rf_gap_ttf)
#--------directions of the cavity can be +1 or -1
directionZ = +1
if (directionZ < 0):
for i_gap in range(len(gap_border_points)):
gap_border_points[i_gap] = -gap_border_points[i_gap]
for i_gap in range(len(gap_positions)):
gap_positions[i_gap] = -gap_positions[i_gap]
rf_gap_ttf = RfGapTTF()
while (1 < 2):
count += 1
rf_gap_ttf = RfGapTTF()
polyT = Polynomial(4)
polyT.coefficient(2, 2.0)
T_ttf = rf_gap_ttf.setT_TTF(polyT)
polyS = Polynomial(5)
polyS.coefficient(3, 3.0)
S_ttf = rf_gap_ttf.setS_TTF(polyS)
polyT = Polynomial(4)
polyT.coefficient(2, 2.0)
T_ttf = rf_gap_ttf.setTp_TTF(polyT)
polyS = Polynomial(5)
polyS.coefficient(3, 3.0)
S_ttf = rf_gap_ttf.setSp_TTF(polyS)
T_ttf = rf_gap_ttf.getT_TTF()
S_ttf = rf_gap_ttf.getS_TTF()
Tp_ttf = rf_gap_ttf.getTp_TTF()
Sp_ttf = rf_gap_ttf.getSp_TTF()
if (count % 100000 == 0): print "count=", count