skySim = SkySim()
# Instantiate the telescope
# There is no camera correction for the comcam at this moment
tele = TeleFacade(M1M3=M1M3, M2=M2, phoSimCommu=phoSimCommu)
# Set the configuration file path
tele.setConfigFile(configFilePath)
# Set the subsystem directory
tele.setSubSysConfigFile(M1M3dataDir=M1M3dataDir,
M2dataDir=M2dataDir,
phosimDir=phosimDir)
# Set the comcam camera
tele.setInstName("comcam10")
# Add the star
skySim.addStarByRaDecInDeg(0, 0, 0, 17)
skySim.addStarByRaDecInDeg(1, 0.0033, 0.0033, 18)
# Set the degree of freedom for intra- and extra-focal images
# Double check the direction with Bo
instFileNameList = {"-1": "starExtra.inst", "1": "starIntra.inst"}
logFileNameList = {"-1": "starExtraPhoSim.log", "1": "starIntraPhoSim.log"}
outputImgDirIntra = os.path.join(outputImgDir, "Intra")
outputImgDirExtra = os.path.join(outputImgDir, "Extra")
outputImgDirList = {"-1": outputImgDirExtra, "1": outputImgDirIntra}
argStringList = []
for ii in (-1, 1):
rotSkyPos=cameraRotation,
mjd=mjdTime)
# Instantiate the subsystems
phoSimCommu = PhosimCommu()
skySim = SkySim()
# Instantiate the telescope
# There is no camera correction for the comcam at this moment
tele = TeleFacade(phoSimCommu=phoSimCommu)
# Set the subsystem directory
tele.setSubSysConfigFile(phosimDir=phosimDir)
# Set the comcam camera
tele.setInstName("lsst")
# Add the star
sensorName = ["R22_S11", "R22_S10"]
starMag = [15, 16.5, 15, 15]
xInpixelInCam = [500, 500, 800, 950]
yInPixelInCam = [1000, 1080, 1300, 1500]
starId = 23
for sensor in sensorName:
for ii in range(len(xInpixelInCam)):
skySim.addStarByChipPos(camera, obs, sensor, starId,
xInpixelInCam[ii], yInPixelInCam[ii],
starMag[ii], folderPath2FocalPlane)
starId += 1