def main(visit=131634,
ccd=None,
threshold=45000,
nSample=1,
showCoeffs=True,
fixXTalk=True,
plot=False,
title=None):
if ccd is None:
visitList = list(range(nSample))
ccdList = [
"simulated",
]
else:
ccdList = makeList(ccd)
visitList = makeList(visit)
coeffs = estimateCoeffs(visitList,
ccdList,
threshold=45000,
plot=plot,
title=title)
if showCoeffs:
printCoeffs(coeffs)
mi = readImage(visitList[0], ccdList[0])
if fixXTalk:
subtractXTalk(mi, coeffs, threshold)
return mi, coeffs
def main(visit=131634, ccd=None, threshold=45000, nSample=1, showCoeffs=True, fixXTalk=True,
plot=False, title=None):
if ccd is None:
visitList = list(range(nSample))
ccdList = ["simulated", ]
else:
ccdList = makeList(ccd)
visitList = makeList(visit)
coeffs = estimateCoeffs(visitList, ccdList, threshold=45000, plot=plot, title=title)
if showCoeffs:
printCoeffs(coeffs)
mi = readImage(visitList[0], ccdList[0])
if fixXTalk:
subtractXTalk(mi, coeffs, threshold)
return mi, coeffs