def main():
inFileName = sys.argv[1] # input .mco file
T = float(sys.argv[2]) # truncation thresholf for ISP
N = int(sys.argv[3]) # truncation parameter for FB expansion
R0 = float(sys.argv[4]) # width of the top-hat part
A0 = float(sys.argv[5]) # decay width of Gaussian
R1 = float(sys.argv[6]) # width of the top-hat part
A1 = float(sys.argv[7]) # decay width of Gaussian
# FETCH IMPULSE RESPONSE FROM mco. FILE
r,z,Arz = io.fetchRawData(inFileName)
# CORRECT AXES TO REFER TO BIN-CENTERS: REF [1], EQS. (4.11), (4.19)
r = r + 0.5*(r[1]-r[0])
z = z + 0.5*(z[1]-z[0])
# CORRECT FOR ACCUMULUATED VOID WALKERS: REF [2], EQ. (22)
Arz[-1:,:]=0
# SET CUSTOM IRRADIATION SOURCE PROFILE FROM OSG LIBRARY
myIsp = lambda x: donut(x, (R0,R1,A0,A1))
# CONVOLVE IMPULSE RESPONSE USING CUSTOM IRRADIATION SOURCE PROFILE
Wrz,recErr = conv.FiskJohnsonConvolvedResponse(r, z, Arz, myIsp, (T,N))
print "# (iz) (zi) (rj) (W(rj,zi))"
for idx,iz in enumerate([int(zi/(z[1]-z[0])) for zi in [0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8]]):
for ir in range(r.size):
print "z%d"%(idx),z[iz],r[ir],Wrz[ir,iz]
print
print "# (ir) (zi) (rj) (W(rj,zi))"
for idx,ir in enumerate([int(ri/(r[1]-r[0])) for ri in [0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8]]):
for iz in range(z.size):
print "r%d"%(idx),z[iz],r[ir],Wrz[ir,iz]
print
def main():
inFileName = sys.argv[1] # input .mco file
outFileName = sys.argv[8] # output .prof file
T = float(sys.argv[2]) # truncation thresholf for ISP
N = int(sys.argv[3]) # truncation parameter for FB expansion
R0 = float(sys.argv[4]) # width of the top-hat part
A0 = float(sys.argv[5]) # decay width of Gaussian
R1 = float(sys.argv[6]) # width of the top-hat part
A1 = float(sys.argv[7]) # decay width of Gaussian
# FETCH IMPULSE RESPONSE FROM mco. FILE
r, z, Arz = io.fetchRawData(inFileName)
# CORRECT AXES TO REFER TO BIN-CENTERS: REF [1], EQS. (4.11), (4.19)
r = r + 0.5 * (r[1] - r[0])
z = z + 0.5 * (z[1] - z[0])
# CORRECT FOR ACCUMULUATED VOID WALKERS: REF [2], EQ. (22)
Arz[-1:, :] = 0
# SET CUSTOM IRRADIATION SOURCE PROFILE FROM OSG LIBRARY
myIsp = lambda x: donut(x, (R0, R1, A0, A1))
# CONVOLVE IMPULSE RESPONSE USING CUSTOM IRRADIATION SOURCE PROFILE
Wrz, recErr = conv.FiskJohnsonConvolvedResponse(r, z, Arz, myIsp, (T, N))
# SAVE DATA IN GNUPLOT FORMAT
io.writeGpl(r, z, Wrz, outFileName)
f = open(outFileName, 'a')
f.write("# ISPRecErr: %lf\n" % (recErr))
f.write("# maxW: %lf\n" % (Wrz.max()))
f.close()
def main():
inFileName = sys.argv[1] # input .mco file
outFileName = sys.argv[2] # output .prof file
# FETCH IMPULSE RESPONSE FROM mco. FILE
r, z, Arz = io.fetchRawData(inFileName)
# ACCOUNT FOR VOID-WALKER
Arz[-1, :] = 0.
# SAVE DATA IN GNUPLOT FORMAT
io.writeGpl(r, z, Arz, outFileName)
def main():
inFileName = sys.argv[1]
outFileName = "testFJ"
T = None # Truncation radius for beam profile
N = 32 # number of samples to consider during FB transform
# FETCH IMPULSE RESPONSE FROM mco. FILE
r, z, Arz = io.fetchRawData(inFileName)
# ACCOUNT FOR VOID-WALKER
Arz[-1, :] = 0.
# SET CUSTOM IRRADIATION SOURCE PROFILE FROM OSG LIBRARY
myIsp = lambda x: isp.flatTop(x, (0.05, 0.05))
# CONVOLVE IMPULSE RESPONSE USING CUSTOM IRRADIATION SOURCE PROFILE
Wrz = conv.FiskJohnsonConvolvedResponse(r, z, Arz, myIsp, (T, N))
# SAVE DATA IN GNUPLOT FORMAT
io.writeGpl(r, z, Wrz, outFileName)